Molecular data and ecological niche modelling reveal a highly dynamic evolutionary history of the East Asian Tertiary relict Cercidiphyllum (Cercidiphyllaceae)

Which contributes more to the relict flora distribution pattern in East Asia, geographical processes or climate change? New evidence from the phylogeography of Rehderodendron kwangtungense

BACKGROUND

Relict species are important for enhancing the understanding of modern biogeographic distribution patterns. Although both geological and climatic changes since the Cenozoic have affected the relict flora in East Asia, the contributions of geographical processes remain unclear. In this study, we employed restriction-site associated DNA sequencing (RAD-seq) and shallow genome sequencing data, in conjunction with ecological niche modeling (ENM), to investigate the spatial genetic patterns and population differentiation history of the relict species Rehderodendron kwangtungense Chun.

RESULTS

A total of 138 individuals from 16 populations were collected, largely covering the natural distribution of R. kwangtungense. The genetic diversity within the R. kwangtungense populations was extremely low (HO = 0.048 ± 0.019; HE = 0.033 ± 0.011). Mantel tests revealed isolation-by-distance pattern (R2 = 0.38, P < 0.001), and AMOVA analysis showed that the genetic variation of R. kwangtungense occurs mainly between populations (86.88%, K = 7). Between 23 and 21 Ma, R. kwangtungense underwent a period of rapid differentiation that coincided with the rise of the Himalayas and the establishment of the East Asian monsoon. According to ENM and population demographic history, the suitable area and effective population size of R. kwangtungense decreased sharply during the glacial period and expanded after the last glacial maximum (LGM).

CONCLUSION

Our study shows that the distribution pattern of southern China mountain relict flora may have developed during the panplain stage between the middle Oligocene and the early Miocene. Then, the flora later fragmented under the force of orogenesis, including intermittent uplift during the Cenozoic Himalayan orogeny and the formation of abundant rainfall associated with the East Asian monsoon. The findings emphasized the predominant role of geographical processes in shaping relict plant distribution patterns.

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.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

Cryptic divergences and repeated hybridizations within the endangered "living fossil" dove tree (Davidia involucrata) revealed by whole genome resequencing

The identification and understanding of cryptic intraspecific evolutionary units (lineages) are crucial for planning effective conservation strategies aimed at preserving genetic diversity in endangered species. However, the factors driving the evolution and maintenance of these intraspecific lineages in most endangered species remain poorly understood. In this study, we conducted resequencing of 77 individuals from 22 natural populations of Davidia involucrata, a "living fossil" dove tree endemic to central and southwest China. Our analysis revealed the presence of three distinct local lineages within this endangered species, which emerged approximately 3.09 and 0.32 million years ago. These divergence events align well with the geographic and climatic oscillations that occurred across the distributional range. Additionally, we observed frequent hybridization events between the three lineages, resulting in the formation of hybrid populations in their adjacent as well as disjunct regions. These hybridizations likely arose from climate-driven population expansion and/or long-distance gene flow. Furthermore, we identified numerous environment-correlated gene variants across the total and many other genes that exhibited signals of positive evolution during the maintenance of two major local lineages. Our findings shed light on the highly dynamic evolution underlying the remarkably similar phenotype of this endangered species. Importantly, these results not only provide guidance for the development of conservation plans but also enhance our understanding of evolutionary past for this and other endangered species with similar histories.

Meta-analysis provides insights into the origin and evolution of East Asian evergreen broad-leaved forests

Evergreen broad-leaved forests (EBLFs) are dominated by a monsoon climate and form a distinct biome in East Asia with notably high biodiversity. However, the origin and evolution of East Asian EBLFs (EAEBLFs) remain elusive despite the estimation of divergence times for various representative lineages. Using 72 selected generic-level characteristic lineages, we constructed an integrated lineage accumulation rate (LAR) curve based on their crown ages. According to the crown-based LAR, the EAEBLF origin was identified at least as the early Oligocene (c. 31.8 million years ago (Ma)). The accumulation rate of the characteristic genera peaked at 25.2 and 6.4 Ma, coinciding with the two intensification periods of the Asian monsoon at the Oligocene - Miocene and the Miocene - Pliocene boundaries, respectively. Moreover, the LAR was highly correlated with precipitation in the EAEBLF region and negatively to global temperature, as revealed through time-lag cross-correlation analyses. An early Oligocene origin is suggested for EAEBLFs, bridging the gap between paleobotanical and molecular dating studies and solving conflicts among previous estimates based on individual representative lineages. The strong correlation between the crown-based LAR and the precipitation brought about by the Asian monsoon emphasizes its irreplaceable role in the origin and development of EAEBLFs.

Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub

Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome-wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid-Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub-central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post-glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

Evaluation of Climate Change Impacts on the Potential Distribution of Wild Radish in East Asia

Climate change can exert a considerable influence on the geographic distribution of many taxa, including coastal plants and populations of some plant species closely related to those used as agricultural crops. East Asian wild radish, Raphanus raphanistrum subsp. sativus, is an annual coastal plant that is a wild relative of the cultivated radish (R. sativus). It has served as source of genetic material that has been helpful to develop and improve the quality and yield of radish crops. To assess the impact of climate change on wild radish in East Asia, we analyzed its distribution at different periods using the maximum entropy model (MaxEnt). The results indicated that the precipitation of the driest month (bio14) and precipitation seasonality (bio15) were the two most dominant environmental factors that affected the geographical distribution of wild radish in East Asia. The total potential area suitable for wild radish is 102.5574 × 104 km2, mainly located along the seacoasts of southern China, Korea, and the Japanese archipelago. Compared with its current distribution regions, the potentially suitable areas for wild radish in the 2070s will further increase and expand northwards in Japan, especially on the sand beach habitats of Hokkaido. This research reveals the spatiotemporal changes for the coastal plant wild radish under global warming and simultaneously provides a vital scientific basis for effective utilization and germplasm innovation for radish cultivars to achieve sustainable agriculture development.

Population structure and phylogeography of three closely related tree peonies

Paeonia decomposita, Paeonia rotundiloba, and Paeonia rockii are three closely related species of Sect. Moutan is distributed in the montane area of the Eastern Hengduan Mountain region. Understanding the population history of these three tree peony species could contribute to unraveling the evolutionary patterns of undergrowth species in this hotspot area. We used one nuclear DNA marker (internal transcribed spacer region, ITS) and two chloroplast DNA markers (matK, ycf1) to reconstruct the phylogeographic pattern of the populations. In total, 228 individuals from 17 populations of the three species were analyzed in this study. Three nuclear clades (Clade I - Clade III) and four maternal clades (Clade A - Clade D) were reconstructed. Molecular dating suggested that young lineages diverged during the late Pliocene and early Pleistocene, younger than the uplift of the Hengduan Mountains but older than the last glacial maximum (LGM). Significant population and phylogeographic structures were detected at both markers. Furthermore, the populations of these tree peonies were overall at equilibrium during the climatic oscillations of the Pleistocene. The simulated palaeoranges of the three species during the LGM period mostly overlapped, which could have led to cross-breeding events. We propose an evolutionary scenario in which mountain orogenesis around the Hengduan Mountain area triggered parapatric isolation between maternal lineages of tree peonies. Subsequent climatic fluctuations drove migration and range recontact of these populations along the valleys. This detailed evolutionary history provides new insights into the phylogeographic pattern of species from mountain-valley systems.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

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.

A clue to the evolutionary history of modern East Asian flora: insights from phylogeography and diterpenoid alkaloid distribution pattern of the Spiraea japonica complex

