Phylogeography of the South China field mouse (Apodemus draco) on the southeastern Tibetan Plateau reveals high genetic diversity and glacial refugia

The southeastern margin of the Tibetan Plateau (SEMTP) is a particularly interesting region due to its topographic complexity and unique geologic history, but phylogeographic studies that focus on this region are rare. In this study, we investigated the phylogeography of the South China field mouse, Apodemus draco, in order to assess the role of geologic and climatic events on the Tibetan Plateau in shaping its genetic structure. We sequenced mitochondrial cytochrome b (cyt b) sequences in 103 individuals from 47 sampling sites. In addition, 23 cyt b sequences were collected from GenBank for analyses. Phylogenetic, demographic and landscape genetic methods were conducted. Seventy-six cyt b haplotypes were found and the genetic diversity was extremely high (π = 0.0368; h = 0.989). Five major evolutionary clades, based on geographic locations, were identified. Demographic analyses implied subclade 1A and subclade 1B experienced population expansions at about 0.052-0.013 Mya and 0.014-0.004 Mya, respectively. The divergence time analysis showed that the split between clade 1 and clade 2 occurred 0.26 Mya, which fell into the extensive glacial period (EGP, 0.5-0.17 Mya). The divergence times of other main clades (2.20-0.55 Mya) were congruent with the periods of the Qingzang Movement (3.6-1.7 Mya) and the Kun-Huang Movement (1.2-0.6 Mya), which were known as the most intense uplift events in the Tibetan Plateau. Our study supported the hypothesis that the SEMTP was a large late Pleistocene refugium, and further inferred that the Gongga Mountain Region and Hongya County were glacial refugia for A. draco in clade 1. We hypothesize that the evolutionary history of A. draco in the SEMTP primarily occurred in two stages. First, an initial divergence would have been shaped by uplift events of the Tibetan Plateau. Then, major glaciations in the Pleistocene added complexity to its demographic history and genetic structure.     

Allopatric divergence and phylogeographic structure of the plateau zokor (Eospalax baileyi), a fossorial rodent endemic to the Qinghai–Tibetan Plateau

Aim Most species of temperate regions are believed to have shifted to lower latitudes or elevations during the glacial periods of the Quaternary. In this study we test whether this phylogeographic assumption is also true for the plateau zokor (Eospalax baileyi), a fossorial rodent endemic to the climate-sensitive Qinghai–Tibetan Plateau (QTP), which ranges in elevation from 2600 to 4600 m. Location The QTP of western China. Methods Phylogeographic analyses were conducted based on the mitochondrial cytochrome b gene sequences of 193 individuals from 20 populations over the entire range of the species. Results A total of 54 haplotypes identified in the present study clustered into four geographically correlated clades located in the interior of the QTP (clade A) and at the plateau edge (B, C and D). Molecular calibrations suggest that the interior plateau (A) and plateau-edge (B–D) clades diverged at 1.2 Ma and that the three plateau-edge clades diverged between 0.85 and 0.80 Ma. These estimates are concordant with diastrophism and glaciation events in the QTP. Coalescent tests rejected both the hypothesis that all current populations originated from a single refugium at a low elevation during the Last Glacial Maximum (LGM) and the hypothesis that the two lineages diverged during the LGM. The tests instead supported the hypothesis that there were four refugia during the LGM, and that the four clades diverged prior to the late Pleistocene. Main conclusions Our results suggest that Quaternary diastrophisms and glaciations repeatedly promoted allopatric divergence of the plateau zokor into geographical clades, and that these regional clades subsequently persisted at high elevations, rather than migrating to the low-elevation plateau edge during subsequent glacial ages.     

Pleistocene persistence and expansion in tarantulas on the Colorado Plateau and the effects of missing data on phylogeographical inferences from RADseq

