Phylogeography and conservation genetics of the Columbia spotted frog (Rana luteiventris; Amphibia, Ranidae)

Low-level pathogen infection and geographic location correlate with the skin microbiomes of Columbia spotted frogs (Rana luteiventris) in a montane landscape

The skin microbiome of amphibians can influence host susceptibility towards the fungal pathogen Batrachochytrium dendrobatidis (Bd), while simultaneously having the potential to be altered by Bd. Severe Bd infections are known to alter the amphibian skin microbiome; however, little is known about microbiome interactions in amphibians with low infection intensity. In addition to disease dynamics, environmental factors may influence the microbiome. To test for patterns in bacterial diversity based on pathogen infection and environmental factors, 399 Columbia spotted frogs (Rana luteiventris) were sampled throughout northern Idaho and northeastern Washington across two years. Bd prevalence and intensity were measured in 376 frogs, revealing a prevalence of 69%, but generally low infection intensity (Mean = 127 Bd zoospore equivalents among infected frogs). Skin bacterial communities were characterized in 92 frogs using 16S rRNA gene amplicon sequencing. Our results indicated correlations of decreasing Shannon diversity and evenness as infection intensity increased. Latitude was correlated with bacterial richness and Faith's Phylogenetic Diversity measures, indicating increased diversity in northern locations. Beta diversity (UniFrac) analyses revealed that skin microbiomes were distinct between infected and uninfected frogs, and infection intensity had a significant effect on microbiome composition. Site explained the majority of microbiome variation (weighted UniFrac: 57.5%), suggesting a combination of local habitat conditions explain variation, as only small proportions of variation could be explained by year, month, temperature, elevation, and latitude individually. Bacterial genera with potential for Bd-inhibitory properties were found with differential relative abundance in infected and uninfected frogs, with higher Stenotrophomonas and lower Pseudomonas relative abundance observed in infected frogs. Further study may indicate if Bd inhibition by members of the skin microbiome is an influence behind the low infection intensities observed and whether low Bd infection intensities are capable of altering skin microbiome composition.

Assay validation of saliva glucocorticoids in Columbia spotted frogs and effects of handling and marking

Non-invasive methods are important to the field of conservation physiology to reduce negative effects on organisms being studied. Glucocorticoid (GC) hormones are often used to assess health of individuals, but collection methods can be invasive. Many amphibians are imperiled worldwide, and saliva is a non- or semi-invasive matrix to measure GCs that has been partially validated for only four amphibian species. Validation ensures that assays are reliable and can detect changes in saliva corticosterone (sCORT) after exposure to stressors, but it is also necessary to ensure sCORT concentrations are correlated with plasma concentrations. To help validate the use of saliva in assessing CORT responses in amphibians, we captured uniquely marked Columbia spotted frogs (Rana luteiventris) on sequential days and collected baseline and stress-induced (after handling) samples. For a subset of individuals, we collected and quantified CORT in both saliva and blood samples, which have not been compared for amphibians. We tested several aspects of CORT responses and, by collecting across separate days, measured repeatability of CORT responses across days. We also evaluated whether methods common to amphibian conservation, such as handling alone or handling, clipping a toe and tagging elevated sCORT. Similar to previous studies, we show that sCORT is reliable concerning parallelism, recovery, precision and sensitivity. sCORT was weakly correlated with plasma CORT (R2 = 0.21), and we detected elevations in sCORT after handling, demonstrating biological validation. Toe clipping and tagging did not increase sCORT over handling alone, but repeated handling elevated sCORT for ~72 hours. However, sCORT responses were highly variable and repeatability was low within individuals and among capture sessions, contrary to previous studies with urinary and waterborne CORT. sCORT is a semi-invasive and rapid technique that could be useful to assess effects of anthropogenic change and conservation efforts, but will require careful study design and future validation.

Effects of changing climate on aquatic habitat and connectivity for remnant populations of a wide-ranging frog species in an arid landscape

Amphibian species persisting in isolated streams and wetlands in desert environments can be susceptible to low connectivity, genetic isolation, and climate changes. We evaluated the past (1900–1930), recent (1981–2010), and future (2071–2100) climate suitability of the arid Great Basin (USA) for the Columbia spotted frog (Rana luteiventris) and assessed whether changes in surface water may affect connectivity for remaining populations. We developed a predictive model of current climate suitability and used it to predict the historic and future distribution of suitable climates. We then modeled changes in surface water availability at each time period. Finally, we quantified connectivity among existing populations on the basis of hydrology and correlated it with interpopulation genetic distance. We found that the area of the Great Basin with suitable climate conditions has declined by approximately 49% over the last century and will likely continue to decline under future climate scenarios. Climate conditions at currently occupied locations have been relatively stable over the last century, which may explain persistence at these sites. However, future climates at these currently occupied locations are predicted to become warmer throughout the year and drier during the frog's activity period (May – September). Fall and winter precipitation may increase, but as rain instead of snow. Earlier runoff and lower summer base flows may reduce connectivity between neighboring populations, which is already limited. Many of these changes could have negative effects on remaining populations over the next 50–80 years, but milder winters, longer growing seasons, and wetter falls might positively affect survival and dispersal. Collectively, however, seasonal shifts in temperature, precipitation, and stream flow patterns could reduce habitat suitability and connectivity for frogs and possibly other aquatic species inhabiting streams in this arid region.

