The role of East Asian monsoon system in shaping population divergence and dynamics of a constructive desert shrub Reaumuria soongarica

RADseq-based population genomic analysis and environmental adaptation of rare and endangered recretohalophyte Reaumuria trigyna

Genetic diversity reflects the survival potential, history, and population dynamics of an organism. It underlies the adaptive potential of populations and their response to environmental change. Reaumuria trigyna is an endemic species in the Eastern Alxa and West Ordos desert regions in China. The species has been considered a good candidate to explore the unique survival strategies of plants that inhabit this area. In this study, we performed population genomic analyses based on restriction-site associated DNA sequencing to understand the genetic diversity, population genetic structure, and differentiation of the species. Analyses of 92,719 high-quality single-nucleotide polymorphisms (SNPs) indicated that overall genetic diversity of R. trigyna was low (HO = 0.249 and HE = 0.208). No significant genetic differentiation was observed among the investigated populations. However, a subtle population genetic structure was detected. We suggest that this might be explained by adaptive diversification reinforced by the geographical isolation of populations. Overall, 3513 outlier SNPs were located in 243 gene-coding sequences in the R. trigyna transcriptome. Potential sites under diversifying selection occurred in genes (e.g., AP2/EREBP, E3 ubiquitin-protein ligase, FLS, and 4CL) related to phytohormone regulation and synthesis of secondary metabolites which have roles in adaptation of species. Our genetic analyses provide scientific criteria for evaluating the evolutionary capacity of R. trigyna and the discovery of unique adaptions. Our findings extend knowledge of refugia, environmental adaption, and evolution of germplasm resources that survive in the Ordos area.

Demographic patterns of two related desert shrubs with overlapping distributions in response to past climate changes

Numerous studies have revealed that past geological events and climatic fluctuations had profoundly affected the genetic structure and demographic patterns of species. However, related species with overlapping ranges may have responded to such environmental changes in different ways. In this study, we compared the genetic structure and population dynamics of two typical desert shrubs with overlapping distributions in northern China, Nitraria tangutorum and Nitraria sphaerocarpa, based on chloroplast DNA (cpDNA) variations and species distribution models. We sequenced two cpDNA fragments (trnH-trnA and atpH-atpI) in 633 individuals sampled from 52 natural populations. Twenty-four chlorotypes, including eight rare chlorotypes, were identified, and a single dominant haplotype (H4) widely occurred in the entire geographical ranges of the two species. There were also a few distinctive chlorotypes fixed in different geographical regions. Population structure analyses suggested that the two species had significantly different levels of total genetic diversity and interpopulation differentiation, which was highly likely correlated with the special habitat preferences of the two species. A clear phylogeographic structure was identified to exist among populations of N. sphaerocarpa, but not exist for N. tangutorum. The neutral tests, together with the distribution of pairwise differences revealed that N. tangutorum experienced a sudden demographic expansion, and its expansion approximately occurred between 21 and 7 Kya before present, while a rapid range expansion was not identified for N. sphaerocarpa. The ecological niche modeling (ENM) analysis indicated that the potential ranges of two species apparently fluctuated during the past and present periods, with obvious contraction in the Last Glacial Maximum (LGM) and recolonization in the present, respectively, comparing to the Last Interglacial (LIG). These findings suggest that the two species extensively occurred in the Northwest of China before the Quaternary, and the current populations of them originated from a few separated glacial refugia following their habitat fragmentation in the Quarternary. Our results provide new insights on the impact of past geological and climatic fluctuations on the population dynamics of desert plants in northwestern China, and further enforce the hypothesis that there were several independent glacial refugia for these species during the Quaternary glaciations.

Gene Profiling of the Ascorbate Oxidase Family Genes under Osmotic and Cold Stress Reveals the Role of AnAO5 in Cold Adaptation in Ammopiptanthus nanus

The uplift of the Qinghai Tibet Plateau has led to a drastic change in the climate in Central Asia, from warm and rainy, to dry and less rainfall. Ammopiptanthus nanus, a rare evergreen broad-leaved shrub distributed in the temperate desert region of Central Asia, has survived the drastic climate change in Central Asia caused by the uplift of the Qinghai-Tibet Plateau. Ascorbate oxidase (AO) regulates the redox status of the apoplast by catalyzing the oxidation of ascorbate acid to dehydroascorbic acid, and plays a key role in the adaptation of plants to environmental changes. Analyzing the evolution, environmental response, and biological functions of the AO family of A. nanus is helpful for understanding how plant genome evolution responds to climate change in Central Asia. A total of 16 AOs were identified in A. nanus, all of which contained the ascorbate oxidase domain, most of which contained transmembrane domain, and many were predicted to be localized in the apoplast. Segmental duplication and tandem duplication are the main factors driving the gene amplification of the AO gene family in A. nanus. Gene expression analysis based on transcriptome data and fluorescence quantitative PCR, as well as enzyme activity measurements, showed that the expression levels of AO genes and total enzyme activity decreased under short-term osmotic stress and low-temperature stress, but the expression of some AO genes (AnAO5, AnAO13, and AnAO16) and total enzyme activity increased under 7 days of cold stress. AnAO5 and AnAO11 are targeted by miR4415. Further functional studies on AnAO5 showed that AnAO5 protein was localized in the apoplast. The expression of AnAO5 in yeast cells and the transient expression in tobacco enhanced the tolerance of yeast and tobacco to low-temperature stress, and the overexpression of AnAO5 enhanced the tolerance of Arabidopsis seedlings to cold stress. Our research provides important data for understanding the role of AOs in plant adaptation to environmental change.

