Cryptic Neogene vicariance and Quaternary dispersal of the red-spotted toad (Bufo punctatus): insights on the evolution of North American warm desert biotas

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".

Hybridization and massive mtDNA unidirectional introgression between the closely related Neotropical toads Rhinella marina and R. schneideri inferred from mtDNA and nuclear markers

Background

The classical perspective that interspecific hybridization in animals is rare has been changing due to a growing list of empirical examples showing the occurrence of gene flow between closely related species. Using sequence data from cyt b mitochondrial gene and three intron nuclear genes (RPL9, c-myc, and RPL3) we investigated patterns of nucleotide polymorphism and divergence between two closely related toad species R. marina and R. schneideri. By comparing levels of differentiation at nuclear and mtDNA levels we were able to describe patterns of introgression and infer the history of hybridization between these species.

Results

All nuclear loci are essentially concordant in revealing two well differentiated groups of haplotypes, corresponding to the morphologically-defined species R. marina and R. schneideri. Mitochondrial DNA analysis also revealed two well-differentiated groups of haplotypes but, in stark contrast with the nuclear genealogies, all R. schneideri sequences are clustered with sequences of R. marina from the right Amazon bank (RAB), while R. marina sequences from the left Amazon bank (LAB) are monophyletic. An Isolation-with-Migration (IM) analysis using nuclear data showed that R. marina and R. schneideri diverged at ≈ 1.69 Myr (early Pleistocene), while R. marina populations from LAB and RAB diverged at ≈ 0.33 Myr (middle Pleistocene). This time of divergence is not consistent with the split between LAB and RAB populations obtained with mtDNA data (≈ 1.59 Myr), which is notably similar to the estimate obtained with nuclear genes between R. marina and R. schneideri. Coalescent simulations of mtDNA phylogeny under the speciation history inferred from nuclear genes rejected the hypothesis of incomplete lineage sorting to explain the conflicting signal between mtDNA and nuclear-based phylogenies.

Conclusions

The cytonuclear discordance seems to reflect the occurrence of interspecific hybridization between these two closely related toad species. Overall, our results suggest a phenomenon of extensive mtDNA unidirectional introgression from the previously occurring R. schneideri into the invading R. marina. We hypothesize that climatic-induced range shifts during the Pleistocene/Holocene may have played an important role in the observed patterns of introgression.

Elucidation of cryptic diversity in a widespread nearctic treefrog reveals episodes of mitochondrial gene capture as frogs diversified across a dynamic landscape

We investigate the evolutionary history of the wide-ranging Nearctic treefrog Hyla arenicolor through the integration of extensive range-wide sampling, phylogenetic analyses of multilocus genetic data, and divergence dating. Previous phylogeographic studies of this frog documented a potential signature of introgressive hybridization from an ecologically and morphologically divergent sister species. Based on our Bayesian phylogenetic analyses of mitochondrial DNA, we inferred strong phylogeographic structure in H. arenicolor as indicated by seven well-supported clades, five of which correspond to well-defined biogeographic regions. Clades from the Balsas Basin and southwestern Central Mexican Plateau in Mexico, and the Grand Canyon of Arizona, group with the morphologically, behaviorally, and ecologically divergent mountain treefrogs in the H. eximia group, rendering H. arenicolor as paraphyletic. The phylogenetic position of at least two of these three H. arenicolor clades within the H. eximia group, however, is most likely the result of several episodes of introgressive hybridization and subsequent mitochondrial gene capture separated in time and space, as supported by evidence from the nuclear genes. Hyla arenicolor from the Balsas Basin appear to be deeply divergent from other H. arenicolor and represent a distinctly different species. Results suggests that introgressive hybridization events, both ancient and contemporary, coupled with late Neogene vicariance and Pleistocene climate-driven range shifts, have all played a role in the historical diversification of H. arenicolor.

