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 FindingsHere, 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/SignificanceIn 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).

DNA sequence variation has been widely reported as the genetic basis for adaptation, in both humans and other animals, to the hypoxic environment experienced at high altitudes. However, little is known about the patterns of gene expression underlying such hypoxic adaptations. In this study, we examined the differences in the transcriptomes of four organs (heart, kidney, liver and lung) between yak and cattle, a pair of closely related species distributed at high and low altitudes respectively. Of the four organs examined, heart shows the greatest differentiation between the two species in terms of gene expression profiles. Detailed analyses demonstrated that some genes associated with the oxygen supply system and the defense systems that respond to threats of hypoxia are differentially expressed. In addition, genes with significantly differentiated patterns of expression in all organs exhibited an unexpected uniformity of regulation along with an elevated frequency of nonsynonymous substitutions. This co-evolution of protein sequences and gene expression patterns is likely to be correlated with the optimization of the yak metabolic system to resist hypoxia.

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'.

Investigating genome-wide variation patterns along a speciation continuum is of central importance to understand the evolutionary processes contributing to lineage diversification. To identify which forces have shaped the genomic landscapes in Populus, we resequenced 201 whole-genomes from eight closely related species, with pairs of species at different stages along the speciation continuum. Using population structure and identity by descent analyses, we show extensive introgression between some species pairs, especially those with parapatric distributions. Inference of historical changes in effective population sizes support species-specific demographic trajectories, including recent population expansions in species characterized by broad present-day distributions. We observe highly conserved genomic landscapes, either focusing on within-species (genetic diversity: π and recombination rate: ρ) or among-species variation (relative divergence: FST and absolute divergence: DXY ). Independent of the stage across the divergence continuum, we recovered positive correlations between the pair π and ρ and the pair DXY and ρ across all species pairs, which is consistent with a substantial contribution of linked selection in shaping these genomic landscapes. However, the positive correlations between π and DXY became weaker as the overall divergence level (da) increased, suggesting that background selection is not the only factor at play. Positive correlations between FST and DXY in all species pairs, regardless of the rate of gene flow, also indicates the high FST could be due to divergent sorting of ancient polymorphism before speciation. Our study showcases the importance of investigating genomic patterns on multiple species across the speciation continuum to better understand the genomic landscapes of diversity and differentiation.

All Qinghai-Tibetan Plateau (QTP) endemic species are assumed to have originated recently, although very rare species most likely diverged early. These ancient species provide an excellent model to examine the origin and evolution of QTP endemic plants in response to the QTP uplifts and the climate changes that followed in this high altitude region. In this study, we examined these hypotheses by employing sequence variation from multiple nuclear and chloroplast DNA of 239 individuals of Juniperus microsperma and its five congeners. Both phylogenetic and population genetic analyses revealed that J. microsperma diverged from its sister clade comprising two species with long isolation around the Early Miocene, which corresponds to early QTP uplift. Demographic modeling and coalescent tests suggest that J. microsperma experienced an obvious bottleneck event during the Quaternary when the global climate greatly oscillated. The results presented here support the hypotheses that the QTP uplifts and Quaternary climate changes played important roles in shaping the evolutionary history of this rare juniper.

No abstract

Leyland cypress [×Hesperotropsis leylandii (A.B. Jacks. & Dallim.) Garland & G. Moore, Cupressaceae] is a well-known horticultural evergreen conifer in the United Kingdom, United States, Australia, New Zealand, and other countries. As demonstrated by previous studies, this taxon is a hybrid between alaska (nootka) cypress [Callitropsis nootkatensis (D. Don) Oerst. ex D.P. Little] and monterey cypress [Hesperocyparis macrocarpa (Hartw. ex Gordon) Bartel]. However, the genetic background of leyland cypress cultivars is unclear. Are they F1 or F2 hybrids or backcrosses? In this study, six individuals that represent major leyland cypress cultivars and two individuals each of its two putative parental species were collected, and three nuclear DNA regions (internal transcribed spacer, leafy and needly), three mitochondrial (mt) DNA regions (coxI, atpA, and rps3), and two chloroplast (cp) DNA regions (matK and rbcL) were sequenced and analyzed. Sequencing results of nuclear DNA regions revealed that leyland cypress cultivars consist of putative F1 and F2 hybrids as well as backcrosses. Analysis of the cp and mt DNA from six cultivars of leyland cypress revealed that their cytoplasmic (cp and mt) genomes came from alaska cypress. Our findings will provide important instructions and background knowledge on the management of these major leyland cypress cultivars as well as future studies. Meanwhile, alaska cypress and monterey cypress may have diverged with each other ≈46 million years ago. The fact that they can produce fertile hybrids indicates that hybridization events may have played an important role in the evolutionary history of the cypress family (Cupressaceae).

No abstract