Montane species endemic to the “sky islands” of the North American southwest were significantly impacted by changing climates during the Pleistocene. We combined mitochondrial and genomic data with species distribution modelling to determine whether Aphonopelma marxi, a large tarantula from the nearby Colorado Plateau, was similarly impacted by glacial climates. Genetic analyses revealed that the species comprises three main clades that diverged in the Pleistocene. A clade distributed along the Mogollon Rim appears to have persisted in place during glacial conditions, whereas the other two clades probably colonized central and northeastern portions of the species' range from refugia in canyons. Climate models support this hypothesis for the Mogollon Rim, but late glacial climate data appear too coarse to detect suitable areas in canyons. Locations of canyon refugia could not be inferred from genomic analyses due to missing data, encouraging us to explore the effect of missing loci in phylogeographical inferences using RADseq. Results from analyses with varying amounts of missing data suggest that samples with large amounts of missing data can still improve inferences, and the specific loci that are missing matters more than the number of missing loci. This study highlights the profound impact of Pleistocene climates on tarantulas endemic to the Colorado Plateau, as well as the mixed nature of the region's fauna. Some animals recently colonized from nearby deserts as glacial climates receded, whereas others, like tarantulas, appear to have persisted on the Mogollon Rim and in refugia associated with the region's famous river‐cut canyons.     

Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the case of the bicoloured shrew (Crocidura leucodon, Soricidae)

We sequenced 1077 bp of the mitochondrial cytochrome b gene and 511 bp of the nuclear Apolipoprotein B gene in bicoloured shrew (Crocidura leucodon, Soricidae) populations ranging from France to Georgia. The aims of the study were to identify the main genetic clades within this species and the influence of Pleistocene climatic variations on the respective clades. The mitochondrial analyses revealed a European clade distributed from France eastwards to north-western Turkey and a Near East clade distributed from Georgia to Romania; the two clades separated during the Middle Pleistocene. We clearly identified a population expansion after a bottleneck for the European clade based on mitochondrial and nuclear sequencing data; this expansion was not observed for the eastern clade. We hypothesize that the western population was confined to a small Italo-Balkanic refugium, whereas the eastern population subsisted in several refugia along the southern coast of the Black Sea.     

Phylogeography and post-glacial colonization patterns of the rainbow darter, Etheostoma caeruleum (Teleostei: Percidae)

Aim To examine the effects of historical climate change and drainage isolation on the distribution of mitochondrial DNA cytochrome b genetic variation within the rainbow darter, Etheostoma caeruleum (Percidae: Etheostomatinae). Location Eastern North American streams including tributaries to the Mississippi River, Great Lakes, Potomac River and Hudson Bay drainages. Methods Parsimony analyses, Bayesian analyses and haplotype networks of mitochondrial DNA sequences. Results Four major clades were recovered from sampled populations of E. caeruleum. Three of four clades are distributed in the western portion of the species’ range (primarily west of the Mississippi River). Samples from this region do not form a monophyletic group, and sequences often vary greatly between samples from adjacent stream systems (up to 7.2% divergence). A basal clade includes samples from the White River system in the Ozark Highlands. The northern Ozarks–upper Midwest clade includes samples from Missouri River tributaries and the upper Midwest (Hudson Bay, upper Mississippi River, and western Lake Michigan drainage). The eastern clade is composed of individuals from the Ohio River, Great Lakes and Potomac River. The Mississippi River corridor clade includes samples from middle and lower Mississippi River tributaries. Main conclusions The four major clades of E. caeruleum are deep allopatric lineages with well‐defined boundaries and have additional phylogeographical structure within each clade. The Ozark Highlands have the greatest levels of diversity relative to distributional area, with marked cytochrome b subdivisions between adjacent stream systems. Samples from previously glaciated areas do not have a subset of the cytochrome b diversity found in unglaciated areas, but four separate source areas are identified based on phylogenetic analyses. Dispersal into previously glaciated areas followed several known glacial outlets and, based on sequence divergence between populations, may have occurred during different glacial or interglacial stages. The disjunct distribution and cytochrome b pattern of E. caeruleum in the Mississippi River corridor clade is consistent with late Pleistocene and Recent changes in the course and characteristics of the middle and lower Mississippi River. Phylogeographical boundaries between clades of E. caeruleum correspond to independent sources of biogeographical information and provide insight into historical stream drainage relationships, post‐glacial colonization and drainage isolation patterns.     

Postglacial colonisation patterns and the role of isolation and expansion in driving diversification in a passerine bird

Pleistocene glacial cycles play a major role in diversification and speciation, although the relative importance of isolation and expansion in driving diversification remains debated. We analysed mitochondrial DNA sequence data from 15 great reed warbler (Acrocephalus arundinaceus) populations distributed over the vast Eurasian breeding range of the species, and revealed unexpected postglacial expansion patterns from two glacial refugia. There were 58 different haplotypes forming two major clades, A and B. Clade A dominated in Western Europe with declining frequencies towards Eastern Europe and the Middle East, but showed a surprising increase in frequency in Western and Central Asia. Clade B dominated in the Middle East, with declining frequencies towards north in Central and Eastern Europe and was absent from Western Europe and Central Asia. A parsimonious explanation for these patterns is independent postglacial expansions from two isolated refugia, and mismatch distribution analyses confirmed this suggestion. Gene flow analyses showed that clade A colonised both Europe and Asia from a refugium in Europe, and that clade B expanded much later and colonised parts of Europe from a refugium in the Middle East. Great reed warblers in the eastern parts of the range have slightly paler plumage than western birds (sometimes treated as separate subspecies; A. a. zarudnyi and A. a. arundinaceus, respectively) and our results suggest that the plumage diversification took place during the easterly expansion of clade A. This supports the postglacial expansion hypothesis proposing that postglacial expansions drive diversification in comparatively short time periods. However, there is no indication of any (strong) reproductive isolation between clades and our data show that the refugia populations became separated during the last glaciation. This is in line with the Pleistocene speciation hypothesis invoking that much longer periods of time in isolation are needed for speciation to occur.     

Further Analysis of Intraspecific Sequence Variation of Chloroplast DNA in Primula cuneifolia Ledeb. (Primulaceae): Implications for Biogeography of the Japanese Alpine Flora

Primula cuneifolia Ledeb. (Primulaceae), we analyzed intraspecific variation of the nucleotide sequences of non-coding regions of chloroplast DNA: the intergenic spacers between trnT (UGU) and the trnL (UAA) 5′exon, the trnL (UAA) 3′exon and trnF (GAA), and atpB and rbcL. In 20 populations of P. cuneifolia, 22 nucleotide substitutions and five insertions/deletions were inferred, and their genetic distances ranged from 0.001 to 0.008. Eight distinct haplotypes could be recognized and each haplotype was found to be geographically structured. Three major clades (the Northern, Hokkaido and Southern clades) were revealed in phylogenetic analyses of the haplotypes. The haplotypes of the Northern clade had a wider distribution area in the populations of Mt. Rausu and Rishiri Island of eastern and northern Hokkaido in Japan, northward to Unalaska Island in the Aleutians, and those of the Hokkaido clade were distributed in the populations of central Hokkaido and Mt. Iwaki of the northern Honshu in Japan; in addition, those of the Southern clade were observed only in the populations of the central Honshu. It was shown that the genetic diversifications of the Southern clade were higher than those of the Northern and Hokkaido clades. Furthermore, it was shown that the topology within the Southern clade was hierarchical, and the haplotypes of the Southern populations in the clade were derivative. From these results, we concluded that the cpDNA haplotypes of the three clades in P. cuneifolia arose and assumed the present distribution areas through several cycles of glacial advance and retreat in the Pleistocene. Received 24 June 1998/ Accepted in revised form 28 December 1998     

Surviving glacial ages within the Biotic Gap: phylogeography of the New Zealand cicada Maoricicada campbelli

Aim New Zealand is an ideal location in which to investigate the roles of landscape and climate change on speciation and biogeography. An earlier study of the widespread endemic cicada Maoricicada campbelli (Myers) found two phylogeographically distinguishable major clades – northern South Island plus North Island (northern‐SI + NI) and Otago. These two clades appeared to have diverged on either side of an area of the South Island known as the Biotic Gap. We sampled more intensively to test competing theories for this divergence. We aimed to discover if M. campbelli had survived within the Biotic Gap during recent glacial maxima, and if predicted areas of secondary contact between the two major clades existed. Location New Zealand. Methods We analysed mitochondrial DNA sequences (1520 bp; 212 individuals; 91 populations) using phylogenetic (maximum likelihood, Bayesian), population genetic (analysis of molecular variance) and molecular dating methods (Bayesian relaxed clock with improved priors). Results We found strong geographical structuring of genetic variation. Our dating analyses suggest that M. campbelli originated 1.83–2.58 Ma, and split into the two major clades 1.45–2.09 Ma. The main subclades in the northern‐SI + NI clade arose almost simultaneously at 0.69–1.03 Ma. Most subclades are supported by long internal branches and began to diversify 0.40–0.78 Ma. We found four narrow areas of secondary contact between the two major clades. We also found a difference between calling songs of the Otago vs. northern‐SI + NI clades. Main conclusions Phylogeographical patterns within M. campbelli indicate an early Pleistocene split into two major clades, followed by late Pleistocene range expansion and in situ population differentiation of subclades. The northern‐SI + NI clade diversified so rapidly that the main subclade relationships cannot be resolved, and we now have little evidence for a disjunction across the Biotic Gap. Structure within the main subclades indicates rapid divergence after a common bottlenecking event, perhaps attributable to an extremely cold glacial maximum at c. 0.43 Ma. Clade structure and dating analyses indicate that M. campbelli survived in many refugia during recent glacial maxima, including within the Biotic Gap. The narrow overlap between the two major clades is attributed to recent contact during the current interglacial and slow gene diffusion. The two major clades appear to be in the early stages of speciation based on genetic and behavioural differences.     

Coalescent-based hypothesis testing supports multiple Pleistocene refugia in the Pacific Northwest for the Pacific giant salamander (Dicamptodon tenebrosus)

Phylogeographic patterns of many taxa are explained by Pleistocene glaciation. The temperate rainforests within the Pacific Northwest of North America provide an excellent example of this phenomenon, and competing phylogenetic hypotheses exist regarding the number of Pleistocene refugia influencing genetic variation of endemic organisms. One such endemic is the Pacific giant salamander, Dicamptodon tenebrosus. In this study, we estimate this species' phylogeny and use a coalescent modeling approach to test five hypotheses concerning the number, location and divergence times of purported Pleistocene refugia. Single refugium hypotheses include: a northern refugium in the Columbia River Valley and a southern refugium in the Klamath-Siskiyou Mountains. Dual refugia hypotheses include these same refugia but separated at varying times: last glacial maximum (20,000 years ago), mid-Pleistocene (800,000 years ago) and early Pleistocene (1.7 million years ago). Phylogenetic analyses and inferences from nested clade analysis reveal distinct northern and southern lineages expanding from the Columbia River Valley and the Klamath-Siskiyou Mountains, respectively. Results of coalescent simulations reject both single refugium hypotheses and the hypothesis of dual refugia with a separation date in the late Pleistocene but not hypotheses predicting dual refugia with separation in early or mid-Pleistocene. Estimates of time since divergence between northern and southern lineages also indicate separation since early to mid-Pleistocene. Tests for expanding populations using mismatch distributions and 'g' distributions reveal demographic growth in the northern and southern lineages. The combination of these results provides strong evidence that this species was restricted into, and subsequently expanded from, at least two Pleistocene refugia in the Pacific Northwest.     

The Balkans and the colonization of Europe: the post‐glacial range expansion of the wild boar,Sus scrofa

Aim We focus on the biogeographical role of the Balkan Peninsula as a glacial refugium and source of northward post‐glacial dispersal for many European taxa. Specifically, we analysed the genetic structure and variation of wild boar (Sus scrofa) samples primarily from Greece, a region that has repeatedly served as a glacial refugium within the Balkan Peninsula. Location Continental Greece, the Aegean island of Samos and Bulgaria. Methods We analysed wild boar samples from 18 localities. Samples from common domestic breeds were also examined to take into account interactions between wild and domesticated animals. Phylogenetic analyses were carried out on a 637‐bp fragment of the mitochondrial DNA control region in 200 wild boar and 27 domestic pigs. The sequences were also compared with 791 Eurasian wild boar and domestic pig D‐loop sequences obtained from GenBank. Results Ninety‐four haplotypes were identified in the European wild boar data set, of which 68 were found in the Balkan samples and assigned to two previously described clades: the E1 European and Near Eastern clades. All of the continental samples clustered in the E1 clade and the samples from Samos fell into the Near Eastern clade, consistent with the island’s proximity to Asia Minor. Intriguingly, 62 novel haplotypes were identified and are found exclusively in the Balkans. Only six haplotypes were shared between wild boar and domestic pigs. Main conclusions Our data reveal numerous novel and geographically restricted haplotypes in wild boar populations, suggesting the presence of separate refugia in the Balkans. Our analyses support the hypothesis of a post‐glacial wild boar expansion consistent with the leading edge model, north and west from modern day Greece, and suggest little maternal introgression of Near Eastern and domestic haplotypes into wild Balkan populations.     

Persistent genetic signatures of colonization in Brachionus manjavacas rotifers in the Iberian Peninsula

Recent phylogeographical assessments have consistently shown that continental zooplankton display high levels of population subdivision, despite the high dispersal capacity of their diapausing propagules. As such, there is an apparent paradox between observed cosmopolitanism in the zooplankton that is associated with long-distance dispersal, and strong phylogeographical structures at a regional scale. Such population dynamics, far from migration-drift equilibrium, have been shown in the rotifer species complex Brachionus plicatilis, a group of over a dozen species inhabiting salt lakes and coastal lagoons worldwide. Here we present the mitochondrial DNA phylogeography of one of these species, Brachionus manjavacas, in the Iberian Peninsula, where it often co-occurs with the morphologically similar species B. plicatilis sensu stricto. We obtained sequences from 233 individuals from diapausing eggs and clonal cultures from 16 lakes in the Iberian Peninsula, and a Tunisian lake. Two strongly supported deep mitochondrial DNA clades were found (A and B). Phylogenetic and nested clade analysis showed that clade A has a strong phylogeographical structure, with a strong similarity of phylogeographical patterns between B. manjavacas clade A and B. plicatilis s.s. These include (i) signatures of allopatric fragmentation between central and southern populations, and (ii) range expansions in the Iberian Peninsula, both likely to have occurred during the Pleistocene. We find evidence for a glacial refugium in the Guadiana basin. Clades A and B co-occurred in several of these lakes because of range expansion and secondary contact between both clades. The co-occurrence between B. plicatilis s.s. and B. manjavacas is not recent, and both species might have experienced similar environmental challenges during the Pleistocene. The strong correlation of genetic and geographical distance found suggests that historical events can lead to such correlation, mirroring the effects of 'isolation by distance' in equilibrium populations.     

Pleistocene climatic fluctuations explain the disjunct distribution and complex phylogeographic structure of the Southern Red-backed Salamander, <Emphasis Type="Italic">Plethodon serratus</Emphasis>

The southeastern United States (U.S.) has experienced dynamic climatic changes over the past several million years that have impacted species distributions. In many cases, contiguous ranges were fragmented and a lack of gene flow between allopatric populations led to genetic divergence and speciation. The Southern Red-backed Salamander, Plethodon serratus, inhabits four widely disjunct regions of the southeastern U.S.: the southern Appalachian Mountains, the Ozark Plateau, the Ouachita Mountains, and the Southern Tertiary Uplands of central Louisiana. We integrated phylogenetic analysis of mitochondrial DNA sequences (1399 base pairs) with ecological niche modeling to test the hypothesis that climate fluctuations during the Pleistocene drove the isolation and divergence of disjunct populations of P. serratus. Appalachian, Ozark, and Louisiana populations each formed well-supported clades in our phylogeny. Ouachita Mountain populations sorted into two geographically distinct clades; one Ouachita clade was sister to the Louisiana clade whereas the other Ouachita clade grouped with the Appalachian and Ozark clades but relationships were unresolved. Plethodon serratus diverged from its sister taxon, P. sherando, ~5.4 million years ago (Ma), and lineage diversification within P. serratus occurred ~1.9–0.6 Ma (Pleistocene). Ecological niche models showed that the four geographic isolates of P. serratus are currently separated by unsuitable habitat, but the species was likely more continuously distributed during the colder climates of the Pleistocene. Our results support the hypothesis that climate-induced environmental changes during the Pleistocene played a dominant role in driving isolation and divergence of disjunct populations of P. serratus.     

Phylogeography of the mountain chickadee (Poecile gambeli): diversification, introgression, and expansion in response to Quaternary climate change

Since the late 1990s, molecular techniques have fuelled debate about the role of Pleistocene glacial cycles in structuring contemporary avian diversity in North America. The debate is still heated; however, there is widespread agreement that the Pleistocene glacial cycles forced the repeated contraction, fragmentation, and expansion of the North American biota. These demographic processes should leave genetic 'footprints' in modern descendants, suggesting that detailed population genetic studies of contemporary species provide the key to elucidating the impact of the late Quaternary (late Pleistocene-Holocene). We present an analysis of mitochondrial DNA (mtDNA) variation in the mountain chickadee (Poecile gambeli) in an attempt to examine the genetic evidence of the impact of the late Quaternary glacial cycles. Phylogenetic analyses reveal two strongly supported clades of P. gambeli: an Eastern Clade (Rocky Mountains and Great Basin) and a Western Clade (Sierra Nevada and Cascades). Post-glacial introgression is apparent between these two clades in the Mono Lake region of Central California. Within the Eastern Clade there is evidence of isolation-by-distance in the Rocky Mountain populations, and of limited gene flow into and around the Great Basin. Coalescent analysis of genetic variation in the Western Clade indicates that northern (Sierra Nevada/Cascades) and southern (Transverse/Peninsular Ranges) populations have been isolated and evolving independently for nearly 60,000 years.     

Population genomic evidence for multiple Pliocene refugia in a montane‐restricted harvestman (Arachnida,Opiliones, Sclerobunus robustus) from the southwestern United States

The integration of ecological niche modelling into phylogeographic analyses has allowed for the identification and testing of potential refugia under a hypothesis‐based framework, where the expected patterns of higher genetic diversity in refugial populations and evidence of range expansion of nonrefugial populations are corroborated with empirical data. In this study, we focus on a montane‐restricted cryophilic harvestman, Sclerobunus robustus, distributed throughout the heterogeneous Southern Rocky Mountains and Intermontane Plateau of southwestern North America. We identified hypothetical refugia using ecological niche models (ENMs) across three time periods, corroborated these refugia with population genetic methods using double‐digest RAD‐seq data and conducted population‐level phylogenetic and divergence dating analyses. ENMs identify two large temporally persistent regions in the mid‐latitude highlands. Genetic patterns support these two hypothesized refugia with higher genetic diversity within refugial populations and evidence for range expansion in populations found outside hypothesized refugia. Phylogenetic analyses identify five to six genetically divergent, geographically cohesive clades of S. robustus. Divergence dating analyses suggest that these separate refugia date to the Pliocene and that divergence between clades pre‐dates the late Pleistocene glacial cycles, while diversification within clades was likely driven by these cycles. Population genetic analyses reveal effects of both isolation by distance (IBD) and isolation by environment (IBE), with IBD more important in the continuous mountainous portion of the distribution, while IBE was stronger in the populations inhabiting the isolated sky islands of the south. Using model‐based coalescent approaches, we find support for postdivergence migration between clades from separate refugia.     

Phylogenetic relationships in genus Niviventer (Rodentia: Muridae) in China inferred from complete mitochondrial cytochrome b gene

Chinese species of the genus Niviventer, predominantly distributed in the southeastern Tibetan Plateau and in Taiwan, are a diverse group and have not yet received a thorough molecular phylogenetic analysis. Here, we reconstructed the phylogenetic relationships of 32 specimens representing nine Chinese species of Niviventer, based on sequences of the complete mitochondrial cytochrome b gene. Maximum parsimony, maximum likelihood and Bayesian analysis resulted in three consistent trees, each supported by high bootstrap values. The results showed that the Niviventer species included here are monophyletic. The nine species were classified into three distinct clades: clade A with Niviventer brahma, N. confucianus, N. coxingi, N. culturatus, N. eha and N. fulvescens; clade B with N. andersoni and N. excelsior; clade C with N. cremoriventer. Our results also suggested that N. culturatus should be a valid species rather than a subspecies of N. confucianus. Divergence times among species were calibrated according to the middle-late Pleistocene (1.2-0.13 Mya) fossil records of N. confucianus. The results demonstrated that the first radiation event of the genus Niviventer occurred in early Pleistocene (about 1.66 Mya), followed by the divergence of clades A and B at about 1.46 Mya. Most of the extant Niviventer species appeared during early to middle Pleistocene (about 1.29-0.67 Mya). These divergence times are coincidental with the last uplift events of the Tibetan Plateau, Kun-Huang movement, Pleistocene glaciations and the vicariant formation of Taiwan Strait. Consequently geographical events and Pleistocene glaciations have played a great role in the diversification of Niviventer.     

Geographical differentiation of the Euchiloglanis fish complex (Teleostei: Siluriformes) in the Hengduan Mountain Region,China: Phylogeographic evidence of altered drainage patterns

The uplift of the Tibetan Plateau caused significant ecogeographical changes that had a major impact on the exchange and isolation of regional fauna and flora. Furthermore, Pleistocene glacial oscillations were linked to temporal large‐scale landmass and drainage system reconfigurations near the Hengduan Mountain Region and might have facilitated speciation and promoted biodiversity in southwestern China. However, strong biotic evidence supporting this role is lacking. Here, we use the Euchiloglanis fish species complex as a model to demonstrate the compound effects of the Tibetan Plateau uplift and Pleistocene glacial oscillations on species formation in this region. The genetic structure and geographical differentiation of the Euchiloglanis complex in four river systems within the Hengduan Mountain Region were deduced using the cytochrome b (cyt b) gene and 10 microsatellite loci from 360 to 192 individuals, respectively. The results indicated that the populations were divided into four independently evolving lineages, in which the populations from the Qingyi River and Jinsha River formed two sub‐lineages. Phylogenetic relationships were structured by geographical isolation, especially near drainage systems. Divergence time estimation analyses showed that the Euchiloglanis complex diverged from its sister clade Pareuchiloglanis sinensis at around 1.3 Million years ago (Ma). Within the Euchiloglanis complex, the divergence time between the Dadu–Yalong and Jinsha–Qingyi River populations occurred at 1.0 Ma. This divergence time was in concordance with recent geological events, including the Kun‐Huang Movement (1.2–0.6 Ma) and the lag time (<2.0 Ma) of river incision in the Hengduan Mountain Region. Population expansion signals were detected from mismatched distribution analyses, and the expansion times were concurrent with Pleistocene glacier fluctuations. Therefore, current phylogeographic patterns of the Euchiloglanis fish complex in the Hengduan Mountain Region were influenced by the uplift event of the Tibetan Plateau and were subsequently altered by paleo‐river transitions during the late Pleistocene glacial oscillations.     

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.     

Intraspecific sequence variation of chloroplast DNA inPedicularis chamissonis Steven (Scrophulariaceae) and geographic structuring of the Japanese “Alpine” plants

In order to clarify evolutionary patterns and processes of intraspecific diversification ofPedicularis chamissonis Steven, we analyzed intraspecific variation of the nucleotide sequences of non-coding regions of chloroplast DNA: the intergenic spacers betweentrnT (UGU) andtrnL (UAA) 5′exon,trnL (UAA) 3′exon andtrnF (GAA), andatpB andrbcL. In 24 populations ofP. chamissonis, 33 nucleotide substitutions and 12 insertions/deletions were inferred, and their genetic distances ranged from 0.001 to 0.014. Seventeen distinct cpDNA haplotypes could be recognized and each haplotype was found to be geographically structured. Two major clades (the Northern and Southern clades) were revealed in phylogenetic analyses of cpDNA haplotypes. The haplotypes of the Northern clade had a wider distribution area in the populations of Mts. lide of central Honshu in Japan, northward to Unalaska Island in the Aleutians. Relationships among most haplotypes were unresolved polytomies. On the other hand, the haplotypes of the Southern clade occurred from the populations of Mt. Gassan southwards to Mt. Arakawa of central Honshu. Within this clade, three subclades were clearly recognized. From these results, we concluded that the haplotypes of the Northern and Southern clades inP. chamissonis might have traveled down to Japanese Archipelago from the north in not a single glacial period.     

Timing of clade divergence and discordant estimates of genetic and morphological diversity in the Slender Madtom,Noturus exilis (Ictaluridae)

Phylogeographic relationships, the timing of clade diversification, and the potential for cryptic diversity in the Slender Madtom, Noturus exilis, was investigated using mitochondrial Cyt b, nuclear RAG2, shape analysis, and meristic and pigmentation data. Three well-supported and deeply divergent clades were recovered from analyses of genetic data: Little Red River (White River drainage) clade, Arkansas + Red River (Mississippi River) clade, and a large clade of populations from the rest of the range of the species. Recovered clades showed little to no diagnostic morphological differences, supporting previous hypotheses of morphological conservatism in catfishes, and indicating morphology may commonly underestimate diversity in this group of fishes. The Little Red River clade is the most distinct lineage of N. exilis with 11 POM pores (vs. 10 in other populations) and unique Cyt b haplotypes and RAG 2 alleles. However, treating it as a species separate from N. exilis would imply that the other major clades of N. exilis are more closely related to one another than they are to the Little Red River clade, which was not supported.The UCLN age estimate for Noturus was 23.9 mya (95% HPD: 13.49, 35.43), indicating a late Oligocene to early Miocene origin. The age of N. exilis was estimated as late Miocene at 9.7 mya (95% HPD: 5.32, 14.93). Diversification within the species spanned the late Miocene to mid-Pleistocene. The largest clade of N. exilis, which dates to the late Miocene, includes populations from the unglaciated Eastern and Interior Highlands as well as the previously glaciated Central Lowlands. Diversification of this clade coincides with a drastic drop in sea-level and diversification of other groups of Central Highlands fishes (Centrarchidae and Cyprinidae). Sub-clades dating to the Pleistocene show that northern populations occurring in previously glaciated regions resulted from dispersal from populations in the Ozarks up the Mississippi River following retreat of the Pleistocene glaciers. Pre-Pleistocene vicariance, such as drainage pattern changes of the Mississippi River, also played a prominent role in the history of the species. The incorporation of a temporal estimate of clade diversification revealed that in some instances, phylogeographic breaks shared with other aquatic species were best explained by different or persistent vicariant events through time, rather than a single shared event.     

Phylogeny and diversification of the cloud forest Morpho sulkowskyi group (Lepidoptera,Nymphalidae) in the evolving Andes

The monophyletic Morpho sulkowskyi butterfly group, endemic of Andean cloud forests, was studied to test the respective contributions of Mio‐Pliocene intense uplift period and Pleistocene glacial cycles on Andean biodiversity. We sampled nine taxa covering the whole geographical range of the group. Two mitochondrial and two nuclear genes were analysed using a Bayesian method. We established a dated phylogeny of the group using a relaxed clock method and a wide‐outgroup approach. To discriminate between two hypotheses, we used a biogeographical probabilistic method. Results suggest that the ancestor of the M. sulkowskyi group originated during the Middle–Late Miocene uplift of the Eastern Cordillera in northern Peru. Biogeographical inference suggests that the M. sulkowskyi and Morpho lympharis clades diverged in the northern Peruvian Andes. The subsequent divergences, from the Late Miocene to the Late Pliocene, should have resulted from a dispersal towards the Northern Andes (M. sulkowskyi clade), after the closure of the West Andean Portal separating the Central and Northern Andes, and a southwards dispersal along the Peruvian and Bolivian Eastern Cordilleras (M. lympharis clade). Only a few divergences occurred at the very end of the Pliocene or during the Pleistocene, a period when the more recent uplifts interfered with Pleistocene glacial cycles.