Evolution of the nocturnal Nearctic Sphaeropthalminae velvet ants (Hymenoptera: Mutillidae) driven by Neogene orogeny and Pleistocene glaciation

The influence of historical climatic and geological changes on patterns of species diversification was investigated in a widely distributed group of North American nocturnal mutillids (Hymenoptera: Mutillidae: Sphaeropthalminae), with particular focus on Pleistocene glacial cycles and earlier patterns of Neogene mountain building. We collected molecular data from two nuclear intergenic regions (internal transcribed spacer regions 1 and 2; approximately 2600 bp in total) to produce the first phylogeny of nocturnal Nearctic mutillids. Bayesian inference of the combined data returned a well-resolved tree with posterior probabilities of over 95% for most nodes. This tree suggested the monophyly of the nocturnal, primarily Nearctic, Sphaeropthalminae genera, but the paraphyly of the three largest genera (Odontophotopsis, Photomorphus and Sphaeropthalma). Dates of species divergences were obtained using r8s (PL and NPRS) and BEAST with the date of Dominican amber set at three different dates (15 Ma, 20 Ma, and 45 Ma) to account for uncertainty in the fossil age. The derived dates ranged from the Pleistocene to the Middle Miocene, suggesting that both recent Pleistocene glaciation cycles and older orogenic events, albeit to a somewhat greater extent, were both causes of major diversification in western North America. Examination of other phylogeographical studies using North American desert taxa indicated that diversification patterns are explained by either mountain building or Pleistocene climate change, depending on the taxa in question.

Conservation phylogenetics of helodermatid lizards using multiple molecular markers and a supertree approach

We analyzed both mitochondrial (mt-) and nuclear (n) DNAs in a conservation phylogenetic framework to examine deep and shallow histories of the Beaded Lizard (Heloderma horridum) and Gila Monster (H. suspectum) throughout their geographic ranges in North and Central America. Both mtDNA and intron markers clearly partitioned each species. One intron and mtDNA further subdivided H. horridum into its four recognized subspecies (H. n. alvarezi, charlesbogerti,exasperatum, and horridum). However, the two subspecies of H. suspectum (H. s. suspectum and H. s. cinctum) were undefined. A supertree approach sustained these relationships. Overall, the Helodermatidae is reaffirmed as an ancient and conserved group. Its most recent common ancestor (MRCA) was Lower Eocene [35.4 million years ago (mya)], with a approximately 25 my period of stasis before the MRCA of H. horridum diversified in Lower Miocene. Another approximately 5 my passed before H. h. exasperatum and H. h. horridum diverged, followed by approximately 1.5 my before H. h. alvarezi and H. h. charlesbogerti separated. Heloderma suspectum reflects an even longer period of stasis (approximately 30 my) before diversifying from its MRCA. Both H. suspectum (México) and H. h. alvarezi also revealed evidence of historic range expansion following a recent bottleneck. Our conservation phylogenetic approach emphasizes the origin and diversification of this group, yields information on the manner by which past environmental variance may have impacted its populations and, in turn, allows us to disentangle historic from contemporary impacts that might threaten its long-term persistence. The value of helodermatid conservation resides in natural services and medicinal products, particularly venom constituents, and these are only now being realized.

Phylogeny, divergence times and species limits of spiny lizards (Sceloporus magister species group) in western North American deserts and Baja California

