Genetic differentiation and delimitation between ecologically diverged Populus euphratica and P. pruinosa

Coalescent-based species delimitation in North American pinyon pines using low-copy nuclear genes and plastomes

PREMISE

Accurate species delimitation is essential for evolutionary biology, conservation, and biodiversity management. We studied species delimitation in North American pinyon pines, Pinus subsection Cembroides, a natural group with high levels of incomplete lineage sorting.

METHODS

We used coalescent-based methods and multivariate analyses of low-copy number nuclear genes and nearly complete high-copy number plastomes generated with the Hyb-Seq method. The three coalescent-based species delimitation methods evaluated were the Generalized Mixed Yule Coalescent (GMYC), Poisson Tree Process (PTP), and Trinomial Distribution of Triplets (Tr2). We also measured admixture in populations with possible introgression.

RESULTS

Our results show inconsistencies among GMYC, PTP, and Tr2. The single-locus based GMYC analysis of plastid DNA recovered a higher number of species (up to 24 entities, including singleton lineages and clusters) than PTP and the multi-locus coalescent approach. The PTP analysis identified 10 species whereas Tr2 recovered 13, which agreed closely with taxonomic treatments.

CONCLUSIONS

We found that PTP and GMYC identified species with low levels of ILS and high morphological divergence (P. maximartinezii, P. pinceana, and P. rzedowskii). However, GMYC method oversplit species by identification of more divergent samples as singletons. Moreover, both PTP and GMYC were incapable of identifying some species that are readily identified morphologically. We suggest that the divergence times between lineages within North American pinyon pines are so disparate that GMYC results are unreliable. Results of the Tr2 method coincided well with previous delimitations based on morphology, DNA, geography, and secondary chemistry.

Phylogeography Reveals Geographic and Environmental Factors Driving Genetic Differentiation of Populus sect. Turanga in Northwest China

Populus sect. Turanga (hereafter referred to as "Populus"), including Populus euphratica and Populus pruinosa, are the predominant tree species in desert riparian forests in northwestern China. These trees play key roles in maintaining ecosystem balance, curbing desertification, and protecting biodiversity. However, the distribution area of Populus forests has been severely diminished and degraded in recent years due to increased habitat destruction and human activity. Understanding the genetic diversity among Populus individuals and populations is essential for designing conservation strategies, but comprehensive studies of their genetic diversity in northwest China are lacking. Here, we assessed the population structures and genetic diversity of 1,620 samples from 85 natural populations of Populus (59 P. euphratica and 26 P. pruinosa populations) covering all of northwestern China using 120 single nucleotide polymorphism (SNP) markers. Analysis of population structure revealed significant differentiation between these two sister species and indicated that strong geographical distribution patterns, a geographical barrier, and environmental heterogeneity shaped the extant genetic patterns of Populus. Both P. euphratica and P. pruinosa populations in southern Xinjiang had higher genetic diversity than populations in other clades, perhaps contributing to local geographic structure and strong gene flow. Analysis of molecular variance (AMOVA) identified 15% variance among and 85% variance within subpopulations. Mantel tests suggested that the genetic variation among P. euphratica and P. pruinosa populations could be explained by both geographical and environmental distance. The genetic diversity of P. euphratica showed a significant negative correlation with latitude and longitude and a positive correlation with various environmental factors, such as precipitation of warmest quarter and driest month, temperature seasonality, and annual mean temperature. These findings provide insights into how the genetic differentiation of endangered Populus species was driven by geographical and environmental factors, which should be helpful for designing strategies to protect these genetic resources in the future.

Demographic history and genetic differentiation of an endemic and endangered Ulmus lamellosa (Ulmus)

BACKGROUND

Ulmus lamellosa (one of the ancient species of Ulmus) is an endemic and endangered plant that has undergone climatic oscillations and geographical changes. The elucidation of its demographic history and genetic differentiation is critical for understanding the evolutionary process and ecological adaption to forests in Northern China.

RESULTS

