Landscape Features and Climatic Forces Shape the Genetic Structure and Evolutionary History of an Oak Species (Quercus chenii) in East China

Population differentiation and dynamics of five pioneer species of Gaultheria from the secondary forests in subtropical China

BACKGROUND

The influence of native secondary succession associated with anthropogenic disturbance on the biodiversity of the forests in subtropical China remains uncertain. In particular, the evolutionary response of small understory shrubs, particularly pioneer species inhabiting continuously disturbed habitats, to topographic heterogeneity and climate change is poorly understood. This study aimed to address this knowledge gap by focusing on the Gaultheria crenulata group, a clade of small pioneer shrubs in subtropical China.

RESULTS

We examined the genetic structure and demographic history of all five species of the G. crenulata group with two maternally inherited chloroplast DNA (cpDNA) fragments and two biparentally inherited low-copy nuclear genes (LCG) over 89 natural populations. We found that the genetic differentiation of this group was influenced by the geomorphological boundary between different regions of China in association with Quaternary climatic events. Despite low overall genetic diversity, we observed an isolation-by-distance (IBD) pattern at a regional scale, rather than isolation-by-environment (IBE), which was attributed to ongoing human disturbance in the region.

CONCLUSION

Our findings suggest that the genetic structure of the G. crenulata group reflects the interplay of geological topography, historical climates, and anthropogenic disturbance during the Pliocene-Pleistocene-Holocene periods in subtropical China. The observed IBD pattern, particularly prominent in western China, highlights the role of limited dispersal and gene flow, possibly influenced by physical barriers or decreased connectivity over geographic distance. Furthermore, the east-to-west trend of gene flow, potentially facilitated by the East Asian monsoon system, underscores the complex interplay of biotic and abiotic factors shaping the genetic dynamics of pioneer species in subtropical China's secondary forests. These findings can be used to assess the impact of environmental changes on the adaptation and persistence of biodiversity in subtropical forest ecosystems.

Ensemble species distribution modeling and multilocus phylogeography provide insight into the spatial genetic patterns and distribution dynamics of a keystone forest species, Quercus glauca

BACKGROUND

Forests are essential for maintaining species diversity, stabilizing local and global climate, and providing ecosystem services. Exploring the impact of paleogeographic events and climate change on the genetic structure and distribution dynamics of forest keystone species could help predict responses to future climate change. In this study, we combined an ensemble species distribution model (eSDM) and multilocus phylogeography to investigate the spatial genetic patterns and distribution change of Quercus glauca Thunb, a keystone of East Asian subtropical evergreen broad-leaved forest.

RESULTS

A total of 781 samples were collected from 77 populations, largely covering the natural distribution of Q. glauca. The eSDM showed that the suitable habitat experienced a significant expansion after the last glacial maximum (LGM) but will recede in the future under a general climate warming scenario. The distribution centroid will migrate toward the northeast as the climate warms. Using nuclear SSR data, two distinct lineages split between east and west were detected. Within-group genetic differentiation was higher in the West than in the East. Based on the identified 58 haplotypes, no clear phylogeographic structure was found. Populations in the Nanling Mountains, Wuyi Mountains, and the southwest region were found to have high genetic diversity.

CONCLUSIONS

A significant negative correlation between habitat stability and heterozygosity might be explained by the mixing of different lineages in the expansion region after LGM and/or hybridization between Q. glauca and closely related species. The Nanling Mountains may be important for organisms as a dispersal corridor in the west-east direction and as a refugium during the glacial period. This study provided new insights into spatial genetic patterns and distribution dynamics of Q. glauca.

Scale-dependent effects of species diversity on aboveground biomass and productivity in a subtropical broadleaved forest on Mt. Huangshan

