Evolutionary history of a widespread tree species Acer mono in East Asia

The Phylogeography of Deciduous Tree Ulmus macrocarpa (Ulmaceae) in Northern China

Disentangling how climate oscillations and geographical events significantly influence plants' genetic architecture and demographic history is a central topic in phytogeography. The deciduous ancient tree species Ulmus macrocarpa is primarily distributed throughout Northern China and has timber and horticultural value. In the current study, we studied the phylogenic architecture and demographical history of U. macrocarpa using chloroplast DNA with ecological niche modeling. The results indicated that the populations' genetic differentiation coefficient (NST) value was significantly greater than the haplotype frequency (GST) (p < 0.05), suggesting that U. macrocarpa had a clear phylogeographical structure. Phylogenetic inference showed that the putative chloroplast haplotypes could be divided into three groups, in which the group Ⅰ was considered to be ancestral. Despite significant genetic differentiation among these groups, gene flow was detected. The common ancestor of all haplotypes was inferred to originate in the middle-late Miocene, followed by the haplotype overwhelming diversification that occurred in the Quaternary. Combined with demography pattern and ecological niche modeling, we speculated that the surrounding areas of Shanxi and Inner Mongolia were potential refugia for U. macrocarpa during the glacial period in Northern China. Our results illuminated the demography pattern of U. macrocarpa and provided clues and references for further population genetics investigations of precious tree species distributed in Northern China.

Wind-dispersed seeds blur phylogeographic breaks: The complex evolutionary history of Populus lasiocarpa around the Sichuan Basin

The strength of phylogeographic breaks can vary among species in the same area despite being subject to the same geological and climate history due to differences in biological traits. Several important phylogeographic breaks exist around the Sichuan Basin in Southwest China but few studies have focused on wind-dispersed plants. Here, we investigated the phylogeographic patterns and the evolutionary history of Populus lasiocarpa, a wind-pollinated and wind-dispersed tree species with a circum-Sichuan Basin distribution in southwest China. We sequenced and analyzed three plastid DNA fragments (ptDNA) and eight nuclear microsatellites (nSSRs) of 265 individuals of P. lasiocarpa from 21 populations spanning the entire distribution range. Distribution patterns based on nSSR data revealed that there are three genetic groups in P. lasiocarpa. This is consistent with the three phylogeographic breaks (Sichuan Basin, the Kaiyong Line and the 105°E line), where the Sichuan basin acts as the main barrier to gene flow between western and eastern groups. However, the distribution pattern based on ptDNA haplotypes poorly matched the phylogeographic breaks, and wind-dispersed seeds may be one of the main contributing factors. Species distribution modelling suggested a larger potential distribution in the last glacial maximum with a severe bottleneck during the last interglacial. A DIYABC model also suggested a population contraction and expansion for both western and eastern lineages. These results indicate that biological traits are likely to affect the evolutionary history of plants, and that nuclear molecular markers, which experience higher levels of gene flow, might be better indicators of phylogeographic breaks.

Genetic structure and trait variation within a maple hybrid zone underscore North China as an overlooked diversity hotspot

Tertiary relict flora in East Asia can be divided into northern and southern regions. North China is a diversity hotspot because it can be the secondary contact zone of ancient lineages from the two regions. To test the extent of ancient lineages hybridization and distinguish between the putative species pair Acer pictum subsp. mono and Acer truncatum, we conducted genetic and ecological studies within a maple hybrid zone in North China. Our results suggest that the two lineages of Acer coexist in the hybrid zone and that adult and offspring populations show typical bimodal genetic patterns. Hybrid individuals are established at intermediate altitudes between the two parental lineages. Flowering phenology is divergent between lineages, whereas the complex sexual system of Acer may ensure pollination among lineages. Leaf and fruit morphologies are different between the northern and southern origin lineages, corresponding to A. pictum subsp. mono and A. truncatum, respectively. Reduced gene flow between lineages suggests that they should be considered as two species. However, large morphological variations within each species and the existence of hybrids offer low reliability of species identification based solely on morphological traits. Our study underscores North China as an overlooked diversity hotspot that requires further study in the future.

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

Genome-scale phylogeography resolves the native population structure of the Asian longhorned beetle, Anoplophora glabripennis (Motschulsky)

Human-assisted movement has allowed the Asian longhorned beetle (ALB, Anoplophora glabripennis (Motschulsky)) to spread beyond its native range and become a globally regulated invasive pest. Within its native range of China and the Korean peninsula, human-mediated dispersal has also caused cryptic translocation of insects, resulting in population structure complexity. Previous studies used genetic methods to detangle this complexity but were unable to clearly delimit native populations which is needed to develop downstream biosurveillance tools. We used genome-wide markers to define historical population structure in native ALB populations and contemporary movement between regions. We used genotyping-by-sequencing to generate 6102 single-nucleotide polymorphisms (SNPs) and amplicon sequencing to genotype 53 microsatellites. In total, we genotyped 712 individuals from ALB's native distribution. We observed six distinct population clusters among native ALB populations, with a clear delineation between northern and southern groups. Most of the individuals from South Korea were distinct from populations in China. Our results also indicate historical divergence among populations and suggest limited large-scale admixture, but we did identify a restricted number of cases of contemporary movement between regions. We identified SNPs under selection and describe a clinal allele frequency pattern in a missense variant associated with glycerol kinase, an important enzyme in the utilization of an insect cryoprotectant. We further demonstrate that small numbers of SNPs can assign individuals to geographic regions with high probability, paving the way for novel ALB biosurveillance tools.