Each subkingdom of East Asian flora (EAF) has a unique evolutionary history, but which has rarely been described based on phylogeographic studies of EAF species. The Spiraea japonica L. complex, which is widespread in East Asia (EA), has received considerable attention because of the presence of diterpenoid alkaloids (DAs). It provides a proxy for understanding the genetic diversity and DA distribution patterns of species under various environmental conditions associated with the geological background in EA. In the present study, the plastome and chloroplast/nuclear DNA of 71 populations belonging to the S. japonica complex and its congeners were sequenced, combined with DA identification, environmental analyses, and ecological niche modelling, to investigate their phylogenetic relationships, genetic and DAs distribution patterns, biogeography, and demographic dynamics. An "ampliative" S. japonica complex was put forward, comprising all species of Sect. Calospira Ser. Japonicae, of which three evolutionary units carrying their respective unique types of DAs were identified and associated with the regionalization of EAF (referring to the Hengduan Mountains, central China, and east China). Moreover, a transition belt in central China with its biogeographic significance was revealed by genetic and DA distribution patterns from the perspective of ecological adaptation. The origin and onset differentiation of the "ampliative" S. japonica complex was estimated in the early Miocene (22.01/19.44 Ma). The formation of Japanese populations (6.75 Ma) was facilitated by the land bridge, which subsequently had a fairly stable demographic history. The populations in east China have undergone a founder effect after the Last Glacial Maximum, which may have been promoted by the expansion potential of polyploidization. Overall, the in-situ origin and diversification of the "ampliative" S. japonica complex since the early Miocene is a vertical section of the formation and development of modern EAF and was shaped by the geological history of each subkingdom.

Allopatric Lineage Divergence of the East Asian Endemic Herb Conandron ramondioides Inferred from Low-Copy Nuclear and Plastid Markers

The evolutionary histories of ornamental plants have been receiving only limited attention. We examined the origin and divergence processes of an East Asian endemic ornamental plant, Conandron ramondioides. C. ramondioides is an understory herb occurring in primary forests, which has been grouped into two varieties. We reconstructed the evolutionary and population demography history of C. ramondioides to infer its divergence process. Nuclear and chloroplast DNA sequences were obtained from 21 Conandron populations on both sides of the East China Sea (ECS) to explore its genetic diversity, structure, and population differentiation. Interestingly, the reconstructed phylogeny indicated that the populations should be classified into three clades corresponding to geographical regions: the Japan (Honshu+Shikoku) clade, the Taiwan-Iriomote clade, and the Southeast China clade. Lineage divergence between the Japan clade and the Taiwan-Iriomote and Southeast China clades occured 1.14 MYA (95% HPD: 0.82-3.86), followed by divergence between the Taiwan-Iriomote and Southeast China clades approximately 0.75 MYA (95% HPD: 0.45-1.3). Furthermore, corolla traits (floral lobe length to tube length ratios) correlated with geographical distributions. Moreover, restricted gene flow was detected among clades. Lastly, the lack of potential dispersal routes across an exposed ECS seafloor during the last glacial maximum suggests that migration among the Conandron clades was unlikely. In summary, the extant Conandron exhibits a disjunct distribution pattern as a result of vicariance rather than long-distance dispersal. We propose that allopatric divergence has occurred in C. ramondioides since the Pleistocene. Our findings highlight the critical influence of species' biological characteristics on shaping lineage diversification of East Asian relic herb species during climate oscillations since the Quaternary.

Comparative Chloroplast Genome Analyses of the Winter-Blooming Eastern Asian Endemic Genus Chimonanthus (Calycanthaceae) With Implications For Its Phylogeny and Diversification

Chimonanthus of Calycanthaceae is a small endemic genus in China, with unusual winter-blooming sweet flowers widely cultivated for ornamentals and medicinal uses. The evolution of Chimonanthus plastomes and its phylogenetic relationships remain unresolved due to limited availability of genetic resources. Here, we report fully assembled and annotated chloroplast genomes of five Chimonanthus species. The chloroplast genomes of the genus (size range 153,010 – 153,299 bp) reveal high similarities in gene content, gene order, GC content, codon usage, amino acid frequency, simple sequence repeats, oligonucleotide repeats, synonymous and non-synonymous substitutions, and transition and transversion substitutions. Signatures of positive selection are detected in atpF and rpoB genes in C. campanulatus. The correlations among substitutions, InDels, and oligonucleotide repeats reveal weak to strong correlations in distantly related species at the intergeneric levels, and very weak to weak correlations among closely related Chimonanthus species. Chloroplast genomes are used to reconstruct a well-resolved phylogenetic tree, which supports the monophyly of Chimonanthus. Within Chimonanthus, C. praecox and C. campanulatus form one clade, while C. grammatus, C. salicifolius, C. zhejiangensis, and C. nitens constitute another clade. Chimonanthus nitens appears paraphyletic and is closely related to C. salicifolius and C. zhejiangensis, suggesting the need to reevaluate the species delimitation of C. nitens. Chimonanthus and Calycanthus diverged in mid-Oligocene; the radiation of extant Chimonanthus species was dated to the mid-Miocene, while C. grammatus diverged from other Chimonanthus species in the late Miocene. C. salicifolius, C. nitens(a), and C. zhejiangensis are inferred to have diverged in the Pleistocene of the Quaternary period, suggesting recent speciation of a relict lineage in the subtropical forest regions in eastern China. This study provides important insights into the chloroplast genome features and evolutionary history of Chimonanthus and family Calycanthaceae.

Spatial dynamics of Chinese Muntjac related to past and future climate fluctuations

Climate fluctuations in the past and in the future are likely to result in population expansions, shifts, or the contraction of the ecological niche of many species, and potentially leading to the changes in their geographical distributions. Prediction of suitable habitats has been developed as a useful tool for the assessment of habitat suitability and resource conservation to protect wildlife. Here, we model the ancestral demographic history of the extant modern Chinese Muntjac Muntiacus reevesi populations using approximate Bayesian computation (ABC) and used the maximum entropy model to simulate the past and predict the future spatial dynamics of the species under climate oscillations. Our results indicated that the suitable habitats for the M. reevesi shifted to the Southeast and contracted during the Last Glacial Maximum, whereas they covered a broader and more northern position in the Middle Holocene. The ABC analyses revealed that the modern M. reevesi populations diverged in the Middle Holocene coinciding with the significant contraction of the highly suitable habitat areas. Furthermore, our predictions suggest that the potentially suitable environment distribution for the species will expand under all future climate scenarios. These results indicated that the M. reevesi diverged in the recent time after the glacial period and simultaneously as its habitat’s expanded in the Middle Holocene. Furthermore, the past and future climate fluctuation triggered the change of Chinese muntjac spatial distribution, which has great influence on the Chinese muntjac’s population demographic history.

Cytotype distribution and chloroplast phylogeography of the Actinidia chinensis complex

BACKGROUND

Plant phylogeographic studies of species in subtropical China have mainly focused on rare and endangered species, whereas few studies have been conducted on taxa with relatively wide distribution, especially polyploid species. We investigated the cytotype and haplotype distribution pattern of the Actinidia chinensis complex, a widespread geographically woody liana with variable ploidy in subtropical China comprising two varieties, with three chloroplast fragments DNA (ndhF-rpl132, rps16-trnQ and trnE-trnT). Macroevolutionary, microevolutionary and niche modeling tools were also combined to disentangle the origin and the demographic history of the species or cytotypes.

RESULTS

