Evidence for multiple Pleistocene refugia in the postglacial expansion of the eastern tiger salamander, Ambystoma tigrinum tigrinum

Depauperate major histocompatibility complex variation in the endangered reticulated flatwoods salamander (Ambystoma bishopi)

Reticulated flatwoods salamander (Ambystoma bishopi) populations began decreasing dramatically in the 1900s. Contemporary populations are small, isolated, and may be susceptible to inbreeding and reduced adaptive potential because of low genetic variation. Genetic variation at immune genes is especially important as it influences disease susceptibility and adaptation to emerging infectious pathogens, a central conservation concern for declining amphibians. We collected samples from across the extant range of this salamander to examine genetic variation at major histocompatibility complex (MHC) class Iα and IIβ exons as well as the mitochondrial control region. We screened tail or toe tissue for ranavirus, a pathogen associated with amphibian declines worldwide. Overall, we found low MHC variation when compared to other amphibian species and did not detect ranavirus at any site. MHC class Iα sequencing revealed only three alleles with a nucleotide diversity of 0.001, while MHC class IIβ had five alleles with a with nucleotide diversity of 0.004. However, unique variation still exists across this species' range with private alleles at three sites. Unlike MHC diversity, mitochondrial variation was comparable to levels estimated for other amphibians with nine haplotypes observed, including one haplotype shared across all sites. We hypothesize that a combination of a historic disease outbreak and a population bottleneck may have contributed to low MHC diversity while maintaining higher levels of mitochondrial DNA variation. Ultimately, MHC data indicated that the reticulated flatwoods salamander may be at an elevated risk from infectious diseases due to low levels of immunogenetic variation necessary to combat novel pathogens.

Hybridizing salamanders experience accelerated diversification

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow-spotted mountain newt, Neurergus derjugini (Caudata: Salamandridae) in the mid-Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.

Post-Pleistocene differentiation in a Central Interior Highlands endemic salamander

AIM

For many endemic species with limited dispersal capacities, the relationship between landscape changes and species distributions is still unclear. We characterized the population structure of the endemic ringed salamander (Ambystoma annulatum) across its distribution in the Central Interior Highlands (CIH) of North America, an area of high species endemism, to infer the ecological and evolutionary history of the species.

METHODS

We sampled 498 individuals across the species distribution and characterized the population genetic structure using nuclear microsatellite and mitochondrial DNA (mtDNA) markers.

RESULTS

Ambystoma annulatum exist in two strongly supported nuclear genetic clusters across the CIH that correspond to a northern cluster that includes the Missouri Ozark populations and a southern cluster that includes the Arkansas and Oklahoma Ozarks and the Ouachita Mountains. Our demographic models estimated that these populations diverged approximately 2,700 years ago. Pairwise estimates of genetic differentiation at microsatellite and mtDNA markers indicated limited contemporary gene flow and suggest that genetic differentiation was primarily influenced by changes in the post-Pleistocene landscape of the CIH.

MAIN CONCLUSIONS

Both the geologic history and post-European settlement history of the CIH have influenced the population genetic structure of A. annulatum. The low mtDNA diversity suggests a retraction into and expansion out of refugial areas in the south-central Ozarks, during temperature fluctuations of the Pleistocene and Holocene epochs. Similarly, the estimated divergence time for the two nuclear clusters corresponds to changes in the post-Pleistocene landscape. More recently, decreased A. annulatum gene flow may be a result of increased habitat fragmentation and alteration post-European settlement.

Phylogeographic structure in three North American tent caterpillar species (Lepidoptera: Lasiocampidae): Malacosoma americana, M. californica, and M. disstria

