A striking lack of genetic diversity across the wide-ranging amphibian Gastrophryne carolinensis (Anura: Microhylidae)

Revisiting the comparative phylogeography of unglaciated eastern North America: 15 years of patterns and progress

In a landmark comparative phylogeographic study, “Comparative phylogeography of unglaciated eastern North America,” Soltis et al. (Molecular Ecology, 2006, 15, 4261) identified geographic discontinuities in genetic variation shared across taxa occupying unglaciated eastern North America and proposed several common biogeographical discontinuities related to past climate fluctuations and geographic barriers. Since 2006, researchers have published many phylogeographical studies and achieved many advances in genotyping and analytical techniques; however, it is unknown how this work has changed our understanding of the factors shaping the phylogeography of eastern North American taxa. We analyzed 184 phylogeographical studies of eastern North American taxa published between 2007 and 2019 to evaluate: (1) the taxonomic focus of studies and whether a previously detected taxonomic bias towards studies focused on vertebrates has changed over time, (2) the extent to which studies have adopted genotyping technologies that improve the resolution of genetic groups (i.e., NGS DNA sequencing) and analytical approaches that facilitate hypothesis‐testing (i.e., divergence time estimation and niche modeling), and (3) whether new studies support the hypothesized biogeographic discontinuities proposed by Soltis et al. (Molecular Ecology, 2006, 15, 4261) or instead support new, previously undetected discontinuities. We observed little change in taxonomic focus over time, with studies still biased toward vertebrates. Although many technological and analytical advances became available during the period, uptake was slow and they were employed in only a small proportion of studies. We found variable support for previously identified discontinuities and identified one new recurrent discontinuity. However, the limited resolution and taxonomic breadth of many studies hindered our ability to clarify the most important climatological or geographical factors affecting taxa in the region. Broadening the taxonomic focus to include more non‐vertebrate taxa, employing technologies that improve genetic resolution, and using analytical approaches that improve hypothesis testing are necessary to strengthen our inference of the forces shaping the phylogeography of eastern North America.

Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time.

Conservation genetics of two threatened frogs from the Mambilla highlands, Nigeria

Amphibians are the vertebrate group with the highest number of species threatened with extinction, and habitat loss and fragmentation are considered to be among the leading causes of their declines and extinctions. Little is known of the population biology of amphibian species inhabiting montane forests in Central and West Africa, where anthropogenic activities such as farming and cattle raising are major threats to native biodiversity. We used Amplified Fragment Length Polymorphisms (AFLPs) to assess the population genetic structure of two poorly known species, Cardioglossa schioetzi and Leptodactylodon bicolor (both in the Arthroleptidae), in and around Ngel Nyaki Forest Reserve on the Mambilla Plateau in eastern Nigeria. The landscape comprises continuous forest on steep slopes and small riparian forest fragments in a grassland matrix. While increased fragmentation is well documented for these and other forests in the mountains of Cameroon and Nigeria over the past century, there are no previous assessments of the impact of forest fragmentation on montane amphibian populations in this region. Our estimates of genetic diversity are similar across populations within each species with levels of heterozygosity values consistent with local population declines. Except for a pair of populations (C. schioetzi) we did not observe genetic differentiation between forest and riparian forest fragment populations, nor across sites within continuous forest (L. bicolor). Our results demonstrate recent gene flow between forest fragments and the adjacent protected forests and suggest that small forest corridors connecting these may lessen the genetic consequences of at least 30 years of intense and severe fragmentation in Ngel Nyaki.

High genetic diversity but low population structure in the frog Pseudopaludicola falcipes (Hensel, 1867) (Amphibia, Anura) from the Pampas of South America