The ploidy levels of 3338 individuals from 128 populations sampled throughout the species distribution range were estimated with flow cytometry. The widespread cytotypes were diploids followed by tetraploids and hexaploids, whereas triploids and octoploids occurred in a few populations. Thirty-one chloroplast haplotypes were detected. The genetic diversity and genetic structure were found to be high between varieties (or ploidy races) chinensis and deliciosa. Our results revealed that these two varieties inhabit significantly different climatic niche spaces. Ecological niche models (ENMs) indicate that all varieties' ranges contracted during the Last Inter Glacial (LIG), and expanded eastward or northward during the Last Glacial Maximum (LGM).

CONCLUSIONS

Pliocene and Plio-Pleistocene climatic fluctuations and vicariance appear to have played key roles in shaping current population structure and historical demography in the A. chinensis complex. The polyploidization process also appears to have played an important role in the historical demography of the complex through improving their adaptability to environmental changes.

Inferring historical survivals of climate relicts: the effects of climate changes, geography, and population-specific factors on herbaceous hydrangeas

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7-9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137-0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.

Climatic Refugia and Geographical Isolation Contribute to the Speciation and Genetic Divergence in Himalayan-Hengduan Tree Peonies (Paeonia delavayi and Paeonia ludlowii)

Himalaya and Hengduan Mountains (HHM) is a biodiversity hotspot, and very rich in endemic species. Previous phylogeographical studies proposed different hypotheses (vicariance and climate-driven speciation) in explaining diversification and the observed pattern of extant biodiversity, but it is likely that taxa are forming in this area in species-specific ways. Here, we reexplored the phylogenetic relationship and tested the corresponding hypotheses within Paeonia subsect. Delavayanae composed of one widespread species (Paeonia delavayi) and the other geographically confined species (Paeonia ludlowii). We gathered genetic variation data at three chloroplast DNA fragments and one nuclear gene from 335 individuals of 34 populations sampled from HHM. We performed a combination of population genetic summary statistics, isolation-with-migration divergence models, isolation by environment, and demographic history analyses. We found evidence for the current taxonomic treatment that P. ludlowii and P. delavayi are two different species with significant genetic differentiation. The significant isolation by environment was revealed within all sampled populations but genetic distances only explained by geographical distances within P. delavayi populations. The results of population divergence models and demographic history analyses indicated a progenitor–derivative relationship and the Late Quaternary divergence without gene flow between them. The coalescence of all sampled cpDNA haplotypes could date to the Late Miocene, and P. delavayi populations probably underwent a severe bottleneck in population size during the last glacial period. Genetic variation in Paeonia subsect. Delavayanae is associated with geographical and environmental distances. These findings point to the importance of geological and climatic changes as causes of the speciation event and lineage diversification within Paeonia subsect. Delavayanae.

Adaptive strategies and driving factors of a montane riparian tree: Trait-specific mechanisms across latitude

Investigating the drivers of phenotypic and genetic divergence can reveal the underlying processes and strategies that species adopt in rapidly changing environments. However, knowledge of adaptive strategies and the underlying mechanisms is lacking for the majority of taxa, especially those living in habitats sensitive to climate change. Here, we investigated 20 populations of a Tertiary-relict tree species, Euptelea pleiospermum (Eupteleaceae), scattered in a mountain riparian habitat in China. We integrated genetic, growth, and reproductive traits, and evaluated the relative contributions of climatic and soil factors on genetic and functional trait divergence. The E. pleiospermum populations were divided into south and north genetic clusters, and there were significant differences in leaf density and seed mass of adult trees between the two. The spatial pattern of genetic divergence resulted from effects of both isolation by distance (IBD) and isolation by environment (IBE), whereas the divergence of growth and reproductive traits resulted solely from IBE effects. Spatial distance and selection by temperature and soils played dominant roles in genetic divergence. Precipitation drove the spatial divergence of sprouting. Both divergence of leaf density and seed mass were prominently induced by genetic divergence, and the influences might be enhanced by temperature and soil nutrients. We infer that E. pleiospermum populations adopt a resource-conservative strategy with low growth rates and higher sprouting under flooding disturbance, with larger seeds for improved seedling recruitment at lower latitudes. In contrast, high growth rate and sexual reproduction with small seeds are strategies adopted by populations at higher latitudes. We conclude that sprouting reflects a plastic response to precipitation, and leaf density and seed mass reflect local adaption under selection by temperature and soil factors. The underlying mechanisms of species adaptation strategies were trait-specific. Temperature and soil conditions are likely the main ecological factors shaping plant divergence in montane riparian regions.

The Tetracentron genome provides insight into the early evolution of eudicots and the formation of vessel elements

Background

Tetracentron sinense is an endemic and endangered deciduous tree. It belongs to the Trochodendrales, one of four early diverging lineages of eudicots known for having vesselless secondary wood. Sequencing and resequencing of the T. sinense genome will help us understand eudicot evolution, the genetic basis of tracheary element development, and the genetic diversity of this relict species.

Results

Here, we report a chromosome-scale assembly of the T. sinense genome. We assemble the 1.07 Gb genome sequence into 24 chromosomes and annotate 32,690 protein-coding genes. Phylogenomic analyses verify that the Trochodendrales and core eudicots are sister lineages and showed that two whole-genome duplications occurred in the Trochodendrales approximately 82 and 59 million years ago. Synteny analyses suggest that the γ event, resulting in paleohexaploidy, may have only happened in core eudicots. Interestingly, we find that vessel elements are present in T. sinense, which has two orthologs of AtVND7, the master regulator of vessel formation. T. sinense also has several key genes regulated by or regulating TsVND7.2 and their regulatory relationship resembles that in Arabidopsis thaliana. Resequencing and population genomics reveals high levels of genetic diversity of T. sinense and identifies four refugia in China.

Conclusions

The T. sinense genome provides a unique reference for inferring the early evolution of eudicots and the mechanisms underlying vessel element formation. Population genomics analysis of T. sinense reveals its genetic diversity and geographic structure with implications for conservation.

Supplementary information

The online version contains supplementary material available at 10.1186/s13059-020-02198-7.

Genomic insights on the contribution of balancing selection and local adaptation to the long-term survival of a widespread living fossil tree, Cercidiphyllum japonicum

'Living fossils' are testimonies of long-term sustained ecological success, but how demographic history and natural selection contributed to their survival, resilience, and persistence in the face of Quaternary climate fluctuations remains unclear. To better understand the interplay between demographic history and selection in shaping genomic diversity and evolution of such organisms, we assembled the whole genome of Cercidiphyllum japonicum, a widespread East Asian Tertiary relict tree, and resequenced 99 individuals of C. japonicum and its sister species, Cercidiphyllum magnificum (Central Japan). We dated this speciation event to the mid-Miocene, and the intraspecific lineage divergence of C. japonicum (China vs Japan) to the Early Pliocene. Throughout climatic upheavals of the late Tertiary/Quaternary, population bottlenecks greatly reduced the genetic diversity of C. japonicum. However, this polymorphism loss was likely counteracted by, first, long-term balancing selection at multiple chromosomal and heterozygous gene regions, potentially reflecting overdominance, and, second, selective sweeps at stress response and growth-related genes likely involved in local adaptation. Our findings contribute to a better understanding of how living fossils have survived climatic upheaval and maintained an extensive geographic range; that is, both types of selection could be major factors contributing to the species' survival, resilience, and persistence.

Phylogeny, Taxonomy, and Biogeography of Pterocarya (Juglandaceae)