The relationship between species diversity and biomass/productivity is a major scientific question in ecology. Exploring this relationship is essential to understanding the mechanisms underpinning the maintenance of biodiversity. Positive, negative, and neutral relationships have been identified in controlled experiments and observational research. However, increasing evidence suggests that the effects of species diversity on aboveground biomass and productivity are influenced by biotic and abiotic factors, but it remains unclear whether scale-dependent effects affect aboveground biomass and productivity. Herein, we used a generalized linear regression model and a structural equation model to explore relationships between species diversity and productivity/aboveground biomass under different scales and to investigate the effects of topographical factors and species diversity on ecosystem functioning. The results revealed a positive relationship between biodiversity and ecosystem functioning based on species diversity and aboveground biomass. Different sampling scales may impact the relationship between species diversity and ecosystem functioning. A positive relationship was found between species richness and productivity at medium and large scales; however, ambiguous relationships were found in productivity and other species diversity indices. Elevation was a key factor affecting both biomass and productivity. These results suggest that species diversity is not the only factor affecting biomass and productivity, and the positive correlation between species diversity and ecosystem functioning is mediated by abiotic factors.

New insights on patterns of genetic admixture and phylogeographic history in Iberian high mountain populations of midwife toads

Multiple Quaternary glacial refugia in the Iberian Peninsula, commonly known as "refugia within refugia", allowed diverging populations to come into contact and admix, potentially boosting substantial mito-nuclear discordances. In this study, we employ a comprehensive set of mitochondrial and nuclear markers to shed light onto the drivers of geographical differentiation in Iberian high mountain populations of the midwife toads Alytes obstetricans and A. almogavarii from the Pyrenees, Picos de Europa and Guadarrama Mountains. In the three analysed mountain regions, we detected evidence of extensive mito-nuclear discordances and/or admixture between taxa. Clustering analyses identified three major divergent lineages in the Pyrenees (corresponding to the eastern, central and central-western Pyrenees), which possibly recurrently expanded and admixed during the succession of glacial-interglacial periods that characterised the Late Pleistocene, and that currently follow a ring-shaped diversification pattern. On the other hand, populations from the Picos de Europa mountains (NW Iberian Peninsula) showed a mitochondrial affinity to central-western Pyrenean populations and a nuclear affinity to populations from the central Iberian Peninsula, suggesting a likely admixed origin for Picos de Europa populations. Finally, populations from the Guadarrama Mountain Range (central Iberian Peninsula) were depleted of genetic diversity, possibly as a consequence of a recent epidemic of chytridiomycosis. This work highlights the complex evolutionary history that shaped the current genetic composition of high mountain populations, and underscores the importance of using a multilocus approach to better infer the dynamics of population divergence.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, “dehesa.” It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors’ group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

Extensive sharing of chloroplast haplotypes among East Asian Cerris oaks: The imprints of shared ancestral polymorphism and introgression

Shared ancestral polymorphism and introgression are two main causes of chloroplast DNA (cpDNA) haplotype sharing among closely related angiosperms. In this study, we explored the roles of these two processes in shaping the phylogeographic patterns of East Asian Cerris oaks by examining the geographic distributions of randomly and locally distributed shared haplotypes, which coincide with the expectations of shared ancestry and introgression, respectively. We sequenced 1340 bp of non‐coding cpDNA from Quercus acutissima (n = 418) and Q. chenii (n = 183) and compiled previously published sequence data of Q. variabilis (n = 439). The phylogenetic relationships among haplotypes were examined using a median‐joining network. The geographic patterns of interspecifically shared haplotypes were assessed to test whether nearby populations have a higher degree of interspecific cpDNA sharing than distant ones. We identified a total of 27 haplotypes that were grouped into three non‐species‐specific lineages with overlapping distributions. Ancestral haplotypes were extensively shared and randomly distributed across populations of the three species. Some young haplotypes were locally shared in mountainous areas that may have been shared refugia. The local exchange of cpDNA resulted in an excess of similar haplotypes between nearby populations. Our study demonstrated that the haplotype sharing pattern among East Asian Cerris oaks reflected the imprints of both shared ancestral polymorphism and introgression. This pattern was also associated with the relatively stable climates and complex landscapes in East Asia, which not only allowed the long‐term persistence of ancestral lineages but also connected the survived populations across refugia.

