Next Article in Journal
Study on Farmers’ Willingness to Maintain the Sloping Land Conversion Program in Ethnic Minority Areas under the Background of Subsidy Expiration
Previous Article in Journal
Bat Assemblages Are Shaped by Land Cover Types and Forest Age: A Case Study from Eastern Ukraine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phenotypic Variation and Diversity in Fruit, Leaf, Fatty Acid, and Their Relationships to Geoclimatic Factors in Seven Natural Populations of Malania oleifera Chun et S.K. Lee

1
Experimental Center of Tropical Forestry, Chinese Academy of Forestry, Pingxiang 532600, China
2
Beijing Floriculture Engineering Technology Research Centre, Beijing Botanical Garden, Beijing 100093, China
3
Guangxi Institute of Chinese Medicine & Pharmaceutical Science, Nanning 530022, China
4
State Key Laboratory of Tree Genetics and Breeding, Key Laboratory of Tree Breeding and Cultivation of State Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Forests 2022, 13(10), 1733; https://doi.org/10.3390/f13101733
Submission received: 13 September 2022 / Revised: 13 October 2022 / Accepted: 17 October 2022 / Published: 20 October 2022
(This article belongs to the Section Forest Biodiversity)

Abstract

:
Malania oleifera Chun et S.K. Lee is an endangered woody oil tree that produces nervonic acid, which is closely related to brain development. The phenotypic variation and diversity in some representative traits in natural populations are poorly understood. To investigate phenotypic variation, diversity and their relationship with geoclimatic factors, we evaluated 21 traits of 97 sampled individuals from seven natural populations. Coefficient of variation for 21 traits widely ranged from 0.05 (shape index of fruit) to 0.62 (water content of dry kernel), with an average of 0.18. The Shannon–Wiener indices ranged from 1.91 (water content of dry kernel) to 2.06 (leaf length and leaf width), with an average of 1.94. The variation of 20 traits was significantly influenced by the population effect (p < 0.05). The 21 traits were reduced to six principal components, which collectively accounted for 86.64% of the total phenotypic variation. The seven population was mainly divided into two groups: Group I, which included four populations distributed throughout the southern region and having wider and shorter leaves, larger fruit transverse diameters, and higher content of fatty acid and nervonic acid in the kernel, and Group II included the left populations with relative opposite characteristics to those in Group I. Trees with shorter leaves tend to produce oblate shape and heavier dried kernel, yielding more fatty acid and nervonic acids, and population-level pairwise kernel morphological distance were positively correlated with fatty acid distance (Mantel r = 0.39, p < 0.05). Annual average humidity was strongly positively correlated with leaf index (length/width, r = 0.78) but had a negative correlation with the kernel dry weight (r = −0.84). The shape of the kernel gradually altered from oblate to globose as longitude and rainfall increased. The environment became colder as latitude increased, resulting in a steady decrease in fatty acid content. No significant isolation by climatic and geographical distance was found for morphological and fatty acid variability. These results can have an important implication for genetic improvement, diversity conservation and resource management of the species in the future.

1. Introduction

Malania oleifera Chun et S.K. Lee, a monotypic species within the Olacaceae family, is an ancient evergreen tree species endemic to SW China (at the junction of Yunnan, Guangxi, and Guizhou provinces) [1]. It is naturally distributed with an extremely small population in a narrow, rocky desertification area and listed as Vulnerable (VU) on the IUCN Red List [2,3] and a Class II key protected species in China. M. oleifera is a famous woody oil plant for its kernels containing high oil content (60.5%), of which nervonic acid (C24H46O2, hereafter referred to as NAC) was the highest, accounting for 55.7–67.10% [4,5]. NAC is a monounsaturated fatty acid and has attracted considerable attention for its close relationship with brain development. There are merely about ten woody plants in China containing NAC content of more than 2% in kernels, with M. oleifera yielding the highest NAC and having the greatest commercial potential [6]. Given its significant economic value, so far, extensive effort in M. oleifera has been devoted to examining the content and composition of fatty acids [5,7], identifying the genes related to NAC biosynthesis through molecular biology methods [8], transcriptome analysis [9], genome sequencing [10], etc. In addition, genetic improvement of key traits is a conventional and promising approach in crops and trees. Furthermore, common gardens and trials established in a genetic improvement program are essential for ex situ conservation. However, genetic improvement of the levels of NAC in M. oleifera kernels is still in its infancy.
Characterization of the diversity and variation of valuable traits is crucial for genetic improvement and reasonable conservation. Generally, the diversity of germplasm resources is the expression of inherent genetic information in a variety of environments, reflecting the richness of genes within a species and its response to environmental changes [11]. In natural populations, forest trees generally harbor great diversity and variation due to inherent genetic-based differences, various environmental conditions and their interactions in adaptation across a range of habitats [12]. The inter- and intra-population diversity and variation have been extensively reported in tree improvement studies [13,14,15,16,17,18]. In most cases, geoclimatic factors of different populations significantly influenced the diversity level and have gained much attention. For instance, the phenotypic variation of leaves was influenced by climatic factors such as temperature and rainfall and changed along gradients of geographical factors like longitude, latitude, and altitude in Tetracentron sinense Oliv. [19]. Seed morphological characteristics such as mass and size were altered along latitudinal and longitudinal gradients in Euptelea pleiosperma J.D. Hooker et Thomson [17]. The fatty acid profile of Xanthoceras sorbifolium Bunge was affected by latitude, mean annual temperature and longitude [18]. In M. oleifera, Wang et al. investigated the morphological characterization of fruits and the concentration of seed oil and its fatty acid composition using fifty-three trees from two populations [20]. However, more comprehensive information on phenotypic variation, diversity and their relationship with geoclimatic factors at the population level is needed.
In this study, we sampled 97 individuals from seven natural populations, based on a survey of the main distribution regions of M. oleifera in China, and measured 21 phenotypic traits of fruits, kernels, leaves, fatty acid and nervonic acid content. Our main objectives were to (1) quantify the phenotypic variation and diversity; (2) test the relationships among phenotype traits to find if there were indicators of fatty acid and nervonic acid content; (3) examine the influence of geoclimatic factors on the phenotypic variability. The results will provide the foundations for subsequent genetic improvement, diversity conservation, afforestation and resource management of the species.

2. Materials and Methods

2.1. Plant Materials

