No abstract

Hybridization and introgression can play an important role in speciation. Here, we examine their roles in the origin and evolution of Picea purpurea, a diploid spruce species occurring on the Qinghai-Tibet Plateau (QTP). Phylogenetic relationships and ecological differences between this species and its relatives, P. schrenkiana, P. likiangensis and P. wilsonii, are unclear. To clarify them, we surveyed sequence variation within and between them for 11 nuclear loci, three chloroplast (cp) and two mitochondrial (mt) DNA fragments, and examined their ecological requirements using ecological niche modelling. Initial analyses based on 11 nuclear loci rejected a close relationship between P. schrenkiana and P. purpurea. BP&P tests and ecological niche modelling indicated substantial divergence between the remaining three species and supported the species status of P. purpurea, which contained many private alleles as expected for a well-established species. Sequence variation for cpDNA and mtDNA suggested a close relationship between P. purpurea and P. wilsonii, while variation at the nuclear se1364 gene suggested P. purpurea was more closely related to P. likiangensis. Analyses of genetic divergence, Bayesian clustering and model comparison using approximate Bayesian computation (ABC) of nuclear (nr) DNA variation all supported the hypothesis that P. purpurea originated by homoploid hybrid speciation from P. wilsonii and P. likiangensis. The ABC analysis dated the origin of P. purpurea at the Pleistocene, and the estimated hybrid parameter indicated that 69% of its nuclear composition was contributed by P. likiangensis and 31% by P. wilsonii. Our results further suggested that during or immediately following its formation, P. purpurea was subject to organelle DNA introgression from P. wilsonii such that it came to possess both mtDNA and cpDNA of P. wilsonii. The estimated parameters indicated that following its origin, P. purpurea underwent an expansion during/after the largest Pleistocene glaciation recorded for the QTP.

A knowledge of intraspecific divergence and range dynamics of dominant forest trees in response to past geological and climate change is of major importance to an understanding of their recent evolution and demography. Such knowledge is informative of how forests were affected by environmental factors in the past and may provide pointers to their response to future environmental change. However, genetic signatures of such historical events are often weak at individual loci due to large effective population sizes and long generation times of forest trees. This problem can be overcome by analysing genetic variation across multiple loci. We used this approach to examine intraspecific divergence and past range dynamics in the conifer Picea likiangensis, a dominant tree of forests occurring in eastern and southern areas of the Qinghai-Tibet Plateau (QTP). We sequenced 13 nuclear loci, two mitochondrial DNA regions and three plastid (chloroplast) DNA regions in 177 individuals sampled from 22 natural populations of this species, and tested the hypothesis that its evolutionary history was markedly affected by Pliocene QTP uplifts and Quaternary climatic oscillations. Consistent with the taxonomic delimitation of the three morphologically divergent varieties examined, all individuals clustered into three genetic groups with intervariety admixture detected in regions of geographical overlap. Divergence between varieties was estimated to have occurred within the Pliocene and ecological niche modelling based on 20 ecological variables suggested that niche differentiation was high. Furthermore, modelling of population-genetic data indicated that two of the varieties (var. rubescens and var. linzhiensis) expanded their population sizes after the largest Quaternary glaciation in the QTP, while expansion of the third variety (var. likiangensis) began prior to this, probably following the Pliocene QTP uplift. These findings point to the importance of geological and climatic changes during the Pliocene and Pleistocene as causes of intraspecific diversification and range shifts of dominant tree species in the QTP biodiversity hot spot region.

Population genetics data based on multiple nuclear loci provide invaluable information to understand demographic, selective, and divergence histories of the current species. We studied nucleotide variation at 13 nuclear loci in 53 populations distributed among four closely related, but morphologically distinct juniper species of the Qinghai-Tibetan Plateau (QTP). We used a novel approach combining Approximate Bayesian Computation and a recently developed neutrality test based on the maximum frequency of derived mutations to examine the demographic and selective histories of individual species, and isolation-with-migration analyses to study the joint history of the species and detect gene flow between them. We found that (1) the four species, which diverged in response to the extensive QTP uplifts, have different demographic histories; (2) two loci, Pgi and CC0822, depart significantly from neutrality in one species and Pgi, is also marginally significant in another; and (3) shared polymorphisms are common, indicating both incomplete lineage sorting and gene flow after species divergence. In addition, the detected unidirectional gene flow provides indirect support for the theoretical prediction that introgression should mostly take place from local to invading species. Our results, together with previous studies, underscore complex evolutionary histories of plant diversification in the biodiversity-hotspot QTP.

Aim To address whether climate changes have altered the distributional ranges of cold‐preferring conifer species, and whether this has influenced interspecific divergence. Location Qinghai‐Tibet Plateau (QTP), Asia. Methods We obtained DNA sequences at 13 nuclear loci in 229 individuals collected from 32 populations across the range of two cold‐preferring alpine species, Picea likiangensis and Picea wilsonii. Bayesian skyline plots (BSPs) were used to examine historical changes in effective population sizes. Species distributions were reconstructed based on presence‐only data at 134 locations and climate variables from the WorldClim database, and the reconstructed distributions were compared with projected palaeodistributions at the Last Glacial Maximum (LGM, about 20 ka). Four models of interspecific gene flow during the Pleistocene and Holocene were compared using approximate Bayesian computation (ABC). Results The effective population sizes of P. likiangensis and P. wilsonii increased continuously from 5 to 0.2 Ma, but decreased after 0.2 Ma. The distributional ranges during the LGM were wider than at present, confirming the trend of a recent decrease. However, the overlapping area became greater due to their shifted distributions, and these shifts altered the major direction of continuing interspecific gene flow. Main conclusions The early increase in the distributions of these two species might have been triggered by enlargements of the available cold habitats brought about by QTP uplift and Pleistocene glaciations. However, the distributional ranges decreased in response to climatic warming since the end of the largest glaciation. These distributional changes further retarded the process of differentiation by altering the pattern of gene flow. Our results shed light on the evolutionary complexity of alpine species in response to climate changes in this climate‐sensitive biodiversity hotspot.