Species delimitation and biogeography of two fir species (Abies) in central China: cytoplasmic DNA variation

Distributional Response of the Rare and Endangered Tree Species Abies chensiensis to Climate Change in East Asia

Globally, increasing temperatures due to climate change have severely affected natural ecosystems in several regions of the world; however, the impact on the alpine plant may be particularly profound, further raising the risk of extinction for rare and endangered alpine plants. To identify how alpine species have responded to past climate change and to predict the potential geographic distribution of species under future climate change, we investigated the distribution records of A. chensiensis, an endangered alpine plant in the Qinling Mountains listed in the Red List. In this study, the optimized MaxEnt model was used to analyse the key environmental variables related to the distribution of A. chensiensis based on 93 wild distribution records and six environmental variables. The potential distribution areas of A. chensiensis in the last interglacial (LIG), the last glacial maximum (LGM), the current period, and the 2050s and 2070s were simulated. Our results showed that temperature is critical to the distribution of A. chensiensis, with the mean temperature of the coldest quarter being the most important climatic factor affecting the distribution of this species. In addition, ecological niche modeling analysis showed that the A. chensiensis distribution area in the last interglacial experiencing population expansion and, during the last glacial maximum occurring, a population contraction. Under the emission scenarios in the 2050s and 2070s, the suitable distribution area would contract significantly, and the migration routes of the centroids tended to migrate toward the southern high-altitude mountains, suggesting a strong response from the A. chensiensis distribution to climate change. Collectively, the results of this study provide a comprehensive and multidimensional perspective on the geographic distribution pattern and history of population dynamics for the endemic, rare, and endangered species, A. chensiensis, and it underscores the significant impact of geological and climatic changes on the geographic pattern of alpine species populations.

Glacial Expansion or Interglacial Expansion? Contrasting Demographic Models of Four Cold-Adapted Fir Species in North America and East Asia

Understanding and forecasting species’ response to climate change is a critical need for future conservation and management. Two expansion hypotheses, the glacial expansion versus the interglacial expansion, have been proposed to interpret how cold-adapted organisms in the northern hemisphere respond to Quaternary climatic fluctuations. To test these two hypotheses, we originally used two pairs of high-low elevation firs from North America (Abies lasiocarpa and Abies balsamea) and East Asia (Abies chensiensis and Abies nephrolepis). Abies lasiocarpa and Abies chensiensis are widely distributed in high-elevation regions of western North America and central China. Abies balsamea and Abies nephrolepis occur in central North America and northeast China, with much lower elevations. These fir species are typical cold-adapted species and sensitive to climate fluctuations. Here, we integrated the mtDNA and cpDNA polymorphisms involving 44 populations and 585 individuals. Based on phylogeographic analyses, recent historical range expansions were indicated in two high-elevation firs (Abies lasiocarpa and Abies chensiensis) during the last glaciation (43.8–28.4 or 21.9–14.2 kya, 53.1–34.5 or 26.6–17.2 kya). Such glacial expansions in high-elevation firs were further confirmed by the evidence of species distribution modelling, geographic-driven genetic patterns, palynological records, and current distribution patterns. Unlike the north American firs, the SDM models indicated unremarkable expansion or contraction in East Asia firs for its much more stable conditions during different historical periods. Taken together, our findings highly supported that high-elevation firs experienced glacial expansion during the Quaternary climate change in East Asia and North America, as interglacial expansion within low-elevation firs. Under this situation, the critically endangered fir species distributed in high elevation would have no enough higher elevational space to migrate. Facing the increasing global warming, thus we proposed ex-situ conservation of defining conservation units as the most meaningful strategy.

Phylogeny and Species Delimitation of Chinese Medicago (Leguminosae) and Its Relatives Based on Molecular and Morphological Evidence