Stepped Geomorphology Shaped the Phylogeographic Structure of a Widespread Tree Species (Toxicodendron vernicifluum, Anacardiaceae) in East Asia

Species' phylogeographic patterns reflect the interplay between landscape features, climatic forces, and evolutionary processes. Here, we used two chloroplast DNA (cpDNA) markers (trnL and trnL-F) to explore the role of stepped geomorphology in shaping the phylogeographic structure of Toxicodendron vernicifluum, an economically important tree species widely distributed in East Asia. The range-wide pattern of sequence variation was analyzed based on a dataset including 357 individuals from China, together with published sequences of 92 individuals mainly from Japan and South Korea. We identified five chloroplast haplotypes based on seven substitutions across the 717-bp alignment. A clear east-west phylogeographic break was recovered according to the stepped landforms of mainland China. The wild trees of the western clade were found to be geographically restricted to the "middle step", which is characterized by high mountains and plateaus, while those of the eastern clade were confined to the "low step", which is mainly made up of hills and plains. The two major clades were estimated to have diverged during the Early Pleistocene, suggesting that the cool glacial climate may have caused the ancestral population to retreat to at least two glacial refugia, leading to allopatric divergence in response to long-term geographic isolation. Migration vector analyses based on the outputs of ecological niche models (ENMs) supported a gradual range expansion since the Last Interglacial. Mountain ranges in western China and the East China Sea land bridge were inferred to be dispersal corridors in the western and eastern distributions of T. vernicifluum, respectively. Overall, our study provides solid evidence for the role of stepped geomorphology in shaping the phylogeographic patterns of T. vernicifluum. The resulting east-west genetic discontinuities could persist for a long time, and could occur at a much larger scale than previously reported, extending from subtropical (e.g., the Xuefeng Mountain) to warm-temperate China (e.g., the Taihang Mountain).

Determinants of Genetic Structure in a Highly Heterogeneous Landscape in Southwest China

Intra-specific genetic diversity is a fundamental component of biodiversity, and is key to species adaptation and persistence. However, significant knowledge gaps still exist in our understanding of the patterns of genetic diversity and their key determinants. Most previous investigations mainly utilized single-species and/or a limited number of explanatory variables; so here we mapped the patterns of plastid genetic diversity within 15 plant species, and explored the key determinants shaping these patterns using a wide range of variables. Population-level cpDNA sequence data for 15 plant species from the Longitudinal Range Gorge Region (LRGR), southwest China, were retrieved from literature and used to estimate haplotype diversity (H _D) and population pairwise genetic differentiation (F _ST) indices. Genetic diversity and divergence landscape surfaces were then generated based on the H _D and F _ST, respectively, to clarify the patterns of genetic structure in the region. Subsequently, we analyzed the relationships between plastid genetic diversity and 16 explanatory variables (classified as anthropogenic, climatic, and topographic). We found that the highest genetic diversity occurred in the Yulong Mountain region, with a significant proportion (~74.81%) of the high diversity land area being located outside of protected areas. The highest genetic divergence was observed approximately along the 25°N latitudinal line, with notable peaks in the western and eastern edges of the LRGR. Genetic diversity (H _D) was weakly but significantly positively correlated with both Latitude (lat) and Annual Mean Wet Day Frequency (wet), yet significantly negatively correlated with all of Longitude (long), Annual Mean Cloud Cover Percent (cld), Annual Mean Anthropogenic Flux (ahf), and Human Footprint Index (hfp). A combination of climatic, topographic, and anthropogenic factors explained a significant proportion (78%) of genetic variation, with topographic factors (lat and long) being the best predictors. Our analysis identified areas of high genetic diversity (genetic diversity "hotspots") and divergence in the region, and these should be prioritized for conservation. This study contributes to a better understanding of the features that shape the distribution of plastid genetic diversity in the LRGR and thus would inform conservation management efforts in this species-rich, but vulnerable region.