Relict species play an important role in understanding the biogeography of intercontinental disjunctions. Pterocarya (a relict genus) is the valuable model taxon for studying the biogeography of East Asian versus southern European/West Asian disjunct patterns. This disjunction has not been as well studied as others (e.g., between Eastern Asia and North America). Several phylogenetic studies on Pterocarya have been conducted, but none have provided a satisfactory phylogenetic resolution. Here, we report the first well-resolved phylogeny of Pterocarya using restriction site-associated DNA sequencing data based on the sampling of all taxa across the entire distribution area of the genus. Taxonomic treatments were also clarified by combining morphological traits. Furthermore, fossil-calibrated phylogeny was used to explore the biogeography of Pterocarya. Our results support the existence of two sections in Pterocarya, which is in accordance with morphological taxonomy. Section Platyptera comprises three species: P. rhoifolia, P. macroptera, and P. delavayi. Section Pterocarya also comprises three species: P. fraxinifolia, P. hupehensis, and P. stenoptera. The divergence between the two sections took place during the early Miocene (20.5 Ma). The formation of the Gobi Desert and climate cooling of northern Siberia in the Middle Miocene (15.7 Ma) might have caused the split of the continuous distribution of this genus and the formation of the East Asian versus southern European/West Asian disjunct pattern. Lastly, the divergence between P. hupehensis and P. stenoptera as well as between P. rhoifolia and P. macroptera/P. delavayi (10.0 Ma) supports the late Miocene diversification hypothesis in East Asia.

Genomic analysis reveals the genetic diversity, population structure, evolutionary history and relationships of Chinese pepper

Chinese pepper, mainly including Zanthoxylum bungeanum and Zanthoxylum armatum, is an economically important crop popular in Asian countries due to its unique taste characteristics and potential medical uses. Numerous cultivars of Chinese pepper have been developed in China through long-term domestication. To better understand the population structure, demographic history, and speciation of Chinese pepper, we performed a comprehensive analysis at a genome-wide level by analyzing 38,395 genomic SNPs that were identified in 112 cultivated and wild accessions using a high-throughput genome-wide genotyping-by-sequencing (GBS) approach. Our analysis provides genetic evidence of multiple splitting events occurring between and within species, resulting in at least four clades in Z. bungeanum and two clades in Z. armatum. Despite no evidence of recent admixture between species, we detected substantial gene flow within species. Estimates of demographic dynamics and species distribution modeling suggest that climatic oscillations during the Pleistocene (including the Penultimate Glaciation and the Last Glacial Maximum) and recent domestication events together shaped the demography and evolution of Chinese pepper. Our analyses also suggest that southeastern Gansu province is the most likely origin of Z. bungeanum in China. These findings provide comprehensive insights into genetic diversity, population structure, demography, and adaptation in Zanthoxylum.

Genomic insights into historical population dynamics, local adaptation, and climate change vulnerability of the East Asian Tertiary relict Euptelea (Eupteleaceae)

The warm-temperate and subtropical climate zones of East Asia are a hotspot of plant species richness and endemism, including a noticeable number of species-poor Tertiary relict tree genera. However, little is understood about when East Asian Tertiary relict plants diversified, how they responded demographically to past environmental change, and to what extent their current genomic composition (and adaptive capacity) might mitigate the effects of global warming. Here, we obtained genomic (RAD-SNP) data for 171 samples from two extant species of Euptelea in China (24 E. pleiosperma populations) and Japan (11 E. polyandra populations) to elucidate their divergence and demographic histories, genome-wide associations with current environmental variables, and genomic vulnerability to future climate change. Our results indicate that Late Miocene changes in climate and/or sea level promoted species divergence, whereas Late Pliocene uplifting in southwest China likely fostered lineage divergence within E. pleiosperma. Its subsequent range expansion into central/east (CE) China bears genomic signatures of climate-driven selection, yet extant CE populations are predicted to be most vulnerable to future climate change. For E. polyandra, geography was the only significant predictor of genomic variation. Our findings indicate a profound impact of Late Neogene geological and climate change on the evolutionary history of Euptelea, with much stronger signals of local adaptation left in China than in Japan. This study deepens our understanding of the complex evolutionary forces that influence the distribution of genetic variation of Tertiary relict trees, and provides insights into their susceptibility to global change and potential for adaptive responses. Our results lay the groundwork for future conservation and restoration programs for Euptelea.

Dynamics behind disjunct distribution, hotspot-edge refugia, and discordant RADseq/mtDNA variability: insights from the Emei mustache toad

Background

The distribution of genetic diversity and the underlying processes are important for conservation planning but are unknown for most species and have not been well studied in many regions. In East Asia, the Sichuan Basin and surrounding mountains constitute an understudied region that exhibits a “ring” of high species richness overlapping the eastern edge of the global biodiversity hotspot Mountains of Southwest China. We examine the distributional history and genetic diversification of the Emei mustache toad Leptobrachium boringii, a typical “ring” element characterized by disjunct ranges in the mountains, by integrating time-calibrated gene tree, genetic variability, individual-level clustering, inference of population splitting and mixing from allele frequencies, and paleoclimatic suitability modeling.

Results

The results reveal extensive range dynamics, including secondary contact after long-term isolation via westward dispersal accompanied by variability loss. They allow the proposal of a model that combines recurrent contractions caused by Quaternary climatic changes and some failed expansions under suitable conditions for explaining the shared disjunct distribution pattern. Providing exceptional low-elevation habitats in the hotspot area, the eastern edge harbors both long-term refugial and young immigrant populations. This finding and a synthesis of evidence from other taxa demonstrate that a certain contributor to biodiversity, one that preserves and receives low-elevation elements of the east in this case, can be significant for only a particular part of a hotspot. By clarifying the low variability of these refugial populations, we show that discordant mitochondrial estimates of diversity can be obtained for populations that experienced admixture, which would have unlikely left proportional immigrant alleles for each locus.

Conclusions

Dispersal after long-term isolation can explain much of the spatial distribution of genetic diversity in this species, while secondary contact and long-term persistence do not guarantee a large variation. The model for the formation of disjunct ranges may apply to many other taxa isolated in the mountains surrounding the Sichuan Basin. Furthermore, this study provides insights into the heterogeneous nature of hotspots and discordant variability obtained from genome-wide and mitochondrial data.

Miocene diversification of a golden-thread nanmu tree species (Phoebe zhennan, Lauraceae) around the Sichuan Basin shaped by the East Asian monsoon

Understanding the role of climate changes and geography as drivers of population divergence and speciation is a long‐standing goal of evolutionary biology and can inform conservation. In this study, we used restriction site‐associated DNA sequencing (RAD‐seq) to evaluate genetic diversity, population structure, and infer demographic history of the endangered tree, Phoebe zhennan which is distributed around the Sichuan Basin. Genomic patterns revealed two distinct clusters, each largely confined to the West and East. Despite sympatry of the two genomic clusters at some sites, individuals show little or no evidence of genomic introgression. Demographic modeling supported an initial divergence time between the West and East lineages at ~15.08 Ma with further diversification within the West lineage at ~7.12 Ma. These times largely coincide with the two independent intensifications of the East Asian monsoon that were initiated during the middle (Langhian) and late Miocene (Messinian), respectively. These results suggest that the Miocene intensification phases of the East Asian monsoon played a pivotal role in shaping the current landscape‐level patterns of genetic diversity within P. zhennan, as has been found for the interspecific divergence of other subtropical Chinese plants. Based on isolation‐by‐distance and species distribution modeling, we hypothesize that P. zhennan followed a ring diversification which was facilitated by the Sichuan Basin acting as barrier to gene flow. In situ and ex situ conservation management plans should consider the results obtained in this study to help secure the future of this beautiful and culturally significant endangered tree.

Independent allopatric polyploidizations shaped the geographical structure and initial stage of reproductive isolation in an allotetraploid fern, Lepisorus nigripes (Polypodiaceae)