Medicago and its relatives, Trigonella and Melilotus comprise the most important forage resources globally. The alfalfa selected from the wild relatives has been cultivated worldwide as the forage queen. In the Flora of China, 15 Medicago, eight Trigonella, and four Melilotus species are recorded, of which six Medicago and two Trigonella species are introduced. Although several studies have been conducted to investigate the phylogenetic relationship within the three genera, many Chinese naturally distributed or endemic species are not included in those studies. Therefore, the taxonomic identity and phylogenetic relationship of these species remains unclear. In this study, we collected samples representing 18 out of 19 Chinese naturally distributed species of these three genera and three introduced Medicago species, and applied an integrative approach by combining evidences from population-based morphological clusters and molecular data to investigate species boundaries. A total of 186 individuals selected from 156 populations and 454 individuals from 124 populations were collected for genetic and morphological analyses, respectively. We sequenced three commonly used DNA barcodes (trnH-psbA, trnK-matK, and ITS) and one nuclear marker (GA3ox1) for phylogenetic analyses. We found that 16 out of 21 species could be well delimited based on phylogenetic analyses and morphological clusters. Two Trigonella species may be merged as one species or treated as two subspecies, and Medicago falcata should be treated as a subspecies of the M. sativa complex. We further found that major incongruences between the chloroplast and nuclear trees mainly occurred among the deep diverging lineages, which may be resulted from hybridization, incomplete lineage sorting and/or sampling errors. Further studies involving a finer sampling of species associated with large scale genomic data should be employed to better understand the species delimitation of these three genera.

Elevation Shift in Abies Mill. (Pinaceae) of Subtropical and Temperate China and Vietnam-Corroborative Evidence from Cytoplasmic DNA and Ecological Niche Modeling

The "elevational shift" scenario has been proposed as a model to explain the response of cold-adapted organisms to Quaternary climatic fluctuations in Europe and North America. However, the elevational shift model has not been well-explored in eastern Asia, which is more topographically complex than the other Northern Hemisphere biogeographic regions. Here, we evaluated the role of elevational shift in the closely related firs, or Abies Mill., of subtropical and temperate China. These firs are typical alpine trees with sensitivity to climate change. We tested the elevational shift hypothesis in firs of China using phylogeographic methods and ecological niche models. Our phylogeographic analyses comprised mitochondrial and chloroplast polymorphisms surveyed across 479 individuals from 43 populations representing 11 species. M1 of the 11 mitotypes and C1 of the 25 chlorotypes were inferred as the ancestral haplotype, and they had the widest distribution. The results of our phylogeographic survey revealed multiple centers of genetic diversity in distinct geographic regions and no latitudinal trend. Moreover, our results showed range expansions for seven taxa during the last glacial (64.9-18.2 or 32.5-9.1 kya), and this was consistent with the Quaternary fossil record of Abies in China. Taken together, our findings support a historical biogeographic pattern in firs of glacial expansions, probably through corridors at lower elevation, and interglacial fragmentations, through isolation at higher elevation peaks. Therefore, Abies in China probably undergoes elevational shift in response to climate change. Facing the forecasting global warming, the risk of several critically endangered firs was further enhanced as these species would have little escape space in situ to higher altitudes. According to our ENMs, we proposed an ex situ conservation strategy in the southern Hengduan Mountains region of south western China.

Molecular phylogenetic relationships reveal contrasting evolutionary patterns in Gorgoniidae (Octocorallia) in the Eastern Pacific

The description and delimitation of species in an evolutionary framework is essential for understanding patterns of biodiversity and distribution, and in the assessment of conservation strategies for natural resources. This study seeks to clarify the evolutionary history and genetic variation within and between closely related octocoral species that are fundamental to benthic marine ecosystems for harbouring a high diversity of associated fauna. For our study system, we focused on members of the Gorgoniidae family in the Eastern Pacific, particularly of the Ecuadorian littoral, a less studied marine ecosystem. According to our results, the diagnosis of the genus Pacifigorgia is here amended to include species previously considered in the genus Leptogorgia. The genera Leptogorgia and Eugorgia are included within a single clade, and neither are recovered as monophyletic. In this case, according to the priority rule of the International Code of Zoological Nomenclature (ICZN), our proposal is to include the species considered in these two genera in Leptogorgia. In addition, we found evidence of interesting speciation patterns: morphological differentiation with no apparent genetic differentiation (in Pacifigorgia), and inconsistencies between mitochondrial and nuclear data that suggest a hybridisation phenomenon (in Leptogorgia). In the first case, recent radiation, ancient hybridisation, sympatric speciation, and in the second, reticulate evolution may have contributed to the evolutionary history of the studied taxa. Therefore, incongruences observed between morphological and molecular evidences in these octocorals, and in corals in general, may reveal the types of events/patterns that have influenced their evolution.

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.