While phylogeographic structure has been examined in many North American vertebrate species, insects have received much less attention despite their central ecological roles. The moth genus Malacosoma (Hübner, 1820), is an important group of forestry pests responsible for large-scale defoliation across much of the Nearctic and Palearctic. The present study uses sequence variation in the mitochondrial cytochrome c oxidase 1 (COI) gene to examine the population genetic structure of the three widespread Malacosoma species (M. americana, M. californica, and M. disstria). Populations of all three species showed highest diversity in the south, suggesting that modern populations derived from southern refugia with loss of variation as these lineages dispersed northwards. However, despite similar life histories and dispersal abilities, the extent of regional variation varied among the taxa. M. americana, a species restricted to eastern North America, showed much less genetic structure than the western M. californica or the widespread M. disstria. The regional differentiation in the latter reflects the likely derivation of modern lineages from several refugia, as well as taxonomic uncertainty in M. californica. In these respects, the three species of Malacosoma share phylogeographic patterns similar to those detected in vertebrates which are characterised by greater phylogeographic breaks in the western half of the continent and limited structure in the east.

Cyto-nuclear discordance suggests complex evolutionary history in the cave-dwelling salamander, Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.

Strong genetic structure corresponds to small-scale geographic breaks in the Australian alpine grasshopper Kosciuscola tristis

Background

Mountain landscapes are topographically complex, creating discontinuous ‘islands’ of alpine and sub-alpine habitat with a dynamic history. Changing climatic conditions drive their expansion and contraction, leaving signatures on the genetic structure of their flora and fauna. Australia’s high country covers a small, highly fragmented area. Although the area is thought to have experienced periods of relative continuity during Pleistocene glacial periods, small-scale studies suggest deep lineage divergence across low-elevation gaps. Using both DNA sequence data and microsatellite markers, we tested the hypothesis that genetic partitioning reflects observable geographic structuring across Australia’s mainland high country, in the widespread alpine grasshopper Kosciuscola tristis (Sjösted).

Results

We found broadly congruent patterns of regional structure between the DNA sequence and microsatellite datasets, corresponding to strong divergence among isolated mountain regions. Small and isolated mountains in the south of the range were particularly distinct, with well-supported divergence corresponding to climate cycles during the late Pliocene and Pleistocene. We found mixed support, however, for divergence among other mountain regions. Interestingly, within areas of largely contiguous alpine and sub-alpine habitat around Mt Kosciuszko, microsatellite data suggested significant population structure, accompanied by a strong signature of isolation-by-distance.

Conclusions

Consistent patterns of strong lineage divergence among different molecular datasets indicate genetic breaks between populations inhabiting geographically distinct mountain regions. Three primary phylogeographic groups were evident in the highly fragmented Victorian high country, while within-region structure detected with microsatellites may reflect more recent population isolation. Despite the small area of Australia’s alpine and sub-alpine habitats, their low topographic relief and lack of extensive glaciation, divergence among populations was on the same scale as that detected in much more extensive Northern hemisphere mountain systems. The processes driving divergence in the Australian mountains might therefore differ from their Northern hemisphere counterparts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0204-1) contains supplementary material, which is available to authorized users.

Integrated fossil and molecular data reveal the biogeographic diversification of the eastern Asian-eastern North American disjunct hickory genus (Carya Nutt.)

The hickory genus (Carya) contains ca. 17 species distributed in subtropical and tropical regions of eastern Asia and subtropical to temperate regions of eastern North America. Previously, the phylogenetic relationships between eastern Asian and eastern North American species of Carya were not fully confirmed even with an extensive sampling, biogeographic and diversification patterns had thus never been investigated in a phylogenetic context. We sampled 17 species of Carya and 15 species representing all other genera of the Juglandaceae as outgroups, with eight nuclear and plastid loci to reconstruct the phylogeny of Carya. The phylogenetic positions of seven extinct genera of the Juglandaceae were inferred using morphological characters and the molecular phylogeny as a backbone constraint. Divergence times within Carya were estimated with relaxed Bayesian dating. Biogeographic analyses were performed in DIVA and LAGRANGE. Diversification rates were inferred by LASER and APE packages. Our results support two major clades within Carya, corresponding to the lineages of eastern Asia and eastern North America. The split between the two disjunct clades is estimated to be 21.58 (95% HPD 11.07-35.51) Ma. Genus-level DIVA and LAGRANGE analyses incorporating both extant and extinct genera of the Juglandaceae suggested that Carya originated in North America, and migrated to Eurasia during the early Tertiary via the North Atlantic land bridge. Fragmentation of the distribution caused by global cooling in the late Tertiary resulted in the current disjunction. The diversification rate of hickories in eastern North America appeared to be higher than that in eastern Asia, which is ascribed to greater ecological opportunities, key morphological innovations, and polyploidy.

