Phylogeographic breaks within Asian butternuts indicate the existence of a phytogeographic divide 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.

Responses of Three Pedicularis Species to Geological and Climatic Changes in the Qinling Mountains and Adjacent Areas in East Asia

The Qinling Mountains in East Asia serve as the geographical boundary between the north and south of China and are also indicative of climatic differences, resulting in rich ecological and species diversity. However, few studies have focused on the responses of plants to geological and climatic changes in the Qinling Mountains and adjacent regions. Therefore, we investigated the evolutionary origins and phylogenetic relationships of three Pedicularis species in there to provide molecular evidence for the origin and evolution of plant species. Ecological niche modeling was used to predict the geographic distributions of three Pedicularis species during the last interglacial period, the last glacial maximum period, and current and future periods, respectively. Furthermore, the distribution patterns of climate fluctuations and the niche dynamics framework were used to assess the equivalence or difference of niches among three Pedicularis species. The results revealed that the divergence of three Pedicularis species took place in the Miocene and Holocene periods, which was significantly associated with the large-scale uplifts of the Qinling Mountains and adjacent regions. In addition, the geographic distributions of three Pedicularis species have undergone a northward migration from the past to the future. The most important environmental variables affecting the geographic distributions of species were the mean diurnal range and annual mean temperature range. The niche divergence analysis suggested that the three Pedicularis species have similar ecological niches. Among them, P. giraldiana showed the highest niche breadth, covering nearly all of the climatic niche spaces of P. dissecta and P. bicolor. In summary, this study provides novel insights into the divergence and origins of three Pedicularis species and their responses to climate and geological changes in the Qinling Mountains and adjacent regions. The findings have also provided new perspectives for the conservation and management of Pedicularis species.

Intraspecific differentiation of Lindera obtusiloba as revealed by comparative plastomic and evolutionary analyses

Lindera obtusiloba Blume is the northernmost tree species in the family Lauraceae, and it is a key species in understanding the evolutionary history of this family. The species of L. obtusiloba in East Asia has diverged into the Northern and Southern populations, which are geographically separated by an arid belt. Though the morphological differences between populations have been observed and well documented, intraspecific variations at the plastomic level have not been systematically investigated to date. Here, ten chloroplast genomes of L. obtusiloba individuals were sequenced and analyzed along with three publicly available plastomes. Comparative plastomic analysis suggests that both the Northern and the Southern populations share similar overall structure, gene order, and GC content in their plastomes although the size of the plasome and the level of intraspecific variability do vary between the two populations. The Northern have relatively larger plastomes while the Southern population possesses higher intraspecific variability, which could be attributed to the complexity of the geological environments in the South. Phylogenomic analyses also support the split of the Northern and Southern clades among L. obtusiloba individuals. However, there is no obvious species boundary between var. obtusiloba and var. heterophylla in the Southern population, indicating that gene flow could still occur between these two varieties, and this could be used as a good example of reticulate evolution. It is also found that a few photosynthesis-related genes are under positive selection, which is mainly related to the geological and environmental differences between the Northern and the Southern regions. Our results provide a reference for phylogenetic analysis within species and suggest that phylogenomic analyses with a sufficient number of nuclear and chloroplast genomic target loci from widely distributed individuals could provide a deeper understanding of the population evolution of the widespread species.

Low but highly geographically structured genomic diversity of East Asian Eurasian otters and its conservation implications

Populations of Eurasian otters Lutra lutra, one of the most widely distributed apex predators in Eurasia, have been depleted mainly since the 1950s. However, a lack of information about their genomic diversity and how they are organized geographically in East Asia severely impedes our ability to monitor and conserve them in particular management units. Here, we re-sequenced and analyzed 20 otter genomes spanning continental East Asia, including a population at Kinmen, a small island off the Fujian coast, China. The otters form three genetic clusters (one of L. l. lutra in the north and two of L. l. chinensis in the south), which have diverged in the Holocene. These three clusters should be recognized as three conservation management units to monitor and manage independently. The heterozygosity of the East Asian otters is as low as that of the threatened carnivores sequenced. Historical effective population size trajectories inferred from genomic variations suggest that their low genomic diversity could be partially attributed to changes in the climate since the mid-Pleistocene and anthropogenic intervention since the Holocene. However, no evidence of genetic erosion, mutation load, or high level of inbreeding was detected in the presumably isolated Kinmen Island population. Any future in situ conservation efforts should consider this information for the conservation management units.

Impacts of Quaternary glaciation, geological history and geography on animal species history in continental East Asia: A phylogeographic review

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

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.

Phylogeny and Historical Biogeography of the East Asian Clematis Group, Sect. Tubulosae, Inferred from Phylogenomic Data

The evolutionary history of Clematis section Tubulosae, an East Asian endemic lineage, has not been comprehensively studied. In this study, we reconstruct the phylogeny of this section with a complete sampling using a phylogenomic approach. The genome skimming method was applied to obtain the complete plastome sequence, the nuclear ribosomal DNA (nrDNA), and the nuclear SNPs data for phylogenetic reconstruction. Using a Bayesian molecular clock approach and ancestral range reconstruction, we reconstruct biogeographical history and discuss the biotic and abiotic factors that may have shaped the distribution patterns of the section. Both nuclear datasets better resolved the phylogeny of the sect. Tubulosae than the plastome sequence. Sect. Tubulosae was resolved as a monophyletic group sister to a clade mainly containing species from the sect. Clematis and sect. Aspidanthera. Within sect. Tubulosae, two major clades were resolved by both nuclear datasets. Two continental taxa, C. heracleifolia and C. tubulosa var. ichangensis, formed one clade. One continental taxon, C. tubulosa, and all the other species from Taiwan island, the Korean peninsula, and the Japanese archipelago formed the other clade. Molecular dating results showed that sect. Tubulosae diverged from its sister clade in the Pliocene, and all the current species diversified during the Pleistocene. Our biogeographical reconstruction suggested that sect. Tubulosae evolved and began species diversification, most likely in mainland China, then dispersed to the Korean peninsula, and then expanded its range through the Japanese archipelago to Taiwan island. Island species diversity may arise through allopatric speciation by vicariance events following the range fragmentation triggered by the climatic oscillation and sea level change during the Pleistocene epoch. Our results highlight the importance of climatic oscillation during the Pleistocene to the spatial-temporal diversification patterns of the sect. Tubulosae.

