- Open Access
A novel method to characterize silica bodies in grasses
Plant Methodsvolume 12, Article number: 3 (2016)
The deposition of silicon into epidermal cells of grass species is thought to be an important mechanism that plants use as a defense against pests and environmental stresses. There are a number of techniques available to study the size, density and distribution pattern of silica bodies in grass leaves. However, none of those techniques can provide a high-throughput analysis, especially for a great number of samples.
We developed a method utilizing the autofluorescence of silica bodies to investigate their size and distribution, along with the number of carbon inclusions within the silica bodies of perennial grass species Koeleria macrantha. Fluorescence images were analyzed by image software Adobe Photoshop CS5 or ImageJ that remarkably facilitated the quantification of silica bodies in the dry ash. We observed three types of silica bodies or silica body related mineral structures. Silica bodies were detected on both abaxial and adaxial epidermis of K. macrantha leaves, although their sizes, density, and distribution patterns were different. No auto-fluorescence was detected from carbon inclusions.
The combination of fluorescence microscopy and image processing software displayed efficient utilization in the identification and quantification of silica bodies in K. macrantha leaf tissues, which should applicable to biological, ecological and geological studies of grasses including forage, turf grasses and cereal crops.
The epidermal cells of grasses (Poaceae) are arranged in parallel rows with combinations of diverse cell types . Some of these cells are specific for silicon (Si) deposition and are called silica cells. The Silicon accumulated in the silica cells develops into the mineral structures of amorphous hydrated silica (SiO2·nH2O) having various shapes and properties called silica bodies, silica phytoliths, or plant opal [2–4].
Silica bodies are one of the most durable structures in grass tissues that remain as particles in the soil even after all other organic parts of plant have naturally decayed or degraded. These particles in the soil and ash can be very important research tools for systematic botanists [2, 5], environmental biologists , archeologists [4, 7, 8], paleontologists/paleobotanists [9–14], and geologists [15–17].
The amounts of silica in plant tissues suggest that silicon has a very important role in growth and development. For example, in rice (Oryza sativa L.), several fold more Si can be detected in shoots compared with the amounts of nitrogen, phosphorus, or potassium , reaching up to ten percent of its dry mass [19, 20].
Functional analyses of plant silica have shown that silicon is critical for mitigating stressors such as fungal infection [21, 22], herbivory [23, 24], wear [25, 26], and drought [27–30]. Mature silica bodies have been found to deter herbivory and increase the abrasiveness of grass leaf blades [31–33]. In addition, ample silica bodies have been associated with photosynthetic activities [29, 34, 35], although the mechanism for this response remains unclear .
Because we are interested in improving stress tolerance response in turf grasses, we wanted to develop a method to efficiently identify and quantify silica bodies in perennial grasses. Such a method could also be extended to other grass species, such as important forage grasses and cereals. In searching for an easy, economical, and fast method to study the morphology and distributional patterns of silica bodies in turf grasses and other plants, we found a number of available techniques. These include dry ash method, wet oxidation method, scanning electron microscopy (SEM) method, and X-ray image analysis. Among which, dry ash-imaging is one of the most commonly used methods for studying silica bodies in modern plants. To study grass leaves, ash imaging has been a method-of-choice to many researchers; however, this method is extremely labor intensive when analyzing the size, density, and distribution patterns using brightfield light microscopy and researchers have to manually measure a great number of silica bodies in order to perform a statistically meaningful analysis [2, 36]. This method can be accomplished by placing samples in porcelain crucibles and into a muffle furnace, or an oven, for 1–2 h at 500 °C, but some morphological changes might occur to certain, lightly silicified phytoliths when the temperature exceeds 600 °C [2, 4, 36, 37]. The wet oxidation method was developed to examine the isolated silica bodies and is suitable for measuring the abundance of silica bodies in plant tissues, but does not work well for analyzing the distribution patterns of silica bodies [2, 4, 38]. In comparison to the dry ash method, the wet oxidation method results in less damaged silica bodies, especially when the samples are exposed in an environment of 600 °C or higher . Due to the limitation of applying light microscopy to examine surface morphology at extra high magnification, scanning electron microscopy (SEM) can also be used to study silica bodies [40, 41]. The SEM method can be combined with X-ray analysis to provide information on surface structure and composition of silica bodies [42–44]. Here we report a method to study silica bodies using fluorescence microscopy to visualize green autofluorescence in combination with the dry ash-imaging technique. This method was developed using a perennial grass species, Koeleria macrantha, commonly known as junegrass, and has potential to be used in all fields of paleobotany and modern plant sciences on silica body research.
Results and discussion
Combination of fluorescence microscopy and image software provides an opportunity to study plant silica bodies with high efficiency
Silicon-containing structures in plants often display auto- or inducible fluorescence emission that can be examined by fluorescence microscopy . Comparing to brightfield microscopy, the fluorescence emission created ideal conditions for image analysis in high throughput studies, because noisy ash background was eliminated and samples demonstrated clear shape, number, and the distribution pattern of silica bodies (Fig. 1a). Autofluorescence micrographs of an ash-image sample can be easily processed using image software such as Adobe Photoshop and Image J for high-throughput analysis, which include occupancy rate of the leaf epidermis, size and number of silica bodies per unit of leaf surface (Fig. 1a–e). Results can be obtained from a single image or from a random combination of sample images. For example, by analyzing a single image of K. macrantha ‘Ireland’, we observed 11 silica bodies in an area of 4909 µm2 (Fig. 1), which converted to 2240 silica bodies in an area of 1 square millimeter (2240 sb/mm2) abaxial leaf epidermis; the silica bodies occupied 8.1 % of the leaf surface (abaxial). The average size of each silica body was 36 µm2 with a standard deviation of 7.94.
With the same approach we analyzed other accessions using 10 randomly selected images per accession and found that KM-MN and KM-CO exhibited averages of 13,676 sb/mm2 and 13,568 sb/mm2, respectively, which is approximately 6.1 times more silica bodies per square millimeter comparing to the Ireland accession. The sizes of the silica bodies also differed significantly among accessions; for example, we observed the largest silica bodies with an average size of 52.8 μm2 in ‘Barkoel’, whereas the smallest silica bodies with an average size of 26.7 μm2 in ‘Canada’ (Fig. 2) [46, 47].
Three types of silica bodies or silica body related mineral structures were detected by brightfield microscopy (BM) and fluorescence microscopy (FM). Under brightfield microscopy, structures remaining in the ash-imaging process are considered silica bodies [2, 36]. In this study, we compared images from the same ash sample under brightfield and fluorescence microscopy and observed three different types of silica bodies or silica body related structures: Type I silica bodies were developed in the short silica cells and detected under both brightfield and fluorescence microscopy (Fig. 3). Type II silica bodies were also developed in the short silica cells, but only seen under the BM (Fig. 3a), and could not be detected under FM (Fig. 3b). Type III silica bodies were only detected under FM, which were not developed in silica cells (short cells), instead, Type III silica bodies were likely developed in the silicon-enriched long-cells and trichomes (Fig. 3c, d). We frequently observed a number of structures that appeared to look like silica bodies under BM (Fig. 3a), which likely contained less silica and more other chemicals, such as carbon seen in the type I silica bodies as carbon inclusions. However, we were not able to utilize correlative fluorescence microscopy  with X-ray analysis in combination with the ash-imaging method to confirm that the intensity of autofluorescence corresponded to the level of silica deposition. Evidence from previous studies using scanning electron microscopy and X-ray analysis suggested that structures developed in silica cells could be classified into two different types based on the silica profile [49, 50]. The majority of silica cells demonstrated significant silica deposition; only a small portion of the short silica cells (~4 %) exhibited silica body-like structures with barely detected silica signal .
Silica bodies presented in both abaxial and adaxial epidermis of K. macrantha
Most of ash-imaging analysis demonstrated a single layer of silica bodies in grass leaves [2, 36]. By carefully examining the ash-image samples of K. macrantha, we found that silica bodies were presented in both abaxial and adaxial leaf epidermis. We observed two layers of silica bodies in some accessions of junegrass at different focus depths of the dry ash slides. Results from the samples that have only one side of epidermal cell layers, either abaxial or adaxial with removal of the opposite epidermis have confirmed that both the abaxial and adaxial epidermis developed silica bodies (Fig. 4). The density, shape, and size of silica bodies in the adaxial epidermis (Fig. 4a), however, differed from those of the abaxial epidermis (Fig. 4b).
No autofluorescence was detected from carbon inclusions. Carbon inclusions are the residues of organic matter entrapped in silica bodies during the process of silicification [51–53] which can occupy approximately 0.85 % of the total volume of silica bodies [51–53]. Recent studies suggested that the carbon inclusions in plant silica bodies contributed to the enhancement of long-term soil carbon sequestration in agro-ecosystem [54–56]. Nevertheless, we did not observe autofluorescence emitted from the carbon inclusions inside silica bodies (Fig. 5b). The number and distribution pattern of carbon inclusions in silica bodies varied in different K. macrantha accessions and cultivars. There was a statistically significant difference in the number of carbon inclusions within silica bodies among accessions: the lowest average was found in ‘Barkoel’ with 2.08 carbon inclusions per silica body and the highest average was found in KM-CO with 5.56 carbon inclusions per silica body (Fig. 5d). Since the number of carbon inclusions in each accession is stable, and the difference between accessions can be used to distinguish the uniqueness in each accession. In addition, the number and size of carbon inclusions can be used to measure the process of cell silicification and the correlation between silica deposition and carbon accumulation [2, 36].
We propose a method of combining fluorescence microscopy and image processing software for the quantification of silica bodies in Koeleria macrantha leaf tissue, which can be applied to biological, ecological and geological studies of grass species. We observed differences between junegrass accessions for both size and density of silica bodies in leaf epidermis. In addition, we identified differences between accessions for carbon inclusions. This study outlines a means to investigate silica bodies in grass models utilizing a novel high throughput method.
Plant material and sample collection
To examine the structure and properties of silica bodies, mature Koeleria macrantha leaf blades were collected from plants grown in the greenhouse in a 2:1 mixture of Sunshine MVP (Sungro Horticulture) and MVP:Turface (PROFILE Products LLC) substrates with no additional fertilizer. Plant material was derived from: (a) populations from the University of Minnesota turfgrass breeding program derived from material collected in either Colorado (KM-CO), Nebraska (KM-NE) or Minnesota (KM-MN) [55–57]; (b) the cultivar ‘Barkoel’; (c) several accessions from the United States Department of Agriculture National Plant Germplasm System including PI 430287 (Ireland), PI 387927 (Canada), W6 33040 (Russia Federation), PI 207489 (Afghanistan), W6 13043 (China), and PI 302912 (Spain).
The middle section of each leaf blade was cut along the transverse plane and used for ash-imaging sample preparation with removal of the tip and leaf base (Fig. 6b). Each leaf blade was cut into two pieces, one of which was placed on the adaxial and the other on the abaxial side of the leaf onto a glass microslide (VWR Micro Slides were used in this study). On the same microslide, one other leaf from a different plant of the same accession was sampled as a biological replicate (Fig. 6b). A total of 24 plants (thus 12 slides) per accession were analyzed.
Dry ash-imaging sample preparation
A microslide with the leaf samples on was then covered with another glass microslide in an effort to not disturb the placement of the leaves and to add appropriate weight to keep the ash sample intact (Fig. 6b). The slides were then heated on either a Corning Hot Plate Stirrer PC-351 (Fig. 6a) or a Tek-Pro Heat-Stir 36 H2397-1 (not shown) placed in a fume hood. The hot plate temperature was gradually increased every 5 min up to 320 °C. The temperature was approximated using an infrared thermometer Ryobi IR001 (CW0938) read at >608.2 °F/320 °C. The ash process usually took 2–3 h depending on the accession. Grass leaf samples first turned dark brown or black (Fig. 6c), and then gray (Fig. 6d) when the ashing process completed. To end the heating process, we turned off the hot plate and kept the slides on the plate for 1 h or longer to slowly cool them down. (Note: The hot plate stirrer and glass slides are extremely hot while preparing the ash samples; do not move the slides directly to a cooler place while they are still hot, which often result in broken slides.) Slow heating and cooling down prevents the microslides from cracking. For the dry ash samples preparation, we recommend to use a hot plate instead of a coiled electric stove, where the slides often break due to its uneven heating surface.
After the slides were cooled down for at least 1 h, the top microslide was then carefully removed and discarded. A 1 ml plastic transfer pipette (Fisher Scientific, Pittsburgh, PA, USA) was cut at the 0.5 ml measurement to make an ease cedar wood oil application, and a single drop of cedar wood oil (Electron Microscopy Sciences, Hatfield, PA, USA) was applied. A cover slip was then placed on the microscopy slide without disrupting the sample. Cedar wood oil was allowed to diffuse fully under the cover slip with slightly warming the slide on an alcohol burner. The slides were then imaged using an Ernst Leitz Wetzlar 307143.004 microscope (Wetzlar, Germany) and photographed with a SPOT Insight 4 Camera (Diagnostic Instruments, USA). The auto-fluorescence was detected using a Green Fluorescence Protein filter cube (SN: 31001, excitation at 480 nm, beamsplitter at 505 nm, and emission at 535 nm) that was manufactured by Chroma (Chroma Technology Corp, Bellows Falls, VT, USA). The camera interference program used to take the images was SPOT Basic v4.6. Up to 5 sets of brightfield and fluorescent images per object on the slide were taken at 200× and 800× magnification rate. Duplicate fluorescent images were analyzed using Photoshop CS5 (Adobe Systems Incorporated, San Jose, CA, USA) and/or Image J (imagej.net) for silica body occupancy rate of leaf surface (percentage of surface area), the size of silica bodies, and the pattern of silica body distribution (number of silica bodies per unit) on epidermis. Granule-like structures in silica bodies are carbon inclusions, which were counted and recorded using ImageJ, Analyzing Particles tool. The data collected from 20 silica bodies per accession, each with the number of carbon inclusions and their spatial distribution patterns, such as tightly clustered, loosely clustered, randomly distributed, or single granule.
Image and data analysis
First, a random selected fluorescence image was imported into Adobe Photoshop CS5 (Fig. 1a). Second, we used the Photoshop Magic Wand Tool to select the dark image background without silica bodies. Third, the Select-Inverse tool was utilized to select all the silica bodies (Fig. 1b). Fourth, we used the Window-Histogram function to read the pixel counts of the entire image (Fig. 1c) and the counts of selected silica bodies (Fig. 1d). To count the number of silica bodies in an image, we used the Window->Tools->Count Tool in the Photoshop CS5, the number of silica bodies was automatically shown. For those who do not have a licensed Adobe Photoshop CS5 or advanced version, the freeware ImageJ is available (http://imagej.nih.gov/ij/download.html) with tutorials videos on YouTube (youtube.com). Using the Analyzing Particles tool in imageJ we could count the number of silica bodies automatically as well. All statistically significant differences were tested at P < 0.05 level. The results were analyzed by ANOVA using R 3.1.2. . The Tukey multiple comparison test was used to test the significant differences of silica bodies among all accessions studied .
Evert RF. Esau’s plant anatomy. Hoboken: John Wiley & Sons, Inc.; 2006.
Chen C, Feng Z, Chen Y, Chao X. Study on fundamental types and classified system of silica bodies in leaf-blade of Poaceae. J East China Normal Univ (Nat Sci). 1996;1996(1):10.
Piperno DR, Sues HD. Paleontology. Dinosaurs dined on grass. Science. 2005;310(5751):1126–8. doi:10.1126/science.1121020.
Piperno DR. Phytolith analysis: an archaeological and geological perspective. San Diego: Academic Press; 1988.
Prychid CJ, Rudall PJ, Gregory M. Systematis and biology of silica bodies in monocotyledons. Bot Rev. 2004;69(4):377–440. doi:10.1038/ncomms1482.
Sangster AG, Hodson MJ. Silica in higher plants. Ciba Found Symp. 1986;121:90–107.
Rapp GR, Mulholland SC. Society for archaeological sciences (U.S.). Phytolith systematics: emerging issues. In: Advances in archaeological and museum science. New York: Plenum Press; 1992, vol. 1.
Piperno DR, Flannery KV. The earliest archaeological maize (Zea mays L.) from highland Mexico: new accelerator mass spectrometry dates and their implications. Proc Natl Acad Sci USA. 2001;98(4):2101–3. doi:10.1073/pnas.98.4.2101.
Wilding LP. Radiocarbon dating of biogenetic opal. Science. 1967;156(3771):66–7. doi:10.1126/science.156.3771.66.
Pearsall DM. Phytolith analysis of archeological soils: evidence for maize cultivation in formative ecuador. Science. 1978;199(4325):177–8. doi:10.1126/science.199.4325.177.
Thomasson JR, Nelson ME, Zakrzewski RJ. A fossil grass (gramineae: chloridoideae) from the miocene with kranz anatomy. Science. 1986;233(4766):876–8. doi:10.1126/science.233.4766.876.
Ciochon RL, Piperno DR, Thompson RG. Opal phytoliths found on the teeth of the extinct ape Gigantopithecus blacki: implications for paleodietary studies. Proc Natl Acad Sci USA. 1990;87(20):8120–4. doi:10.1073/pnas.87.20.8120.
Prasad V, Stromberg CA, Leache AD, Samant B, Patnaik R, Tang L, et al. Late Cretaceous origin of the rice tribe provides evidence for early diversification in Poaceae. Nat Commun. 2011;2:480. doi:10.1038/ncomms1482.
Lu H, Zhang J, Wu N, Liu KB, Xu D, Li Q. Phytoliths analysis for the discrimination of Foxtail millet (Setaria italica) and Common millet (Panicum miliaceum). PLoS One. 2009;4(2):e4448. doi:10.1371/journal.pone.0004448.
Folger DW, Burckle LH, Heezen BC. Opal phytoliths in a north atlantic dust fall. Science. 1967;155(3767):1243–4. doi:10.1126/science.155.3767.1243.
Zhang J, Lu H, Gu W, Wu N, Zhou K, Hu Y, et al. Early mixed farming of millet and rice 7800 years ago in the Middle Yellow River region, China. PLoS One. 2012;7(12):e52146. doi:10.1371/journal.pone.0052146.
Zhang J, Lu H, Wu N, Yang X, Diao X. Phytolith analysis for differentiating between foxtail millet (Setaria italica) and green foxtail (Setaria viridis). PLoS One. 2011;6(5):e19726. doi:10.1371/journal.pone.0019726.
Savant NK, Snyder GH, Datnoff LE. Silicon management and sustainable rice production. Adv Agron. 1997;58:151–99. doi:10.1015/S0065-2113(08)60255-2.
Epstein E. The anomaly of silicon in plant biology. Proc Natl Acad Sci USA. 1994;91(1):11–7.
Ma JF, Tamai K, Ichii M, Wu GF. A rice mutant defective in Si uptake. Plant Physiol. 2002;130(4):2111–7. doi:10.1104/pp.010348.
Chen X, Li H, Chan WF, Wu C, Wu F, Wu S, et al. Arsenite transporters expression in rice (Oryza sativa L.) associated with arbuscular mycorrhizal fungi (AMF) colonization under different levels of arsenite stress. Chemosphere. 2012;89(10):1248–54. doi:10.1016/j.chemosphere.2012.07.054.
Park JJ, Kim KW, Park TJ, Park EW, Kim Y. Solid-state NMR spectroscopy of silicon-treated rice with enhanced host resistance against blast. Anal Sci. 2006;22(4):645–8. doi:10.2116/analsci/22.645.
Hunt JW, Dean AP, Webster RE, Johnson GN, Ennos AR. A novel mechanism by which silica defends grasses against herbivory. Ann Bot. 2008;102(4):653–6. doi:10.1093/aob/mcn130.
Garbuzov M, Reidinger S, Hartley SE. Interactive effects of plant-available soil silicon and herbivory on competition between two grass species. Ann Bot. 2011;108(7):1355–63. doi:10.1093/aob/mcr230.
Peters CR. Electron-optical microscopic study of incipient dental microdamage from experimental seed and bone crushing. Am J Phys Anthropol. 1982;57(3):283–301. doi:10.1002/ajpa.1330570306.
Massey FP, Hartley SE. Physical defences wear you down: progressive and irreversible impacts of silica on insect herbivores. J Anim Ecol. 2009;78(1):281–91. doi:10.1111/j.1365-2656.2008.01472.x.
Schaller J, Brackhage C, Paasch S, Brunner E, Baucker E, Dudel EG. Silica uptake from nanoparticles and silica condensation state in different tissues of Phragmites australis. Sci Total Environ. 2013;442:6–9. doi:10.1016/j.scitotenv.2012.10.016.
Ma JF. Silicon transporters in higher plants. Adv Exp Med Biol. 2010;679:99–109.
Chen W, Yao X, Cai K, Chen J. Silicon alleviates drought stress of rice plants by improving plant water status, photosynthesis and mineral nutrient absorption. Biol Trace Elem Res. 2011;142(1):67–76. doi:10.1007/s12011-010-8742-x.
Richmond KE, Sussman M. Got silicon? The non-essential beneficial plant nutrient. Curr Opin Plant Biol. 2003;6(3):268–72. doi:10.1016/S1369-5266(03)00041-4.
Vicari M, Bazely DR. Do grasses fight back? The case for antiherbivore defences. Trends Ecol Evol. 1993;8:137–41. doi:10.1016/0169-5347(93)90026-L.
Massey FP, Ennos AR, Hartley SE. Herbivore specific induction of silica-based plant defences. Oecologia. 2007;152(4):677–83. doi:10.1007/s00442-007-0703-5.
Massey FP, Ennos AR, Hartley SE. Grasses and the resource availability hypothesis: the importance of silica-based defences. J Ecol. 2007;95:414–24. doi:10.1111/J.1365-2745.2007.01223.x.
Detmann KC, Araujo WL, Martins SC, Sanglard LM, Reis JV, Detmann E, et al. Silicon nutrition increases grain yield, which, in turn, exerts a feed-forward stimulation of photosynthetic rates via enhanced mesophyll conductance and alters primary metabolism in rice. New Phytol. 2012;196(3):752–62. doi:10.1111/j.1469-8137.2012.04299.x.
Agarie S, Agata W, Uchida H, Kubota F, Kaufman PB. Function of silica bodies in the epidermal system of rice (Oryza sativa L.): testing the window hypothesis. J Exp Bot. 1996;47(298):655–60. doi:10.1093/jxb/47.5.655.
Chen C, Ni B, Xu S, Ma Q. On morphology and morphogenesis of silica bodies in Bambusa multiplex var. nana. Bull Bot Res. 2003;23(4):5.
Wu Y, Wang C, Hill DV. The transformation of phytolith morphology as the result of their exposure to high temperature. Microsc Res Tech. 2012;75(7):852–5. doi:10.1002/jemt.22004.
Gallego L, Distel RA. Phytolith assemblages in grasses native to central Argentina. Ann Bot. 2004;94(6):865–74. doi:10.1093/aob/mch214.
Sun X, Wu Y, Wang C, Hill DV. Comparing dry ashing and wet oxidation methods. The case of the rice husk (Oryza sativa L.). Microsc Res Tech. 2012;75(9):1272–6. doi:10.1002/jemt.22060.
Sangster AG. Studies of opaline silica deposits in the leaf of Sieglingia decumbens L. ‘Bernh’., using the scanning electron microscope. Ann Bot. 1968;32(2):237–40.
Theunissen JD. A method for isolating and preparing silica bodies in grasses for scanning electron microscopy. Biotech Histochem. 1994;69(5):291–4.
Kaufmian PB, Bigelow WC, Petering LB, Drogosz FB. Silica in developing epidermal cells of Avena internodes: electron microprobe analysis. Science. 1969;166(3908):1015–7. doi:10.1126/science.166.3908.1015.
Takeoka Y, Matsumura O, Kaufman PB. Studies on silicification of epidermal tissues of grasses as investigated by soft X-ray image analysis. Jpn J Crop Sci. 1983;52(4):544–50. doi:10.1626/jcs.52.544.
Lichtenberger O, Neumann D. Analytical electron microscopy as a powerful tool in plant cell biology: examples using electron energy loss spectroscopy and X-ray microanalysis. Eur J Cell Biol. 1997;73(4):378–86.
Soukup M, Martinka M, Cigan M, Ravaszova F, Lux A. New method for visualization of silica phytoliths in Sorghum bicolor roots by fluorescence microscopy revealed silicate concentration-dependent phytolith formation. Planta. 2014;240(6):1365–72. doi:10.1007/s00425-014-2179-y.
Dixon JM. Koeleria macrantha (Ledeb.) Schultes (K. alpigena Domin, K. cristata (L.) Pers. pro parte, K. gracilis Pers., K. albescens auct. non DC.). J Ecol. 2001;88(4):709–26. doi:10.1007/s00216-004-2521-5.
Piperno DR, Pearsall DM. The silica bodies of tropical American grasses: morphology, taxonomy, and implications for grass systematics and fossil phytolith identification. In: Smithsonian contributions to botany. Washington, D.C.: Smithsonian Institution Press; 1998, vol. 85.
Jahn K, Barton D, Braet F. Correlative fluorescence- and scanning, transmission electron microscopy for biomolecular investigation. In: Mendez-Vilas A, Diaz JD, editors. Modern research and educational topics in microscopy. Extremadura: Formatiex; 2007. p. 203–11.
Kaufman PB, Yakeoka Y, Bigelow WC. Scanning electron microscopy and X-ray microanalysis of silica in the leaf sheath pulvinus and internodal intercalary meristem of rice. Jpn J Crop Sci. 1979;48(Sup.1):187–8.
Kaufman PB, Dayanandan P, Takeoka Y, Bigelow WC, Jones JD, Iler R. Silica in shoots of higher plants. In: Simpson TL, Volcani BE, editors. Silicon and siliceous structures in biological systems. New York: Springer-Verlag New York, Inc; 1978. p. 409–49.
Gallagher KL, Alfonso-Garcia A, Sanchez J, Potma EO, Santos GM. Plant growth conditions alter phytolith carbon. Front Plant Sci. 2015;6:753. doi:10.3389/fpls.2015.00753.
Rowlett RM, Pearsall DM. Archaeological age determinations derived from opal phytoliths by thermoluminescence. In: Pearsall DM, Piperno DR, editors. Current research in phytolith analysis: applications in archaeology and paleoecology. Philadelphia: University of Pennsylvania; 1993. p. 25–30.
Krull ES, Skjemstad JO, Graetz D, Grice K, Dunning W, Cook G, et al. 13C-depleted charcoal from C4 grasses and the role of occluded carbon in phytoliths. Org Geochem. 2003;34:1337–52. doi:10.1016/S0146-6380(03)00100-1.
Jones RL, Beavers AH. Division S-5—soil genesis, morphology, and classification—aspects of catenary and depth distribution of opal phytoliths in Illinois soils. Soil Sci Soc Proc. 1964;28:413–6.
Currie HA, Perry CC. Silica in plants: biological, biochemical and chemical studies. Ann Bot. 2007;100(7):1383–9. doi:10.1093/aob/mcm247.
Rajendiran S, Vassanda Coumar M, Kundu Ajay S, Dotaniya ML, Subba Rao A. Role of phytolith occluded carbon of crop plants for enhancing soil carbon sequestration in agro-ecosystems. Curr Sci. 2012;103(8):911–20.
Yang J, Wu J, Jiang P, Xu Q, Zhao P, He S. A study of Phytolith-occluded Carbon Stock in Monopodial Bamboo in China. Sci Rep. 2015;5:13292. doi:10.1038/srep13292.
R Development Core team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2011.
Shaffer JP. Multiple hypothesis testing. Annu Rev Phychol. 1995;46:561–84. doi:10.1146/annurev.ps.46.020195.003021.
CD, EW, and CC designed the experiments and wrote the manuscript; CD and JO performed the experiments. All authors read and approved the final manuscript.
We would like to thank Dr. Stefanie Dukowic-Schulze for technical assistance, and Andrew Hollman for plant care. We deeply appreciated two reviewers’ comments, which helped us to improve the protocol significantly. This project was sponsored by funds from the United States Golf Association (2007-16-357) and the Minnesota Agricultural Experimental Station (MIN-21-041, MIN-21-031). The authors are grateful to the University of Minnesota Libraries for funds to support this open access publication.
The authors declare that they have no competing interests.