Multi-locus phylogeographic and population genetic analysis of Anolis carolinensis: historical demography of a genomic model species

The green anole (Anolis carolinensis) has been widely used as an animal model in physiology and neurobiology but has recently emerged as an important genomic model. The recent sequencing of its genome has shed new light on the evolution of vertebrate genomes and on the process that govern species diversification. Surprisingly, the patterns of genetic diversity within natural populations of this widespread and abundant North American lizard remain relatively unknown. In the present study, we use 10 novel nuclear DNA sequence loci (N = 62 to 152) and one mitochondrial locus (N = 226) to delimit green anole populations and infer their historical demography. We uncovered four evolutionarily distinct and geographically restricted lineages of green anoles using phylogenetics, bayesian clustering, and genetic distance methods. Molecular dating indicates that these lineages last shared a common ancestor ∼2 million years ago. Summary statistics and analysis of the frequency distributions of DNA polymorphisms strongly suggest range-wide expansions in population size. Using Bayesian Skyline Plots, we inferred the timing of population size expansions, which differ across lineages, and found evidence for a relatively recent and rapid westward expansion of green anoles across the Gulf Coastal Plain during the mid-Pleistocene. One surprising result is that the distribution of genetic diversity is not consistent with a latitudinal shift caused by climatic oscillations as is observed for many co-distributed taxa. This suggests that the most recent Pleistocene glacial cycles had a limited impact on the geographic distribution of the green anole at the northern limits of its range.

Ancient DNA assessment of tiger salamander population in Yellowstone National Park

Recent data indicates that blotched tiger salamanders (Ambystoma tigrinum melanostictum) in northern regions of Yellowstone National Park are declining due to climate-related habitat changes. In this study, we used ancient and modern mitochondrial haplotype diversity to model the effective size of this amphibian population through recent geological time and to assess past responses to climatic changes in the region. Using subfossils collected from a cave in northern Yellowstone, we analyzed >700 base pairs of mitochondrial sequence from 16 samples ranging in age from 100 to 3300 years old and found that all shared an identical haplotype. Although mitochondrial diversity was extremely low within the living population, we still were able to detect geographic subdivision within the local area. Using serial coalescent modelling with Bayesian priors from both modern and ancient genetic data we simulated a range of probable population sizes and mutation rates through time. Our simulations suggest that regional mitochondrial diversity has remained relatively constant even through climatic fluctuations of recent millennia.

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

BACKGROUND

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal "barcoding gene". An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

METHODS AND FINDINGS

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a "barcode gap" using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

CONCLUSIONS

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a "species".

Uncoordinated phylogeography of Borrelia burgdorferi and its tick vector, Ixodes scapularis