Genomic Data Reveals Profound Genetic Structure and Multiple Glacial Refugia in Lonicera oblata (Caprifoliaceae), a Threatened Montane Shrub Endemic to North China

Characterizing genetic diversity and structure and identifying conservation units are both crucial for the conservation and management of threatened species. The development of high-throughput sequencing technology provides exciting opportunities for conservation genetics. Here, we employed the powerful SuperGBS method to identify 33, 758 high-quality single-nucleotide polymorphisms (SNP) from 134 individuals of a critically endangered montane shrub endemic to North China, Lonicera oblata. A low level of genetic diversity and a high degree of genetic differentiation among populations were observed based on the SNP data. Both principal component and phylogenetic analyses detected seven clusters, which correspond exactly to the seven geographic populations. Under the optimal K = 7, Admixture suggested the combination of the two small and geographically neighboring populations in the Taihang Mountains, Dongling Mountains, and Lijiazhuang, while the division of the big population of Jiankou Great Wall in the Yan Mountains into two clusters. High population genetic diversity and a large number of private alleles were detected in the four large populations, while low diversity and non-private alleles were observed for the remaining three small populations, implying the importance of these large populations as conservation units in priority. Demographic history inference suggested two drastic contractions of population size events that occurred after the Middle Pleistocene Transition and the Last Glacial Maximum, respectively. Combining our previous ecological niche modeling results with the present genomic data, there was a possible presence of glacial refugia in the Taihang and Yan Mountains, North China. This study provides valuable data for the conservation and management of L. oblata and broadens the understanding of the high biodiversity in the Taihang and Yan Mountains.

Different rates of pollen and seed gene flow cause branch-length and geographic cytonuclear discordance within Asian butternuts

Topological cytonuclear discordance is commonly observed in plant phylogenetic and phylogeographic studies, yet few studies have attempted to detect two other forms of cytonuclear discordance (branch length and geographical) and to uncover the causes of the discordance. We used the whole nuclear and chloroplast genome data from 80 individual Asian butternuts to reveal the pattern and processes of cytonuclear discordance. Our findings indicate that the chloroplast genome had substantially deeper divergence (branch-length discordance) and a steeper cline in the contact zone (geographic discordance) compared with the nuclear genome. After various hypothesis have been tested, the results suggest that incomplete lineage sorting, positive selection and cytonuclear incompatibility are probably insufficient to explain this pattern. However, isolation-by-distance analysis and gene flow estimation point to a much higher level of gene flow by pollen compared with by seeds, which may have slowed down lineage divergence and mediated wider contact for nuclear genome compared with the chloroplast genome. Altogether, this study highlights a critical role of sex-biased dispersal in causing discordance between the nuclear and plastid genome of Asian butternuts. Given its ubiquity among plants, asymmetric gene flow should be given a high priority in future studies of cytonuclear discordance.

Distinct late Pleistocene subtropical-tropical divergence revealed by fifteen low-copy nuclear genes in a dominant species in South-East China

In East Asia, genetic divergence is usually considered to be correlated to different floristic regions, however, subtropical-tropical divergence is largely ignored, compared to widely explored temperate-subtropical divergence. Lindera aggregata (Lauraceae), a dominant species in South-East China was selected to address this issue. Fifteen low-copy nuclear genes (LCGs) and four chloroplast DNA (cpDNA) fragments were used to detect its evolutionary history. In LCGs, STRUCTURE and dated Bayesian phylogeny analyses detect distinct subtropical-tropical divergence since late Pleistocene. Approximate Bayesian calculation (ABC) further supports the distinct subtropical-tropical divergence, and close related Taiwan and South China populations are diverged at the last interglacial. Isolation by distance, isolation by environment and isolation by resistance analyses suggest the current climatic difference rather than geographical distance contributes to the genetic differentiation. Principle component analysis shows populations of tropical cluster occur in warmer area with higher precipitation. Ancestral area reconstruction based on Bayesian phylogeny indicates that ancestral L. aggregata populations are distributed in tropical region. In cpDNA, although unique haplotypes are found in tropical region, distinct subtropical-tropical divergence is absent. In conclusion, distinct late Pleistocene subtropical-tropical divergence of L. aggregata is triggered by climate. It is likely that L. aggregata is originated in Southwest-South China and experienced hierarchical dispersal from south to north. The South China Sea land bridge has dual role in connecting or isolating Taiwan and mainland populations since the last glaciation.