Although polyploidy is pervasive and its evolutionary significance has been recognized, it remains unclear how newly formed polyploid species become established. In particular, the impact of multiple origins on genetic differentiation among populations of a polyploid species and whether lineages of independent origins have different evolutionary potentials remain open questions. We used population genetic and phylogenetic approaches to identify genetic differentiation between lineages with independent origins within an allotetraploid fern, Lepisorus nigripes. A total of 352 individuals from 51 populations were collected throughout the distribution range. To examine the genetic structure, multilocus genotyping, Bayesian population structure analysis, and neighbor-net analysis were carried out using single-copy nuclear genes. Phylogenetic trees were constructed to detect recurrent polyploid origins. Proportions of abortive spores were analysed as the measure of postzygotic reproductive isolation. Two genetically distinct lineages, the East-type and the West-type, were distributed mainly in the eastern and western parts, respectively, of the Japanese archipelago. Phylogenetic analyses indicated independent origins of these types and detected additional independent origins within each type. We also revealed limited genetic recombination between both types, even in their sympatric regions. F₁ hybrids between the East- and West-types showed a reduction in fertility. It is likely that the East- and West-types formed independently in the eastern and western parts of Japan, respectively. The limited genetic recombination and reduced fertility of hybrids suggest that the two types are at an incipient stage of speciation. Two polyploid lineages with independent geographic origins could develop reproductive isolation barrier(s).

Sky islands as foci for divergence of fig trees and their pollinators in southwest China

The dynamics of populations and their divergence over time have shaped current levels of biodiversity and in the case of the "sky islands" of mountainous southwest (SW) China have resulted in an area of exceptional botanical diversity. Ficus tikoua is a prostrate fig tree subendemic to the area that displays unique intraspecific diversity, producing figs typical of different pollination modes in different parts of its range. By combining climate models, genetic variation in populations of the tree's obligate fig wasp pollinators and distributions of the different plant phenotypes, we examined how this unusual situation may have developed. We identified three genetically distinct groups of a single Ceratosolen pollinator species that have largely parapatric distributions. The complex topography of the region contributed to genetic divergence among the pollinators by facilitating geographical isolation and providing refugia. Migration along elevations in response to climate oscillations further enhanced genetic differentiation of the three pollinator groups. Their distributions loosely correspond to the distributions of the functionally significant morphological differences in the male figs of their host plants, but postglacial expansion of one group has not been matched by spread of its associated plant phenotype, possibly due to a major river barrier. The results highlight how interplay between the complex topography of the "sky island" complex and climate change has shaped intraspecies differentiation and relationships between the plant and its pollinator. Similar processes may explain the exceptional botanical diversity of SW China.

Evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae), in south-east China during the late Neogene: old lineage, young populations

BACKGROUND AND AIMS

Many monotypic gymnosperm lineages in south-east China paradoxically remain in relict status despite long evolutionary histories and ample opportunities for allopatric speciation, but this paradox has received little attention and has yet to be resolved. Here, we address this issue by investigating the evolutionary history of a relict conifer, Pseudotaxus chienii (Taxaceae).

METHODS

DNA sequences from two chloroplast regions and 14 nuclear loci were obtained for 134 samples. The demographic history was inferred and the contribution of isolation by environment (IBE) in patterning genetic divergence was compared with that of isolation by distance (IBD).

KEY RESULTS

Three genetic clusters were identified. Approximate Bayesian computation analyses showed that the three clusters diverged in the late Pliocene (~3.68 Ma) and two admixture events were detected. Asymmetric gene flow and similar population divergence times (~ 3.74 Ma) were characterized using the isolation with migration model. Neither IBD nor IBE contributed significantly to genetic divergence, and the contribution of IBE was much smaller than that of IBD.

CONCLUSIONS

These results suggest that several monotypic relict gymnosperm lineages like P. chienii in south-east China did not remain in situ and undiversified for millions of years. On the contrary, they have been evolving and the extant populations have become established more recently, having insufficient time to speciate. Our findings provide a new perspective for understanding the formation and evolution of the relict gymnosperm flora of China as well as of the Sino-Japanese Flora.

Contrasting patterns of genetic structure and phylogeography in the marine agarophytes Gelidiophycus divaricatus and G. freshwateri (Gelidiales, Rhodophyta) from East Asia

The evolutionary and population demographic history of marine red algae in East Asia is poorly understood. Here, we reconstructed the phylogeographies of two upper intertidal species endemic to East Asia, Gelidiophycus divaricatus and G. freshwateri. Phylogenetic and phylogeographic inferences of 393 mitochondrial cox1, 128 plastid rbcL, and 342 nuclear ITS2 sequences were complemented with ecological niche models. Gelidiophycus divaricatus, a southern species adapted to warm water, is characterized by a high genetic diversity and a strong geographical population structure, characteristic of stable population sizes and sudden reduction to recent expansion. In contrast, G. freshwateri, a northern species adapted to cold temperate conditions, is genetically relatively homogeneous with a shallow population structure resulting from steady population growth and recent equilibrium. The overlap zone of the two species roughly matches summer and winter isotherms, indicating that surface seawater temperature is a key feature influencing species range. Unidirectional genetic introgression was detected at two sites on Jeju Island where G. divaricatus was rare while G. freshwateri was common, suggesting the occurrence of asymmetric natural hybrids, a rarely reported event for rhodophytes. Our results illustrate that Quaternary climate oscillations have left strong imprints on the current day genetic structure and highlight the importance of seawater temperature and sea level change in driving speciation in upper intertidal seaweed species.

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.

Biodiversity can benefit from climate stabilization despite adverse side effects of land-based mitigation

Limiting the magnitude of climate change via stringent greenhouse gas (GHG) mitigation is necessary to prevent further biodiversity loss. However, some strategies to mitigate GHG emission involve greater land-based mitigation efforts, which may cause biodiversity loss from land-use changes. Here we estimate how climate and land-based mitigation efforts interact with global biodiversity by using an integrated assessment model framework to project potential habitat for five major taxonomic groups. We find that stringent GHG mitigation can generally bring a net benefit to global biodiversity even if land-based mitigation is adopted. This trend is strengthened in the latter half of this century. In contrast, some regions projected to experience much growth in land-based mitigation efforts (i.e., Europe and Oceania) are expected to suffer biodiversity loss. Our results support the enactment of stringent GHG mitigation policies in terms of biodiversity. To conserve local biodiversity, however, these policies must be carefully designed in conjunction with land-use regulations and societal transformation in order to minimize the conversion of natural habitats.

Greenhouse gas mitigation can involve land-use changes that alter the habitat available for wildlife. Here, Ohashi et al. perform an integrated assessment showing that climate mitigation can be beneficial for global biodiversity but may entail local biodiversity losses where land-based mitigation is implemented.

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.

Multiple origins and the population genetic structure of Rubus takesimensis (Rosaceae) on Ulleung Island: Implications for the genetic consequences of anagenetic speciation

To determine the origin and genetic consequences of anagenesis in Rubus takesimensis on Ulleung Island, Korea, we compared the genetic diversity and population structure of R. takesimensis with those of its continental progenitor R. crataegifolius. We broadly sampled a total of 315 accessions in 35 populations and sequenced five noncoding regions of chloroplast DNA. Rubus takesimensis emerged as nonmonophyletic and several geographically diverse continental populations were likely responsible for the origin of R. takesimensis; the majority of R. takesimensis accessions were sisters to the clade containing accessions of R. crataegifolius, primarily from the Korean peninsula, while rare accessions from three populations shared common ancestors with the ones from the southern part of the Korean peninsula, Jeju Island, and Japan. A few accessions from the Chusan population originated independently from the Korean peninsula. Of 129 haplotypes, 81 and 48 were found exclusively in R. crataegifolius and R. takesimensis, respectively. We found unusually high genetic diversity in two regions on Ulleung Island and no geographic population structure. For R. crataegifolius, two major haplotype groups were found; one for the northern mainland Korean peninsula, and the other for the southern Korean peninsula and the Japanese archipelago. Compared with populations of R. crataegifolius sampled from Japan, much higher haplotype diversity was found in populations from the Korean peninsula. The patterns of genetic consequences in R. takesimensis need to be verified for other endemic species based on chloroplast DNA and independent nuclear markers to synthesize emerging patterns of anagenetic speciation on Ulleung Island.

Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil

As Charles Darwin anticipated, living fossils provide excellent opportunities to study evolutionary questions related to extinction, competition, and adaptation. Ginkgo (Ginkgo biloba L.) is one of the oldest living plants and a fascinating example of how people have saved a species from extinction and assisted its resurgence. By resequencing 545 genomes of ginkgo trees sampled from 51 populations across the world, we identify three refugia in China and detect multiple cycles of population expansion and reduction along with glacial admixture between relict populations in the southwestern and southern refugia. We demonstrate multiple anthropogenic introductions of ginkgo from eastern China into different continents. Further analyses reveal bioclimatic variables that have affected the geographic distribution of ginkgo and the role of natural selection in ginkgo’s adaptation and resilience. These investigations provide insights into the evolutionary history of ginkgo trees and valuable genomic resources for further addressing various questions involving living fossil species.

Ginkgo is one of the living fossils from the plant kingdom. Here, authors conduct population genomics analyses to reveal its refugia and demographic history, and provide evidence of multiple anthropogenic introductions of ginkgo from eastern China into different continents.

Phylogeographic Analysis and Genetic Structure of an Endemic Sino-Japanese Disjunctive Genus Diabelia (Caprifoliaceae)

The Sino-Japanese Floristic Region (SJFR) is a key area for plant phylogeographical research, due to its very high species diversity and disjunct distributions of a large number of species and genera. At present, the root cause and temporal origin of the discontinuous distribution of many plants in the Sino-Japanese flora are still unclear. Diabelia (Caprifoliaceae; Linnaeoideae) is a genus endemic to Asia, mostly in Japan, but two recent discoveries in China raised questions over the role of the East China Sea (ECS) in these species' disjunctions. Chloroplast DNA sequence data were generated from 402 population samples for two regions (rpl32-trnL, and trnH-psbA) and 11 nuclear microsatellite loci were screened for 549 individuals. Haplotype, population-level structure, combined analyses of ecological niche modeling, and reconstruction of ancestral state in phylogenies were also performed. During the Last Glacial Maximum (LGM) period after the Tertiary, Diabelia was potentially widely distributed in southeastern China, the continental shelf of the East China Sea and Japan (excluding Hokkaido). After LGM, all populations in China have disappeared except those in Zhejiang which may represent a Glacial refuge. Populations of Diabelia in Japan have not experienced significant bottleneck effects, and populations have maintained a relatively stable state. The observed discontinuous distribution of Diabelia species between China and Japan are interpreted as the result of relatively ancient divergence. The phylogenetic tree of chloroplast fragments shows the characteristics of multi-origin evolution (except for D. sanguinea). STRUCTURE analysis of nuclear Simple Sequence Repeat (nSSR) showed that the plants of the Diabelia were divided into five gene pools: D. serrata, D. spathulata, D. sanguinea, D. ionostachya (D. spathulata var. spathulata-Korea), and populations of D. ionostachya var. ionostachya in Yamagata prefecture, northern Japan. Molecular evidence provides new insights of Diabelia into biogeography, a potential glacial refuge, and population-level genetic structure within species. In the process of species differentiation, ECS acts as a corridor for two-way migration of animals and plants between China and Japan during glacial maxima, providing the possibility of secondary contact for discontinuously distributed species between China and Japan, or as a filter (creating isolation) during glacial minima. The influence of the ECS in speciation and biogeography of Diabelia in the Tertiary remains unresolved in this study. Understanding origins, evolutionary histories, and speciation will provide a framework for the conservation and cultivation of Diabelia.

Evolutionary History of Rhus chinensis (Anacardiaceae) From the Temperate and Subtropical Zones of China Based on cpDNA and Nuclear DNA Sequences and Ecological Niche Model

To explore the origin and evolution of local flora and vegetation, we examined the evolutionary history of Rhus chinensis, which is widely distributed in China’s temperate and subtropical zones, by sequencing three maternally inherited chloroplast DNAs (cpDNA: trnL-trnF, psbA-trnH, and rbcL) and the biparentally inherited nuclear DNA (nuDNA: LEAFY) from 19 natural populations of R. chinensis as well as the ecological niche modeling. In all, 23 chloroplast haplotypes (M1–M23) and 15 nuclear alleles (N1–N15) were detected. The estimation of divergence time showed that the most recent common ancestor dated at 4.2 ± 2.5 million years ago (Mya) from cpDNA, and the initial divergence of genotypes occurred at 4.8 ± 3.6 Mya for the nuDNA. Meanwhile, the multimodality mismatch distribution curves and positive Tajima’s D values indicated that R. chinensis did not experience population expansion after the last glacial maximum. Besides, our study was also consistent with the hypothesis that most refugia in the temperate and subtropical zones of China were in situ during the glaciation.

New insight into the phylogeographic pattern of Liriodendron chinense (Magnoliaceae) revealed by chloroplast DNA: east-west lineage split and genetic mixture within western subtropical China

Background

Subtropical China is a global center of biodiversity and one of the most important refugia worldwide. Mountains play an important role in conserving the genetic resources of species. Liriodendron chinense is a Tertiary relict tree largely endemic to subtropical China. In this study, we aimed to achieve a better understanding of the phylogeographical pattern of L. chinense and to explore the role of mountains in the conservation of L. chinense genetic resources.

Methods

Three chloroplast regions (psbJ-petA, rpl32-ndhF, and trnK5’-matK) were sequenced in 40 populations of L. chinense for phylogeographical analyses. Relationships among chloroplast DNA (cpDNA) haplotypes were determined using median-joining networks, and genetic structure was examined by spatial analysis of molecular variance (SAMOVA). The ancestral area of the species was reconstructed using the Bayesian binary Markov Chain Monte Carlo (BBM) method according to its geographic distribution and a maximum parsimony (MP) tree based on Bayesian methods.

Results

Obvious phylogeographic structure was found in L. chinense. SAMOVA revealed seven groups matching the major landscape features of the L. chinense distribution area. The haplotype network showed three clades distributed in the eastern, southwestern, and northwestern regions. Separate northern and southern refugia were found in the Wu Mountains and Yungui Plateau, with genetic admixture in the Dalou Mountains and Wuling Mountains. BBM revealed a more ancient origin of L. chinense in the eastern region, with a west–east split most likely having occurred during the Mindel glacial stage.

Discussion

The clear geographical distributions of haplotypes suggested multiple mountainous refugia of L. chinense. The east–west lineage split was most likely a process of gradual genetic isolation and allopatric lineage divergence when the Nanling corridor was frequently occupied by evergreen or coniferous forest during Late Quaternary oscillations. Hotspots of haplotype diversity in the Dalou Mountains and Wuling Mountains likely benefited from gene flow from the Wu Mountains and Yungui Plateau. Collectively, these results indicate that mountain regions should be the main units for conserving and collecting genetic resources of L. chinense and other similar species in subtropical China.

Evolutionary Hotspots of Seed Plants in Subtropical China: A Comparison With Species Diversity Hotspots of Woody Seed Plants