Phylogeography and Population Genetics Analyses Reveal Evolutionary History of the Desert Resource Plant Lycium ruthenicum (Solanaceae)

Climactic oscillations during the Quaternary played a significant role in the formation of genetic diversity and historical demography of numerous plant species in northwestern China. In this study, we used 11 simple sequence repeats derived from expressed sequence tag (EST-SSR), two chloroplast DNA (cpDNA) fragments, and ecological niche modeling (ENM) to investigate the population structure and the phylogeographic history of Lycium ruthenicum, a plant species adapted to the climate in northwestern China. We identified 20 chloroplast haplotypes of which two were dominant and widely distributed in almost all populations. The species has high haplotype diversity and low nucleotide diversity based on the cpDNA data. The EST-SSR results showed a high percentage of total genetic variation within populations. Both the cpDNA and EST-SSR results indicated no significant differentiation among populations. By combining the evidence from ENM and demographic analysis, we confirmed that both the last interglacial (LIG) and late-glacial maximum (LGM) climatic fluctuations, aridification might have substantially narrowed the distribution range of this desert species, the southern parts of the Junggar Basin, the Tarim Basin, and the eastern Pamir Plateau were the potential glacial refugia for L. ruthenicum during the late middle Pleistocene to late Pleistocene Period. During the early Holocene, the warm, and humid climate promoted its demographic expansion in northwestern China. This work may provide new insights into the mechanism of formation of plant diversity in this arid region.

Evolutionary history of a desert perennial Arnebia szechenyi (Boraginaceae): Intraspecific divergence, regional expansion and asymmetric gene flow

The complex interactions of historical, geological and climatic events on plant evolution have been an important research focus for many years. However, the role of desert formation and expansion in shaping the genetic structures and demographic histories of plants occurring in arid areas has not been well explored. In the present study, we investigated the phylogeography of Arnebia szechenyi, a desert herb showing a near-circular distribution surrounding the Tengger Desert in Northwest China. We measured genetic diversity of populations using three maternally inherited chloroplast DNA (cpDNA) fragments and seven bi-paternally inherited nuclear DNA (nDNA) loci that were sequenced from individuals collected from 16 natural populations across its range and modelled current and historical potential habitats of the species. Our data indicated a considerably high level of genetic variation within A. szechenyi and noteworthy asymmetry in historical migration from the east to the west. Moreover, two nuclear genetic groups of populations were revealed, corresponding to the two geographic regions separated by the Tengger Desert. However, analysis of cpDNA data did not show significant geographic structure. The most plausible explanation for the discrepancy between our findings based on cpDNA and nDNA data is that A. szechenyi populations experienced long periods of geographic isolation followed by range expansion, which would have promoted generalized recombination of the nuclear genome. Our findings further highlight the important role that the Tengger Desert, together with the Helan Mountains, has played in the evolution of desert plants and the preservation of biodiversity in arid Northwest China.

Genomic Adaptive Evolution of Sand Rice (Agriophyllum squarrosum) and Its Implications for Desert Ecosystem Restoration