Vector-borne microbes necessarily co-occur with their hosts and vectors, but the degree to which they share common evolutionary or biogeographic histories remains unexplored. We examine the congruity of the evolutionary and biogeographic histories of the bacterium and vector of the Lyme disease system, the most prevalent vector-borne disease in North America. In the eastern and midwestern US, Ixodes scapularis ticks are the primary vectors of Borrelia burgdorferi, the bacterium that causes Lyme disease. Our phylogeographic and demographic analyses of the 16S mitochondrial rDNA suggest that northern I. scapularis populations originated from very few migrants from the southeastern US that expanded rapidly in the Northeast and subsequently in the Midwest after the recession of the Pleistocene ice sheets. Despite this historical gene flow, current tick migration is restricted even between proximal sites within regions. In contrast, B. burgdorferi suffers no barriers to gene flow within the northeastern and midwestern regions but shows clear interregional migration barriers. Despite the intimate association of B. burgdorferi and I. scapularis, the population structure, evolutionary history, and historical biogeography of the pathogen are all contrary to its arthropod vector. In the case of Lyme disease, movements of infected vertebrate hosts may play a larger role in the contemporary expansion and homogenization of the pathogen than the movement of tick vectors whose populations continue to bear the historical signature of climate-induced range shifts.

Pleistocene glacial refugia across the Appalachian Mountains and coastal plain in the millipede genus Narceus: evidence from population genetic, phylogeographic, and paleoclimatic data

Background

Species that are widespread throughout historically glaciated and currently non-glaciated areas provide excellent opportunities to investigate the role of Pleistocene climatic change on the distribution of North American biodiversity. Many studies indicate that northern animal populations exhibit low levels of genetic diversity over geographically widespread areas whereas southern populations exhibit relatively high levels. Recently, paleoclimatic data have been combined with niche-based distribution modeling to locate possible refugia during the Last Glacial Maximum. Using phylogeographic, population, and paleoclimatic data, we show that the distribution and mitochondrial data for the millipede genus Narceus are consistent with classical examples of Pleistocene refugia and subsequent post-glacial population expansion seen in other organismal groups.

Results

The phylogeographic structure of Narceus reveals a complex evolutionary history with signatures of multiple refugia in southeastern North America followed by two major northern expansions. Evidence for refugial populations were found in the southern Appalachian Mountains and in the coastal plain. The northern expansions appear to have radiated from two separate refugia, one from the Gulf Coastal Plain area and the other from the mid-Atlantic coastal region. Distributional models of Narceus during the Last Glacial Maximum show a dramatic reduction from the current distribution, with suitable ecological zones concentrated along the Gulf and Atlantic coastal plain. We found a strong correlation between these zones of ecological suitability inferred from our paleo-model with levels of genetic diversity derived from phylogenetic and population estimates of genetic structuring.

Conclusion

The signature of climatic change, during and after the Pleistocene, on the distribution of the millipede genus Narceus is evident in the genetic data presented. Niche-based historical distribution modeling strengthens the conclusions drawn from the genetic data and proves useful in identifying probable refugia. Such interdisciplinary biogeographic studies provide a comprehensive approach to understanding these processes that generate and maintain biodiversity as well as the framework necessary to explore questions regarding evolutionary diversification of taxa.

Controlling for the effects of history and nonequilibrium conditions in gene flow estimates in northern bullfrog (Rana catesbeiana) populations

Nonequilibrium conditions due to either allopatry followed by secondary contact or recent range expansion can confound measurements of gene flow among populations in previously glaciated regions. We determined the scale at which gene flow can be estimated among breeding aggregations of bullfrogs (Rana catesbeiana) at the northern limit of their range in Ontario, Canada, using seven highly polymorphic DNA microsatellite loci. We first identified breeding aggregations that likely share a common history, determined from the pattern of allelic richness, factorial correspondence analysis, and a previously published mtDNA phylogeography, and then tested for regional equilibrium by evaluating the association between pairwise F(ST) and geographic distance. Regional breeding aggregations in eastern Ontario separated by <100 km were determined to be at or near equilibrium. High levels of gene flow were measured using traditional F-statistics and likelihood estimates of Nm. Similarly high levels of recent migration (past one to three generations) were estimated among the breeding aggregations using nonequilibrium methods. We also show that, in many cases, breeding aggregations separated by up to tens of kilometers are not genetically distinct enough to be considered separate genetic populations. These results have important implications both for the identification of independent "populations" and in assessing the effect of scale in detecting patterns of genetic equilibrium and gene flow.