Insights Into the Significance of the Chinense Loess Plateau for Preserving Biodiversity From the Phylogeography of Speranskia tuberculata (Euphorbiaceae)

The significance of the Chinese Loess Plateau (CLP) in maintaining biodiversity for northern China has rarely been shown, as previous phylogeographic studies are mostly woody species and they have revealed that Quaternary refugia are mainly located in mountain regions. We selected a drought-enduring endemic herb, Speranskia tuberculata (Euphorbiaceae), to determine its glacial refugia and postglacial demographic history. To this end, we sampled 423 individuals from 38 populations covering its entire geographic distribution. Three chloroplast DNA (cpDNA) fragments, two low-copy nuclear genes, and six nuclear microsatellites (nSSRs) were used and supplemented with ecological niche modeling (ENM) to infer the phylogeographic history of this species. Populations with private haplotypes and high haplotype diversity of cpDNA are mainly located in the CLP or scattered around northeastern China and the coastal region. Spatial expansion, detected using a neutrality test and mismatch distribution, may have resulted in a widely distributed ancestral cpDNA haplotype, especially outside of the CLP. For nuclear DNA, private haplotypes are also distributed mainly in the CLP. In nSSRs, STRUCTURE clustering identified two genetic clusters, which are distributed in the west (western cluster) and east (eastern cluster), respectively. Many populations belonged, with little to no admixture, to the western cluster while (hardly) pure populations of the eastern cluster were barely found. Genetic differentiation is significantly correlated with geographic distance, although genetic diversity is uniformly distributed. ENM suggests that the distribution of S. tuberculata has recently expanded northwards from the southern CLP, whereas it has experienced habitat loss in the south. Thus, S. tuberculata populations probably survived the last glacial maximum (LGM) in the southern CLP and experienced post-glacial expansion. Wind-dispersed pollen could bring the majority of genotypes to the front during spatial expansion, resulting in uniformly distributed genetic diversity. Based on evidence from molecular data and vegetation and climate changes since the LGM, we conclude that drought-enduring species, especially herbaceous species, are likely to have persisted in the CLP during the LGM and to have experienced expansion to other regions in northern China.

Evolutionary history of endangered and relict tree species Dipteronia sinensis in response to geological and climatic events in the Qinling Mountains and adjacent areas

Geological and climatic events are considered to profoundly affect the evolution and lineage divergence of plant species. However, the evolutionary histories of tree species that have responded to past geological and climate oscillations in central China's mountainous areas remain mostly unknown. In this study, we assessed the evolutionary history of the endangered and relict tree species Dipteronia sinensis in the Qinling Mountains (QM) and adjacent areas in East Asia based on variations in the complete chloroplast genomes (cpDNA) and reduced-genomic scale single nucleotide polymorphisms (SNPs). Population structure and phylogenetic analysis based on the cpDNA variations suggested that D. sinensis could be divided into two intraspecific genetic lineages in the eastern and western sides of the QM (EQM and WQM, respectively) in East Asia. Molecular dating suggested that the intraspecific divergence of D. sinensis occurred approximately 39.2 million years ago during the later Paleogene. It was significantly correlated with the orogeny of the QM, where the formation of this significant geographic barrier in the region may have led to the divergence of independent lineages. Bayesian clustering and demographic analysis showed that intraspecific gene flow was restricted between the EQM and WQM lineages. Isolation-with-migration analysis indicated that the two genetic lineages experienced significant demographic expansions after the Pleistocene ice ages. However, the genetic admixture was determined in some populations between the two lineages by the large scale of SNP variations due to DNA incompatibility, the large significant population size, and rapid gene flow of nuclear DNA markers. Our results suggest that two different conservation and management units should be constructed for D. sinensis in the EQM and WQM areas. These findings provide novel insights into the unprecedented effects of tectonic changes and climatic oscillations on lineage divergence and plant population evolution in the QM and adjacent areas in East Asia.

Biogeography of land snail genus Acusta (Gastropoda: Camaenidae): Diversification on East Asian islands

East Asia has highly diverse and endemic biota due to its complex geological and climatic history and its diversified topography. The continental and insular distributions of land snail genus Acusta in East Asia provide a good opportunity to compare the evolutionary processes in this group under different biogeographical conditions. In this study, we inferred the evolutionary history of the land snail genus Acusta by a molecular phylogeny and investigated how the palaeogeographic events shaped species diversity and the distribution of the Acusta genus within the island arc. A concatenated dataset generated from sequences of one nuclear (ITS2) and two mitochondrial (16S, COI) gene fragments, include most of nominal taxa of the genus, four related species and one outgroup. We constructed the phylogeny and the evolutionary history of the genus through maximum-likelihood and Bayesian inference methods, using a Bayesian molecular clock and ancestral range estimation. Our results suggested that currently recognized species in Acusta are polyphyletic. The traditionally accepted concept of the affinity of Acusta and Bradybaena is not supported. The hypothesis of colonization via land bridges during the Pleistocene glaciations for the biota of East Asian islands is not supported. Instead, the origin and diversification of the genus Acusta was dated to the late Miocene-Pliocene from an area around North and Northeast China to South China and East Asian islands Three major evolutionary lineages were identified. Two of the major lineages demonstrate distinct evolutionary histories, as sympatric speciation is the major speciation process for the continental clade, while the insular clade originated from founder events. Taiwan functioned as an important source of diversification for species on the East Asian islands possibly through passive dispersal of different mechanisms. The sea level fluctuations caused by the Pleistocene glacial cycles play a role in the subsequent dispersion and diversification of species of the continental clade, such as the more recent range expansion of A. redfieldi from South China to Taiwan and Japan.