The broad distribution of the Sceloporus magister species group (squamata: phrynosomatidae) throughout western North America provides an appropriate model for testing biogeographical hypotheses explaining the timing and origins of diversity across mainland deserts and the Baja California Peninsula. We inferred concordant phylogenetic trees describing the higher-level relationships within the magister group using 1.6 kb of mitochondrial DNA (mtDNA) and 1.7 kb of nuclear DNA data. These data provide strong support for the parallel divergence of lineages endemic to the Baja California Peninsula (S. zosteromus and the orcutti complex) in the form of two sequential divergence events at the base of the magister group phylogeny. A relaxed phylogenetic analysis of the mtDNA data using one fossil and one biogeographical constraint provides a chronology of these divergence events and evidence that further diversification within the Baja California clades occurred simultaneously, although patterns of geographical variation and speciation between clades differ. We resolved four major phylogeographical clades within S. magister that (i) provide a novel phylogenetic placement of the Chihuahuan Desert populations sister to the Mojave Desert; (ii) illustrate a mixed history for the Colorado Plateau that includes Mojave and Sonoran Desert components; and (iii) identify an area of overlap between the Mojave and Sonoran Desert clades near Yuma, Arizona. Estimates of bidirectional migration rates among populations of S. magister using four nuclear loci support strong asymmetries in gene flow among the major mtDNA clades. Based on the nonexclusivity of mtDNA haplotypes, nuclear gene flow among populations and wide zones of phenotypic intergradation, S. magister appears to represent a single geographically variable and widespread species.

Multi-locus DNA sequence data reveal a history of deep cryptic vicariance and habitat-driven convergence in the desert night lizardXantusia vigilisspecies complex (Squamata: Xantusiidae)

The lizard genus Xantusia of southwestern North America has received recent attention in relation to delimiting species. Using more than 500 lizards from 156 localities, we further test hypothesized species boundaries and clarify phylogeographical patterns, particularly in regions of potential secondary contact. We sequenced the entire mitochondrial cytochrome b gene for every lizard in the study, plus a second mitochondrial DNA (mtDNA) region and two nuclear introns for subsets of the total sample. Phylogenetic analyses of the mtDNA recover a well-resolved, novel hypothesis for species in the Xantusia vigilis complex. The nuclear DNA (nDNA) data provide independent support for the recognition of X. arizonae, X. bezyi and X. wigginsi. Differences between the respective mtDNA and nDNA topologies result from either the effects of lineage sorting or ancient introgression. Nuclear data confirm the inference that some populations of X. vigilis in northwestern Arizona converged on rock-crevice-dwelling morphology and are not X. arizonae with an introgressed X. vigilis mtDNA genome. The historical independence of ancient cryptic lineages of Xantusia in southern California is also corroborated, though limited introgression is detected. Our proposed biogeographical scenario indicates that diversification of this group was driven by vicariance beginning in the late Miocene. Additionally, Pleistocene climatical changes influenced Xantusia distribution, and the now inhospitable Colorado Desert previously supported night lizard presence. The current taxonomy of the group likely underestimates species diversity within the group, and our results collectively show that while convergence on the rock-crevice-dwelling morphology is one hallmark of Xantusia evolution, morphological stasis is paradoxically another.

Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Nearctic-Neotropical transition

The Western Lyresnake (Trimorphodon biscutatus) is a widespread, polytypic taxon inhabiting arid regions from the warm deserts of the southwestern United States southward along the Pacific versant of Mexico to the tropical deciduous forests of Mesoamerica. This broadly distributed species provides a unique opportunity to evaluate a priori biogeographical hypotheses spanning two major distinct biogeographical realms (the Nearctic and Neotropical) that are usually treated separately in phylogeographical analyses. I investigated the phylogeography of T. biscutatus using maximum likelihood and Bayesian phylogenetic analysis of mitochondrial DNA (mtDNA) from across this species' range. Phylogenetic analyses recovered five well-supported clades whose boundaries are concordant with existing geographical barriers, a pattern consistent with a model of vicariant allopatric divergence. Assuming a vicariance model, divergence times between mitochondrial lineages were estimated using Bayesian relaxed molecular clock methods calibrated using geological information from putative vicariant events. Divergence time point estimates were bounded by broad confidence intervals, and thus these highly conservative estimates should be considered tentative hypotheses at best. Comparison of mtDNA lineages and taxa traditionally recognized as subspecies based on morphology suggest this taxon is comprised of multiple independent lineages at various stages of divergence, ranging from putative secondary contact and hybridization to sympatry of 'subspecies'.

