Anna E. Savage

MHC genotypes associate with resistance to a frog-killing fungus

The emerging amphibian disease chytridiomycosis is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian populations and species differ in susceptibility to Bd, yet we know surprisingly little about the genetic basis of this natural variation. MHC loci encode peptides that initiate acquired immunity in vertebrates, making them likely candidates for determining disease susceptibility. However, MHC genes have never been characterized in the context of chytridiomycosis. Here, we performed experimental Bd infections in laboratory-reared frogs collected from five populations that show natural variation in Bd susceptibility. We found that alleles of an expressed MHC class IIB locus associate with survival following Bd infection. Across populations, MHC heterozygosity was a significant predictor of survival. Within populations, MHC heterozygotes and individuals bearing MHC allele Q had a significantly reduced risk of death, and we detected a significant signal of positive selection along the evolutionary lineage leading to allele Q. Our findings demonstrate that immunogenetic variation affects chytridiomycosis survival under controlled experimental conditions, confirming that host genetic polymorphisms contribute to chytridiomycosis resistance.

Toward Immunogenetic Studies of Amphibian Chytridiomycosis: Linking Innate and Acquired Immunity

Recent declines in amphibian diversity and abundance have contributed significantly to the global loss of biodiversity. The fungal disease chytridiomycosis is widely considered to be a primary cause of these declines, yet the critical question of why amphibian species differ in susceptibility remains unanswered. Considerable evidence links environmental conditions and interspecific variability of the innate immune system to differential infection responses, but other sources of individual, population, or species-typical variation may also be important. In this article we review the preliminary evidence supporting a role for acquired immune defenses against chytridiomycosis, and advocate for targeted investigation of genes controlling acquired responses, as well as those that functionally bridge the innate and acquired immune systems. Immunogenetic data promise to answer key questions about chytridiomycosis susceptibility and host-pathogen coevolution, and will draw much needed attention to the importance of considering evolutionary processes in amphibian conservation management and practice.

Composition of symbiotic bacteria predicts survival in Panamanian golden frogs infected with a lethal fungus.

Symbiotic microbes can dramatically impact host health and fitness, and recent research in a diversity of systems suggests that different symbiont community structures may result in distinct outcomes for the host. In amphibians, some symbiotic skin bacteria produce metabolites that inhibit the growth of Batrachochytrium dendrobatidis (Bd), a cutaneous fungal pathogen that has caused many amphibian population declines and extinctions. Treatment with beneficial bacteria (probiotics) prevents Bd infection in some amphibian species and creates optimism for conservation of species that are highly susceptible to chytridiomycosis, the disease caused by Bd. In a laboratory experiment, we used Bd-inhibitory bacteria from Bd-tolerant Panamanian amphibians in a probiotic development trial with Panamanian golden frogs, Atelopus zeteki, a species currently surviving only in captive assurance colonies. Approximately 30% of infected golden frogs survived Bd exposure by either clearing infection or maintaining low Bd loads, but this was not associated with probiotic treatment. Survival was instead related to initial composition of the skin bacterial community and metabolites present on the skin. These results suggest a strong link between the structure of these symbiotic microbial communities and amphibian host health in the face of Bd exposure and also suggest a new approach for developing amphibian probiotics.

Fighting a Losing Battle: Vigorous Immune Response Countered by Pathogen Suppression of Host Defenses in the Chytridiomycosis-Susceptible Frog Atelopus zeteki

The emergence of the disease chytridiomycosis caused by the chytrid fungus Batrachochytrium dendrobatidis (Bd) has been implicated in dramatic global amphibian declines. Although many species have undergone catastrophic declines and/or extinctions, others appear to be unaffected or persist at reduced frequencies after Bd outbreaks. The reasons behind this variance in disease outcomes are poorly understood: differences in host immune responses have been proposed, yet previous studies suggest a lack of robust immune responses to Bd in susceptible species. Here, we sequenced transcriptomes from clutch-mates of a highly susceptible amphibian, Atelopus zeteki, with different infection histories. We found significant changes in expression of numerous genes involved in innate and inflammatory responses in infected frogs despite high susceptibility to chytridiomycosis. We show evidence of acquired immune responses generated against Bd, including increased expression of immunoglobulins and major histocompatibility complex genes. In addition, fungal-killing genes had significantly greater expression in frogs previously exposed to Bd compared with Bd-naïve frogs, including chitinase and serine-type proteases. However, our results appear to confirm recent in vitro evidence of immune suppression by Bd, demonstrated by decreased expression of lymphocyte genes in the spleen of infected compared with control frogs. We propose susceptibility to chytridiomycosis is not due to lack of Bd-specific immune responses but instead is caused by failure of those responses to be effective. Ineffective immune pathway activation and timing of antibody production are discussed as potential mechanisms. However, in light of our findings, suppression of key immune responses by Bd is likely an important factor in the lethality of this fungus.

Adaptive tolerance to a pathogenic fungus drives major histocompatibility complex evolution in natural amphibian populations.

Amphibians have been affected globally by the disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), and we are just now beginning to understand how immunogenetic variability contributes to disease susceptibility. Lineages of an expressed major histocompatibility complex (MHC) class II locus involved in acquired immunity are associated with chytridiomycosis susceptibility in controlled laboratory challenge assays. Here, we extend these findings to natural populations that vary both in exposure and response to Bd. We find that MHC alleles and supertypes associated with Bd survival in the field show a molecular signal of positive selection, while those associated with susceptibility do not, supporting the hypothesis that heritable Bd tolerance is rapidly evolving. We compare MHC supertypes to neutral loci to demonstrate where selection versus demography is shaping MHC variability. One population with Bd tolerance in nature shows a significant signal of directional selection for the same allele (allele Q) that was significantly associated with survival in an earlier laboratory study. Our findings indicate that selective pressure for Bd survival drives rapid immunogenetic adaptation in some natural populations, despite differences in environment and demography. Our field-based analysis of immunogenetic variation confirms that natural amphibian populations have the evolutionary potential to adapt to chytridiomycosis.