Rear-edge, low-diversity, and haplotypic uniformity in cold-adapted Bupleurum euphorbioides interglacial refugia populations

The high genetic diversity of rear-edge refugia populations is predicted to have resulted from species repeatedly migrating to low latitudes during glacial periods over the course of Quaternary climate change. However, several recent empirical studies of cold-tolerant plants revealed the opposite pattern. We investigated whether current habitats of the cold-adapted and range-restricted Bupleurum euphorbioides in the Baekdudaegan, South Korea, and North Korea could be interglacial refugia, and documented how their rear-edge populations differ genetically from those of typical temperate species. Phylogeographic analysis and ecological niche modeling (ENM) were used. The genetic structure was analyzed using microsatellite markers and chloroplast DNA sequences. The congener B. longiradiatum was included as a typical temperate plant species. Despite having almost identical life history traits, these congeneric species exhibited contrasting patterns of genetic diversity. ENM revealed an apparent maximum range contraction during the last interglacial. In contrast, its range expanded northward to the Russian Far East (Primorsky) during the Last Glacial Maximum. Thus, we hypothesize that B. euphorbioides retreated to its current refugia during interglacial periods. Unlike populations in the central region, the rear-edge populations were genetically impoverished and uniform, both within populations and in pooled regional populations. The rear-edge B. euphorbioides survived at least one past interglacial, contributing to the species' genetic diversity. We believe that such genetic variation in the cold-adapted B. euphorbioides gives the species the necessary adaptations to survive an upcoming favorable environment (the next glacial), unless there is artificial environmental change.

A Genomewide Scan for Genetic Structure and Demographic History of Two Closely Related Species, Rhododendron dauricum and R. mucronulatum (Rhododendron, Ericaceae)

Understanding the processes of divergence and speciation is an important task for evolutionary research, and climate oscillations play a pivotal role. We estimated the genetic structure and demographic history of two closely related species of Rhododendron, R. dauricum, and R. mucronulatum, distributed in northeastern China using 664,406 single nucleotide polymorphic loci of specific-locus amplified fragment sequencing (SLAF-seq) and 4 chloroplast DNA (cpDNA) fragments, sampling 376 individuals from 39 populations of these two species across their geographic distributions. The geographical distribution of cpDNA haplotypes revealed that R. dauricum and R. mucronulatum have different spatial genetic structures and haplotype diversity. Analysis of molecular variance (AMOVA) results showed that these two species have significant genetic differentiation and that the phylogeny demonstrates that these two species clustered a monophyletic group based on SLAF data, respectively, but not in cpDNA data. The evidence of significant gene flow was also detected from R. mucronulatum to R. dauricum. A deep divergence between the two species was observed and occurred during the early Oligocene. The niche models showed that the two species have different demographic histories. Thus, our results imply that geography and climate changes played important roles in the evolutionary process of R. dauricum and R. mucronulatum, and although there was an interspecific gene flow, the divergence was maintained by natural selection.

Plastome comparative genomics in maples resolves the infrageneric backbone relationships

Maples (Acer) are among the most diverse and ecologically important tree genera of the north-temperate forests. They include species highly valued as ornamentals and as a source of timber and sugar products. Previous phylogenetic studies employing plastid markers have not provided sufficient resolution, particularly at deeper nodes, leaving the backbone of the maple plastid tree essentially unresolved. We provide the plastid genome sequences of 16 species of maples spanning the sectional diversity of the genus and explore the utility of these sequences as a source of information for genetic and phylogenetic studies in this group. We analyzed the distribution of different types of repeated sequences and the pattern of codon usage, and identified variable regions across the plastome. Maximum likelihood and Bayesian analyses using two partitioning strategies were performed with these and previously published sequences. The plastomes ranged in size from 155,212 to 157,023 bp and had structure and gene content except for Acer palmatum (sect. Palmata), which had longer inverted repeats and an additional copy of the rps19 gene. Two genes, rps2 and rpl22, were found to be truncated at different positions and might be non-functional in several species. Most dispersed repeats, SSRs, and overall variation were detected in the non-coding sequences of the LSC and SSC regions. Fifteen loci, most of which have not been used before in the genus, were identified as the most variable and potentially useful as molecular markers for barcoding and genetic studies. Both ML and Bayesian analyses produced similar results irrespective of the partitioning strategy used. The plastome-based tree largely supported the topology inferred in previous studies using cp markers while providing resolution to the backbone relationships but was highly incongruous with a recently published nuclear tree presenting an opportunity for further research to investigate the causes of discordance, and particularly the role of hybridization in the diversification of the genus. Plastome sequences are valuable tools to resolve deep-level relationships within Acer. The variable loci and SSRs identified in this study will facilitate the development of markers for ecological and evolutionary studies in the genus. This study underscores the potential of plastid genome sequences to improve our understanding of the evolution of maples.