Polymorphic haplotypes were detected in most populations of U. lamellosa via DNA sequencing. All haplotypes were divided into three phylogeographic clades fundamentally corresponding to their geographical distribution, namely THM (Taihang Mountains), YM (Yinshan Mountains), and YSM (Yanshan Mountains) groups. The YSM group, which is regarded as ancestral, possessed higher genetic diversity and significant genetic variability in contrast to the YSM and YM groups. Meanwhile, the divergence time of intraspecies haplotypes occurred during the Miocene-Pliocene, which was associated with major Tertiary geological and/or climatic events. Different degrees of gene exchanges were identified between the three groups. During glaciation, the YSM and THM regions might have served as refugia for U. lamellosa. Based on ITS data, range expansion was not expected through evolutionary processes, except for the THM group. A series of mountain uplifts (e.g., Yanshan Mountains and Taihang Mountains) following the Miocene-Pliocene, and subsequently quaternary climatic oscillations in Northern China, further promoted divergence between U. lamellosa populations.

CONCLUSIONS

Geographical topology and climate change in Northern China played a critical role in establishing the current phylogeographic structural patterns of U. lamellosa. These results provide important data and clues that facilitate the demographic study of tree species in Northern China.

Change in Spatial Distribution Patterns and Regeneration of Populus euphratica under Different Surface Soil Salinity Conditions

Ecological conservation and restoration have increasingly captured attention worldwide due to the degradation of ecosystems. As one of the most ecologically fragile areas, the Tarim River Basin, of Xinjiang, China, encountered serious decline of desert riparian forests. The Chinese government has implemented the "Ecological Water Conveyance Project" (EWCP) since 2000, protecting and restoring the dominant natural species of the desert riparian forests, i.e., Populus euphratica Oliv. The regenerative effect after the water conveyance was noteworthy. For the purpose of clarifying the mechanism of P. euphratica forest regeneration to find a better prescription for the ecological restoration works in the Tarim River Basin, we investigated the relationship between the distribution of P. euphratica and soil salinity. Experimentally evaluated the effects of surface soil salinity on P. euphratica seed germination and the influence of river flooding on the salinity of surface soils. The results showed that (1) P. euphratica trees mainly spread along the river channel within 2 km; with increasing vertical distance to the channel, the number of trees declined significantly; (2) where the salinity of the surface soil is high, there are less living P. euphratica trees; (3) the germination of P. euphratica seeds decreases with increased soil conductivity; when the soil conductivity was higher than 7 ms/cm, the germination of P. euphratica seeds was severely curtailed. (4) Flooding regimes were a pre-condition of P. euphratica restoration; they had profound effects on improving the germination of the seeds via ameliorating water conditions and reducing salinity. Our results point out that the most efficient ecological prescription for restoring and protecting desert riparian forests is to induce flooding twice yearly during June to August with 10- to 15-day durations each time. Such a plan (especially in the Tarim River Basin) should prioritize the protection of seedlings.

Transcriptomic Analysis of Seed Germination Under Salt Stress in Two Desert Sister Species (Populus euphratica and P. pruinosa)

As a major abiotic stress, soil salinity limits seed germination and plant growth, development and production. Seed germination is highly related not only to the seedlings survival rate but also subsequent vegetative growth. Populus euphratica and P. pruinosa are closely related species that show a distinguished adaptability to salinity stress. In this study, we performed an integrative transcriptome analyses of three seed germination phases from P. euphratica and P. pruinosa under salt stress. A two-dimensional data set of this study provides a comprehensive view of the dynamic biochemical processes that underpin seed germination and salt tolerance. Our analysis identified 12831 differentially expressed genes (DEGs) for seed germination processes and 8071 DEGs for salt tolerance in the two species. Furthermore, we identified the expression profiles and main pathways in each growth phase. For seed germination, a large number of DEGs, including those involved in energy production and hormonal regulation pathways, were transiently and specifically induced in the late phase. In the comparison of salt tolerance between the two species, the flavonoid and brassinosteroid pathways were significantly enriched. More specifically, in the flavonoid pathway, FLS and F3'5'H exhibited significant differential expression. In the brassinosteroid pathway, the expression levels of DWF4, BR6OX2 and ROT3 were notably higher in P. pruinosa than in P. euphratica. Our results describe transcript dynamics and highlight secondary metabolite pathways involved in the response to salt stress during the seed germination of two desert poplars.