Genomic signals of local adaptation and hybridization in Asian white birch

Disentangling the numerous processes that affect patterns of genome-wide diversity in widespread tree species has important implications for taxonomy, conservation, and forestry. Here, we investigate the population genomic structure of Asian white birch (Betula platyphylla) in China and seek to explain it in terms of hybridization, demography and adaptation. We generate whole genome sequence data from 83 individuals across the species range in China. Combining this with an existing data set for 79 European and Russian white birches, we show a clear distinction between B. pendula and B. platyphylla, which have sometimes been lumped taxonomically. Genomic diversity of B. platyphylla in north-western China and Central Russia is affected greatly by hybridization with B. pendula. Excluding these hybridized populations, B. platyphylla in China has a linear distribution from north-eastern to south-western China, along the edge of the inland mountainous region. Within this distribution, three genetic clusters are found, which we model as long diverged with subsequent episodes of gene flow. Patterns of covariation between allele frequencies and environmental variables in B. platyphylla suggest the role of natural selection in the distribution of diversity at 7609 SNPs of which 3767 were significantly differentiated among the genetic clusters. The putative adaptive SNPs are distributed throughout the genome and span 1633 genic regions. Of these genic regions, 87 were previously identified as candidates for selective sweeps in Eurasian B. pendula. We use the 7609 environmentally associated SNPs to estimate the risk of nonadaptedness for each sequenced B. platyphylla individual under a scenario of future climate change, highlighting areas where populations may be under future threat from rising temperatures.

Monsoon intensification in East Asia triggered the evolution of its flora

INTRODUCTION

East Asia (EA), which falls within the region of the Asian monsoon that is composed of the East Asia monsoon (EAM) and the Indian monsoon (IM), is known for its high species diversity and endemism. This has been attributed to extreme physiographical heterogeneity in conjunction with climate and sea-level changes during the Pleistocene, this hypothesis has been widely proven by phylogeographic studies. Recently, dated phylogenies have indicated that the origins (stem age) of the flora occurred after the Oligocene-Miocene boundary and are related to the establishment of the EAM.

METHODS

Hence, this study further examined whether the strengthening of the monsoons triggered floral evolution via a meta-analysis of the tempo-spatial pattern of evolutionary radiation dates (crown ages) of 101 endemic seed plant genera.

RESULTS

Taxonomic diversification began during the late Eocene, whereas the accumulated number of diversifications did not significantly accelerate until the late Miocene. The distribution of the weighted mean and the average divergence times in the EAM, IM, or transitional regions all fall within the mid-late Miocene. Fossils of the Tertiary relict genera are mostly and widely distributed outside EA and only half of the earliest fossils in the EA region are not older than Miocene, while their divergence times are mostly after the late Miocene. The pattern of divergence time of monotypic and polytypic taxa suggest the climatic changes after the late Pliocene exert more influence on monotypic taxa.

DISCUSSION

The two key stages of floral evolution coincide with the intensifications of the EAM and IM, especially the summer monsoon which brings a humid climate. An integrated review of previous studies concerning flora, genus, and species levels further supports our suggestion that monsoon intensification in EA triggered the evolution of its flora.

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.

Population-genomic analyses reveal bottlenecks and asymmetric introgression from Persian into iron walnut during domestication

Background

Persian walnut, Juglans regia, occurs naturally from Greece to western China, while its closest relative, the iron walnut, Juglans sigillata, is endemic in southwest China; both species are cultivated for their nuts and wood. Here, we infer their demographic histories and the time and direction of possible hybridization and introgression between them.

Results

We use whole-genome resequencing data, different population-genetic approaches (PSMC and GONE), and isolation-with-migration models (IMa3) on individuals from Europe, Iran, Kazakhstan, Pakistan, and China. IMa3 analyses indicate that the two species diverged from each other by 0.85 million years ago, with unidirectional gene flow from eastern J. regia and its ancestor into J. sigillata, including the shell-thickness gene. Within J. regia, a western group, located from Europe to Iran, and an eastern group with individuals from northern China, experienced dramatically declining population sizes about 80 generations ago (roughly 2400 to 4000 years), followed by an expansion at about 40 generations, while J. sigillata had a constant population size from about 100 to 20 generations ago, followed by a rapid decline.

Conclusions

Both J. regia and J. sigillata appear to have suffered sudden population declines during their domestication, suggesting that the bottleneck scenario of plant domestication may well apply in at least some perennial crop species. Introgression from introduced J. regia appears to have played a role in the domestication of J. sigillata.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02720-z.

The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis

BACKGROUND

Manchurian walnut (Juglans mandshurica Maxim.) is a tree with multiple industrial uses and medicinal properties in the Juglandaceae family (walnuts and hickories). J. mandshurica produces juglone, which is a toxic allelopathic agent and has potential utilization value. Furthermore, the seed of J. mandshurica is rich in various unsaturated fatty acids and has high nutritive value.