Natural selection is a significant driver of population divergence and speciation of plants. Due to local adaptation to geographic regions with ecological gradients, plant populations harbored a wide range of adaptive genetic variation to enable them to survive the heterogeneous habitats. This is all the more necessary for desert plants, as they must tolerant more striking gradients of abiotic stresses. However, the genomic mechanism by which desert plants adapt to ecological heterogeneity remains unclear, which could help to guide the sustainability of desert ecosystems. Here, using restriction-site-associated DNA sequencing in 38 natural populations, we investigated the genomic divergence and environmental adaptation of sand rice, Agriophyllum squarrosum, an annual pioneer species that covers sand dunes in northern China. Population genetic structure analyses showed that sand rice could be divided into three geographically distinct lineages, namely, Northwest, Central, and East. Phylogeographic analyses revealed that the plant might originate locally in Bergen County and further differentiated into the East lineage and then the Central lineage. Ecological niche modeling found that different lineages occupied distinct ecological niches, suggesting that the ecological gradient would have triggered genomic differentiation among sand rice lineages. Ecological association study supported that the three SNPs under divergent selection were closely correlated with precipitation gradients, indicating that precipitation might be the most important stress trigger for lineage diversity in sand rice. These adaptive SNPs could be used to genotype suitable germplasms for the ecological restoration of specific desertified lands. Further analyses found that genetic structure could significantly overestimate the signals for balancing selection. Within the Central lineage, we still found that 175 SNPs could be subject to balancing selection, which could be the means by which sand rice maintains genetic diversity and adapts to multiple stresses across heterogeneous deserts and sandy lands. From a genomic point of view, this study highlighted the local and global adaptation patterns of a desert plant to extreme and heterogeneous habitats. Our data provide molecular guidance for the restoration of desertified lands in the arid and semi-arid regions of China and could facilitate the marker assistant breeding of this potential crop to mitigate climate change.

Leaf size variations in a dominant desert shrub, Reaumuria soongarica, adapted to heterogeneous environments

The climate in arid Central Asia (ACA) has changed rapidly in recent decades, but the ecological consequences of this are far from clear. To predict the impacts of climate change on ecosystem functioning, greater attention should be given to the relationships between leaf functional traits and environmental heterogeneity. As a dominant constructive shrub widely distributed in ACA, Reaumuria soongarica provided us with an ideal model to understand how leaf functional traits of desert ecosystems responded to the heterogeneous environments of ACA. Here, to determine the influences of genetic and ecological factors, we characterized species-wide variations in leaf traits among 30 wild populations of R. soongarica and 16 populations grown in a common garden. We found that the leaf length, width, and leaf length to width ratio (L/W) of the northern lineage were significantly larger than those of other genetic lineages, and principal component analysis based on the in situ environmental factors distinguished the northern lineage from the other lineages studied. With increasing latitude, leaf length, width, and L/W in the wild populations increased significantly. Leaf length and L/W were negatively correlated with altitude, and first increased and then decreased with increasing mean annual temperature (MAT) and mean annual precipitation (MAP). Stepwise regression analyses further indicated that leaf length variation was mainly affected by latitude. However, leaf width was uncorrelated with altitude, MAT, or MAP. The common garden trial showed that leaf width variation among the eastern populations was caused by both local adaptation and phenotypic plasticity. Our findings suggest that R. soongarica preferentially changes leaf length to adjust leaf size to cope with environmental change. We also reveal phenotypic evidence for ecological speciation of R. soongarica. These results will help us better understand and predict the consequences of climate change for desert ecosystem functioning.

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.

Dispersal and genetic differentiation of Syntrichia caninervis populations across different desert regions in China

The moss Syntrichia caninervis is widely distributed in cool temperate and cold deserts where environmental pressures create a dependence on asexual reproduction (fragment reproduction). However, when compared to sporophyte-producing mosses, there is a lack of evidence to support the capacity of drought-tolerant mosses that predominantly fragment and produce protonema to disperse over long distances. We used 20 microsatellite loci to study genetic variation and structure in six populations (five natural populations and one population from a regeneration site) in three contrasting and widely separated regions of China. The genetic diversity and expected heterozygosity were lower in populations from the Tengger Desert than in populations from the other regions. Using PCoA, UPGMA and Structure analysis, the genetic grouping divided the three regions into three distinct groups. This may indicate that in regions where S. caninervis reproduces predominantly asexually, propagules are spread mainly by short-distance dispersal. The genetic diversity of the population from the regeneration site in the Tengger Desert was slightly higher than that of the nearby, naturally occurring population, and included some input from the Pamir Plateau almost 2,300 km to the west, suggesting long-distance dispersal of S. caninervis propagules across the region. Predominantly asexually reproducing populations of S. caninervis are mainly dependent on short-distance dispersal. Long-distance dispersal of S. caninervis propagules across the region is difficult. Establishment of populations with dominant asexual reproduction will eventually result in genetic differentiation.

Phylogeographic patterns of the desert poplar in Northwest China shaped by both geology and climatic oscillations

Background

The effects of historical geology and climatic events on the evolution of plants around the Qinghai-Tibetan Plateau region have been at the center of debate for years. To identify the influence of the uplift of the Tianshan Mountains and/or climatic oscillations on the evolution of plants in arid northwest China, we investigated the phylogeography of the Euphrates poplar (Populus euphratica) using chloroplast DNA (cpDNA) sequences and nuclear microsatellites, and estimated its historical distribution using Ecological Niche Modeling (ENM).

Results

