Whole exome sequencing identifies the potential for genetic rescue in iconic and critically endangered Panamanian harlequin frogs

A roadmap for harlequin frog systematics, with a partial revision of Amazonian species related to Atelopus spumarius

Harlequin frogs, genus Atelopus, are a species-rich group of bufonid anurans from the Neotropics with more than 100 species. For nearly four decades now, this group has suffered from massive population declines. Almost all species are threatened with extinction, and many populations and several species are considered extinct or possibly extinct. This results in a limited sampling available for studies on harlequin frog systematics, especially in terms of molecular genetic information. However, efficient conservation of harlequin frogs requires an improved taxonomy. This is further complicated through the circumstance that many Atelopus species are relatively poor in external morphological characters combined with a high level of intra-specific character variation (e.g. coloration and body size). At the same time, cryptic diversity exists with well differentiated species (supported by osteology and molecular genetics) almost indistinguishable by external morphology. We compiled the largest dataset to date for mitochondrial (12S, 16S, cyt b) and nuclear (POMC, RAG1) markers and present a phylogeny (likelihood and Bayesian inference methods) including 152 samples from 104 populations scattered over the entire geographic range of the genus. Four allo- or parapatric main clades are distinguished: I. Sierra Nevada; II. Venezuelan-Andean; III. Andean-Chocó-Central American (with the ignescens and the varius-longirostris clades); and IV Amazonian (containing the tricolor and the flavescens-spumarius clades). The phylogenetic relationships within these clades remain to be resolved. Taxonomic implications included both splitting and lumping, but taxonomic action is here only taken for populations related to A. spumarius from western Amazonia. Besides redescriptions of A. spumarius sensu stricto and A. colomai, we describe two new species based on morphology, skull osteology and bioacoustics. Additional yet understudied populations from Amazonia may be allocated to these species or may represent additional undescribed taxa.

Four new genome sequences of the Pallas's cat (Otocolobus manul): an insight into the patterns of within-species variability

Manul (Otocolobus manul) is the only representative of the genus Otocolobus, which makes up the Leopard Cat lineage along with the genus Prionailurus. Their habitat is characterized by harsh environmental conditions. Although their populations are probably more stable than previously thought, it is still the case that their population size is declining. Conservation programs exist to protect manuls, but those based on captive breeding are often unsuccessful due to their increased susceptibility to diseases. The manul is therefore a suitable model species for evolutionary and diversity studies as well as for studying mechanisms of adaptation to harsh environment and mechanisms of susceptibility to diseases. Recently, the genome of the O. manul based on nanopore long-range sequencing has been published. Aiming to better understand inter- and intraspecific variation of the species, we obtained information on genome sequences of four other manuls, based on whole genome resequencing via the Illumina platform. On average, we detected a total of 3,636,571 polymorphic variants. Information on different types of structural variants and on the extent of SNP homozygosity, not available from the reference genome, was retrieved. The average whole-genome heterozygosity was almost identical to that found in the O. manul reference genome. In this context, we performed a more detailed analysis of the candidate gene EPAS1 potentially related to adaptation to the hypoxic environment. This analysis revealed both inter- and intraspecific variation, confirmed the presence of a previously described non-synonymous substitution in exon 15 unique to manuls and identified three additional unique non-synonymous substitutions located in so far not analyzed EPAS1 exonic sequences. The analysis of lncRNA located in the intron 7 of EPAS1 revealed interspecific variability and monomorphic nature of the sequence among analyzed manuls. The data obtained will allow more detailed analyses of the manul genome, focusing on genes and pathways involved in their adaptation to the environment and in susceptibility to diseases. This information can be helpful for optimizing conservation programs for this understudied species.

Outbreeding management offers the promise of genetic rescue for an endangered lizard