Footprints of divergent evolution in two Na+/H+ type antiporter gene families (NHX and SOS1) in the genus Populus

Populus, a deciduous tree species of major economic and ecological value, grows across the range in which trees are distributed in the Northern Hemisphere. Patterns of DNA variation are often used to identify the evolutionary forces shaping the genotypes of distinctive species lineages. Sodium/hydrogen (Na+/H+) antiporter genes have been shown to play a central role in plant salt tolerance. Here, we analyzed DNA nucleotide polymorphisms in the Na+/H+ antiporter (NHX and SOS1) gene families across 30 different Populus species using several methods of phylogenetic analysis and functional verification. NHX and SOS1 gene families in the genus Populus have expanded from the state in their common ancestors by duplication events, and their distinct lineages have been retained. Signals of positive selection at certain amino acid sites in different members of the Na/H antiporter gene families show that the dynamics that drive the evolution of each gene vary. SOS1 has undergone duplication in Populus euphratica and been subjected to adaptive evolution in section Turanga; the paralog of PeSOS1 (PeSOS1.2) can complement the Escherichia coli mutant EP432; and the expression pattern of PeSOS1.2 is different from that of PeSOS1, a fact which may have been beneficial for P. euphratica, conferring a fitness advantage in saline habitats. The divergent evolution of the individual members of the NHX and SOS1 gene families is likely to have been influenced by the varied ecological and environmental niches occupied by the different poplar species, giving rise to evolutionary footprints that reflect the specific functions and subcellular localizations of the proteins encoded by these genes.

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.

Identification of processed Chinese medicinal materials using DNA mini-barcoding

Most of Chinese medicinal herbs are subjected to traditional processing procedures, including stir-frying, charring, steaming, boiling, and calcining before they are released into dispensaries. The marketing and identification of processed medicinal materials is a growing issue in the marketplace. However, conventional methods of identification have limitations, while DNA mini-barcoding, based on the sequencing of a short-standardized region, has received considerable attention as a new potential means to identify processed medicinal materials. In the present study, six DNA barcode loci including ITS2, psbA-trnH, rbcL, matK, trnL (UAA) intron and its P6 loop, were employed for the authentication of 45 processed samples belonging to 15 species. We evaluated the amplification efficiency of each locus. We also examined the identification accuracy of the potential mini-barcode locus, of trnL (UAA) intron P6 loop. Our results showed that the five primary barcode loci were successfully amplified in only 8.89%-20% of the processed samples, while the amplification rates of the trnL (UAA) intron P6 loop were higher, at 75.56% successful amplification. We compared the mini-barcode sequences with Genbank using the Blast program. The analysis showed that 45.23% samples could be identified to genus level, while only one sample could be identified to the species level. We conclude that trnL (UAA) p6 loop is a candidate mini-barcode that has shown its potential and may become a universal mini-barcode as complementary barcode for authenticity testing and will play an important role in medicinal materials control.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

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

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

De novo transcriptome assembly, development of EST-SSR markers and population genetic analyses for the desert biomass willow, Salix psammophila

Salix psammophila, a sandy shrub known as desert willow, is regarded as a potential biomass feedstock and plays an important role in maintaining local ecosystems. However, a lack of genomic data and efficient molecular markers limit the study of its population evolution and genetic breeding. In this study, chromosome counts, flow cytometry and SSR analyses indicated that S. psammophila is tetraploid. A total of 6,346 EST-SSRs were detected based on 71,458 de novo assembled unigenes from transcriptome data. Twenty-seven EST-SSR markers were developed to evaluate the genetic diversity and population structure of S. psammophila from eight natural populations in Northern China. High levels of genetic diversity (mean 10.63 alleles per locus; mean HE 0.689) were dectected in S. psammophila. The weak population structure and little genetic differentiation (pairwise FST = 0.006-0.016) were found among Population 1-Population 7 (Pop1-Pop7; Inner Mongolia and Shaanxi), but Pop8 (Ningxia) was clearly separated from Pop1-Pop7 and moderate differentiation (pairwise FST = 0.045-0.055) was detected between them, which may be influenced by local habitat conditions. Molecular variance analyses indicated that most of the genetic variation (94.27%) existed within populations. These results provide valuable genetic informations for natural resource conservation and breeding programme optimisation of S. psammophila.

