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

Across two phylogeographic breaks: Quaternary evolutionary history of a mountain aspen (Populus rotundifolia) in the Hengduan Mountains

Biogeographical barriers to gene flow are central to plant phylogeography. In East Asia, plant distribution is greatly influenced by two phylogeographic breaks, the Mekong-Salween Divide and Tanaka-Kaiyong Line, however, few studies have investigated how these barriers affect the genetic diversity of species that are distributed across both. Here we used 14 microsatellite loci and four chloroplast DNA fragments to examine genetic diversity and distribution patterns of 49 populations of Populus rotundifolia, a species that spans both the Mekong-Salween Divide and the Tanaka-Kaiyong Line in southwestern China. Demographic and migration hypotheses were tested using coalescent-based approaches. Limited historical gene flow was observed between the western and eastern groups of P. rotundifolia, but substantial flow occurred across both the Mekong-Salween Divide and Tanaka-Kaiyong Line, manifesting in clear admixture and high genetic diversity in the central group. Wind-borne pollen and seeds may have facilitated the dispersal of P. rotundifolia following prevalent northwest winds in the spring. We also found that the Hengduan Mountains, where multiple genetic barriers were detected, acted on the whole as a barrier between the western and eastern groups of P. rotundifolia. Ecological niche modeling suggested that P. rotundifolia has undergone range expansion since the last glacial maximum, and demographic reconstruction indicated an earlier population expansion around 600 Ka. The phylogeographic pattern of P. rotundifolia reflects the interplay of biological traits, wind patterns, barriers, niche differentiation, and Quaternary climate history. This study emphasizes the need for multiple lines of evidence in understanding the Quaternary evolution of plants in topographically complex areas.

Hybrid de novo genome assembly of the sexually dimorphic Lady Amherst's pheasant

Pheasants are an important group of birds, valued for their economic benefit as poultry birds, game birds, and as ornamental species for their plumage. Lady Amherst's pheasant Chrysolophus amherstiae is an ornamental species, valued for its elaborate and beautiful plumage. In this study, we present a high-quality de novo hybrid genome assembly of C. amherstiae. Previous attempts to sequence the genome of this species resulted in draft-level assemblies, which are not available in the public domain. Using a combination of Illumina short reads and Oxford Nanopore's long-reads, we assembled a high-quality genome of N50 ~3.9 Mb and near complete BUSCO assessment. We observed a correlation between effective population size and past climatic conditions, with an increase in population size during the warm interglacial periods. We further observed significant fluctuations in genes involved with the immune system and visual perception. C. amherstiae is a highly dimorphic species, and significant fluctuations in gene families involved in immune response, visual perception, among others, suggesting a role of mate choice and sexual selection in the evolution and maintenance of exaggerated traits in the males.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Geographic isolation and environmental heterogeneity contribute to genetic differentiation in Cephalotaxus oliveri

Evaluating the contributions of geographic distance and environmental heterogeneity to the genetic divergence can inform the demographic history and responses to environmental change of natural populations. The isolation-by-distance (IBD) reveals that genetic differentiation among populations increases with geographic distance, while the isolation-by-environment (IBE) assumes a linear relationship between genetic variation and environmental differences among populations. Here, we sampled and genotyped 330 individuals from 18 natural populations of Cephalotaxus oliveri throughout the species' distribution. Twenty-eight EST-SSR markers were applied to analyze population genetics, for the investigation of the driving factors that shaped spatial structure. In addition, we identified the outlier loci under positive selection and tested their association with environmental factors. The results showed a moderate genetic diversity in C. oliveri and high genetic differentiation among populations. Population structure analyses indicated that 18 populations were clustered into two major groups. We observed that the genetic diversity of central populations decreased and the genetic differentiation increased towards the marginal populations. Additionally, the signatures of IBD and IBE were detected in C. oliveri, and IBE provided a better contribution to genetic differentiation. Six outlier loci under positive selection were demonstrated to be closely correlated with environmental variables, among which bio8 was associated with the greatest number of loci. Genetic evidence suggests the consistency of the central-marginal hypothesis (CMH) for C. oliveri. Furthermore, our results suggest that temperature-related variables played an important role in shaping genetic differentiation.