We found that the Euphrates poplar differed from another desert poplar, P. pruinosa, in both nuclear and chloroplast DNA. The low clonal diversity in both populations reflected the low regeneration rate by seed/seedlings in many locations. Both cpDNA and nuclear markers demonstrated a clear divergence between the Euphrates poplar populations from northern and southern Xinjiang regions. The divergence time was estimated to be early Pleistocene based on cpDNA, and late Pleistocene using an Approximate Bayesian Computation analysis based on microsatellites. Estimated gene flow was low between these two regions, and the limited gene flow occurred mainly via dispersal from eastern regions. ENM analysis supported a wider distribution of the Euphrates poplar at 3 Ma, but a more constricted distribution during both the glacial period and the interglacial period.

Conclusions

These results indicate that the deformation of the Tianshan Mountains has impeded gene flow of the Euphrates poplar populations from northern and southern Xinjiang, and the distribution constriction due to climatic oscillations further accelerated the divergence of populations from these regions. To protect the desert poplars, more effort is needed to encourage seed germination and seedling establishment, and to conserve endemic gene resources in the northern Xinjiang region.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1194-1) contains supplementary material, which is available to authorized users.

Population dynamics of Agriophyllum squarrosum, a pioneer annual plant endemic to mobile sand dunes, in response to global climate change

Climate change plays an important role in the transition of ecosystems. Stratigraphic investigations have suggested that the Asian interior experienced frequent transitions between grassland and desert ecosystems as a consequence of global climate change. Using maternally and bi-parentally inherited markers, we investigated the population dynamics of Agriophyllum squarrosum (Chenopodiaceae), an annual pioneer plant endemic to mobile sand dunes. Phylogeographic analysis revealed that A. squarrosum could originate from Gurbantunggut desert since ~1.6 Ma, and subsequently underwent three waves of colonisation into other deserts and sandy lands corresponding to several glaciations. The rapid population expansion and distribution range shifts of A. squarrosum from monsoonal climate zones suggested that the development of the monsoonal climate significantly enhanced the population growth and gene flow of A. squarrosum. These data also suggested that desertification of the fragile grassland ecosystems in the Qinghai-Tibetan Plateau was more ancient than previously suggested and will be aggravated under global warming in the future. This study provides new molecular phylogeographic insights into how pioneer annual plant species in desert ecosystems respond to global climate change, and facilitates evaluation of the ecological potential and genetic resources of future crops for non-arable dry lands to mitigate climate change.

Different Sets of Post-Embryonic Development Genes Are Conserved or Lost in Two Caryophyllales Species (Reaumuria soongorica and Agriophyllum squarrosum)

Reaumuria soongorica and sand rice (Agriophyllum squarrosum) belong to the clade of Caryophyllales and are widely distributed in the desert regions of north China. Both plants have evolved many specific traits and adaptation strategies to cope with recurring environmental threats. However, the genetic basis that underpins their unique traits and adaptation remains unknown. In this study, the transcriptome data of R. soongorica and sand rice were compared with three other species with previously sequenced genomes (Arabidopsis thaliana, Oryza sativa, and Beta vulgaris). Four different gene sets were identified, namely, the genes conserved in both species, those lost in both species, those conserved in R. soongorica only, and those conserved in sand rice only. Gene ontology showed that post-embryonic development genes (PEDGs) were enriched in all gene sets, and different sets of PEDGs were conserved or lost in both the R. soongorica and sand rice genomes. Expression profiles of Arabidopsis orthologs further provided some clues to the function of the species-specific conserved PEDGs. Such orthologs included LEAFY PETIOLE, which could be a candidate gene involved in the development of branch priority in sand rice.

Different Sets of Post-Embryonic Development Genes Are Conserved or Lost in Two Caryophyllales Species (Reaumuria soongorica and Agriophyllum squarrosum)

Reaumuria soongorica and sand rice (Agriophyllum squarrosum) belong to the clade of Caryophyllales and are widely distributed in the desert regions of north China. Both plants have evolved many specific traits and adaptation strategies to cope with recurring environmental threats. However, the genetic basis that underpins their unique traits and adaptation remains unknown. In this study, the transcriptome data of R. soongorica and sand rice were compared with three other species with previously sequenced genomes (Arabidopsis thaliana, Oryza sativa, and Beta vulgaris). Four different gene sets were identified, namely, the genes conserved in both species, those lost in both species, those conserved in R. soongorica only, and those conserved in sand rice only. Gene ontology showed that post-embryonic development genes (PEDGs) were enriched in all gene sets, and different sets of PEDGs were conserved or lost in both the R. soongorica and sand rice genomes. Expression profiles of Arabidopsis orthologs further provided some clues to the function of the species-specific conserved PEDGs. Such orthologs included LEAFY PETIOLE, which could be a candidate gene involved in the development of branch priority in sand rice.