FINDINGS

Here, we present a high-quality chromosome-scale reference genome assembly and annotation for J. mandshurica (n = 16) with a contig N50 of 21.4 Mb by combining PacBio high-fidelity reads with high-throughput chromosome conformation capture data. The assembled genome has an estimated sequence size of 548.7 Mb and consists of 657 contigs, 623 scaffolds, and 40,453 protein-coding genes. In total, 60.99% of the assembled genome consists of repetitive sequences. Sixteen super-scaffolds corresponding to the 16 chromosomes were assembled, with a scaffold N50 length of 33.7 Mb and a BUSCO complete gene percentage of 98.3%. J. mandshurica displays a close sequence relationship with Juglans cathayensis, with a divergence time of 13.8 million years ago. Combining the high-quality genome, transcriptome, and metabolomics data, we constructed a gene-to-metabolite network and identified 566 core and conserved differentially expressed genes, which may be involved in juglone biosynthesis. Five CYP450 genes were found that may contribute to juglone accumulation. NAC, bZip, NF-YA, and NF-YC are positively correlated with the juglone content. Some candidate regulators (e.g., FUS3, ABI3, LEC2, and WRI1 transcription factors) involved in the regulation of lipid biosynthesis were also identified.

CONCLUSIONS

Our genomic data provide new insights into the evolution of the walnut genome and create a new platform for accelerating molecular breeding and improving the comprehensive utilization of these economically important tree species.

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.

Phylogeography and Population Genetics of Rosa chinensis var. spontanea and R. lucidissima Complex, the Important Ancestor of Modern Roses

Rosa chinensis var. spontanea and R. lucidissima complex are the morphologically very similar key ancestors of modern roses with high importance in rose research and breeding. Although widely distributed in subtropical central and southwestern China, these two taxa are highly endangered. We sampled a total of 221 specimens and 330 DNA samples from 25 populations across the two taxa's whole range. Leaf morphological traits were compared. Two chloroplast DNA intergenic spacers (trnG-trnS, petL-psbE) and ITS were used for population genetics and phylogenetic study to delimit the boundary between the two taxa, assess the genetic variation, uncover the possible evolutionary mechanism responsible for the differentiation within the complex, and make the conservation recommendations. The complex exhibited high levels of genetic variation (h _TcpDNA = 0.768, h _TITS = 0.726) and high population differentiation even over small geographic distance. We suggest R. chinensis var. spontanea and R. lucidissma be treated as independent taxa, and the northern populations around and within the Sichuan Basin being R. chinensis var. spontanea, having broader leaflets and paler full-blooming flowers, while those in the middle and southern Yunnan-Guizhou Plateau and the adjacent regions being R. lucidissma, having narrower leaflets and darker full-blooming flowers. Transitional areas between the southeastern Sichuan Basin and northeastern Guizhou are the contact or the hybridization zone of the two taxa. Ancestral haplotypes of the complex (R. lucidissma) evolved at about 1.21-0.86 Mya in southeastern Yunnan-Guizhou Plateau and its adjacent regions and survived there during the Quaternary Oscillation. Ancestral haplotypes of R. chinensis var. spontanea deviated from R. lucidissma at about 0.022-0.031 Mya at the transitional areas (Daloushan and Wulingshan Mountains) between the northeastern edge of Yunnan-Guizhou Plaeteau and the southeastern border of Sichuan Basin, where they survived the LGM. The evolution of the complex included spatial isolation and inter-species hybridization. The complex's endangered status might be the result of over-exploitation for its ornamental and medical value, or due to reforestation of some originally open habitats. We provide specific recommendations for the two taxa's in situ and ex situ conservation.

Dead-End Hybridization in Walnut Trees Revealed by Large-Scale Genomic Sequence Data

Although hybridization plays a large role in speciation, some unknown fraction of hybrid individuals never reproduces, instead remaining as genetic dead-ends. We investigated a morphologically distinct and culturally important Chinese walnut, Juglans hopeiensis, suspected to have arisen from hybridization of Persian walnut (J. regia) with Asian butternuts (J. cathayensis, J. mandshurica, and hybrids between J. cathayensis and J. mandshurica). Based on 151 whole-genome sequences of the relevant taxa, we discovered that all J. hopeiensis individuals are first-generation hybrids, with the time for the onset of gene flow estimated as 370,000 years, implying both strong postzygotic barriers and the presence of J. regia in China by that time. Six inversion regions enriched for genes associated with pollen germination and pollen tube growth may be involved in the postzygotic barriers that prevent sexual reproduction in the hybrids. Despite its long-recurrent origination and distinct traits, J. hopeiensis does not appear on the way to speciation.

Spatial phylogenetics of the Chinese angiosperm flora provides insights into endemism and conservation

The flora of China is well known for its high diversity and endemism. Identifying centers of endemism and designating conservation priorities are essential goals for biodiversity studies. However, there is no comprehensive study from a rigorous phylogenetic perspective to understand patterns of diversity and endemism and to guide biodiversity conservation in China. We conducted a spatial phylogenetic analysis of the Chinese angiosperm flora at the generic level to identify centers of neo- and paleo-endemism. Our results indicate that: (i) the majority of grid cells in China with significantly high phylogenetic endemism (PE) were located in the mountainous regions; (ii) four of the nine centers of endemism recognized, located in northern and western China, were recognized for the first time; (iii) arid and semiarid regions in Northwest China were commonly linked to significant PE, consistent with other spatial phylogenetic studies worldwide; and (iv) six high-priority conservation gaps were detected by overlaying the boundaries of China's nature reserves on all significant PE cells. Overall, we conclude that the mountains of southern and northern China contain both paleo-endemics (ancient relictual lineages) and neo-endemics (recently diverged lineages). The areas we highlight as conservation priorities are important for broad-scale planning, especially in the context of evolutionary history preservation.

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.