Genomic signatures and evolutionary history of the endangered blue-crowned laughingthrush and other Garrulax species

BACKGROUND

The blue-crowned laughingthrush (Garrulax courtoisi) is a critically endangered songbird endemic to Wuyuan, China, with population of ~323 individuals. It has attracted widespread attention, but the lack of a published genome has limited research and species protection.

RESULTS

We report two laughingthrush genome assemblies and reveal the taxonomic status of laughingthrush species among 25 common avian species according to the comparative genomic analysis. The blue-crowned laughingthrush, black-throated laughingthrush, masked laughingthrush, white-browed laughingthrush, and rusty laughingthrush showed a close genetic relationship, and they diverged from a common ancestor between ~2.81 and 12.31 million years ago estimated by the population structure and divergence analysis using 66 whole-genome sequencing birds from eight laughingthrush species and one out group (Cyanopica cyanus). Population inference revealed that the laughingthrush species experienced a rapid population decline during the last ice age and a serious bottleneck caused by a cold wave during the Chinese Song Dynasty (960-1279 AD). The blue-crowned laughingthrush is still in a bottleneck, which may be the result of a cold wave together with human exploitation. Interestingly, the existing blue-crowned laughingthrush exhibits extremely rich genetic diversity compared to other laughingthrushes. These genetic characteristics and demographic inference patterns suggest a genetic heritage of population abundance in the blue-crowned laughingthrush. The results also suggest that fewer deleterious mutations in the blue-crowned laughingthrush genomes have allowed them to thrive even with a small population size. We believe that cooperative breeding behavior and a long reproduction period may enable the blue-crowned laughingthrush to maintain genetic diversity and avoid inbreeding depression. We identified 43 short tandem repeats that can be used as markers to identify the sex of the blue-crowned laughingthrush and aid in its genetic conservation.

CONCLUSIONS

This study supplies the missing reference genome of laughingthrush, provides insight into the genetic variability, evolutionary potential, and molecular ecology of laughingthrush and provides a genomic resource for future research and conservation.

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

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

Evolutionary history of two rare endemic conifer species from the eastern Qinghai-Tibet Plateau

BACKGROUND AND AIMS

Understanding the population genetics and evolutionary history of endangered species is urgently needed in an era of accelerated biodiversity loss. This knowledge is most important for regions with high endemism that are ecologically vulnerable, such as the Qinghai-Tibet Plateau (QTP).

METHODS

The genetic variation of 84 juniper trees from six populations of Juniperus microsperma and one population of Juniperus erectopatens, two narrow-endemic junipers from the QTP that are sister to each other, was surveyed using RNA-sequencing data. Coalescent-based analyses were used to test speciation, migration and demographic scenarios. Furthermore, positively selected and climate-associated genes were identified, and the genetic load was assessed for both species.

KEY RESULTS

Analyses of 149 052 single nucleotide polymorphisms showed that the two species are well differentiated and monophyletic. They diverged around the late Pliocene, but interspecific gene flow continued until the Last Glacial Maximum. Demographic reconstruction by Stairway Plot detected two severe bottlenecks for J. microsperma but only one for J. erectopatens. The identified positively selected genes and climate-associated genes revealed habitat adaptation of the two species. Furthermore, although J. microsperma had a much wider geographical distribution than J. erectopatens, the former possesses lower genetic diversity and a higher genetic load than the latter.

CONCLUSIONS

This study sheds light on the evolution of two endemic juniper species from the QTP and their responses to Quaternary climate fluctuations. Our findings emphasize the importance of speciation and demographic history reconstructions in understanding the current distribution pattern and genetic diversity of threatened species in mountainous regions.

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.

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

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

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.

Evolution of strong reproductive isolation in plants: broad-scale patterns and lessons from a perennial model group