Genetic diversity is a fundamental level of biodiversity. However, it is frequently neglected in conservation prioritization because intraspecific genetic diversity is difficult to measure at large scales. In this study, we synthesized population genetic or phylogeographic datasets of 33 seed plants in subtropical China into multi-species genetic landscapes. The genetic landscapes identified 18 evolutionary hotspots with high within-population genetic diversity (WGD), and among-population genetic diversity (AGD), or both. The western subtropical China is rich in AGD (possessing four major AGD hotspots), deserving a high conservation priority. We found that WGD was positively correlated with longitude, with most WGD hotspots locating in east subtropical China. The results showed that the locations of 12 of 18 evolutionary hotspots corresponded approximately to those of previously identified species diversity (SD) hotspots, however, a positive and significant correlation existed only between AGD and SD, not between WGD and SD. Therefore, spatial patterns of species richness in plants in subtropical China cannot generally be used as surrogate for their intraspecific diversity. This study identified multi-species evolutionary hotspots and correlated multi-species genetic diversity with SD across subtropical China for the first time, providing profound implications for the conservation of biodiversity in this important ecoregion.

Genetic diversity and population structure of the endangered basal angiosperm Brasenia schreberi (Cabombaceae) in China

Brasenia schreberi J.F. Gmelin (Cabombaceae), an aquatic herb that occurs in fragmented locations in China, is rare and endangered. Understanding its genetic diversity and structure is crucial for its conservation and management. In this study, 12 microsatellite markers were used to estimate the genetic diversity and variation in 21 populations of B. schreberi in China. A total of 61 alleles were found; assessment of allelic richness (Ar = 1.92) and observed and expected heterozygosity (HO = 0.200, HE = 0.256) suggest lower genetic diversity compared to some endangered species, and higher variation was observed within populations (58.68%) rather than among populations (41.32%). No significant correlation between geographical and genetic distance among populations was detected (Mantel test, r = 0.0694; P = 0.7985), which may have likely resulted from barriers to gene flow (Nm = 0.361) that were produced by habitat fragmentation. However, Bayesian and neighbor-joining cluster analyses suggest a population genetic structure consisting of two clusters (I and II) or four subclusters (I-1, 2 and II-1, 2). The genetic structure and distribution of B. schreberi in China may have involved glacial refugia that underwent range expansions, introgression, and habitat fragmentation. The findings of the present study emphasize the importance for both in situ and ex situ conservation efforts.

Molecular phylogeography of East Asian Boea clarkeana (Gesneriaceae) in relation to habitat restriction

Subfamily Cyrtandroideae (Gesneriaceae) comprises a broadly distributed group of rocky-slope herbs, with China being the center of its distributional range. The normal growth of many species within the family is particularly dependent on special habitats. Due to the paucity of molecular studies, very little is known regarding East Asian herb phylogeographic pattern. Here, we investigate the molecular phylogeography of Boea clarkeana Hemsl., a unique resurrection herb endemic to China, focusing on geographically restrictive effects of habitat distribution on evolutionary history. Variation in three chloroplast DNA (cpDNA) intergenic spacers (psbA-trnH, rps12-rpl20, and trnL-trnF), the ribosomal internal transcribed spacer (ITS) and simple sequence repeats in expressed sequence tags (EST-SSRs) was investigated across 18 populations to assess genetic diversity, genetic structure and historical dynamics. Genetic diversity was low within populations (cpDNA, hS = 0.03, πS×10(3) = 0.17; ITS, hS = 0.16, πS×10(3) = 0.43) but high for species (cpDNA, hT = 0.82, πT×10(3) = 3.12; ITS, hT = 0.88, πT×10(3) = 6.39); 76 alleles were detected in this highly inbred species (FIS = 0.22), with a significantly low average of 1.34 alleles per locus. No cpDNA or ITS haplotypes were shared between regions. Based on cpDNA results, the Mt. Huangshan-Tianmu and Mt. Qinling-Daba haplotypes are ancestral; these two regions represent potential refugia. Although no evidence of significant retreat-migration phenomena during glacial cycles was detected, interglacial range expansion from northern Mt. Qinling-Daba was identified (121,457 yr BP). Rapid agricultural growth caused bottlenecks in many populations, especially on Mt. Huang-Tianmu. Habitat restriction and fragmentation, weak seed and pollen dispersal abilities, and long-term isolation caused by human-induced or environmental changes are considered the main causes of extinction of several populations and low genetic diversity within populations and regions. These analyses clarify the effects of habitat restriction on B. clarkeana, representing an evolutionary reference for similar gesneriads, and enrich our understanding of the molecular phylogeography of East Asian rocky-slope herbs.

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.

Ginkgo biloba's footprint of dynamic Pleistocene history dates back only 390,000 years ago

BACKGROUND

At the end of the Pliocene and the beginning of Pleistocene glaciation and deglaciation cycles Ginkgo biloba went extinct all over the world, and only few populations remained in China in relict areas serving as sanctuary for Tertiary relict trees. Yet the status of these regions as refuge areas with naturally existing populations has been proven not earlier than one decade ago. Herein we elaborated the hypothesis that during the Pleistocene cooling periods G. biloba expanded its distribution range in China repeatedly. Whole plastid genomes were sequenced, assembled and annotated, and sequence data was analyzed in a phylogenetic framework of the entire gymnosperms to establish a robust spatio-temporal framework for gymnosperms and in particular for G. biloba Pleistocene evolutionary history.

RESULTS

Using a phylogenetic approach, we identified that Ginkgoatae stem group age is about 325 million years, whereas crown group radiation of extant Ginkgo started not earlier than 390,000 years ago. During repeated warming phases, Gingko populations were separated and isolated by contraction of distribution range and retreated into mountainous regions serving as refuge for warm-temperate deciduous forests. Diversification and phylogenetic splits correlate with the onset of cooling phases when Ginkgo expanded its distribution range and gene pools merged.

CONCLUSIONS

Analysis of whole plastid genome sequence data representing the entire spatio-temporal genetic variation of wild extant Ginkgo populations revealed the deepest temporal footprint dating back to approximately 390,000 years ago. Present-day directional West-East admixture of genetic diversity is shown to be the result of pronounced effects of the last cooling period. Our evolutionary framework will serve as a conceptual roadmap for forthcoming genomic sequence data, which can then provide deep insights into the demographic history of Ginkgo.

Resolving relationships and phylogeographic history of the Nyssa sylvatica complex using data from RAD-seq and species distribution modeling

Nyssa sylvatica complex consists of several woody taxa occurring in eastern North America. These taxa were recognized as two or three species including three or four varieties by different authors. Due to high morphological similarities and complexity of morphological variation, classification and delineation of taxa in the group have been difficult and controversial. Here we employ data from RAD-seq to elucidate the genetic structure and phylogenetic relationships within the group. Using the genetic evidence, we evaluate previous classifications and delineate species. We also employ Species Distribution Modeling (SDM) to evaluate impacts of climatic changes on the ranges of the taxa and to gain insights into the relevant refugia in eastern North America. Results from Molecular Variance Analysis (AMOVA), STRUCTURE, phylogenetic analyses using Maximum likelihood, Bayesian Inference, and Splittree methods of RAD-seq data strongly support a two-clade pattern, largely separating samples of N. sylvatica from those of N. biflora-N. ursina mix. Divergence time analysis with BEAST suggests the two clades diverged in the mid Miocene. The ancestor of the present trees of N. sylvatica was suggested to be in the Pliocene and that of N. biflora-N. ursina mix in the end of the Miocene. Results from SDM predicted a smaller range in the southern part of the species present range of each clade during the Last Glacial Maximum (LGM). A northward expansion of the ranges during interglacial period and a northward shift of the ranges in the future under a model of global warming were also predicted. Our results support the recognition of two species in the complex, N. sylvatica and N. biflora, following the phylogenetic species concept. We found no genetic evidence supporting recognitions of intraspecific taxa. However, we propose subsp. ursina and subsp. biflora within N. biflora due to their distinction in habits, distributions, and habitats. Our results further support movements of trees in eastern North America in response to climatic changes. Finally, our study demonstrates that RAD-seq data and a combination of population genomics and SDM are valuable in resolving relationship and biogeographic history of closely related species that are taxonomically difficult.

