The Quaternary evolutionary history,potential distribution dynamics,and conservation implications for a Qinghai–Tibet Plateau endemic herbaceous perennial,Anisodus tanguticus (Solanaceae)

Various hypotheses have been proposed about the Quaternary evolutionary history of plant species on the Qinghai–Tibet Plateau (QTP), yet only a handful of studies have considered both population genetics and ecological niche context. In this study, we proposed and compared climate refugia hypotheses based on the phylogeographic pattern of Anisodus tanguticus (three plastid DNA fragments and nuclear internal transcribed spacer regions from 32 populations) and present and past species distribution models (SDMs). We detected six plastid haplotypes in two well‐differentiated lineages. Although all haplotypes could be found in its western (sampling) area, only haplotypes from one lineage occurred in its eastern area. Meanwhile, most genetic variations existed between populations (F_ST = 0.822). The SDMs during the last glacial maximum and last interglacial periods showed range fragmentation in the western area and significant range contraction in the eastern area, respectively, in comparison with current potential distribution. This species may have undergone intraspecific divergence during the early Quaternary, which may have been caused by survival in different refugia during the earliest known glacial in the QTP, rather than geological isolation due to orogenesis events. Subsequently, climate oscillations during the Quaternary resulted in a dynamic distribution range for this species as well as the distribution pattern of its plastid haplotypes and nuclear genotypes. The interglacial periods may have had a greater effect on A. tanguticus than the glacial periods. Most importantly, neither genetic data nor SDM alone can fully reveal the climate refugia history of this species. We also discuss the conservation implications for this important Tibetan folk medicine plant in light of these findings and SDMs under future climate models. Together, our results underline the necessity to combine phylogeographic and SDM approaches in future investigations of the Quaternary evolutionary history of species in topographically complex areas, such as the QTP.     

Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae)

Local adaptation to different environments has the potential to maintain divergence between populations despite recurrent gene flow and is an important driver for generating biological diversity. In this study, we investigate the role of adaptation in the maintenance of two parapatric varieties of a forest tree. We used sequence variation of chloroplastic DNA and restriction site‐associated DNA to investigate the genetic structure of two varieties of Castanopsis carlesii in subtropical China and relate it to climatic variation. We used niche reconstruction methods to investigate niche differentiation between the two varieties and to estimate the past distribution of this species. A deep divergence was observed between the two varieties, but evidence of introgression and genetic admixture was detected in two phenotypically and geographically intermediate populations. Niche reconstruction suggests that the distribution of the two varieties was disjunct during periods of global cooling and that the two varieties occupy significantly different niches. The genetic structure was mainly driven by environmental factors, and 13 outlier loci under divergent selection were correlated with climatic variation. These results suggest that the two varieties evolved in allopatry and came back into secondary contact after the last glacial maximum and that they are an evolutionary example of divergence maintained by climatic selection despite recurrent gene flow.     

Extensive population expansion of Pedicularis longiflora (Orobanchaceae) on the Qinghai‐Tibetan Plateau and its correlation with the Quaternary climate change

The Qinghai‐Tibetan Plateau (QTP) is thought to be more strongly affected by the Quaternary glaciations than most other regions of the same latitude. It would be of great interest to investigate the population genetic structure of organisms distributed on the platform and its correlation with the Quaternary climatic oscillations. Here we used the chloroplast (cp)DNA trnT‐trnF sequence to study genetic variation and phylogeography of Pedicularis longiflora, an alpine herb with extensive distribution on the QTP. Based on a range‐wide sampling comprising 41 populations and 910 individuals, we detected 30 cpDNA haplotypes that were divided into five clades by phylogenetic and network analyses and a strong phylogeographical structure. All haplotypes but one in the three basal clades occur exclusively in the southeast QTP, whereas haplotypes in the young clade V occupy almost the whole species range. In particular, the young haplotype H18 occurs in 420 individuals, even at a frequency of 100% in some QTP platform populations and the Altai population. The haplotype distribution pattern, together with molecular clock estimation and mismatch distribution analysis, suggests that the southeast QTP was either a refuge for P. longiflora during the Quaternary climatic change or is the place of origin of the species. The present wide distribution of the species on the QTP platform has resulted from recent population expansions which could be dated back to 120 000–17 000 years ago, a period mostly before the last glacial maximum. The possible relationships among geographic genetic structure, climatic change and species diversification in Pedicularis are also discussed.     