Many recent studies have addressed the mechanisms operating during the early stages of speciation, but surprisingly few studies have tested theoretical predictions on the evolution of strong reproductive isolation (RI). To help address this gap, we first undertook a quantitative review of the hybrid zone literature for flowering plants in relation to reproductive barriers. Then, using Populus as an exemplary model group, we analysed genome-wide variation for phylogenetic tree topologies in both early- and late-stage speciation taxa to determine how these patterns may be related to the genomic architecture of RI. Our plant literature survey revealed variation in barrier complexity and an association between barrier number and introgressive gene flow. Focusing on Populus, our genome-wide analysis of tree topologies in speciating poplar taxa points to unusually complex genomic architectures of RI, consistent with earlier genome-wide association studies. These architectures appear to facilitate the 'escape' of introgressed genome segments from polygenic barriers even with strong RI, thus affecting their relationships with recombination rates. Placed within the context of the broader literature, our data illustrate how phylogenomic approaches hold great promise for addressing the evolution and temporary breakdown of RI during late stages of speciation. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Large-Scale Hybridisation as an Extinction Threat to the Suweon Treefrog (Hylidae: Dryophytes suweonensis)

Simple Summary

A large number of amphibian species are now endangered, mostly because of human activities. An example is land modification, which may bring species that were previously isolated in contact, and allows them to hybridise. Here, we assessed the presence of hybrid individuals between the endangered Suweon treefrog (Dryophytes suweonensis) and the widespread Japanese treefrog (Dryophytes japonicus). We found hybrids to be relatively widespread and present at all populations where the Suweon treefrog occurred. This is important, as it results in an additional threat to the Suweon treefrog.

Abstract

Amphibians are in the midst of a sixth mass extinction, and human activities play a major role in pushing species towards extinction. Landscape anthropisation has impacts that indirectly threaten species, in addition to the obvious destruction of natural habitats. For instance, land modification may bring human-commensal species in contact with sister-clades from which they were previously isolated. The species in these new contact zones are then able to hybridise to the point of reaching lineage fusion, through which the gene pool of the two species merges and one of the parental lineages becomes extirpated. Here, we documented the patterns of hybridisation between the spatially restricted D. suweonensis and the widespread D. japonicus. On the basis of the analysis of Cytochrome c oxidase subunit I mitochondrial DNA sequences (404 individuals from 35 sites) and six polymorphic microsatellites (381 individuals from 34 sites), we revealed a generalised, bi-directional, and geographically widespread hybridisation between the two species. Evidence of fertile back-crosses is provided by relatively high numbers of individuals in cyto-nuclear disequilibrium, as well as the presence of hybrid individuals further south than the species distribution limit, determined on the basis of call properties. Hybridisation is an additional threat to the endangered D. suweonensis.

Genetic Diversity and Population Structure of Rhododendron rex Subsp. rex Inferred from Microsatellite Markers and Chloroplast DNA Sequences

Genetic diversity is vital to the sustainable utilization and conservation of plant species. Rhododendron rex subsp. rex Lévl. is an endangered species endemic to the southwest of China. Although the natural populations of this species are facing continuous decline due to the high frequency of anthropogenic disturbance, the genetic information of R. rex subsp. rex is not yet elucidated. In the present study, 10 pairs of microsatellite markers (nSSRs) and three pairs of chloroplast DNA (cpDNAs) were used in the elucidation of the genetic diversity, population structure, and demographic history of 11 R. rex subsp. rex populations. A total of 236 alleles and 12 haplotypes were found. A moderate genetic diversity within populations (HE = 0.540 for nSSRs, Hd = 0.788 for cpDNA markers), high historical and low contemporary gene flows, and moderate genetic differentiation (nSSR: FST = 0.165***; cpDNA: FST = 0.841***) were detected among the R. rex subsp. rex populations. Genetic and geographic distances showed significant correlation (p < 0.05) determined by the Mantel test. The species exhibited a conspicuous phylogeographical structure among the populations. Using the Bayesian skyline plot and species distribution models, we found that R. rex subsp. rex underwent a population demography contraction approximately 50,000-100,000 years ago. However, the species did not experience a recent population expansion event. Thus, habitat loss and destruction, which result in a population decline and species inbreeding depression, should be considered in the management and conservation of R. rex subsp. rex.