Different processes shape the patterns of divergence in the nuclear and chloroplast genomes of a relict tree species in East Asia

Isolation by spatial distance (IBD), environment (IBE), and historical climatic instability (IBI) are three common processes assessed in phylogeographic and/or landscape genetic studies. However, the relative contributions of these three processes with respect to spatial genetic patterns have seldom been compared. Moreover, whether the relative contribution differs in different regions or when assessed using different genetic markers has rarely been reported. Lindera obtusiloba has been found to have two sister genetic clades of chloroplast (cpDNA) and nuclear microsatellite (nSSR), both of which show discontinuous distribution in northern and southern East Asia. In this study, we used the Mantel test and multiple matrix regression with randomization (MMRR) to determine the relative contributions of IBD, IBE, and IBI with respect to L. obtusiloba populations. Independent Mantel tests and MMRR calculations were conducted for two genetic data sets (cpDNA and nSSR) and for different regions (the overall species range, and northern and southern subregions of the range). We found a significant IBI pattern in nSSR divergence for all assessed regions, whereas no clear IBI pattern was detected with respect to cpDNA. In contrast, significant (or marginal) divergent IBD patterns were detected for cpDNA in all regions, whereas although a significant IBE was apparent with respect to the overall range, the effect was not detected in the two subregions. The differences identified in nSSR and cpDNA population divergence may be related to differences in the heredity and ploidy of the markers. Compared with the southern region, the northern region showed less significant correlation patterns, which may be related to the shorter population history and restricted population range. The findings of this study serve to illustrate that comparing between markers or regions can contribute to gaining a better understanding the population histories of different genomes or within different regions of a species' range.

Comparative Plastome Analyses and Phylogenetic Applications of the Acer Section Platanoidea

The Acer L. (Sapindaceae) is one of the most diverse and widespread genera in the Northern Hemisphere. Section Platanoidea harbours high genetic and morphological diversity and shows the phylogenetic conflict between A. catalpifolium and A. amplum. Chloroplast (cp) genome sequencing is efficient for the enhancement of the understanding of phylogenetic relationships and taxonomic revision. Here, we report complete cp genomes of five species of Acer sect. Platanoidea. The length of Acer sect. Platanoidea cp genomes ranged from 156,262 bp to 157,349 bp and detected the structural variation in the inverted repeats (IRs) boundaries. By conducting a sliding window analysis, we found that five relatively high variable regions (trnH-psbA, psbN-trnD, psaA-ycf3, petA-psbJ and ndhA intron) had a high potential for developing effective genetic markers. Moreover, with an addition of eight plastomes collected from GenBank, we displayed a robust phylogenetic tree of the Acer sect. Platanoidea, with high resolutions for nearly all identified nodes, suggests a promising opportunity to resolve infrasectional relationships of the most species-rich section Platanoidea of Acer.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Population genetic structures of two ecologically distinct species Betula platyphylla and B. ermanii inferred based on nuclear and chloroplast DNA markers

Climatic oscillations during the last glacial maximum (LGM) significantly affected the distribution patterns and genetic structure of extant plants. Northeast China (NEC) is a major biodiversity center in East Asia, and the influence of historical climate change on NEC populations is critical for understanding species responses to future climate change. However, only a few phylogeographic studies of cool temperate deciduous tree species have been conducted in the area, and results are inconsistent for species with different niches or distribution areas. We employed multiple chloroplast and nuclear markers to investigate the genetic structure of two ecologically contrasting species, Betula platyphylla and B. ermanii, in NEC. Rare haplotypes were identified in the chloroplast genome of these species, and both exhibited high levels of nucleotide diversity based on a fragment of the nuclear gene G3PDH and microsatellites. Moreover, significant phylogeographic structure was detected for B. platyphylla, suggesting that these populations had recolonized from independent glacial refuges, whereas no genetic structure was found for B. ermanii.

OPEN RESEARCH BADGES

The nSSR datasets used in the current study and the table of pairwise FST (below diagonal) and its standardized F'ST (above diagonal) among 25 populations based on seven SSRs are available from the Dryad (DOI: https://doi.org/10.5061/dryad.230d176). Sequences generated from this study were deposited in GenBank under Accession nos. KY199568-KY200162 and MK819541-MK819970.

Phylogeography and paleodistribution models of a widespread birch (Betula platyphylla Suk.) across East Asia: Multiple refugia, multidirectional expansion, and heterogeneous genetic pattern

