Rangewide phylogeography and landscape genetics of the Western U.S. endemic frog Rana boylii (Ranidae): implications for the conservation of frogs and rivers

From eDNA to decisions using a multi-method approach to restoration planning in streams

Reintroduction efforts are increasingly used to mitigate biodiversity losses, but are frequently challenged by inadequate planning and uncertainty. High quality information about population status and threats can be used to prioritize reintroduction and restoration efforts and can transform ad hoc approaches into opportunities for improving conservation outcomes at a landscape scale. We conducted comprehensive environmental DNA (eDNA) and visual encounter surveys to determine the distribution of native and non-native aquatic species in two high-priority watersheds to address key uncertainties-such as the distribution of threats and the status of existing populations-inherent in restoration planning. We then used these occurrence data to develop a menu of potential conservation actions and a decision framework to benefit an endangered vertebrate (foothill yellow-legged frog, Rana boylii) in dynamic stream systems. Our framework combines the strengths of multiple methods, allowing managers and conservation scientists to incorporate conservation science and site-specific knowledge into the planning process to increase the likelihood of achieving conservation goals.

Biogeographic inferences across spatial and evolutionary scales

The field of biogeography unites landscape genetics and phylogeography under a common conceptual framework. Landscape genetics traditionally focuses on recent-time, population-based, spatial genetics processes at small geographical scales, while phylogeography typically investigates deep past, lineage- and species-based processes at large geographical scales. Here, we evaluate the link between landscape genetics and phylogeographical methods using the western fence lizard (Sceloporus occidentalis) as a model species. First, we conducted replicated landscape genetics studies across several geographical scales to investigate how population genetics inferences change depending on the spatial extent of the study area. Then, we carried out a phylogeographical study of population structure at two evolutionary scales informed by inferences derived from landscape genetics results to identify concordance and conflict between these sets of methods. We found significant concordance in landscape genetics processes at all but the largest geographical scale. Phylogeographical results indicate major clades are restricted to distinct river drainages or distinct hydrological regions. At a more recent timescale, we find minor clades are restricted to single river canyons in the majority of cases, while the remainder of river canyons include samples from at most two clades. Overall, the broad-scale pattern implicating stream and river valleys as key features linking populations in the landscape genetics results, and high degree of clade specificity within major topographic subdivisions in the phylogeographical results, is consistent. As landscape genetics and phylogeography share many of the same objectives, synthesizing theory, models and methods between these fields will help bring about a better understanding of ecological and evolutionary processes structuring genetic variation across space and time.

Stepping into the past to conserve the future: Archived skin swabs from extant and extirpated populations inform genetic management of an endangered amphibian

Moving animals on a landscape through translocations and reintroductions is an important management tool used in the recovery of endangered species, particularly for the maintenance of population genetic diversity and structure. Management of imperiled amphibian species rely heavily on translocations and reintroductions, especially for species that have been brought to the brink of extinction by habitat loss, introduced species, and disease. One striking example of amphibian declines and associated management efforts is in California's Sequoia and Kings Canyon National Parks with the mountain yellow-legged frog species complex (Rana sierrae/muscosa). Mountain yellow-legged frogs have been extirpated from more than 93% of their historic range, and limited knowledge of their population genetics has made long-term conservation planning difficult. To address this, we used 598 archived skin swabs from both extant and extirpated populations across 48 lake basins to generate a robust Illumina-based nuclear amplicon data set. We found that samples grouped into three main genetic clusters, concordant with watershed boundaries. We also found evidence for historical gene flow across watershed boundaries with a north-to-south axis of migration. Finally, our results indicate that genetic diversity is not significantly different between populations with different disease histories. Our study offers specific management recommendations for imperiled mountain yellow-legged frogs and, more broadly, provides a population genetic framework for leveraging minimally invasive samples for the conservation of threatened species.

Spatial and temporal patterns of environmental DNA detection to inform sampling protocols in lentic and lotic systems

The development of efficient sampling protocols for the capture of environmental DNA (eDNA) could greatly help improve accuracy of occupancy monitoring for species that are difficult to detect. However, the process of developing a protocol in situ is complicated for rare species by the fact that animal locations are often unknown. We tested sampling designs in lake and stream systems to determine the most effective eDNA sampling protocols for two rare species: the Sierra Nevada yellow-legged frog (Rana sierrae) and the foothill yellow-legged frog (Rana boylii). We varied water volume, spatial sampling, and seasonal timing in lakes and streams; in lakes we also tested multiple filter types. We found that filtering 2 L versus 1 L increased the odds of detection in streams 5.42X (95% CI: 3.2-9.19X) in our protocol, from a probability of 0.51-0.85 per technical replicate. Lake sample volumes were limited by filter clogging, and we found no effect of volume or filter type. Sampling later in the season increased the odds of detection in streams by 1.96X for every 30 days (95% CI: 1.3-2.97X) but there was no effect for lakes. Spatial autocorrelation of the quantity of yellow-legged frog eDNA captured in streams ceased between 100 and 200 m, indicating that sampling at close intervals is important.