Evolution of rattlesnakes (Viperidae; Crotalus) in the warm deserts of western North America shaped by Neogene vicariance and Quaternary climate change

During Pleistocene, the Laurentide ice sheet rearranged and diversified biotic distributions in eastern North America, yet had minimal physical impact in western North America where lineage diversification is instead hypothesized to result from climatic changes. If Pleistocene climatic fluctuations impacted desert species, the latter would reflect patterns of restricted gene flow concomitant with indications of demographic bottlenecks. Accordingly, molecular evidence for refugia should be present within these distributions and for subsequent range expansions as conditions improved. We sought answers to these questions by evaluating mitochondrial DNA (mtDNA) sequences from four species of rattlesnakes [Crotalus mitchellii (speckled rattlesnake), Crotalus cerastes (sidewinder), Crotalus tigris (tiger rattlesnake), Crotalus ruber (red diamond rattlesnake)] with distributions restricted to desert regions of southwestern North America. We inferred relationships using parsimony and maximum likelihood, tested intraspecific clades for population expansions, applied an isolation-with-migration model to determine bi-directional migration rates (m) among regions, and inferred divergence times for species and clades by applying a semiparametric penalized likelihood approach to our molecular data. Evidence for significant range expansion was present in two of eight regions in two species (Crotalus mitchellii pyrrhus, C. tigris region north). Two species (C. cerastes, C. mitchellii) showed a distribution concomitant with northward displacement of Baja California from mainland México, followed by vicariant separation into subclades. Effects of Pleistocene climate fluctuations were found in the distributions of all four species. Three regional diversification patterns were identified: (i) shallow genetic diversity that resulted from Pleistocene climatic events (C. tigris, C. ruber); (ii) deep Pleistocene divisions indicating allopatric segregation of subclades within refugia (C. mitchellii, C. cerastes); and (iii) lineage diversifications that extended to Pliocene or Late Miocene (C. mitchellii, C. cerastes). Clade-diversifying and clade-constraining effects impacted the four species of rattlesnakes unequally. We found relatively high levels of molecular diversification in the two most broadly distributed species (C. mitchellii, C. cerastes), and lower levels of genetic diversification in the two species (C. tigris, C. ruber) whose ranges are relatively more restricted. Furthermore, in several cases, the distributions of subspecies were not congruent with our molecular information. We suggest regional conservation perspectives for southwestern deserts cannot rely upon subspecies as biodiversity surrogates, but must instead employ a molecular and deep historical perspective as a primary mechanism to frame biodiversity reserves within this region.

Phylogeography of Pseudacris regilla (Anura: Hylidae) in western North America, with a proposal for a new taxonomic rearrangement

The Baja California populations of Pseudacris regilla, a widespread species in Western North America ranging from British Columbia to southern Baja California, are characterized by extensive geographic fragmentation. We performed phylogeographic and historical demographic analyses on 609 bp of the cytochrome b mitochondrial gene of 110 individuals representing 28 populations to determine the relative influences of current and historical processes in shaping the present distribution of genetic diversity on the Baja California Peninsula. Haplotypes from this area were nested in a clade with three well-differentiated groups. Two of these groups are from Baja California Sur and another is from California and Baja California. The estimated date for the split of these groups, between 0.9-1 Ma, fits with previously proposed hypotheses of vicariance due to different transpeninsular seaways, although successive population fragmentation and expansion due to climatic oscillations during Pleistocene glaciations cannot be discarded. Historical demographic analyses detected signs of past population expansions, especially in the southernmost group. With respect to populations north of this region, two older clades were identified, one with haplotypes mainly distributed in central California, and the other corresponding to the northern half of the species range, in what apparently is a recurrent pattern in the Pacific coast of North America. Based on the concordance between mt-DNA and available allozyme data indicating that these species have a long independent evolutionary history, we propose to consider the three major clades as distinct species: P. regilla, P. pacifica, and P. hypochondriaca.