青藏高原岷县龙胆的遗传分化与种群动态历史分析

青藏高原是全球生物多样性中心之一,是研究物种形成和适应性进化的热点地区。岷县龙胆（Gentiana algida var. purdomii）是青藏高原地区特有植物,具有很高的药用价值。本研究以岷县龙胆种群为对象,基于不同遗传方式的分子标记,通过种群遗传结构、遗传分化、种群动态历史和物种分布模型等分析,探讨其在青藏高原的分化历史。结果表明：岷县龙胆具有高的遗传多样性,叶绿体数据表明种内遗传分化程度很高（F_ST=0.452）,但核基因数据表明种内遗传分化程度很低（F_ST=0.022）,揭示岷县龙胆种内存在强烈的基因流。种群动态历史分析表明岷县龙胆种群大小经历了近期扩张,物种分布模型的结果表明岷县龙胆在末次盛冰期以来分布范围略有扩张。分子钟模型的结果表明岷县龙胆种内的遗传分化主要发生在第四纪冰期以来。以上结果一致表明,第四纪气候波动是岷县龙胆分化的主要外部因素。本研究为龙胆属及青藏高原植物类群的物种分化、适应性研究以及资源开发与保护提供了参考。     

Phylogeography of the genus Dasiphora (Rosaceae) in the Qinghai‐Tibetan Plateau: divergence blurred by expansion

Plateau uprisings and climatic oscillations are considered to have caused extensive allopatric divergences that account for the rich species diversity of the Qinghai‐Tibetan Plateau (QTP). However, secondary contact during range shifts in the Quaternary glacial cycles or inter‐uplift stages may have restored the gene flow between species and so counteracted these divergences, particularly in rapidly‐adapting dominant elements. We tested this hypothesis by determining the phylogeographical history of Dasiphora (Rosaceae), a genus of two species that are widely distributed on the QTP and co‐exist in numerous localities. We sequenced two chloroplast DNA fragments (rbcL, trnT‐L) for 559 individuals from 87 populations. Bayesian methods were used to identify phylogenetic relationships and to estimate divergence times. Demographic histories were inferred using neutrality tests, mismatch distribution analysis, and coalescent simulation. A total of 112 haplotypes that clustered into three major groups were identified. The formation of these groups and their subgroups was dated to between the Pliocene and the late Pleistocene. In addition, we found that some groups underwent multiple extensive expansions. Species‐specific haplotypes were identified for each species, although these haplotypes phylogenetically intermixed. These results suggest that recent plateau uplifts and climatic oscillations might have caused the deep divergences observed within this genus. However, later range expansions probably blurred these divergences and possible species boundaries. Our results shed new light on the complex evolutionary history of the QTP alpine plants. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 777–788.     

Late Pleistocene origin of the entire circumarctic range of the arctic‐alpine plant Kalmia procumbens