Distinct patterns of natural selection in Na+/H+ antiporter genes in Populus euphratica and Populus pruinosa

Salt tolerance genes constitute an important class of loci in plant genomes. Little is known about the extent to which natural selection in saline environments has acted upon these loci, and what types of nucleotide diversity such selection has given rise to. Here, we surveyed genetic diversity in three types of Na⁺/H⁺ antiporter gene (SOS, NhaD, and NHX, belonging to the cation/proton antiporter 1 family), which have well‐characterized essential roles in plant salt tolerance. Ten Na⁺/H⁺ antiporter genes and 16 neutral loci randomly selected as controls were sequenced from 17 accessions of two closely related members of the genus Populus, Populus euphratica and Populus pruinosa, section Turanga, which are native to northwest China. The results show that salt tolerance genes are common targets of natural selection in P. euphratica and P. pruinosa. Moreover, the patterns of nucleotide variation across the three types of Na⁺/H⁺ antiporter gene are distinctly different in these two closely related Populus species, and gene flow from P. pruinosa to P. euphratica is highly restricted. Our results suggest that natural selection played an important role in shaping the current distinct patterns of Na⁺/H⁺ antiporter genes, resulting in adaptive evolution in P. euphratica and P. pruinosa.

Phylogeography of East Asian Lespedeza buergeri (Fabaceae) based on chloroplast and nuclear ribosomal DNA sequence variations

The dynamic changes in land configuration during the Quaternary that were accompanied by climatic oscillations have significantly influenced the current distribution and genetic structure of warm-temperate forests in East Asia. Although recent surveys have been conducted, the historical migration of forest species via land bridges and, especially, the origins of Korean populations remains conjectural. Here, we reveal the genetic structure of Lespedeza buergeri, a warm-temperate shrub that is disjunctively distributed around the East China Sea (ECS) at China, Korea, and Japan. Two non-coding regions (rpl32-trnL, psbA-trnH) of chloroplast DNA (cpDNA) and the internal transcribed spacer of nuclear ribosomal DNA (nrITS) were analyzed for 188 individuals from 16 populations, which covered almost all of its distribution. The nrITS data demonstrated a genetic structure that followed geographic boundaries. This examination utilized AMOVA, comparisons of genetic differentiation based on haplotype frequency/genetic mutations among haplotypes, and Mantel tests. However, the cpDNA data showed contrasting genetic pattern, implying that this difference was due to a slower mutation rate in cpDNA than in nrITS. These results indicated frequent migration by this species via an ECS land bridge during the early Pleistocene that then tapered gradually toward the late Pleistocene. A genetic isolation between western and eastern Japan coincided with broad consensus that was suggested by the presence of other warm-temperate plants in that country. For Korean populations, high genetic diversity indicated the existence of refugia during the Last Glacial Maximum on the Korean Peninsula. However, their closeness with western Japanese populations at the level of haplotype clade implied that gene flow from western Japanese refugia was possible until post-glacial processing occurred through the Korea/Tsushima Strait land bridge.

Application of Microsatellite Loci for Molecular Identification of Elite Genotypes, Analysis of Clonality, and Genetic Diversity in Aspen Populus tremula L. (Salicaceae)

Testing systems for molecular identification of micropropagated elite aspen (Populus tremula L.) genotypes were developed on the base on microsatellite (SSR) loci. Out of 33 tested microsatellite loci, 14 were selected due to sustainable PCR amplification and substantial variability in elite clones of aspen aimed for establishment of fast-rotated forest plantations. All eight tested clones had different multilocus genotypes. Among 114 trees from three reference native stands located near the established plantations, 80 haplotypes were identified while some repeated genotypes were attributed to natural clones which appeared as a result of sprouting. The selected set of SSR markers showed reliable individual identification with low probability of appearance of identical aspen genotypes (a minimum of 4.8 · 10(-10) and 1 × 10(-4) for unrelated and related individuals, resp.). Case studies demonstrating practical applications of the test system are described including analysis of clonal structure and levels of genetic diversity in three natural aspen stands growing in the regions where plantations made of elite clones were established.

