- Research article
- Open Access
Roles of metabolic regulation in developing Quercus variabilis acorns at contrasting geologically-derived phosphorus sites in subtropical China
BMC Plant Biology volume 20, Article number: 389 (2020)
Phosphorus (P) -rich soils develop in phosphorite residing areas while P-deficient soils are ubiquitous in subtropical regions. Little has been reported that how metabolites participate in the seed development and the processes involved in their coping with contrasting-nutrient environments.
Here we quantified the metabolites of Quercus variabilis acorns in the early (July), middle (August), late (September) development stages, and determined element (C, H, O, N, P, K, Ca, Mg, S, Fe, Al, Mn, Na, Zn, and Cu) concentrations of acorns in the late stage, at geologically-derived contrasting-P sites in subtropical China. The primary metabolic pathways included sugar metabolism, the TCA cycle, and amino acid metabolism. Most metabolites (especially C- and N-containing metabolites) increased and then decreased from July to September. Acorns between the two sites were significantly discriminated at the three stages, respectively, by metabolites (predominantly sugars and organic acids). Concentrations of P, orthophosphoric acid and most sugars were higher; erythrose was lower in late-stage acorns at P-rich sites than those at P-deficient sites. No significant differences existed in the size and dry mass of individual acorns between oak populations at the two sites.
Oak acorns at the two sites formed distinct metabolic phenotypes related to their distinct geologically-derived soil conditions, and the late-stage acorns tended to increase P-use-efficiency in the material synthesis process at P-deficient sites, relative to those at P-rich sites.
Subtropical soils are often characterized by phosphorus (P), calcium (Ca), and magnesium (Mg) deficiencies, with iron (Fe) and aluminum (Al) enriched [1, 2]. However, there are often some P-rich ores mixed in the P-deficient sites in some subtropical areas of China, which leads to significant changes in P and other elements across the region [3, 4]. These contrasting-P sites generally give rise to plants with different stoichiometry characteristics [3, 5,6,7,8], which can affect the metabolism and formation of compounds within cells and organisms [9, 10]. The metabolome represents all of the small molecules in organisms at a given moment [11, 12]. In ecological and physiological studies, metabolomics is primarily employed to explore the physiological status of organisms in response to variable environments [13,14,15]. To date, non-target metabolomics has been applied in many field experiments to elucidate the effects of environmental changes on the composition of metabolites [16,17,18].
Since the concept of ionome is first presented in 2003 , it has been used in the analysis of the responses to environmental factors in different plants , such as Arabidopsis (Arabidopsis thaliana) , tea (Camellia sinensis) , rice (Oryza sativa) . With response to the variation of the genetic, developmental and environmental factors, elements are always bound to biological molecules [24,25,26], which are involved in lots of metabolisms, including carbon and nitrogen metabolism [27, 28]. Hence, the integration of metabolomics with element variations has enhanced the understanding of phenotypic responses involving both physiological and molecular mechanisms to different spatiotemporal scales [16, 29].
Serving as a sink tissue, seeds can mobilize stored nutrients to support the nascent processes of seed germination and seedling establishment in the plant life cycle [30,31,32]. Based on metabolomics, under a certain environment condition, researches have been conducted on the seed development of different plants, such as maize (Zea mays) [33, 34], barley (Hordeum vulgare) , soybean (Glycine max) [35,36,37] and lotus (Nelumbo nucifera) . However, the adaption of seed metabolomes during development to variable soil nutrients in the field experiments has been rarely studied, which would deepen our knowledge of the adaption mechanism of the in situ plant to the environment.
Quercus variabilis is an important deciduous broadleaved oak for ecology, economy, and culture in Eastern Asia. In previous studies, there were distinct stoichiometry characteristics in both the leaves and seeds of Q. variabilis populations growing at P-rich and P-deficient sites [3, 39]. Moreover, strong correlations between metabolites and stoichiometry characteristics of leaves were exhibited in Q. variabilis populations at the two sites . Therefore, for this study, using developing Q. variabilis seeds as samples, the questions we aimed to address were: 1) How do the metabolite profiles of developing plant seeds differ between P-rich and P-deficient sites? 2) Which specific metabolites play key roles in the metabolic regulation of seeds, in coping with variable nutrients at P-rich and P-deficient sites? This study aims to fill the gaps in our understanding as relates to how in situ plants metabolically adapt to sites with contrasting nutrient availability.
Morphological characteristics of developing acorns at P-rich and P-deficient sites
The morphological characteristics of developing acorns at two contrasting-P sites were shown in Fig. 1: Acorn length significantly increased from July to October (p < 0.05) at both the P-rich and P-deficient sites (Fig. 1a); the width and dry mass of individual acorns significantly increased during developing, until reached the maximum values in September (Fig. 1b, c). Further, there were no significant differences in the acorn length and width, or the dry mass of single acorns between the P-rich and P-deficient sites at the four developmental stages (p > 0.05) (Fig. 1).
Metabolomic changes in developing acorns
A total of 100 metabolites were annotated and quantified at the three development stages, as shown in Table S1 (Additional file 4), which could be categorized into six groups based on the molecular structure: amino acids, sugars, organic acids, alcohols, amines, and esters. Additionally, a metabolic map of the process was developed based on the results of pathway analysis, which encompassed all of the metabolites of the three developmental stages (Fig. 2). The identified metabolites of the different developmental stages were different (Fig. 2). The central metabolic pathways of the acorns, including sugar metabolism, the TCA cycle, and amino acid metabolism, were observed in July, August, and September (Fig. 2; Additional file 5: Table S2). The target pathways were those with a pathway impact (PI) of > 0.1 (Additional file 5: Table S2). As shown in Table S2 (Additional file 5), β-alanine metabolism, pantothenate, and CoA biosynthesis were found in the acorns in July only.
To better compare the metabolites of the acorns, we compared the concentrations of 47 metabolites (predominately sugars and organic acids), which were identified in acorns of all the three developmental stages (Additional file 2: Figure S2). Almost all of these 47 metabolites from the two sites of contrasting P availability were highest in August, except for glycine, d-mannopyranose, allose, glutaric acid, and ethanolamine (Additional file 2: Figure S2). The orthophosphoric acid, d-fructose, sucrose, and myo-inositol concentrations revealed obvious changes within the developing acorn at P-rich and P-deficient sites (Fig. 3; Additional file 4: Table S1): the concentrations of orthophosphoric acid and sucrose in acorns differed significantly among the three developmental stages at both the P-rich and P-deficient sites (p < 0.05); the concentrations of d-fructose and myo-inositol differed between the acorns in July and August, as well as acorns in August and September at the two sites (p < 0.05). Concretely, the concentrations of orthophosphoric acid, myo-inositol and d-fructose in the acorns in August were above 20 fold higher than those both in July and September at two sites; sucrose concentration in August were more than 7 fold higher than in July, and more than 25 fold higher than in September at the two sites (Fig. 3).
Differential metabolites and discriminating elements in acorns between P-deficient and P-rich sites
Variance analysis of the metabolite concentration in developing Q. variabilis acorns at two contrasting-P sites revealed that the predominant metabolites of significant difference in the acorns of the two sites were glyceric acid, l-lactic acid, d-lyxose, and glycine in July, sugars (e.g., erythrose, allose, sucrose) and organic acids (e.g., glycine, linoleic acid, tartaric acid) in August, and sugars (e.g., erythrose, cellobiose, floridoside) in September (Fig. 5).
According to the significant OPLS-DA model with R2 value of > 0.7 and Q2 value of > 0.5 (Additional file 3: Fig. S3A), the acorns of the two contrasting-P sites could be distinguished by metabolites (Fig. 4a) in July, with major contributions from four metabolites (glyceric acid, l-lactic acid, d-lyxose, and glycine). The concentrations of these metabolites were higher at the P-rich sites in contrast to the P-deficient sites (Fig. 5a; Additional file 4: Table S1). Similarly, based on significant OPLS-DA model with R2 value of > 0.7 and Q2 value of > 0.5 (Additional file 3: Figures. S3B, S3C), the acorns at the two contrasting-P sites were significantly discriminated by metabolites, with the major contributions from 26 metabolites in August (Fig. 4b, Fig. 5b), and 16 metabolites in September (Fig. 4c, Fig. 5c). In August, among the 26 metabolites, except for glycine and erythrose, other metabolites exhibited significantly higher concentrations at the P-rich sites in contrast to the P-deficient sites (p < 0.05), including sucrose and orthophosphoric acid (Fig. 5b; Additional file 4: Table S1). In particular, erythose concentration was higher at P-deficient sites, relative to the P-rich sites (Fig. 5b; Additional file 4: Table S1). In September, the 16 differential metabolites at the P-deficient sites were more abundant in contrast to those at the P-rich sites, including orthophosphoric acid. In particular, the acorn sugar (floridoside, glyceryl-glycoside, and β-d-glucopyranose) and organic acid (ribonic acid) concentrations at the P-deficient sites were two-fold higher compared with the P-rich sites (Fig. 5c; Additional file 4: Table S1).
Considering the physiological maturity of Q. variabilis acorns at the two sites was in September, when the dry mass of acorns reached stability during developing (Fig. 1), we also focused on metabolites-elements relationships of acorns in September. Based on OPLS-DA (Additional file 3: Fig. S3), the acorns from the P-rich and P-deficient sites were significantly discriminated via the ionomic data in September (Fig. 4c), with the major contributions from five elements (N, P, S, Mn, and Cu) (Fig. 6). Besides, concentrations of acorn N, P, S and Cu was higher, and acorn Mn was lower at P-rich sites than those at P-deficient sites (Fig. 6).
The relationships among these discriminating elements were complex (Additional file 7: Table S4): there were significantly positive correlations between N and P, between N and S, between P and S, and between Cu and S, and significantly negative correlations between N and Mn, and between Cu and Mn. Moreover, in September, acorn N and Mn showed complex correlations with acorn differential metabolites (Additional file 8: Table S5): acorn N was positively correlated with glycine (amino acid) and orthophosphoric acid, but negatively related to some sugars (e.g., l-sorbofuranose, β-d-glucopyranose, floridoside and glyceryl-glycoside); Mn was positively correlated with some sugars (e.g., β-d-glucopyranose, floridoside and glyceryl-glycoside), but negatively correlated with orthophosphoric acid (p < 0.05) (Fig. 6). Besides, P and S were positively related to glycine, negatively related to some sugars (e.g., floridoside and glyceryl-glycoside), and Cu was negatively related to some sugars (e.g., Lactulose, β-d-glucopyranose, floridoside and glyceryl-glycoside) in September (Additional file 8: Table S5; Fig. 6).
To the best of our knowledge, this study was the first to demonstrate how in situ plants generate seeds by adapting to variable nutrient environments, in terms of the metabolome of developing seeds in subtropical regions, where P was limiting for plant growth and development. Based on our results, the discussions focused on differences in the metabolic profiles of acorns, and the identification of metabolites that played key roles in metabolic regulation at the contrasting P sites.
Sequential dynamic in C- and N-containing metabolites during acorn development
Our results clearly revealed that, for the oak populations at the two sites of contrasting P availability, from July to September, concentrations of most sugars, amino acids, and organic acids in acorns initially increased and then decreased (Additional file 4: Table S1; Fig. 2, Additional file 2: Fig. S2), with PI values of the same metabolomic pathways in the three developmental stages different (Additional file 5: Table S2). These results suggested that, during developing, metabolic states of acorns changed based on the sequentially dynamic of C- and N- containing metabolites. This was consistent with the results of Wang et al. , who found that the metabolites and proteins involved in the development of Lotus (Nelumbo nucifera) seeds were sequentially dynamic. On one hand, this might be caused by the vital function of these metabolites in seed development. From July to August, acorns synthesized many C- and N- containing metabolites (e.g., amino acids, sugars, and organic acids) to synthesize the storage substances in September, when acorns reached the physiological maturation with the dry mass being relatively stable (Fig. 1) . As reported by others, key metabolites (e.g., amino acids, sugars, and organic acids) were significantly altered in developing maize , lotus , and wheat , even in developing fruits, such as, medlar (Mespilus germanica) , navel oranges (Citrus sinensis) , and ponkan (Citrus reticulata) .
On the other hand, the sequentially dynamic of metabolites in developing acorns at the contrasting-P sites might be related to the ordered regulation of genes in metabolite synthesis during plant development. The significant reduction of C- and N- containing metabolites was regulated by gene expression . Although the most active genes in seeds were shared throughout development, each developmental stage utilized a set of genes that was up-regulated compared with other stages, or special at the GeneChip level . Similarly, the transcriptome research of developing acorns by Miguel et al.  revealed that, carbohydrate metabolism, including sugar metabolism, the TCA cycle, and amino acid metabolism, was most represented during acorn development, but every stage had special up-regulated metabolisms, such as carbohydrates and energy metabolism for the middle stage. Similarly, even though some special metabolisms appeared in July, sugar metabolism, the TCA cycle, and amino acid metabolism comprised the core metabolic pathways in developing acorns (Fig. 2; Additional file 5: Table S2), with the size and dry mass increasing significantly from July to September at both P-rich and P-deficient sites (p < 0.05) (Fig. 1).
Besides, as a storage organ, acorns primarily accumulated starch (46.3–68.6%), followed by fats and proteins . Therefore, during acorn development at the two contrasting-P sites, noticeable changes occurred in the concentrations of sugars, sugar alcohols, and amino acids, and the changes were related to their functions in the synthesis of starch and proteins (Fig. 2). Firstly, as a form of energy for plants, starch steadily accumulated during acorn development , resulting in the remarkable increase of acorn dry mass (Fig. 1). During this process, small sugar molecules (particularly sucrose and d-fructose) were converted to large sugar molecules (e.g., raffinose and starch) [38, 41], coupled with the consumption of sugar alcohols (e.g., myo-inositol and galactitol) (Fig. 2) , or acted as vital regulators [51, 52]. The accumulated large sugar molecules would serve as energy reserves (Fig. 2) , which could provide energy materials in the establishment of seedlings. Hence, at the two sites, during maturation process (from August to September), concentrations of sucrose and d-fructose, as well as myo-inositol and galactitol decreased significantly (Figs. 3, Additional file 2: Figure S2; Additional file 4: Table S1). Secondly, amino acids, which were generally synthesized in leaves, could be transported into seeds for further protein synthesis , and acted as intermediates for glycolysis and the TCA cycle in energy metabolism (Fig. 2) [55,56,57]. Hence, the decrease of most amino acids (particularly l-glutamic acid) during maturation (Additional file 2: Figure S2; Additional file 4: Table S1) might have been caused by the incorporation of amino acids in the synthesis of storage proteins.
Differential C- and N-containing metabolites in oak acorn populations at P-rich and P-deficient sites
Our results clearly showed that metabolites, predominantly sugars and organic acids, differed significantly in acorns of oak populations at the P-rich and P-deficient sites (Fig. 5). This suggested that these types of sugars and organic acids were involved in the responses of in situ plant seeds to the variable supplies of soil nutrients due to geologically-derived environments. This had been proved by Ji et al. , who found that sugars and organic acids of Q. variabilis leaves played predominant roles in clearly discriminating Q. variabilis trees at the contrasting geologic-P sites. The vital functions of sugars and organic acids were also showed in other studies on the response of plant leaves and seeds under different nutritional conditions. Yan et al.  revealed that the added P could increase concentrations of sugars in oilseed flax seeds. P deficiencies  and Zn stress  predominantly altered concentrations of sugars and organic acids in tea (Camellia sinensis) leaves.
As the reproductive organ of Q. variabilis trees, acorns were characterized by the enrichment of starch during the development process . Carbohydrate metabolism, including glycolysis, gluconeogenesis, amino sugar and nucleotide sugar metabolism, served as the essential biochemical process (Fig. 2) . As two vital forms of metabolites, sugars and organic acids hold a larger percentage in plants than other metabolites, and sugars (e.g., fructose, sucrose, glucose) and organic acids (e.g., shikimate, glycine, glutamate) were the primary intermediates that participated in the carbon and nitrogen metabolism of plant seeds (Fig. 2) . Similarly, the important functions of sugars and organic acids were observed in developing fruits, such as ponkan (Citrus reticulata) , navel oranges (Citrus sinensis) , and medlar (Mespilus germanica) .
Increased efficiency of P use in the late stage acorns at P-deficient sites over P-rich sites
In September, the dry mass of Q. variabilis acorns reached stable (Fig. 1). We found that there were distinct metabolic shifts in acorns of this stage at such contrasting geologic-P sites (Fig. 6). The distinct metabolic shift of plants under different nutrient supply situation has been reported in other studies. Christian and Oliver  showed that Chlamydomonas reinhardtii developed highly distinctive metabolite profiles under N, P, S, or Fe deficiency conditions. Gargallogarriga et al.  revealed that increased P availability could lead to shifts in the metabolome through higher investments in the protection mechanisms of plants.
Concretely, our results indicated that the utilization of P was different to guarantee the same dry mass in acorns at P-rich and P-deficient sites with different soil nutrient supply (Fig. 6). Erythrose could be synthesized by erythrose-4-phosphate, which was an important intermediate product of the pentose phosphate pathway [63, 64], and the pentose phosphate pathway consumed less inorganic P than glycolysis . In our study, the concentrations of most sugars, most organic acids, orthophosphoric acid, and P were lower, and erythrose was higher in the late-stage acorns at P-deficient sites than those at P-rich sites. But the dry mass and C concentration showed no significant differences in late-stage acorns at the two sites (Fig. 6; Additional file 6: Table S3). Unlikely, seeds of some Proteaceae species were considered to accumulate comparatively higher P concentration in seeds to adapt to low P environments , and produce small seeds in such P-impoverished habitats . P addition could enhance the dry mass of maize (Zea mays) . Besides, shifts in P usage were associated with variable P supplies to plants [69,70,71]. Hence, late-stage acorns at P-deficient sites might boost the pentose phosphate pathway to decrease P quotas through bypassing P consumption in glycolysis reactions, which increased the efficiency of P use in acorns at P-deficient sites (Fig. 6).
Further, similar to previous reports concerning the element concentrations of other organs [25, 72], the protein/nucleic-related (e.g., N, P, and S) and enzyme-related (e.g., Mn, Cu) elements of late-stage acorns were more highly impacted, in contrast to other types of elements at P-rich and P-deficient sites in the present study (Fig. 6). However, unlike our results, supply of P could enhance the concentrations of macro-elements (e.g., N, P, and S) in crops, such as urdbean (Vigna mungo) , wheat (Triticum aestivum), maize (Zea mays), and faba bean (Vicia faba) , but decrease the concentrations of micro-elements (e.g., Mn and Cu) in rice (Oryza sativa) . Considering the strongly close relationships between these elements (e.g., N, P, S, Mn, and Cu) and metabolites in plants (Fig. 6) , and the vital effects of these elements on the synthesis of storage substances in plant seeds and the yield of plant seeds (Fig. 6) [75, 77, 78], it also could prove the distinct metabolic shifts in the late-stage acorns at the two sites.
Metabolic regulation is one of the essential mechanisms for plants to generate formal seeds under the in situ conditions, by adapting to stresses due to nutrient deficiencies in soils, particularly in geologically-derived nutrient limited areas. In the present study, our results characterized the roles of metabolic regulation in the development of woody plant seeds at contrasting P sites in subtropical regions where soils were featured with deficient P, Ca and Mg. First, sugar metabolism, the TCA cycle, and amino acid metabolism were vital for the development of oak acorns due to their involvement in material synthesis. Moreover, the metabolites involved were sequentially dynamic, from July to September, with the most significant changes observed in the concentrations of orthophosphoric acid, d-fructose, sucrose, and myo-inositol. This demonstrated the characteristics of energy and nutrient accumulation in oak seeds. Second, there were significant differences in the acorns of oaks between the P-rich and P-deficient sites in terms of metabolites (predominantly sugars and organic acids) and elements, while no differences were observed in the size and the dry weight of individual acorns between P-rich and P-deficient sites. This verified the inconsistent between morphologies and chemical compositions of acorns in oak populations of the two site types. Third, in the late-stage (September) acorns, the concentrations of P, orthophosphoric acid and most sugars (particularly D-fructose and sucrose) were significantly higher, and erythrose was lower, at the P-rich sites in contrast to those at the P-deficient sites. This suggested that the late-stage acorns tended to increase the efficiency of P use in the process of material synthesis at the P-deficient sites relative to these at the P-rich sites. Our findings deepen our understanding of how in situ plants adapt to variable nutrient environments through metabolic mechanisms in subtropical regions.
This study was based on a field experiment. The P-rich sites develop on phosphate rocks with high P content, locating in Kunming City (latitude 24°58′54.38″N, longitude 102°26′47.58″E, altitude 1869 m a.s.l), whereas the P-deficient sites develop on non-phosphate rocks with low P content, locating at Mouding, in Chuxiong City (latitude 25°14′48.25″N, longitude 101°32′35.62″E, altitude 1846 m a.s.l), of Yunnan Province, China. The stands at the P-rich (Kunming) and P-deficient (Mouding) sites are natural forests, which have been left anthropogenically untouched through history . So, the trees we selected were the same age. Further, the chemical compositions of these phosphate and non-phosphate rocks have been revealed by Ji . The level of P2O5 in the phosphate rocks is significantly higher over the non-phosphate rocks. The total P concentration are 2.23 mg g− 1 at P-rich sites and 0.3 mg g− 1 at P-deficient sites, while the available P concentration are 0.045 mg g− 1 at P-rich sites and 0.007 mg/g at P-deficient sites .
The climates of the two sites are similar, with the average annual temperatures of 16.0 °C at the two sites and the mean annual precipitation of 909.40 mm at Mouding, and 978.00 mm at Kunming in 2017. Figure S1 (Additional file 1) depicts the average monthly temperatures and precipitation of Mouding and Kunming, from 1981 to 2010, as well as the precipitation of the two sites, which occurred primarily during July and August (Additional file 1: Figure S1).
The chemical compositions of the soil from P-rich sites at Kunming and P-deficient sites at Mouding have been revealed by Wen et al.  and Ji et al. . The pH of the soils at the two sites was ~ 5.22, with the soil C, N, P, K, Mg, and Mn concentration being much higher at the P-rich sites, in contrast to the P-deficient sites [16, 80]. Contingent on the location and extent of the Q. variabilis stands, six P-rich sites (20 m × 20 m) and six P-deficient sites (20 m × 20 m) were established in July. In August and September, eight P-rich sites (20 m × 20 m) and 12 P-deficient sites (20 m × 20 m) were established, respectively, where each site contained three plots.
Wild acorns were sampled from Quercus variabilis populations at P-rich (Kunming) and P-deficient (Mouding) sites in central Yunnan, China, where Q. variabilis trees were typically distributed. No permission was required to collect these acorn samples, and the study had no impact on the biological diversity of central Yunnan. Prof. Chunjiang Liu undertook the formal identification of samples at the two sites. In the P-rich site stand, the height of the Q. variabilis trees ranged from 17.67 m to 31.67 m, and the DBH (diameter at breast height) spanned 21.74 cm to 48.41 cm. In the P-deficient site stand, the height of Q. variabilis trees ranged from 22.50 m to 25.00 m, with the DBH spanning 16.03 cm to 28.66 cm. Four stages of acorn development, early stage (July), middle stage (August), late stage (September) and the last stage (October), were established, based on variations of developmental characteristics and the content of major substance (starch). Acorns were completely covered by the cupule in July, which was visible and in green in August, but became brown in September, and fell in October. Besides, in July and August, the accumulation rate of starch was much faster than other developmental stages, and in September, the starch content reached the maximum [49, 81]. So, morphological characteristics of Q. variabilis acorns were analyzed at the four developmental stages, but the metabolite profiling in acorns of July, August, and September were analyzed. In 2017, six acorns were randomly collected from three trees at each oak plot in the middle of July, August, September, and October, respectively. Parts of cotyledon from nine non-parasitized and peeled acorn samples from a given site were formed a composite sample in July. Parts of cotyledon from three non-parasitized and peeled acorn samples from a given plot were formed a composite sample in August and September. Voucher specimens were deposited in a public herbarium in School of Agriculture and Biology of Shanghai Jiao Tong University. To conduct GC-MS analysis, the acorns in July, August and September were immediately frozen in dry ice and stored at − 80 °C. The length, width, and dry weight of individual acorns at the four developmental stages were determined.
Extraction, derivation, and analysis of acorn metabolites
Metabolite profiling analysis was performed using gas chromatography-mass spectrometry (7890A-5975C, Agilent, USA). The metabolite extraction technique was modified from Wu et al.  and Du et al. . After composite acorn samples were defrosted, a 100 mg quantity of the fresh composite sample was introduced into a 2.0 mL centrifuge tube. Subsequently, 0.8 mL methanol-chloroform (3:1, v/v) and 30.0 μL of ribitol (2.0 mg mL− 1 stock in water) was introduced into each tube as an internal quantitative standard. The mixture was ground at 60.0 HZ for 80 s and centrifuged at 1200 r min− 1. for 10 min. 400.0 μL of the polar phase sample (aqueous and organic) was collected independently into 1.5 mL HPLC glass vials and dried in a benchtop centrifugal concentrator (Labconco Corporation, Kansas City, MI) for 3.0 h. Once the upper polar phase was thoroughly dried, methoximation, incubating the dried fraction at 37 °C for 1.5 h with 80.0 μL 15 mg mL− 1 methoxyamine hydrochloride, were carried out at first, and then trimethylsilylation, incubating the dried fraction at 70 °C for 1.0 h with 80.0 μL TMCS (BSTFA: TMCS = 99:1), was conducted.
The GC-MS analysis was carried out according to Du et al. . The derivatization of samples was carried using a PerkinElmer gas chromatograph and TurboMass-Autosystem XL mass spectrometer (PerkinElmer lnc., Waltham, MS). A 1 μL aliquot of each sample was injected into a DB-5MS capillary column (30 mm × 0.25 mm × 250 mm) (Agilent JW Scientific, Folsom, CA). Following a solvent delay of 5 min, the GC oven temperature was adjusted to 60 °C; after injection for 1 min, the temperature of oven was raised from 60 °C to 300 °C at 5 °C min− 1 for 27 min. The injector and ion source temperatures were adjusted to 280 °C and 230 °C, respectively. Helium was applied as the carrier gas at a constant rate of 1.0 mL min− 1. Measurements were achieved with an electron impact at 70 e and at full scan mode, with a mass scan range of 33–600 m z− 1.
The metabolite levels were determined by the mass spectral library (2011) in MSD ChemStation (version E.02.02.1431; Agilent Inc., CA, USA), the National Institute of Standards and Technology (NIST), and ChromaTOF (version 22.214.171.124; Leco Corporation, MI, USA). The raw files were converted to NetCDF format, and then sequentially processed by ChromaTOF (version 126.96.36.199; Leco Corporation, MI, USA) for correcting the baseline, searching peaks, deconvoluting spectrum, and aligning compounds of the different samples. A number of artificial peaks were removed from the dataset based on the mass spectra and retention index (RI) comparison to the NIST library. The metabolites were expressed in the peak area, then normalized to the area of the internal standard ribitol for further analysis.
In September, the oven-dried composite sample from a given plot was ground for element analysis, where 5 mg samples were used to analyze the total C, N, and H concentrations using an elemental analysis-stable isotope-ratio mass spectrometer (Vario ELIII; Elementar, Germany); 3 mg samples were used for total O concentration analysis via an Elemental Analyzer (Vario EL Cube; Germany). Total concentrations of P, K, Ca, Mg, S, Fe, Al, Mn, Na, Zn, and Cu were analyzed via acid digestion under high temperature. Briefly, 100 mg samples were introduced into 50 mL beakers, then added 3 mL of nitric acid and 0.5 mL of perchloric acid in them. After left undisturbed overnight (~ 12 h), they were digested on an electric heating plate until the composite liquid becoming clear (~ 6 h). The digested portions were placed into 15-mL flasks, to which distilled water (10 mL) was added. Subsequently, the plasma optical emission spectrometer (ICP-OES) (Iris Advantage 1000; Thermo Jarrell Ash, Franklin, MA) was used to determine concentrations of these elements at the Instrumental Analysis Center, Shanghai Jiao Tong University.
The pathway analysis was carried out with the online software MetaboAnalyst 4.0 (http://www.metaboanalyst.ca/faces/ModuleView.xhtml). According to the results of the pathway analysis, the potential metabolic target pathways, with the value of pathway impact value (PI) > 0.1, were filtered out from the pathway topology. Analysis of variance  was performed to evaluate the differences of metabolites levels between groups with SPSS 20.0 (SPSS Inc., USA). The SIMCA-P version 14.1 (Umetrics, Sweden) was utilized to assess the effects of different types of geological P on the metabolic data and ionomic data by conducting supervised orthogonal partial least squares discrimination analysis (OPLS-DA). Permutation tests (200) or ANOVA of the cross-validated residuals (CV-ANOVA) were conducted to validate the OPLS-DA model. A R2 value of > 0.7 and a Q2 value of > 0.5 in OPLS-DA denoted that the models were highly significant . Furthermore, the variable importance in projection (VIP) was critical for explaining the data obtained with OPLS-DA. Metabolites with a VIP of above 1.0 and a p value of below 0.05, were selected as differential metabolites and discriminating elements (biomarkers), which played a larger role in distinguishing acorns growing at two sites having P-rich and P-deficient soils. Pearson’s correlation test was performed to analyze the relationships among discriminating elements, and relationships of discriminating elements and differential metabolites of acorns in September with SPSS 20.0 (SPSS Inc., USA). All data were log10-transformed to improve normality prior to analysis.
Availability of data and materials
The original dataset of metabolites is available from the Metabolights database (https://www.ebi.ac.uk/metabolights) with study identifier MTBLS888.
Zhu L, Zhu J, Chen Q, Zhu A, Shi Q. The distribution law and controlling factors of rich phosphorite ore in Anning phosphorous deposit, Yunnan Province. Geology Resour. 2008;17:40–4 (In Chinese).
Tian H, Chen G, Zhang C, Melillo JM, Hall CAS. Pattern and variation of C:N:P ratios in China’s soils: a synthesis of observational data. Biogeochemistry. 2010;98:139–51.
Ji H, Du B, Liu C. Elemental stoichiometry and compositions of weevil larvae and two acorn hosts under natural phosphorus variation. Sci Rep. 2017;7:45810.
Zhou X, Sun X, Du B, Yin S, Liu C. Multielement stoichiometry in Quercus variabilis under natural phosphorus variation in subtropical China. Sci Rep. 2015;5:7839.
Reich PB, Oleksyn J. Global patterns of plant leaf N and P in relation to temperature and latitude. PNAS. 2004;101:11001–6.
Sun X, Kang H, Kattge J, Gao Y, Liu C. Biogeographic patterns of multi-element stoichiometry of Quercus varia. Can J For Res. 2015;45:1827–34.
Sun X, Kang H, Chen HY, Berg B, Bartels SF, Liu C. Biogeographic patterns of nutrient resorption from Quercus variabilis Blume leaves across China. Plant Biol. 2016;18:505–13.
Wen J, Ji H, Sun N, Tao H, Du B, Hui D, Liu C. Imbalanced plant stoichiometry at contrasting geologic-derived phosphorus sites in subtropics: the role of microelements and plant functional group. Plant Soil. 2018;430:1–13.
Baxter I. Ionomics: studying the social network of mineral nutrients. Curr Opin Plant Biol. 2009;12:381–6.
Rivas-Ubach A, Peñuelas J. Strong relationship between elemental stoichiometry and metabolome in plants. PNAS. 2012;109:4181–6.
Fan WM. Metabolite profiling by one- and two-dimensional NMR analysis of complex mixtures. Prog Nucl Magn Reson Spectrosc. 1996;28:161–219.
Saito K, Matsuda F. Metabolomics for functional genomics, systems biology, and biotechnology. Annu Rev Plant Biol. 2010;61:463–89.
Iles RA, Hind AJ, Chalmers RA. Use of proton nuclear magnetic resonance spectroscopy in detection and study of organic acidurias. Clin Chem. 1985;31:1795.
Fiehn O. Metabolomics-the link between genotypes and phenotypes. Plant Mol Biol. 2002;48:155–71.
Kociakowska M. Growth and metabolism of cucumber in phosphate-deficient conditions. J Plant Nutr. 2002;25:1–1.
Ji H, Ossipov V, Du B, Wen J, Liu C. Differences in the relationship between metabolomic and ionomic traits of Quercus variabilis growing at contrasting geologic-phosphorus sites in subtropics. Plant Soil. 2019;439:339–55.
Fester T, Merbach I, Schulz E, Härtig C. Metabolic response of Medicago sativa to severe nutrient imbalances and disturbances under field conditions. J Plant Nutr Soil Sci. 2014;177:245–59.
Pitts-Singer TL, Hagen MM, Helm BR, Highland S, Buckner JS, Kemp WP. Comparison of the chemical compositions of the cuticle and dufour’s gland of two solitary bee species from laboratory and field conditions. J Chem Ecol. 2017;43:451–68.
Lahner B, Gong JM, Smith EL, Abid KB, Rogers EE. Genomic scale profiling of nutrient and trace elements in Arabidopsis thaliana. Nat Biotechnol. 2003;21:1215–21.
Baxter I. Should we treat the ionome as a combination of individual elements, or should we be deriving novel combined traits? J Exp Bot. 2015;66:2127–31.
Campos A, Kruijer W, Alexander R, Akkers R, Danku J, Salt D, Aarts M. Natural variation in Arabidopsis thaliana reveals shoot ionome, biomass, and gene expression changes as biomarkers for zinc deficiency tolerance. J Exp Bot. 2017;68:3643–56.
Ding Z, Jia S, Wang Y, Xiao J, Zhang Y. Phosphate stresses affect ionome and metabolome in tea plants. Plant Physiol Bioch Ppb. 2017;120:30.
Yan B, Bing W, Gao YH, Wu JM, Niu JY, Xie YP, Cui ZZ, Zhang ZK. Effects of nitrogen and phosphorus on the regulation of nonstructural carbohydrate accumulation, translocation and the yield formation of oilseed flax. Field Crop Res. 2018;219:229–41.
Bernillon, Stéphane, Benoît, Deborde, Catherine, Maucourt, Mickaël, Cabasson, Cécile. Metabolomic and elemental profiling of melon fruit quality as affected by genotype and environment. Metabolomics 2013, 9: 57–77.
Wu D, Shen Q, Cai S, Chen ZH, Dai F, Zhang G. Ionomic responses and correlations between elements and metabolites under salt stress in wild and cultivated barley. Plant Cell Physiol. 2013;54:1976–88.
Zubek S, Rola K, Szewczyk A, Majewska ML, Turnau K. Enhanced concentrations of elements and secondary metabolites in Viola tricolor L. induced by arbuscular mycorrhizal fungi. Plant Soil. 2015;390:129–42.
Yruela I. Transition metals in plant photosynthesis. Metallomics. 2013;5:1090.
Sanchez D, Pieckenstain F, Escaray F, Erban A, Kraemer U, Udvardi M, Kopka J. Comparative ionomics and metabolomics in extremophile and glycophytic Lotus species under salt stress challenge the metabolic pre-adaptation hypothesis. Plant Cell Environ. 2015;34:605–17.
Jia S, Wang Y, Hu J, Ding Z, Liang Q, Zhang Y, Wang H. Mineral and metabolic profiles in tea leaves and flowers during flower development. Plant Physiol Bioch Ppb. 2016;106:316.
Tan H, Yang X, Zhang F, Zheng X, Qu C, Mu J, Fu F, Li J, Guan R, Zhang H. Enhanced seed oil production in canola by conditional expression of Brassica napus LEAFY COTYLEDON1 and LEC1-LIKE in developing seeds. Plant Physiol. 2011;156:1577–88.
Bewley JD. Seed germination and dormancy. Plant Cell. 1997;9:1055–66.
Aparicio N, Villegas D, Araus JL, Blanco R, Royo C. Seedling development and biomass as affected by seed size and morphology in durum wheat. J Agr Sci. 2002;139:143–50.
Skogerson K, Harrigan GG, Reynolds TL, Halls SC, Ruebelt M, Iandolino A, Pandravada A, Glenn KC, Fiehn O. Impact of genetics and environment on the metabolite composition of maize grain. J Agric Food Chem. 2010;58:3600–10.
Méchin V, Thévenot C, Le GM, Prioul JL, Damerval C. Developmental analysis of maize endosperm proteome suggests a pivotal role for pyruvate orthophosphate dikinase. Plant Physiol. 2007;143:1203.
Neuberger T, Sreenivasulu N, Rokitta M, Rolletschek H, Göbel C, Rutten T, Radchuk V, Feussner I, Wobus U, Jakob P, Webb A, Borisjuk L. Quantitative imaging of oil storage in developing crop seeds. Plant Biotechnol J. 2007;6:31–45.
Hajduch M, Ganapathy A, Stein JW, Thelen JJ. A systematic proteomic study of seed filling in soybean. Establishment of high-resolution two-dimensional reference maps, expression profiles, and an interactive proteome database. Plant Physiol. 2005;137:1397–419.
Nguyen QT, Kisiala A, Andreas P, Emery RJN, Narine S. Soybean seed development: fatty acid and phytohormone metabolism and their interactions. Curr Genomics. 2016;17:241–60.
Wang L, Fu J, Li M, Lena F, Wolfram W, Yang P. Metabolomic and proteomic profiles reveal the dynamics of primary metabolism during seed development of lotus (Nelumbo nucifera). Front Plant Sci. 2016;7:750.
Ji H, Wen J, Du B, Sun N, Berg B, Liu C. Comparison of the nutrient resorption stoichiometry of Quercus variabilis Blume growing in two sites contrasting in soil phosphorus content. Ann Forest Sci. 2018;75:59.
Wang Y. Study on seed developing physiological charavteristics of Tilia miqueiana. China: Nanjing Forestry University; 2010. (In Chinese).
Reynolds TL, Nemeth MA, Glenn KC, Ridley WP, Astwood JD. Natural variability of metabolites in maize grain: differences due to genetic background. J Agric Food Chem. 2005;53:10061–7.
Schupp N, Ziegler P. The relation of starch phosphorylases to starch metabolism in wheat. Plant Cell Physiol. 2004;45:1471–84.
Obenland D, Collin S, Mackey B, Sievert J, Fjeld K, Arpaia ML. Determinants of flavor acceptability during the maturation of navel oranges. Postharvest Biol Tech. 2009;52:156–63.
Cui S. Study on the characteristics of growth and nutrient accumulation of different peanut varieties. China: Hebei Agricultural University; 2015. (In Chinese).
Fait A, Angelovici R, Less H, Ohad I, Urbanczykwochniak E, Fernie AR, Galili G. Arabidopsis seed development and germination is associated with temporally distinct metabolic switches. Plant Physiol. 2006;142:839.
Le BH, Cheng C, Bui AQ, Wagmaister JA, Henry KF, Pelletier J, Kwong L, Belmonte M, Kirkbride R, Horvath S, Drews GN, Fischer RL, Okamuro JK, Harada JJ, Goldberg RB. Global analysis of gene activity during Arabidopsis seed development and identification of seed-specific transcription factors. PANS. 2010;107:8083–70.
Miguel A, Vegabartol JD, Marum L, Chaves I, Santo T, Leitão J, Varela MC, Miguel CM. Characterization of the cork oak transcriptome dynamics during acorn development. BMC Plant Biol. 2015;15:158.
Liu Y, Yang H, Chai H. Studies on the chemical constituents of acorns of some Chinese oaks (1) deciduous species. J Beijing Forestry College. 1984;4 (In Chinese).
Huang Z. Development and component analysis on fruits of energy plants Quercus variabilis and Quercus wutaishanica. China: Beijing Forestry University; 2014. (In Chinese).
Joët T, Laffargue A, Salmona J, Doulbeau S, Descroix F, Bertrand B, Kochko AD, Dussert S. Metabolic pathways in tropical dicotyledonous albuminous seeds: Coffea arabica as a case study. New Phytol. 2009;182:146–62.
Koch KE. Carbohydrate-modulated gene expression in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996;47:509–40.
Wobus U, Weber H. Sugars as signal molecules in plant seed development. Biol Chem. 1999;380:937–44.
Thomas P, Andreas R. Biochemistry and physiology of rafnose family oligosaccharides and galactosyl cyclitols in seeds. Seed Sci Res. 2001;11:185–97.
Chang T, Zhu X. Source-sink interaction: a century old concept under the light of modern molecular systems biology. J Exp Bot. 2017;68:4417–31.
Hildebrandt T, Nunes N, Araújo W, Braun H. Amino acid catabolism in plants. Mol Plant. 2015;8:1563–79.
Galili G, Avin-Wittenberg T, Angelovici R, Fernie A. The role of photosynthesis and amino acid metabolism in the energy status during seed development. Front Plant Sci. 2014;5:447.
Yuan J, Sun N, Du H, Muhammad U, Kang H, Du B, Yin S, Liu C. Correlated metabolic and elemental variations between the leaves and seeds of oak trees at contrasting geologically derived phosphorus sites. Sci Total Environ. 2019;691:178–86.
Zhang Y, Wang Y, Ding Z, Wang H, Song L, Jia S, Ma D. Zinc stress affects ionome and metabolome in tea plants. Plant Physiol Bioch Ppb. 2017;111:318.
Lin Q, Wang C, Dong W, Jiang Q, Wang D, Li S, Chen M, Liu C, Sun C, Chen K. Transcriptome and metabolome analyses of sugar and organic acid metabolism in Ponkan (Citrus reticulata) fruit during fruit maturation. Gene. 2015;554:64–74.
Glew RH, Ayaz FA, Sanz C, Vanderjagt DJ, Huang HS, Chuang LT, Strnad M. Changes in sugars, organic acids and amino acids in medlar (Mespilus germanica L.) during fruit development and maturation. Food Chem. 2003;83:363–9.
Christian B, Oliver F. Metabolite profiling of Chlamydomonas reinhardtii under nutrient deprivation. Plant Physiol. 2005;139:1995–2005.
Gargallogarriga A, Wright SJ, Sardans J. Long-term fertilization determines different metabolomic profiles and responses in saplings of three rainforest tree species with different adult canopy position. PLoS One. 2017;12:e0177030.
Kruger NJ, Schaewen AV. The oxidative pentose phosphate pathway: structure and organisation. Curr Opin Plant Biol. 2003;6:236–46.
Shetty K. Role of proline-linked pentose phosphate pathway in biosynthesis of plant phenolics for functional food and environmental applications: a review. Process Biochem. 2004;39:789–804.
Pan R, Wang X, Li N. Plant physiology. China: Higher Education Press; 2012. (In Chinese).
Lambers H, Finnegan P, Jost R, Plaxton W, Shane M, Stitt M. Phosphorus nutrition in Proteaceae and beyond. Nat Plants. 2015;1:15109.
Groom PK, Lamont BB. Phosphorus accumulation in Proteaceae seeds: a synthesis. Plant Soil. 2010;334:61–72.
Rashid M, Iqbal M. Effect of phosphorus fertilizer on the yield and quality of maize (Zea mays L.) fodder on clay loam soil. J Anim Plant Sci. 2012;22:199–203.
Dyhrman ST, Jenkins BD, Rynearson TA, Saito MA, Mercier ML, Alexander H, Whitney LP, Drzewianowski A, Bulygin VV, Bertrand EM. The transcriptome and proteome of the diatom Thalassiosira pseudonana reveal a diverse phosphorus stress response. PLoS One. 2012;7:1643–6.
Sulpice R, Ishihara H, Schlereth A, Cawthray G, Encke B, Giavalisco P, Ivakov A, Arrivault S, Jost R, Krohn N, et al. Low levels of ribosomal RNA partly account for the very high photosynthetic phosphorus-use efficiency of Proteaceae species. Plant Cell Environ. 2014;37:1276–98.
Lambers H, Cawthray G, Giavalisco P, Kuo J, Laliberté E, Pearse S, Scheible W, Stitt M, Teste F, Turner B. Proteaceae from severely phosphorus-impoverished soils extensively replace phospholipids with galactolipids and sulfolipids during leaf development to achieve a high photosynthetic phosphorus-use-efficiency. New Phytol. 2012;196:1098–108.
Lebzien P. Nutrient elements in grassland: soil–plant–animal relationships. In: Whitehead DC, editor. . Wallingford: CABI Publishing; 2000. p. 209.
Rathore DS, Purohit HS, Yadav BL. Integrated phosphorus management on yield and nutrient uptake of urdbean under rainfed conditions of southern Rajasthan. J Food Legumes. 2010;23:128–31.
Li W, Li L, Sun J, Zhang F, Christie P. Effects of nitrogen and phosphorus fertilizers and intercropping on uptake of nitrogen and phosphorus by wheat, maize, and faba bean. Jl Plant Nutr. 2003;26:629–42.
E S, Yuan J, Ding Z, Yao F, Yu X, Lou F. Effect of N, P, K fertilizers on Fe, Zn, cu, Mn, Ca and mg contents and yields in rice. Chinese J Rice Sci. 2005;19:434–40 (In Chinese).
Lott JNA, Greenwood JS, Batten GD. Mechanisms and regulation of mineral nutrient storage during seed development. In: Kigel J, Galili G, editors. Seed development and germination. New York: Marcel Dekker; 1995. p. 215–35.
Hänsch R, Mendel RR. Physiological functions of mineral micronutrients (cu, Zn, Mn, Fe, Ni, Mo, B, cl). Curr Opin Plant Biol. 2009;12:259–66.
Portsmouth GB. The effect of manganese on carbon assimilation in the potato plant as determined by a modified half-leaf method. Ann Bot. 1949;13:113–33.
Wu CY, Zhu YC, Jiang HC. The vegetation of Yunnan. Beijing: Science Press; 1987. (In Chinese).
Ji H. Stoichiometry, metabolome and genetic variation of Quercus variabilis at contrasting geological phosphorus sites in Central Yunnan. China: Shanghai Jiao Tong University; 2019. (In Chinese).
Merouani H, Apolinario LM, Almeida MH, Pereira JS. Morphological and physiological maturation of acorns of cork oak (Quercus suber L.). Seed Sci Technol. 2003;31:111–24.
Du HM, Wang ZL, Yu WJ, Liu YM, Huang BR. Differential metabolic responses of perennial grass Cynodon transvaalensis×Cynodon dactylon (C4) and Poa Pratensis (C3) to heat stress. Physiol Plantarum. 2011;141:251–64.
Doncheva S, Georgieva K, Vassileva V, Stoyanova Z, Popov NGI. Effects of succinate on manganese toxicity in pea plants. J Plant Nutr. 2005;28:47–62.
Bjerrum JT, Nielsen OH, Hao F, Tang H, Nicholson JK, Wang Y, Olsen J. Metabonomics in ulcerative colitis: diagnostics, biomarker identification, and insight into the pathophysiology. J Proteome Res. 2010;9:954.
We thanks to Wenjuan Yu, Xiaoyu Gao, Jiahui Xu and Shi Xu (Instrumental Analysis Center of Shanghai Jiao tong University) for help with the sample analysis.
This work was funded by the National Key R&D Program of China (2017YFC0505501) and the National Natural Science Foundation of China (31670626, 31270640 and 31070532). The funders supported this research project but played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript. We appreciated the support of the Instrumental Analysis Center, Shanghai Jiao Tong University in the interpretation of data.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Mean monthly temperature and mean monthly precipitation of Mouding and Kunming, from 1981 to 2010.
Variation of the same metabolites in acorns during three developmental stages at P-rich and P-deficient sites. A, Overlap of all metabolites identified in acorns during three developmental stages. R&D, at P-rich and P-deficient sites. B, Heatmap analysis of the same metabolites identified in acorns during three developmental stages. Values were normalized with the Z-transformation. PR, P-rich sites; PD, P-deficient sites; J, July; A, August; S, September.
Over-fitting of OPLS-DA model validation based on metabolomic and ionomic data of developing acorns.
All metabolites identified in acorns at P-rich and P-deficient sites in July, August and September (mg g− 1). --, not determined. All data are mean ± SE (in July, n = 4 of P-rich sites, n = 5 of P-deficient sites; in August and September, n = 8 of P-rich sites, of P-deficient sites, n = 12).
Results of pathway analysis involving all identified metabolites in acorns during three developmental stages. All pathways shown in the table are potential target metabolic pathways with pathway impacts of above 0.1.
Concentrations of acorn elements in September at P-rich and P-deficient sites (mg/g). All data are mean ± SE (n = 8 at P-rich sites; n = 12, at P-deficient sites). Discriminating elements, with bold text for p value and VIP (variable importance projection plot) value, were selected based on the p < 0.05 from the t test (Significance, p < 0.05) and VIP > 1 from OPLS-DA models.
Correlations (r, Pearson’s correlation coefficient) between acorn discriminating elements between P-rich and P-deficient sites in September. Significance, p value < 0.05: *, p < 0.05; **, p < 0.01.
Correlations (r, Pearson’s correlation coefficient) between differential metabolites and discriminating elements in acorns between P-rich and P-deficient sites in September. Significance, p value < 0.05: *, p < 0.05; **, p < 0.01.
About this article
Cite this article
Yuan, J., Sun, N., Du, H. et al. Roles of metabolic regulation in developing Quercus variabilis acorns at contrasting geologically-derived phosphorus sites in subtropical China. BMC Plant Biol 20, 389 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-020-02605-y
- Quercus variabilis
- Phospharite areas
- Acorn development
- Subtropical soils