ORIGINAL RESEARCH article

Front. Plant Sci., 05 August 2022

Sec. Plant Metabolism and Chemodiversity

Volume 13 - 2022 | https://doi.org/10.3389/fpls.2022.945897

Association Analysis Between Genotype and Environment: Differentiation Between Cyclocarya paliurus Resources That Accumulate Triterpenoids

  • 1. College of Forestry, Nanjing Forestry University, Nanjing, China

  • 2. Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China

Article metrics

View details

8

Citations

1,8k

Views

767

Downloads

Abstract

Cyclocarya paliurus is mainly distributed in subtropical areas of China. Its leaves are rich in beneficial triterpenoids that have bioactivities against human diseases, including hyperlipemia, diabetes, and hypertension. In this study, data on the genetic diversity, distributing environment, and triterpenoids of C. paliurus samples were collected from a wide area in China. The data covered 316 C. paliurus germplasms collected from 26 distinct populations. Association analysis between genotype and triterpenoids was carried out to describe triterpenoids accumulation pattern. Based on our analyses, we identified the important trend that genotypes with higher triterpenoid contents belonged to a unique genotype subgroup. The results showed that pterocaryoside B and pterocaryoside A significantly vary among the genotypic subgroups. In addition, the different genotypic subgroups showed distinct geographical distributing areas. These findings provide information about the relationship between genetic and environmental factors and how this affects triterpenoids accumulation. This information will be valuable for targeted breeding and for further germplasm selection of C. paliurus resources.

Key Message

An association analysis based on genotypic, biochemical, geographical, and environmental data revealed that there are high-triterpenoid Cyclocarya paliurus subgroups. Each subgroup responds differently to environmental and climatic factors.

Introduction

Cyclocarya paliurus is the only extant species in the Cyclocarya genus, which is distributed only in subtropical areas of China (Manchester et al., 2009; Sun et al., 2021). The leaves are rich in triterpenoids, which can improve insulin resistance and metabolism disorder of glucose and lipid, and have the effects of reducing the lipids and pressure of blood (Kurihara et al., 2003; Xie et al., 2006). However, the contents of triterpenoids in leaves of C. paliurus vary greatly, which hindrances the development of the C. paliurus industry. Previous studies have analyzed the composition of compounds in leaves of C. paliurus, including six main triterpene compounds (e.g., arjunolic acid, cyclocaric acid B, pterocaryoside B, pterocaryoside A, hederagenin, and oleanolic acid) (Shu et al., 1995; Wright et al., 2014; Shang et al., 2015). While the triterpenoids variation patterns of different C. paliurus genotypes remain unclear.

Currently, considerable efforts are focused on the selection of C. paliurus resources according to their patterns of secondary metabolite accumulation (Cao et al., 2017; Liu et al., 2018; Sun et al., 2020a,b). Previous studies have estimated the genetic diversity of C. paliurus populations (Li X. C. et al., 2017; Qu et al., 2021a,b), the responses of secondary metabolites to cultivation measures (Deng et al., 2015, 2018; Yang et al., 2017; Liu et al., 2018; Wang et al., 2021), and secondary metabolites variations among populations (Fang et al., 2011; Cao et al., 2017, 2018). Genotyping tests of the C. paliurus germplasm have been carried out at the population level, and have shown that different genotypes respond differently to abiotic factors (Deng et al., 2015, 2017, 2019; Liu et al., 2016). In previous studies, variations in secondary metabolite accumulation could only be attributed to the specific genotype being studied, and no general patterns were identified. Cao et al. (2018) noted this and considered that there were huge variations in secondary metabolite accumulation both at the population level and at the individual level. Thus, the degree to which the genotype explains variations in secondary metabolite accumulation and the balance among different metabolites is unknown. Clarifying these relationships may be the key to understanding the strong effect of the genotype on the quality of C. paliurus resources, especially when it is hard to construct clonal propagation systems.

In this research, based on previous studies on the genotype, geographic variation, and biochemistry of C. paliurus, data from multiple C. paliurus germplasms were collected and analyzed to investigate the patterns of triterpenoids accumulation. The objectives were to identify the effects of genotype on the triterpenoids accumulation in C. paliurus. Based on our analyses, we tried to identify general patterns of triterpenoids accumulation at the genotype level. These results will be valuable for targeted breeding and for further germplasm selection.

Materials and Methods

Plant Materials

Leaves were collected from C. paliurus trees in China in 2014 (Li X. C. et al., 2017; Cao et al., 2018). The leaves were collected from 316 individual trees at 26 different sites (Table 1). The sites were distributed across 12 provinces in China, covering a range of 5.4° latitude, 12.6° longitude, and 1,370 m elevation (Table 1). Mature leaves were collected in October to ensure that leaf secondary metabolites were fully developed. The leaves were mainly from dominant or subdominant branches of trees that were >20 years old. All samples were dried, sliced, ground, and stored at room temperature until analysis.

Table 1