Relative to South America's ecoregions, the temperate grasslands of the Pampas have been poorly studied from a phylogeographic perspective. Based on an intermediate biogeographic setting between subtropical forest (Atlantic Forest) and arid ecosystems (Chaco and Patagonia), Pampean species are expected to show unstable demographic histories due to the Quaternary climatic oscillations. Herein, we investigate the phylogenetic relatedness and phylogeographic history of Pseudopaludicola falcipes, a small and common frog that is widely distributed across the Pampean grasslands. First, we use molecular data to assess if P. falcipes represents a single or multiple, separately evolving cryptic lineages. Because P. falcipes is a small-size species (<20mm) with extensive coloration and morphological variation, we suspected that it might represent a complex of cryptic species. In addition, we expected strong genetic and geographic structuring within Pseudopaludicola falcipes due to its large geographic distribution, potentially short dispersal distances, and multiple riverine barriers. We found that P. falcipes is a single evolutionary lineage with poor geographic structuring. Furthermore, current populations of P. falcipes have a large effective population size, maintain ancestral polymorphisms, and have a complex network of gene flow. We conclude that the demographic history of P. falcipes, combined with its ecological attributes and the landscape features of the Pampas, favored a unique combination among anurans of small body size, large population size, high genetic variability, but high cohesiveness of populations over a wide geographic distribution.

Cryptic genetic diversity and complex phylogeography of the boreal North American scorpion, Paruroctonus boreus (Vaejovidae)

Diverse studies in western North America have revealed the role of topography for dynamically shaping genetic diversity within species though vicariance, dispersal and range expansion. We examined patterns of phylogeographical diversity in the widespread but poorly studied North American vaejovid scorpion, Paruroctonus boreus Girard 1854. We used mitochondrial sequence data and parsimony, likelihood, and Bayesian inference to reconstruct phylogenetic relationships across the distributional range of P. boreus, focusing on intermontane western North America. Additionally, we developed a species distribution model to predict its present and historical distributions during the Last Glacial Maximum and the Last Interglacial Maximum. Our results documented complex phylogeographic relationships within P. boreus, with multiple, well-supported crown clades that are either geographically-circumscribed or widespread and separated by short, poorly supported internodes. We also observed subtle variation in predicted habitat suitability, especially at the northern, eastern and southern edges of the predicted distributional range under past climatic conditions. The complex phylogenetic relationships of P. boreus suggests that historical isolation and expansion of populations may have occurred. Variation in the predicted distributional range over time may implicate past climatic fluctuations in generating the patterns of genetic diversity observed in P. boreus. These findings highlight both the potential for cryptic biodiversity in widespread North American scorpion species and the importance of phylogeographical studies for understanding the factors responsible for generating the biodiversity of western North America.

Low phylogeographic structure of Rhinella arunco (Anura: Bufonidae), an endemic amphibian from the Chilean Mediterranean hotspot

Background

The Mediterranean zone of central Chile (30° to 38°S) is one of the 25 diversity hotspots in the world. However, there are few phylogeographic studies which identify the factors that have influenced population diversification in the fauna of this area. In this study, we investigated the phylogeographic structure of Rhinella arunco, an anuran endemic to Mediterranean Chile, using 160 individuals from 23 localities representative of its entire distribution (32° to 38°S).

Results

The haplotype network revealed four haplogroups, three of which overlap geographically and only one of which has an exclusive geographic distribution. An analysis of molecular variance indicated that neither watershed limits nor the main rivers in the current distribution of this species have been important geographic barriers against the dispersion of individuals. Finally, the Geneland analysis showed three population units, one of which concurs with one of the haplogroups found in the haplotype network. Together, these analyses indicated a low level of phylogeographic structure for this species. On the other hand, the highest levels of intrapopulational genetic variation were concentrated in the central part of the distribution (33° to 34°S), which may indicate an effect of Pleistocene glaciations on the genetic diversity of the populations in the extreme south of its range.

Conclusions

The low phylogeographic structure observed in R. arunco is a rarely documented pattern for amphibians and contrasts with the phylogeographic studies of other vertebrates which inhabit the same zone. This result may be attributed to a series of attributes of bufonids, related mainly to water retention and their reproductive biology, which have allowed them to disperse and colonize an enormous variety of environments.

Testing the role of meander cutoff in promoting gene flow across a riverine barrier in ground skinks (Scincella lateralis)