The circumarctic ranges of arctic‐alpine plants are thought to have been established in the late Pliocene/early Pleistocene, when the modern arctic tundra was formed in response to climate cooling. Previous findings of range‐wide genetic structure in arctic‐alpine plants have been thought to support this hypothesis, but few studies have explicitly addressed the temporal framework of the genetic structure. Here, we estimated the demographic history of the genetic structure in the circumarctic Kalmia procumbens using sequences of multiple nuclear loci and examined whether its genetic structure reflects prolonged isolation throughout the Pleistocene. Both Bayesian clustering and phylogenetic analyses revealed genetic distinction between alpine and arctic regions, whereas detailed groupings were somewhat discordant between the analyses. By assuming a population grouping based on the phylogenetic analyses, which likely reflects a deeper intraspecific divergence, we conducted model‐based analyses and demonstrated that the intraspecific genetic divergence in K. procumbens likely originated during the last glacial period. Thus, there is no need to postulate range separation throughout the Pleistocene to explain the current genetic structure in this species. This study demonstrates that range‐wide genetic structure in arctic‐alpine plants does not necessarily result from the late Pliocene/early Pleistocene origin of their circumarctic ranges and emphasizes the importance of a temporal framework of the current genetic structure for understanding the biogeographic history of the arctic flora.     

Intraspecific differentiation of Pleurospermum hookeri (Apiaceae), and its interspecific relationships with two close relatives in the genus Pleurospermum

In order to clarify the interspecific relationships of a lineage in Pleurospermum, P. hookeri C. B. Clarke, P. yunnanense Franch., and P. giraldii Diels, and to understand intraspecific divergence of P. hookeri, a phylogeographic study was carried out based on 198 individuals from 24 populations. Three chloroplast DNA regions, ndhF-rpl32, trnL-trnF, and trnQ-rps16, were sequenced in the present study. The genetic relationship between P. hookeri and P. giraldii is not as close as previously assumed. Pleurospermum hookeri and P. giraldii may originate from an unknown ancestor located in the Qinling region. Pleurospermum yunnanensewas found to be the closest relative of P. hookeri in all the species included in the phylogenetic analysis. The two haplotypes identified from P. yunnanense are shared with P. hookeri, which is potentially a result of both incomplete linkage sorting and introgression. Three large divergences within P. hookeri were identified, located at the northeastern edge, southeastern edge, and platform of the Qinghai–Tibet Plateau (QTP), respectively. Long-term history can explain the deep intraspecific divergence of P. hookeri. The uplift of the QTP played a key role in that divergence, and then were the climatic changes in the Quaternary. In addition, we found one refugium at the northeastern edge of the QTP, one at the southeastern edge, and at least one in the Hengduan Mountains region on the platform of the QTP.     

Pliocene–Pleistocene ecological niche evolution shapes the phylogeography of a Mediterranean plant group

Estimating species ability to adapt to environmental changes is crucial to understand their past and future response to climate change. The Mediterranean Basin has experienced remarkable climatic changes since the Miocene, which have greatly influenced the evolution of the Mediterranean flora. Here, we examine the evolutionary history and biogeographic patterns of two sedge sister species (Carex, Cyperaceae) restricted to the western Mediterranean Basin, but with Pliocene fossil record in central Europe. In particular, we estimated the evolution of climatic niches through time and its influence in lineage differentiation. We carried out a dated phylogenetic–phylogeographic study based on seven DNA regions (nDNA and ptDNA) and fingerprinting data (AFLPs), and modelled ecological niches and species distributions for the Pliocene, Pleistocene and present. Phylogenetic and divergence time analyses revealed that both species form a monophyletic lineage originated in the late Pliocene–early Pleistocene. We detected clear genetic differentiation between both species with distinct genetic clusters in disjunct areas, indicating the predominant role of geographic barriers limiting gene flow. We found a remarkable shift in the climatic requirements between Pliocene and extant populations, although the niche seems to have been relatively conserved since the Pleistocene split of both species. This study highlights how an integrative approach combining different data sources and analyses, including fossils, allows solid and robust inferences about the evolutionary history of a plant group since the Pliocene.     

Molecular phylogeography and intraspecific divergence of Spiraea alpina (Rosaceae) distributed in the Qinghai-Tibetan Plateau and adjacent regions inferred from nrDNA