Population codeLocationLatitude (°N)Longitude (°E)Altitude (m)Annual mean temperature (°C)Annual mean rainfall (mm)Annual sunlight (h)Frost-free days (day)Sample size
JQ6Jixi County, Anhui30.15118.8973012.11,7261,92023311
QJ4Jingde County, Anhui30.23118.4561013.81,6571784.12403
QJ4Qimen County, Anhui30.02117.5340013.91,7751,6692325
SW5Shucheng County, Anhui31.02116.5477012.31,6061,96922412
FS13Ningbo City, Zhejiang29.76121.22715141,5271,8502328
LF14Longquan City, Zhejiang27.91119.191,20011.82,1191849.82637
LF14Wencheng County, Zhejiang27.88119.7991513.91,9131764.42855
PK11Pucheng County, Fujian27.93118.7693015.81,9981,90025410
SM12Mingxi County, Fujian26.57116.9356417.41,8211788.626112
H29Lueyang County, Shanxi33.36105.871,200137141558.32369
JY2Yongshun County, Hunan128.88110.3367014.81,4921,3062869
SS1Shaoyang City, Hunan26.37110.131,00016.71,3201348.930416
JY3Yongshun County, Hunan224.92112.0384517.41,5021,75830810
NH22Hefeng County, Hubei29.88110.421,12511.31,4991,3422456
YW21Wufeng County, Hubei30.19110.996912.61,6501,53324015
JX17Xiushui County, Jiangxi28.16114.5282713.51,6551600.424715
JJ18Jinggangshan County, Jiangxi26.51114.196713.21,8161,51124110
FD16Fenyi County, Jiangxi27.63114.5356515.31,6911,25127014
BL23Baise City, Guangxi124.46106.341,44815.61,3641906.635712
BT24Baise City, Guangxi224.61104.951,77016.71,2121569.335711
GL25Longsheng County, Guangxi25.62109.8960614.51,6291,30931411
GZ26Ziyuan County, Guangxi25.92110.3885015.51,5801,27530010
QL9Liping County, Guizhou26.34109.2472716.11,3111317.927711
QJ10Jianhe County, Guizhou26.37108.381,24015.11,2651236.330015
TS7Shiqian County, Guizhou27.35108.111,23913.31,2561232.931616
TL8Yinjiang County, Guizhou27.74108.511,03214.51,2311,2553007
M20Muchuan County, Sichuan28.97103.781,20017.31,33296833234
GQ19Guangyuan City, Sichuan32.42104.861,57014.11,0271,29224312

Collection locations and associated environmental factors.

Determination of Triterpenoids

Triterpenoids were identified and quantified in our previous studies (Cao et al., 2017, 2018). Briefly, powdered leaves (75 g from each location) were extracted with 80% ethanol in a water bath at 90°C for 2 h, then shaken for 15 min. The mixture was centrifuged at 8,000 r/min for 15 min, and a cooling centrifuge Sigma 3-18KS instrument was applied in the present study (3-18KS; Sigma, Aachen, Germany). Supernatants were concentrated at 40°C to obtain extracts. Triterpenoid contents were analyzed by high-performance liquid chromatography (HPLC) according to Cao et al. (2017, 2018), with a C18 column (X-Bridge, 250 × 4.6 mm; Waters, Milford, MA, USA) connected to an Agilent 1200 series HPLC system (Agilent, Waldbronn, Germany). A stepwise elution program was applied according to Cao et al. (2017, 2018). The characteristic common peaks of individual triterpenoids, which were sufficient to evaluate the quality of C. paliurus leaves, were identified and quantified by comparison with external standards (Liu et al., 2018). The reference standards of hederagenin, arjunolic acid, and oleanolic acid (purity >98%) were purchased from Shanghai Yuanye Biotechnology Co., Ltd. (Shanghai, China), and cyclocaric acid B, pterocaryoside A, and pterocaryoside B (purity >98%) were isolated and purified from the laboratory of China Pharmaceutical University (Nanjing, China) (Cao et al., 2017) (Figure A1; Table A1).

Environmental Factors

Geographical data (longitude, latitude, and altitude) were collected at each survey site with a GPS unit (JUNO® SCSD, Trimble, Sunnyvale, CA, USA). Climate data were obtained from ClimateAP for the historical period 1991–2014 according to Liu et al. (2018) and included annual mean temperature, annual mean rainfall, annual sunlight, and the number of frost-free days in each region (http://climateap.net/).

Extraction of Genomic DNA and PCR Amplification

The DNA extraction and PCR amplification were conducted in a previous study (Li X. C. et al., 2017) (Table 2). Briefly, DNA was isolated from dried leaves and purified using the Ezup Column Plant Genomic DNA Purification Kit (Sangon Biotech, Shanghai, China). Nine inter simple sequence repeat (ISSR) and 6 simple sequence repeat (SSR) selected primers were selected based on the reproducibility and band polymorphism (Li X. C. et al., 2017). The products were separated by 8% non-denatured polyacrylamide gel electrophoresis (PAGE) and silver-stained. The genetic similarity among populations was assessed by PCA analysis by Gen Al Ex V.6.5 Mxcomp.

Table 2

TriterpenoidsMinimum (mg/g)Maximum (mg/g)Variable coefficientMedian (Quartile) (mg/g)
Arjunolic acid0.1911.470.602.97 (1.94–4.42)
Cyclocaric acid B0.296.40.701.25 (0.77–1.94)
Pterocaryoside B019.171.081.97 (0.41–4.02)
Pterocaryoside A0.1216.410.753.04 (1.68–4.75)
Hederagenin0.024.680.710.84 (0.55–1.32)
Oleanolic acid0.071.510.620.35 (0.23–0.49)
Content of total six triterpenoids1.6844.50.5611.39 (7.51–16.24)

Triterpenoid concentrations in leaves of 316 C. paliurus individuals.

Table 3

Extract coefficientPCA1PCA2
Arjunolic acid0.690.81−0.18
Cyclocaric acid B0.720.76−0.37
Pterocaryoside B0.900.470.82
Pterocaryoside A0.820.810.40
Hederagenin0.630.69−0.39
Oleanolic acid0.680.68−0.47
Content of total six triterpenoids1.000.940.33
Total variance of interpretation (%)56.3921.24

Results of PCA analysis of triterpenoids of C. paliurus.

Table 4

Genotype groupLocation (sample size)No. of samples
1Guangyuan City, Sichuan (12)12
2Ningbo City, Zhejiang (8), Longquan City, Zhejiang (7), Pucheng County, Fujian (10), Mingxi County, Fujian (12), Xiushui County, Jiangxi (15), Jinggangshan County, Jiangxi (10), Fenyi County, Jiangxi (14)76
3Jixi County, Anhui (11), Jingde County, Anhui (3), Qimen County, Anhui (5), Shucheng County, Anhui (12), Wencheng County, Zhejiang (5), Lueyang County, Shanxi (9), Yongshun County, Hunan1 (9), Shaoyang City, Hunan (16), Yongshun County, Hunan2 (10), Hefeng County, Hubei (6), Wufeng County, Hubei (15), Baise City, Guangxi1 (12), Baise City, Guangxi2 (11), Longsheng County, Guangxi (11), Ziyuan County, Guangxi (10), Liping County, Guizhou (11), Jianhe County, Guizhou (15), Shiqian County, Guizhou (16), Yinjiang County, Guizhou (7), Muchuan County, Sichuan (34)228

C. paliurus genotype structure differentiation (Li X. C. et al., 2017).

Association Analyses

We combined the genetic diversity data with triterpenoids data and environmental data for the association analysis. The secondary metabolite data concluded triterpenoids compounds in 316 individuals. The genetic diversity data were for 26 populations at the population level. The environmental data obtained as described above (sections Environmental Factors and 2.4) were collected at the population level.

A principal component analysis (PCA) was constructed to describe the accumulation characteristics of triterpenoids. According to the genotype division, another PCA was conducted to detect differences in patterns of genotypic subgroups distributing characteristics related to environmental variables. Secondary metabolite data were marked with specific genotypic subgroups for the association analysis between genotype and secondary metabolite contents. Multiple comparisons, non-parametric tests, and PCA were conducted using SPSS software, version 19.0.0 (SPSS Inc., Chicago, IL, USA).

Results

Triterpenoid Accumulation in C. paliurus Leaves

The C. paliurus leaves were collected from 316 individual trees of 26 different population (Table 2), and the total triterpenoid content ranged from 1.68 to 44.50 mg/g (the median, 11.39 mg/g) and was the sum of six compounds: arjunolic acid (2.97 mg/g), pterocaryoside B (1.25 mg/g), pterocaryoside A (1.97 mg/g), cyclocaric acid B (3.04 mg/g), hederagenin (0.84 mg/g), and oleanolic acid (0.35 mg/g). Among them, the variable coefficient ranged from 0.56 to 1.08, cyclocaric acid B (0.70), pterocaryoside B (1.08), pterocaryoside A (0.75), and hederagenin (0.71) exhibit higher variable coefficient, while arjunolic acid (0.60), oleanolic acid (0.62), and total six triterpenoids (0.56) exhibit lower.

To explore the triterpenoids variation pattern, a principal component analysis of 7 triterpenoid indicators was constructed (Figure 1). In general, the first principal components have explained 77.6% of total variance with PCA 1 explaining 56.39% and PCA 2 explaining 21.24%. Most triterpenoid information was included, especially for total triterpenoid (1), cyclocaric acid B (0.72), pterocaryoside B (0.90), and pterocaryoside A (0.82) (Table 1). As shown in Figure 1, different types of triterpenoids showed almost the same variation direction, indicating that they changed in basically the same pattern.

Figure 1

Genetic Structure of Natural Populations

In a previous study on the genetic diversity of C. paliurus resources, no significant differences in coefficients were detected in a UPGMA (Li X. C. et al., 2017). Despite long-term geographical isolation, the 26 populations separated into five groups at the 0.91 level in the ISSR and SSR marker-based UPGMA dendrogram: GQ19, TS7, and SS1 each formed separate groups (Li X. C. et al., 2017); group 1 consisted of PK11, SM12, LF14, FS13, FD16, JX17, and JJ18; and group 2 consisted of JY2, M20, H29, YW21, NH22, QJ4, SW5, JQ6, QL9, TL8, QJ10, JY3, BL23, BT24, GL25, and GZ26. However, the UPGMA method seemed to produce ambiguous results that were rather limited for association analyses of population division.

In this study, instead of conducting analyses at the population or individual levels as was done in previous studies (Deng et al., 2015, 2017, 2019; Liu et al., 2016, 2018; Cao et al., 2017, 2018; Zhou et al., 2019), we conducted a PCA to analyze the phylogenetic relationships (Li X. C. et al., 2017) in the C. paliurus population (Figure 2). The 26 populations were divided into three clear subgroups according to the first two components. Subgroup I included population GQ19, subgroup II included populations JJ18, FS13, SM12, JX17, PK11, FD16, and LF14, and subgroup III contained the populations NH22, M20, QL9, JQ6, YW21, SW5, JY2, QJ10, TS7, JY3, GL25, QJ4, TL8, BL23, SS1, GZ26, H29, and BT24. This partly agreed with the results of the previous UPGMA clustering analysis (Li X. C. et al., 2017), but the groups were more clearly defined for further analyses.

Figure 2

According to the PCA results, 26 different C. paliurus population has been divided into three genotype subgroups. The first genotype subgroup included 12 individuals of the C. paliurus population located in Guangyuan City, Sichuan. And the second genotype subgroup included 76 individuals of 7 C. paliurus populations from Zhejiang province, Fujian province, and Jiangxi province. The third genotype subgroup included 228 individuals of 20 C. paliurus populations from Anhui province, Zhejiang province, Shanxi province, Hunan province, Hubei province, Guangxi province, Guizhou province, and Sichuan province. According to the clustering result, most germplasm resources belong to genotype subgroup 3, and subgroup 1 only contained 1 population.

Genotype and Environment Association Analysis

To determine whether the different genotype subgroups of C. paliurus showed different responses to environmental factors, a PCA was carried out to analyze the distributing environmental data of different genotype C. paliurus subgroups (Figure 3). The first two PCAs (Table 5), analyzed 7 environmental variables, namely, latitude, longitude, altitude, annual mean temperature, annual mean rainfall, annual sunlight, and frost-free days, explained 74.95% of the total variation (PC 1 explained 44.95% and PC 2 explained 30%). Most environmental information was included, especially for latitude (0.83), longitude (0.95), annual mean rainfall (0.80), and frost-free days (0.90). According to the results, the three genotype subgroups showed a significant difference in geographical distributing patterns. As Figure 3 showed, individuals from genotype subgroup 2 were almost all distributed in the first quartile, with higher longitude, annual mean rainfall, and annual sunlight. While individuals from genotype subgroup 3 were almost all distributed in the second and fourth quartile, with higher latitude, annual mean temperature, and frost-free days. And genotype subgroup 1 was distributed only in the third quartile with higher altitudes.

Figure 3

Table 5

Extract coefficientPCA1PCA2
Latitude0.830.34−0.85
Longitude0.950.930.28
Altitude0.58−0.67−0.35
Annual mean temperature0.70−0.620.56
Annual precipitation0.800.720.54
Annual sunlight0.500.690.15
Annual frost-free days0.90−0.860.41
Total variance of interpretation (%)44.9530.00

Results of PCA analysis of environmental factors of C. paliurus.

To further analyze the geographical variation pattern, a Kruskal–Wallis test was made to find the environment factors with significant differences between the three genotype subgroups (Table 6). According to the result, it was found that the geographical distribution of different genotype subgroups exhibits significant difference along with longitude (P < 0.05) and annual mean rainfall (P < 0.01). Multiple comparison results also showed that genotype subgroup 2 was almost distributed in the environment with the median longitude being 114.53° (quartile 114.52°-118.76°), and genotype subgroup 3 almost distributed at lower longitude being 109.24° (quartile 106.34°-110.9°). Besides, genotype subgroup 2 is almost all distributed in the environment with the median annual mean rainfall being 1,816 mm (quartile 1,655–1,821 mm), much higher than genotype subgroup 3 with the median being 1,332 mm and the quartile being 1,265–1,606 mm. In general, environment gradient, including longitude and annual mean rainfall, has been certified significantly correlated with the genotype geographical differentiation of C. paliurus germplasm.

Table 6

Environment factorsSignificanceGenotype 1Genotype 2Genotype 3
Latitude0.347
Longitude0.011104.86Aa114.53 (114.52–118.76)Cc109.24 (106.34–110.9)Bb
Altitude0.189
Annual mean temperature0.997
Annual mean rainfall0.0011027Aa1,816 (1,655–1,821)Cc1,332 (1,265–1,606)Bb
Annual sunlight0.429
Frost-free days0.548

Kruskal–Wallis test and multiple comparisons between genotypes and environmental factors of C. paliurus.

Lower case letters indicate that the environmental factors have a significant difference at the P < 0.05 level in genotypes; upper case letters indicate that different genotypes have a significant difference at the P < 0.01 level.

Genotype Subgroups Associated With Particular Patterns of Triterpenoids Accumulation

According to the secondary triterpenoids variation pattern of C. paliurus individuals, no significant trend was explored. To further study, the triterpenoids PCA plot was combined with genotypic subgroup data, and it was found that almost all the high-triterpenoids type germplasms were in genotype subgroup 2, but the genotype subgroup 1 showed much lower content of triterpenoids than the others (Figure 4). Despite most individuals, more samples in subgroup 3 showed relatively lower high triterpenoids content than genotype subgroup 2. And this may imply that genotype subgroup clustering influenced C. paliurus leaf triterpenoids accumulation to some extent.

Figure 4

To further analyze the association between genotype subgroups and triterpenoids, a Kruskal–Wallis test was made between seven triterpenoid indicators and two genotype subgroups (Table 7), and genotype subgroup 1 was not evaluated because of the small sample size. According to the Kruskal–Wallis test result, it was found that five triterpenoid indicators exhibit a significant difference between genotype group 2 and genotype group 3 at P < 0.01 level. Content of a total of six triterpenoids, pterocaryoside B, pterocaryoside A, hederagenin, and oleanolic acid all showed a significant difference between genotype groups, while arjunolic acid and cyclocaric acid B were not influenced between genotype subgroup 2 and genotype subgroup 3.

Table 7

TriterpenoidsSignificanceGenotype 1Genotype 2Genotype 3
Arjunolic acid0.071.1 (1–1.34)Aa3.2 (2.16–5.35)Bb2.98 (2–4.23)Bb
Cyclocaric acid B0.080.41 (0.3–0.51)Aa1.49 (0.88–2.14)Bb1.24 (0.79–1.86)Bb
Pterocaryoside B00.44 (0.25–0.7)Aa3.35 (1.78–5.7)Cc1.68 (0.19–3.51)Bb
Pterocaryoside A00.57 (0.25–0.73)Aa5.01 (3.42–7.44)Cc2.65 (1.51–4.06)Bb
Hederagenin01 (0.7–1.36)ABab1.18 (0.73–1.63)Bb0.73 (0.49–1.13)Aa
Oleanolic acid00.64 (0.51–0.68)Bb0.44 (0.34–0.62)Bb0.3 (0.2–0.42)Aa
Content of total six triterpenoids04.18 (3.63–4.81)Aa15 (11.27–22.24)Cc10.34 (6.97–14.99)Bb

Kruskal–Wallis test and multiple comparisons between genotypes and triterpenoids of C. paliurus.

Lower case letters indicate that the triterpenoid individual has a significant difference at the P < 0.05 level in genotypes; upper case letters indicate that the triterpenoid individuals of different genotypes have a significant difference at the P < 0.01 level.

For further analysis, multiple comparisons were carried out to deeply identify the triterpenoids variation pattern (Table 7). Genotype subgroup 2 has shown the highest total triterpenoid content with the median being 15 mg/g and the quartile being 11.27–22.24 mg/g, significantly higher (P < 0.01) than genotype subgroup 3 with the median being 10.34 mg/g and the quartile being 6.97–14.99 mg/g. And genotype subgroup 1 has shown the lowest total triterpenoid content with the median being only 4.18 mg/g and the quartile being 3.63–4.81 mg/g. Besides, pterocaryoside B (3.35, 1.78–5.7 mg/g) and pterocaryoside A (5.01, 3.42–7.44 mg/g) of genotype subgroup 2 showed higher content than genotype subgroup 3 in P < 0.01 level (pterocaryoside B, 1.68 mg/g, 0.19–3.51 mg/g; pterocaryoside A, 2.65 mg/g, 1.51–4.06 mg/g). And genotype subgroup 1 has also showed the lowest pterocaryoside B (0.44, 0.25–0.7 mg/g) and pterocaryoside A (0.57, 0.25–0.73 mg/g). For hederagenin, genotype subgroup 2 showed significant higher content (1.18, 0.73–1.63 mg/g) than genotype subgroup 3 (0.73, 0.49–1.13 mg/g) at P < 0.01 level. And hederagenin content of genotype subgroup 1 showed no significant difference from the other two. To the arjunolic acid content, genotype subgroup 2 (3.2, 2.16–5.35 mg/g) and genotype subgroup 3 (2.98, 2–4.23 mg/g) showed no significant difference, while significant higher than genotype subgroup 1 (0.41, 0.3–0.51 mg/g) in P < 0.01 level. In general, most triterpenoid individuals, including total triterpenoid content, have been certified significantly correlated with the genotype differentiation of C. paliurus germplasm.

Discussion

In a previous study on the genetic diversity of C. paliurus resources, no significant differences in coefficients were detected in a UPGMA (Li X. C. et al., 2017). In that study, the 26 C. paliurus populations were divided into subgroups based on genetic diversity detected using nine ISSR primers and six SSR primers (Li X. C. et al., 2017). Despite long-term geographical isolation, the 26 populations were separated into five groups at the 0.91 level in the ISSR and SSR marker-based UPGMA dendrogram. However, the UPGMA method seemed to produce ambiguous results that were rather limited for association analyses of population division. In this study, instead of conducting analyses at the population or individual levels as was done in previous studies (Cao et al., 2018; Liu et al., 2018; Deng et al., 2019; Zhou et al., 2019), we conducted a PCA to analyze the phylogenetic relationships in the C. paliurus population. The 26 populations were divided into three clear subgroups according to the first two components, this partly agreed with the results of the previous UPGMA clustering analysis (Li X. C. et al., 2017), but the groups were more clearly defined for further analyses.

The genotype × environment interaction is considered to play an important role in secondary metabolite accumulation (Deng et al., 2019). In this research, it was found that genotype subgroups influenced secondary metabolite accumulation, consistent with the results of Deng et al. (2015, 2017). Deepak et al. (2018) found that genotype had a stronger effect than provenance on variations in secondary metabolite accumulation in leaves of silver birch. Genotypes also have been shown to differ in their responses to environmental factors (Keinänen et al., 1999; Laitinen et al., 2005). Accordingly, the results of this research showed that genotype generally had a stronger effect than environmental factors on triterpenoids differentiation in leaves of C. paliurus.

In a previous study, environmental factors have been considered to be the main reason for influencing plant triterpenoids accumulation. Yeh et al. (2016) found that nitrogen significantly promoted triterpenoid production in an endophytic fungus. Park et al. (2005) and Schlag and McIntosh (2006) found that the triterpenoid contents and types differed among the same plant species grown in different environments. Sun et al. (2018, 2019) and Fan et al. (2013) found that Sapindus mukorossi growing at sites with higher temperatures had higher contents of the triterpenoid saponin in fruits. While fewer researches have focused on triterpenoid differentiation between genotypes. Deng et al. (2018) concluded that the environment had a stronger influence than the genotype on total triterpenoid content, but the genotype had a stronger influence than the environment on the accumulation of individual triterpenoids.

Our previous research has considered that altitude and longitude were significantly correlated with C. paliurus triterpenoid (Sun et al., 2021), while the conclusion did not take the genotype influence into account. In this research, according to the statistical results, total triterpenoids have shown a significant difference between genotypes, which can improve insulin resistance, improve glucolipid metabolism disorder, and have the effects of lowering blood sugar and blood lipid. Besides, the response of different triterpenoid individuals to genotype differences was not consistent. Tetracyclic triterpenes and pentacyclic triterpenes have the same synthetic precursors, but correspond to different synthetic paths (Deng et al., 2017). Therefore, the response of these two different types of triterpenoids to genotype differentiation may also be different, which is agreed with Deng et al. (2018). In this research, arjunolic acid, cyclocaric acid B, hederagenin, and oleanolic acid were pentacyclic triterpenoid acids, while pterocaryoside B and pterocaryoside A were tetracyclic triterpenoid saponins, which exhibit a significant difference between genotype subgroups. Pterocaryoside B and pterocaryoside A are both specific to the species and are highly sweet secodammarane glycosides (Kennelly et al., 1995; Kinghorn and Soejarto, 2002). The genotype difference in pterocaryoside B and pterocaryoside A accumulation may lead to different sweetness and other potential bioactivities.

More previous studies have investigated genotypic variation at population or family levels, and the conclusions are drawn only apply to those specific populations or families. However, in this study, we detected a general pattern of triterpenoids accumulation among C. paliurus resources and found that triterpenoid accumulation responses differed among genotype subgroups. According to the genotype subgroups defined based on molecular marker data, the C. paliurus germplasm resources could be divided into a high-triterpenoid group and others. Previous studies have suggested that environmental factors have stronger effects than the genetic background on the accumulation of triterpenoids (Deng et al., 2018) in the leaves of C. paliurus. However, the wide variations in secondary metabolite contents detected in this study implied that genotype may also be an important factor. Therefore, the extent to which the variations among populations were influenced by genotype and environmental factors was evaluated.

Conclusion

In this study, three genotype subgroups exhibit significant distributing characteristics to the environmental gradient. Different genotype subgroups of C. paliurus showed different patterns of total triterpenoids accumulation. The germplasm could be divided into the high-triterpenoids group and normal groups, which will be useful for further selection. Among triterpenoids, different triterpenoid individuals' variation pattern to genotype differences was not consistent, pterocaryoside B and pterocaryoside A significantly accumulated in the high-triterpenoids genotype group. This information will be useful for the targeted breeding of C. paliurus for particular purposes, and for optimizing cultivation conditions to favor the accumulation of certain groups of metabolites. More secondary metabolites should be considered in further association analysis research to better understand the effect of the genotype × environment interaction on secondary metabolite accumulation in C. paliurus leaves.

Funding

This work was supported by the National Natural Science Foundation of China (31971642), the Key Research and Development Program of Jiangsu Province (BE2019388).

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

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

CS and SF: methodology, experimental design, and writing—review and editing. WY, HD, YC, and XL: data curation. CS: data analysis and writing—original draft preparation. XS: funding acquisition. All authors have read and agreed to the published version of the manuscript.

Conflict of interest

The 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.

    Abbreviations

  • ISSR

    intersimple sequence repeat

  • PAGE

    polyacrylamide gel electrophoresis

  • PCA

    principal component analysis

  • PCR

    polymerase chain reaction

  • SSR

    simple sequence repeat

  • UPGMA

    unweighted pair group method with arithmetic mean.

References

  • 1

    CaoY. N.DengB.FangS. Z.ShangX. L.FuX. X.YangW. X. (2018). Genotypic variation in tree growth and selected flavonoids in leaves of Cyclocarya paliurus. Molecules80, 6774. 10.2989/20702620.2016.1274862

  • 2

    CaoY. N.FangS. Z.YinZ. Q.FuX. X.ShangX. L.YangW. X.et al. (2017). Chemical fingerprint and multicomponent quantitative analysis for the quality evaluation of Cyclocarya paliurus leaves by HPLC-Q-TQF-MS. Molecules22, 19271942. 10.3390/molecules22111927

  • 3

    DeepakM.LihavainenJ.Keski-SaariS.Kontunen-SoppelaS.SalojarviJ.TenkanenA.et al. (2018). Genotype- and provenance-related variation in the leaf surface secondary metabolites of silver birch. Can. J. Forest Res. 48, 494505. 10.1139/cjfr-2017-0456

  • 4

    DengB.CaoY. N.FangS. Z.ShangX. L.YangW. X.QianC. Y. (2015). Variation and stability of growth and leaf flavonoid content in Cyclocarya paliurus across environments. Ind. Crop Prod. 76, 386393. 10.1016/j.indcrop.2015.07.011

  • 5

    DengB.FangS.ShangX.FuX.LiY. (2017). Influence of provenance and shade on biomass production and triterpenoid accumulation in Cyclocarya paliurus. Agroforest Syst. 1, 110. 10.1007/s10457-017-0138-x

  • 6

    DengB.FangS.ShangX.FuX.YangW. (2018). Influence of genotypes and environmental factors on leaf triterpenoid content and growth of Cyclocarya paliurus. J. For. Res. 4, 110.

  • 7

    DengB.FangS. Z.ShangX. L.FuX. X.YangW. X. (2019). Influence of genotypes and environmental factors on leaf triterpenoid content and growth of Cyclocarya paliurus. J. For. Res.30, 789798. 10.1007/s11676-018-0680-z

  • 8

    FanH. H.YaoX. M.ZhangT. Y.LinX. Q.TangX. H.MaL. (2013). The selection of high-saponin-yield plus tree of Sapindus mukorossi. J. For. Sci. Technol. 40, 6973. 10.3969/j.issn.1002-7351.2013.04.14

  • 9

    FangS. Z.YangW. X.ChuX. L.ShangX. L.SheC. Q.FuX. X. (2011). Provenance and temporal variations in selected flavonoids in leaves of Cyclocarya paliurus. Food Chem. 124, 13821386. 10.1016/j.foodchem.2010.07.095

  • 10

    KeinänenM.Julkunen-TiittoR.MutikainenP.WallsM.OvaskaJ.VapaavuoriE. (1999). Trade-offs in phenolic metabolism of silver birch: effects of fertilization, defoliation, and genotype. Ecology80, 19701986. 10.1890/0012-9658(1999)080[1970:TOIPMO]2.0.CO;2

  • 11

    KennellyE.CaiL.LongL.ShamonL.ZawK.ZhouB.et al. (1995). Novel Highly Sweet secodammarane glycosides from Pterocarya paliurus. J. Agric. Food Chem.43, 26022607. 10.1021/jf00058a009

  • 12

    KinghornA. D.SoejartoD. D. (2002). Discovery of terpenoid and phenolic sweeteners from plants. Pure Appl. Chem.74, 11691179. 10.1351/pac200274071169

  • 13

    KuriharaH.AsamiS.ShibataH.FukamiH.TanakaT. (2003). Hypolipemic effect of Cyclocarya paliurus (Batal) Iljinskaja in lipid-loaded mice. Biol. Pharm. Bull. 26, 383385. 10.1248/bpb.26.383

  • 14

    LaitinenM. L.Julkunen-TiittoR.TahvanainenJ.HeinonenJ.RousiM. (2005). Variation in birch (Betula pendula) shoot secondary chemistry due to genotype, environment, and ontogeny. J. Chem. Ecol. 31, 697717. 10.1007/s10886-005-3539-7

  • 15

    LiX. C.FuX. X.ShangX. L.YangW. X.FangS. Z. (2017). Natural population structure and genetic differentiation for heterodicogamous plant: Cyclocarya paliurus (Batal.) Iljinskaja (Juglandaceae). Tree Genet. Genomes13, 80. 10.1007/s11295-017-1157-5

  • 16

    LiuY.CaoY. N.FangS. Z.WangT. L.YinZ. Q.ShangX. L.et al. (2018). Antidiabetic effect of Cyclocarya paliurus leaves depends on the contents of antihyperglycemic flavonoids and antihyperlipidemic triterpenoids. Molecules23, 10421058. 10.3390/molecules23051042

  • 17

    LiuY.QianC. Y.DingS. H.ShangX. L.YangW. X.FangS. Z. (2016). Effect of light regime and provenance on leaf characteristics, growth and flavonoid accumulation in Cyclocarya paliurus (Batal) Iljinskaja coppices. Bot. Stud. 57, 28. 10.1186/s40529-016-0145-7

  • 18

    ManchesterS. R.ChenZ. D.LuA. M.UemuraK. (2009). Eastern Asian endemic seed plant genera and their paleogeographic history throughout the Northern Hemisphere. J. Syst. Evol.47, 142. 10.1111/j.1759-6831.2009.00001.x

  • 19

    ParkJ. D.RheeD. K.LeeY. H. (2005). Biological activities and chemistry of saponins from Panax ginseng C. A. Meyer. Phytochem. Rev. 4, 159175. 10.1007/s11101-005-2835-8

  • 20

    QuY. Q.KongW. L.WangQ.FuX. X. (2021a). Genome-wide identification MIKC-Type MADS-box gene family and their roles during development of floral buds in wheel wingnut (Cyclocarya paliurus). Int. J. Mol. Sci. 22, 114. 10.3390/ijms221810128

  • 21

    QuY. Q.WangQ.YuY. H.FuX. X. (2021b). Comparative analysis of the complete chloroplast genome between tetraploidy and diploidy of Cyclocarya paliurus (Batal.) Iljinskaja. Mitochondrial DNA Part B Resourc. 6, 26692671. 10.1080/23802359.2021.1964402

  • 22

    SchlagE. M.McIntoshM. S. (2006). Ginsenoside content and variation among and within American ginseng (Panax quinquefolius L.) populations. Phytochemistry67, 15101519. 10.1016/j.phytochem.2006.05.028

  • 23

    ShangX. L.WuZ. F.YinZ. Q.ZhangJ.LiuZ. J.FangS. Z. (2015). Simultaneous determination of flavonoids and triterpenoids in Cyclocarya paliurus leaves using high-performance liquid chromatography. Afr. J. Tradit. Complem.12, 125134. 10.4314/ajtcam.v12i3.16

  • 24

    ShuR. G.XuC. R.LiL. N. (1995). Studies on the sweet principals from the leaves of Cyclocarya paliurus (Batal.) Iljinsk. Acta Pharmaceut. 30, 757761.

  • 25

    SunC.FangS.ShangX. (2020a). Triterpenoids biosynthesis regulation for leaf coloring of wheel wingnut (Cyclocarya paliurus). Forests12, 111. 10.3390/f12121733

  • 26

    SunC.ShangX.DingH.CaoY.FangS. (2020b). Natural variations in flavonoids and triterpenoids of Cyclocarya paliurus leaves. J. For. Res.31, 805814. 10.1007/s11676-020-01139-1

  • 27

    SunC.ZhouY.FangS.ShangX. (2021). Ecological gradient analysis and environmental interpretation of Cyclocarya paliurus communities. Forests12, 146. 10.3390/f12020146

  • 28

    SunC. W.JiaL. M.XiB. Y.LiuJ. M.WangL.WengX. H. (2019). Genetic diversity and association analyses of fruit traits with microsatellite ISSRs in Sapindus. J. For. Res.30, 193203. 10.1007/s11676-017-0580-7

  • 29

    SunC. W.WangL. C.LiuJ. M.ZhaoG. C.GaoS. L.XiB. Y.et al. (2018). Genetic structure and biogeographic divergence among Sapindus species: an inter-simple sequence repeat-based study of germplasms in China. Ind. Crop. Prod. 118, 110. 10.1016/j.indcrop.2018.03.029

  • 30

    WangZ. K.XuZ. H.ChenZ. Y.KowalchukG. A.FuX. X.KuramaeE. E. (2021). Microbial inoculants modulate growth traits, nutrients acquisition and bioactive compounds accumulation of Cyclocarya paliurus (Batal.) Iljinskaja under degraded field condition. For. Ecol. Manage. 482, 111. 10.1016/j.foreco.2020.118897

  • 31

    WrightM.ByrdJ.GaoY.StubblefieldJ.ParkH.DunlapN. (2014). Isolation and structural clarification of triterpenes from Cyclocarya paliurus: cyclocaric acid A and B. Planta Med. 80, 139143. 10.1055/s-0034-1382440

  • 32

    XieM. Y.LiL.NieS. P.WangX. R.LeeF. S. C. (2006). Determination of speciation of elements related to blood sugar in bioactive extracts from Cyclocarya paliurus leaves by FIA-ICP-MS. Eur. Food Res. Technol. 223, 202209. 10.1007/s00217-005-0173-0

  • 33

    YangW.LiuY.FangS.DingH.ZhouM.ShangX. (2017). Variation in growth, photosynthesis and water-soluble polysaccharide of Cyclocarya paliurus under different light regimes. IForest10, 468474. 10.3832/ifor2185-010

  • 34

    YehC. W.KanS. C.LinC. C.ShiehC. J.LiuY. C. (2016). Polyhydroxylated steroids and triterpenoids from an entophytic fungus, Hypocreales sp. NCHU01 isolated from Tuber magnatum. J. Taiwan Inst. Chem. E64, 2230. 10.1016/j.jtice.2016.03.049

  • 35

    ZhouM. M.LinY.FangS. Z.LiuY.ShangX. L. (2019). Phytochemical content and antioxidant activity in aqueous extracts of Cyclocarya paliurus leaves collected from different populations. Peer J.7, e6492. 10.7717/peerj.6492

Appendix

Figure A1

Table A1

Peak No.tR (min)[M-H]-MS/MS fragment ion (m/z)FormulaIdentification
1247.1487.343445.2942, 401.3056, 389.2698C30H48O5Arjunolic acid
1349.7485.328/C30H46O5Cyclocaric acid B
1454.2621.400521.3107, 489.3571C35H58O9Pterocaryoside B
1560.0635.416535.3265, 489.3573C36H60O9Pterocaryoside A
1660.5471.348145.0286C30H48O4Hederagenin
2083.5455.355/C30H48O3Oleanolic acid

Identification of 6 triterpenoid compounds from leaves of C. paliurus by developed HPLC-Q-TOF-MS.

Summary

Keywords

triterpenoids, genotype, genetic resource, Cyclocarya paliurus, association analysis

Citation

Sun C, Shang X, Fang S, Yang W, Cao Y, Ding H and Li X (2022) Association Analysis Between Genotype and Environment: Differentiation Between Cyclocarya paliurus Resources That Accumulate Triterpenoids. Front. Plant Sci. 13:945897. doi: 10.3389/fpls.2022.945897

Received

17 May 2022

Accepted

20 June 2022

Published

05 August 2022

Volume

13 - 2022

Edited by

Luigi Lucini, Catholic University of the Sacred Heart, Italy

Reviewed by

Hiroshi Noguchi, Nihon Pharmaceutical University, Japan; Leilei Zhang, Catholic University of the Sacred Heart, Italy

Updates

Copyright

*Correspondence: Xulan Shang Shengzuo Fang ;

This article was submitted to Plant Metabolism and Chemodiversity, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics