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

: 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. Coefﬁcient 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 signiﬁcantly inﬂuenced 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 signiﬁcant 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. FSI, acid The second component (PC2) explained 21.80% of observed variability and mostly deﬁned the fruit and kernel-related traits. The third principal component, in addition to predominantly deﬁning leaf-related traits, also contributed to the FAC. The fourth principal component mainly de-scribed FVD, PT, and KSI. The ﬁfth and sixth principal components had lower contributions, which primarily characterized LI, KDW, and FAC, respectively.


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 (C 24 H 46 O 2 , 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.

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

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

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 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].
where p i 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].

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.

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.

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

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

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.

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.