The Qinghai-Tibetan plateau (QTP) uplift had a decisive effect on climatic and eco-environmental evolution in East Asia during the Quaternary. In the current study phylogeographic structure and diversification history of Spiraea alpina across the QTP were investigated for the first time based on nuclear internal transcribed spacer. The nuclear internal transcribed spacers (ITS_1a–ITS₄) were generated for a total of 284 individuals distributed within 31 natural populations. A clear phylogeographic structure was found for S. alpina. The results showed that this species colonized in three different glacial refugia during the Quaternary extensive glaciation and expanded during the Interglacial period. Analysis of molecular variance (AMOVA) showed 74.13% genetic diversity among populations and 25.87% genetic variation within populations with distinct phylogeographic structure (F_ST = 0.741*). The estimated divergence time revealed that the main lineages of S. alpina diversified during the Quaternary 1.2–0.6 million years ago. The study concluded that severe climatic oscillations during Quaternary and the uplift of QTP had a profound effect on intraspecific divergence of S. alpina.     

Evolutionary history of Purple cone spruce (Picea purpurea) in the Qinghai–Tibet Plateau: homoploid hybrid origin and Pleistocene expansion

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.     

Initial diversification,glacial survival,and continuous range expansion of Gentiana straminea (Gentianaceae) in the Qinghai–Tibet Plateau

The Qinghai–Tibet Plateau (QTP) has been considered as one of the most sensitive regions to climate change on Earth, and the growth and distribution of alpine species on this plateau have been suggested to depend greatly on their ability to survive within a small range of temperatures. However, the responses of most species in the QTP to the Quaternary climatic oscillation remain largely unknown. We sequenced two cpDNA fragments and nrITS to examine genetic variations in 22 natural populations across the range of distribution in this region to investigate the phylogeographical distributional pattern of Gentiana straminea (Gentianaceae) in the QTP. The high haplotype diversity from populations on the platform suggested the existence of intraspecific diversification. Molecular dating estimated that all haplotypes have differentiated before the Last Glacial Maximum (LGM). Moreover, the haplotype distribution map based on both cpDNA and nrDNA data suggested expansions from QTP to its outer edges. Finally, ecological niche modeling further demonstrated the glacial survival of this species on the platform and continuous expansion to the platform edge. These findings imply that G. straminea should have experienced initial diversification, glacial survival on the platform, and continuous expansion to the QTP edge during the glacial period.     

Genetic variation and phylogeographic history of <Emphasis Type="Italic">Picea likiangensis</Emphasis> revealed by RAPD markers

Repeated cycles of retreat and recolonization during the Quaternary ice ages are thought to have greatly influenced current species distributions and their genetic diversity. It remains unclear how this climatic oscillation has affected the distribution of genetic diversity between populations of wind-pollinated conifers in the Qinghai-Tibetan region. In this study, we investigated the within-species genetic diversity and phylogenetic relationships of Picea likiangensis, a dominant forest species in this region using polymorphic DNA (RAPD) markers. Our results suggest that this species has high overall genetic diversity, with 85.42% of loci being polymorphic and an average expected heterozygosity (H _E) of 0.239. However, there were relatively low levels of polymorphism at population levels and the differences between populations were not significant, with percentages of polymorphic bands (PPB) ranging from 46.88 to 69.76%, Nei’s gene diversity (H _E) from 0.179 to 0.289 and Shannon’s indices (Hpop) from 0.267 to 0.421. In accordance with our proposed hypothesis, a high level of genetic differentiation among populations was detected based on Nei’s genetic diversity (G _ST = 0.256) and AMOVA analysis (Phi _st = 0.236). Gene flow between populations was found to be limited (Nm = 1.4532) and far lower than reported for other conifer species with wide distribution ranges from other regions. No clusters corresponding to three morphological varieties found in the south, north and west, respectively, were detected in either UPGMA or PCO analyses. Our results suggest that this species may have had different refugia during the glacial stages in the southern region and that the northern variety may have multiple origins from these different refugia.     