From September to October 2021, we surveyed the primary distribution regions of M. oleifera. The mature fruits, branches, and leaves were collected from each sampled tree and brought back for measurement. The sample sizes differed due to the population size, the fruit availability and the sampling difficulty. A total of 97 individuals were sampled from seven populations, including Bama, Fengshan, Leye, Lingyun, Tianlin in Guangxi province, and Funing, Guangnan in Yunnan province (Figure 1). The altitude, latitude, and longitude of each population were documented by the global positioning systems (GPS). Annual mean temperature, annual mean rainfall, and relative humidity of each population were obtained from the China Meteorological Data Network (http://data.cma.cn/; 1981–2010, accessed on 1 April 2022). Basic information on the populations, including sample sizes and ages, as well as climate and geographic factors, is presented in Table 1.

2.2. Trait Measurements

Morphological traits of the fruit, kernel, and leaf were measured with calipers (in mm). For each tree, twelve fruits were measured and then the measurements were averaged for each trait. The thickness of the pericarp is the difference between the transverse diameter of the fruit and kernel. The shape index of the fruit or kernel was calculated as the ratio of the vertical diameter to the transverse diameter of the fruit or kernel. The weight of the fruit and the kernel was measured with a portable electronic balance (in g). The water content of fresh kernels was determined after drying in an oven at 60 °C for 72 h to constant weight. The length, width, area, and perimeter of leaves were measured using the WSeen LA-S System (WSeen Testing Technology Co., Ltd., Hangzhou, China). For each tree, 20 leaves were measured and then the measurements were averaged for each trait. The leaf index was the length-to-width ratio of a leaf.
Quality traits of dried kernels consist of water, fatty, and nervonic acid content. The dried kernels were ground to powder in liquid nitrogen and subsequently analyzed by the Agricultural Products Quality Supervision, Inspection and Testing Center (Zhengzhou) of the Ministry of Agriculture and Rural Affairs in accordance with the following relevant national technical standards. The water content was determined by the direct drying method following the National food safety standard–Determination of moisture in foods (GB5009.3–2016). The determination of fatty acid and nervonic acid content was carried out using the A1-3 6890N gas chromatograph by the alkaline hydrolysis method following the National food safety standard–Determination of fatty acids in foods (GB5009.168–2016).

2.3. Statistical Analysis

All statistical analyses were performed using R software (version 4.0.2) [21]. The coefficients of phenotypic variation (CV) indicate the magnitude of variation and were calculated as
C V = σ / x
where σ and x are standard deviation and arithmetic mean.
The Shannon–Wiener index (H) was used to evaluate the population diversity for each trait. The phenotypic trait was converted to a ten-grade variable: the values < x 2 σ were assigned to the first grade; the remaining grades were set in steps of 0.5 σ [22].
H = p i ( l n p i )
where pi indicates the percentage of each grade.
Variance analysis combined with F-test was used to test the effects of population on each trait and was implemented using the R function ‘aov’. Subsequently, multiple comparisons and the statistical significance at p < 0.05 were performed using the ‘duncan.test’ function in the R package agricolae (version 1.3-5) [23].
Principal component analysis (PCA) was used to reduce the variability and explore the trait structure and relationships among populations. Data standardization using the Z-score method was carried out prior to PCA. PCA was performed using the R function ‘prcompt’.
To investigate the relationships among phenotype traits and their correlations with geoclimatic factors, firstly, Pearson’s correlation coefficients and the statistical significance were obtained using the R function ‘cor.test’. Next, the Mantel test [24] was carried out to test the correlations between fruit, kernel, and leaf morphology, chemical traits, as well as their variation patterns of isolation by isolation by climatic and geographical distance based on six different matrices of distances among the studied populations. The first five matrices were obtained using the ‘vegdist’ function in the R package vegan (version 2.6-2) [25]. Squared Mahalanobis distances between the studied populations were calculated to obtain the matrices of morphometric (fruit, kernel and leaf) and chemical distances. Climatic differences were calculated as the Euclidian distance for three climatic factors of the site of sample collection (Table 1). Geographic distances in the last matrix were calculated from the latitude and longitude of the site of sample collection using the ‘distm’ function with the Haversine method in the R package geosphere (version 1.5-14) [26]. The Mantel test was implemented in the R package vegan (version 2.6-2) [25] using the ‘mantel’ function. The Canberra distance, a weighted version of the Manhattan distance, was produced by the R function ‘dist’ and then was subjected to a clustering procedure (UPGMA method) using the ‘upgma’ function in the R package phangorn (version 2.7.0) [27].

3. Results

3.1. Phenotypic Variation and Diversity

The phenotypic variation for 21 traits of M. oleifera is shown in Table 2. CVs of the 21 traits ranged from 0.05 (shape index of fruit, FSI) to 0.62 (water content of dry kernel, WCF), with a mean value of 0.18. For fruit- and kernel-related traits, CVs of diameter and shape index ranged from 0.05 to 0.07, which was lower than the CV of weight (0.16–0.27). For leaf-related traits, CVs ranged from moderate (0.12–0.18) to high (0.33). For fruit quality traits, CVs were moderate, at 0.18 and 0.11 for fatty acid content (FAC) and nervonic acid content (NAC), respectively.
Shannon–Wiener index (H) of the 21 traits had a mean value of 1.94 and varied from 1.51 (water content of dry kernel, WCD) to 2.06 (leaf length, LL, and leaf width, LW). For fruit-related traits, the pattern of H was typically consistent with CV. Pericarp thickness (PT) had the highest H of 2.02, whereas FSI had the lowest at 1.74. For kernel-related traits, however, the opposite trend was observed for H related to CV. The largest H (2.05) was found in the kernel transverse diameter (KTD, while the lowest (1.51) was found in the WCD. For leaf-related traits, the maximum H was 2.06 for LL and LW, and the lowest H was 1.98 for LSI. For fruit quality traits, the H of FAC and NAC were low, at 1.87 and 1.5, respectively.

3.2. Differences among Seven Natural Populations

The population significantly affected the variation of all the traits except for FTD (Table 2). Table 3 shows the means, standard deviations, and multiple comparisons of the 21 traits at the population level. P1 (Bama) had globose-shaped fruits (FSI = 0.94, KSI = 0.96) and had the greatest NAC (48.52%). P2 (Fengshan) had narrow and long leaves, of which the LW was the smallest at 4.9 cm, and the LI was the largest at 3.6. P3 (Leye) had the highest WCF and KDW but the least FAC and NAC. P4 (Lingyun) had the largest fruit, kernel, and leaves. The measurements for all the traits of P5 (Tianlin) ranged from moderate to high. P6 (Funing) had oblate-shaped fruits (FSI = 0.86, KSI = 0.89) and had the greatest FAC (61.94%). P7 (Guangnan) yields the heaviest fruit (KFW = 16.34), with a similar fruit shape to P6 and the highest FAC and NAC. Therefore, P6, P7 and P1 were more advantageous in selecting trees with higher fatty acid and nervonic acid content for future genetic improvement, while the remaining populations were better suited to preserving genetic diversity.
In general, the two populations in Yunnan province had oblate-shaped fruit and kernel characteristics (the smallest shape index of fruit and kernel) and high fatty acid content; the five populations in Guangxi province had relatively rounder shapes of fruits and kernels and a lower fatty acid content than the Yunnan population.

3.3. Principal Components and Cluster Analysis of the 21 Traits

The principal components of the 21 traits are shown in Table 4. The proportion of the cumulative contribution of the principal components with six eigenvalues (eigenvalue > 1) was 86.64%. The first principal component (PC1) explained 28.44% of observed variability and was strongly positively with the leaf-related traits, FSI, KSI, WCD, and strongly negatively correlated with the fatty acid contents (Figure 2A). The second principal component (PC2) explained 21.80% of observed variability and mostly defined the fruit and kernel-related traits. The third principal component, in addition to predominantly defining leaf-related traits, also contributed to the FAC. The fourth principal component mainly described FVD, PT, and KSI. The fifth and sixth principal components had lower contributions, which primarily characterized LI, KDW, and FAC, respectively.
The seven populations were not well clustered by the first two principal components except P7 (Guangnan; Figure 2A). However, they could be divided into two groups based on cluster analysis with UPGMA and Canberra distance (Figure 2B). Group I included P6 (Funing), P7, P4 (Lingyun) and P1 (Bama). The Group I population, which consisted of four populations distributed in the southern region, had wider and shorter leaves, larger FTD, and higher fatty acid and nervonic acid content in kernels. Group II comprised P2 (Fengshan), P3 (Leye) and P5 (Tianlin), which had higher WCF and WCD, larger LI, longer and narrower leaves, and lower nervonic acid content.

3.4. Correlations among the 21 Traits and Their Relationship with Geoclimatic Factors

Figure 3 depicts the pairwise trait correlations and their variation patterns of isolation by isolation by climatic and geographical distance. In general, within the four trait types, fruit, kernel, leaf, and fatty acid, traits were all strongly correlated with each other; however, among the four trait types, there was no relationship between leaf-related traits with fruit and kernel traits. FAC and NAC were positively correlated with KDW (r = 0.67 and 0.42) and negatively correlated with the shape index of fruit and kernel (r = −0.41 and −0.36) and leaf-related traits except for KS and LW. PL and LL were positively correlated with FSI and KSI and negatively correlated with KDW. It is indicated that trees with shorter leaves tend to produce oblate-shaped and heavier dry kernels, which yield more fatty acid and nervonic acid. Population-level pairwise morphological distance of kernel were positively correlated with fatty acid distance (Mantel r = 0.39, p < 0.05) and morphological distance of fruit (Mantel r = 0.47, p < 0.05).
Most of the relationships between the geographical and climate factors of the seven populations and the 21 traits were not significant (Figure 3). There was no correlation between AAT and the 21 traits. AAH was strongly positively correlated with LI (r = 0.78) and had a negative correlation with the KDW (r = −0.84), indicating the higher the humidity, the longer and narrower leaves and the smaller kernel dry weight. AAP and LON were highly positively correlated with FSI and KSI. From Yunnan to Guangxi province, with the increase in longitude and rainfall, the shape of the kernel gradually changed from oblate to globose. There was a negative correlation between ALT and PL (r = −0.80). LAT was found to be positively correlated with WCF (r = 0.86) and negatively correlated with FAC (r = −0.79), indicating that as latitude increased, the colder environment might reduce the fatty acid content. No significant isolation by climatic and geographical distance was found for morphological and fatty acid variability.

4. Discussion

Phenotypic diversity is affected by genetic and environmental factors and their interactions [28]. Phenotyping of desirable traits is the simplest and most direct method to survey and evaluate the diversity of forest germplasm resources, which is crucial for their reasonable conservation and genetic improvement. There is likely considerable phenotypic variation and diversity in forest tree species for the underlying complex environmental conditions, long-term geographic isolation, natural selection, and also phenotypic plasticity [29,30]. In the present study, we comprehensively evaluated the phenotypic diversity and variation in 21 phenotypic traits of fruits, kernels, leaves, fatty acid and nervonic acid content of M. oleifera from seven natural populations. We also investigated the relationships between these traits and geoclimatic factors, as well as the influence of these factors on phenotypic variability.
The coefficient of variation (CV) reflects the variation discovered in M. oleifera germplasm resources for each trait. At the individual level, great differences in CV were observed across the 21 traits, ranging from 0.05 (shape index of fruit, FSI) to 0.62 (water content of dry kernel, WCF) (Table 2). The size (e.g., diameter) of kernel/seed typically has limited and stable variation, for instance, at 0.06 in this study, in agreement with a previous report of M. oleifera (0.04–0.10) [20], Torreya grandis Fort. ex Lindl. (0.10) [15], Quercus variabilis Blume (around 0.10) [17]. The CVs of fruit size-related traits were slightly higher than those of kernel/seed, for instance, 0.07 in this study, 0.10–0.21 in Sorbus domestica L. [31], 0.12–0.18 in Actinidia eriantha Benth. [32] and 0.16–0.19 in Camellia oleifera Abel. [33]. The CVs of weight, as well as water content, had higher fluctuations, 0.16–0.62 in this study, 0.37 in Sweet Chestnut (Castanea sativa Mill.) [34], 0.06–0.40 in Sapindus mukorossi Gaertn. [35], which could be affected by fruit maturity and maternal factors, such as the position of the fruit in the tree. The CVs for fatty acid content (FAC) and nervonic acid content (NAC) were 0.18 and 0.11, respectively (Table 2), lower or similar to those reported for other tree species, like 0.20 for crude fat in S. domestica L. [31], 0.25 for crude fat in Castanea sativa Mill. [34], 0.41 for oil content in Acer truncatum Bunge [36]. In forest tree species, traits with lower CV generally have higher heritability [37], indicating that they were mainly controlled by the additive genetic effect, and it was capable of capturing higher genetic gain through traditional breeding approaches [38] or dissecting the underlying genetic architecture using contemporary association studies. On the other hand, for the traits with high CVs, such as the water content of fresh and dry kernels, to a large extent, the variation was caused by the environment and genetic-environment interaction effects. It requires a multi-year measurement in natural populations or multi-site provenance-progeny trials to obtain the genetic components in the variation in the forward studies.
The Shannon–Wiener index (H) reflects the diversity or evenness of the traits. At the individual level, the H of the 21 traits ranged from 1.51 (water content of dry kernels) to 2.06 (leaf width), with an average value of 1.94. The average level of diversity was comparable to that of Eucommia ulmoides Oliv. (H = 1.92) [39], slightly lower than that of Erythrophleum fordii Oliv. (H = 2.03) [22], both of which are Class II key protected species in China. It was similar to that of Populus cathayana Rehd. (H = 1.97) [40] and Catalpa fargesii Bur. (H = 1.92) [41]. The estimates of H had no significant relationship with those corresponding CVs. The lower the H, the fewer grades concentratedly distributed by the measurements. This phenomenon can be explained by the measurements from a small number of individuals showing some departure from those taken by the majority. In this study, the low H of the fatty acid content (FAC, 1.87) and nervonic acid content (NAC, 1.5), along with the moderate CV, demonstrated that M. oleifera has a great deal of breeding potential for genetic improvement.
It is crucial to identify selection indicators for excellent resources in forest tree species. For example, high-quality feeds, leaf biomass and crude protein content can be considered indicators in the selection of excellent provenance in Broussonetia papyrifera (L.) L’Hér. ex Vent. [42]. In this study, the following relationship between the traits was found: the larger the fruit, especially in terms of the transverse diameter, the higher the content of nervonic acid; meanwhile, small and wide short leaves and petioles tend to correspond to higher fat and nervonic acid content. As a result, small leaves, short petioles, and large fruits can be used as key indicators for the selection and breeding of varieties of M. oleifera with high FAC and NAC. In addition, this conclusion was also further supported by the result that population-level pairwise kernel morphological distance were positively correlated with fatty acid distance (Mantel r = 0.39, p < 0.05) and morphological distance of fruit (Mantel r = 0.47, p < 0.05). Leaves play a key role in effective water balance, light capture, and temperature regulation [43], as well as the petiole, which plays an important role in the adjustment of foliage and its inclination angles for optimal light capture [44]. Thus, leaf-related traits, such as leaf area development and leaf size, were robust indicators of biomass production in trees [45]. Similarly, seed size could be used as an indicator for the selection of seed sources [46]. Therefore, germplasm resources with high FAC and NAC could be selected simply based on the measurement of leaf and fruit characteristics to circumvent the time-consuming and costly chemical assays.
Longitude, latitude, and altitude, as key geographical factors, mainly affect plant traits through temperature and precipitation. With the increase in longitude and rainfall from west to east, the longitudinal diameters of fruits and kernels gradually increased, and their shape gradually changed from oblate to globose (Figure 3). With the increased humidity, the leaves tended to become narrower, and the dry kernel weight tended to decrease (Figure 3). The petiole tended to shorten with the gradual increase in altitude. As the latitude increased, the climate became colder, and the fresh kernel tended to become smaller with higher water content and lower FAC. Similar results have been documented in other species. In Eucommia ulmoides Oliver, latitude and precipitation were the key environmental factors that impacted the seed fatty acid composition [47]. The oil content was negatively correlated with latitude and annual rainfall in A. truncatum [36]. Latitude was the main effect factor in the variation of growth traits and nutritional traits [42]. The fruit length-to-width ratio showed a concave variation trend with longitude [17]. Therefore, annual average humidity and latitude should be considered when we introduce M. oleifera to the northern regions to expand its planting area.
With the increase of longitude and precipitation, the distribution area of M. oleifera expanded, which was consistent with the prediction that the distribution range would extend to the east and north under global climate change [48]. Historical research revealed that M. oleifera was dispersed in 22 counties (19 in Guangxi province and three in Yunnan province). Based on the climate data of 22 counties, Xie et al. [49] reported that the extreme minimum temperature for M. oleifera survival was between −5.8 °C and −0.3 °C, and the extreme maximum temperature was between 33.4 °C and 42.5 °C; heat and precipitation were the main factors affecting its distribution, and it was hard to adapt to high temperature and dry climate. According to Gong et al. [48], the coefficient of seasonal variation in temperature was the most important factor affecting the distribution and population scale of M. oleifera. Therefore, annual average humidity, latitude and temperature were the main limiting factors of M. oleifera distribution.
The clustering relationships and geographical distribution characteristics of the seven groups on 21 traits reflected, to some extent, the pattern of geographical variation dominated by latitudinal changes. The four populations distributed in the southern region had wider and shorter leaves, larger FTD, and higher fatty acid and nervonic acid content in kernels. At the same time, the three northern populations had higher WCF and WCD, larger LI, longer and narrower leaves, and lower NAC. The seven populations were not clustered strictly according to geographic distances, although the two populations in Yunnan province (Guangnan and Funing) were clustered together on the dendrogram (Figure 2B), indicating that geographic distance was not a critical factor in the differentiation of M. oleifera. Furthermore, the results of the Mantel test also suggested that morphological and fatty acid variability was not significantly associated with climatic and geographical distance (Figure 3). There was no significant correlation between the AAT and the 21 traits of the sample population. However, for M. oleifera, a small population with narrow distribution, temperature undoubtedly affects its species distribution and population size. Moreover, Gong et al. [48] used the maximum entropy model (MaxEnt) and 20 main environmental factors to predict the potential suitable range of M. oleifera: the longitude and latitude were 104°–107° E and 22°–26° N, respectively; the top four environmental factors were all temperature-related factors. The negligible contribution of AAT on population morphological diversity can be attributed to the little temperature difference resulting from the relatively narrow research region with a longitude spread of 2.2° and a latitude span of 1.22°.

5. Conclusions

There was abundant variation and diversity in 21 traits in seven natural populations of Malania oleifera. There is considerable breeding potential for genetic gain improvement of the fatty acid content and nervonic acid content in M. oleifera. Leaf length could be an indicator trait for superior tree selection because trees with shorter leaves tend to produce oblate shapes and heavier dry kernels, which yields more fatty acid and nervonic acid. The Yunnan province provenances had higher fatty acid content. Annual average humidity and latitude are key geographical factors affecting fatty acid content. Annual average humidity, latitude and temperature were the main limiting factors of M. oleifera distribution. Climatic and geographic distances were not critical factors in the morphological and fatty acid differentiation of M. oleifera.

Author Contributions

Conceptualization, H.L. and L.D.; investigation, H.L., R.D., H.F. and Y.L.; data curation, R.D., H.F. and Y.L.; formal analysis, H.L. and R.W.; writing—original draft preparation, H.L., R.W. and L.D.; writing—review and editing, Z.T., J.X. and W.S.; visualization, Y.Z.; funding acquisition, H.L. and L.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Natural Science Foundation of Guangxi Province (2020GXNSFAA297261) and the Natural Science Foundation of Beijing (6224061).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, S.K. Malania, a new genus of oil-yielding plant. Bull. Bot. Lab. North East. For. Inst. 1980, 1, 67–72. [Google Scholar]
  2. Sun, W. Malania oleifera. IUCN Red List Threat. Species 1998, e.T32361A9701100. [Google Scholar] [CrossRef]
  3. Ma, Y.; Chen, G.; Edward Grumbine, R.; Dao, Z.; Sun, W.; Guo, H. Conserving plant species with extremely small populations (PSESP) in China. Biodivers. Conserv. 2013, 22, 803–809. [Google Scholar] [CrossRef]
  4. Wang, X.; Wang, S. Research status and application perspective of nervonic acid. China Oils Fats 2010, 35, 1–5. [Google Scholar]
  5. Tang, T.-F.; Liu, X.-M.; Ling, M.; Lai, F.; Zhang, L.; Zhou, Y.-H.; Sun, R.-R. Constituents of the essential oil and fatty acid from Malania oleifera. Ind. Crops Prod. 2013, 43, 1–5. [Google Scholar] [CrossRef]
  6. China Oleaginous Plants Writing Committee. China Grease Plant; Science Press: Beijing, China, 1987. [Google Scholar]
  7. Zhou, Y.H.; Li, W.G.; Yi, F.P.; Liu, X.M. Determination of fatty acids in Malania oleifera oil by gas chromatography–mass spectrometry. Chin. J. Chromatogr. 2001, 19, 147–148. [Google Scholar]
  8. Li, Z.; Ma, S.; Song, H.; Yang, Z.; Zhao, C.; Taylor, D.; Zhang, M. A 3-ketoacyl-CoA synthase 11 (KCS11) homolog from Malania oleifera synthesizes nervonic acid in plants rich in 11Z-eicosenoic acid. Tree Physiol. 2021, 41, 331–342. [Google Scholar] [CrossRef]
  9. Yang, T.; Yu, Q.; Xu, W.; Li, D.Z.; Chen, F.; Liu, A. Transcriptome analysis reveals crucial genes involved in the biosynthesis of nervonic acid in woody Malania oleifera oilseeds. BMC Plant Biol. 2018, 18, 247. [Google Scholar] [CrossRef] [Green Version]
  10. Xu, C.Q.; Liu, H.; Zhou, S.S.; Zhang, D.X.; Zhao, W.; Wang, S.; Chen, F.; Sun, Y.Q.; Nie, S.; Jia, K.H.; et al. Genome sequence of Malania oleifera, a tree with great value for nervonic acid production. Gigascience 2019, 8, giy164. [Google Scholar] [CrossRef] [Green Version]
  11. Ramanatha Rao, V.; Hodgkin, T. Genetic diversity and conservation and utilization of plant genetic resources. Plant Cell Tissue Organ Cult. 2002, 68, 1–19. [Google Scholar] [CrossRef]
  12. Zhang, Y.; Dong, L.; Xie, Y.; Chen, D.; Sun, X. Altitude shape genetic and phenotypic variations in growth curve parameters of Larix kaempferi. J. For. Res. 2022. [Google Scholar] [CrossRef]
  13. Mckown, A.D.; Guy, R.D.; Klápště, J.; Geraldes, A.; Friedmann, M.; Cronk, Q.C.B.; El-Kassaby, Y.A.; Mansfield, S.D.; Douglas, C.J. Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa. New Phytol. 2014, 201, 1263–1276. [Google Scholar] [CrossRef] [PubMed]
  14. Zheng, Y.L.; Sun, W.B.; Zhou, Y.; Coombs, D. Variation in seed and seedling traits among natural populations of Trigonobalanus doichangensis (A. Camus) Forman (Fagaceae), a rare and endangered plant in southwest China. New For. 2009, 37, 285–294. [Google Scholar] [CrossRef]
  15. Dong, L.; Zeng, Y.; Wu, Y.; Huang, Y. Variations in phenotypic traits and chemical compositions of seeds from a natural population in Torreya grandis. J. Zhejiang A F Univ. 2014, 31, 224–230. [Google Scholar]
  16. Wu, H.; Meng, H.; Wang, S.; Wei, X.; Jiang, M. Geographic patterns and environmental drivers of seed traits of a relict tree species. For. Ecol. Manag. 2018, 422, 59–68. [Google Scholar] [CrossRef]
  17. Gao, S.; Ren, Y.; Masabni, J.; Zou, F.; Xiong, H.; Zhu, J. Influence of Geographical and Climatic Factors on Quercus variabilis Blume Fruit Phenotypic Diversity. Diversity 2021, 13, 329. [Google Scholar] [CrossRef]
  18. Ma, Y.; Bi, Q.; Li, G.; Liu, X.; Fu, G.; Zhao, Y.; Wang, L. Provenance variations in kernel oil content, fatty acid profile and biodiesel properties of Xanthoceras sorbifolium Bunge in northern China. Ind. Crops Prod. 2020, 151, 112487. [Google Scholar] [CrossRef]
  19. Li, Y.; Li, S.; Lu, X.; Wang, Q.; Han, H.; Zhang, X.; Ma, Y.; Gan, X. Leaf phenotypic variation of endangered plant Tetracentron sinense Oliv. and influence of geographical and climatic factors. J. For. Res. 2021, 32, 623–636. [Google Scholar] [CrossRef]
  20. Wang, S.H.; Chen, J.; Yang, W.; Hua, M.; Ma, Y.P. Fruiting character variability in wild individuals of Malania oleifera, a highly valued endemic species. Sci. Rep. 2021, 11, 23605. [Google Scholar] [CrossRef]
  21. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2014. [Google Scholar]
  22. Li, H.; Chen, D.; Xu, J.; Liu, G.; Pang, X.; Ye, J.; Mo, X.; Chen, H. Phenotypic diversity and variation in natural populations of Erythrophleum fordii, an endangered plant species. Sci. Silvae Sin. 2019, 55, 69–83. [Google Scholar]
  23. de Mendiburu, F. Agricolae: Statistical Procedures for Agricultural Research. Available online: https://cran.r-project.org/web/packages/agricolae/index.html (accessed on 1 April 2022).
  24. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
  25. Oksanen, J.; Simpson, G.L.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Solymos, P.; Stevens, M.H.H.; Szoecs, E.; et al. Vegan: Community Ecology Package. Available online: https://cran.r-project.org/web/packages/vegan/index.html (accessed on 1 April 2022).
  26. Hijmans, R.J. Geosphere: Spherical Trigonometry. Available online: https://cran.r-project.org/web/packages/geosphere/index.html (accessed on 1 April 2022).
  27. Schliep, K.; Potts, A.J.; Morrison, D.A.; Grimm, G.W. Intertwining phylogenetic trees and networks. Methods Ecol. Evol. 2017, 8, 1212–1220. [Google Scholar] [CrossRef] [Green Version]
  28. Falconer, D.S.; Mackay, T.F.C. Introduction to Quantitative Genetics, 4th ed.; Longman Group: Harlow, UK, 1996; ISBN 9780582243026. [Google Scholar]
  29. White, T.L.; Adams, W.T.; Neale, D.B. Forest Genetics, 1st ed.; CABI: Cambridge, MA, USA, 2007; Volume 167, ISBN 9781845932855. [Google Scholar]
  30. Poljak, I.; Kajba, D.; Ljubić, I.; Idžojtić, M. Morphological variability of leaves of Sorbus domestica L. in Croatia. Acta Soc. Bot. Pol. 2015, 84, 249–259. [Google Scholar] [CrossRef] [Green Version]
  31. Poljak, I.; Vahčić, N.; Liber, Z.; Tumpa, K.; Pintar, V.; Zegnal, I.; Vidaković, A.; Valković, B.; Kajba, D.; Idžojtić, M. Morphological and chemical diversity and antioxidant capacity of the service tree (Sorbus domestica L.) fruits from two eco-geographical regions. Plants 2021, 10, 1691. [Google Scholar] [CrossRef]
  32. Tao, J.; Wu, M.; Zhong, W.; Jiao, X.; Chen, S.; Jia, H.; Jia, D.; Huang, C. Changes in Fruit Quality and Sugar Components of Wild Actinidia eriantha of Different Varieties (Lines) at the Ripening Stage. Horticulturae 2022, 8, 824. [Google Scholar] [CrossRef]
  33. Gao, S.; Wang, B.; Liu, F.; Zhao, J.; Yuan, J.; Xiao, S.; Masabni, J.; Zou, F.; Yuan, D. Variation in Fruit Morphology and Seed Oil Fatty Acid Composition of Camellia oleifera Collected from Diverse Regions in Southern China. Horticulturae 2022, 8, 818. [Google Scholar] [CrossRef]
  34. Poljak, I.; Vahčić, N.; Liber, Z.; Šatović, Z.; Idžojtić, M. Morphological and Chemical Variation of Wild Sweet Chestnut (Castanea sativa Mill.) Populations. Forests 2022, 13, 55. [Google Scholar] [CrossRef]
  35. Sun, C.; Wang, J.; Duan, J.; Zhao, G.; Weng, X.; Jia, L. Association of Fruit and Seed Traits of Sapindus mukorossi Germplasm with Environmental Factors in Southern China. Forests 2017, 8, 491. [Google Scholar] [CrossRef] [Green Version]
  36. Chang, P.; Ma, J.; Xin, H.; Wang, S.; Chen, Z.; Hong, X.; Zhang, B.; Li, L. Comparative Study of the Fatty Acid Composition of the Acer truncatum Bunge from Different Producing Areas. Forests 2022, 13, 1409. [Google Scholar] [CrossRef]
  37. Dong, L.; Xie, Y.; Sun, X. Full-diallel-based analysis of genetic parameters for growth traits in Japanese larch (Larix kaempferi). New For. 2020, 51, 261–271. [Google Scholar] [CrossRef]
  38. Lai, M.; Dong, L.; Yi, M.; Sun, S.; Zhang, Y.; Fu, L.; Xu, Z.; Lei, L.; Leng, C.; Zhang, L. Genetic Variation, Heritability and Genotype × Environment Interactions of Resin Yield, Growth Traits and Morphologic Traits for Pinus elliottii at Three Progeny Trials. Forests 2017, 8, 409. [Google Scholar] [CrossRef] [Green Version]
  39. Li, H.; Du, H.; Jia, H.; Chen, H.; Xu, J.; Du, Q. Establishment of Male Core Collection of Eucommia ulmoides Based on Phenotypic Traits. Mol. Plant Breed. 2018, 16, 591–601. [Google Scholar]
  40. Cao, D.; Zhang, Y.; Cheng, X.; Xiang, X.; Zhang, L.; Hu, J. Genetic Variation of Leaf Phenotypic Traits in Different Populations of Populus cathayana. Sci. Silvae Sin. 2021, 57, 56–67. [Google Scholar]
  41. Xiao, Y.; Ma, W.-J.; Yi, F.; Yang, G.-J.; Wang, Q.-X.; Wang, P.; Wang, J.-H. Genetic Variation of Growth Traits and Genetic Diversity of Phenotypic Traits in Catalpa fargesii f. duclouxii Germplasm. Bull. Bot. Res. 2018, 38, 843–851. [Google Scholar]
  42. Zhou, P.; Que, Q.; Ke, B.; Cui, S.; Chen, X.; Zhou, W. Geographical Variation in the Growth and Nutritional Traits of Leaf Powder from Broussonetia papyrifera (L.) L’Hér. ex Vent. from Different Provenances. Forests 2022, 13, 868. [Google Scholar] [CrossRef]
  43. Kidner, C.A.; Umbreen, S. Why is Leaf Shape so Variable? Int. J. Plant Dev. Biol. 2010, 4, 64–75. [Google Scholar]
  44. Niinemets, Ü. Adjustment of foliage structure and function to a canopy light gradient in two co-existing deciduous trees. Variability in leaf inclination angles in relation to petiole morphology. Trees-Struct. Funct. 1998, 12, 446–451. [Google Scholar] [CrossRef]
  45. Ren, J.; Ji, X.; Wang, C.; Hu, J.; Nervo, G.; Li, J. Variation and genetic parameters of leaf morphological traits of eight families from Populus simonii × P. nigra. Forests 2020, 11, 1319. [Google Scholar] [CrossRef]
  46. Mamo, N.; Mihretu, M.; Fekadu, M.; Tigabu, M.; Teketay, D. Variation in seed and germination characteristics among Juniperus procera populations in Ethiopia. For. Ecol. Manag. 2006, 225, 320–327. [Google Scholar] [CrossRef]
  47. Du, Q.; Wang, L.; Liu, P.; Qing, J.; Sun, C.; Sun, Z.; Du, H. Fatty acids variation in seed of Eucommia ulmoides populations collected from different regions in China. Forests 2018, 9, 505. [Google Scholar] [CrossRef]
  48. Gong, M.; Wang, J.; Fu, X.; Kou, W.; Lu, N.; Wang, Q.; Lai, H. Suitable regions forecasting and environmental influencing factors of Malania oleifera in Yunnan and Guangxi. J. Nanjing For. Univ. Sci. Ed. 2022, 46, 44–52. [Google Scholar]
  49. Xie, W.; Chen, J.; Lai, J.; Shi, H.; Huang, K.; Liu, J.; Li, X. Analysis on Relationship between Geographic Distribution of Malania oleifera and Hydro-thermal Factors. J. Trop. Subtrop. Bot. 2009, 17, 388–394. [Google Scholar]
Figure 1. Locations of the seven populations sampled of M. oleifera. Upper panel: A map of China with grey indicating the studied area. Lower panels: The enlarged studied area, where solid black dots indicate the studied populations.
Figure 1. Locations of the seven populations sampled of M. oleifera. Upper panel: A map of China with grey indicating the studied area. Lower panels: The enlarged studied area, where solid black dots indicate the studied populations.
Forests 13 01733 g001
Figure 2. Biplot of the principal component analysis (A) and clustering relationship of the seven populations by UPGMA (B) based on 21 traits.
Figure 2. Biplot of the principal component analysis (A) and clustering relationship of the seven populations by UPGMA (B) based on 21 traits.
Forests 13 01733 g002
Figure 3. Correlations among the 21 traits and geoclimate factors and association between population and geographic distance.
Figure 3. Correlations among the 21 traits and geoclimate factors and association between population and geographic distance.
Forests 13 01733 g003
Table 1. Population, sample size, age, and geoclimatic factors for seven natural populations of M. oleifera.
Table 1. Population, sample size, age, and geoclimatic factors for seven natural populations of M. oleifera.
PopulationSample SizeSample AgeAnnual Average Temperature (°C, AAT)Annual Average Humidity (%, AAH)Average Annual Precipitation (mm, AAP)Latitude (° N, LAT)Longitude (° E, LON)Altitude (m, ALT)
Bama (P1)622–4620.7801505.924.12107.12532.6–541.2
Fengshan (P2)825–6019.4791509.324.47106.83698.2–953.2
Leye (P3)2030–8016.8821327.224.84106.40986.1–1027.8
Lingyun (P4)330–10020.4771669.624.33106.77705.6–777.8
Tianlin (P5)1035–10020.8811193.324.39106.32600.8–955.0
Funing (P6)1330–8019.8781103.523.65105.73819.9–878.2
Guangnan (P7)3730–8017.078994.724.24104.921217.1–1314.6
Table 2. Variation and diversity of the 21 traits in seven natural populations of M. oleifera (* and ** denote significance for p < 0.05 and p < 0.01, respectively.).
Table 2. Variation and diversity of the 21 traits in seven natural populations of M. oleifera (* and ** denote significance for p < 0.05 and p < 0.01, respectively.).
Traits CodeMeanSDMaxMinRangeCVHF Value
FruitFruit vertical diameter (mm)FVD36.7744.1626.282.5817.880.071.965.81 **
Fruit transverse diameter (mm)FTD41.3745.8125.852.7519.960.071.862.18
Shape index of fruitFSI0.891.060.700.050.360.051.7413.15 **
Fresh fruit weight (g)FFW36.1946.6610.215.8936.450.161.932.44 *
Pericarp thickness (mm)PT9.4213.344.051.639.290.172.023.19 **
KernelKernel vertical diameter (mm)KVD29.6234.5624.271.9010.290.062.029.05 **
Kernel transverse diameter (mm)KTD32.1236.7226.601.9510.120.062.054.14 **
Shape index of kernelKSI0.921.150.810.060.340.071.967.55 **
Fresh kernel weight (g)KFW15.1421.137.172.3613.960.162.006.55 **
Dry kernel weight (g)KDW8.8413.422.022.3911.400.272.0020.82 **
Water content of fresh kernel (g)WCF6.3014.082.732.2211.350.351.9225.50 **
Water content of dry kernel (%)WCD7.7528.203.544.8424.660.621.517.66 **
LeafPetiole length (mm)PL12.3317.158.412.088.740.172.0011.03 **
Knot spacing (mm)KS1.222.440.490.411.950.332.045.93 **
Leaf length (mm)LL19.1227.8010.503.4517.300.182.0610.88 **
Leaf width (mm)LW5.818.063.850.964.210.171.985.71 **
Leaf indexLI3.324.282.350.401.930.122.068.33 **
Leaf area (mm2)LA70.67128.9227.7423.35101.180.332.049.99 **
Leaf perimeter (mm)LP44.5664.0726.387.6137.690.172.0510.38 **
Fatty acidFatty acid content (%)FAC53.8567.8028.609.8139.200.181.8727.53 **
Nervonic acid content (%)NAC45.0251.4021.704.8729.700.111.693.91 **
Mean (SD) -----0.18 (0.13)1.94 (0.14)
Table 3. Mean value, standard deviation, and multiple comparison of 21 traits (Means with the same letter are not significantly different at the p < 0.05 level.).
Table 3. Mean value, standard deviation, and multiple comparison of 21 traits (Means with the same letter are not significantly different at the p < 0.05 level.).
Traits P1P2P3P4P5P6P7
FruitFVD38.22 ± 3.69 abc36.68 ± 2.03 bcd37.18 ± 2.63 bc40.49 ± 0.59 a38.98 ± 2.46 ab34.61 ± 3.28 d36.19 ± 1.2 cd
FTD40.91 ± 3.67 b40.05 ± 3.78 b40.53 ± 3.72 b43.9 ± 2.32 a42.12 ± 2.05 ab40.28 ± 2.35 b42.17 ± 1.64 ab
FSI0.94 ± 0.03 a0.92 ± 0.05 a0.92 ± 0.04 a0.92 ± 0.04 a0.92 ± 0.03 a0.86 ± 0.06 b0.86 ± 0.02 b
FFW35.76 ± 9.81 b32.68 ± 8.23 b34.15 ± 6.85 b42.54 ± 4.38 a37.97 ± 5.35 ab34.37 ± 5.64 b37.77 ± 3.19 ab
PT11.15 ± 1.37 a9.26 ± 2.57 b8.88 ± 1.26 b11 ± 0.77 a10.35 ± 1.04 ab8.82 ± 1.85 b9.29 ± 1.45 b
KernelKVD28.64 ± 2.73 cd29.09 ± 1.33 cd31.35 ± 1.54 ab32.34 ± 0.35 a30.01 ± 1.42 bc28.04 ± 1.92 d29.18 ± 1.29 cd
KTD29.76 ± 2.51 c30.79 ± 2.35 bc32.42 ± 1.59 ab32.9 ± 3.03 a31.77 ± 1.18 ab31.46 ± 2.44 abc32.88 ± 1.37 a
KSI0.96 ± 0.04 a0.95 ± 0.06 ab0.97 ± 0.05 a0.99 ± 0.1 a0.94 ± 0.04 ab0.89 ± 0.05 b0.89 ± 0.06 b
KFW11.66 ± 3.42 c14.03 ± 1.96 b15.22 ± 1.68 ab16.28 ± 3.39 a14.96 ± 1.76 ab13.77 ± 2.82 b16.34 ± 1.6 a
KDW7.7 ± 2.54 c8.03 ± 0.99 bc5.72 ± 1.73 d10.54 ± 1.62 a9.01 ± 1.67 abc9.47 ± 1.92 ab10.49 ± 1.31 a
WCF3.96 ± 0.92 c6 ± 1.26 b9.5 ± 2.13 a5.74 ± 2.26 b5.95 ± 1.18 b4.3 ± 1.24 c5.85 ± 0.98 b
WCD7.61 ± 4.02 abc9.14 ± 5.14 abc11.8 ± 5.64 a8.39 ± 0.18 abc10.53 ± 7.5 ab7.08 ± 2.77 bc4.73 ± 0.96 c
LeafPL15.02 ± 1.87 a11.62 ± 1.79 cd13.38 ± 1.93 abc13.1 ± 1.38 bc13.69 ± 1.82 ab12.97 ± 2.32 bc10.83 ± 0.95 d
KS1.41 ± 0.5 ab1.21 ± 0.21 bc1.21 ± 0.43 bc1.68 ± 0.35 a1.31 ± 0.22 abc1.6 ± 0.43 ab1 ± 0.31 c
LL19.63 ± 2.14 ab17.25 ± 2.8 bc21.99 ± 3.33 a20.3 ± 2.9 a21.15 ± 1.29 a20.67 ± 2.36 a16.7 ± 2.81 c
LI3.2 ± 0.25 bc3.6 ± 0.41 a3.56 ± 0.3 ab2.9 ± 0.06 c3.57 ± 0.3 ab3.19 ± 0.54 bc3.15 ± 0.33 c
LW6.15 ± 0.35 bc4.9 ± 1.07 d6.21 ± 0.91 ab7.03 ± 1.12 a5.97 ± 0.5 bc6.6 ± 0.99 ab5.32 ± 0.69 cd
LA74.02 ± 10.37 ab53.93 ± 20.12 c85.24 ± 24.65 a92.72 ± 28.93 a77.37 ± 9.51 a89.2 ± 19.55 a55.76 ± 16.45 bc
LP45.68 ± 4.7 ab39.81 ± 6.32 b50.6 ± 7.49 a47.59 ± 6.51 a48.45 ± 2.86 a48.65 ± 5.52 a39.41 ± 6.1 b
Fatty acidFAC55.08 ± 9.12 bc50.52 ± 7.54 c41.28 ± 6.87 d41.7 ± 0.6 d51.13 ± 7.99 c61.94 ± 4.65 a60.04 ± 4.55 ab
NAC48.52 ± 2.22 a42.59 ± 7.76 b41.89 ± 3.83 b48.4 ± 0.72 a44.46 ± 2.48 ab44.31 ± 7.01 ab46.79 ± 3.56 ab
Table 4. Principal component analysis of the 21 traits.
Table 4. Principal component analysis of the 21 traits.
Traits 123456
FruitFVD0.050.75−0.280.500.00−0.02
FTD−0.400.83−0.020.170.21−0.13
FSI0.59−0.10−0.380.44−0.320.13
FFW−0.340.870.010.200.12−0.01
PT−0.210.250.010.820.21−0.21
KernelKVD0.300.61−0.590.11−0.220.16
KTD−0.250.73−0.21−0.55−0.05−0.01
KSI0.51−0.09−0.380.62−0.200.17
KFW−0.290.73−0.30−0.40−0.190.20
KDW−0.670.410.300.02−0.230.30
WCD0.410.34−0.64−0.440.04−0.10
LeafPL0.660.070.210.180.100.29
KS0.490.340.540.05−0.10−0.14
LL0.820.340.34−0.090.200.12
LI0.41−0.04−0.26−0.080.710.44
LW0.600.400.57−0.05−0.28−0.20
LA0.730.400.50−0.10−0.07−0.06
LP0.800.360.39−0.100.150.10
WCF0.59−0.14−0.36−0.120.03−0.46
Fatty acidFAC−0.66−0.050.500.11−0.120.39
NAC−0.580.270.200.180.30−0.33
Eigen value5.974.583.012.371.171.10
Contribution28.4421.8014.3311.275.555.26
Cumulative contribution28.4450.2464.5775.8481.3986.64
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, H.; Wang, R.; Tian, Z.; Xu, J.; Sun, W.; Duan, R.; Fu, H.; Li, Y.; Zhang, Y.; Dong, L. Phenotypic Variation and Diversity in Fruit, Leaf, Fatty Acid, and Their Relationships to Geoclimatic Factors in Seven Natural Populations of Malania oleifera Chun et S.K. Lee. Forests 2022, 13, 1733. https://doi.org/10.3390/f13101733

AMA Style

Li H, Wang R, Tian Z, Xu J, Sun W, Duan R, Fu H, Li Y, Zhang Y, Dong L. Phenotypic Variation and Diversity in Fruit, Leaf, Fatty Acid, and Their Relationships to Geoclimatic Factors in Seven Natural Populations of Malania oleifera Chun et S.K. Lee. Forests. 2022; 13(10):1733. https://doi.org/10.3390/f13101733

Chicago/Turabian Style

Li, Hongguo, Ruizhen Wang, Zuwei Tian, Jihuang Xu, Wensheng Sun, Runmei Duan, Hao Fu, Yunmu Li, Yalin Zhang, and Leiming Dong. 2022. "Phenotypic Variation and Diversity in Fruit, Leaf, Fatty Acid, and Their Relationships to Geoclimatic Factors in Seven Natural Populations of Malania oleifera Chun et S.K. Lee" Forests 13, no. 10: 1733. https://doi.org/10.3390/f13101733

APA Style

Li, H., Wang, R., Tian, Z., Xu, J., Sun, W., Duan, R., Fu, H., Li, Y., Zhang, Y., & Dong, L. (2022). Phenotypic Variation and Diversity in Fruit, Leaf, Fatty Acid, and Their Relationships to Geoclimatic Factors in Seven Natural Populations of Malania oleifera Chun et S.K. Lee. Forests, 13(10), 1733. https://doi.org/10.3390/f13101733

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop