Next Article in Journal
Ethnolichenology—The Use of Lichens in the Himalayas and Southwestern Parts of China
Previous Article in Journal
Bacterial Number and Genetic Diversity in a Permafrost Peatland (Western Siberia): Testing a Link with Organic Matter Quality and Elementary Composition of a Peat Soil Profile
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diversity
Volume 13
Issue 7
10.3390/d13070329
diversity-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Geographical and Climatic Factors on Quercus variabilis Blume Fruit Phenotypic Diversity

by
Shuang Gao
1,
Yue Ren
2,3,
Joseph Masabni
4,
Feng Zou
1,
Huan Xiong
1,* and
Jingle Zhu
2,*
1
Key Laboratory of Cultivation and Protection for Non-Wood Forest Trees, Ministry of Education, Central South University of Forestry and Technology, Changsha 410004, China
2
Paulownia Research and Development Center, State Forestry Administration, Zhengzhou 450003, China
3
Forestry College, Henan Agriculture University, Zhengzhou 450002, China
4
Texas A & M Agri Life Research and Extension Center, 17360 Coit Rd, Dallas, TX 75252, USA
*
Authors to whom correspondence should be addressed.
Diversity 2021, 13(7), 329; https://0-doi-org.brum.beds.ac.uk/10.3390/d13070329
Submission received: 7 June 2021 / Revised: 15 July 2021 / Accepted: 15 July 2021 / Published: 18 July 2021
Browse Figures
Versions Notes

Abstract

:
Quercus variabilis Blume is one of the most ecologically valuable tree species in China and is known to have adaptive mechanisms to climate change. Our objective was to quantify the variation pattern in the fruit morphology of Q. variabilis. Fruit samples were collected from 43 natural populations in autumn of 2019. Our results indicated that the coefficient of variation (CV) of the fruit length (FL) and fruit width (FW) were 10.08% and 11.21%, respectively. There were significant differences in the FL, FW, and fruit length-to-width ratios (FL/FW) among the studied populations. Also, there was a significant positive correlation between the FW and FL. The FL decreased with increasing precipitation in the wettest quarter (PWQ). A concave trend was observed in the variations in FL with the equivalent latitude (ELAT), longitude (LON), annual mean air temperature (MAT), and annual precipitation (AP). A similar concave trend was observed for the FL/FW with LON, MAT, and AP. A positive correlation was observed between the FW, FL and FL/FW, and the ELAT. The cluster analysis revealed five groups of the 43 natural populations. Our study findings suggests that Q. variabilis has high levels of phenotypic plasticity for geographical and climatic factors.

1. Introduction

Phenotypic variation among a population of the same species usually reflects the plant’s response to different environmental pressures [1,2]. The phenotypic diversitycommonly refers to leaf, flower, fruit, and seed morphology [3,4,5,6,7,8]. In general, seed size is considered a component of a co-evolving complex of variables, including plant biomass, dispersal, niche specialization, seed dormancy, and competitiveness [9]. Seed size also affects plant dissemination, settlement, and individual development [10,11]. Seed phenotype is malleable to changes in the environment that allow the plant to be better fitted to environmental conditions [1].
Studying the variation pattern of seed morphology and its relationship with environmental factors helps in the understanding of plant plasticity in response to changes in environmental factors and reveals the role of genetics and the environment in plant ecological adaptation [12]. The diversity of fruit morphology often indicates adaptation to different dispersal modes [13]. Recent research has focused on fruit size differences [14,15] and the relationship between the environment and fruit characteristics [16,17,18]. Maranz and Wiesman [19] found that temperature and rainfall had significant effects on fruit size. Wu et al. [20] indicated that fruit size (length and width) decreased with longitude from west to east. Also, Liu et al. [21] noted an increase in fruit width and fruit weight with longitude, and a gradual decrease with increasing latitude. Leal-Sáenz et al. [22] suggested that populations found in warmer and wetter climates had larger fruits. Soil properties can influence fruit size as well. For instance, the sweet chestnut populations from acidophilic and thermophilic forests were characterized with smaller fruits. In contrast, populations from forests of mesophilic character were characterized with larger fruits [23].
More than 400 species of the genus Quercus L. (Fagaceae) are distributed in the northern hemisphere, and there are 51 species of Quercus in China [24]. Quercus variabilis Blume is an important deciduous broadleaved oak of benefit to the ecology, economy, and culture in China [25,26]. At present, research on Q. variabilis has focused mainly on biological characteristics [27,28], population ecosystems [25], and physiological characteristics [29,30]. In this study, we analyzed the morphological characteristics of fruit from 43 natural populations of Q. variabilis. To understand the morphological variation of Q. variabilis fruits due to environmental influences, we asked: what is the degree of fruit phenotypic diversity within and among populations of Q. variabilis? What are the factors playing key roles in the fruit morphological variation of the Q. variabilis population? Thus, the aim of this study was to quantify the phenotypic diversity of Q. variabilis fruits and the response to climate factors.

2. Materials and Methods

2.1. Study Area

In the geographical distribution area of Q. variabilis in China [31], the climate may be temperate, warm temperate, north subtropical, and south subtropical. The soil types in these regions may be dark brown, cinnamon, yellow, red, and other main zonal soil types [12]. Forty-three typical natural secondary forests of Q. variabilis were selected in seventeen Chinese provinces (Figure S1; Table 1). The geographical span of the selected region ranged between 23.33° N (Honghe Prefecture, Yunnan) and 43.83° N (Urumqi city, Xinjiang), and between 87.62° E (Urumqi city, Xinjiang) and 123.30° E (Liaoyang city, Liaoning). The annual precipitation ranged between 231 mm (Urumqi city, Xinjiang) and 1556 mm (Nanchang city, Jiangxi). The wettest season precipitation ranged from 87 mm in Urumqi, Xinjiang, to 761 mm in Honghe Prefecture, Yunnan. The average temperature was 7.06 °C in Urumqi, Xinjiang, and 17.78 °C in Changsha city, Hunan. The highest monthly temperature ranged from 20.16 °C in Linzhi city, Tibet, to 33.42 °C in Nanchang city, Jiangxi (Table 1).

2.2. Sample Collection

In each forest stand, we selected a study area with healthy Q. variabilis trees and with little to no human interference. We established a 20 m × 20 m plot and recorded the latitude, longitude, and altitude of the test plot [12]. At the peak of fruit maturity in autumn 2019, at least 100 fully developed, disease-free fruits were collected from each plot. Fruits were sterilized by soaking in 5% chlorine bleach for ten minutes [32]. The samples were placed in separate nylon bags and air dried, then stored at 2 °C until the measurements were done [33].

2.3. Geographic Information and Climate Data of the Sample Plot