Landscape genetics reveal broad and fine-scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east-west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east-west split occurred ~13-18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5-11 kya during a period of extreme drought cycles. A range-wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid-Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.

Population genomic data reveal extreme geographic subdivision and novel conservation actions for the declining foothill yellow-legged frog

Genomic data have the potential to inform high resolution landscape genetic and biological conservation studies that go far beyond recent mitochondrial and microsatellite analyses. We characterize the relationships of populations of the foothill yellow-legged frog, Rana boylii, a declining, "sentinel" species for stream ecosystems throughout its range in California and Oregon. We generated RADseq data and applied phylogenetic methods, hierarchical Bayesian clustering, PCA and population differentiation with admixture analyses to characterize spatial genetic structure across the species range. To facilitate direct comparison with previous analyses, we included many localities and individuals from our earlier work based on mitochondrial DNA. The results are striking, and emphasize the power of our landscape genomic approach. We recovered five extremely differentiated primary clades that indicate that R. boylii may be the most genetically differentiated anuran yet studied. Our results provide better resolution and more spatially consistent patterns than our earlier work, confirming the increased resolving power of genomic data compared to single-locus studies. Genomic structure is not equal across the species distribution. Approximately half the range of R. boylii consists of a single, relatively uniform population, while Sierra Nevada and coastal California clades are deeply, hierarchically substructured with biogeographic breaks observed in other codistributed taxa. Our results indicate that clades should serve as management units for R. boylii rather than previously suggested watershed boundaries, and that the near-extinct population from southwestern California is particularly diverged, exhibits the lowest genetic diversity, and is a critical conservation target for species recovery.

Rapid extirpation of a North American frog coincides with an increase in fungal pathogen prevalence: Historical analysis and implications for reintroduction

As extinctions continue across the globe, conservation biologists are turning to species reintroduction programs as one optimistic tool for addressing the biodiversity crisis. For repatriation to become a viable strategy, fundamental prerequisites include determining the causes of declines and assessing whether the causes persist in the environment. Invasive species-especially pathogens-are an increasingly significant factor contributing to biodiversity loss. We hypothesized that Batrachochytrium dendrobatidis (Bd), the causative agent of the deadly amphibian disease chytridiomycosis, was important in the rapid (<10 years) localized extirpation of a North American frog (Rana boylii) and that Bd remains widespread among extant amphibians in the region of extirpation. We used an interdisciplinary approach, combining interviews with herpetological experts, analysis of archived field notes and museum specimen collections, and field sampling of the extant amphibian assemblage to examine (1) historical relative abundance of R. boylii; (2) potential causes of R. boylii declines; and (3) historical and contemporary prevalence of Bd. We found that R. boylii were relatively abundant prior to their rapid extirpation, and an increase in Bd prevalence coincided with R. boylii declines during a time of rapid change in the region, wherein backcountry recreation, urban development, and the amphibian pet trade were all on the rise. In addition, extreme flooding during the winter of 1969 coincided with localized extirpations in R. boylii populations observed by interview respondents. We conclude that Bd likely played an important role in the rapid extirpation of R. boylii from southern California and that multiple natural and anthropogenic factors may have worked in concert to make this possible in a relatively short period of time. This study emphasizes the importance of recognizing historical ecological contexts in making future management and reintroduction decisions.

Mitochondrial genome of Pachytriton feii (Urodela: Salamandridae)

In this paper, the complete mitochondrial genome of Pachytriton feii is sequenced and reported for the first time. In our result, the complete mitochondrial genome of P. feii is 16 293 bp in length. Similar to the typical mtDNA of vertebrates, it contained 37 genes (13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes) and 1 non-coding region (D-loop). All the genes in P. feii were distributed on the H-strand, except for the ND6 subunit gene and eight tRNA genes which were encoded on the L-strand. The phylogenetic tree was reconstructed using the Bayesian analysis method and contained two major clades. Among Pachytriton (P. feii, Pachytriton labiatus, and Pachytriton brevipes), Paramesotriton laoensis make up the first lineage.

Extremely low genetic diversity of mtDNA control region and remarkable population differentiation of Ichthyophis bannanicus (Amphibia: Gymnophiona)

It is widely accepted that the mitochondrial DNA (mtDNA) control region (CR) gene evolves more quickly than protein-encoding genes, such as ND2 and Cyt b, with few exceptions. However, some species have a different evolution pattern. In the present study, we sequenced the mtDNA CR partial sequences (454 base pairs) of 142 individuals from five sampling sites of Ichthyophis bannanicus and compared the genetic diversity and structure with the information from the NADH dehydrogenase subunit 2 (ND2) and Cytochrome b (Cyt b) genes within this species. Extremely low genetic diversity was found in the mtDNA CR compared with those of the ND2 and Cyt b genes. These results showed that the relatively mean clock rate of the CR was broadly lower than those of the ND2 (about 2.55 times) and Cyt b (about 3.14 times) genes. Despite the extremely low genetic diversity of CR, the population structure analysis identified two groups, Xishuangbanna and Northern Vietnam-Yulin-Yangchun-Deqing, which indicated that the Red River systems may have acted as gene-flow barriers for I. bannanicus.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Mitochondrial genome of Emberiza elegans (Emberizidae: Emberiza)