Demographic history and local adaptation of Myripnois dioica (Asteraceae) provide insight on plant evolution in northern China flora

The flora of northern China forms the main part of the Sino-Japanese floristic region and is located in a south-north vegetative transect in East Asia. Phylogeographic studies have demonstrated that an arid belt in this region has promoted divergence of plants in East Asia. However, little is known about how plants that are restricted to the arid belt of flora in northern China respond to climatic oscillation and environmental change. Here, we used genomic-level data of Myripnois dioica across its distribution as a representative of northern China flora to reconstruct plant demographic history, examine local adaptation related to environmental disequilibrium, and investigate the factors related to effective population size change. Our results indicate M. dioica originated from the northern area and expanded to the southern area, with the Taihang Mountains serving as a physical barrier promoting population divergence. Genome-wide evidence found strong correlation between genomic variation and environmental factors, specifically signatures associated with local adaptation to drought stress in heterogeneous environments. Multiple linear regression analyses revealed joint effects of population age, mean temperature of coldest quarter, and precipitation of wettest month on effective population size (Ne). Our current study uses M. dioica as a case for providing new insights into the evolutionary history and local adaptation of northern China flora and provides qualitative strategies for plant conservation.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Improved de novo chromosome-level genome assembly of the vulnerable walnut tree Juglans mandshurica reveals gene family evolution and possible genome basis of resistance to lesion nematode

Manchurian walnut (Juglans mandshurica Maxim.) is a synonym of J. cathayensis, a diploid, vulnerable, temperate deciduous tree valued for its wood and nut. It is also valued as a rootstock for Juglans regia because of its reported tolerance of lesion nematode. Reference genomes are available for several Juglans species, our goal was to produce a de novo, chromosome-level assembly of the J. mandshurica genome. Here, we reported an improved assembly of J. mandshurica with a contig N50 size of 6.49 Mb and a scaffold N50 size of 36.1 Mb. The total genome size was 548 Mb encoding 29,032 protein coding genes which were annotated. The collinearity analysis showed that J. mandshurica and J. regia originated from a common ancestor, with both species undergoing two WGD events. A genomic comparison showed that J. mandshurica was missing 1657 genes found in J. regia, and J. mandshurica includes 2827 genes not found in of the J. regia genome. The J. mandshurica contained 1440 unique paralogues that were highly enriched for flavonoid biosynthesis, phenylpropanoid biosynthesis, and plant-pathogen interaction. Four gene families related to disease resistance notable contraction (rapidly evolving; LEA, WAK, PPR, and PR) in J. mandshurica compared to eight species. JmaPR10 and JmaPR8 contained three orthologous gene pairs with J. regia that were highly expressed in root bark. JmaPR10 is a strong candidate gene for lesion nematodes resistance in J. mandshurica. The J. mandshurica genome should be a useful resource for study of the evolution, breeding, and genetic variation in walnuts (Juglans).

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.

Subspecies divergence and pronounced phylogenetic incongruence in the East-Asia-endemic shrub Magnolia sieboldii

BACKGROUND AND AIMS

The biogeographic patterns of the East-Asia-endemic shrub Magnolia sieboldii, in which the range of the subsp. sieboldii is interposed with the disjunct distribution of subsp. japonica, implies a complex evolutionary history, involving rapid speciation and hybridization. Here, we aim to reveal the evolutionary and phylogeographic histories of the species with a particular focus on the time of subspecies divergence, the hypothesis of secondary hybridization and the Pleistocene survival of each subspecies, using a combination of genetic analyses and ecological niche modelling.

METHODS

Genetic variation, genetic structures and phylogenetic relationships were elucidated based on nuclear low-copy genes, chloroplast DNA, and nuclear simple sequence repeats (SSRs). A scenario selection analysis and divergence time estimation were performed using coalescent simulation in DIYABC and *BEAST. Ecological niche modelling and a test of niche differentiation were performed using Maxent and ENMTools.

KEY RESULTS

All marker types showed deep, but pronouncedly incongruent, west-east genetic divergences, with the subspecies being delineated only by the nuclear low-copy genes. Phylogenetic tree topologies suggested that ancient hybridization and introgression were likely to have occurred; however, this scenario did not receive significant support in the DIYABC analysis. The subspecies differentiated their niches, but both showed a dependence on high humidity and were predicted to have persisted during the last glacial cycle by maintaining a stable latitudinal distribution via migration to lower altitudes.

CONCLUSIONS

We found a deep genetic divergence and a pronounced phylogenetic incongruence among the two subspecies of M. sieboldii, which may have been driven by major paleogeographic and paleoclimatic events that have occurred since the Neogene in East Asia, including global cooling, climate oscillations and the formation of land bridges. Both subspecies were, however, considered to persist in situ in stable climatic conditions during the late Pleistocene.

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.

The East Asian Winter Monsoon Acts as a Major Selective Factor in the Intraspecific Differentiation of Drought-Tolerant Nitraria tangutorum in Northwest China