The climate data for all sites, including the annual mean air temperature (MAT), annual precipitation (AP), maximum temperature of warmest month (MTW), and precipitation of wettest quarter (PWQ), were obtained from the Global Climate and Weather Data website (http://www.worldclim.org (accessed on 15 December 2020)) [34] and are presented in Table 1. We used ArcMap10.8 to extract climate data from the WorldClim data website according to Slave [35]. Due to the large differences in altitude among sites, we converted all latitudes into equivalent latitudes (ELAT), as proposed by Alena et al. [36], in order to determine the true effect of latitude and eliminate the influence of altitude factors. The conversion formulas used are as follows:
equivalent latitude = latitude + (altitude −300)/200 (for altitude lower than 300 m)
equivalent latitude = latitude + (altitude −300)/140 (for altitude greater than 300 m)

2.4. Determination of Fruit Morphological Characters

The fruit width and fruit length of at least 20 fruits in each population were measured by vernier caliper with an accuracy of 0.01 mm [30]. Fruit width (FW) refers to the distance between the widest points on the left and right side of the fruit, and fruit length (FL) refers to the distance from the bottom to the top of the fruit. All measurements were accurate to two decimal places. Each fruit was measured three times and the data was averaged. The fruit length-to-width ratio (FL/FW) of each fruit was also calculated [17].

2.5. Statistical Analysis

First, the main statistics (arithmetic means, maximum value, minimum value, standard deviation, and coefficient of variation) were calculated for the studied characteristics and populations. Second, a oneway analysis of variance (ANOVA), as implemented in SPSS 25.0 software (IBM company, New York, NY, USA), was used to test the differences between the studied populations [37]. Significant differences among means were assessed using Duncan’s multiple comparison at p ≤ 0.05. Origin 2019b software (Origin Lab company, Northampton, MA, USA) was used to quantify the variability in phenotypic traits as influenced by regional climatic factors [38]. To eliminate dimension and the order of magnitude, the original data were standardized, converted into dimensionless data with a mean of 0 and a variance of 1, and then principal component analysis (PC) was performed. Cluster analysis was carried out by the between groups linkage method and measured according to squared Euclidean distance. Canoco5 (Microcomputer Power, New York, NY, USA) was used in the principal component analysis (PCA) of Q. variabilis fruits morphology [39]. These statistical analyses were performed using SPSS 25.0 [1]. Correlation between phenotypic traits and environmental factors were analyzed using Pearson correlation analysis by SPSS 25.0 [40]. In addition, we computed and tested the correlations between: (1) the matrix of geographical distances between pairs of populations and the matrix of morphological differences among populations—the isolation-by-distance pattern [41], and between (2) the matrix of environmental distances and the matrix of morphological differences among populations—the isolation-by-environmental distance [42]. The significance level was assessed after 9,999 permutations, and the simple Mantel test was performed with the R package “Vegan” [43].

3. Results

3.1. Fruit Morphological Characters and Variation Characteristics

The width, length, and length-to-width ratio of fruits were significantly different among the various populations (Table 2 and Table 3). The variation in fruit width among groups ranged between 7.96 mm and 27.17 mm, with an average of 17.35 mm. Fruits from the SXH site had the largest width (20.93 mm), while the smallest fruit width was observed in fruits from ESH (11.70 mm). Fruit length ranged between 9.96 mm and 35.04 mm, with an average of 19.97 mm. The average fruit length was largest in UQX (31.69 mm), and shortest in fruits came from BDH (17.17 mm). Fruit length-to-width ratio of fruits ranged between 0.70 and 2.31, with an average value of 1.17. The mean variation coefficient of fruit width, fruit length, and fruit length-to-width ratio among the different populations were 11.21%, 10.08%, and 14.48%, respectively (Table 2). In addition, we revealed that there was a significant correlation between Q. variabilis fruit width and fruit length of Q. variabilis (R2 = 0.17, p ≤ 0.0001) (Figure S2).

3.2. The Relationship in the Variation Pattern between Fruit Morphology and Geographical and Environmental Factors

Results of the correlation analysis between the various morphological indicators of Q. variabilis fruits and the geographical and ecological factors are presented in Figure 1. Fruit length and ELAT had a concave variation trend (R2 = 0.12; p = 0.029; y = 0.016x2 − 1.1x + 37.94), with a minimum value at 33°–35° N (Figure 1A). Fruit length also had a concave variation trend with longitude (R2 = 0.43; p < 0.0001; y = 0.014x2 − 3.19x + 198.26), with a minimum value at 111°–113° E (Figure 1B). The fruit length-to-width ratio did not correlate with ELAT, but showed a concave variation trend with longitude (R2 = 0.12; p = 0.027; y = 0.0001x2 − 0.15x + 9.55), with a minimum value at 111°–113° E (Figure 1B).
Fruit length and annual mean air temperature showed a concave variation trend (R2 = 0.16; p = 0.013; y = 0.096x2 − 2.65x + 37.58) and reached a minimum value at 13.46–14.65 °C (Figure 1C). Fruit length-to-width ratio and annual mean air temperature also showed a concave variation trend (R2 = 0.17; p = 0.01; y = 0.009x2 − 0.24x + 2.69) and reached a minimum value at 12.93–13.72 °C (Figure 1C). However, fruit length and length-to-width ratio did not correlate with the maximum temperature of the warmest month (Figure 1D).
Fruit length and annual precipitation showed a concave variation trend (R2 = 0.36; p < 0.0001; y = 0.0001x2 − 0.024x + 31.29) and reached a minimum value at 974–1071 mm (Figure 1E). The fruit length-to-width ratio and annual precipitation showed a concave variation trend (R2 = 0.34; p < 0.0001; y = 0.0001x2 − 0.002x + 2.02) and reached a minimum value at 939–1004 mm (Figure 1C). Fruit length was negatively correlated with precipitation in the wettest quarter (R2 = 0.10; p = 0.029; y = −0.005x + 22.53) (Figure 1F). But there was no correlation between the fruit length-to-width ratio and precipitation in the wettest quarter.
Fruit width, fruit length, and fruit length-to-width ratio were positively correlated with the equivalent latitude (ELAT) (Table 4). However, fruit length was significantly negatively correlated with precipitation of the wettest quarter and longitude.
The Simple Mantel test identified significant correlations between the morphological, geographic, and climatic distance matrices (Table 5). Correlations were higher between morphological and climatic matrices (r = 0.988, p < 0.001), and weaker, but significant, correlations were observed between the morphological and geographical matrices (r = 0.633, p < 0.001).

3.3. Principal Component Analysis and Cluster Analysis

Based on the results of principal component analysis (Figure 2, Table 6), 99.60% of the total variance was explained by the first two principal components. The principal component score, associated to each variable on the three principal components, identifies the variables that mostly defined them (Table 6). The PC1, showing 57.57% of total variability, was positively correlated with following variables: FW and FL/FW. The PC2 indicated 42.03% of cumulative variance was positively correlated with FL and FW.
The 43 natural populations of Q. variabilis were clustered into five groups with a Euclidean distance of five as the threshold (Figure 3). TBH, SXH, HNA, ZSH, BKH, JGH, JSH, CZA, JSA, WHH, YHZ, LAZ, BDH, NCJ, BSG, MYB, PGB, DLL, YTS, LYL, MTB, TSG, XCS, LSH, SNS, DFH, SYS, XYS, WNS, NZH, JYH, and BJS populations were clustered into the first group, ALG, HHY, XGG, ANY, KMY, and ZYY populations into the second, XPH, CSH, and ESH into the third, LZT was in the fourth, and UQX was in the fifth.

4. Discussion

Indeed, fruit morphology can vary significantly among locations, but plants of the same genus can maintain similar fruit size. In this paper, the variation in the average fruit length-to-width ratios of Q. variabilis ranged between 0.96 and 1.78. Chen et al. [44] reported that the variation range of the fruit length-to-width ratio of Q. virginiana Mill. was between 1.19 and 1.93. Li et al. [45] indicated that the fruit length-to-width ratio of Q. mongolica Fisch. ex Ledeb. was between 1.21 and 1.49. Chang et al. [46] observed that the fruit length-to-width ratio of Q. cocciferoides Hand.-Mazz. ranged between 0.86 and 1.48. These results are similar with our data. Therefore, the variation in the fruit length-to-width ratios in the Quercus genus is relatively moderate. In addition, in the present study, the fruit width and fruit length of Q. variabilis showed significant differences among the various populations, with a relatively low intrapopulation variability. These results are consistent with the results published by Zhou et al. [12].
The variation of fruit phenotypic diversity is closely associated with geographical area [1,47]. The geographic isolation from other regions leads to fruit phenotypic variation [20]. Fu et al. [48] showed that the fruit morphological characteristics of Cornus officinalis Sieb. have obvious geographical characteristics. Similar trends were observed in our study. Fruit morphological traits of Q. variabilis exhibited geographic patterns. The hierarchical tree dendrogram of the 43 populations of Q. variabilis showed that populations with close geographic distance, such as the SYS, SNS, WNS, and XYS populations, and the XPH and CSH populations are clustered together. This indicates the regional and continuous pattern of geographic variation. However, a few distant populations were found to be clustered together, such as MYB and LAZ populations, demonstrating random variation patterns. These results indicated that the fruit phenotypic diversity of Q. variabilis represented three geographical variation patterns: continuous variation, regional variation, and random variation. In addition, the populations of Q. variabilis were characterized by significant phylogeographic structure. The results of previous studies showed morphological variability of seeds along the latitudinal gradient [21,49]. Rewicz et al. [49] found that there was a trend of an increase in seed morphology towards the higher latitude areas. Equally, the results of our study showed that the fruit length, fruit width, and fruit length-to-width ratio of Q. variabilis had a positive correlation with the equivalent latitude. This correlation may reflect the adaptability and high phenotypic plasticity of Q. variabilis.
In addition to geographical factors, climatic factors also affected the fruit morphology. Fruit phenotypic diversity provides useful results for understanding the adaptation mechanism of plants under different climatic conditions [1,49]. The phenotypic variation was the synthetic effect of multiple climatic factors [50]. The variation trend of Q. variabilis fruit width and fruit length, as affected by climatic factors, showed regular variations in fruit morphological characteristics (the fruit width, fruit length, and fruit length-to-width ratio). The fruit length-to-width ratio showed a concave variation trend as affected by annual mean air temperature and annual precipitation, and was similar to fruit length. This showed that the fruit length-to-width ratio of Q. variabilis had a synergistic variation in different climatic conditions. Moreover, our results indicated that the fruit width of Q. variabilis was negatively correlated with temperature and precipitation, similar to the responses of Pinus tabuliformis Carr [50].
Early research on Quercus showed that phenotypic diversity was tightly related to local climatic factors with respect to temperature or rainfall [30,51,52]. Thus, the phenotypic characters of Q. variabilis should be variable in response to the changes of climatic factors because of its large distribution area. Generally, in a climate of low temperature and low annual precipitation, fruit length is negatively correlated with temperature and rainfall. On the contrary, in a climate of high temperature and high annual precipitation, fruit length is positively correlated with temperature and rainfall [12,53,54,55]. Similarly, we found various correlation patterns between those phenotypic traits and the climatic factors. In our study, the fruit length of Q. variabilis was negatively correlated with precipitation in the wettest quarter. During the development period of fruits between July and September, precipitation in the wettest quarter significantly affected the morphology of Q. variabilis fruits. In addition, when the annual precipitation was below 1000 mm, the fruit length was negatively correlated with the annual precipitation. Besides, we found that when the annual precipitation was greater than 1000 mm, the fruit length increased with increasing annual precipitation, and when the temperature was above 13 °C, the fruit length increased with the annual mean air temperature. In addition, we observed that when the temperature was less than 13 °C, the fruit length was negatively correlated with annual mean air temperature. Those various correlation patterns between phenotypic traits and the climatic factors have resulted in a strong adaptability to climate changes.
According to the results of the simple Mantel tests, it was determined that both geographical and climate variables affect the structure of Q. variabilis populations [56,57]. Similarly, Poljak et al. [43] suggested that phenotypic divergence of Alnus incana (L.) Moench subsp. incana populations was the result of a significant level of isolation both by distance and by the environment. Furthermore, DeWoody et al. [42] revealed that adaptive differentiation and persistent isolation by colonization (IBC) acted in combination to produce the genetic and morphological patterns observed in Populus nigra L. populations. In our study, we found that acorn morphological variation was significantly correlated with both climatic heterogeneity and geographic distance.

5. Conclusions

Our study revealed that the phenotypic diversity of Q. variabilis fruit is related to both geographical and climatic factors. In terms of geographical factors, the acorn morphology of Q. variabilis was significantly different among and within the 43 natural populations, with gradations in different geographical distribution areas. Moreover, the hierarchical tree dendrogram of the five major groups reflects diversified geographic variation patterns of Q. variabilis fruit phenotype, namely, random variation, continuous variation, and regional variation. In terms of climatic factors, Q. variabilis fruit phenotypic diversity was tightly related to local climate factors with respect to temperature and rainfall. In summary, the fruits of Q. variabilis have high phenotypic plasticity, which is conducive to the expansion of its adaptation range.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d13070329/s1, Figure S1: Location of the 43 study sites, covering the geographical distribution of Quercus variabilis in China. Figure S2: Correlation analysis between fruit width (FW) and fruit length (FL) of 43 populations of Quercus variabilis.

Author Contributions

S.G. contributed the data analysis and wrote the manuscript. H.X., F.Z., Y.R., and J.Z. collected the samples. J.M., H.X., J.Z., and F.Z. revised the manuscript. J.Z. and H.X. secured funding for the research. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the basic scientific research business expenses of the central public welfare scientific research institutes of the Chinese Academy of Forestry (CAFYBB2018ZB001-11), the key research and development program of Hunan Province (2020NK2017), and Innovation and Entrepreneurship training program of Central South University of Forestry and Technology ([2021] no.34).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data in this study are available in presented Tables.

Conflicts of Interest

The authors have no conflict of interest.

References

  1. Li, Y.; Li, S.; Lu, X.H.; Wang, Q.Q.; Han, H.Y.; Zhang, X.M.; Ma, Y.H.; Gan, X.H. Leaf phenotypic variation of endangered plant Tetracentron sinense Oliv. And influence of geographical and climatic factors. J. For. Res. 2020, 32, 623–636. [Google Scholar] [CrossRef]
  2. Alcántara-Ayala, O.; Oyama, K.; Ríos-Muñoz, C.; Rivas, G.; Ramirez-Barahona, S.; Luna-Vega, I. Morphological variation of leaf traits in the Ternstroemia lineata species complex (Ericales: Penthaphylacaceae) in response to geographic and climatic variation. PeerJ 2020, 8, e830. [Google Scholar] [CrossRef] [Green Version]
  3. Thomas, C.L.; Alcock, T.D.; Graham, N.S.; Hayden, R.; Matterson, S.; Wilson, L.; Young, S.D.; Dupuy, L.X.; White, P.J.; Hammond, J.P.; et al. Root morphology and seed and leaf ionomic traits in a Brassica napus L. diversity panel show wide phenotypic variation and are characteristic of crop habit. BMC Plant Biol. 2016, 16, 214. [Google Scholar] [CrossRef] [Green Version]
  4. Jasińska, A.K.; Rucińska, B.; Kozlowski, G.; Bétrisey, S.; Safarov, H.; Boratyńska, K.; Boratyński, A. Morphological differentiation of leaves in the relict tree Zelkova carpinifolia (Ulmaceae). Dendrobiology 2015, 74, 109–122. [Google Scholar] [CrossRef] [Green Version]
  5. Miljković, D.; Čortan, D. Morphometric and morphological analysis of Populus nigra L. leaves in flooded regions. Sumar. List 2020, 144, 139–147. [Google Scholar] [CrossRef]
  6. 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]
  7. Gómez, J.M.; Perfectti, F.; Armas, C.; Narbona, E.; González-Megías, A.; Navarro, L.; DeSoto, L.; Torices, R. Within-individual phenotypic plasticity in flowers fosters pollination niche shift. Nat. Commun. 2020, 11, 4019. [Google Scholar] [CrossRef] [PubMed]
  8. Daničić, V.; Kovačević, B.; Ballian, D. Variability in fruit morphology of European sweet chestnut (Castanea sativa Mill.) in natural populations in Bosnia and Herzegovina. Sumar. List 2018, 142, 517–528. [Google Scholar] [CrossRef] [Green Version]
  9. Pluess, A.R.; Schütz, W.; Stocklin, J. Seed weight increases with altitude in the Swiss Alps between related species but not populations of individual species. Oecologia 2005, 144, 55–61. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Yi, F.Y.; Wang, Z.R.; Baskin, C.C.; Baskin, J.M.; Ye, R.H.; Sun, H.L.; Zhang, Y.Y.; Ye, X.H.; Liu, G.F.; Yang, X.J.; et al. Seed germination responses to seasonal temperature and drought stress are species-specific but not related to seed size in a desert steppe: Implications for effect of climate change on community structure. Ecol. Evol. 2019, 9, 2149–2159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Wirth, L.R.; Graf, R.; Gugerli, F.; Landergott, U.; Holdergger, R. Between-year variation in seed weights across altitudes in the high-alpine plant Eritrichium nanum. Plant Ecol. 2010, 207, 227–231. [Google Scholar] [CrossRef] [Green Version]
  12. Zhou, X.; He, Z.B.; Kang, H.Z.; Sun, X.; Liu, C.J. Variations of seed morphology related to climate for Quercus variabilis across temperate subtropical China. Chin. J. Plant Ecol. 2013, 37, 481–491. [Google Scholar] [CrossRef]
  13. Liu, M.L.; Yu, W.B.; Li, D.Z.; Mill, R.; Wang, H. Seed morphological diversity of Pedicularis (Orobanchaceae) and its taxonomic significance. Plant Syst. Evol. 2013, 299, 1645–1657. [Google Scholar] [CrossRef]
  14. Kaliniewicz, Z.; Markowski, P.; Anders, A.; Jadwisieńczak, K.; Żuk, Z.; Krzysiak, Z. Physical properties of seeds of eleven fir species. Forests 2019, 10, 142. [Google Scholar] [CrossRef] [Green Version]
  15. Velázquez-Rosas, N.; Ruiz-Guerra, B.; Sánchez-Coronado, M.E.; Buen, A.G.; Orozco-Segovia, A. Morphological variation in fruits and seeds of Ceiba aesculifolia and its relationship with germination and seedling biomass. Bot. Sci. 2017, 95, 81–91. [Google Scholar] [CrossRef] [Green Version]
  16. Kijowska-Oberc, J.; Staszak, A.M.; Wawrzyniak, M.K.; Ratajczak, E. Changes in proline levels during seed development of orthodox and recalcitrant seeds of genus Ace in a climate change scenario. Forests 2020, 11, 1362. [Google Scholar] [CrossRef]
  17. Liu, G.F.; Zang, R.G.; Liu, H.; Bai, Z.Q.; Guo, Z.J.; Ding, Y. Geographic variation of seed morphological traits of Picea schrenkiana var. tianschanica in Tianshan Mountains, Xinjiang of Northwest China. Chin. J. Appl. Ecol. 2012, 23, 1455–1461. [Google Scholar]
  18. Wang, Y.J.; Wang, J.J.; Lai, L.M.; Jiang, L.H.; Zhuang, P.; Zhang, L.H.; Zheng, Y.R.; Jerry, M.; Baskin, J.M.; Baskin, C.C. Geographic variation in seed traits within and among forty-two species of Rhododendron (Ericaceae) on the Tibetan plateau: Relationships with altitude, habitat, plant height, and phylogeny. Ecol. Evol. 2014, 4, 1913–1923. [Google Scholar] [CrossRef]
  19. Maranz, S.; Wiesman, Z. Evidence for indigenous selection and distribution of the shea tree, Vitellaria paradoxa and its potential significance to prevailing parkland savanna tree patterns in sub-Saharan Africa north of the equator. J. Biogeogr. 2003, 30, 1505–1516. [Google Scholar] [CrossRef]
  20. Wu, H.; Meng, H.J.; Wang, S.T.; Wei, X.Z.; Jiang, M.X. Geographic patterns and environmental drivers of seed traits of a relict tree species. For. Ecol. Manag. 2018, 422, 59–68. [Google Scholar] [CrossRef]
  21. Liu, Z.L.; Yu, M.K.; Ma, Y.; Tang, L.Z.; Fang, S.Z. A trend surface analysis of geographic variation in the traits of seeds and seedlings from different Quercus acutissima provenances. Acta Ecol. Sin. 2011, 31, 6796–6804. [Google Scholar]
  22. Leal-Sáenz, A.; Waring, K.M.; Menon, M.; Cushman, S.A.; Eckert, A.; Flores-Rentería, L.; Hernández-Díaz, J.C.; López-Sánchez, C.A.; Martínez-Guerrero, J.H.; Wehenkel, C. Morphological differences in Pinus strobiformis across latitudinal and elevational gradients. Front. Plant Sci. 2020, 11, 559–697. [Google Scholar] [CrossRef] [PubMed]
  23. Poljak, I.; Idžojtić, M.; Zebec, M.; Perković, N. The variability of European sweet chestnut (Castanea sativa Mill.) in the region of northwest Croatia according to morphology of fruits. Sumar. List 2012, 136, 479–488. [Google Scholar]
  24. Liu, M.S.; Hong, B.G. The distribution of Fagaceae in China and its relationship with climatic and geographic characters. Acta Phytoecol. Sin. 1998, 22, 41–50. [Google Scholar]
  25. Gao, W.Q.; Liu, J.F.; Xue, M.Z.; Zhang, Y.T.; Gao, Z.H.; Ni, Y.Y.; Wang, X.F.; Jiang, Z.P. Geographical patterns and drivers of growth dynamics of Quercus variabilis. For. Ecol. Manag. 2018, 429, 256–266. [Google Scholar] [CrossRef]
  26. Yuan, J.; Sun, N.X.; Du, H.M.; Yin, S.; Kang, H.Z.; Umair, M.; Liu, C.J. Roles of metabolic regulation in developing Quercus variabilis acorns at contrasting geologically-derived phosphorus sites in subtropical China. BMC Plant Biol. 2020, 20, 389. [Google Scholar] [CrossRef]
  27. Shi, W.H.; Villar-Salvador, P.; Li, G.L.; Jiang, X.X. Acorn size is more important than nursery fertilization for outplanting performance of Quercus variabilis container seedlings. Ann. For. Sci. 2019, 76, 1–12. [Google Scholar] [CrossRef] [Green Version]
  28. Sun, X.; Kang, H.; Chen, H.Y.H.; Björn, B.; Samuel, B.F.; Liu, C. Biogeographic patterns of nutrient resorption from Quercus variabilis Blume leaves across China. Plant Biol. 2016, 18, 505–513. [Google Scholar] [CrossRef] [PubMed]
  29. Zhang, Y.W.; Yang, L.N.; Wang, D.M.; Li, D.W. Structure elucidation and properties of different lignins isolated from acorn shell of Quercus variabilis Bl. Int. J. Biol. Macromol. 2018, 107, 1193–1202. [Google Scholar] [CrossRef]
  30. Sun, Q.Y.; Dumroese, R.K.; Liu, Y. Container volume and subirrigation schedule influence Quercus variabilis seedling growth and nutrient status in the nursery and field. Scand. J. For. Res. 2018, 33, 560–567. [Google Scholar] [CrossRef]
  31. Liu, J.F.; Deng, Y.P.; Wang, X.F.; Ni, Y.Y.; Wang, Q.; Xiao, W.F.; Lei, J.P.; Jiang, Z.P.; Li, M.H. The concentration of non-structural carbohydrates, N, and P in Quercus variabilis does not decline toward its northernmost distribution range along a 1500 km transect in China. Front. Plant Sci. 2018, 9, 1444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Llanderal-Mendoza, J.; Gugger, P.F.; Oyama, K.; Uribe-Salas, D.; González-Rodríguez, A. Climatic determinants of acorn size and germination percentage of Quercus rugosa (Fagaceae) along a latitudinal gradient in Mexico. Bot. Sci. 2017, 95, 37–45. [Google Scholar] [CrossRef] [Green Version]
  33. Long, T.J.; Jones, R.H. Seedling growth strategies and seed size effects in fourteen oak species native to different soil moisture habitat. Trees Struct. Funct. 1996, 11, 1–8. [Google Scholar] [CrossRef]
  34. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  35. Slave, C. Using arcmap to create a cadastral database. Case study. J. Inf. Syst. Oper. Manag. 2018, 12, 180–190. [Google Scholar]
  36. Alena, J.; Gösta, E.; Ingegerd, D.; Jan, I. Studies on frost hardiness of Pinus contorta Dougl. seedlings grown in climate chambers. Studia For. Suec. 1981, 157, 4–47. [Google Scholar]
  37. Wagstaff, K.; Cardie, C.; Rogers, S.; Schroedl, S. Constrained K-means clustering with background knowledge. Int. Conf. Mach. Learn. 2001, 1, 577–584. [Google Scholar]
  38. Smith, I.M.; Forbes, A.B. Algorithms and software for fitting polynomial functions constrained to pass through the origin. J. Phys. Conf. Ser. 2018, 1065, 212022. [Google Scholar] [CrossRef]
  39. Lepš, J.; Šmilauer, P. Multivariate analysis of ecological data using canoco 5. Bull. Ecol. Soc. Am. 2014, 85, 5. [Google Scholar]
  40. Li, H.C.; Gan, X.H.; Zhang, Z.P.; Zhang, C.; Song, L. Effects of altitudes and the DBH of seed trees on biological characteristics of Tetracentron sinense (Tetracentraceae) seeds. Plant Divers. Resour. 2015, 37, 177–183. [Google Scholar]
  41. Wright, S. Isolation by distance. Genetics 1943, 28, 139–156. [Google Scholar] [CrossRef] [PubMed]
  42. DeWoody, J.; Trewin, H.; Taylor, G. Genetic and morphological differentiation in Populus nigra L. Isolation by colonization or isolation by adaptation? Mol. Ecol. 2015, 24, 2641–2655. [Google Scholar] [CrossRef] [Green Version]
  43. Poljak, I.; Idžojtić, M.; Šapić, I.; Korijan, P.; Vukelić, J. Diversity and structure of Croatian continental and Alpine-Dinaric populations of grey alder (Alnus incana/L./Moench subsp. incana): Isolation by distance and environment explains phenotypic divergence. Sumar. List 2018, 142, 19–31. [Google Scholar] [CrossRef]
  44. Chen, Y.T.; Wang, S.F.; Chen, Y.C.; Yuan, J.; Wang, T.; Chen, B.; Qian, Y.N. Variation and stability of seed yield, fall-off process, seed size and form in live oak. For. Res. 2015, 28, 524–530. [Google Scholar]
  45. Li, Y.Q.; Li, Y.C.; Wu, Z.Z. Variation in phenotype characters and starch content of Quercus mongolica Fisch seed from different provenances. For. Res. 2013, 26, 528–532. [Google Scholar]
  46. Chang, E.F.; Zhang, Q.; Xiao, G.Y.; Li, P.R.; Li, Y.; Ding, Y.X.; Huang, C.L.; Jing, Y.B. Morphological characteristics and variation analysis of seeds from different provenances and families of Quercus cocciferoides. Seed 2020, 39, 53–58. [Google Scholar]
  47. Rao, G.R.; Shanker, A.K.; Srinivas, I.; Korwar, G.R.; Venkateswarlu, B. Diversity and variability in seed characters and growth of Pongamia pinnata (L.) Pierre accessions. Trees 2011, 25, 725–734. [Google Scholar] [CrossRef]
  48. Fu, X.X.; Liu, H.N.; Zhou, X.D.; Hong, X.T. Morphological variation of Cornus Officinalis seeds in relation to environmental factors. Chin. J. Ecol. 2013, 32, 27–32. [Google Scholar]
  49. Rewicz, A.; Bomanowska, A.; Magda, J.; Rewicz, T. Morphological variability of Consolida regalis seeds of south-eastern and central Europe. Syst. Biodivers. 2016, 15, 25–34. [Google Scholar] [CrossRef]
  50. Ji, M.F.; Zhang, X.W.; Wang, Z.Q.; Zhang, Q.; Deng, J.M. Intra- versus inter-population variation of cone and seed morphological traits of Pinus tabulaeformis Carr in northern China: Impact of climate-related conditions. Pol. J. Ecol. 2011, 59, 717–727. [Google Scholar]
  51. Rodríguez-Gómez, F.; Oyama, K.; Ochoa-Orozco, M.; Mendoza-Cuenca, L.; Gaytán-Legaria, R.; González-Rodríguez, A. Phylogeography and climate-associated morphological variation in the endemic white oak Quercus deserticola (Fagaceae) along the trans-Mexican volcanic belt. Botany 2018, 96, 121–133. [Google Scholar] [CrossRef] [Green Version]
  52. Du, B.M.; Zhu, Y.H.; Kang, H.Z.; Liu, C.J. Spatial variations in stomatal traits and their coordination with leaf traits in Quercus variabilis across eastern Asia. Sci. Total Environ. 2021, 789, 147757. [Google Scholar] [CrossRef] [PubMed]
  53. Zhang, J.B.; Dong, X.; Xin, Z.M.; Liu, M.H.; Zhang, R.H.; Huang, Y.R.; Sun, F. Effects of artificial simulated precipitation on seed characters and germination of Nitraria tangutorum. Southwest China J. Agric. Sci. 2019, 5, 1181–1186. [Google Scholar]
  54. Shu, X.; Yang, Z.L.; Yang, X.; Duan, H.P.; Yu, H.H.; Huang, J.C.; Li, S.C. Variation in seed characters of Magnolia officinalis from different locations. For. Res. 2010, 23, 457–461. [Google Scholar]
  55. Gao, Z.Y.; Zhang, H.F.; Chen, G.P.; Feng, X.M.; Zhao, T.J.; Gao, X.; Shi, F.C. Fruit stone morphology and geographic variation in Juglans mandshurica populations. Chin. J. Appl. Environ. Biol. 2017, 23, 609–615. [Google Scholar]
  56. Scrivanti, L.R.; Mestre, L.; Anton, A.M. Phenotypical variation and taxonomic correlates of five closely related andean species of Poa (Poaceae) along geographic and climatic gradients. Phytotaxa 2014, 183, 121–144. [Google Scholar] [CrossRef]
  57. Malanson, G.P.; Cheney, A.B.; Kinney, M. Climatic and geographic relations of alpine tundra floras in western North America. Alp. Bot. 2015, 125, 21–29. [Google Scholar] [CrossRef]
