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.

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.

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