The influence of Quaternary climate fluctuation on the geographical structure and genetic diversity of species distributed in the regions of the Qinghai–Tibet Plateau (QTP) has been well established. However, the underlying role of the East Asian monsoon system (EAMS) in shaping the genetic structure of the population and the demography of plants located in the arid northwest of China has not been explored. In the present study, Nitraria tangutorum, a drought-tolerant desert shrub that is distributed in the EAMS zone and has substantial ecological and economic value, was profiled to better understand the influence of EAMS evolution on its biogeographical patterns and demographic history. Thus, the phylogeographical structure and historical dynamics of this plant species were elucidated using its five chloroplast DNA (cpDNA) fragments. Hierarchical structure analysis revealed three distinct, divergent lineages: West, East-A, and East-B. The molecular dating was carried out using a Bayesian approach to estimate the time of intraspecies divergence. Notably, the eastern region, which included East-A and East-B lineages, was revealed to be the original center of distribution and was characterized by a high level of genetic diversity, with the intraspecific divergence time dated to be around 2.53 million years ago (Ma). These findings, combined with the data obtained by ecological niche modeling analysis, indicated that the East lineages have undergone population expansion and differentiation, which were closely correlated with the development of the EAMS, especially the East Asian winter monsoon (EAWM). The West lineage appears to have originated from the migration of N. tangutorum across the Hexi corridor at around 1.85 Ma, and subsequent colonization of the western region. These results suggest that the EAWM accelerated the population expansion of N. tangutorum and subsequent intraspecific differentiation. These findings collectively provide new information on the impact of the evolution of the EAMS on intraspecific diversification and population demography of drought-tolerant plant species in northwest China.

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.

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.

Population genetic structure and adaptive differentiation of iron walnut Juglans regia subsp. sigillata in southwestern China

Southwestern (SW) China is an area of active tectonism and erosion, yielding a dynamic, deeply eroded landscape that influences the genetic structure of the resident populations of plants and animals. Iron walnut (Juglans regia subsp. sigillata) is a deciduous tree species endemic to this region of China and cultivated there for its edible nuts. We sampled 36 iron walnut populations from locations throughout the species' range in SW China and genotyped a total of 765 individuals at five chloroplast DNA regions and 22 nuclear microsatellite loci. Species distribution models were produced to predict the evolution and historical biogeography of iron walnut and to estimate the impacts of climate oscillations and orographic environments on the species' demography. Our results indicated that J. regia subsp. sigillata had relatively low genetic diversity, high interpopulation genetic differentiation, and asymmetric interpopulation gene flow. Based on DIYABC analysis, we identified two lineages of J. sigillata in southwestern China. The lineages (subpopulations) diverge during the last glacial period (~1.34 Ma). Southwestern China was a glacial refuge during the last glacial period, but increasingly colder and arid climates might have fostered the fragmentation of J. regia subsp. sigillata within this refugium. Finally, we found that recent habitat fragmentation has led to a reduction in population connectivity and increased genetic differentiation by genetic drift in isolated populations. Our results support a conclusion that geological and climatic factors since the Miocene triggered the differentiation, evolutionary origin, and range shifts of J. sigillata in the studied region.

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.

Molecular Phylogeography Analysis Reveals Population Dynamics and Genetic Divergence of a Widespread Tree Pterocarya stenoptera in China

The geological events, past climatic fluctuations, and river systems played key roles in the spatial distribution, population dynamics, and genetic differentiation of species. In this work, we selected Pterocarya stenoptera, a widespread tree species in China, to test the roles of these factors. Four noncoding spacers, eight microsatellite (simple sequence repeat) markers, and species distribution modeling were used to examine the phylogeographical pattern of P. stenoptera. Based on chloroplast DNA data, populations of P. stenoptera were clearly clustered into three groups. The divergence time of these groups fell into the stage of the Qinghai–Tibet Movement, 1.7–2.6 Ma. For simple sequence repeat data, only one western marginal population YNYB could be separated from other populations, whereas other populations were mixed together. Our results indicated that the environmental heterogeneity resulting from the Qinghai–Tibet movement might be response for this genetic divergence. The climatic fluctuations in the Pleistocene did not cause the substantial range shift of P. stenoptera, while the fluctuations affected its population size. Moreover, we also confirmed the river systems did not act as channels or barrier of dispersal for P. stenoptera.

Gene Introgression among Closely Related Species in Sympatric Populations: A Case Study of Three Walnut (Juglans) Species

Gene introgression usually results from natural hybridization occurring among closely related species in sympatric populations. In this study, we discussed two rare and frequent gene flow phenomena between three species of Juglans plants and analyzed the possible causes for the difference. We collected 656 individuals from 40 populations of Persian walnut (Juglans regia L.), Chinese walnut (J. cathayensis Dode), and Iron walnut (J. sigillata Dode) that were genotyped at 17 expressed sequence tag simple sequence repeat (EST-SSR) loci to analyze the introgressions between J. regia and J. cathayensis, and J. regia and J. sigillata. Our study compared the spatial patterns of expected heterozygosity (HE), allelic richness (Rs), and private allele richness (PAR) so as to vividly infer the biogeographic history of related species of Juglans in the two regions. The results of the PCoA, UPGMA, and STRUCTURE analyses showed that all J. regia and J. sigillata populations clustered into one group, and the J. cathayensis populations clustered into the other group. The results of the historical gene flow analysis indicated that J. regia and J. sigillata have no genetic barriers, and the directional gene flow is mainly from J. regia to J. sigillata. For the three species of Juglans, all the above results indicated that gene flow was common among the same group of Juglans, and only rare and low-level gene flow appeared in distinct groups. Therefore, our study revealed multiple phenomena of gene flow and introgression among closely related species in sympatric populations, thereby providing a theoretical basis for the genetic evolution of the genus Juglans.

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.