Diversity 13 00329 g001 550
Figure 1. Relationship between fruit width (FW), fruit length (FL), and FL/FW with geographic and climatic factors in 43 populations of Q. variabilis. (A) relationship between morphological factors with ELAT; (B) relationship between morphological factors with LON; (C) relationship between morphological factors with MAT; (D) relationship between morphological factors with MTW; (E) relationship between morphological factors with AP; (F) relationship between morphological factors with PWQ. AP, annual precipitation; ELAT, equivalent latitude; LON, longitude; MAT, annual mean air temperature; MTW, maximum temperature of warmest month; and PWQ, precipitation of wettest quarter.
Figure 1. Relationship between fruit width (FW), fruit length (FL), and FL/FW with geographic and climatic factors in 43 populations of Q. variabilis. (A) relationship between morphological factors with ELAT; (B) relationship between morphological factors with LON; (C) relationship between morphological factors with MAT; (D) relationship between morphological factors with MTW; (E) relationship between morphological factors with AP; (F) relationship between morphological factors with PWQ. AP, annual precipitation; ELAT, equivalent latitude; LON, longitude; MAT, annual mean air temperature; MTW, maximum temperature of warmest month; and PWQ, precipitation of wettest quarter.
Diversity 13 00329 g001
Diversity 13 00329 g002 550
Figure 2. Principal components analysis (PCA) ordination diagram based on morphological variables.
Figure 2. Principal components analysis (PCA) ordination diagram based on morphological variables.
Diversity 13 00329 g002
Diversity 13 00329 g003 550
Figure 3. Hierarchical tree dendrogram of 43 populations of Q. variabilis based on fruit morphological characters and environmental factors.
Figure 3. Hierarchical tree dendrogram of 43 populations of Q. variabilis based on fruit morphological characters and environmental factors.
Diversity 13 00329 g003
Table
Table 1. Geographical locations and climatic conditions of the 43 study sites of Q. variabilis in China.
Table 1. Geographical locations and climatic conditions of the 43 study sites of Q. variabilis in China.
Site NumberSite LocationCodeLON
(° E)		LAT
(° N)		ALT (m)MAT (°C)MTW (°C)AP (mm)PWQ (mm)ELAT (°)
1Chuzhou City, Anhui ProvinceCZA117.9732.358015.2830.8993944831.25
2Huainan County, Anhui ProvinceHNA117.0032.634915.6431.6489743131.37
3Mentougou District, BeijingMTB116.0939.9621311.7030.2555242239.53
4Miyun County, BeijingMYB117.0740.503579.7528.7651438240.90
5Pinggu District, BeijingPGB117.1340.283539.9028.5953039940.66
6Tianshui City, Gansu ProvinceTSG106.5534.47118910.8126.9164234140.82
7Baise City, Guangxi ProvinceBSG106.5524.7714216.5527.41127068223.98
8Anlong county, Guizhou ProvinceALG104.7024.85169714.6624.19118363634.82
9Xingren County, Guizhou ProvinceXGG104.9525.25129815.9426.16126567632.38
10Dengfeng City, Henan ProvinceDFH113.0534.4537113.4629.7667336534.96
11Jiyuan City, Henan ProvinceJYH112.6035.0715514.6031.9556732634.34
12Lushi County, Henan ProvinceLSH111.0534.0588013.3330.5867133538.20
13Nanzhao County, Henan ProvinceNZH112.4333.4625115.1131.2480438233.21
14Tongbai County, Henan ProvinceTBH113.6832.5317015.1330.8697444431.88
15Badong County, Hubei ProvinceBDH110.3431.0459815.1430.07121654033.17
16Baokang County, Hubei ProvinceBKH111.2631.8868013.7229.35105847234.59
17Enshi City, Hubei ProvinceESH109.4930.2849116.2631.48146865231.65
18Jianshi County, Hubei ProvinceJSH109.7330.6073014.8429.81138360033.67
19Jingshan City, Hubei ProvinceJSA113.1231.0210316.0831.17107146730.03
20Suixian County, Hubei ProvinceSXH112.9831.5328715.0530.20103344731.46
21Jiangui County, Hubei ProvinceJGH110.9830.8361015.6230.49116753633.04
22Wuhan City, Hubei ProvinceWHH114.3130.592717.2633.00126556129.23
23Zhushan County, Hubei ProvinceZSH110.2332.2241815.2831.21100445433.07
24Xiangxi Prefecture, Hunan ProvinceXPH109.7428.3127717.3332.57133960428.20
25Changsha City, Hunan ProvinceCSH112.9428.236117.7833.33140359727.03
26Nanchang City, Jiangxi ProvinceNCJ115.8328.763717.7033.42155672627.44
27Dalian City, Liaoning ProvinceDLL121.7939.1013710.2226.6964640238.28
28Liaoyang City, Liaoning ProvinceLYL123.3041.081718.0427.2574247340.43
29Yantai City, Shandong ProvinceYTS121.7437.2622211.1326.3972143436.87
30Xia County, Shanxi ProvinceXCS111.3735.0111859.9526.4961133341.33
31Baoji City, Shaanxi ProvinceBJS107.1434.3768013.1130.2568135837.08
32Shanyang County, Shaanxi ProvinceSYS109.8833.5372612.9329.0377837536.58
33Shangnan County, Shaanxi ProvinceSNS110.8833.5382614.6931.0577437337.29
34Weinan County, Shaanxi ProvinceWNS109.5034.5053613.1630.9460629736.19
35Xianyang City, Shaanxi ProvinceXYS108.0834.2748612.7629.7864832635.60
36Linzhi City, TibetLZT94.3629.6531647.9020.1665135950.11
37Urumqi, XinjiangUQX87.6243.838997.0630.222318748.10
38Anning City, Yunnan ProvinceANY102.4524.99185215.2524.8189849636.07
39Honghe Prefecture, Yunnan ProvinceHHY103.6123.33165514.7122.98136776133.01
40Kunming, Yunnan ProvinceKMY102.7525.14205114.3023.7092150937.65
41Zhanyi County, Yunnan ProvinceZYY103.5525.59221413.8423.2393851539.27
42Lin’an District, Zhejiang ProvinceLAZ119.4430.3332014.7830.13139955430.47
43Yuhang District, Zhejiang ProvinceYHZ120.3030.42916.5132.29126247228.96
ELAT, equivalent latitude; LON, longitude; LAT, latitude; ALT, altitude; MAT, annual mean air temperature; MTW, maximum temperature of warmest month; AP, annual precipitation; and PWQ, precipitation of wettest quarter. The climate data for all sites were obtained from the Global Climate and Weather Data website (http://www.worldclim.org (accessed on 15 December 2020)).
Table
Table 2. Statistical data of fruit length (FL), fruit width (FW), and FL/FW of 43 populations of Q. variabilis.
Table 2. Statistical data of fruit length (FL), fruit width (FW), and FL/FW of 43 populations of Q. variabilis.
SiteFLFWFL/FWSample Size
MeanMaxMinSDCV (%)MeanMaxMinSDCV (%)MeanMaxMinSDCV (%)
119.51 defghijk23.5615.021.718.7418.21 ijk26.7813.131.588.701.07 bcde1.580.700.109.08100
220.42 hijkl24.2215.881.135.5319.99 lmn27.1713.821.507.511.04 abc1.320.760.066.21120
319.20 cdefghij22.9415.491.497.7914.82 bc18.4610.871.268.501.30 i1.601.050.075.5680
420.52 ijkl22.7518.680.974.7418.88 klm21.9016.331.437.561.09 bcde1.360.940.109.2620
520.40 hijkl24.2716.691.899.2617.42 fghijk20.9715.001.418.091.18 efgh1.500.860.1411.5520
617.81 abc21.5213.362.1211.8815.75 cde19.5212.441.7711.261.14 cdefg1.400.940.1210.5520
720.31 hijkl24.4914.142.1310.4816.53 defghi20.6913.751.579.521.23 fghi1.521.000.1310.4630
819.53 defghijk22.8013.451.929.8218.82 klm24.5614.262.0710.991.04 abc1.240.850.109.3730
921.15 klmn23.8418.201.426.7120.26 mn22.8917.971.025.031.04 abc1.230.900.076.8930
1019.49 defghijk26.8013.971.608.2317.42 fghijk22.0812.571.488.481.12 bcdef1.610.820.108.72438
1121.69 lmn26.1017.081.617.4417.68 fghijk21.5713.971.236.941.23 fghi1.500.980.096.9460
1219.89 efghijk20.7418.920.462.2917.07 efghij18.8314.471.277.411.17 defgh1.361.000.097.8520
1318.65 abcdefg23.0013.901.638.7616.20 cdef21.5010.451.539.431.16 cdefgh1.760.850.119.2580
1420.36 hijkl24.4515.701.296.3520.20 mn23.4216.61.497.351.01 ab1.200.860.055.1560
1517.17 a24.35 9.961.438.3314.95 bc20.089.781.268.441.14 cdefg1.440.800.1210.8260
1619.71 fghijkl22.8214.621.326.7017.82 ghijk20.4014.261.055.891.11 cdefg1.410.900.098.5540
1717.22 a20.7911.561.649.5411.70 a13.799.280.917.741.48 j1.991.030.1610.6260
1820.27 ghijkl23.8917.581.919.4118.19 ijk21.9915.311.669.141.12 bcde1.410.990.108.8730
1920.05 fghijkl23.5817.021.055.2517.29 efghijk21.6313.551.166.691.17 defgh1.620.940.097.3160
2020.08 fghijkl22.0216.651.216.0420.93 n23.1318.491.125.360.96 a1.120.820.077.2530
2118.90 bcdefghi22.8214.621.9510.3417.72 fghijk19.8815.301.196.701.07 abcde1.360.900.1312.2220
2219.51 cdefghijk24.4913.791.286.5717.73 fghijk23.3312.721.206.751.10 bcde1.450.810.076.7899
2320.95 jklm22.7319.190.904.3019.70 lmn21.7217.221.065.381.06 abcde1.130.960.054.4920
2419.99 efghijk23.9315.121.507.5213.58 b18.209.461.349.841.50 j2.011.080.138.4749
2519.58 defghijk26.7912.871.356.9011.89 a20.177.960.857.111.68 k2.251.010.148.26180
2621.12 klmn23.6117.461.517.1419.99 lmn22.3416.791.547.681.06 abcde1.220.950.065.9520
2719.93 efghijk22.4516.911.497.4916.16 cdefg18.2112.741.348.271.24 ghi1.440.960.118.6220
2819.00 bcdefghi24.7314.302.0610.8516.96 defghij21.3213.192.0812.271.13 bcdefg1.420.930.1210.8530
2919.79 efghijk26.6915.661.698.5515.73 cde19.0312.221.278.081.26 hi1.861.040.118.8280
3018.42 abcdef20.9415.681.427.7117.34 efghijk21.4214.321.538.811.07 abcde1.220.910.076.6520
3122.46 n24.8620.481.295.7415.45 cd17.2313.681.066.861.46 j1.581.340.064.3820
3218.84 abcde21.7814.621.437.5617.54 efghijk21.4212.071.448.191.08 abcde1.350.840.087.6960
3319.83 efghijk21.7817.021.015.1217.88 ghijk21.4214.161.588.811.11 bcde1.300.910.087.5520
3418.04 abcd22.6214.081.789.8917.06 efghij21.3714.091.468.541.06 abcde1.280.810.1110.0960
3520.55 ijkl24.0416.241.919.3017.93 hijk20.8813.721.538.541.16 cdefg1.440.860.1513.2530
3618.63 abcdefgh22.5013.841.9210.3217.78 fghijk21.3413.451.8710.511.05 abcd1.250.860.098.7630
3731.69 o35.0424.652.437.6717.86 ghijk20.9714.531.769.861.78 l2.311.560.168.9730
3818.46 a20.6415.271.367.3817.29 defgh19.4815.521.046.021.07 abc1.250.890.066.0250
3920.67 ijkl23.3318.000.934.5218.69 jkl21.0915.771.407.501.11 bcde1.340.960.109.1630
4017.54 ab20.4211.781.8010.2916.65 defghi20.9112.281.519.091.06 abcde1.240.810.098.9347
4122.16 mn25.4817.781.506.7718.95 klm20.8315.671.166.141.17 defgh1.421.030.097.8530
4219.62 defghijk24.7014.762.0210.2917.69 fghijk21.6115.531.448.151.11 bcde1.260.910.087.5430
4319.48 defghijk21.7916.011.306.6818.14 hijk21.4613.251.749.601.08 bcde1.420.910.1110.1630
Total19.9723.6515.771.537.7217.3521.0913.771.408.111.171.460.940.108.41-
CV betweenPopulations(%)10.8011.2114.48
Different letters in a column indicate significant differences at p ≤ 0.05. Refer to Table 1 for details of site locations.
Table
Table 3. ANOVA of fruit width, fruit length, and fruit length-to-width ratio of 43 populations of Q. variabilis.
Table 3. ANOVA of fruit width, fruit length, and fruit length-to-width ratio of 43 populations of Q. variabilis.
Variance SourcedfSSMSFp
Fruit LengthInter-group718342171.01539.28<0.001
Intra-group10,23223504.354
Total17,4152392
Fruit WidthInter-group11,45242272.66767.70<0.001
Intra-group946423504.027
Total20,9162392
Fruit Length-to-Width RatioInter-group84422.009105.42<0.001
Intra-group4523500.019
Total1292392
Table
Table 4. Correlation between geoclimatic factors and phenotypic traits of Q. variabilis populations.
Table 4. Correlation between geoclimatic factors and phenotypic traits of Q. variabilis populations.
LON (°)ELAT (°)ALT (m)MAT (°C)MTW (°C)AP (mm)PWQ (mm)
FL−0.402 **0.243−0.040−0.2610.103−0.282−0.333 *
FW−0.0380.0420.126−0.019−0.155−0.055−0.042
FL/FW−0.2440.071−0.174−0.1030.263−0.088−0.150
AP, annual precipitation; ELAT, equivalent latitude; LON, longitude; MAT, annual mean air temperature; MTW, maximum temperature of warmest month; PWQ, precipitation of wettest quarter; and ALT, altitude. ** p ≤ 0.01, * p ≤ 0.05.
Table
Table 5. Correlations between morphological (FL, FW), climatic (MAT, MTW, AP, and PWQ) and geographic (LON, LAT, ALT, and ELAT) distance matrices.
Table 5. Correlations between morphological (FL, FW), climatic (MAT, MTW, AP, and PWQ) and geographic (LON, LAT, ALT, and ELAT) distance matrices.
Comparisonrp-Value
Morphological, Geographic0.633<0.001
Morphological, Climate0.988<0.001
Table
Table 6. Principal component analysis loadings and percentage variance contributions by principal components.
Table 6. Principal component analysis loadings and percentage variance contributions by principal components.
Principal
Component		PC1PC2
FL0.380.93
FW−0.780.62
FL/FW0.990.14
Eigenvalue1.731.26
Variance (%)57.5742.03
% Total Variance57.5799.60
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://0-doi-org.brum.beds.ac.uk/10.3390/d13070329

AMA Style

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(7):329. https://0-doi-org.brum.beds.ac.uk/10.3390/d13070329

Chicago/Turabian Style

Gao, Shuang, Yue Ren, Joseph Masabni, Feng Zou, Huan Xiong, and Jingle Zhu. 2021. "Influence of Geographical and Climatic Factors on Quercus variabilis Blume Fruit Phenotypic Diversity" Diversity 13, no. 7: 329. https://0-doi-org.brum.beds.ac.uk/10.3390/d13070329

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