Persistence at distributional edges: Columbia spotted frog habitat in the arid Great Basin, USA

A common challenge in the conservation of broadly distributed, yet imperiled species is understanding which factors facilitate persistence at distributional edges, locations where populations are often vulnerable to extirpation due to changes in climate, land use, or distributions of other species. For Columbia spotted frogs (Rana luteiventris) in the Great Basin (USA), a genetically distinct population segment of conservation concern, we approached this problem by examining (1) landscape-scale habitat availability and distribution, (2) water body-scale habitat associations, and (3) resource management-identified threats to persistence. We found that areas with perennial aquatic habitat and suitable climate are extremely limited in the southern portion of the species’ range. Within these suitable areas, native and non-native predators (trout and American bullfrogs [Lithobates catesbeianus]) are widespread and may further limit habitat availability in upper- and lower-elevation areas, respectively. At the water body scale, spotted frog occupancy was associated with deeper sites containing abundant emergent vegetation and nontrout fish species. Streams with American beaver (Castor canadensis) frequently had these structural characteristics and were significantly more likely to be occupied than ponds, lakes, streams without beaver, or streams with inactive beaver ponds, highlighting the importance of active manipulation of stream environments by beaver. Native and non-native trout reduced the likelihood of spotted frog occupancy, especially where emergent vegetation cover was sparse. Intensive livestock grazing, low aquatic connectivity, and ephemeral hydroperiods were also negatively associated with spotted frog occupancy. We conclude that persistence of this species at the arid end of its range has been largely facilitated by habitat stability (i.e., permanent hydroperiod), connectivity, predator-free refugia, and a commensalistic interaction with an ecosystem engineer. Beaver-induced changes to habitat quality, stability, and connectivity may increase spotted frog population resistance and resilience to seasonal drought, grazing, non-native predators, and climate change, factors which threaten local or regional persistence.

Genetic structure and diversity of the endangered growling grass frog in a rapidly urbanizing region

Two pervasive and fundamental impacts of urbanization are the loss and fragmentation of natural habitats. From a genetic perspective, these impacts manifest as reduced genetic diversity and ultimately reduced genetic viability. The growling grass frog (Litoria raniformis) is listed as vulnerable to extinction in Australia, and endangered in the state of Victoria. Remaining populations of this species in and around the city of Melbourne are threatened by habitat loss, degradation and fragmentation due to urban expansion. We used mitochondrial DNA (mtDNA) and microsatellites to study the genetic structure and diversity of L. raniformis across Melbourne's urban fringe, and also screened four nuclear gene regions (POMC, RAG-1, Rhod and CRYBA1). The mtDNA and nuclear DNA sequences revealed low levels of genetic diversity throughout remnant populations of L. raniformis. However, one of the four regions studied, Cardinia, exhibited relatively high genetic diversity and several unique haplotypes, suggesting this region should be recognized as a separate Management Unit. We discuss the implications of these results for the conservation of L. raniformis in urbanizing landscapes, particularly the potential risks and benefits of translocation, which remains a contentious management approach for this species.

Morphometric and genetic differentiation among populations of Eupemphix nattereri (Amphibia, Anura, Leiuperidae) from central Brazil

To assess genetic structure and phenotypic diversity of Eupemphix nattereri Steindachner, 1863, morphometric and molecular analyses were carried out for nine populations from the State of Goiás. A total of 11 morphometric traits were evaluated and genetic information was estimated using RAPD markers. Genetic and phenotypic distances were determined as a function of geographical origin. Correlation among genetic, morphometric, micro, and macroenviromental were analyzed by the Mantel test. Genetic data indicated high levels of genetic diversity (Φst= 0.3) among the nine populations. Mantel tests did not reveal a significant positive correlation between genetic and geographical distances, indicating that locally geographical populations were not genetically similar, even in distances smaller than 50 km. Discriminant analysis on 11 morphometric measurements showed a high divergence among the nine populations. However, a marginally significant correlation (P=0.08) between genetic and morphometric distances was found. The observed correlation was not causal in terms of the relationship between phenotype and genotype, but indicated common spatial structures. Thus, our results suggest that isolation-by-distance processes may explain population divergence in Eupemphix nattereri.