Landscape genetics reveals that adaptive genetic divergence in Pinus bungeana (Pinaceae) is driven by environmental variables relating to ecological habitats

Background

Understanding the genetic basis of local adaptation has long been the concern of biologists. Identifying these adaptive genetic variabilities is crucial not only to improve our knowledge of the genetic mechanism of local adaptation but also to explore the adaptation potential of species.

Results

Using 10 natural populations and 12 start codon targeted (SCoT) markers, a total of 430 unambiguous loci were yielded. The Bayesian analysis of population structure clearly demonstrated that the 10 populations of P. bungeana could be subdivided into three groups. Redundancy analysis showed that this genetic divergence was caused by divergence selection from environmental variables related to the ecological habitats of “avoidance of flooding” and “avoidance of high temperature and humidity.” LFMM results indicated that Bio1, Bio5, Bio8, Bio12, Bio14, and Bio16, which are related to the ecological habitat of P. bungeana, were correlated with the highest numbers of environment-associated loci (EAL).

Conclusions

The results of EAL characterization in P. bungeana clearly supported the hypothesis that environmental variations related to the ecological habitat of species are the key drivers of species adaptive divergence. Moreover, a method to calculate the species landscape adaptation index and quantify the adaptation potential of species was proposed and verified using ecological niche modeling. This model could estimate climatically suitable areas of species spatial distribution. Taking the results together, this study improves the current understanding on the genetic basis of local adaptation.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1489-x) contains supplementary material, which is available to authorized users.

Integrating Phylogeographic Analysis and Geospatial Methods to Infer Historical Dispersal Routes and Glacial Refugia of Liriodendron chinense

Liriodendron chinense (Hemsl.), a Tertiary relic tree, is mainly distributed in subtropical China. The causes of the geographical distribution pattern of this species are poorly understood. In this study, we inferred historical dispersal routes and glacial refugia of this species by combining genetic data (chloroplast DNA (cpDNA), nuclear ribosomal DNA (nrDNA), and nuclear DNA (nDNA)) and geospatial data (climate and geology) with the methods of landscape genetics. Additionally, based on sequence variation at multiple loci, we employed GenGIS and Barrier software to analyze L. chinense population genetic structure. Dispersal corridors and historical gene flow between the eastern and western populations were detected, and they were located in mountainous regions. Based on species distribution model (SDMs), the distribution patterns in paleoclimatic periods were consistent with the current pattern, suggesting the presence of multiple refuges in multiple mountainous regions in China. The genetic structure analysis clustered most eastern populations into a clade separated from the western populations. Additionally, a genetic barrier was detected between the eastern and western populations. The dispersal corridors and historical gene flow detected here suggested that the mountains acted as a bridge, facilitating gene flow between the eastern and western populations. Due to Quaternary climatic fluctuations, the habitats and dispersal corridors were frequently inhabited by warm-temperate evergreen forests, which may have fragmented L. chinense habitats and exacerbated the differentiation of eastern and western populations. Ultimately, populations retreated to multiple isolated mountainous refugia, shaping the current geographical distribution pattern. These dispersal corridors and montane refugia suggested that the mountains in subtropical China play a crucial role in the conservation of genetic resources and migration of subspecies or related species in this region.

Phylogenomics Reveals an Ancient Hybrid Origin of the Persian Walnut

Persian walnut (Juglans regia) is cultivated worldwide for its high-quality wood and nuts, but its origin has remained mysterious because in phylogenies it occupies an unresolved position between American black walnuts and Asian butternuts. Equally unclear is the origin of the only American butternut, J. cinerea. We resequenced the whole genome of 80 individuals from 19 of the 22 species of Juglans and assembled the genome of its relatives Pterocarya stenoptera and Platycarya strobilacea. Using phylogenetic-network analysis of single-copy nuclear genes, genome-wide site pattern probabilities, and Approximate Bayesian Computation, we discovered that J. regia (and its landrace J. sigillata) arose as a hybrid between the American and the Asian lineages and that J. cinerea resulted from massive introgression from an immigrating Asian butternut into the genome of an American black walnut. Approximate Bayesian Computation modeling placed the hybrid origin in the late Pliocene, ∼3.45 My, with both parental lineages since having gone extinct in Europe.

Phylogeographic Analyses of the East Asian Endemic Genus Prinsepia and the Role of the East Asian Monsoon System in Shaping a North-South Divergence Pattern in China

Prinsepia Royle (Rosaceae) is a genus native to China and the Himalayan region. In order to explain its current fragmented distribution pattern and to compare the impact of relatively recent climate changes on the genetic structure of Prinsepia species in different regions of China, a total of 66 populations and 617 individuals of four species of Prinsepia were genotyped, using three cpDNA markers. Meanwhile, phylogenetic reconstructions and divergence dating were conducted using the cpDNA haplotypes dataset and the nuclear ribosomal internal transcribed spacer (ITS) dataset, respectively. Ecological niche modeling (ENM) was performed to predict the potential distribution of each species of Prinsepia at present and during the Last Glacial Maximum. Both ITS and cpDNA gene trees support a north-south divergence of Prinsepia species in China. The divergence time of the northern and southern Clades occurred around the late Oligocene epoch. Combining the present distribution of Prinsepia species and their habitats, we inferred that the transition to a monsoon climate system in East Asia during the late Oligocene epoch, created a humid forest vegetation zone from central to East China, which potentially gave rise to the north-south divergence of Prinsepia species. Both regional climates and allopatric divergence may have played an important role in the speciation of P. sinensis and P. uniflora. P. sinensis had the lowest genetic diversity and a putative northward post-glacial colonization. The distribution range of P. uniflora was also extremely sensitive to interglacial-glacial cycles. P. utilis from southwestern China preserved more haplotypes than P. sinensis and P. uniflora due to its multiple and isolated refugia.