Combining the Effects of Global Warming, Land Use Change and Dispersal Limitations to Predict the Future Distributions of East Asian Cerris Oaks (Quercus Section Cerris, Fagaceae) in China

Species shift their ranges in response to climate change (CC). However, they may not be able to track optimal conditions as soon as possible, due to limited dispersal ability or habitat fragmentation, caused by land use and land cover change (LULC). This study aimed to explore the combined impacts of CC, LULC and dispersal limitations on the future range dynamics of Quercus acutissima Carruth., Q. variabilis Blume and Q. chenii Nakai, three dominant Cerris oak tree species in warm-temperate and subtropical deciduous forests of China. We used the Maximum Entropy (Maxent) algorithm to predict the suitable habitats for the years 2050 and 2070, under three representative concentration pathways (RCPs). Habitat fragmentation patterns were examined to assess the influence of LULC. Two migration scenarios (full- and partial-migration) were compared to evaluate the effect of dispersal limitations. We found that annual precipitation (AP), minimum temperature in the coldest month (MTCM) and temperature seasonality (TS) play a key role in determining the present distributions of Q. chenii, while AP, MTCM and annual mean temperature (AMT) contribute the most to the distribution models of Q. variabilis and Q. acutissima. For all the three species, LULC will increase the level of habitat fragmentation and lead to the loss of core areas, while limited dispersal ability will restrict the accessibility of future potentially suitable habitats. Under the scenarios of CC and LULC, the suitable areas of Q. chenii will decrease sharply, while those of Q. variabilis in South China will become unsuitable. Our findings highlight the importance of considering dispersal ability, as well as land use and land cover change, for modeling species’ range shifts in the face of global warming. Our study also provides vital information for guiding the management of East Asian Cerris oaks in China; Q. chenii should be listed as a species requiring priority protection, and the threatened habitats of Q. variabilis should be protected to buffer the impacts of CC and LULC.

Genetic Evaluation of Juniperus sabina L. (Cupressaceae) in Arid and Semi-Arid Regions of China Based on SSR Markers

Juniperus sabina L., a shrub distributed in patches in arid and semi-arid areas of the northern hemisphere, plays an important role in preventing land desertification and maintaining ecosystems. However, few studies have reported genetic diversity and genetic structure of widely distributed populations of J. sabina in northwest China. Here, we evaluated the genetic diversity and genetic structure and predicted the isolation barriers among 11 populations based on 20 simple sequence repeats (SSRs). A total of 134 alleles were generated and the average number of alleles per locus was 6.70. The Shannon diversity index ranged from 0.659 to 0.951, with an average of 0.825. Population structure analysis revealed that the populations were assigned into two genetic groups. The analysis of molecular variance (AMOVA) indicated that 88% of genetic variation existed within populations. Moderate population differentiation was occurred with FST value of 0.090. Finally, we concluded that geographic isolation is the main factor affecting the genetic structure of J. sabina populations. The results of this study provide a foundation for the strategies for J. sabina genetic conservation and management.

Genetic variation of the relict maple Acer miyabei: uncovering its history of disjunct occurrence and the role of mountain refugia in shaping genetic diversity

PREMISE

Relict species provide valuable insights into the origin and formation of extant vegetation. Here, we aimed to elucidate the genetic structure and diversity of a riparian relic, Acer miyabei, in Japan. Once widely distributed, it now occurs in three isolated regions. The most northern regional group is located at low elevation on Hokkaido Island, whereas the southernmost group in central Honshu Island is at high elevation in a mountainous landscape. This contrastive distribution enables us to examine the effects of climate oscillations on genetic diversity in relation to topographic variation.

METHODS

We collected 604 individuals of A. miyabei from 43 sites. Their genetic structure and diversity were analyzed using 12 microsatellite markers and cpDNA sequences.

RESULTS

According to structure analyses, ∆K was lowest at K = 2; the clustering essentially separated many of the individuals in the most northern regional group from the others. In contrast, the two southern groups were not clearly differentiated from each other, despite their geographic discontinuity. The proportion of private alleles was high in populations from the mountain terrain in the southern group although the number of extant populations is limited.

CONCLUSIONS

Genetic clustering of A. miyabei is not perfectly congruent with the current patterns of geographic distribution. We infer that disjunction of the two southern groups occurred more recently than that between these groups and the northern group. The mountainous landscape in the most southern region likely provided multiple refugia and contributed to the retention of distinctive genetic variation.

Mismatch Between Specific and Genetic Diversity in an Evergreen Broadleaf Forest in Southeast China: A Study Case of 10.24 ha Forest Dynamics Plot of Huangshan

For a long time, forestry management has often focused on the protection of species diversity, and mistakenly believed that protecting species diversity protects genetic diversity. Therefore, research that integrates community ecology and population genetics has become important because it can help elucidate whether the targets for protecting specific and genetic diversity are congruent. In this study, we have emphasized the impact of the community on the population because no previous studies have considered the community composition of a place a priori. Based on the Huangshan 10.24 ha dynamics forest plot, we a priori considered the community composition in the plot to test species-genetic diversity among the tree layers. Firstly, a redundancy analysis (RDA) found that Castanopsis eyrei and Pinus massoniana were the dominant species. Secondly, specific and genetic diversity are not congruent in Huang Shan. Finally, the structural equation model (SEM) showed that the different degrees of response by community composition and population structure to environmental heterogeneity are the main reasons for the mismatch between species diversity and genetic diversity. The results suggest that we must focus on genetic diversity, as well as on protecting species diversity.

Influence of Pliocene and Pleistocene climates on hybridization patterns between two closely related oak species in China

BACKGROUND AND AIMS

Contemporary patterns of genetic admixture reflect imprints of both ancient and recent gene flow, which can provide us with valuable information on hybridization history in response to palaeoclimate change. Here, we examine the relationships between present admixture patterns and past climatic niche suitability of two East Asian Cerris oaks (Quercus acutissima and Q. chenii) to test the hypothesis that the mid-Pliocene warm climate promoted while the Pleistocene cool climate limited hybridization among local closely related taxa.

METHODS

We analyse genetic variation at seven nuclear microsatellites (1111 individuals) and three chloroplast intergenic spacers (576 individuals) to determine the present admixture pattern and ancient hybridization history. We apply an information-theoretic model selection approach to explore the associations of genetic admixture degree with past climatic niche suitability at multiple spatial scales.

KEY RESULTS

More than 70 % of the hybrids determined by Bayesian clustering analysis and more than 90 % of the individuals with locally shared chloroplast haplotypes are concentrated within a mid-Pliocene contact zone between ~30°N and 35°N. Climatic niche suitabilities for Q. chenii during the mid-Pliocene Warm Period [mPWP, ~3.264-3.025 million years ago (mya)] and during the Last Glacial Maximum (LGM, ~0.022 mya) best explain the admixture patterns across all Q. acutissima populations and across those within the ancient contact zone, respectively.

CONCLUSIONS

Our results highlight that palaeoclimate change shapes present admixture patterns by influencing the extent of historical range overlap. Specifically, the mid-Pliocene warm climate promoted ancient contact, allowing widespread hybridization throughout central China. In contrast, the Pleistocene cool climate caused the local extinction of Q. chenii, reducing the probability of interspecific gene flow in most areas except those sites having a high level of ecological stability.

Linked selection shapes the landscape of genomic variation in three oak species