Late Pleistocene climate change promoted divergence between Picea asperata and P. crassifolia on the Qinghai–Tibet Plateau through recent bottlenecks

Divergence during the early stage of speciation can be driven by a population bottleneck via reduced gene flow and enhanced lineage sorting. In this study, we aimed to examine whether such bottlenecks occurred during the initial speciation of two closely related spruce species Picea asperata and P. crassifolia occurring on the Qinghai–Tibet Plateau (QTP). We analyzed sequences of three chloroplast, two mitochondrial DNA fragments and a further 13 nuclear loci from 216 individuals of the two species. Both species showed a low level of genetic diversity in contrast to other congeners occurring in the QTP and adjacent regions. The estimated population sizes of P. asperata and P. crassifolia are less than the ancestral population size before splitting. These results together with multiple statistical tests (Tajima's D, Fu and Li's D* and F*) suggest that these two species underwent recent bottlenecks. Based on approximate Bayesian computation (ABC), we also determined that the period of the population shrinkage was consistent with the interspecific divergence during the late Pleistocene. The reduced population sizes and the divergent selection may together have triggered the initial divergence under high gene flow between these two species. Our results therefore highlight the importance of climatic oscillations during the late Pleistocene in promoting speciation through changing demographic sizes of the ancestral species on the QTP and in adjacent regions.     

Genetic divergence of the Mesoamerican azure‐crowned hummingbird (Amazilia cyanocephala,Trochilidae) across the Motagua‐Polochic‐Jocotán fault system

The cloud forests of Mesoamerica are highly endangered habitats and the existence of narrowly distributed cryptic endemics will increase the number of taxa at potential risk of extinction. Here, we investigate genetic divergence between populations of the azure‐crowned hummingbird (Amazilia cyanocephala), a species complex of endemic hummingbirds to the montane forests of Mesoamerica, by analysing DNA sequences of four mitochondrial markers, morphological data and ecological niche modelling. Our results revealed the presence of two mtDNA lineages corresponding to subspecies A. c. cyanocephala distributed from Tamaulipas to Chiapas in Mexico and Amazilia c. guatemalensis distributed from southern Chiapas to Guatemala. The lineage split can be explained as a consequence of relative isolation of the populations in the different mountain ranges separated by the Motagua‐Polochic‐Jocotán fault system and corresponds to differences in morphology and to the lack of overlap in environmental space between subspecies. The divergence time estimates do not support the proposed model of a highly constrained temporal window at the end of the Pliocene as divergence at this barrier between cyanocephala and guatemalensis and splits of other bird taxa occurred during the Pleistocene.     

Phylogeography and evolution of two closely related oak species (<Emphasis Type="Italic">Quercus</Emphasis>) from north and northeast China

Orographic and climatic influences during the Pleistocene have had a crucial role on interspecific divergence and population demography during speciation. However, associations between demographic histories of closely related species and related climatic events, especially in north and northeast China, are still underexplored. Genetic analyses with four chloroplast DNA and two nuclear genes and species distribution modeling were used for two closely related oak species (Quercus liaotungensis and Quercus mongolica) to test if their interspecific divergence and phylogeographical histories were possibly related to the Pleistocene-era climatic events. Potential divergence of the two oak species was estimated at about 0.92–2.15 Ma. Species distribution models and genetic data showed varying phylogeographical histories and spatial population structures between the two oaks, leading to different patterns of interspecific gene flow between the chloroplast and nuclear genes. The results indicate that speciation event between the two species is recent and may have been triggered by geological and climatic fluctuations linked to the upheavals of the Qinghai-Tibetan Plateau at the Pliocene/Pleistocene boundary. The two closely related oaks possess varying population demography during the interglacial-glacial climatic oscillations of the Quaternary, probably due to the various niche adaptations among different distribution ranges across their species trajectories.     

Ancient east-west divergence,recent admixture,and multiple marginal refugia shape genetic structure of a widespread oak species (<Emphasis Type="Italic">Quercus acutissima</Emphasis>) in China