Climatic Change Can Influence Species Diversity Patterns and Potential Habitats of Salicaceae Plants in China

Salicaceae is a family of temperate woody plants in the Northern Hemisphere that are highly valued, both ecologically and economically. China contains the highest species diversity of these plants. Despite their widespread human use, how the species diversity patterns of Salicaceae plants formed remains mostly unknown, and these may be significantly affected by global climate warming. Using past, present, and future environmental data and 2673 georeferenced specimen records, we first simulated the dynamic changes in suitable habitats and population structures of Salicaceae. Based on this, we next identified those areas at high risk of habitat loss and population declines under different climate change scenarios/years. We also mapped the patterns of species diversity by constructing niche models for 215 Salicaceae species, and assessed the driving factors affecting their current diversity patterns. The niche models showed Salicaceae family underwent extensive population expansion during the Last Inter Glacial period but retreated to lower latitudes during and since the period of the Last Glacial Maximum. Looking ahead, as climate warming intensifies, suitable habitats will shift to higher latitudes and those at lower latitudes will become less abundant. Finally, the western regions of China harbor the greatest endemism and species diversity of Salicaceae, which are significantly influenced by annual precipitation and mean temperature, ultraviolet-B (UV-B) radiation, and the anomaly of precipitation seasonality. From these results, we infer water–energy dynamic equilibrium and historical climate change are both the main factors likely regulating contemporary species diversity and distribution patterns. Nevertheless, this work also suggests that other, possibly interacting, factors (ambient energy, disturbance history, soil condition) influence the large-scale pattern of Salicaceae species diversity in China, making a simple explanation for it unlikely. Because Southwest China likely served as a refuge for Salicaceae species during the Last Glacial Maximum, it is a current hotspot for endemisms. Under predicted climate change, Salicaceae plants may well face higher risks to their persistence in southwest China, so efforts to support their in-situ conservation there are urgently needed.

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.

PHENETIC ANALYSIS OF POPULUS NIGRA, P. LAURIFOLIA AND P. × JRTYSCHENSIS IN NATURAL HYBRIDIZATION ZONE

The wide spread of hybridization in the genus Populus, including spontaneous hybridization, caused by cultivars, requires studying the variability and inheritance of morphological traits by hybrids for initial tracking of these processes. The analysis of endogenous, intra- and inter-population variability was performed on 533 individual poplar trees in seven populations of P. nigra, seven populations of P. laurifolia and four populations of P. × jrtyschensis in the Tom river basin. On each specimen, 15 leaves from short mid-crown branches were collected to determine the shape of the leaf blade, the shape of its tip and base, as well as the branch morphotype. Some biometric indicators were proposed for geometric assessment of the leaf blade shapes of poplar species. The analysis showed that of all the traits examined only the leaf blade shape did not meet the criterion for “phene”, since it is usually represented by several forms in the crown of one and the same tree. In all the species studied, the level of their intra-population diversity was found to be much higher than the inter-population one. According to the increase of intra-population variability of qualitative traits, the taxa could be ranked as P. nigra < P. laurifolia < P. × jrtyschensis. The share of inter-population diversity differed among the species studied, accounting for 21.5% in P. laurifolia, 13.8% in P. nigra and 8.0% in P. × jrtyschensis. The P. laurifolia populations showed the greatest inter-population differentiation, most likely because of a disjunct distribution due to narrow specialization to the montane river environment. The lower differentiation in P. nigra is probably due to the facts that this species dominates the poplar stands of the Tom River basin; its populations are not fragmented and are linked by vast gene flow. In P. nigra, an increase in the diversity of qualitative characteristics and phenotypes was observed in populations confined to hybridization centers. Natural selection is most likely the factor governing the inter-population differentiation in P. × jrtyschensis, leading to the predominance of F1 hybrids in populations and hence to a sharp decrease in inter-population variability.