Natural selection shapes genome-wide patterns of diversity within species and divergence between species. However, quantifying the efficacy of selection and elucidating the relative importance of different types of selection in shaping genomic variation remain challenging. We sequenced whole genomes of 101 individuals of three closely related oak species to track the divergence history, and to dissect the impacts of selective sweeps and background selection on patterns of genomic variation. We estimated that the three species diverged around the late Neogene and experienced a bottleneck during the Pleistocene. We detected genomic regions with elevated relative differentiation ('F_ST -islands'). Population genetic inferences from the site frequency spectrum and ancestral recombination graph indicated that F_ST -islands were formed by selective sweeps. We also found extensive positive selection; the fixation of adaptive mutations and reduction neutral diversity around substitutions generated a signature of selective sweeps. Prevalent negative selection and background selection have reduced genetic diversity in both genic and intergenic regions, and contributed substantially to the baseline variation in genetic diversity. Our results demonstrate the importance of linked selection in shaping genomic variation, and illustrate how the extent and strength of different selection models vary across the genome.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

A Meta-Analysis Indicates Positive Correlation between Genetic Diversity and Species Diversity

Simple Summary

Understanding species and genetic correlations (SGDCs) is essential to establish community composition. In this study, 295 observations from 39 studies explored the SGDCs and the underlying drivers through conducting a global meta-analysis. A positive correlation was found, suggesting that parallel processes (environmental heterogeneity, area, and connectivity etc.) have effects on two diversities. As current biodiversity hotspots have mainly been identified based on high species diversity and high endemism of taxon, the understanding of SGDC will substantially help us to determine whether and how genetic diversity can be used in identifying biodiversity hotspots, as well as in developing conservation practices and policies for biodiversity.

Abstract

Species diversity (SD) and genetic diversity (GD) are the two basic levels of biodiversity. In general, according to the consensus view, the parallel effects of environmental heterogeneity, area, and connectivity on two levels, can drive a positive correlation between GD and SD. Conversely, a negative correlation or no correlation would be expected if these effects are not parallel. Our understanding of the relationships between SD and GD among different ecosystems, sampling methods, species, and under climate change remains incomplete. In the present study, we conducted a hierarchical meta-analysis based on 295 observations from 39 studies and found a positive correlation between genetic diversity and species diversity (95% confidence interval, 7.6–22.64%). However, significant relationships were not found in some ecosystems when we conducted species–genetic diversity correlation analysis based on a single ecosystem. Moreover, the magnitudes of the correlations generally decreased with the number of sampling units and the annual average the temperature of sampling units. Our results highlight the positive correlation between GD and SD, thereby indicating that protecting SD involves protecting GD in conservation practice. Furthermore, our results also suggest that global increases in temperature during the 21st century will have significant impacts on global biodiversity.

Genetic Differentiation and Demographic History of Three Cerris Oak Species in China Based on Nuclear Microsatellite Makers

Knowledge of interspecific divergence and population expansions/contractions of dominant forest trees in response to geological events and climatic oscillations is of major importance to understand their evolution and demography. However, the interspecific patterns of genetic differentiation and spatiotemporal population dynamics of three deciduous Cerris oak species (Q. acutissima, Q. variabilis and Q. chenii) that are widely distributed in China remain poorly understood. In this study, we genotyped 16 nuclear loci in 759 individuals sampled from 44 natural populations of these three sibling species to evaluate the plausible demographical scenarios of the closely related species. We also tested the hypothesis that macro- and microevolutionary processes of the three species had been triggered and molded by Miocene–Pliocene geological events and Quaternary climatic change. The Bayesian cluster analysis showed that Q. acutissima and Q. chenii were clustered in the same group, whereas Q. variabilis formed a different genetic cluster. An approximate Bayesian computation (ABC) analyses suggested that Q. variabilis and Q. acutissima diverged from their most common ancestor around 19.84 Ma, and subsequently Q. chenii diverged from Q. acutissima at about 9.6 Ma, which was significantly associated with the episodes of the Qinghai–Tibetan Plateau (QTP). In addition, ecological niche modeling and population history analysis showed that these three Cerris oak species repeatedly underwent considerable ‘expansion–contraction’ during the interglacial and glacial periods of the Pleistocene, although they have varying degrees of tolerance for the climatic change. Overall, these findings indicated geological and climatic changes during the Miocene–Pliocene and Pleistocene as causes of species divergence and range shifts of dominant tree species in the subtropical and warm temperature areas in China.

