Shanshan Zhou1,2
Xia Yan1,3*
Jian Yang4
Chaoju Qian1
Xiaoyue Yin1,2
Xingke Fan1,2
Tingzhou Fang1,2
Yuan Gao5
Yuxiao Chang6
Weimin Liu7
Xiao-Fei Ma1,3*- 1Key Laboratory of Stress Physiology and Ecology in Cold and Arid Regions of Gansu Province, Department of Ecology and Agriculture Research, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, China
- 2College of Resources and Environment, University of Chinese Academy of Sciences, Beijing, China
- 3School of Life Sciences, Nantong University, Nantong, China
- 4State Key Laboratory Breeding Base of Dao-di Herbs, National Resource Center for Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, China
- 5Faculty of Materials and Chemical Engineering, Yibin University, Yibin, China
- 6Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, China
- 7Zhongnong Haidao (Shenzhen) Biotech Co., Ltd., Shenzhen, China
Agriophyllum squarrosum (L.) Moq., a pioneer plant endemic to the temperate deserts of Asia, could be domesticated into an ideal crop with outstanding ecological and medicinal characteristics. A previous study showed differential flavonoid accumulation between two in situ altitudinal ecotypes. To verify whether this accumulation was determined by environmental or genetic factors, we conducted flavonoid-targeted metabolic profiling among 14 populations of A. squarrosum collected from regions with different altitudes based on a common garden experiment. Results showed that the most abundant flavonoid in A. squarrosum was isorhamnetin (48.40%, 557.45 μg/g), followed by quercetin (13.04%, 150.15 μg/g), tricin (11.17%, 128.70 μg/g), isoquercitrin (7.59%, 87.42 μg/g), isovitexin (7.20%, 82.94 μg/g), and rutin (7.00%, 80.62 μg/g). However, based on a common garden at middle-altitude environment, almost none of the flavonoids was enriched in the high-altitude populations, and even some flavonoids, such as quercetin, tricin, and rutin, were significantly enriched in low-altitude populations. This phenomenon indicated that the accumulation of flavonoids was not a result of local adaptation to high altitude. Furthermore, association analysis with in situ environmental variables showed that the contents of quercetin, tricin, and rutin were strongly positively correlated with latitude, longitude, and precipitation gradients and negatively correlated with temperature gradients. Thus, we could conclude that the accumulations of flavonoids in A. squarrosum were more likely as a result of local adaption to environmental heterogeneity combined with precipitation and temperature other than high altitude. This study not only provides an example to understand the molecular ecological basis of pharmacognosy, but also supplies methodologies for developing a new industrial crop with ecological and agricultural importance.
Introduction
As sessile organisms, plants have evolved to produce a variety of chemical metabolites in the process of adaptation to the changing environment, and each metabolite plays a vital role in responding to abiotic stresses. On the other hand, the long-term environmental heterogeneity could exert local selection pressure on plants which increases the fitness of individuals in specific environments (Kawecki and Ebert, 2004). When populations inhabit different environments, divergent selection pressure can result in phenotypic differentiation by metabolite accumulation to confer a local fitness advantage (Kawecki and Ebert, 2004). Normally, when exhibit better fitness in a local population, metabolites, particularly secondary metabolites, should be variated along with environmental gradients or multiple stresses, such as drought, high and low temperature, UV-B, pathogens, and herbivores (Vanwallendael et al., 2019). For example, climatic conditions at a high altitude with higher exposure to UV radiation, lower temperature, and stressed conditions that induce the production of secondary metabolites (Barnes et al., 1987; Turunen and Latola, 2005; Zidorn et al., 2005; Spitaler et al., 2008). Many studies have shown that flavonoid content increased with the increasing altitude (Zidorn et al., 2005; Rieger et al., 2008; Kalim et al., 2010; Macar and Kalefetoğlu Macar, 2018; Rana et al., 2020). However, to date, the contribution of a secondary metabolite involved in the high-altitudinal adaptation based on the common garden experiment has seldom been addressed. In addition, breeding of elite germplasms with local adaptation is an important strategy for mitigating the negative impacts of climate change on ecology and agriculture (Howden et al., 2007; Takeda and Matsuoka, 2008). Thus, it is necessary to elucidate the potential role of secondary metabolites for a plant species to adapt along an environmental gradient.
Agriophyllum squarrosum (L.) Moq., also called sand rice, is a pioneer annual psammophyte of the Amaranthaceae family. It is widely distributed in mobile sand and semi-fixed dunes across all the arid and semi-arid regions in Central Asia, the Caucasus, Mongolia, and Siberia, including the Qinghai-Tibet Plateau (Wu et al., 2003). It can survive in extremely high temperatures and drought conditions and can tolerate sand burial (Chen et al., 2014; Li et al., 2015; Zhao et al., 2017). Since the withered plants can reduce wind velocity by more than 90%, they also acted as nutrient providers and reservoirs with a rich source of carbon and nitrogen in infertile soil environments. A. squarrosum can effectively prevent wind and fix sand, and plays a crucial role in maintaining and restoring the fragile desert ecosystems (Ma et al., 2008; Chen et al., 2009). Although A. squarrosum grows in infertile and sandy soils, its seeds are highly nutritious. Compared with Chenopodium quinoa, a world-recognized crop relative with total nutrients (Jacobsen, 2006; Chen et al., 2014), A. squarrosum seeds contain lesser quantities of carbohydrates and more protein and crude fiber (Gao et al., 1991; Ren et al., 2005; Yin et al., 2021), and they are also rich in minerals and trace elements, such as calcium, magnesium, iron, zinc, selenium, and potassium (Zhang et al., 2006; Wang et al., 2009). Besides the high nutritious content in its seeds, the young fresh leaves and stems of A. squarrosum are suitable for fodder and medical usage (Gao, 2002). A Supplement to the Compendium of Materia Medica references that A. squarrosum seeds are healthy for the stomach, spleen, and large intestine (Zhao, 1765), and its aboveground parts can be used as antidiabetes, diuresis, analgesic, and antipyretic in Mongolian folk medicine (Gong et al., 2012). Therefore, as a promising crop candidate with pharmacological and agricultural importance, A. squarrosum provides an ideal model to understand the molecular ecological basis of pharmacognosy.
Previous studies have indicated significant genetic divergence and phenotypic variations among the populations of A. squarrosum (Qian et al., 2016; Yin et al., 2016a,b). Based on common garden trials, the in situ environmental heterogeneity was significantly correlated with phenotypic traits and geographical distribution of A. squarrosum (Yin et al., 2016a,b). In addition, non-targeted metabonomics analysis showed that A. squarrosum is rich in flavonoids (Yin et al., 2018). Flavonoids are ubiquitous plant secondary metabolites that received and continue to receive a great deal of attention in science and medicine, due to the great variety of biological activities and the vast array of biological functions (Winkel-Shirley, 2002). Flavonoids contribute to the adaptation of the plant to environmental stresses, including resistance to diverse abiotic stresses (Liu et al., 2013), such as drought, extreme temperature fluctuations, high soil salinity, and increased exposure to UV radiation (Albert et al., 2018). Moreover, flavonoids play an important role in various plant physiological or developmental functions, as well as in the agronomic and industrial qualities of plant products (Lepiniec et al., 2006; Li et al., 2018). Numerous studies have found that flavonoids have antioxidant, anti-inflammatory, antimicrobial, anti-fungal, antiviral, anti-tumor properties, and also help to scavenge free radicals and stabilize the reactive oxygen species; these keep the body safe from chronic diseases, such as cancer, diabetes, and cardiovascular and heart disease, enhance the immune system, for example, suppress pro-inflammatory immune response, in particular auto-reactive T cells (Havsteen, 2002; Verbeek et al., 2004; Falcone Ferreyra et al., 2012; Mouradov and Spangenberg, 2014; Lei et al., 2019; Sarbu et al., 2019).
Previous studies showed that there was no significant divergence among the altitudinal wild populations in neutral molecular genetic markers (nrITS and chloroplast DNA; Qian et al., 2016, 2021). However, non-target metabolomics analysis supported that the accumulations of a couple of flavonoids, such as hesperetin, quercetin, and apigenin, were significantly enriched in the ecotype with high altitude (unpublished data). Thus, we proposed that the enrichment of flavonoids in the high-altitude populations of A. squarrosum was a consequence of physiological response to environmental stresses or of genetic differentiation involved in local adaptation to the high altitude. To test the above hypothesis, we conducted flavonoid-targeted metabolic profiling among 14 populations of A. squarrosum based on a common garden experiment, a classic approach for studying local adaptation (Anderson et al., 2010; Fournier-Level et al., 2011; Savolainen et al., 2013). Thus, we aimed to address the following questions: (1) Whether flavonoids are still enriched in the high-altitude populations of A. squarrosum at the common garden site. If the results obey the trend, we can conclude that the flavonoid accumulated due to the genetic basis of local adaptation to the high altitude. (2) If not, what other kinds of in situ environmental variables influenced the differences in flavonoid metabolites among the populations. Under this circumstance, the differentiation of flavonoids could be determined by the heterogeneity of precipitation and temperature among the populations, which could also be a result of local adaptation for the nature populations. This study could not only shed light on the molecular ecological basis for A. squarrosum adapted to the desert heterogeneity, but also provide instructive guidelines for the development and utilization of the wild plant in food, pharmaceutical, and cosmetic industries under global climate change.
Materials and Methods
Plant Materials and Common Garden Experiment
Seeds of 14 wild A. squarrosum populations, representing the different altitude regions of A. squarrosum in China, were sampled in 2016 (Figure 1; Table 1). According to the altitude, the 14 A. squarrosum populations were divided into three groups: low-altitude (~0–1,000 m), middle-altitude (~1,000–2,000 m), and high-altitude (~2,000–4,000 m) groups. To identify the accumulation of flavonoid metabolites in these A. squarrosum populations, a common garden experiment was conducted in natural environmental conditions in Wuwei City, Gansu Province (WW, 37°54′10.98″N, 102°54′4.2″E, 1530 m) at the southern edge of the Tengger Desert (Figure 1). Seeds of A. squarrosum were sown on April 20, 2018, and seedlings initially raised in a nursery bed were transplanted into the common garden after 10 days. Twenty randomly selected seedlings of each population were planted at 2-m intervals. Fresh tissues of A. squarrosum, including leaves, stems, and spikes, were randomly collected from six individuals of each population in early August 2018 at the flowering phase. All the materials were frozen immediately in liquid nitrogen and then stored at −80°C for further use.
Figure 1. Geographical distribution of 14 Agriophyllum squarrosum populations in China. The red dot indicates the planting place of the common garden experiment.
Table 1. Details of the locality information of the sampled populations.
Flavonoid-Targeted Metabolic Profiling
The samples of A. squarrosum were dried at 60°C in a drying oven after washing. For each sample, 50 mg homogenized sample powder was suspended in 1.5 ml 70% v/v ethanol overnight. Then, the extract solution was extracted by ultrasonication for 30 min at a 40 KHz SCIENT ultrasonic processor (Ningbo, China). The extract was centrifuged for 10 min at 12,000 × g and the supernatant was subsequently filtered through a 0.22-μm syringe filter for UPLC-QQQ-MS/MS analysis. Chromatographic analysis was performed using the Acquity UPLCT I-Class system (Waters, Milford, MA, United States). The column used was Acquity UPLC BEH C18 column (2.1 mm × 100 mm, 1.8 μm), and the column temperature was maintained at 40°C. The binary gradient consisted of solvent system A (formic acid/water, 0.1:99.9, v/v) and solvent system B (formic acid/acetonitrile, 0.1:99.9, v/v). The chromatographic conditions were as follows: 0 min, 5% B; 1.0 min, 25% B; 3.5 min, 40% B; and 4.5 min, 60% B. The injection volume was 1.0 μl, and the flow rate was 0.60 ml/min. Tandem mass spectrometry (MS/MS) was performed using QTRAP 6500 system (AB SCIEX, Los Angeles, CA, United States) equipped with an electrostatic ionization (ESI) source (AB SCIEX). The MS spectra were acquired in negative ion mode, which was carried out by optimization of the product ion obtained from the fragment of the isolated precursor ion for each analyte. The ion spray potential was −4,500 V, and the source temperature was set at 550°C. After the product ions were chosen, the multiple reaction monitoring conditions for each standard were further optimized to achieve maximum sensitivity. The retention times (tR), quantitative ion pairs, declustering potential (DP), collision energy (CE), and cell exit potential (CXP) are listed in Table 2.
Table 2. Mass spectra properties of flavonoids evaluated in A. squarrosum.
Environmental Variables
To test the correlations of different flavonoid metabolites with environmental variables, we extracted 22 environmental variables of the original locations of the 14 A. squarrosum populations from the world climate database1 (Supplementary Table S1), including altitude, annual mean temperature (AMT), mean monthly temperature range (MMTR), isothermality, temperature seasonality (TS), maximum temperature of the warmest month (MATWM), minimum temperature of the coldest month (MITCM), temperature annual range (TAR), mean temperature of the wettest quarter (MTWEQ), mean temperature of the driest quarter (MTDQ), mean temperature of the warmest quarter (MTWAQ), mean temperature of the coldest quarter (MTCQ), annual precipitation (AP), precipitation of the wettest month (PWM), precipitation of the driest month (PDM), precipitation seasonality (PS), precipitation of the wettest quarter (PWEQ), precipitation of the driest quarter (PDQ), precipitation of the warmest quarter (PWAQ), precipitation of the coldest quarter (PCQ), longitude, and latitude.
Statistical Analyses
To determine the content differences in flavonoids among the examined groups, one-way ANOVA, Tukey honest significant difference (HSD) test, and pairwise Pearson’s correlations were performed using IBM SPSS Statistics (version 25; IBM Corporation, Armonk, NY, United States) and OriginPro2020 (OriginLab Corporation, Northampton, MA, United States). To discern the similarities and differences in the profiles of the flavonoids among these A. squarrosum populations, principal component analysis (PCA) was performed using SIMCA software (version 13.0; Umetrics AB, Umeå, Sweden). To maximize the separation between the different groups, the supervised models were subsequently formulated by orthogonal partial least-squares discrimination analysis (OPLS-DA). To show the variation in the flavonoid metabolites in different A. squarrosum populations, a cluster heatmap was implemented by the heatmap package in R (version 1.1.463; RStudio, Boston, MA, United States). To investigate the relationship between environmental variables and flavonoid metabolites in A. squarrosum, univariate linear regression analysis was performed using OriginPro2020.
Results
Composition and Content of Flavonoid Metabolites in A. squarrosum
As shown in Figure 2 and Supplementary Table S2, 20 flavonoid metabolites were identified from 14 A. squarrosum populations. Figures 2A,B show the averages from 14 A. squarrosum populations. The flavonoid metabolites content in A. squarrosum varied from 0.04 μg/g (luteoloside, 0.003%) to 557.45 μg/g (isorhamnetin, 48.40%). Quercetin (13.04%, 150.15 μg/g), tricin (11.17%, 128.70 μg/g), isoquercitrin (7.59%, 87.42 μg/g), isovitexin (7.20%, 82.94 μg/g), rutin (7.00%, 80.62 μg/g), and isorhamnetin-3-glucoside (3.73%, 43.01 μg/g) were highly accumulated in A. squarrosum (Figures 2A,B; Supplementary Table S2b).
Figure 2. Flavonoid content in A. squarrosum. The composition (A) and content (B) of flavonoids in A. squarrosum. Values are Mean ± SE (n = 84). (C) The composition and content of flavonoids in 14 A. squarrosum populations. Values are Mean ± SE (n = 6). (D) The composition and content of flavonoids in different altitude populations. Values are Mean ± SE (n = 12, 54, 18). Values marked with the different or same letters indicate statistically significant and not significant differences, respectively, using Tukey honest significant difference (HSD) test at the 5% significance level.
Among the different populations, the total flavonoid content ranged from 754.74 μg/g [population in Tiegaixiang, Qaidam (TGX)] to 1548.62 μg/g [population in Xiangshawan, Kubuqi (XSW)], averaging 1151.70 μg/g (Figure 2C; Supplementary Table S2c). The level of quercetin, tricin, and rutin stood out with the higher content in low-altitude group than in the high- and middle-altitude groups (Figure 2D; Supplementary Table S2d).
Differences in Flavonoid Metabolites Among A. squarrosum Populations
After a comprehensive view of the similarities and differences in the profiles of the flavonoids among A. squarrosum populations, PCA and OPLS-DA were performed. The unsupervised PCA model was also performed to observe the trends of flavonoid metabolites in A. squarrosum (Figure 3A). The principal components 1 (PC1, 22.9%) and 2 (PC2, 17.7%) explain 45.2% of the metabolic variance among samples. As shown in Figure 3A, PCA did not show a clear separation among the three altitude groups. To identify the variables contributing to the separation between the two sample groups, a supervised OPLS-DA was performed. OPLS-DA between the high-altitude and low-altitude groups gave a strong model with R2 = 0.714 and Q2 = 0.777, and the plot showed complete separation of the two groups (Figure 3B). The total variation in the high-altitude group vs. the low-altitude group was 35.6%, of which 13.3% was related to the separation between the high-altitude and low-altitude groups, while 22.1% was within group variation. OPLS-DA between the high-altitude and middle-altitude groups gave a good model with R2 = 0.422 and Q2 = 0.586, and the plot showed definite groupings with a small overlap (Figure 3C). Class separation was not as clear as in the model middle-altitude vs. low-altitude sample (Figure 3D). These results indicated that the accumulation of 20 flavonoids was found to be significantly different between the high-altitude group and the other two altitude groups, while not significantly different between the middle-altitude and low-altitude groups.
Figure 3. Score plots of flavonoid metabolites in A. squarrosum. (A) Principal component analysis score plot of total samples (R2 = 0.559, Q2 = 0.213). (B) Orthogonal partial least-squares discrimination analysis (OPLS-DA) score plot of high-altitude vs. low-altitude sample (R2X = 0.714, R2Y = 0.917, Q2 = 0.777). (C) OPLS-DA score plot of high-altitude vs. middle-altitude sample (R2X = 0.422, R2Y = 0.586, Q2 = 0.397). (D) OPLS-DA score plot of middle-altitude vs. low-altitude sample (R2X = 0.316, R2Y = 0.404, Q2 = 0.227). Green, blue, and red circles indicate sample sites of low-altitude (~0–1,000 m), middle-altitude (~1,000–2,000 m), and high-altitude (~2,000–4,000 m), respectively.
The level of isorhamnetin was more than 2-fold higher than that of the other 19 flavonoid metabolites, and this prevented the cluster heatmap from showing the differences in flavonoid metabolites in the different populations. So, we produced the cluster heatmap of 19 flavonoid metabolites except for isorhamnetin among 14 A. squarrosum populations. Figure 4 shows that the content of flavonoid metabolites varied among different populations. As shown in Figure 4, the 19 flavonoid metabolites were clustered into three distinct metabolite groups. The first group contained high levels of three flavonoid metabolites, including quercetin, tricin, and rutin. The second group contained medium levels of two flavonoid metabolites, including isoquercitrin and isovitexin. The last group contained low levels of 14 flavonoids, such as isorhamnetin-3-glucoside, astragaline, and so on.
Figure 4. Cluster heatmap of flavonoid metabolites in A. squarrosum. Metabolites and populations two-way clustering heatmap. Color depth represents average intensity of metabolite contents in 14 populations of A. squarrosum. Red represents highest content and blue represents lowest content. Altitude categories are encoded using three colors: green, blue, and red indicates sample sites of low-altitude (~0–1,000 m), middle-altitude (~1,000–2,000 m), and high-altitude (~2,000–4,000 m), respectively.
As shown in Figure 4, populations in Dulan, Qaidam, and TGX of the high-altitude group had lower values in most of the measured flavonoids. The middle-altitude group was characterized by medium levels and the low-altitude group had high values in most of the measured flavonoids. These results indicated that most flavonoid metabolites, especially quercetin, tricin, and rutin, were more in the low-altitude populations than in the high- and middle-altitude groups (Figures 2D, 4).
Relationship Between Environmental Variables and Flavonoid Metabolites in A. squarrosum
Univariate linear regression was implemented to identify the major in situ environmental variables determining the flavonoid metabolites in the A. squarrosum (Figure 5; Supplementary Figures S1–S22). We found that quercetin, tricin, and rutin were strongly positively correlated with latitude, longitude, and precipitation gradients, such as AP, PWAQ, PWEQ, and PWM. While rutin was significantly negatively correlated with temperature gradients, such as MTDQ and mean MTCQ, quercetin was significantly negatively correlated with isothermality. Naringenin, isovitexin, and isorhamnetin-3-glucoside were mainly positively correlated with temperature gradients, such as AMT, MITCM, MTDQ, and MTCQ. Isorhamnetin and hesperitin had significantly positive correlationship with MITCM. Eriodictyol had a significantly positive correlationship with AMT and MTWEQ. These results indicated that naringenin, isovitexin, isorhamnetin-3-glucoside, hesperitin, and eriodictyol were mainly positively correlated with temperature, while rutin and quercetin were significantly negatively correlated with temperature.
Figure 5. Cluster heatmap based on correlation between flavonoid metabolites and the environmental variables. Red represents positive correlation, blue represents negative correlation. *Indicates statistical significance (p < 0.05).
Discussion
Flavonoids play an important role in the agronomic and industrial qualities of plant products with tremendous medicinal value (Lepiniec et al., 2006; Li et al., 2018; Sarbu et al., 2019). Our results showed that isorhamnetin was the most abundant flavonoid among 20 identified flavonoid metabolites in A. squarrosum, accounted for more than 48.40%, averaging 0.56 mg/g in our study (Figures 2A,B), which is equivalent to the isorhamnetin content in another medicinal plant kale (Schmidt et al., 2010). Isorhamnetin has a wide range of pharmacological effects on cardiovascular, hyperuricemia, and pulmonary fibrosis diseases and a variety of tumors, and possesses the potential of preventing neurodegenerative diseases such as Alzheimer’s (Gong et al., 2020). As an annual medicinal psammophyte with large aboveground biomass, A. squarrosum provides ecological benefits of windbreak and sand fixation, making good use of barren sandy land. So, A. squarrosum possesses great potential for promotion and application, and its high isorhamnetin content makes it attractive to the pharmaceutical and food industries.
As sessile organisms, plants are often exposed to various environmental stress factors. Among them, abiotic stresses (such as cold, heat, drought, salinity, and UV radiation) are major constraints that affect plant development and growth and pose serious threats to plants’ life. Hence, plants must regulate their growth and development in response to abiotic stresses through primary and secondary metabolism (Cetinkaya et al., 2016). Some secondary metabolite accumulation was positively influenced by abiotic stress conditions, as a consequence of plant adaptive strategies concerned with the establishment of some changes allowing the plant to sustain its life under abiotic stress conditions. Among them, flavonoids could respond to the different stressors, which seems a universal mechanism that flavonoid accumulation is enhanced under abiotic stress (Vanderauwera et al., 2005; Walia et al., 2005; Sun et al., 2020). On the other hand, long-term selection can lead to the development of morphological and physiological adaptations to the local environment, generating genetic differentiation in important traits, including specific metabolites (Kawecki and Ebert, 2004; Savolainen et al., 2007). Individuals within a species may diversify by orders of magnitude in chemical constituents, sizes, and reproductive and growth rates (Callaway et al., 2003). The processes involved in the metabolism and accumulation of chemical constituents in plants are under both genetic and environmental control. Therefore, it is necessary to conduct a common garden experiment across the populations with heterogeneity, which could minimize the environmental effect on the accumulation of the metabolites.
Agriophyllum squarrosum populations distributed across a wide range of environmental conditions are subjected to differential selective pressures. Our recent study showed that flavonoid metabolites were significantly enriched in the in situ high-altitude ecotype than in the in situ middle-altitude ecotypes of A. squarrosum (unpublished data). However, in this common garden experiment, we did not find that flavonoid metabolites were significantly enriched in the high-altitude populations, and even some flavonoids, such as quercetin, tricin, and rutin, were significantly enriched in the low-altitude populations (Figures 2D, 4). Considering that the common garden experiment was performed in the middle-altitude environment where the environmental stress for the high-altitude population was decreased, but environmental stress for the low-altitude populations was increased, accordingly, flavonoids were significantly enriched in the low-altitude populations than in the high- and middle-altitude populations. Thus, we proposed that the accumulation of flavonoids was not related to local adaptation to high altitude.
Based on stress treatment, some studies showed that low temperature would increase the accumulation of quercetin (Albert et al., 2009; Jaakola and Hohtola, 2010). However, another study found that high temperature can increase rutin concentration (Lumingkewas et al., 2015). Some field sampling studies have found that the concentration of quercetin derivates correlated positively with the in situ latitude (Stark et al., 2008; Sun et al., 2020). Theoretically, the environment of parents may influence characters, such as phenotype and physiological activity, in the offspring (Galloway, 2001). However, without common garden experiments across populations, it is hard to eliminate the environmental imprint. In this study, through investigating the correlations of 20 flavonoid metabolites with environmental variables from the original locations of 14 A. squarrosum populations, we detected that quercetin, tricin, and rutin were strongly positively correlated with latitude, longitude, and precipitation gradients, such as AP, PWAQ, PWEQ, and PWM. While rutin was significantly negatively correlated with temperature gradients, such as MTDQ and MTCQ; quercetin was significantly negatively correlated with isothermality (Figure 5). This phenomenon suggested that flavonoids accumulated among the populations of A. squarrosum as a result of natural selection of environmental heterogeneity. Of course, more treatment experiments should be conducted to verify this point. These findings could guide us to screen the germplasm and planting field with proper precipitation and temperatures to guarantee maximized yield of flavonoid in A. squarrosum.
Conclusion
This study showed a high variation of the 20 measured flavonoid metabolites among 14 A. squarrosum populations from different altitudinal regions based on a common garden experiment. Among them, the contents of quercetin, tricin, and rutin were significantly enriched in low-altitude populations than in middle- and high-altitude, which were positively correlated to precipitation of the in situ habitat. This study could shed light on guidelines for future breeding and field management of A. squarrosum.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
SZ: writing–original draft, software, and visualization. JY, XYi, XF, TF, and YG: investigation. CQ, YC, and WL: writing–review and editing. X-FM and XYa: conceptualization, methodology, supervision, and writing–review and editing. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Key Research and Development Program of China (grant nos. 2016YFC0500903 and 2017YF100200); the National Natural Science Foundation of China (grant nos. 31901079 and 31770416); the Key Research and Development Program of Gansu Province (grant no. 18YF1NA123); the Xinjiang Production and Construction Crops Key Laboratory of Protection and Utilization of Biological Resources in Tarim Basin (grant no. BRZD1809); the Foundation for Excellent Youth Scholars of NIEER, CAS (grant no. Y92972); and the Special Funds for Scientific and Technological Innovation and Development in Gansu Province (Research and Development Center of Desert Nutrition and Medicinal Plant Germplasm Innovation).
Conflict of Interest
WL was employed by the company Zhongnong Haidao (Shenzhen) Biotech Co., Ltd.The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.683265/full#supplementary-material
Footnotes
References
Albert, A., Sareedenchai, V., Heller, W., Seidlitz, H. K., and Zidorn, C. (2009). Temperature is the key to altitudinal variation of phenolics in Arnica montana L. cv. ARBO. Oecologia 160, 1–8. doi: 10.1007/s00442-009-1277-1
Albert, N. W., Thrimawithana, A. H., McGhie, T. K., Clayton, W. A., Deroles, S. C., Schwinn, K. E., et al. (2018). Genetic analysis of the liverwort Marchantia polymorpha reveals that R2R3MYB activation of flavonoid production in response to abiotic stress is an ancient character in land plants. New Phytol. 218, 554–566. doi: 10.1111/nph.15002
Anderson, C. D., Epperson, B. K., Fortin, M.-J., Holderegger, R., James, P. M. A., Rosenberg, M. S., et al. (2010). Considering spatial and temporal scale in landscape-genetic studies of gene flow. Mol. Ecol. 19, 3565–3575. doi: 10.1111/j.1365-294X.2010.04757.x
Barnes, P. W., Flint, S. D., and Caldwell, M. M. (1987). Photosynthesis damage and protective pigments in plants from a latitudinal arctic/alpine gradient exposed to supplemental UV-B radiation in the field. Arct. Alp. Res. 19, 21–27. doi: 10.2307/1550996
Callaway, R., Pennings, S., and Richards, C. (2003). Phenotypic plasticity and interactions among plants. Ecology 84, 1115–1128. doi: 10.1890/0012-9658(2003)084[1115:PPAIAP]2.0.CO;2
Cetinkaya, H., Koc, M., and Kulak, M. (2016). Monitoring of mineral and polyphenol content in olive leaves under drought conditions: application chemometric techniques. Ind. Crop. Prod. 88, 78–84. doi: 10.1016/j.indcrop.2016.01.005
Chen, B. M., Wang, G. X., and Peng, S. L. (2009). Role of desert annuals in nutrient flow in arid area of Northwestern China: a nutrient reservoir and provider. J. Plant Ecol. 201, 401–409. doi: 10.1007/s11258-008-9526-7
Chen, G., Zhao, J., Zhao, X., Zhao, P., Duan, R., Nevo, E., et al. (2014). A psammophyte Agriophyllum squarrosum (L.) Moq.: a potential food crop. Genet. Resour. Crop. Evol. 61, 669–676. doi: 10.1007/s10722-014-0083-8
Falcone Ferreyra, M. L., Rius, S. P., and Casati, P. (2012). Flavonoids: biosynthesis, biological functions, and biotechnological applications. Front. Plant Sci. 3:222. doi: 10.3389/fpls.2012.00222
Fournier-Level, A., Korte, A., Cooper, M. D., Nordborg, M., Schmitt, J., and Wilczek, A. M. (2011). A map of local adaptation in Arabidopsis thaliana. Science 334, 86–89. doi: 10.1126/science.1209271
Galloway, L. F. (2001). Parental environmental effects on life history in the herbaceous plant Campanula americana. Ecology 82, 2781–2789. doi: 10.2307/2679960
Gao, Q. (2002). The “grass seed” is Agriophyllum squarrosum in Dunhuang manuscripts. J. Dunhuang Stud. 42, 43–44.
Gao, B., Sun, H., and Liu, H. (1991). The investigation and utilization of nutritional value of Agriophyllum squarrosum seeds. Food Sci. 12, 50–53.
Gong, G., Guan, Y.-Y., Zhang, Z.-L., Rahman, K., Wang, S.-J., Zhou, S., et al. (2020). Isorhamnetin: a review of pharmacological effects. Biomed. Pharmacother. 128:110301. doi: 10.1016/j.biopha.2020.110301
Gong, B., Zhan, K., Zhou, Y., Zhang, L., Hui, Y., and Li, Y. (2012). Separation and identification of chemical constituents from Agriophyllum squarrosum (L.) Moq. Chin. Med. Mat. 14, 7–11. doi: 10.13313/j.issn.1673-4890.2012.10.016
Havsteen, B. H. (2002). The biochemistry and medical significance of the flavonoids. Pharmacol. Ther. 96, 67–202. doi: 10.1016/S0163-7258(02)00298-X
Howden, S. M., Soussana, J.-F., Tubiello, F. N., Chhetri, N., Dunlop, M., and Meinke, H. (2007). Adapting agriculture to climate change. Proc. Natl. Acad. Sci. U. S. A. 104:19691. doi: 10.1073/pnas.0701890104
Jaakola, L., and Hohtola, A. (2010). Effect of latitude on flavonoid biosynthesis in plants. Plant Cell Environ. 33, 1239–1247. doi: 10.1111/j.1365-3040.2010.02154.x
Jacobsen, S.-E. (2006). The Worldwide Potential for Quinoa (Chenopodium quinoa Willd.). Food Rev. Int. 19, 167–177. doi: 10.1081/FRI-120018883
Kalim, M. D., Bhattacharyya, D., Banerjee, A., and Chattopadhyay, S. (2010). Oxidative DNA damage preventive activity and antioxidant potential of plants used in Unani system of medicine. BMC Complement. Altern. Med. 10:77. doi: 10.1186/1472-6882-10-77
Kawecki, T. J., and Ebert, D. (2004). Conceptual issues in local adaptation. Ecol. Lett. 7, 1225–1241. doi: 10.1111/j.1461-0248.2004.00684.x
Lei, Z., Sumner, B. W., Bhatia, A., Sarma, S. J., and Sumner, L. W. (2019). UHPLC-MS analyses of plant flavonoids. Curr. Protoc. Plant Biol. 4:e20085. doi: 10.1002/cppb.20085
Lepiniec, L., Debeaujon, I., Routaboul, J.-M., Baudry, A., Pourcel, L., Nesi, N., et al. (2006). Genetics and biochemistry of seed flavonoids. Annu. Rev. Plant Biol. 57, 405–430. doi: 10.1146/annurev.arplant.57.032905.105252
Li, Y., Fang, J., Qi, X., Lin, M., Zhong, Y., Sun, L., et al. (2018). Combined analysis of the fruit metabolome and transcriptome reveals candidate genes involved in flavonoid biosynthesis in Actinidia arguta. Int. J. Mol. Sci. 19:1471. doi: 10.3390/ijms19051471
Li, J., Qu, H., Zhao, H., Zhou, R., Yun, J., and Pan, C. (2015). Growth and physiological responses of Agriophyllum squarrosum to sand burial stress. J. Arid. Land 7, 94–100. doi: 10.1007/s40333-014-0033-5
Liu, Z., Liu, Y., Pu, Z., Wang, J., Zheng, Y., Li, Y., et al. (2013). Regulation, evolution, and functionality of flavonoids in cereal crops. Biotechnol. Lett. 35, 1765–1780. doi: 10.1007/s10529-013-1277-4
Lumingkewas, A. M. W., Koesmaryono, Y., and Aziz, S., Impron (2015). The influence of temperature to rutin concentration of buckwheat grains in humid tropic. Int. J. Sci. Basic Appl. Res. 20, 1–9.
Ma, Q. L., Wang, J. H., Zhang, J. G., Zhan, K. J., Zhang, D. K., and Chen, F. (2008). Ecological protective function of a pioneer species (Agriophyllum squarrosum) on shifting sand dunes. J. Soil Water Conserv. 22, 140–142.
Macar, O., and Kalefetoğlu Macar, T. (2018). Altitude triggers some biochemical adaptations of Polygonum cognatum Meissn. Plants. Cumhur. Sci. J. 39, 621–627. doi: 10.17776/csj.375050
Mouradov, A., and Spangenberg, G. (2014). Flavonoids: a metabolic network mediating plants adaptation to their real estate. Front. Plant Sci. 5:620. doi: 10.3389/fpls.2014.00620
Qian, C., Yan, X., Fang, T., Yin, X., Zhou, S., Fan, X., et al. (2021). Genomic adaptive evolution of sand rice (Agriophyllum squarrosum) and its implications for desert ecosystem restoration. Front. Genet. 12:656061. doi: 10.3389/fgene.2021.656061
Qian, C., Yin, H., Shi, Y., Zhao, J., Yin, C., Luo, W., et al. (2016). Population dynamics of Agriophyllum squarrosum, a pioneer annual plant endemic to mobile sand dunes, in response to global climate change. Sci. Rep. 6:26613. doi: 10.1038/srep26613
Rana, P. S., Saklani, P., and Chandel, C. (2020). Influence of altitude on secondary metabolites and antioxidant activity of Coleus forskohlii root extracts. Res. J. Med. Plant 14, 43–52. doi: 10.3923/rjmp.2020.43.52
Ren, W., Liu, X., and Ni, C. (2005). An analysis on nutritional composition of natural Agriophyllum squarrosum of Maowusu desert. J. Inner Mong. Inst. Agric. Anim. Husb. 26, 88–90.
Rieger, G., Müller, M., Guttenberger, H., and Bucar, F. (2008). Influence of altitudinal variation on the content of phenolic compounds in wild populations of Calluna vulgaris, Sambucus nigra, and Vaccinium myrtillus. J. Agric. Food Chem. 56, 9080–9086. doi: 10.1021/jf801104e
Sarbu, L., Bahrin, L. G., Babii, C., Stefan, M., and Birsa, L. (2019). Synthetic flavonoids with antimicrobial activity: a review. J. Appl. Microbiol. 127, 1282–1290. doi: 10.1111/jam.14271
Savolainen, O., Lascoux, M., and Merilä, J. (2013). Ecological genomics of local adaptation. Nat. Rev. Genet. 14, 807–820. doi: 10.1038/nrg3522
Savolainen, O., Pyhäjärvi, T., and Knürr, T. (2007). Gene flow and local adaptation in trees. Annu. Rev. Ecol. Evol. Syst. 38, 595–619. doi: 10.1146/annurev.ecolsys.38.091206.095646
Schmidt, S., Zietz, M., Schreiner, M., Rohn, S., Kroh, L. W., and Krumbein, A. (2010). Genotypic and climatic influences on the concentration and composition of flavonoids in kale (Brassica oleracea var. sabellica). Food Chem. 119, 1293–1299. doi: 10.1016/j.foodchem.2009.09.004
Spitaler, R., Winkler, A., Lins, I., Yanar, S., Stuppner, H., and Zidorn, C. (2008). Altitudinal variation of phenolic contents in flowering heads of Arnica montana cv. ARBO: a 3-year comparison. J. Chem. Ecol. 34, 369–375. doi: 10.1007/s10886-007-9407-x
Stark, S., Julkunen-Tiitto, R., Holappa, E., Mikkola, K., and Nikula, A. (2008). Concentrations of foliar quercetin in natural populations of white birch (Betula pubescens) increase with latitude. J. Chem. Ecol. 34, 1382–1391. doi: 10.1007/s10886-008-9554-8
Sun, C., Shang, X., Ding, H., Cao, Y., and Fang, S. (2020). Natural variations in flavonoids and triterpenoids of Cyclocarya paliurus leaves. J. For. Res. 32, 805–814. doi: 10.1007/s11676-020-01139-1
Takeda, S., and Matsuoka, M. (2008). Genetic approaches to crop improvement: responding to environmental and population changes. Nat. Rev. Genet. 9, 444–457. doi: 10.1038/nrg2342
Turunen, M., and Latola, K. (2005). UV-B radiation and acclimation in timberline plants. Environ. Pollut. 137, 390–403. doi: 10.1016/j.envpol.2005.01.030
Vanderauwera, S., Zimmermann, P., Rombauts, S., Vandenabeele, S., Langebartels, C., Gruissem, W., et al. (2005). Genome-wide analysis of hydrogen peroxide-regulated gene expression in Arabidopsis reveals a high light-induced transcriptional cluster involved in anthocyanin biosynthesis. Plant Physiol. 139:806. doi: 10.1104/pp.105.065896
Vanwallendael, A., Soltani, A., Emery, N., Peixoto, M., Olsen, J., and Lowry, D. (2019). A molecular view of plant local adaptation: incorporating stress-response networks. Annu. Rev. Plant Biol. 70, 559–583. doi: 10.1146/annurev-arplant-050718-100114
Verbeek, R., Plomp, A. C., van Tol, E. A. F., and Van Noort, J. M. (2004). The flavones luteolin and apigenin inhibit in vitro antigen-specific proliferation and interferon-gamma production by murine and human autoimmune T cells. Biochem. Pharmacol. 68, 621–629. doi: 10.1016/j.bcp.2004.05.012
Walia, H., Wilson, C., Condamine, P., Liu, X., Ismail, A. M., Zeng, L., et al. (2005). Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage. Plant Physiol. 139:822. doi: 10.1104/pp.105.065961
Wang, Y., Zhao, P., Li, Q., Lin, Y., and Qi, X. (2009). Evaluation on nutritional composition of Agriophyllum squarrosum of Tengger desert. Sci. Technol. Food Ind. 30, 286–288.
Winkel-Shirley, B. (2002). Biosynthesis of flavonoids and effects of stress. Curr. Opin. Plant Biol. 5, 218–223. doi: 10.1016/S1369-5266(02)00256-X
Wu, Z., Raven, P. H., and Hong, D. (2003). Flora of China: Ulmaceae Through Basellaceae. Vol. 5. Beijing and St. Louis: Science Press and Missouri Botanical Garden Press, 367.
Yin, C., Qian, C., Chen, G., Yan, X., and Ma, X. (2016a). The influence of selection of ecological differentiation to the phenotype polymorphism of Agriophyllum squarrosum. J. Desert Res. 36, 364–373.
Yin, X., Wang, W., Qian, C., Fan, X., Yan, X., Yan, G., et al. (2018). Analysis of metabolomics in Agriophyllum squarrosum based on UPLC-MS. Chin. J. Exp. Tradit. Med. Form. 24, 51–56. doi: 10.13422/j.cnki.syfjx.20181503
Yin, X., Yan, X., Qian, C., Zhou, S., Fang, T., Fan, X., et al. (2021). Comparative transcriptome analysis to identify genes involved in terpenoid biosynthesis in Agriophyllum squarrosum, a folk medicinal herb native to Asian temperature deserts. Plant Biotechnol. Rep. doi: 10.1007/s11816-021-00674-5
Yin, C. L., Zhao, J. C., Jinling, H. U., Qian, C. J., Xie, Z. K., Chen, G. X., et al. (2016b). Phenotypic variation of a potential food crop, Agriophyllum squarrosum, impacted by environmental heterogeneity. Sci. Sinica 46, 1324–1335. doi: 10.1360/N052015-00294
Zhang, J., Zhao, J., and Li, J. (2006). Determination and analysis of seed nutrients of Agriophyllum squarrosum. Pratacult. Sci. 23, 77–79.
Zhao, X. (1765). The Supplement to Compendium of Materia Medica. Beijing: People’s medical Publishing House.
Zhao, P., Zhang, J., Qian, C., Zhou, Q., Zhao, X., Chen, G., et al. (2017). SNP discovery and genetic variation of candidate genes relevant to heat tolerance and agronomic traits in natural populations of sand rice (Agriophyllum squarrosum). Front. Plant Sci. 8:536. doi: 10.3389/fpls.2017.00536
Keywords: Agriophyllum squarrosum, flavonoid, targeted metabolic profiling, common garden experiment, local adaptation
Citation: Zhou S, Yan X, Yang J, Qian C, Yin X, Fan X, Fang T, Gao Y, Chang Y, Liu W and Ma X-F (2021) Variations in Flavonoid Metabolites Along Altitudinal Gradient in a Desert Medicinal Plant Agriophyllum squarrosum. Front. Plant Sci. 12:683265. doi: 10.3389/fpls.2021.683265
Edited by:
Guiguo Zhang, Shandong Agricultural University, ChinaReviewed by:
Nutan Kaushik, Amity University, IndiaXuezhao Sun, Jilin Agricultural Science and Technology University, China
Copyright © 2021 Zhou, Yan, Yang, Qian, Yin, Fan, Fang, Gao, Chang, Liu and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiao-Fei Ma, maxiaofei@lzb.ac.cn; Xia Yan, yanxia@lzb.ac.cn