Inbreeding and loss of genetic diversity increase the extinction risk of small isolated populations. Genetic rescue by augmenting gene flow is a powerful means for the restoration of lost genetic variation. In this study, we used multigenerational pedigrees and neutral genetic markers to assess the consequences of outbreeding management in the Chinese and Vietnamese populations of the endangered crocodile lizard, Shinisaurus crocodilurus. Compared with the purebred Chinese population, the outbreeding population exhibited greater molecular genetic variation and a 3-fold larger population size. Moreover, the first-generation hybrids had a longer lifespan than purebreds, suggesting that outbreeding depression did not occur, but the long-term fitness effect of outbreeding needs to be further evaluated. Our study provides valuable insights into the potential for genetic rescue in the endangered crocodile lizard, emphasizing the importance of an evidence-based genetic management approach to address the risks of inbreeding and outbreeding depression in threatened populations.

What Can Frogs Teach Us about Resilience? Adaptive Renewal in Amphibian and Academic Ecosystems

Examples of resilience in nature give us hope amid a growing biodiversity crisis. While resilience has many definitions across disciplines, here I discuss resilience as the ability to continue to adapt and persist. Naturally, as biologists, we seek to uncover the underlying mechanisms that can help us explain the secrets of resilience across scales, from individuals to species to ecosystems and beyond. Perhaps we also ponder what the secrets to resilience are in our own lives, in our own research practices, and academic communities. In this paper, I highlight insights gained through studies of amphibian resilience following a global disease outbreak to uncover shared patterns and processes linked to resilience across amphibian communities. I also reflect on how classical resilience heuristics could be more broadly applied to these processes and to our own academic communities. Focusing on the amphibian systems that I have worked in-the Golden Frogs of Panama (Atelopus zeteki/varius) and the Mountain Yellow-Legged Frogs of California (Rana muscosa/sierrae)-I highlight shared and unique characteristics of resilience across scales and systems and discuss how these relate to adaptive renewal cycles. Reflecting on this work and previous resilience scholarship, I also offer my own thoughts about academia and consider what lessons we could take from mapping our own adaptive trajectories and addressing threats to our own community resilience.

Rapid Evolution of Resistance and Tolerance Leads to Variable Host Recoveries following Disease-Induced Declines

AbstractRecoveries of populations that have suffered severe disease-induced declines are being observed across disparate taxa. Yet we lack theoretical understanding of the drivers and dynamics of recovery in host populations and communities impacted by infectious disease. Motivated by disease-induced declines and nascent recoveries in amphibians, we developed a model to ask the following question: How does the rapid evolution of different host defense strategies affect the transient recovery trajectories of hosts following pathogen invasion and disease-induced declines? We found that while host life history is predictably a major driver of variability in population recovery trajectories (including declines and recoveries), populations that use different host defense strategies (i.e., tolerance, avoidance resistance, and intensity-reduction resistance) experience notably different recoveries. In single-species host populations, populations evolving tolerance recovered on average four times slower than populations evolving resistance. Moreover, while populations using avoidance resistance strategies had the fastest potential recovery rates, these populations could get trapped in long transient states at low abundance prior to recovery. In contrast, the recovery of populations evolving intensity-reduction resistance strategies were more consistent across ecological contexts. Overall, host defense strategies strongly affect the transient dynamics of population recovery and may affect the ultimate fate of real populations recovering from disease-induced declines.

Coevolution of a generalist pathogen with many hosts: the case of the amphibian chytrid Batrachochytrium dendrobatidis

Generalist pathogens maintain infectivity in numerous hosts; how this broad ecological niche impacts host-pathogen coevolution remains to be widely explored. Batrachochytrium dendrobatidis (Bd) is a highly generalist pathogenic fungus that has caused devastating declines in hundreds of amphibian species worldwide. This review examines amphibian chytridiomycosis host-pathogen interactions and available evidence for coevolution between Bd and its numerous hosts. We summarize recent evidence showing that Bd genotypes vary in geographic distribution and virulence, and that amphibian species also vary in Bd susceptibility according to their geographic distribution. How much variation can be explained by phenotypic plasticity or genetic differences remains uncertain. Recent research suggests that Bd genotypes display preferences for specific hosts and that some hosts are undergoing evolution as populations rebound from Bd outbreaks. Taken together, these findings suggest the potential for coevolution to occur and illuminate a path for addressing open questions through integrating historical and contemporary genetic data.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Importance of Genetic-Fitness Correlations for the Conservation of Amphibians