Seasonal and ontogenetic variation of skin microbial communities and relationships to natural disease dynamics in declining amphibians

Recently, microbiologists have focused on characterizing the probiotic role of skin bacteria for amphibians threatened by the fungal disease chytridiomycosis. However, the specific characteristics of microbial diversity required to maintain health or trigger disease are still not well understood in natural populations. We hypothesized that seasonal and developmental transitions affecting susceptibility to chytridiomycosis could also alter the stability of microbial assemblages. To test our hypothesis, we examined patterns of skin bacterial diversity in two species of declining amphibians (Lithobates yavapaiensis and Eleutherodactylus coqui) affected by the pathogenic fungus Batrachochytrium dendrobatidis (Bd). We focused on two important transitions that affect Bd susceptibility: ontogenetic (from juvenile to adult) shifts in E. coqui and seasonal (from summer to winter) shifts in L. yavapaiensis. We used a combination of community-fingerprinting analyses and 16S rRNA amplicon sequencing to quantify changes in bacterial diversity and assemblage composition between seasons and developmental stages, and to investigate the relationship between bacterial diversity and pathogen load. We found that winter-sampled frogs and juveniles, two states associated with increased Bd susceptibility, exhibited higher diversity compared with summer-sampled frogs and adult individuals. Our findings also revealed that hosts harbouring higher bacterial diversity carried lower Bd infections, providing support for the protective role of bacterial communities. Ongoing work to understand skin microbiome resilience after pathogen disturbance has the potential to identify key taxa involved in disease resistance.

Selection, trans-species polymorphism, and locus identification of major histocompatibility complex class IIβ alleles of New World ranid frogs.

Genes encoded by the major histocompatibility complex (MHC) play key roles in the vertebrate immune system. However, our understanding of the evolutionary processes and underlying genetic mechanisms shaping these genes is limited in many taxa, including amphibians, a group currently impacted by emerging infectious diseases. To further elucidate the evolution of the MHC in frogs (anurans) and develop tools for population genetics, we surveyed allelic diversity of the MHC class II β1 domain in both genomic and complementary DNA of seven New World species in the genus Rana (Lithobates). To assign locus affiliation to our alleles, we used a “gene walking” technique to obtain intron 2 sequences that flanked MHC class IIβ exon 2. Two distinct intron sequences were recovered, suggesting the presence of at least two class IIβ loci in Rana. We designed a primer pair that successfully amplified an orthologous locus from all seven Rana species. In total, we recovered 13 alleles and documented trans-species polymorphism for four of the alleles. We also found quantitative evidence of selection acting on amino acid residues that are putatively involved in peptide binding and structural stability of the β1 domain of anurans. Our results indicated that primer mismatch can result in polymerase chain reaction (PCR) bias, which influences the number of alleles that are recovered. Using a single locus may minimize PCR bias caused by primer mismatch, and the gene walking technique was an effective approach for generating single-copy orthologous markers necessary for future studies of MHC allelic variation in natural amphibian populations.

Positive Darwinian Selection at Single Amino Acid Sites Conferring Plant Virus Resistance

Explicit evaluation of the accuracy and power of maximum likelihood and Bayesian methods for detecting site-specific positive Darwinian selection presents a challenge because selective consequences of single amino acid changes are generally unknown. We exploited extensive molecular and functional characterization of amino acid substitutions in the plant gene eIF4E to evaluate the performance of these methods in detecting site-specific positive selection. We documented for the first time a molecular signature of positive selection within a recessive resistance gene in plants. We then used two statistical platforms, Phylogenetic Analysis Using Maximum Likelihood and Hypothesis Testing Using Phylogenies (HyPhy), to look for site-specific positive selection. Their relative power and accuracy are assessed by comparing the sites they identify as being positively selected with those of resistance-determining amino acids. Our results indicate that although both methods are surprisingly accurate in their identification of resistance sites, HyPhy appears to more accurately identify biologically significant amino acids using our data set.

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2–4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.

Linking genetic and environmental factors in amphibian disease risk.

A central question in evolutionary biology is how interactions between organisms and the environment shape genetic differentiation. The pathogen Batrachochytrium dendrobatidis (Bd) has caused variable population declines in the lowland leopard frog (Lithobates yavapaiensis); thus, disease has potentially shaped, or been shaped by, host genetic diversity. Environmental factors can also influence both amphibian immunity and Bd virulence, confounding our ability to assess the genetic effects on disease dynamics. Here, we used genetics, pathogen dynamics, and environmental data to characterize L. yavapaiensis populations, estimate migration, and determine relative contributions of genetic and environmental factors in predicting Bd dynamics. We found that the two uninfected populations belonged to a single genetic deme, whereas each infected population was genetically unique. We detected an outlier locus that deviated from neutral expectations and was significantly correlated with mortality within populations. Across populations, only environmental variables predicted infection intensity, whereas environment and genetics predicted infection prevalence, and genetic diversity alone predicted mortality. At one locality with geothermally elevated water temperatures, migration estimates revealed source–sink dynamics that have likely prevented local adaptation. We conclude that integrating genetic and environmental variation among populations provides a better understanding of Bd spatial epidemiology, generating more effective conservation management strategies for mitigating amphibian declines.