Response to drought and salt stress in leaves of poplar (Populus alba × Populus glandulosa): expression profiling by oligonucleotide microarray analysis

Drought and salt stresses are major environmental constraints on forest productivity. To identify genes responsible for stress tolerance, we conducted a genome-wide analysis in poplar (Populus alba × Populus glandulosa) leaves exposed to drought and salt (NaCl) stresses. We investigated gene expression at the mRNA level using oligonucleotide microarrays containing 44,718 genes from Populus trichocarpa. A total of 1604 and 1042 genes were up-regulated (≥2-fold; P value < 0.05) by drought and salt stresses, respectively, and 765 genes were up-regulated by both stresses. In addition, 2742 and 1685 genes were down-regulated by drought and salt stresses, respectively, and 1564 genes were down-regulated by both stresses. The large number of genes regulated by both stresses suggests that crosstalk occurs between the drought and salt stress responses. Most up-regulated genes were involved in functions such as subcellular localization, signal transduction, metabolism, and transcription. Among the up-regulated genes, we identified 47 signaling proteins, 65 transcription factors, and 43 abiotic stress-related genes. Several genes were modulated by only one of the two stresses. About 25% of the genes significantly regulated by these stresses are of unknown function, suggesting that poplar may provide an opportunity to discover novel stress-related genes.

Transcriptome dynamics of a desert poplar (Populus pruinosa) in response to continuous salinity stress

Using RNA sequencing analysis, we identified 9,216 regulatory and salt-related genes with differential expression and temporal expression trends which provide a clear picture of transcriptomic dynamics in response to continuous salinity stress in a desert poplar, Populus pruinosa. Populus pruinosa Schrenk is native to the desert region of western China and extraordinarily well adapted to the local salt stress. Thus, it is an ideal model for studying plants' adaptation to salt stress, but its transcriptomic responses have not been previously characterized. Thus, we analyzed time- courses of these responses via a series of sequencings. In total, we generated 157.4 million 100 bp paired-end clean reads and identified 9,216 differentially expressed genes (DEGs) between salt-stressed calli and controls. Gene ontology classification analysis revealed that salt stress-related categories--including 'oxidation reduction', 'transcription factor activity', 'membrane' and 'ion channel activity'--were highly enriched among these DEGs. In addition, we grouped the 9,216 DEGs by their expression dynamics into four clusters, and the genes in each cluster showed enrichment for particular functional categories. We also found that most DEGs were activated within 24 h of the stress and their expression stabilized after 48 h. All these findings suggest that gene expression rapidly and coordinately changes during this species' adaptation to salt stress. In addition, the identified DEGs provide critical genetic resources for further functional analyses and indications of potential transgenic modifications for developing salt-tolerant poplars.

Transcriptome differences between two sister desert poplar species under salt stress

BACKGROUND

Populus euphratica Oliv and P. pruinosa Schrenk (Salicaceae) both grow in dry desert areas with high summer temperatures. However, P. euphratica is distributed in dry deserts with deep underground water whereas P. pruinosa occurs in deserts in which there is underground water close to the surface. We therefore hypothesized that these two sister species may have evolved divergent regulatory and metabolic pathways during their interaction with different salt habitats and other stresses. To test this hypothesis, we compared transcriptomes from callus exposed to 24 h of salt stress and control callus samples from both species and identified differentially expressed genes (DEGs) and alternative splicing (AS) events that had occurred under salt stress.

RESULTS

A total of 36,144 transcripts were identified and 1430 genes were found to be differentially expressed in at least one species in response to salt stress. Of these DEGs, 884 and 860 were identified in P. euphratica and P. pruinosa, respectively, while 314 DEGs were common to both species. On the basis of parametric analysis of gene set enrichment, GO enrichment in P. euphratica was found to be significantly different from that in P. pruinosa. Numerous genes involved in hormone biosynthesis, transporters and transcription factors showed clear differences between the two species in response to salt stress. We also identified 26,560 AS events which were mapped to 8380 poplar genomic loci from four libraries. GO enrichments for genes undergoing AS events in P. euphratica differed significantly from those in P. pruinosa.