The postglacial recolonization of Northern Europe by Rana arvalis as revealed by microsatellite and mitochondrial DNA analyses

The postglacial history of the moor frog (Rana arvalis) in Northern Europe was investigated with the aid of eight variable microsatellite loci and a 661 bp sequence of the mitochondrial cytochrome b gene. A division between eastern and western mitochondrial lineages was discovered, supporting two recolonization routes to Fennoscandia since the last glacial maximum. This result was corroborated by the microsatellite data, which revealed a contact zone between the two lineages in Northern Sweden. These findings add to the increasing evidence that an intraspecific genetic biodiversity founded on the existence of eastern and western clades is a common element in Fennoscandian fauna and flora.

Postglacial colonization of the Tibetan plateau inferred from the matrilineal genetic structure of the endemic red-necked snow finch, Pyrgilauda ruficollis

Most phylogeographical studies of postglacial colonization focus on high latitude locations in the Northern Hemisphere. Here, we studied the phylogeographical structure of the red-necked snow finch Pyrgilauda ruficollis, an endemic species of the Tibetan plateau. We analysed 879 base pairs (bp) of the mitochondrial cytochrome b gene and 529 bp of the control region in 41 birds from four regional groups separated by mountain ranges. We detected 34 haplotypes, 31 of which occurred in a single individual and only three of which were shared among sampling sites within regional groups or among regional groups. Haplotype diversity was high (h = 0.94); nucleotide diversity was low (eth = 0.00415) and genetic differentiation was virtually non-existent. Analyses of mismatch distributions and geographically nested clades yielded results consistent with contiguous range expansion, and the expansion times were estimated as 0.07-0.19 million years ago (Ma). Our results suggest that P. ruficollis colonized the Tibetan plateau after the extensive glacial period (0.5-0.175 Ma), expanding from the eastern margin towards the inner plateau. Thus, in contrast to many of the post-glacial phylogeographical structures known at high latitudes, this colonization occurred without matrilineal population structuring. This might be due to the short glacial cycles typical of the Tibetan plateau, adaptation of P. ruficollis to cold conditions, or refugia and colonized habitat being semicontinuous and thus promoting population mixing.

Population structure of Columbia spotted frogs (Rana luteiventris) is strongly affected by the landscape

Landscape features such as mountains, rivers, and ecological gradients may strongly affect patterns of dispersal and gene flow among populations and thereby shape population dynamics and evolutionary trajectories. The landscape may have a particularly strong effect on patterns of dispersal and gene flow in amphibians because amphibians are thought to have poor dispersal abilities. We examined genetic variation at six microsatellite loci in Columbia spotted frogs (Rana luteiventris) from 28 breeding ponds in western Montana and Idaho, USA, in order to investigate the effects of landscape structure on patterns of gene flow. We were particularly interested in addressing three questions: (i) do ridges act as barriers to gene flow? (ii) is gene flow restricted between low and high elevation ponds? (iii) does a pond equal a 'randomly mating population' (a deme)? We found that mountain ridges and elevational differences were associated with increased genetic differentiation among sites, suggesting that gene flow is restricted by ridges and elevation in this species. We also found that populations of Columbia spotted frogs generally include more than a single pond except for very isolated ponds. There was also evidence for surprisingly high levels of gene flow among low elevation sites separated by large distances. Moreover, genetic variation within populations was strongly negatively correlated with elevation, suggesting effective population sizes are much smaller at high elevation than at low elevation. Our results show that landscape features have a profound effect on patterns of genetic variation in Columbia spotted frogs.

Genetic diversity and divergence among freshwater mussel (Anodonta) populations in the Bonneville Basin of Utah

Populations of the freshwater mussel genus Anodonta appear to be in a state of rapid decline in western North America, following a trend that unfortunately seems to be prevalent among these animals (Mollusca: Unionoida). Here we describe the patterns of molecular divergence and diversity among Anodonta populations in the Bonneville Basin, a large sub-basin of the Great Basin in western North America. Using amplified fragment length polymorphism (AFLP) analysis, we found a striking lack of nuclear diversity within some of these populations, along with a high degree of structuring among populations (FST = 0.61), suggesting post-Pleistocene isolation, due either to a long-term loss of hydrologic connectivity among populations or to more recent fish introductions. We also found evidence of recent hybridization in one of these populations, possibly mediated by fish-stocking practices. Using mitochondrial sequence data, we compared the Bonneville Basin populations to Anodonta in several other drainages in western North America. We found a general lack of resolution in these phylogenetic reconstructions, although there was a tendency for the Bonneville Basin Anodonta (tentatively A. californiensis) to cluster with A. oregonensis from the adjacent Lahontan Basin in Nevada. We recommend further investigation of anthropogenic factors that may be contributing to the decline of western Anodonta and a broad-scale analysis and synthesis of genetic and morphological variation among Anodonta in western North America.