The interactions of major topographic features with climate fluctuations over time are key factors that influence the spatial genetic structure of trees species. As an important characteristic of China’s terrain, the boundary between the mountainous West China and the lowland East China has been found to contribute to the divergence of many species. In this study, we examined whether the stepped geomorphology of China has shaped the genetic structure of Quercus acutissima, a widespread oak species native to East Asia, at a much larger geographic scale. Thirty natural populations of the species were sampled across the entire distribution in China and genotyped using ten nuclear microsatellite markers. To determine the genetic structure and demographic history, we conducted Bayesian cluster analysis and approximate Bayesian computation. A significant east-west divergence was detected at the species level, with western populations presenting much stronger differentiation and lower genetic diversity than eastern populations. The divergence between the groups of Southwest China and East China was estimated to have occurred 3.55–2.37 million years ago (Ma), closely associated with the recent uplift of the Qinghai-Tibetan Plateau (QTP) and climate changes during the late Pliocene to the early Pleistocene. The Central China group may arise from the admixture of these two groups, which was dated to 0.83–0.55 Ma, a period covering several interglacial intervals with unusually warm and wet climate in Central China. The spatial genetic structure of Q. acutissima indicates the possibility of central-marginal dispersal dynamics along latitude and existence of multiple refugia at distribution edges. Our findings highlight that the stepped geomorphology in China and climatic changes since the Pliocene play an important role in shaping the intraspecific genetic structure of widespread species in East Asia.     

Population‐genomic approach reveals adaptive floral divergence in discrete populations of a hawk moth‐pollinated violet

Local adaptation to contrasting biotic or abiotic environments is an important evolutionary step that presumably precedes floral diversification at the species level, yet few studies have demonstrated the adaptive nature of intraspecific floral divergence in wild plant populations. We combine a population‐genomic approach with phenotypic information on floral traits to examine whether the differentiation in metric floral traits exhibited by 14 populations of the southern Spanish hawk moth‐pollinated violet Viola cazorlensis reflects adaptive divergence. Screening of many amplified fragment length polymorphism (AFLP) loci using a multiple‐marker‐based neutrality test identified nine outlier loci (2.6% of the total) that departed from neutral expectations and were potentially under selection. Generalized analysis of molecular variance revealed significant relationships between genetic distance and population divergence in three floral traits when genetic distance was based on outlier loci, but not when it was based on neutral ones. Population means of floral traits were closely correlated with population scores on the first principal coordinate axis of the genetic distance matrix using outlier loci, and with the allelic frequencies of four of the outlier loci. Results strongly support the adaptive nature of intraspecific floral divergence exhibited by V. cazorlensis and illustrate the potential of genome scans to identify instances of adaptive divergence when used in combination with phenotypic information.     

Intraspecific divergences of Rhodiola alsia (Crassulaceae) based on plastid DNA and internal transcribed spacer fragments

In the present study, we used two maternally inherited plastid DNA intergenic spacers, rpl20‐rps12 and trnS‐trnG, and the biparentally inherited nuclear ribosomal internal transcribed spacer (ITS) region to explore genetic variation and phylogeographical history of Rhodiola alsia, a herb endemic to the Qinghai‐Tibetan Plateau (QTP). Based on range‐wide sampling (18 populations and 227 individuals), we detected 45 plastid DNA haplotypes and 19 ITS sequence types. Only three plastid DNA haplotypes were widespread; most haplotypes were restricted to single sites or to neighbouring populations. Analysis of molecular variance revealed that most of the genetic variance was found within populations (51.24%) but that populations were also distinct (F_ST = 0.48759). We found three areas with relatively high plastid DNA diversity and these could further be recognized as potentially isolated divergence centres based on the ITS sequence type distribution. These represent three potentially isolated glacial refugia for R. alsia: one of them has long been recognized as an important refugium on the south‐eastern edge of the QTP, whereas the others are new and located in the north and south of the Tanggula Mountains on the plateau platform. Divergence time estimates based on ITS suggest that the main lineages of R. alsia diverged from each other 0.35–0.87 Mya, indicating that climatic oscillations during the Pleistocene may have been an important driver of intraspecific divergence in R. alsia. Rhodiola alsia probably experienced a phylogeographical history of retreat to isolated glacial refugia during Quaternary glaciations that led to different degrees of allopatric intraspecific divergence. © 2011 The Linnean Society of London, Botanical Journal of the Linnean Society, 2012, 168 , 204–215.     