CONCLUSIONS

A number of salt-responsive genes in both P. euphratica and P. pruinosa were identified and candidate genes with potential roles in the salinity adaptation were proposed. Transcriptome comparisons of two sister desert poplar species under salt stress suggest that these two species may have developed different genetic pathways in order to adapt to different desert salt habitats. The DEGs that were found to be common to both species under salt stress may be especially important for future genetic improvement of cultivated poplars or other crops through transgenic approaches in order to increase tolerance of saline soil conditions.

Speciation of two desert poplar species triggered by Pleistocene climatic oscillations

Despite the evidence that the Pleistocene climatic fluctuations have seriously affected the distribution of intraspecific diversity, less is known on its impact on interspecific divergence. In this study, we aimed to test the hypothesis that the divergence of two desert poplar species Populus euphratica Oliv. and P. pruinosa Schrenk. occurred during the Pleistocene. We sequenced 11 nuclear loci in 60 individuals from the two species to estimate the divergence time between them and to test whether gene flow occurred after species separation. Divergence time between the two species was estimated to be 0.66-1.37 million years ago (Ma), a time at which glaciation was at its maximum in China and deserts developed widely in central Asia. Isolation-with-Migration model also indicated that the two species had diverged in the presence of gene flow. We also detected evidence of selection at GO in P. euphratica and to a lesser extent at PhyB2. Together, these results underscore the importance of Pleistocene climate oscillations in triggering plant speciation as a result of habitats divergence.

Barcoding poplars (Populus L.) from western China

Background

Populus is an ecologically and economically important genus of trees, but distinguishing between wild species is relatively difficult due to extensive interspecific hybridization and introgression, and the high level of intraspecific morphological variation. The DNA barcoding approach is a potential solution to this problem.

Methodology/Principal Findings

Here, we tested the discrimination power of five chloroplast barcodes and one nuclear barcode (ITS) among 95 trees that represent 21 Populus species from western China. Among all single barcode candidates, the discrimination power is highest for the nuclear ITS, progressively lower for chloroplast barcodes matK (M), trnG-psbK (G) and psbK-psbI (P), and trnH-psbA (H) and rbcL (R); the discrimination efficiency of the nuclear ITS (I) is also higher than any two-, three-, or even the five-locus combination of chloroplast barcodes. Among the five combinations of a single chloroplast barcode plus the nuclear ITS, H+I and P+I differentiated the highest and lowest portion of species, respectively. The highest discrimination rate for the barcodes or barcode combinations examined here is 55.0% (H+I), and usually discrimination failures occurred among species from sympatric or parapatric areas.

Conclusions/Significance

In this case study, we showed that when discriminating Populus species from western China, the nuclear ITS region represents a more promising barcode than any maternally inherited chloroplast region or combination of chloroplast regions. Meanwhile, combining the ITS region with chloroplast regions may improve the barcoding success rate and assist in detecting recent interspecific hybridizations. Failure to discriminate among several groups of Populus species from sympatric or parapatric areas may have been the result of incomplete lineage sorting, frequent interspecific hybridizations and introgressions. We agree with a previous proposal for constructing a tiered barcoding system in plants, especially for taxonomic groups that have complex evolutionary histories (e.g. Populus).

Rapidly evolving genes and stress adaptation of two desert poplars, Populus euphratica and P. pruinosa

Understanding which genes have evolved rapidly with the recent tree speciation in arid habitats can provide valuable insights into different adaptation mechanisms. We employed a comparative evolutionary analysis of expressed sequence tags (ESTs) from two desert poplars, Populus pruinosa and P. euphratica, which diverged in the recent past. Following an approach taken previously with P. euphratica, we conducted a deep transcriptomic analysis of P. pruinosa. To maximize representation of conditional transcripts, mRNA was obtained from living tissues of two types of callus and desert-grown trees. De novo assembly generated 114,866 high-quality unique sequences using Solexa sequence data. Following assembly we were able to identify, with high confidence, 2859 orthologous sequence pairs between the two species. Based on the ratio of nonsynonymous (Ka) to synonymous (Ks) substitutions, we identified a total of 84 (2.9%) ortholog pairs exhibiting rapid evolution with signs of strong selection (Ka/Ks>1). Genes homologous to these ortholog pairs in model species are mainly involved in ‘responses to stress’, ‘ubiquitin-dependent protein catabolic processes’, and ‘biological regulation’. Finally, we examined the expression patterns of candidate genes with rapid evolution in response to salt stress. Only one pair of orthologs up-regulated their expression in both species while three and four genes were found to up-regulated in P. pruinosa and in P. euphratica respectively. Our findings together suggest that the genes at the same category or network but with differentiated expressions or functions may have evolved rapidly during adaptive divergence of the two species to differentiated salty desert habitats.

