TY - JOUR AU - Albrechtová, Jana AB - Abstract Chloroplast number per cell is a frequently examined quantitative anatomical parameter, often estimated by counting chloroplast profiles in two-dimensional (2D) sections of mesophyll cells. However, a mesophyll cell is a three-dimensional (3D) structure and this has to be taken into account when quantifying its internal structure. We compared 2D and 3D approaches to chloroplast counting from different points of view: (i) in practical measurements of mesophyll cells of Norway spruce needles, (ii) in a 3D model of a mesophyll cell with chloroplasts, and (iii) using a theoretical analysis. We applied, for the first time, the stereological method of an optical disector based on counting chloroplasts in stacks of spruce needle optical cross-sections acquired by confocal laser-scanning microscopy. This estimate was compared with counting chloroplast profiles in 2D sections from the same stacks of sections. Comparing practical measurements of mesophyll cells, calculations performed in a 3D model of a cell with chloroplasts as well as a theoretical analysis showed that the 2D approach yielded biased results, while the underestimation could be up to 10-fold. We proved that the frequently used method for counting chloroplasts in a mesophyll cell by counting their profiles in 2D sections did not give correct results. We concluded that the present disector method can be efficiently used for unbiased estimation of chloroplast number per mesophyll cell. This should be the method of choice, especially in coniferous needles and leaves with mesophyll cells with lignified cell walls where maceration methods are difficult or impossible to use. Chloroplast counting, confocal microscopy, disector method, mesophyll, coniferous needle structure, Norway spruce (Picea abies L. Karst.), profile counting, stereology. Introduction Chloroplasts are important organelles of plant photosynthesizing cells as loci where the photosynthetic processes take place. In mesophyll cells, chloroplasts are usually located next to the cytoplasmic membrane adjacent to intercellular spaces to decrease the resistance to CO2 diffusion (Terashima et al., 2011). The chloroplast number per cell represents a frequently examined quantitative anatomical parameter, reflecting various leaf internal and external conditions. It can be influenced by internal factors, such as ploidy (Mochizuki and Sueoka, 1955) and environmental factors, e.g. CO2 concentration (Wang et al., 2004; Teng et al., 2006). Well-established methods for estimation of the chloroplast number per cell previously used in two-dimensional (2D) space are based on: (i) counting chloroplasts in stomatal guard cells in epidermal peels, (ii) determining the number of chloroplasts in flattened mesophyll cells after maceration, and (iii) counting profiles of chloroplasts in leaf sections. Chloroplast counting in guard cells of epidermal peels using a light microscope is not problematic, enabling counting of chloroplasts contained in one cell in one focal plane; this has been applied to various plants, e. g. Beta vulgaris L. (Mochizuki and Sueoka, 1955), Solanum tuberosum L. (Frandsen, 1968), and Citrullus lanatus (Thunb.) Matsum. and Nakai. (Sari et al., 1999). One of the most frequently used methods for estimation of chloroplast number per mesophyll cell in herbaceous plants in 2D is counting chloroplasts in separated mesophyll cells obtained by maceration procedures, as described by Possingham and Saurer (1969). Fixation and maceration of Spinacia oleracea L. leaves enabled the separation of mesophyll cells, which were then flattened into a single plane of focus, allowing the chloroplasts to be viewed in a single layer and thus to be easily counted under phase contrast in an optical microscope with an eyepiece graticule. Methods based on these principles were then improved (Possingham and Smith, 1972), and modified protocols have been widely used, e.g. for Triticum aestivum L. (Boffey et al., 1979), Spinacia oleracea L. (Chaly et al., 1980), Pisum sativum L. (Lamppa et al., 1980), Medicago sativa L. (Molin et al., 1982), B. vulgaris L. (Tymms et al., 1983), and Chenopodium album L. (Yamasaki et al., 1996). The maceration process was further adapted and applied to Arabidopsis thaliana L. and counting was accomplished after chloroplast thresholding using an image analysis program (Pyke and Leech, 1991; Marrison et al., 1999; Stettler et al., 2009). Sung and Chen (1989) extracted chloroplasts from macerated leaves of Glycine max L. Merrill and counted chloroplasts per leaf using a haemocytometer—a special microscope slide with an engraved graticule. To our knowledge, a method for chloroplast counting using maceration has not yet been applied to leaves with thick, lignified cell walls, such as coniferous needles. Another method frequently used for the estimation of chloroplast number per mesophyll cell in 2D is based on counting chloroplast profiles in semi-thin (1–4 μm thick) physical sections of a leaf using transmission electron and light microscopy (Boffey et al., 1979; Zechmann et al., 2003; Wang et al., 2004; Jin et al., 2011; Simon et al., 2013). This method is frequently used, but its practical application including sampling design is often insufficiently described (Sam et al., 2003; Hayashida et al., 2005; Teng et al., 2006; Gopi et al., 2008). Moreover, in many recent studies, the chloroplast and/or cell profiles in 2D were counted to get an estimation of chloroplast number in 3D, which is a theoretically incorrect approach. For example, the number of chloroplasts per unit of leaf area was estimated from the number of chloroplast profiles in a 0.8 μm thick section observed by light microscopy (Miyazawa and Terashima, 2001; Oguchi et al., 2005) and the number of mesophyll cells per leaf area was determined from the number of cell profiles in 5 μm thick sections (Adachi et al., 2013). Some further methods have involved analysis of several optical sections from the same specimen, enabling the counting of chloroplasts in 3D directly—either during focusing through the whole thick physical section by a light microscope (e.g. for T. aestivum L. and Triticum monococcum L.; Ellis and Leech, 1985), focusing through protoplasts in a cell suspension of Marchantia polymorpha L. (Bockers et al., 1997), or in a 3D reconstruction created from images of series of optical sections acquired by confocal microscopy. The latter approach was used for determining the chloroplast number per stomatal guard cell in leaves of potato S. tuberosum L. (Mozafari et al., 1997), A. thaliana L. (Dinkins et al., 2001), Eucalyptus saligna Sm. (Xu et al., 2012), and Glycine dolichocarpa Tateishi and H. Ohashi (Coate et al., 2012). Chloroplasts in mesophyll cells of Oryza sativa L. were counted after comparing photographs from three planes of focus located inside one cell (Hassan and Wazuddin, 2000). The importance of using a proper method for counting/sampling chloroplasts is obvious, because incorrect approaches can lead to serious bias in the chloroplast number estimation. However, in many studies, the sampling method is not well specified. In general, it is important to note that an unbiased method for sampling/counting particles must sample any of the particles with the same probability (Sterio, 1984). It should be stressed that profiles of particles in a 2D section of a specimen do not represent an unbiased sample of the particles, as the larger particles are sampled (i.e. sectioned) with a higher probability. Moreover, it is known that anatomical parameters often exhibit gradients along the leaf blade (e.g. Pazourek, 1966; Possingham and Saurer, 1969). Thus, a proper design of sampling leaf segments for an analysis is crucial to obtain unbiased results. For an unbiased estimation of the number of particles in 3D with no assumption about particle shape and size, the optical disector method based on 3D unbiased sampling probe was developed (Sterio, 1984; Gundersen, 1986). This method enables particles of varying shape, such as chloroplasts in a cell and cells within a tissue, to be sampled interactively and counted in an unbiased way. This method has already been used successfully for estimation of mesophyll cell number (Albrechtová and Kubínová, 1991; Kubínová, 1991, 1993, 1994; Kubínová et al., 2002; Albrechtová et al. 2007), but to our knowledge it has never been used for chloroplast number estimation. We hypothesised that the number of chloroplast profiles counted in cell sections in 2D would be different from the number of chloroplasts counted by an unbiased method in 3D. To verify our hypothesis, we compared 2D and 3D approaches to chloroplast counting (i) in practical measurements of mesophyll cells of Norway spruce needles, (ii) in a 3D model of a mesophyll cell with chloroplasts, and (iii) by a theoretical analysis. Materials and methods Norway spruce (Picea abies L. Karst.) current-year needles were collected from crowns of 18-year-old trees planted on an experimental site of the Global Change Research Centre, Academy of Sciences, Czech Republic, the Bílý Kříž in Moravskoslezské Beskydy mountains, in 2004 (Urban et al., 2001). Spruce needles were collected from south- and south-west-facing branches from the middle crown part and were stored in a deep freeze until processing (Lhotáková et al., 2008). The images were captured during 2009 and analysed during 2009–2012. Sampling of 0.2mm thick needle cross-sections cut off using a hand microtome and placed into a drop of water was done in a systematic uniform random (SUR) way at intervals of 3mm (Fig. 1). Fig. 1. Open in new tabDownload slide Systematic uniform random (SUR) sampling of positions of cross-sections along the needle. 0, tip; z, position of the first cross-section: a random integer from an interval 0–5 was chosen using a table of random numbers; this number determined the position of the first cross-section (z) in the needle (0 corresponded to 0.5mm from the tip and 1 corresponded to 1mm, up to 5, which corresponded to 3mm from the tip), T=3mm (the interval between subsequent cross-sections). (This figure is available in colour at JXBV online.) Fig. 1. Open in new tabDownload slide Systematic uniform random (SUR) sampling of positions of cross-sections along the needle. 0, tip; z, position of the first cross-section: a random integer from an interval 0–5 was chosen using a table of random numbers; this number determined the position of the first cross-section (z) in the needle (0 corresponded to 0.5mm from the tip and 1 corresponded to 1mm, up to 5, which corresponded to 3mm from the tip), T=3mm (the interval between subsequent cross-sections). (This figure is available in colour at JXBV online.) The images of the spruce needle cross-sections were captured using a Leica SP2 AOBS confocal laser-scanning microscope (Leica Microsystems, Wetzlar, Germany) with Ar laser excitation of 488nm; phenolics autofluorescence was detected in the green channel (494–577nm) and autofluorescence of chlorophyll in chloroplasts was detected in the red channel (625–710nm) (Fig. 2A). The acquired images were analysed using Ellipse software (ViDiTo, Košice, Slovakia, http://www.ellipse.sk/, last accessed 21 November 2013), which offers special modules for stereological measurements. Fig. 2. Open in new tabDownload slide Norway spruce needle cross-section: image acquisition, sampling, and processing. (A) Anatomical structure of a Norway spruce needle in a cross-section. Autofluorescence of chlorophyll in chloroplasts was detected in the red channel and autofluorescence of phenolics was detected in the green channel. (B) Histochemical lignin detection in cell walls (pink) of mesophyll cells using a phloroglucinol/HCl test (O’Brian and McCully, 1981). (C) Sampling frames were superimposed on a needle transverse section using the Rectangles module in Ellipse software. (D) Stacks of optical sections were acquired at the positions of rectangles in (C) with higher resolution. (E) Estimation of needle cross-sectional area by a point counting method using the Point Grid module in Ellipse software (F) Two frames showing the subsequent acquisition of series for counting mesophyll cells. (A, C–F), confocal microscopy: (B) bright-field light microscopy. Fig. 2. Open in new tabDownload slide Norway spruce needle cross-section: image acquisition, sampling, and processing. (A) Anatomical structure of a Norway spruce needle in a cross-section. Autofluorescence of chlorophyll in chloroplasts was detected in the red channel and autofluorescence of phenolics was detected in the green channel. (B) Histochemical lignin detection in cell walls (pink) of mesophyll cells using a phloroglucinol/HCl test (O’Brian and McCully, 1981). (C) Sampling frames were superimposed on a needle transverse section using the Rectangles module in Ellipse software. (D) Stacks of optical sections were acquired at the positions of rectangles in (C) with higher resolution. (E) Estimation of needle cross-sectional area by a point counting method using the Point Grid module in Ellipse software (F) Two frames showing the subsequent acquisition of series for counting mesophyll cells. (A, C–F), confocal microscopy: (B) bright-field light microscopy. The Ellipse Disector module (Tomori et al., 2001) enables application of a virtual 3D probe to the stack of serial optical sections. It is possible to browse through the successive optical sections and to mark the selected particles (Figs. 3B–D). The optical disector probe (Sterio, 1984, Gundersen, 1986) represents a virtual 3D block with three exclusion planes (Fig. 3A). The bottom look-up plane belongs to exclusion planes together with two side planes. The reference plane on the top does not belong to the exclusion planes. Particles lying within this block or intersecting its planes, except the exclusion ones, are counted (Fig. 3). The particles are sampled without bias if each particle has the same probability of being sampled, i.e. if the 3D space is filled in by the spatial mosaic of the shifted copies of the disector probe and each particle in the space is sampled by one and only one shifted copy of the 3D sampling probe (Fig. 4). Fig. 3. Open in new tabDownload slide Chloroplast counting using the virtual 3D disector probe. (A) Scheme of the disector probe with chloroplasts. The disector probe is a 3D block. Chloroplasts lying within this block or intersecting its planes are counted, except those intersecting the exclusion planes. The exclusion planes in this scheme are represented by the dark planes and the transparent front plane (all bordered by a full line): the bottom KLMN (so-called look-up plane), KNRO, KLPO, and half-planes PLS and RNT. The top rectangle OPQR is called the reference plane and does not belong to the exclusion planes. Chloroplasts that lie fully inside the probe are always counted. In this scheme, eight chloroplasts are counted (ticked) and six chloroplasts are not counted (crossed). (B–D) Three optical sections of the stack of real serial optical sections acquired by confocal microscopy. Within section numbers 16–27, the disector probe was placed using the Disector module in the Ellipse software, section 16 being the reference plane of the probe (B), section 22 inside the disector probe (C), and section 27 the look-up plane of the disector probe (D). Counted particles are those within the probe not intersecting the exclusion planes: point, counted particle; cross, particle is not counted. In this example, eight chloroplasts were counted (chloroplasts numbers 3 and 10–13 are intersecting the exclusion planes). (This figure is available in colour at JXB online.) Fig. 3. Open in new tabDownload slide Chloroplast counting using the virtual 3D disector probe. (A) Scheme of the disector probe with chloroplasts. The disector probe is a 3D block. Chloroplasts lying within this block or intersecting its planes are counted, except those intersecting the exclusion planes. The exclusion planes in this scheme are represented by the dark planes and the transparent front plane (all bordered by a full line): the bottom KLMN (so-called look-up plane), KNRO, KLPO, and half-planes PLS and RNT. The top rectangle OPQR is called the reference plane and does not belong to the exclusion planes. Chloroplasts that lie fully inside the probe are always counted. In this scheme, eight chloroplasts are counted (ticked) and six chloroplasts are not counted (crossed). (B–D) Three optical sections of the stack of real serial optical sections acquired by confocal microscopy. Within section numbers 16–27, the disector probe was placed using the Disector module in the Ellipse software, section 16 being the reference plane of the probe (B), section 22 inside the disector probe (C), and section 27 the look-up plane of the disector probe (D). Counted particles are those within the probe not intersecting the exclusion planes: point, counted particle; cross, particle is not counted. In this example, eight chloroplasts were counted (chloroplasts numbers 3 and 10–13 are intersecting the exclusion planes). (This figure is available in colour at JXB online.) Fig. 4. Open in new tabDownload slide Scheme of the hypothetic spatial mosaic of 3D disector sampling probes illustrating that particles are unambiguously sampled by only one disector sampling probe. The probes a, b, and c are located next to each other. The exclusion planes in this scheme are represented by grey planes. Particles are labelled according to the label of the probe they are counted in. In this scheme, in probe a eight chloroplasts are counted (labelled a), in probe b eight chloroplasts are counted (labelled b), and in probe c nine chloroplasts are counted (labelled c). Chloroplasts labelled by crosses are not counted by any of the three probes shown here, but they would be unambiguously sampled by adjacent hypothetical probes. (This figure is available in colour at JXB online.) Fig. 4. Open in new tabDownload slide Scheme of the hypothetic spatial mosaic of 3D disector sampling probes illustrating that particles are unambiguously sampled by only one disector sampling probe. The probes a, b, and c are located next to each other. The exclusion planes in this scheme are represented by grey planes. Particles are labelled according to the label of the probe they are counted in. In this scheme, in probe a eight chloroplasts are counted (labelled a), in probe b eight chloroplasts are counted (labelled b), and in probe c nine chloroplasts are counted (labelled c). Chloroplasts labelled by crosses are not counted by any of the three probes shown here, but they would be unambiguously sampled by adjacent hypothetical probes. (This figure is available in colour at JXB online.) At first, the image of an entire cross-section was captured using a dry plan apochromatic 10× objective [numeric aperture (NA) 0.4] with resolution 1024×1024 pixels (1500×1500 µm). From these images, the areas of needle cross-sections and the proportions of mesophyll were estimated by the point counting method (Weibel, 1979) using the Ellipse Point Grid software module (Fig. 2E). The spruce needle volume was estimated by Cavalieri principle (Gundersen and Jensen, 1987): the mean of the areas of cross-sections was multiplied by the length of the needle. For mesophyll cell counting, 30 µm thick stacks comprising 16 serial optical sections 2 µm apart were acquired by using a plan apochromatic 20× water-immersion objective (NA 0.7) with resolution 2048×2048 pixels (750×750 µm) (Fig. 2F, Supplementary Videos S1 and Supplementary Data, at JXB online). Several stacks were acquired in each cross-section to cover its whole area; altogether, 157 stacks were acquired. The height of the disector probe was 10 μm, placed in the middle of captured stacks, i.e. within a substack of six optical sections 2 μm apart, the size of the sampling frame was adjusted to cover the whole cross-section area. A double disector (Gundersen, 1986) was used: the probe was used in both directions, i.e. at first, cells were counted by the disector with the reference plane and look-up planes in the positions shown in Fig. 3A, and then additional cells were counted by the disector with the reference plane set to the position of the former look-up plane and vice versa. Thus, the total height of the probe was twice its real height (20 μm, in this case). The needle volume in the stack was determined by multiplying the height of the double disector and the needle area estimated from the cross-section captured by a 10× objective. The mesophyll cell density (number of cells per needle volume) was estimated by the formula (Sterio, 1984; Gundersen, 1986): where estNVneedle(cell) is the estimated number of mesophyll cells per needle volume, is the sum of all sampled cells in all disector probes within a needle, ∑Pi is the sum of all points falling within a needle in all disector probes (for calculation needle volume in each probe) used for cell counting, p is the number of test points in a grid in a sampling frame used for cell counting, a is the area of the disector sampling frame (the base plane of the 3D probe) used for cell counting, and h is the height of the disector probe used for cell counting For chloroplast counting, the positions of captured stacks were selected by random application of rectangle grid in the Rectangles module of the Ellipse software (ViDiTo, Košice, SR) (Fig. 2C) on the images captured by a 10× objective, while the size of the rectangles was set equal to the size of the field of view and the distance between the rectangles was set to capture on average 24 stacks of optical sections per spruce needle. In each position, 20 µm thick stacks comprising 41 serial optical sections 0.5 µm apart were acquired using a plan apochromatic 63× water-immersion objective (NA 1.2) with zoom 2× and resolution 512×512 pixels (119×119 µm) (Fig. 2D). Altogether, 464 stacks were acquired and used for counting chloroplasts by applying the disector method. The dimensions of the disector sampling frame for chloroplast counting in the middle of the stack were set to 500×500 pixels (116×116 μm) and the height of the probe was 5.5 μm, as there were 12 optical sections included, 0.5 μm apart. To estimate the volume of the spruce needle within each disector probe, the needle area in the middle image of the stack was determined by the point-grid method and multiplied by the disector height. The chloroplast density (number of chloroplasts per needle volume) was estimated by the formula (Sterio, 1984; Gundersen, 1986): where estNVneedle(chl) is the estimated number of chloroplasts per needle volume, is the sum of all sampled chloroplasts in all disector probes within a needle, ∑P’i is the sum of all points falling within a needle in all 3D probes (for calculation of needle volume in each probe) used for chloroplast counting, p’ is the number of test points in a grid in a sampling frame used for chloroplast counting, a’ is the area of the disector sampling frame (the base plane of the 3D probe) used for chloroplast counting, and h’ is the height of the disector probe used for chloroplast counting. The estimated number of chloroplasts per mesophyll cell was calculated by the ratio of the chloroplast number and the mesophyll cell number per needle volume: where estNNcell(chl) is the estimated number of chloroplasts per mesophyll cell, estNVneedle(cell) is the estimated number of mesophyll cells per needle volume, and estNVneedle(chl) is the estimated number of chloroplasts per needle volume. The profiles of chloroplasts and cells in 2D images were counted in the same stacks as those used for 3D measurements in the middle optical sections using the 2D unbiased sampling frame (Gundersen, 1977) (Fig. 2D). The number of chloroplasts per cell was calculated as the ratio of the chloroplast profile number per needle area to the cell profile number per needle area, while the needle area was estimated previously by the point counting method (Weibel, 1979). In order to obtain comparable results by 2D and 3D methods, we calculated the chloroplast number per cell in each needle cross-section (68 in total) by both methods. In the 3D method, we summed the number of chloroplasts per volume as estimated by the disector method in each cross-section and divided it by the sum of the estimated number of cells in each cross-section. In the 2D profile counting method, we summed the estimated number of chloroplasts per area in each cross-section and divided it by the sum of the estimated number of cells in each cross-section. Acquired values were compared using statistical analysis (paired t-test) in NCSS 2000 program (Number Cruncher Statistical Systems, Kaysville, Utah, USA). For visualization of mesophyll cell surface (Fig. 5A–E) and 3D chloroplast arrangement in a mesophyll cell (Fig. 5G), the 3D reconstructions were created based on processing series of optical sections 2 µm apart acquired by a confocal microscope using a 20× objective. A model of a simplified mesophyll cell with 210 chloroplasts was made in the Cortona software (Fig. 5F, H, Supplementary Video S3 at JXB online) in order to demonstrate comparison of the 3D disector and the 2D profile counting methods. Fig. 5. Open in new tabDownload slide 3D reconstruction and 3D model of a Norway spruce mesophyll cell. (A–D) Different optical sections of a mesophyll cell in a stack of images captured by a confocal microscope by using an objective 20×. (E) 3D reconstruction of the surface of the mesophyll cell from the confocal image stack shown in (A)–(D) using Ellipse software. (F) Face view of a model of a mesophyll cell with 210 chloroplasts made in Cortona 3D software. (G) 3D reconstruction of the mesophyll cell with chloroplasts from a stack of optical sections from confocal microscopy using Ellipse software. The positions of chloroplasts do not correspond to the in vivo state. (H) Side view of the model from (F) with six lines indicating planes of cross-sections. (I) 2D sections through the model cell shown in (F) and (H). Letters indicate the position of section planes in (H). Numbers indicate the number of chloroplast profiles in each of these sections. Fig. 5. Open in new tabDownload slide 3D reconstruction and 3D model of a Norway spruce mesophyll cell. (A–D) Different optical sections of a mesophyll cell in a stack of images captured by a confocal microscope by using an objective 20×. (E) 3D reconstruction of the surface of the mesophyll cell from the confocal image stack shown in (A)–(D) using Ellipse software. (F) Face view of a model of a mesophyll cell with 210 chloroplasts made in Cortona 3D software. (G) 3D reconstruction of the mesophyll cell with chloroplasts from a stack of optical sections from confocal microscopy using Ellipse software. The positions of chloroplasts do not correspond to the in vivo state. (H) Side view of the model from (F) with six lines indicating planes of cross-sections. (I) 2D sections through the model cell shown in (F) and (H). Letters indicate the position of section planes in (H). Numbers indicate the number of chloroplast profiles in each of these sections. For maceration of needle mesophyll samples, in order to count chloroplasts directly in separated cells, different methods were tested: a 1M aqueous solution of HCl (Possingham and Saurer, 1969); a 3.5% aqueous solution of glutaraldehyde and 0.1M aqueous solution of Na2EDTA (Boffey et al., 1979); and a solution of 10% aqueous CrO3 and 10% aqueous HNO3 (O’Brian and McCully, 1981). A phloroglucinol/HCl test (O’Brian and McCully, 1981) was applied to demonstrate lignification of mesophyll cell walls (Fig. 2B). Results The Norway spruce needle structure (Fig. 2A) is composed of a one-cell-layer epidermis with highly cutinized periclinal cell walls on the needle surface; below the epidermis, there is one-cell layer of sclerenchymatic hypodermis not present below the stomata, followed by mesophyll and a central cylinder with a vascular bundle surrounded by transfusion tissue. Mesophyll cells contain chloroplasts and have an irregular shape with lobed anticlinal cell walls and are prolonged in the direction perpendicular to the needle surface (Fig. 5E), forming layers of tightly connected mesophyll cells surrounded by intercellular spaces (Supplementary Videos S1 and Supplementary Data, Supplementary Fig. S1 at JXB online). Numerous maceration methods were tested in order to count chloroplasts directly in separated spruce mesophyll cells. However, the spruce mesophyll cells kept connected together by the middle lamella and cell walls. We concluded that maceration methods were unsuccessful due to lignification of mesophyll cell walls detected histochemically (Fig. 2B) and observed regularly in our previous study (Soukupová et al., 2000). The average Norway spruce needle volume was 8.17±1.00mm3 ([mean ± standard error (SE)] and the proportion of mesophyll in the needle volume was 72.0±0.8 %. The average number of chloroplasts per mesophyll cell estimated by profile counting in 2D was 20.92±1.32, while the average number of chloroplasts estimated by the disector probe in 3D was 10 times higher, i.e. 209.65±17.44 (Fig. 6). Unsurprisingly, statistical analysis showed a significant difference between results obtained in 2D and 3D. The methods of chloroplast counting in 2D and 3D thus yielded estimates of the chloroplast number per mesophyll cell that were different by one order of magnitude. Fig. 6. Open in new tabDownload slide Comparison of 3D- and 2D-based methods for chloroplast number estimation. The mean number of chloroplasts per mesophyll cell in Norway spruce needle estimated by 3D disector was 209.65±17.44 (mean±SE) (open column), and by profile counting in 2D (filled column) was 20.92±1.32. This result was statistically significantly different (P<0.001). Fig. 6. Open in new tabDownload slide Comparison of 3D- and 2D-based methods for chloroplast number estimation. The mean number of chloroplasts per mesophyll cell in Norway spruce needle estimated by 3D disector was 209.65±17.44 (mean±SE) (open column), and by profile counting in 2D (filled column) was 20.92±1.32. This result was statistically significantly different (P<0.001). To visualize differences in the results yielded by both methods, we used the 3D model of a mesophyll cell (Fig. 5F, H, I, Supplementary Video S3). We made 111 sections through the cell model and counted profiles in each section. The mean number of particle profiles per section was 17.23±0.69 (mean±SE). However, there were 210 particles in the model cell. Our test showed that the average number of chloroplasts estimated by profile counting was more than 10 times underestimated in comparison with the real number of chloroplasts in the model. This is clearly demonstrated in Fig. 6. The difference between both estimators tested can be also explained theoretically as shown below. The number of chloroplasts per mesophyll cell can be calculated by: Instead of direct counting of N, it could be estimated by estimations of densities : obtained by stereological methods. should be estimated with a small error, , because it is not possible to divide by zero, and with an error of about (Marsaglia, 2006). The number of chloroplasts cannot be directly estimated by counting profiles in the section, because the areal density of profiles is related not only to density Nv but also to height h: . If we count the ratio of estimations, we get close to the value which is different from the desired value by the ratio of the height of chloroplasts to the height of mesophyll cells: The chloroplasts are much smaller than mesophyll cells (chloroplast mean height is approximately 4 μm, while the cell mean height is approximately 40–60 μm); therefore, the theoretically expected difference is high, also of one order of magnitude. Discussion Counting chloroplast profiles in 2D yielded values of chloroplast number per cell that were 10 times underestimated in comparison with the mean number of chloroplasts estimated by the unbiased disector method in 3D. We determined that the systematic shift in an estimated value is inversely related to the ratio of the height of chloroplasts to the height of mesophyll cells. The 3D model also showed that neglecting the 3D appearance of the cell led to the underestimation of the number of chloroplasts (Fig. 6). All our practical and theoretical tests thus clearly showed that the frequently used method for chloroplast number estimation by counting profiles of particles from 2D sections (e. g. Boffey et al., 1979; Sam et al., 2003; Wang et al., 2004; Hayashida et al., 2005; Teng et al., 2006; Gopi et al., 2008) yielded biased estimates and that the results may be one order of magnitude different from the real chloroplast numbers. A review of the previously reported results on the number of chloroplasts per mesophyll cell with focus on the method used is presented in Supplementary Table S1 at JXB online. The studies were conducted on various plant species, considering herbaceous species in the majority of cases. For example, the extreme variability can be seen in the number of chloroplasts in Arabidopsis mesophyll cells estimated by various methods. The number of chloroplasts varied from eight to 10 chloroplasts per cell, based on counting chloroplast profiles in thin sections (Teng et al., 2006; Jin et al., 2011), to more than 100 chloroplasts per cell if the number was determined by the maceration method (Pyke and Leech, 1992; Marrison et al., 1999), or even more than 200 (Meyer et al., 2006). It is generally known that the number of chloroplasts depends on plant internal factors—developmental stage, age, and species genotype (reviewed by Possingham, 1980)—or environmental factors, although the methodical bias could conceal effects of these factors. Authors often do not describe the sampling design and the method used for chloroplast counting. Nevertheless, sampling design is crucial for yielding unbiased results, as anatomical gradients within a leaf blade were demonstrated also for the number of chloroplasts per cell (Dean and Leech, 1982). Gradients in leaf anatomical parameters, e.g. the number of chloroplasts per cell, also exist in dependence on distance of a leaf from the root system (Possingham and Saurer, 1969). The chloroplast number per cell in conifers was estimated by profile counting to be six to 16 chloroplasts per cell in different species sampled in August (Maslova et al., 2009) (Supplementary Table S1). These numbers seem to be very low and may not reflect the real chloroplast number in a mesophyll cell. For example, the number of chloroplasts per pair of stomatal guard cells in epidermal peels varies between 10 and 20 (Mochizuki and Sueoka, 1955; Frandsen, 1968; Nicholson, 1981; Qin and Rotino, 1995; Sari et al., 1999) and the mesophyll cells are obviously larger than guard cells. In flat cells, such as the stomatal guard cells, it is possible to observe the chloroplasts practically in one optical section and to count them correctly in contrary to the mesophyll cells, which are larger and thicker, so it is unlikely that profiles of all of the chloroplasts within a cell would be present in a single optical section. Chloroplast counting in separated mesophyll cells after leaf maceration (e.g. Possingham and Saurer, 1969; Stettler et al., 2009) is appropriate in the case when the chloroplasts in the specimen are not overlapping. Counting of chloroplasts in a solution from the macerated leaf segments (Sung and Chen, 1989) does not enable the determination of chloroplast number per cell as the number of mesophyll cells, from which the chloroplasts present in a solution are released, is not known. Applying the maceration method to coniferous needles is problematic or even impossible, as shown in our study, as they contain phenolic compounds and lignin in cell walls (Fig. 2B), as proved in a histochemical study by Soukupová et al. (2000). Counting chloroplasts in cells directly during focusing through the specimen using conventional light microscopy (Ellis and Leech, 1985; Bockers et al., 1997) can be applied if entire cells can be focused through and the chloroplasts are sparsely distributed in cells; however, chloroplasts usually tend to be densely packed along the cytoplasmic membrane. Determining the chloroplast number per cell in 3D reconstructions made from a series of confocal microscope images (Mozafari et al., 1997; Dinkins et al., 2001; Coate et al., 2012; Xu et al., 2012) can yield an unbiased estimate if SUR sampling is applied and a sufficient number of cells is analysed. However, this is a much more time-consuming approach than application of the disector method. In order to get unbiased results, it is also important to apply unbiased sampling of locations within the leaf where the measurement is performed. Many authors of previous studies reporting the chloroplast number per mesophyll cell do not explain the sampling design in a sufficient detail. However, a SUR sampling, along with both assumptions—the unbiased method for counting particles and the sufficient number of measurements—are essential for getting unbiased and accurate results (Sterio, 1984). In the present study, sample preparation and image acquisition were carefully designed according to the principle of SUR sampling (Gundersen and Jensen, 1987). Needles collected from Norway spruce trees were stored frozen before processing; freezing of needles does not influence various geometrical parameters of mesophyll (Lhotáková et al., 2008). Similarly, the number of mesophyll cells and the number of chloroplasts cannot change during the needle storage. Norway spruce chloroplasts kept their chlorophyll fluorescence activity, even after 5 years of storage in a freezer, and were perfectly suitable for confocal microscopy. To eliminate the possible shrinkage of a tissue caused by cutting, the dry stem pith of Sambucus nigra L. was used to fix the sample during cutting. It was found that some cells might have been pulled out of the specimen during cutting (Lhotáková et al., 2008), and this is why the disector probe was placed in the middle of a stack. In general, plant mesophyll cells contain a large central vacuole and, in vivo, chloroplasts are usually located in the peripheral cytoplasm, close to the cell wall. However, the spatial distribution of organelles may be disrupted by freezing. Therefore, we are aware that our 3D models of mesophyll cells (Fig. 5) do not correspond to in vivo chloroplast arrangement; however, they illustrate well the size of Norway spruce mesophyll cells, their variable and irregular shape, as well as the number and size of chloroplasts within a cell. The particular mesophyll cell chosen for 3D reconstruction was a smaller one, so that it was possible to focus through the entire cell. Most mesophyll cells were too high to be focused through; therefore, a method based on direct counting of chloroplasts within the entire cell could not be used. The present 3D model of mesophyll cell (Fig. 5F, H, I) was constructed such that it contained the mean number of chloroplasts obtained in our study. Conclusions Based on the above-presented lines of evidence, we conclude that the presented optical disector method for chloroplast counting, using stacks of confocal microscopic images acquired from thick tissue sections accompanied by SUR sampling, is a very efficient method for estimating chloroplast number per mesophyll cell. While using the optical disector method for chloroplast counting, an important consideration that a mesophyll cell is a 3D structure is taken into account. We propose this method as a universal unbiased one. It should be the method of choice, especially in conifers or other xeromorphic leaves with mesophyll cells with lignified walls, where maceration methods are difficult or impossible to apply. Abbreviations: Abbreviations: 2D two-dimensional 3D three-dimensional SUR systematic uniform random. Acknowledgements We are grateful to Drahomíra Bartáková for technical assistance, Barbora Radochová for assistance with confocal microscopy, and Marianna Dubánková for help with sample collection. This work was supported by the Czech Science Foundation (P501/10/0340), funds provided by the Academy of Sciences of the Czech Republic (AV0Z50110509 and RVO:67985823), and by the Charles University in Prague (SVV 265203). References Adachi S Nakae T Uchida M et al. 2013 . The mesophyll anatomy enhancing CO2 diffusion is a key trait for improving rice photosynthesis . Journal of Experimental Botany 64 , 1061 – 1072 . Google Scholar Crossref Search ADS PubMed WorldCat Albrechtová J Kubínová L . 1991 . Quantitative analysis of the structure of etiolated barley leaf using stereological methods . Journal of Experimental Botany 42 , 1311 – 1314 . Google Scholar Crossref Search ADS WorldCat Albrechtová J Janáček J Lhotáková Z Radochová B Kubínová L . 2007 . Novel efficient methods for measuring mesophyll anatomical characteristics from fresh thick sections using stereology and confocal microscopy: application on acid rain-treated Norway spruce needles . Journal of Experimental Botany 58 , 1451 – 1461 . Google Scholar Crossref Search ADS PubMed WorldCat Bockers M Čapková V Tichá I Schafer C . 1997 . Growth at high CO2 affects the chloroplast number but not the photosynthetic efficiency of photoautotrophic Marchantia polymorpha culture cells . Plant Cell Tissue and Organ Culture 48 , 103 – 110 . Google Scholar Crossref Search ADS WorldCat Boffey SA Ellis JR Sellden G Leech RM . 1979 . Chloroplast division and DNA synthesis in light-grown wheat leaves . Plant Physiology 64 , 502 – 505 . Google Scholar Crossref Search ADS PubMed WorldCat Chaly N Possingham JV Thomson WW . 1980 . Chloroplast division in spinach leaves examined by scanning electron-microscopy and freeze-etching . Journal of Cell Science 46 , 87 – 96 . Google Scholar PubMed OpenURL Placeholder Text WorldCat Coate JE Luciano AK Seralathan V Minchew KJ Owens TG Doyle JJ . 2012 . Anatomical, biochemical, and photosynthetic responses to recent allopolyploidy in Glycine dolichocarpa (Fabaceae) . American Journal of Botany 99 , 55 – 67 . Google Scholar Crossref Search ADS PubMed WorldCat Dean C Leech RM . 1982 . Genome expression during normal leaf development.1. Cellular and chloroplast numbers and DNA, RNA, and protein-levels in tissues of different ages within a 7-day-old wheat leaf . Plant Physiology 69 , 904 – 910 . Google Scholar Crossref Search ADS PubMed WorldCat Dinkins R. Srinivasa Reddy MS Leng M Collins GB . 2001 . Overexpression of the Arabidopsis thaliana MinD1 gene alters chloroplast size and number in transgenic tobacco plants . Planta 214 , 180 – 188 . Google Scholar Crossref Search ADS PubMed WorldCat Ellis JR Leech RM . 1985 . Cell-size and chloroplast size in relation to chloroplast replication in light-grown wheat leaves . Planta 165 , 120 – 125 . Google Scholar Crossref Search ADS PubMed WorldCat Frandsen N . 1968 . Die Plastidenzahl als Merkmal bei der Kartoffel . Theoretical and Applied Genetics 38 , 153 – 167 . Google Scholar Crossref Search ADS PubMed WorldCat Gopi R Jaleel CA Panneerselvam R . 2008 . Leaf anatomical responses of Amorphophallus campanulatus to triazoles fungicides . EurAsian Journal of BioSciences 2 , 46 – 52 . Google Scholar OpenURL Placeholder Text WorldCat Gundersen HJG . 1986 . Stereology of arbitrary particles—a review of unbiased number and size estimators and the presentation of some new ones, in memory of Thompson, William, R . Journal of Microscopy 143 , 3 – 45 . Google Scholar Crossref Search ADS PubMed WorldCat Gundersen HJG Jensen EB . 1987 . The efficiency of systematic sampling in stereology and its prediction . Journal of Microscopy 147 , 229 – 263 . Google Scholar Crossref Search ADS PubMed WorldCat Gundersen HJG . 1977 . Notes on the estimation of the numerical density of arbitrary profiles: the edge effect . Journal of Microscopy 111 , 219 – 223 . Google Scholar Crossref Search ADS WorldCat Hassan L Wazuddin M . 2000 . Colchicine-induced variation of cell size and chloroplast number in leaf mesophyll of rice . Plant Breeding 119 , 531 – 533 . Google Scholar Crossref Search ADS WorldCat Hayashida A Takechi K Sugiyama M Kubo M Itoh RD Takio S Fujita T Hiwatashi Y Hasebe M Takano H . 2005 . Isolation of mutant lines with decreased numbers of chloroplasts per cell from a tagged mutant library of the moss Physcomitrella patens . Plant Biology 7 , 300 – 306 . Google Scholar Crossref Search ADS PubMed WorldCat Jin B Wang L Wang J Jiang KZ Wang Y Jiang XX Ni CY Wang YL Teng NJ . 2011 . The effect of experimental warming on leaf functional traits, leaf structure and leaf biochemistry in Arabidopsis thaliana . BMC Plant Biology 11 , 10 . Google Scholar Crossref Search ADS PubMed WorldCat Kubínová L . 1991 . Stomata and mesophyll characteristics of barley leaf as affected by light—stereological analysis . Journal of Experimental Botany 42 , 995 – 1001 . Google Scholar Crossref Search ADS WorldCat Kubínová L . 1993 . Recent stereological methods for the measurement of leaf anatomical characteristics: estimation of volume density, volume and surface area . Journal of Experimental Botany 44 , 165 – 173 . Google Scholar Crossref Search ADS WorldCat Kubínová L . 1994 . Recent stereological methods for measuring leaf anatomical characteristics: estimation of the number and sizes of stomata and mesophyll cells . Journal of Experimental Botany 45 , 119 – 127 . Google Scholar Crossref Search ADS WorldCat Kubínová L Janáček J Krekule I . 2002 . Stereological methods for estimating geometrical parameters of microscopical structure studied by three-dimensional microscopical techniques . In: Diaspro A, ed. Confocal and two-photon microscopy , New York : Wiley-Liss . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Lamppa GK Elliot LV Bendich AJ . 1980 . Changes in chloroplast number during pea leaf development—an analysis of a protoplast population . Planta 148 , 437 – 443 . Google Scholar Crossref Search ADS PubMed WorldCat Lhotáková Z Albrechtová J Janáček J Kubínová L . 2008 . Advantages and pitfalls of using free-hand sections of frozen needles for three-dimensional analysis of mesophyll by stereology and confocal microscopy . Journal of Microscopy 232 , 56 – 63 . Google Scholar Crossref Search ADS PubMed WorldCat Marrison JL Rutherford SM Robertson EJ Lister C Dean C Leech RM . 1999 . The distinctive roles of five different ARC genes in the chloroplast division process in Arabidopsis . Plant Journal 18 , 651 – 662 . Google Scholar Crossref Search ADS PubMed WorldCat Marsaglia G . 2006 . Ratios of normal variables . Journal of Statistical Software 16 , 1 – 10 . Google Scholar Crossref Search ADS WorldCat Maslova TG Mamushina NS Sherstneva OA Bubolo LS Zubkova EK . 2009 . Seasonal structural and functional changes in the photosynthetic apparatus of evergreen conifers . Russian Journal of Plant Physiology 56 , 607 – 615 . Google Scholar Crossref Search ADS WorldCat Meyer R Yuan J Afzal J Iqbal MJ Zhu MX Garvey G Lightfoot DA . 2006 . Identification of Gsr1 in Arabidopsis thaliana: a locus inferred to regulate gene expression in response to exogenous glutamine . Euphytica 151 , 291 – 302 . Google Scholar Crossref Search ADS WorldCat Miyazawa SI Terashima I . 2001 . Slow development of leaf photosynthesis in an evergreen broad-leaved tree, Castanopsis sieboldii: relationships between leaf anatomical characteristics and photosynthetic rate . Plant, Cell & Environment 24 , 279 – 291 . Google Scholar Crossref Search ADS WorldCat Mochizuki A Sueoka N . 1955 . Genetic studies on the number of plastid in stomata. I. Effect of autopolyploidy in sugar beets . Cytologia 20 , 358 – 366 . Google Scholar Crossref Search ADS WorldCat Molin WT Meyers SP Baer GR Schrader LE . 1982 . Ploidy effects in isogenic populations of alfalfa. 2. Photosynthesis, chloroplast number, ribulose-1,5-bisphosphate carboxylase, chlorophyll, and DNA in protoplasts . Plant Physiology 70 , 1710 – 1714 . Google Scholar Crossref Search ADS PubMed WorldCat Mozafari J Wolyn DJ AliKhan ST . 1997 . Chromosome doubling via tuber disc culture in dihaploid potato as determined by confocal microscopy . Plant Cell Reports 16 , 329 – 333 . Google Scholar OpenURL Placeholder Text WorldCat Nicholson GG . 1981 . The use of chloroplast numbers in guard cells as a means of distinguishing the chromosome races of Ranunculus ficaria L . Annals of Botany 48 , 909 – 913 . Google Scholar Crossref Search ADS WorldCat O’Brian TP McCully ME . 1981 . The study of plant structure: principles and selected methods . Melbourne, Australia : Termarcarphi Pty Ltd . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Oguchi R Hikosaka K Hirose T . 2005 . Leaf anatomy as a constraint for photosynthetic acclimation: differential responses in leaf anatomy to increasing growth irradiance among three deciduous trees . Plant, Cell & Environment 28 , 916 – 927 . Google Scholar Crossref Search ADS WorldCat Pazourek J . 1966 . Anatomical gradients . Acta Universitatis Carolinae—Biologica , Suppl. 1/2 , 19 – 25 . Google Scholar OpenURL Placeholder Text WorldCat Possingham JV Saurer W . 1969 . Changes in chloroplast number per cell during leaf development in spinach (Spinacea oleracea) . Planta 86 , 186 – 194 . Google Scholar Crossref Search ADS PubMed WorldCat Possingham JV Smith JW . 1972 . Factors affecting chloroplast replication in spinach . Journal of Experimental Botany 23 , 1050 – 1059 . Google Scholar Crossref Search ADS WorldCat Possingham JV . 1980 . Plastid replication and development in the life-cycle of higher-plants . Annual Review of Plant Physiology and Plant Molecular Biology 31 , 113 – 129 . Google Scholar Crossref Search ADS WorldCat Pyke KA Leech RM . 1991 . Rapid image analysis screening procedure for identifying chloroplast number mutants in mesophyll cells of Arabidopsis thaliana (L.) Heynh . Plant Physiology 96 , 1193 – 1195 . Google Scholar Crossref Search ADS PubMed WorldCat Pyke KA Leech RM . 1992 . Chloroplast division and expansion is radically altered by nuclear mutations in Arabidopsis thaliana . Plant Physiology 99 , 1005 – 1008 . Google Scholar Crossref Search ADS PubMed WorldCat Qin X Rotino GL . 1995 . Chloroplast number in guard cells as ploidy indicator of in vitro-grown androgenic pepper plantlets . Plant Cell Tissue and Organ Culture 41 , 145 – 149 . Google Scholar Crossref Search ADS WorldCat Sam O Ramirez C Coronado MJ Testillano PS Risueno MC . 2003 . Changes in tomato leaves induced by NaCl stress: leaf organization and cell ultrastructure . Biologia Plantarum 47 , 361 – 366 . Google Scholar Crossref Search ADS WorldCat Sari N Abak K Pitrat M . 1999 . Comparison of ploidy level screening methods in watermelon: Citrullus lanatus (Thunb.) Matsum. and Nakai . Scientia Horticulturae 82 , 265 – 277 . Google Scholar Crossref Search ADS WorldCat Simon UK Polanschütz LM Koffler BE Zechmann B . 2013 . High resolution imaging of temporal and spatial changes of subcellular ascorbate, glutathione and H2O2 Distribution during Botrytis cinerea infection in Arabidopsis . PLoS ONE 8 , 1 – 11 . Google Scholar OpenURL Placeholder Text WorldCat Soukupová J Cvikrová M Albrechtová J Rock BN Eder J . 2000 . Histochemical and biochemical approaches to the study of phenolic compounds and peroxidases in needles of Norway spruce (Picea abies) . New Phytologist 146 , 403 – 414 . Google Scholar Crossref Search ADS WorldCat Sterio DC . 1984 . The unbiased estimation of number and sizes of arbitrary particles using the disector . Journal of Microscopy 134 , 127 – 136 . Google Scholar Crossref Search ADS PubMed WorldCat Stettler M Eicke S Mettler T Messerli G Hortensteiner S Zeeman SC . 2009 . Blocking the metabolism of starch breakdown products in Arabidopsis leaves triggers chloroplast degradation . Molecular Plant 2 , 1233 – 1246 . Google Scholar Crossref Search ADS PubMed WorldCat Sung FJM Chen JJ . 1989 . Changes in photosynthesis and other chloroplast traits in lanceolate leaflet isoline of soybean . Plant Physiology 90 , 773 – 777 . Google Scholar Crossref Search ADS PubMed WorldCat Teng NJ Wang J Chen T Wu XQ Wang YH Lin JX . 2006 . Elevated CO2 induces physiological, biochemical and structural changes in leaves of Arabidopsis thaliana . New Phytologist 172 , 92 – 103 . Google Scholar Crossref Search ADS PubMed WorldCat Terashima I Hanba YT Tholen D Niinemets U . 2011 . Leaf functional anatomy in relation to photosynthesis . Plant Physiology 155 , 108 – 116 . Google Scholar Crossref Search ADS PubMed WorldCat Tomori Z Krekule I Kubínová L . 2001 . DISECTOR program for unbiased estimation of particle number, numerical density and mean volume . Image Analysis & Stereology 20 , 119 – 130 . Google Scholar Crossref Search ADS WorldCat Tymms MJ Scott NS Possingham JV . 1983 . DNA content of Beta vulgaris chloroplasts during leaf cell expansion . Plant Physiology 71 , 785 – 788 . Google Scholar Crossref Search ADS PubMed WorldCat Urban O Janouš D Pokorný R Marková I Pavelka M Fojtík Z Šprtová M Kalina J Marek MV . 2001 . Glass domes with adjustable windows: a novel technique for exposing juvenile forest stands to elevated CO2 concentration . Photosynthetica 39 , 395 – 401 . Google Scholar Crossref Search ADS WorldCat Wang XZ Anderson OR Griffin KL . 2004 . Chloroplast numbers, mitochondrion numbers and carbon assimilation physiology of Nicotiana sylvestris as affected by CO2 concentration . Environmental and Experimental Botany 51 , 21 – 31 . Google Scholar Crossref Search ADS WorldCat Weibel ER . 1979 . Practical methods for biological morphometry . In: Stereological methods , vol. 1 . New York : Academic Press , 1 – 415 . Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC Xu CY Salih A Ghannoum O Tissue DT . 2012 . Leaf structural characteristics are less important than leaf chemical properties in determining the response of leaf mass per area and photosynthesis of Eucalyptus saligna to industrial-age changes in [CO2] and temperature . Journal of Experimental Botany 63 , 5829 – 5841 . Google Scholar Crossref Search ADS PubMed WorldCat Yamasaki T Kudoh T Kamimura Y Katoh S . 1996 . A vertical gradient of the chloroplast abundance among leaves of Chenopodium album . Plant and Cell Physiology 37 , 43 – 48 . Google Scholar Crossref Search ADS WorldCat Zechmann B Müller M Zellnig G . 2003 . Cytological modifications in zucchini yellow mosaic virus (ZYMV)-infected Styrian pumpkin plants . Archives of Virology 148 , 1119 – 1133 . Google Scholar Crossref Search ADS PubMed WorldCat © The Author 2013. Published by Oxford University Press on behalf of the Society for Experimental Biology. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. © The Author 2013. Published by Oxford University Press on behalf of the Society for Experimental Biology. TI - Unbiased estimation of chloroplast number in mesophyll cells: advantage of a genuine three-dimensional approach JF - Journal of Experimental Botany DO - 10.1093/jxb/ert407 DA - 2014-02-01 UR - https://www.deepdyve.com/lp/oxford-university-press/unbiased-estimation-of-chloroplast-number-in-mesophyll-cells-advantage-Z5zrtHX5GR SP - 609 EP - 620 VL - 65 IS - 2 DP - DeepDyve ER -