Phylogeography and conservation genetics of the endangered Tugarinovia mongolica (Asteraceae) from Inner Mongolia, Northwest China

Tugarinovia (Family Asteraceae) is a monotypic genus. It’s sole species, Tugarinovia mongolica Iljin, is distributed in the northern part of Inner Mongolia, with one additional variety, Tugarinovia mongolica var ovatifolia, which is distributed in the southern part of Inner Mongolia. The species has a limited geographical range and declining populations. To understand the phylogeographic structure of T. mongolica, we sequenced two chloroplast DNA regions (psbA-trnH and psbK-psbI) from 219 individuals of 16 populations, and investigated the genetic variation and phylogeographic patterns of T. mongolica. The results identified a total of 17 (H1-H17) chloroplast haplotypes. There were no haplotypes shared between the northern (T. mongolica) and southern groups (T. mongolica var. ovatifolia), and they formed two distinct lineages. The regional split was also supported by AMOVA and BEAST analyses. AMOVA showed the main variation that occurred between the two geographic groups. The time of divergence of the two groups can be dated to the early Pleistocene epoch, when climate fluctuations most likely resulted in the allopatric divergence of T. mongolica. The formation of the desert blocked genetic flow and enhanced the divergence of the northern and southern groups. Our results indicate that the genetic differences between T. mongolica and T. mongolica var. ovatifolia are consistent with previously proposed morphological differences. We speculate that the dry, cold climate and the expansion of the desert during the Quaternary resulted in the currently observed distribution of extant populations of T. mongolica. In the northern group, the populations Chuanjinsumu, Wuliji and Yingen displayed the highest genetic diversity and should be given priority protection. The southern group showed a higher genetic drift (F_ST = 1, G_ST = 1), and the inbreeding load (H_S = 0) required protection for each population. Our results propose that the protection of T. mongolica should be implemented through in situ and ex situ conservation practices to increase the effective population size and genetic diversity.

A Phytogeographic Divide Along the 500 mm Isohyet in the Qinghai-Tibet Plateau: Insights From the Phylogeographic Evidence of Chinese Alliums (Amaryllidaceae)

The Qinghai-Tibet Plateau (QTP) has been biogeographically divided into the eastern monsoonal and the western continental climatic zones along the 500 mm isohyet. However, this biogeographic hypothesis has been rarely tested using a phylogeographic approach. The members of the genus Allium subgenus Cyathophora coincidentally distribute across this biogeographical divide. Intriguingly, Allium fasciculatum of subgenus Amerallium co-occurs in the distribution range of subgenus Cyathophora. To illuminate the role of this biogeographic divide on the genetic divergence, we genotyped 466 individuals of 52 populations of subgenus Cyathophora and 110 individuals of 19 populations of A. fasciculatum using three chloroplast DNA fragments, whole nrITS and nine nuclear microsatellite loci, supplemented with the present environmental space and paleo-distribution modeling. Our phylogeographical evidence recovered the concordant east-west genetic breaks both for subgenus Cyathophora and A. fasciculatum along the 500 mm isohyet. The divergence time estimations and environmental niche differentiations suggested this east-west genetic breaks could have been triggered by the climatic-induced vicariance during the early Pleistocene. Noticeably, this split within subgenus Cyathophora could have been deepened by the morphological vicariance from the eastern umbel to the western spicate, while that within A. fasciculatum could have been obscured due to the pollen flows from the east to west caused by the postglacial expansion. The genetic structures and ecological niche modelings (ENMs) recovered the distinct responses to the Quaternary climatic oscillations for species constricted to different climatic zones, further highlighting the profound effect of the climatic differences and tectonic uplifts on the genetic diversification. Overall, our findings offer strong evidence for the existence of a biogeographic divide between the eastern monsoonal and the west continental climatic zones of the QTP nearly along the 500 mm isohyet.

Multiple refugia from penultimate glaciations in East Asia demonstrated by phylogeography and ecological modelling of an insect pest

BACKGROUND

Refugial populations in Quaternary glaciations are critical to understanding the evolutionary history and climatic interactions of many extant species. Compared with the well-studied areas of Europe and Northern America, refugia of species in eastern Asia remain largely unknown. Here, we investigated the phylogeographic history of a globally important insect pest, the oriental fruit moth Grapholita molesta, in its native range of China.

RESULTS

Genetic structure analyses unveiled three distinct groups and a set of populations with admixture. Approximate Bayesian Computation (ABC) analyses support range expansion of this moth from southwest groups of Yunnan and Sichuan to northern and eastern China. A set of admixed populations was found around these two ancestral groups. This pattern of genetic structure points to two refugia located in the Yunnan region and Sichuan Basin. The split of the two refugia was dated to 329.2 thousand years ago in the penultimate glacial period. One of the lineages was exclusively found around the Sichuan Basin, indicating the formation of endemic populations in this refugium. Ecological niche model analysis suggested a shrinking distribution from the LIG period to the MID period in the Sichuan lineage but a wide and stable distribution in the other lineage.

CONCLUSIONS

Our results for the first time suggest that Yunnan and Sichuan jointly served as two large-scale refugia in eastern Asia in Quaternary glaciations, helping to maintain genetic diversity overall.

The Phytogeographic History of Common Walnut in China