Spatial Genetic Structure and Demographic History of the Dominant Forest Oak Quercus fabri Hance in Subtropical China

Oak trees (Quercus L.) are important models for estimating abiotic impacts on the population structure and demography of long life span tree species. In this study, we generated genetic data for 17 nuclear microsatellite loci in 29 natural populations of Quercus fabri to estimate the population genetic structure. We also integrated approximate Bayesian computation (ABC) and ecological niche analysis to infer the population differentiation processes and demographic history of this oak species. The genetic analyses indicated two genetic clusters across the 29 populations collected, where most approximately corresponded to the intraspecific differentiation among populations from western and eastern China, whereas admixed populations were mainly found in central mountains of China. The best model obtained from hierarchical ABC simulations suggested that the initial intraspecific divergence of Q. fabri potentially occurred during the late Pliocene (ca. 3.99 Ma) to form the two genetic clusters, and the admixed population group might have been generated by genetic admixture of the two differentiated groups at ca. 53.76 ka. Ecological analyses demonstrated clear differentiation among the Q. fabri population structures, and association estimations also indicated significant correlations between geography and climate with the genetic variation in this oak species. Our results suggest abiotic influences, including past climatic changes and ecological factors, might have affected the genetic differentiation and demographic history of Q. fabri in subtropical China.

Distribution Pattern of Endangered Plant Semiliquidambar cathayensis (Hamamelidaceae) in Response to Climate Change after the Last Interglacial Period

Semiliquidambar cathayensis is a special and endangered plant in China, used for traditional Chinese medicine and in landscape applications. Predicting the impact of climate change on the distribution of S. cathayensis is crucial for its protection and the sustainable use of resources. We used the maximum entropy (MaxEnt) model optimized by the ENMeval data packet to analyze the potential geographic distribution changes of S. cathayensis in 12 provinces of Southern China for the different periods since the last interglacial period (LIG, 120–140 ka). Considering the potential geographic distribution changes in the province, and based on the two climate scenarios of Representative Concentration Pathways (RCP) 2.6 and RCP 8.5, the distribution range of S. cathayensis was analyzed and we predicted the range for the 2050s (average for 2041–2060) and 2070s (average for 2061–2080). The area under AUC (Area under the receiver operating characteristic (ROC) curve) is 0.9388 under these parameters, which indicates that the model is very accurate. We speculate that the glacial period refugia were the Nanling and Wuyi Mountains for S. cathayensis, and central and Western Fujian and Taiwan are likely to be the future climate refugia. In the mid-Holocene (MH, 6 ka), the growth habitat was 32.41% larger than the modern habitat; in the 2050s and 2070s (except RCP2.6–2070s), the growth habitat will shrink to varying degrees, so efforts to support its in situ and ex situ conservation are urgently needed. The jackknife test showed that the main factors affecting the geographical distribution of S. cathayensis were annual precipitation, precipitation of the wettest month, and precipitation of the driest month. The annual precipitation may be the key factor restricting the northward distribution of S. cathayensis. In general, the centroid of the distribution of S. cathayensis will move northward. The centroid of the adaptive habitats will move northward with the highest degree of climate abnormality. We think that Hainan Island is the most likely origin of S. cathayensis. These findings provide a theoretical basis for the establishment of genetic resources protection measures, the construction of core germplasm resources, and the study of the formation and evolution of Hamamelidaceae.