Endangered animals suffer from isolation of their habitats. Isolation leads to a reduction in population size as well as a decrease in genetic diversity and a concomitant increase in the risk of extinction. Amphibians are the most endangered vertebrate class. Besides habitat loss, fragmentation and isolation, amphibians are threatened by emerging diseases e.g., chytrid fungus or Ranavirus. By employing experiments, researchers investigate whether changes in genetic diversity within or among isolated populations affect amphibian fitness. While genetic diversity estimates are based on molecular markers, typically microsatellites, fitness is mostly measured as tadpole performance in rearing experiments often under varying environmental conditions. Tadpole performances (e.g., body mass, growth rate and survival) have been found to be negatively affected by low genetic diversity, as several studies have found a positive association between genetic diversity and these fitness traits. Moreover, infection with pathogens also seems to be more likely in individuals or populations with lower genetic diversity. Overall, these genetic-fitness correlations seem to be more pronounced or detectable in smaller, declining populations but not in larger populations. Genomic studies, which sample a larger fraction of the genome, are still scarce in the conservation genetic literature on amphibians. These are likely to increase in upcoming years and may reveal adaptive variants that protect against dangerous pathogens or environmental changes. Altogether, genetic-fitness correlation studies should be a priority in order to develop effective management plans for the genetic rescue of isolated, imperilled amphibian populations.

Genetic mixing in conservation translocations increases diversity of a keystone threatened species, Bettongia lesueur

Translocation programmes are increasingly being informed by genetic data to monitor and enhance conservation outcomes for both natural and established populations. These data provide a window into contemporary patterns of genetic diversity, structure and relatedness that can guide managers in how to best source animals for their translocation programmes. The inclusion of historical samples, where possible, strengthens monitoring by allowing assessment of changes in genetic diversity over time and by providing a benchmark for future improvements in diversity via management practices. Here, we used reduced representation sequencing (ddRADseq) data to report on the current genetic health of three remnant and seven translocated boodie (Bettongia lesueur) populations, now extinct on the Australian mainland. In addition, we used exon capture data from seven historical mainland specimens and a subset of contemporary samples to compare pre-decline and current diversity. Both data sets showed the significant impact of population founder source (whether multiple or single) on the genetic diversity of translocated populations. Populations founded by animals from multiple sources showed significantly higher genetic diversity than the natural remnant and single-source translocation populations, and we show that by mixing the most divergent populations, exon capture heterozygosity was restored to levels close to that observed in pre-decline mainland samples. Relatedness estimates were surprisingly low across all contemporary populations and there was limited evidence of inbreeding. Our results show that a strategy of genetic mixing has led to successful conservation outcomes for the species in terms of increasing genetic diversity and provides strong rationale for mixing as a management strategy.

When Defenses Fail: Atelopus zeteki Skin Secretions Increase Growth of the Pathogen Batrachochytrium dendrobatidis