Widespread tree species cover large geographical areas and play important roles in various vegetation types. Understanding how these species responded to historical climatic changes is important for understanding community assembly mechanisms with evolutionary and conservation implications. However, the location of refugial areas and postglacial history of widespread trees in East Asia remain poorly known. We combined microsatellite data (63 populations, 1756 individuals) and ecological niche modeling to examine the range-wide population diversity, genetic structure, and historical demography of a pioneer tree species, Asian white birch (Betula platyphylla Suk.) across East Asia. We found a north-to-south trend of declining genetic diversity and five clusters, corresponding to geographical regions. Different clusters were inferred to have diverged through Pleistocene climatic oscillations and have different expansion routes, leading to genetic admixture in some populations. Ecological niche models indicated that the distribution of B. platyphylla during the last glacial maximum still had a large latitude span with slight shifts toward southeast, and northern populations had more variable distribution ranges than those in the south during later climatic oscillations. Our results reflect the relatively stable distribution through the last glacial-interglacial cycles and recent multidirectional expansion of B. platyphylla, providing new hypotheses for the response pattern of widespread tree species to climate change. The gradual genetic pattern from northeast to southwest and alternative distribution dynamics possibly resulted from environmental differences caused by latitude and topographic heterogeneity.

Phylogeography of Schisandra chinensis (Magnoliaceae) Reveal Multiple Refugia With Ample Gene Flow in Northeast China