Exploiting the transcriptome of Euphrates Poplar, Populus euphratica (Salicaceae) to develop and characterize new EST-SSR markers and construct an EST-SSR database

BACKGROUND

Microsatellite markers or Simple Sequence Repeats (SSRs) are the most popular markers in population/conservation genetics. However, the development of novel microsatellite markers has been impeded by high costs, a lack of available sequence data and technical difficulties. New species-specific microsatellite markers were required to investigate the evolutionary history of the Euphratica tree, Populus euphratica, the only tree species found in the desert regions of Western China and adjacent Central Asian countries.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 94,090 non-redundant Expressed Sequence Tags (ESTs) from P. euphratica comprising around 63 Mb of sequence data were searched for SSRs. 4,202 SSRs were found in 3,839 ESTs, with 311 ESTs containing multiple SSRs. The most common motif types were trinucleotides (37%) and hexanucleotides (33%) repeats. We developed primer pairs for all of the identified EST-SSRs (eSSRs) and selected 673 of these pairs at random for further validation. 575 pairs (85%) gave successful amplification, of which, 464 (80.7%) were polymorphic in six to 24 individuals from natural populations across Northern China. We also tested the transferability of the polymorphic eSSRs to nine other Populus species. In addition, to facilitate the use of these new eSSR markers by other researchers, we mapped them onto Populus trichocarpa scaffolds in silico and compiled our data into a web-based database (http://202.205.131.253:8080/poplar/resources/static_page/index.html).

CONCLUSIONS

The large set of validated eSSRs identified in this work will have many potential applications in studies on P. euphratica and other poplar species, in fields such as population genetics, comparative genomics, linkage mapping, QTL, and marker-assisted breeding. Their use will be facilitated by their incorporation into a user-friendly web-based database.

Two highly validated SSR multiplexes (8-plex) for Euphrates' poplar, Populus euphratica (Salicaceae)

Multiplex PCR amplification of microsatellites has significantly increased the throughput and decreased the costs of genotyping. We have developed two highly polymorphic microsatellite multiplexes for Populus euphratica, the only tree species found in desert regions of Western China and adjacent Central Asian countries. The first of these multiplex kits comprises an eight-Plex of genomic SSRs (gSSRs) obtained from published databases. The second comprises an eight-plex of newly designed EST-SSRs (eSSRs) based on expressed sequence tags for P. euphratica. Both kits were tested on a sample of 170 individuals from four populations. The gSSRs exhibited slightly more polymorphism than the eSSRs. The new multiplex protocols yielded consistent results in the hands of multiple researchers, demonstrating their robustness. The 16 loci used in the kits exhibited a high transferability rate (82.0%) in eight other poplar species belonging to five different sections of the genus. Both kits should therefore be useful for further investigations of population genetics in P. euphratica and related species. Our results indicate that it is essential to follow recently established recommendations when developing microsatellite markers, including verifying the amplification efficiency, detecting null alleles and carefully measuring error rates.

Plant DNA barcodes and the influence of gene flow

Success of species assignment using DNA barcodes has been shown to vary among plant lineages because of a wide range of different factors. In this study, we confirm the theoretical prediction that gene flow influences species assignment with simulations and a literature survey. We show that the genome experiencing the highest gene flow is, in the majority of the cases, the best suited for species delimitation. Our results clearly suggest that, for most angiosperm groups, plastid markers will not be the most appropriate for use as DNA barcodes. We therefore advocate shifting the focus from plastid to nuclear markers to achieve an overall higher success using DNA barcodes.