Mitochondrial phylogeography of the moor frog, Rana arvalis

The moor frog Rana arvalis is a lowland species with a broad Eurasiatic distribution, from arctic tundra through forest to the steppe zone. Its present-day range suggests that glacial refugia of this frog were located outside southern European peninsulas. We studied the species-wide phylogeographical pattern using sequence variation in a 682 base pairs fragment of mtDNA cytochrome b gene; 223 individuals from 73 localities were analysed. Two main clades, A and B, differing by c. 3.6% sequence divergence were detected. The A clade is further subdivided into two subclades, AI and AII differing by 1.0%. All three lineages are present in the Carpathian Basin (CB), whereas the rest of the species range, including huge expanses of Eurasian lowlands, are inhabited solely by the AI lineage. We infer that AII and B lineages survived several glacial cycles in the CB but did not expand, at least in the present interglacial, to the north. The geographical distribution and genealogical relationships between haplotypes from the AI lineage indicate that this group had two glacial refugia, one located in the eastern part of the CB and the other probably in southern Russia. Populations from both refugia contributed to the colonization of the western part of the range, whereas the eastern part was colonized from the eastern refugium only. The effective population size as evidenced by theta(ML) is an order of magnitude higher in the AI lineage than in the AII and B lineages. Demographic expansion was detected in all three lineages.

Huge populations and old species of Costa Rican and Panamanian dirt frogs inferred from mitochondrial and nuclear gene sequences

Molecular genetic data were used to investigate population sizes and ages of Eleutherodactylus (Anura: Leptodactylidae), a species-rich group of small leaf-litter frogs endemic to Central America. Population genetic structure and divergence was investigated for four closely related species surveyed across nine localities in Costa Rica and Panama. DNA sequence data were collected from a mitochondrial gene (ND2) and a nuclear gene (c-myc). Phylogenetic analyses yielded concordant results between loci, with reciprocal monophyly of mitochondrial DNA haplotypes for all species and of c-myc haplotypes for three of the four species. Estimates of genetic differentiation among populations (FST) based upon mitochondrial data were always higher than nuclear-based FST estimates, even after correcting for the expected fourfold lower effective population size (Ne) of the mitochondrial genome. Comparing within-population variation and the relative mutation rates of the two genes revealed that the Ne of the mitochondrial genome was 15-fold lower than the estimate of the nuclear genome based on c-myc. Nuclear FST estimates were approximately 0 for the most proximal pairs of populations, but ranged from 0.5 to 1.0 for all other pairs, even within the same nominal species. The nuclear locus yielded estimates of Ne within localities on the order of 105. This value is two to three orders of magnitude larger than any previous Ne estimate from frogs, but is nonetheless consistent with published demographic data. Applying a molecular clock model suggested that morphologically indistinguishable populations within one species may be 107 years old. These results demonstrate that even a geologically young and dynamic region of the tropics can support very old lineages that harbour great levels of genetic diversity within populations. The association of high nucleotide diversity within populations, large divergence between populations, and high species diversity is also discussed in light of neutral community models.

Palaeoclimatic changes explain Anatolian mountain frog evolution: a test for alternating vicariance and dispersal events

Holarctic biodiversity has been influenced by climatic fluctuations since the Pliocene. Asia Minor was one of the major corridors for postglacial invasions in the Palearctic. Today this area is characterized by an extraordinarily rich fauna with close affiliation to European, Asian and Indo-African biota. However, exact scenarios of range expansion and contraction are lacking. Using a phylogeographical approach we (i). identify monophyletic lineages among Anatolian mountain frogs and (ii). derive a spatio-temporal hypothesis for the invasion process in Anatolia. We sequenced 540 bp of the mitochondrial 16S rRNA gene from 40 populations of mountain frogs from Anatolia, the Elburz Mountains and the Caucasus. Our samples comprise all known species and subspecies: Rana macrocnemis macrocnemis, R. m. tavasensis, R. m. pseudodalmatina, R. camerani and R. holtzi. They include the type localities of four of these taxa. We used a nested clade analysis (NCA) to infer historical and recurrent events that account for the observed geographical distribution of haplotypes. None of the extant species is monophyletic. Based on a molecular clock calibration using homologous sequences of Western Palearctic water frogs of the same genus, we estimated that a basic radiation into three lineages c. 2 Mya was followed by several dispersal and fragmentation events. The geographical distribution of resident and widespread haplotypes allows us to infer and date scenarios of range expansion and fragmentation that are aligned with dramatic climatic oscillations that have occurred during the last 600000 years. Consequently, Pliocene and Pleistocene climatic oscillations triggered the evolution of Anatolian mountain frogs through an interplay of vicariance and dispersal events.