The impact of early Quaternary climate change on the diversification and population dynamics of a South American cactus species

Aim

Climatic oscillations have been suggested to promote speciation and changes in species distributions, mostly in connection with the Last Glacial Maximum (LGM). However, the LGM is just the most recent of the 20+ glacial‐interglacial periods that characterise the Quaternary. Here, we investigate the role of climatic changes and geomorphological features in shaping the evolution, distribution and population dynamics of the South American cactus Cereus hildmannianus.

Location

South‐eastern South America.

Methods

We built a large fossil‐calibrated phylogeny for cacti (family Cactaceae), comprising 128 species distributed in all subfamilies, using a Bayesian relaxed clock. We used the results to derive a secondary calibration for a population‐level phylogeny in C. hildmannianus. We amplified two plastid (trnQ‐5′rps16 and psbJ‐petA) and one nuclear marker (PhyC) for 24 populations. We estimated population dynamics, ancestral areas, and species distribution models to infer the clade's evolutionary history in time and space.

Results

Our results show a major population divergence of C. hildmannianus at c. 2.60 Ma, which is strikingly coincident with the transition of the Pliocene–Pleistocene and onset of Quaternary glaciations. This was followed by a complex phylogeographic scenario involving population expansions across ecologically diverse regions.

Main conclusions

Contrary to the dominant research focus on the LGM, our study indicates a major impact of the first Quaternary glaciation on the distribution and population divergence of a South American plant species. Further intraspecific events seem related to successive climatic changes and geomorphology, including the development of the coastal plain and its peculiar diversity. We propose that the first Quaternary glaciation acted as a major evolutionary bottleneck, whereby many warm‐adapted lineages succumbed, while those that survived could diversify and better cope with subsequent climatic oscillations.     

Genetic structure and ecological niche segregation of Indian gray mongoose (Urva edwardsii) in Iran

Combining genetic data with ecological niche models is an effective approach for exploring climatic and nonclimatic environmental variables affecting spatial patterns of intraspecific genetic variation. Here, we adopted this combined approach to evaluate genetic structure and ecological niche of the Indian gray mongoose (Urva edwardsii) in Iran, as the most western part of the species range. Using mtDNA, we confirmed the presence of two highly differentiated clades. Then, we incorporated ensemble of small models (ESMs) using climatic and nonclimatic variables with genetic data to assess whether genetic differentiation among clades was coupled with their ecological niche. Climate niche divergence was also examined based on a principal component analysis on climatic factors only. The relative habitat suitability values predicted by the ESMs for both clades revealed their niche separation. Between‐clade climate only niche comparison revealed that climate space occupied by clades is similar to some extent, but the niches that they utilize differ between the distribution ranges of clades. We found that in the absence of evidence for recent genetic exchanges, distribution models suggest the species occurs in different niches and that there are apparent areas of disconnection across the species range. The estimated divergence time between the two Iranian clades (4.9 Mya) coincides with the uplifting of the Zagros Mountains during the Early Pliocene. The Zagros mountain‐building event seems to have prevented the distribution of U. edwardsii populations between the western and eastern parts of the mountains as a result of vicariance events. Our findings indicated that the two U. edwardsii genetic clades in Iran can be considered as two conservation units and can be utilized to develop habitat‐specific and climate change‐integrated management strategies.