Despite considerable attention, the long-term impact of rivers on species diversification remains uncertain. Meander loop cutoff (MLC) is one river phenomenon that may compromise a river's diversifying effects by passively transferring organisms from one side of the river to the other. However, the ability of MLC to promote gene flow across rivers has not been demonstrated empirically. Here, we test several predictions of MLC-mediated gene flow in populations of North American ground skinks (Scincella lateralis) separated by a well-established riverine barrier, the Mississippi River: 1) individuals collected from within meander cutoffs should be more closely related to individuals across the river than on the same side, 2) individuals within meander cutoffs should contain more immigrants than individuals away from meander cutoffs, 3) immigration rates estimated across the river should be highest in the direction of the cutoff event, and 4) the distribution of alleles native to one side of the river should be better predicted by the historical rather than current path of the river. To test these predictions we sampled 13 microsatellite loci and mitochondrial DNA from ground skinks collected near three ancient meander loops. These predictions were generally supported by genetic data, although support was stronger for mtDNA than for microsatellite data. Partial support for genetic divergence of samples within ancient meander loops also provides evidence for the MLC hypothesis. Although a role for MLC-mediated gene flow was supported here for ground skinks, the transient nature of river channels and morphologies may limit the long-term importance of MLC in stemming population divergence across major rivers.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Extremely low genetic diversity indicating the endangered status of Ranodon sibiricus (Amphibia: Caudata) and implications for phylogeography

BACKGROUND

The Siberian salamander (Ranodon sibiricus), distributed in geographically isolated areas of Central Asia, is an ideal alpine species for studies of conservation and phylogeography. However, there are few data regarding the genetic diversity in R. sibiricus populations.

METHODOLOGY/PRINCIPAL FINDINGS

We used two genetic markers (mtDNA and microsatellites) to survey all six populations of R. sibiricus in China. Both of the markers revealed extreme genetic uniformity among these populations. There were only three haplotypes in the mtDNA, and the overall nucleotide diversity in the mtDNA was 0.00064, ranging from 0.00000 to 0.00091 for the six populations. Although we recovered 70 sequences containing microsatellite repeats, there were only two loci that displayed polymorphism. We used the approximate Bayesian computation (ABC) method to study the demographic history of the populations. This analysis suggested that the extant populations diverged from the ancestral population approximately 120 years ago and that the historical population size was much larger than the present population size; i.e., R. sibiricus has experienced dramatic population declines.

CONCLUSION/SIGNIFICANCE

Our findings suggest that the genetic diversity in the R. sibiricus populations is the lowest among all investigated amphibians. We conclude that the isolation of R. sibiricus populations occurred recently and was a result of recent human activity and/or climatic changes. The Pleistocene glaciation oscillations may have facilitated intraspecies genetic homogeneity rather than enhanced divergence. A low genomic evolutionary rate and elevated inbreeding frequency may have also contributed to the low genetic variation observed in this species. Our findings indicate the urgency of implementing a protection plan for this endangered species.

Snapshot of DNA methylation changes associated with hybridization in Xenopus

Hybridization often results in dramatic genome reconfigurations including epigenetic changes that control gene expression. Here we survey methylation patterns of interspecific Xenopus F1 hybrids relative to parental species X. laevis and X. muelleri, using methyl-sensitive amplification polymorphisms (MSAPs). Out of a total of 546 MSAP markers, 364 were effective in elucidating the difference in methylation patterns between the hybrids and the parental species. Principal coordinate analysis of methylated fragments revealed four distinct clusters with the two parental species separate from hybrid males and females. On average, hybrids were characterized by a higher proportion (70.6%) of methylated fragments compared with the parental species (64.5%), and this difference was consistent with previously observed disruptions of hybrid transcriptomes. The proportion of methylated fragments did not correlate with variation in genome size, as measured with flow cytometry. The levels of methylation in sterile hybrid males (73.8%) were higher than in fertile hybrid females (68.6%), but this difference was not statistically significant. A total of 76 methylated fragments (20.9%) were hybrid-unique, presumably originating from methylation alterations in hybrid genomes.

Molecular systematics of the Middle American genus Hypopachus (Anura: Microhylidae)