Emberiza elegans is a common bunting with very wide geographical distribution. In this paper, the complete mitochondrial genome of E. elegans (16,779 bp in length) was analyzed for building the database. The results showed that it contained 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and one control region. The base composition of mtDNA was A (29.4%), G (14.7%), C (32.7%), and T (23.2%), therefore, the percentage of A and T (52.6%) was slightly higher than G and C (47.4%). All the genes in E. elegans were distributed on the H-strand, except for the ND6 subunit gene and ten tRNA genes which were encoded on the L-strand.

The complete mitochondrial genome of the Hynobius maoershanensis (Caudata: hynobiidae)

The Hynobius maoershanensis is a member of hynobiidae, endemic to Mountain Maoer in Guangxi province, China. It was first found and reported in 2006 and so far there is a little molecular research about it. The complete mitochondrial genome of H. maoershanensis has been obtained for the first time in this study. The circular genome (16,412 bp in length) consisted of 37 typical animal mitochondrial genes (13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes) and 1 control region. Overall base composition of the complete mitochondrial DNA was 33% A, 32% T, 21% C, and 14% G with AT (65%).

Mitochondrial genome of the Emberiza rutila (Emberizidae: Emberiza)

Emberiza rutila is a passerine bird of eastern Asia which belongs to the genus Emberiza in the bunting family Emberizidae. The complete mitochondrial genome of E. rutila was obtained for the first time in this study. The circular genome (16,803 bp in length) consists of 37 typical animal mitochondrial genes (13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes) and 1 control region. Overall base composition of the complete mitochondrial DNA was 30% A, 22.8% T, 32.8% C and 14.3% G. Except for 8 tRNA genes and ND6 gene, the relative position and orientation of all the genes were identical to those of most vertebrates.

Evolutionary dynamics of a rapidly receding southern range boundary in the threatened California Red-Legged Frog (Rana draytonii)

Populations forming the edge of a species range are often imperiled by isolation and low genetic diversity, with proximity to human population centers being a major determinant of edge stability in modern landscapes. Since the 1960s, the California red‐legged frog (Rana draytonii) has undergone extensive declines in heavily urbanized southern California, where the range edge has rapidly contracted northward while shifting its cardinal orientation to an east‐west trending axis. We studied the genetic structure and diversity of these frontline populations, tested for signatures of contemporary disturbance, specifically fire, and attempted to disentangle these signals from demographic events extending deeper into the past. Consistent with the genetic expectations of the ‘abundant‐center’ model, we found that diversity, admixture, and opportunity for random mating increases in populations sampled successively further away from the range boundary. Demographic simulations indicate that bottlenecks in peripheral isolates are associated with processes extending tens to a few hundred generations in the past, despite the demographic collapse of some due to recent fire‐flood events. While the effects of recent disturbance have left little genetic imprint on these populations, they likely contribute to an extinction debt that will lead to continued range contraction unless management intervenes to stall or reverse the process.

Mitochondrial genome of the Emberiza pusilla (Emberizidae: Emberiza)

The Emberiza pusilla is a common bunting with very wide geographical range. Here, the complete mitochondrial genome of E. pusilla (16,790 bp in length) has been analyzed for building the database. The results showed that it consisted of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region. The base composition of mtDNA was A (29.6%), G (14.8%), C (32.8%), and T (22.9%), so the percentage of A and T (52.5%) was slightly higher than G and C. All the genes in E. pusilla were distributed on the H-strand, except for the ND6 subunit gene and 10 tRNA genes, which were encoded on the L-strand.

Phylogeography and demographic history of Babina pleuraden (Anura, Ranidae) in southwestern China

Factors that determine genetic structure of species in southwestern China remain largely unknown. In this study, sequences of two mitochondrial genes (COI and cyt b) were determined to investigate the phylogeography and demography of Babina pleuraden, a pond frog endemic to southwestern China. A total of 262 individuals from 22 populations across the entire range of the species were collected. Our results indicate that B. pleuraden comprises five well-supported mitochondrial lineages roughly corresponding to five geographical areas. The phylogeographic structure of B. pleuraden has been shaped primarily by the unique regional responses of the Yunnan Plateau to the rapid uplift of the Qinghai-Tibetan Plateau occurred c. 2.5 Mya (B phrase of Qingzang Movement) and climatic oscillation during middle Pleistocene (c. 0.64–0.36 Mya), rather than by the paleo-drainage systems. The present wide distribution of the species has resulted from recent population expansion (c. 0.053–0.025 Mya) from multiple refugia prior to the Last Glacial Maximum, corresponding to the scenario of “refugia within refugia”.