Importance of incomplete lineage sorting and introgression in the origin of shared genetic variation between two closely related pines with overlapping distributions

Genetic variation shared between closely related species may be due to retention of ancestral polymorphisms because of incomplete lineage sorting (ILS) and/or introgression following secondary contact. It is challenging to distinguish ILS and introgression because they generate similar patterns of shared genetic diversity, but this is nonetheless essential for inferring accurately the history of species with overlapping distributions. To address this issue, we sequenced 33 independent intron loci across the genome of two closely related pine species (Pinus massoniana Lamb. and Pinus hwangshanensis Hisa) from Southeast China. Population structure analyses revealed that the species showed slightly more admixture in parapatric populations than in allopatric populations. Levels of interspecific differentiation were lower in parapatry than in allopatry. Approximate Bayesian computation suggested that the most likely speciation scenario explaining this pattern was a long period of isolation followed by a secondary contact. Ecological niche modeling suggested that a gradual range expansion of P. hwangshanensis during the Pleistocene climatic oscillations could have been the cause of the overlap. Our study therefore suggests that secondary introgression, rather than ILS, explains most of the shared nuclear genomic variation between these two species and demonstrates the complementarity of population genetics and ecological niche modeling in understanding gene flow history. Finally, we discuss the importance of contrasting results from markers with different dynamics of migration, namely nuclear, chloroplast and mitochondrial DNA.

Species delimitation of Chinese hop-hornbeams based on molecular and morphological evidence

Species delimitation through which infers species boundaries is emerging as a major work in modern systematics. Hop-hornbeam species in Ostrya (Betulaceae) are well known for their hard and heavy woods. Five species were described in China and their interspecific delimitations remain unclear. In this study, we firstly explored their distributions in all recorded field sites distributed in China. We then selected 110 samples from 22 natural populations of five species from this genus and one type specimen of O. yunnanensis, for molecular barcoding analyses. We sequenced four chloroplast (cp) DNA fragments (trnH-psbA, trnL-trnF, rps16, and trnG) and the nuclear internal transcribed spacer (ITS) region for all samples. Sequence variations of Ostrya from four cpDNA fragments identified three groups that showed no correspondence to any morphological delimitation because of the incomplete lineage sorting and/or possible interspecific introgression in the history. However, phylogenetic analyses of ITS sequence variations discerned four species, O. japonica, O. rehderiana, O. trichocarpa, and O. multinervis while O. yunnanensis nested within O. multinervis. Morphological clustering also discerned four species and showed the complete consistency with molecular evidence. Moreover, our phylogenetic analyses-based ITS sequence variations suggested that O. trichocarpa comprised an isolated lineage different from the other Eurasian ones. Based on these results, hop-hornbeams in China should be treated as four separate species. Our results further highlight the importance of ITS sequence variations in delimitating and discerning the closely related species in plants.

The extent and meaning of hybridization and introgression between Siberian spruce (Picea obovata) and Norway spruce (Picea abies): cryptic refugia as stepping stones to the west?

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well-separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.

Genetic divergence and phylogeographic history of two closely related species (Leucomeris decora and Nouelia insignis) across the 'Tanaka Line' in Southwest China

BACKGROUND