We present the first phylogenetic study on the widespread Middle American microhylid frog genus Hypopachus. Partial sequences of mitochondrial (12S and 16S ribosomal RNA) and nuclear (rhodopsin) genes (1275 bp total) were analyzed from 43 samples of Hypopachus, three currently recognized species of Gastrophryne, and seven arthroleptid, brevicipitid and microhylid outgroup taxa. Maximum parsimony (PAUP), maximum likelihood (RAxML) and Bayesian inference (MrBayes) optimality criteria were used for phylogenetic analyses, and BEAST was used to estimate divergence dates of major clades. Population-level analyses were conducted with the programs NETWORK and Arlequin. Results confirm the placement of Hypopachus and Gastrophryne as sister taxa, but the latter genus was strongly supported as paraphyletic. The African phrynomerine genus Phrynomantis was recovered as the sister taxon to a monophyletic Chiasmocleis, rendering our well-supported clade of gastrophrynines paraphyletic. Hypopachus barberi was supported as a disjunctly distributed highland species, and we recovered a basal split in lowland populations of Hypopachus variolosus from the Pacific versant of Mexico and elsewhere in the Mesoamerican lowlands. Dating analyses from BEAST estimate speciation within the genus Hypopachus occurred in the late Miocene/early Pliocene for most clades. Previous studies have not found bioacoustic or morphological differences among these lowland clades, and our molecular data support the continued recognition of two species in the genus Hypopachus.

Phylogeographic analysis and environmental niche modeling of the plain-bellied watersnake (Nerodia erythrogaster) reveals low levels of genetic and ecological differentiation

Species that exhibit geographically defined phenotypic variation traditionally have been divided into subspecies. Subspecies based on phenotypic features may not comprise monophyletic groups due to selection, gene flow, and/or convergent evolution. In many taxonomic groups the number of species once designated as widespread is dwindling rapidly, and many workers reject the concept of subspecies altogether. We tested whether currently recognized subspecies in the plain-bellied watersnake Nerodia erythrogaster are concordant with relationships based on mitochondrial markers, and whether it represents a single widespread species. The range of this taxon spans multiple potential biogeographic barriers (especially the Mississippi and Apalachicola Rivers) that correspond with lineage breaks in many species, including other snakes. We sequenced three mitochondrial genes (NADH-II, Cyt-b, Cox-I) from 156 geo-referenced specimens and developed ecological niche models using Maxent and spatially explicit climate data to examine historical and ecological factors affecting variation in N. erythrogaster across its range. Overall, we found little support for the recognized subspecies as either independent evolutionary lineages or geographically circumscribed units and conclude that although some genetic and niche differentiation has occurred, most populations assigned to N. erythrogaster appear to represent a single, widespread species. However, additional sampling and application of nuclear markers are necessary to clarify the status of the easternmost populations.

Marked genetic structuring and extreme dispersal limitation in the Pyrenean brook newt Calotriton asper (Amphibia: Salamandridae) revealed by genome-wide AFLP but not mtDNA

Direct estimation of dispersal rates at large geographic scales can be technically and logistically challenging, especially in small animals of low vagility like amphibians. The use of molecular markers to reveal patterns of genetic structure provides an indirect way to infer dispersal rates and patterns of recent and historical gene flow among populations. Here, we use mitochondrial DNA (mtDNA) sequence data and genome-wide amplified fragment length polymorphism markers to examine population structure in the Pyrenean brook newt (Calotriton asper) across four main drainages in the French Pyrenees. mtDNA sequence data (2040 bp) revealed three phylogroups shallowly differentiated and with low genetic diversity. In sharp contrast, variation in 382 amplified fragment length polymorphism loci was high and revealed a clear pattern of isolation by distance consistent with long-term restriction of gene flow at three spatial scales: (i) among all four main drainages, (ii) between sites within drainages, and (iii) even between adjacent populations separated by less than 4 km. The high pairwise F(ST) values between localities across numerous loci, together with the high frequency of fixed alleles in several populations, suggests a combination of marked geographic isolation, small population sizes and very limited dispersal in C. asper. The contrasting lack of variation detected in mtDNA sequence data is intriguing and underscores the importance of multilocus approaches to detect true patterns of gene flow in natural populations of amphibians.