To combat the threat of emerging infectious diseases in wildlife, ecoimmunologists seek to understand the complex interactions among pathogens, their hosts, and their shared environments. The cutaneous fungal pathogen Batrachochytrium dendrobatidis (Bd), has led to the decline of innumerable amphibian species, including the Panamanian golden frog (Atelopus zeteki). Given that Bd can evade or dampen the acquired immune responses of some amphibians, nonspecific immune defenses are thought to be especially important for amphibian defenses against Bd. In particular, skin secretions constitute a vital component of amphibian innate immunity against skin infections, but their role in protecting A. zeteki from Bd is unknown. We investigated the importance of this innate immune component by reducing the skin secretions from A. zeteki and evaluating their effectiveness against Bd in vitro and in vivo. Following exposure to Bd in a controlled inoculation experiment, we compared key disease characteristics (e.g., changes in body condition, prevalence, pathogen loads, and survival) among groups of frogs that had their skin secretions reduced and control frogs that maintained their skin secretions. Surprisingly, we found that the skin secretions collected from A. zeteki increased Bd growth in vitro. This finding was further supported by infection and survival patterns in the in vivo experiment where frogs with reduced skin secretions tended to have lower pathogen loads and survive longer compared to frogs that maintained their secretions. These results suggest that the skin secretions of A. zeteki are not only ineffective at inhibiting Bd but may enhance Bd growth, possibly leading to greater severity of disease and higher mortality in this highly vulnerable species. These results differ from those of previous studies in other amphibian host species that suggest that skin secretions are a key defense in protecting amphibians from developing severe chytridiomycosis. Therefore, we suggest that the importance of immune components cannot be generalized across all amphibian species or over time. Moreover, the finding that skin secretions may be enhancing Bd growth emphasizes the importance of investigating these immune components in detail, especially for species that are a conservation priority.

Linking pathogen-microbiome-host interactions to explain amphibian population dynamics

Symbiotic interactions can determine the evolutionary trajectories of host species, influencing genetic variation through selection and changes in demography. In the context of strong selective pressures such as those imposed by infectious diseases, symbionts providing defences could contribute to increase host fitness upon pathogen emergence. Here, we generated genome-wide data of an amphibian species to find evidence of evolutionary pressures driven by two skin symbionts: a batrachochytrid fungal pathogen and an antifungal bacterium. Using demographic modelling, we found evidence of decreased effective population size, probably due to pathogen infections. Additionally, we investigated host genetic associations with infection status, antifungal bacterium abundance and overall microbiome diversity using structural equation models. We uncovered relatively lower nucleotide diversity in infected frogs and potential heterozygote advantage to recruit the candidate beneficial symbiont and fight infections. Our models indicate that environmental conditions have indirect effects on symbiont abundance through both host body traits and microbiome diversity. Likewise, we uncovered a potential offsetting effect among host heterozygosity-fitness correlations, plausibly pointing to different ecological and evolutionary processes among the three species due to dynamic interactions. Our findings revealed that evolutionary pressures not only arise from the pathogen but also from the candidate beneficial symbiont, and both interactions shape the genetics of the host. Our results advance knowledge about multipartite symbiotic relationships and provide a framework to model ecological and evolutionary dynamics in wild populations. Finally, our study approach can be applied to inform conservation actions such as bioaugmentation strategies for other imperilled amphibians affected by infectious diseases.

A review of chemical defense in harlequin toads (Bufonidae: Atelopus)

Toads of the genus Atelopus are chemically defended by a unique combination of endogenously synthesized cardiotoxins (bufadienolides) and neurotoxins which may be sequestered (guanidinium alkaloids). Investigation into Atelopus small-molecule chemical defenses has been primarily concerned with identifying and characterizing various forms of these toxins while largely overlooking their ecological roles and evolutionary implications. In addition to describing the extent of knowledge about Atelopus toxin structures, pharmacology, and biological sources, we review the detection, identification, and quantification methods used in studies of Atelopus toxins to date and conclude that many known toxin profiles are unlikely to be comprehensive because of methodological and sampling limitations. Patterns in existing data suggest that both environmental (toxin availability) and genetic (capacity to synthesize or sequester toxins) factors influence toxin profiles. From an ecological and evolutionary perspective, we summarize the possible selective pressures acting on Atelopus toxicity and toxin profiles, including predation, intraspecies communication, disease, and reproductive status. Ultimately, we intend to provide a basis for future ecological, evolutionary, and biochemical research on Atelopus.