Strong population bottleneck and repeated demographic expansions of Populus adenopoda (Salicaceae) in subtropical China

BACKGROUND AND AIMS

Glacial refugia and inter-/postglacial recolonization routes during the Quaternary of tree species in Europe and North America are well understood, but far less is known about those of tree species in subtropical eastern Asia. Thus, we have examined the phylogeographic history of Populus adenopoda (Salicaceae), one of the few poplars that naturally occur in this subtropical area.

METHODS

Genetic variations across the range of the species in subtropical China were surveyed using ten nuclear microsatellite loci and four chloroplast fragments (matK, trnG-psbK, psbK-psbI and ndhC-trnV). Coalescent-based analyses were used to test demographic and migration hypotheses. In addition, species distribution models (SDMs) were constructed to infer past, present and future potential distributions of the species.

KEY RESULTS

Thirteen chloroplast haplotypes were detected, and haplotype-rich populations were found in central and southern parts of the species' range. STRUCTURE analyses of nuclear microsatellite loci suggest obvious lineage admixture, especially in peripheral and northern populations. DIYABC analysis suggests that the species might have experienced two independent rounds of demographic expansions and a strong bottleneck in the late Quaternary. SDMs indicate that the species' range contracted during the Last Glacial Maximum (LGM), and contracted northward but expanded eastward during the Last Interglacial (LIG).

CONCLUSIONS

Chloroplast data and SDMs suggest that P. adenopoda might have survived in multiple glacial refugia in central and southern parts of its range during the LGM. Populations of the Yunnan-Guizhou Plateau in the southern part have high chloroplast DNA diversity, but may have contributed little to the postglacial recolonization of northern and eastern parts. The three major demographic events inferred by DIYABC coincide with the initiation of the LIG, start of the LGM and end of the LGM, respectively. The species may have experienced multiple rounds of range contraction during glacial periods and range expansion during interglacial periods. Our study corroborates the importance of combining multiple lines of evidence when reconstructing Quaternary population evolutionary histories.

Phylogeographic analysis of the East Asian goldenrod (Solidago virgaurea complex, Asteraceae) reveals hidden ecological diversification with recurrent formation of ecotypes

BACKGROUND AND AIMS

The processes and mechanisms underlying lineage diversification are major topics in evolutionary biology. Eurasian goldenrod species of the Solidago virgaurea complex show remarkable morphological and ecological diversity in the Japanese Archipelago, with ecotypic taxa well adapted to specific environments (climate, edaphic conditions and disturbance regimes). The species complex is a suitable model to investigate the evolutionary processes of actively speciating plant groups, due to its ability to evolve in relation to environmental adaptation and its historical population dynamics.

METHODS

Two chloroplast markers, 18 nuclear microsatellite markers and ddRAD-sequencing were used to infer population genetic demography of S. virgaurea complex with its related species/genera.

KEY RESULTS

Our analysis showed that populations in Japan form an evolutionary unit, which was genetically diverged from adjacent continental populations. The phylogenetic structure within the archipelago strongly corresponds to the geography, but interestingly there is no concordance between genetic structure and ecotypic boundaries; neighbouring populations of distinct ecotypes share a genetic background.

CONCLUSIONS

We propose that the traits specific to the ecotypic entities are maintained by natural selection or are very recently generated and have little effect on the genomes, making genome-wide genetic markers unsuitable for detecting ecotypic differentiation. Furthermore, some sporadically distributed taxa (found as rheophytes and alpine plants) were repeatedly generated from a more widespread taxon in geographically distant areas by means of selection. Overall, this study showed that the goldenrod complex has a high ability to evolve, enabling rapid ecological diversification over a recent timeframe.

Ecological and Genetic Divergences with Gene Flow of Two Sister Species (Leucomeris decora and Nouelia insignis) Driving by Climatic Transition in Southwest China

Understanding of the processes of divergence and speciation is a major task for biodiversity researches and may offer clearer insight into mechanisms generating biological diversity. Here, we employ an integrative approach to explore genetic and ecological differentiation of Leucomeris decora and Nouelia insignis distributed allopatrically along the two sides of the biogeographic boundary 'Tanaka Line' in Southwest China. We addressed these questions using ten low-copy nuclear genes and nine plastid DNA regions sequenced among individuals sampled from 28 populations across their geographic ranges in China. Phylogenetic, coalescent-based population genetic analyses, approximate Bayesian computation (ABC) framework and ecological niche models (ENMs) were conducted. We identified a closer phylogenetic relationship in maternal lineage of L. decora with N. insignis than that between L. decora and congeneric Leucomeris spectabilis. A deep divergence between the two species was observed and occurred at the boundary between later Pliocene and early Pleistocene. However, the evidence of significant chloroplast DNA gene flow was also detected between the marginal populations of L. decora and N. insignis. Niche models and statistical analyses showed significant ecological differentiation, and two nuclear loci among the ten nuclear genes may be under divergent selection. These integrative results imply that the role of climatic shift from Pliocene to Pleistocene may be the prominent factor for the divergence of L. decora and N. insignis, and population expansion after divergence may have given rise to chloroplast DNA introgression. The divergence was maintained by differential selection despite in the face of gene flow.

Phylogeography and conservation genetics of the rare and relict Bretschneidera sinensis (Akaniaceae)

Bretschneidera sinensis, a class-I protected wild plant in China, is a relic of the ancient Tertiary tropical flora endemic to Asia. However, little is known about its genetics and phylogeography. To elucidate the current phylogeographic patterns and infer the historical population dynamics of B. sinensis, and to make recommendations for its conservation, three non-coding regions of chloroplast DNA (trnQ-rps16, rps8-rps11, and trnT-trnL) were amplified and sequenced across 256 individuals from 23 populations of B. sinensis, spanning 10 provinces of China. We recognized 13 haplotypes, demonstrating relatively high total haplotype diversity (hT = 0.739). Almost all of the variation existed among populations (98.09%, P < 0.001), but that within populations was low (1.91%, P < 0.001). Strong genetic differentiation was detected among populations (GST = 0.855, P < 0.001) with limited estimations of seed flow (Nm = 0.09), indicating that populations were strongly isolated from one another. According to SAMOVA analysis, populations of B. sinensis in China could be divided into five geographic groups: (1) eastern Yunnan to western Guangxi; (2) Guizhou-Hunan-Hubei; (3) central Guangdong; (4) northwestern Guangdong; and (5) the Luoxiao-Nanling-Wuyi -Yangming Mountain. Network analysis showed that the most ancestral haplotypes were located in the first group, i.e., the eastern Yungui Plateau and in eastern Yunnan, which is regarded as a putative glacial refugia for B. sinensis in China. B. sinensis may have expanded its range eastward from these refugia and experienced bottleneck or founder effects in southeastern China. Populations in Liping (Guizhou Province), Longsheng (Guangxi Province), Huizhou (Guangdong Province), Chongyi (Jiangxi Province), Dong-an (Hunan Province), Pingbian (Yunnan Province) and Xinning (Hunan Province) are proposed as the priority protection units.