Common walnut (Juglans regia L.) is an economically important hardwood tree species cultivated worldwide for its high quality wood and edible nuts. It is generally accepted that after the last glaciation J. regia survived and grew in almost completely isolated stands in Asia, and that ancient humans dispersed walnuts across Asia and into new habitats via trade and cultural expansion. The history of common walnut in China is a matter of debate, however. We estimated the genetic diversity and spatial genetic structure of 31 walnut populations sampled across its Chinese range using 22 microsatellite markers (13 neutral and 9 non-neutral). Using historical data and population genetic analysis, including approximate Bayesian analysis (ABC), we reconstructed the demographic history of J. regia in China. The genetic data indicated the likely presence of J. regia in glacial refugia in the Xinjiang province (Northwest China), Northeastern China (Beijing, Shandong, and Changbai Mountains), Central China (Qinling and Baishan Mountains and Xi'an), and Southwestern China (Tibet, Yunnan, Guizhou, and Sichuan provinces). Based on DIY-ABC analysis, we identified three ancient lineages of J. regia in China. Two lineages (subpopulation A and subpopulation B+C) diverged about 2.79 Mya, while Southwestern China, and Qinling and Baishan Mountains lineages diverged during the Quaternary glaciations (about 1.13 Mya). Remnants of these once-distinct genetic clusters of J. regia may warrant ecological management if they are to be retained as in situ resources. A population size expansion in Northeastern China was detected in the last five centuries. The present distribution of walnut in China resulted from the combined effects of expansion/contraction from multiple refugia after the Last Glacial Maximum and later human exploitation.

Genomic Variation Among and Within Six Juglans Species

Genomic analysis in Juglans (walnuts) is expected to transform the breeding and agricultural production of both nuts and lumber. To that end, we report here the determination of reference sequences for six additional relatives of Juglans regia: Juglans sigillata (also from section Dioscaryon), Juglans nigra, Juglans microcarpa, Juglans hindsii (from section Rhysocaryon), Juglans cathayensis (from section Cardiocaryon), and the closely related Pterocarya stenoptera. While these are ‘draft’ genomes, ranging in size between 640Mbp and 990Mbp, their contiguities and accuracies can support powerful annotations of genomic variation that are often the foundation of new avenues of research and breeding. We annotated nucleotide divergence and synteny by creating complete pairwise alignments of each reference genome to the remaining six. In addition, we have re-sequenced a sample of accessions from four Juglans species (including regia). The variation discovered in these surveys comprises a critical resource for experimentation and breeding, as well as a solid complementary annotation. To demonstrate the potential of these resources the structural and sequence variation in and around the polyphenol oxidase loci, PPO1 and PPO2 were investigated. As reported for other seed crops variation in this gene is implicated in the domestication of walnuts. The apparently Juglandaceae specific PPO1 duplicate shows accelerated divergence and an excess of amino acid replacement on the lineage leading to accessions of the domesticated nut crop species, Juglans regia and sigillata.

Long-distance dispersal or postglacial contraction? Insights into disjunction between Himalaya-Hengduan Mountains and Taiwan in a cold-adapted herbaceous genus, Triplostegia

Current disjunct patterns can result from long-distance dispersal or postglacial contraction. We herein investigate the evolutionary history of Triplostegia to elucidate the disjunction between the Himalaya-Hengduan Mountain region (HHM) and Taiwan (TW). Genetic structure of Triplostegia was investigated for 48 populations using sequences from five chloroplast loci and the ribosomal nuclear internal transcribed spacer. Divergence time estimation, ancestral area reconstruction, and species distribution modeling (SDM) were employed to examine the biogeographic history of Triplostegia. Substantial genetic differentiation among populations from southwestern China (SW), Central China (CC), and TW was detected. Triplostegia was inferred to have originated in SW, and diversification began during the late Miocene; CC was colonized in the mid-Pliocene, and TW was finally colonized in the early Pleistocene. SDM suggested an expansion of climatically suitable areas during the Last Glacial Maximum and range contraction during the Last interglacial in Triplostegia. Disjunction between HHM and TW in Triplostegia is most likely the consequence of topographic isolation and postglacial contraction. The potential climatic suitability areas for Triplostegia by 2070s (2061-2080) are predicted to slightly shrink and move northward. With continued global warming and human-induced deforestation, extinction risk may increase for the cold-adapted species, and appropriate strategies should be employed for ecosystem conservation.

Ice age unfrozen: severe effect of the last interglacial, not glacial, climate change on East Asian avifauna

BACKGROUND

The glacial-interglacial cycles in the Pleistocene caused repeated range expansion and contraction of species in several regions in the world. However, it remains uncertain whether such climate oscillations had similar impact on East Asian biota, despite its widely recognized importance in global biodiversity. Here we use both molecular and ecological niche profiles on 11 East Asian avian species with various elevational ranges to reveal their response to the late Pleistocene climate changes.

RESULTS

The ecological niche models (ENM) consistently showed that these avian species might substantially contract their ranges to the south during the Last Interglacial period (LIG) and expanded their northern range margins through the Last Glacial Maximum (LGM), leading to the LGM ranges observed for all 11 species. Consistently, coalescent simulations based on 25-30 nuclear genes retrieved signatures of significant population growth through the last glacial period across all species studied. Climate statistics suggested that high climatic variability during the LIG and a relatively mild climate at the LGM potentially explained the historical population dynamics of these birds.

CONCLUSIONS

This is the first study based on multiple species and both lines of ecological niche profiles and genetic data to characterize the unique response of East Asian biota to late Pleistocene climate. The present study highlights regional differences in the evolutionary consequence of climate change during the last glacial cycle and implies that global warming might pose a great risk to species in this region given potentially higher climatic variation in the future analogous to that during the LIG.