Temperate conifers and broadleaved mixed forests in northeast China are ideal to investigate the genetic consequences of climate changes during the last glacial maximum (LGM), 29 - 16 kya. As previous studies were focused on tree species with long generation time; here, the evolutionary history of Schisandra chinensis, a climber species with a generation time of five years, was investigated using chloroplast DNA (cpDNA), nuclear single copy gene (nSCG), and nuclear single sequence repeats (nSSRs, i.e., microsatellite) markers, along with ecological niche modeling (ENM), which predicted a suitable habitat in Korea Peninsula (KP) during the LGM. Private haplotypes and high genetic diversity of both cpDNA and nSCG were mainly found in KP and Changbai Mt. (CB). Although no significant phylogeographic structure was detected in the cpDNA and nSCG, three nSSRs clusters roughly distributed in west (CB and KP), east (north China), and north (Xiaoxing'an Range, XR) regions were found in Structure analysis. The approximate Bayesian computation analysis showed the west cluster diverged at 35.45 kya, and the other two clusters at 19.85 kya. The genetic diversity calculated for each of the three markers showed no significant correlation with latitude. Genetic differentiation of nSSRs was also not correlated with geographic distance. Migrate analysis estimated extensive gene flow between almost all genetic cluster pairs and BOTTLENECK analysis showed that few populations experienced severe bottlenecks. Overall, results indicate that S. chinensis survived the LGM in situ in multiple refugia, which likely include two macrorefugia (KP and CB) and two microrefugia (XR and north China). Extensive postglacial gene flow among the three nSSRs clusters led to uniformly distributed genetic diversity and low genetic differentiation.

Combining mitochondrial and nuclear genome analyses to dissect the effects of colonization, environment, and geography on population structure in Pinus tabuliformis

The phylogeographic histories of plants in East Asia are complex and shaped by both past large-scale climatic oscillations and dramatic tectonic events. The impact of these historic events, as well as ecological adaptation, on the distribution of biodiversity remains to be elucidated. Pinus tabuliformis is the dominant coniferous tree in northern China, with a large distribution across wide environmental gradients. We examined genetic variation in this species using genotyping-by-sequencing and mitochondrial (mt) DNA markers. We found population structure on both nuclear and mt genomes with a geographic pattern that corresponds well with the landscape of northern China. To understand the contributions of environment, geography, and colonization history to the observed population structure, we performed ecological niche modeling and partitioned the among-population genomic variance into isolation by environment (IBE), isolation by distance (IBD), and isolation by colonization (IBC). We used mtDNA, which is transmitted by seeds in pine, to reflect colonization. We found little impact of IBE, IBD, and IBC on variation in neutral SNPs, but significant impact of IBE on a group of outlier loci. The lack of IBC illustrates that the maternal history can be quickly eroded from the nuclear genome by high rates of gene flow. Our results suggest that genomic variation in P. tabuliformis is largely affected by neutral and stochastic processes, and the signature of local adaptation is visible only at robust outlier loci. This study enriches our understanding on the complex evolutionary forces that shape the distribution of genetic variation in plant taxa in northern China, and guides breeding, conservation, and reforestation programs for P. tabuliformis.

Recent Fragmentation May Not Alter Genetic Patterns in Endangered Long-Lived Species: Evidence From Taxus cuspidata

Forestland fragmentation caused by overexploitation of forest resources can in principle reduce genetic diversity, limit gene flow and eventually lead to species developing strong genetic structure. However, the genetic consequences of recent anthropogenic fragmentation of tree species remain unclear. Taxus cuspidata, which has extremely small populations distributed mainly in Changbai Mt. in Northeast (NE) China, has recently endured severe habitat fragmentation. Here, we investigate the pattern of genetic diversity and structure, identify risk factors, predict the future distribution and finally provide guidelines for the conservation and management of this species. We used three chloroplast and two mitochondrial DNA fragments, which are both paternally inherited in yews but differ in mutation rates, to genotype a total of 265 individuals from 26 populations covering the distribution of the species in China. Both chloroplast and mitochondrial data showed high degrees of genetic diversity, extensive gene flow over the entire geographical range and historical stability of both effective population size and distribution of the species. However, ecological niche modeling suggests a decrease in suitable areas for this species by the years 2050 and 2070. The maintenance of high genetic diversity and the existence of sufficient gene flow suggest that recent fragmentation has not affected the genetic composition of the long-lived tree T. cuspidata. However, severe impacts of anthropogenic activities are already threatening the species. Conservation and management strategies should be implemented in order to protect the remnant populations.

The Application and Limitation of Universal Chloroplast Markers in Discriminating East Asian Evergreen Oaks

The East Asian subtropics mostly occupied by evergreen broad-leaved forests (EBLFs), is one of the global diversity centers for evergreen oaks. Evergreen oaks are keystone canopy trees in EBLFs with important ecosystem function and crucial significance for regional biodiversity conservation. However, the species composition and diversity of Asian evergreen oaks are poorly understood. Here, we test whether the four chloroplast markers atpI-atpH, matK, psbA-trnH, and ycf1, can discriminate the two evergreen oak sections in Asia - Cyclobalanopsis and Ilex. Two hundred and seventy-two individuals representing 57 species were scanned and 17 species from other oaks sections were included for phylogenetic reconstruction. The genetic diversity of the Quercus sections was also compared. Overall, we found that universal chloroplast DNA (cpDNA) barcoding markers could resolve two clades in Quercus, i.e., subgenus Cerris (Old World Clade) and subgenus Quercus (New World Clade). The chloroplast markers distinguished the main sections, with few exceptions. Each cpDNA region showed no barcoding gap and none of them provided good resolution at the species level. The best species resolution (27.78%) was obtained when three or four markers were combined and analyzed using BLAST. The high conservation of the cpDNA and complicated evolutionary patterns, due to incomplete lineage sorting, interspecific hybridization and introgressions may hinder the ability of cpDNA markers to discriminate different species. When comparing diversification pattern across Quercus sections (Cyclobalanopsis, Ilex, Cerris, Quercus, and Protobalanus), we found that section Ilex was the most genetically diverse, and section Cyclobalanopsis was lower genetically diverse. This diversification pattern may have resulted from the interplay of the Eurasia Cenozoic tectonic movements, climate changes and different niches of their ancestral lineages.

Comparison of genetic variation between northern and southern populations of Lilium cernuum (Liliaceae): Implications for Pleistocene refugia

The so-called "Baekdudaegan" (BDDG), a mountain range that stretches along the Korean Peninsula, has been recently proposed as a major "southern" glacial refugium for boreal or temperate plant species based on palaeoecological and, especially, genetic data. Genetic studies comparing genetic variation between population occurring on the BDDG and more northern ones (i.e. in NE China and/or in Russian Far East) are, however, still too few to draw firm conclusions on the role of the BDDG as a refugium and a source for possible northward post-glacial recolonizations. In order to fill this gap, we selected a boreal/temperate herb, Lilium cernuum, and compared levels of allozyme-based genetic diversity of five populations from NE China with five populations from South Korea (home of its hypothesized refuge areas). As a complementary tool, we used the maximum entropy algorithm implemented in MaxEnt to infer the species' potential distribution for the present time, which was projected to different past climate scenarios for the Last Glacial Maximum (LGM). Permutation tests revealed that Korean populations harbored significantly higher levels of within-population genetic variation than those from NE China (expected heterozygosity = 0.173 vs. 0.095, respectively). Our results suggest that the lowered levels of genetic diversity in NE Chinese populations might be due to founder effects associated with post-glacial migration from southern regions. Congruent with genetic data, past distribution models showed higher probability of occurrence in southern ranges than in northern ones during the LGM. In addition, a positive correlation was detected between the expected heterozygosity and environmental LGM suitability. From a conservation perspective, our results further suggest that the southern populations in South Korea may be particularly worthy of protection.

Species Delimitation and Lineage Separation History of a Species Complex of Aspens in China

Species delimitation in tree species is notoriously challenging due to shared polymorphisms among species. An integrative survey that considers multiple operational criteria is a possible solution, and we aimed to test it in a species complex of aspens in China. Genetic [four chloroplast DNA (cpDNA) fragments and 14 nuclear microsatellite loci (nSSR)] and morphological variations were collected for 76 populations and 53 populations, respectively, covering the major geographic distribution of the Populus davidiana-rotundifolia complex. Bayesian clustering, analysis of molecular variance (AMOVA), Principle Coordinate Analysis (PCoA), ecological niche modeling (ENM), and gene flow (migrants per generation), were employed to detect and test genetic clustering, morphological and habitat differentiation, and gene flow between/among putative species. The nSSR data and ENM suggested that there are two separately evolving meta-population lineages that correspond to P. davidiana (pd) and P. rotundifolia (pr). Furthermore, several lines of evidence supported a subdivision of P. davidiana into Northeastern (NEC) and Central-North (CNC) groups, yet they are still functioning as one species. CpDNA data revealed that five haplotype clades formed a pattern of [pdNEC, ((pdCNC, pr), (pdCNC, pr))], but most haplotypes are species-specific. Meanwhile, PCA based on morphology suggested a closer relationship between the CNC group (P. davidiana) and P. rontundifolia. Discrepancy of nSSR and ENM vs. cpDNA and morphology could have reflected a complex lineage divergence and convergence history. P. davidiana and P. rotundifolia can be regarded as a recently diverged species pair that experienced parapatric speciation due to ecological differentiation in the face of gene flow. Our findings highlight the importance of integrative surveys at population level, as we have undertaken, is an important approach to detect the boundary of a group of species that have experienced complex evolutionary history.

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation-interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.

Phylogeographic breaks within Asian butternuts indicate the existence of a phytogeographic divide in East Asia

East Asia has been hypothesized to be subdivided into two distinct northern and southern areas, separated by a band of dry climate that was far more severe in the early Tertiary but still exists today. However, this biogeographic hypothesis has rarely been tested using a molecular phylogeographic approach. We genotyped 70 populations throughout the distributional range of Asian butternuts (Juglans section Cardiocaryon) using eight chloroplast DNA regions, one single-copy nuclear gene, and 17 nuclear microsatellite loci, supplemented with paleodistribution modeling of the major genetic clades. The genetic data consistently identified two clades, one northern, comprising Juglans mandshurica and Juglans ailantifolia, and one southern, comprising Juglans cathayensis. The two clades diverged through climate-induced vicariance of an ancestral northern range during the mid-Miocene and remained mostly separate thereafter, with geographical isolation of the Japanese Islands and refugial isolation or secondary contacts in the late Pleistocene producing further subdivision within the northern clade. But beyond all that, we also discovered a role of environmental adaptation in maintaining and/or reinforcing the north-south divergence. Asian butternuts offer a strong case for the existence of a biogeographic divide between the northern and southern parts of East Asia during the Neogene and into the Pleistocene.

Efficient Identification of the Forest Tree Species in Aceraceae Using DNA Barcodes

Aceraceae is a large forest tree family that comprises many economically and ecologically important species. However, because interspecific and/or intraspecific morphological variations result from frequent interspecific hybridization and introgression, it is challenging for non-taxonomists to accurately recognize and identify the tree species in Aceraceae based on a traditional approach. DNA barcoding is a powerful tool that has been proposed to accurately distinguish between species. In this study, we assessed the effectiveness of three core standard markers (matK, rbcL and ITS) plus the chloroplast locus trnS-trnG as Aceraceae barcodes. A total of 231 sequences representing 85 species in this forest family were collected. Of these four barcode markers, the discrimination power was highest for the ITS (I) region (50%) and was progressively reduced in the other three chloroplast barcodes matK (M), trnS-trnG (T) and rbcL (R); the discrimination efficiency of the ITS marker was also greater than any two-locus combination of chloroplast barcodes. However, the combinations of ITS plus single or combined chloroplast barcodes could improve species resolution significantly; T+I (90.5% resolution) and R+M+T+I (90.5% resolution) differentiated the highest portion of species in Aceraceae. Our current results show that the nuclear ITS fragment represents a more promising DNA barcode marker than the maternally inherited chloroplast barcodes. The most efficient and economical method to identify tree species in Aceraceae among single or combined DNA barcodes is the combination of T+I (90.5% resolution).

Contributions of multiple refugia during the last glacial period to current mainland populations of Korean pine (Pinus koraiensis)

The northern microrefugia that existed during the Last Glacial Maximum (LGM) are a key factor in the demographic history of species. Pinus koraiensis has a unique distribution in northeast Asia. The Changbai Mountains and the Korean peninsula (CM/KP) are usually considered to be the LGM refugia for P. koraiensis. However, the Xiaoxingan Range (XR), at the northern part of this species' distribution, is another possible refugium. We used chloroplast sequencing and ten nuclear single-copy gene loci to calculate the genetic diversity pattern of P. koraiensis. The probabilities of a single LGM refugium and of multiple LGM refugia were calculated based on approximate Bayesian computation. The effect of the latitudinal gradient on genetic diversity was not significant. However, unique alleles occurred at low frequencies in CM/KP and XR. A conservative estimate of the coalescence time between CM/KP and XR is 0.4 million years ago, a time prior to the LGM. Gene flow between CM/KP and XR was estimated to be more than one in per generation, an amount that may be sufficient to limit genetic divergence between the regions. Our study strongly supports the hypothesis that XR was another LGM refugium in addition to CM/KP.