Leucomeris decora and Nouelia insignis (Asteraceae) are narrowly and allopatrically distributed species, separated by the important biogeographic boundary Tanaka Line in Southwest China. Previous morphological, cytogenetic and molecular studies suggested that L. decora is sister to N. insignis. However, it is less clear how the two species diverged, whether in full isolation or occurring gene flow across the Tanaka Line. Here, we performed a molecular study at the population level to characterize genetic differentiation and decipher phylogeographic history in two closely related species based on variation examined in plastid and nuclear DNAs using a coalescent-based approach.

RESULTS

These morphologically distinct species share plastid DNA (cpDNA) haplotypes. In contrast, Bayesian analysis of nuclear DNA (nDNA) uncovered two distinct clusters corresponding to L. decora and N. insignis. Based on the IMa analysis, no strong indication of migration was detected based on both cpDNA and nDNA sequences. The molecular data pointed to a major west-east split in nuclear DNA between the two species corresponding with the Tanaka Line. The coalescent time estimate for all cpDNA haplotypes dated to the Mid-Late Pleistocene. The estimated demographic parameters showed that the population size of L. decora was similar to that of N. insignis and both experienced limited demographic fluctuations recently.

CONCLUSIONS

The study revealed comprehensive species divergence and phylogeographic histories of N. insignis and L. decora divided by the Tanaka Line. The phylogeographic pattern inferred from cpDNA reflected ancestrally shared polymorphisms without post-divergence gene flow between species. The marked genealogical lineage divergence in nDNA provided some indication of Tanaka Line for its role as a barrier to plant dispersal, and lent support to its importance in promoting strong population structure and allopatric divergence.

Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.

Spatial and temporal population genetic variation and structure of Nothotsuga longibracteata (Pinaceae), a relic conifer species endemic to subtropical China

Nothotsuga longibracteata, a relic and endangered conifer species endemic to subtropical China, was studied for examining the spatial-temporal population genetic variation and structure to understand the historical biogeographical processes underlying the present geographical distribution. Ten populations were sampled over the entire natural range of the species for spatial analysis, while three key populations with large population sizes and varied age structure were selected for temporal analyses using both nuclear microsatellites (nSSR) and chloroplast microsatellites (cpSSR). A recent bottleneck was detected in the natural populations of N. longibracteata. The spatial genetic analysis showed significant population genetic differentiation across its total geographical range. Notwithstanding, the temporal genetic analysis revealed that the level of genetic diversity between different age class subpopulations remained constant over time. Eleven refugia of the Last Glacial Maximum were identified, which deserve particular attention for conservation management.

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

BACKGROUND

The fixed genetic differences between ecologically divergent species were found to change greatly depending on the markers examined. With such species it is difficult to differentiate between shared ancestral polymorphisms and past introgressions between the diverging species. In order to disentangle these possibilities and provide a further case for DNA barcoding of plants, we examine genetic differentiation between two ecologically divergent poplar species, Populus euphratica Oliver and P. pruinosa Schrenk using three different types of genetic marker.

METHODOLOGY/PRINCIPAL FINDINGS

We genotyped 290 individuals from 29 allopatric and sympatric populations, using chloroplast (cp) DNA, nuclear (nr) ITS sequences and eight simple sequence repeat (SSR) loci. Three major cpDNA haplotypes were widely shared between the two species and between-species cpDNA differentiation (F(CT)) was very low, even lower than among single species populations. The average SSR F(CT) values were higher. Bayesian clustering analysis of all loci allowed a clear delineation of the two species. Gene flow, determined by examining all SSR loci, was obvious but only slightly asymmetrical. However, the two species were almost fixed for two different nrITS genotypes that had the highest F(CT), although a few introgressed individuals were detected both in allopatric and sympatric populations.

CONCLUSIONS

The two species shared numerous ancestral polymorphisms at cpDNA and a few SSR loci. Both ITS and a combination of nuclear SSR data could be used to differentiate between the two species. Introgressions and gene flow were obvious between the two species either during or after their divergence. Our findings underscore the complex genetic differentiations between ecologically diverged species and highlight the importance of nuclear DNA (especially ITS) differentiation for delimiting closely related plant species.