Functional variation at an expressed MHC class IIβ locus associates with Ranavirus infection intensity in larval anuran populations

Evolutionary analysis of major histocompatibility complex variants in chytrid-resistant and susceptible amphibians

An amphibian emerging infectious disease (EID), chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), originated in Asia but primarily led to declines and extinctions in amphibian populations outside of Asia. Host major histocompatibility complex (MHC) molecules exhibit high polymorphism, and the evolution of MHC can be influenced by recombination and pathogens. Previous studies have indicated that host MHC class II is associated with Bd resistance. In this study, I conducted recombination and selection tests on functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) and an Australasian Bd-susceptible species (Litoria caerulea). Recombination at the same site was identified in both species, supporting the hypothesis that recombination contributes to MHC IIß1 diversity in amphibians. Positive selection was observed in MHC IIß1 alleles in both species. In L. caerulea, at least four amino acid sites were identified under significant positive selection in the MHC IIß1, whereas these sites were either negatively selected or conserved in B. gargarizans. This suggests these sites might be selected for Bd resistance. Hydrophobicity was detected in certain amino acid sites relating to Bd resistance, suggesting this physicochemical property may be a factor selected to counteract Bd infection. These findings of this study provide an evolutionary basis for understanding how amphibian MHC IIß1 may undergo selection in response to chytrid infection.

Association between the skin microbiome and MHC class II diversity in an amphibian

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

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.

The wood frog (Rana sylvatica): An emerging comparative model for anuran immunity and host-ranavirus interactions

The wood frog (Rana sylvatica) is widely distributed across North America and is the only amphibian found north of the Arctic Circle due to its remarkable ability to tolerate whole-body freezing. Recent mass mortalities attributable to Ranavirus spp. (family Iridoviridae) in wild juvenile wood frogs, coupled with the apparent high susceptibility of wood frogs to experimental infection with frog virus 3 (FV3), the type species of the Ranavirus genus, or FV3-like isolates underscore the serious threat ranaviruses poses to wood frog populations. Despite the ecological relevance and unique life history of wood frogs, our understanding of the wood frog immune system and antiviral response to ranaviral infections is in its infancy. Here we aim to (1) synthesize the limited knowledge of wood frog immune defences, (2) review recent progress in establishing the wood frog as a study system for ranavirus infection, and (3) highlight the future use of wood frogs as a model anuran to provide insight into the evolution of anuran immune systems and antiviral responses.

The amphibian immune system

Amphibians are at the forefront of bridging the evolutionary gap between mammals and more ancient, jawed vertebrates. Currently, several diseases have targeted amphibians and understanding their immune system has importance beyond their use as a research model. The immune system of the African clawed frog, Xenopus laevis, and that of mammals is well conserved. We know that several features of the adaptive and innate immune system are very similar for both, including the existence of B cells, T cells and innate-like T cells. In particular, the study of the immune system at early stages of development is benefitted by studying X. laevis tadpoles. The tadpoles mainly rely on innate immune mechanisms including pre-set or innate-like T cells until after metamorphosis. In this review we lay out what is known about the innate and adaptive immune system of X. laevis including the lymphoid organs as well as how other amphibian immune systems are similar or different. Furthermore, we will describe how the amphibian immune system responds to some viral, bacterial and fungal insults. This article is part of the theme issue 'Amphibian immunity: stress, disease and ecoimmunology'.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Divergent allele advantage in the MHC and amphibian emerging infectious disease

Genetic variation in the major histocompatibility complex (MHC) may be associated with resistance to the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). The pathogen originated in Asia, then spread worldwide, causing amphibian population declines and species extinctions. We compared the expressed MHC IIβ1 alleles of a Bd-resistant toad species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian frog species, Litoria caerulea. We found at least six expressed MHC IIβ1 loci in each of the two species. Amino acid diversity encoded by these MHC alleles was similar between species, but the genetic divergence of those alleles known for broader pathogen-derived peptide binding was greater in the Bd-resistant species. In addition, we found a potentially rare allele in one resistant individual from the Bd-susceptible species. Deep next-generation sequencing recovered approximately triple the genetic resolution accessible from traditional cloning-based genotyping. Targeting more than one MHC IIβ1 expressed locus enables us to better understand how host MHC may adapt to emerging infectious diseases.

Amphibians as a model to study the role of immune cell heterogeneity in host and mycobacterial interactions

Mycobacterial infections represent major concerns for aquatic and terrestrial vertebrates including humans. Although our current knowledge is mostly restricted to Mycobacterium tuberculosis and mammalian host interactions, increasing evidence suggests common features in endo- and ectothermic animals infected with non-tuberculous mycobacteria (NTMs) like those described for M. tuberculosis. Importantly, most of the pathogenic and non-pathogenic NTMs detected in amphibians from wild, farmed, and research facilities represent, in addition to the potential economic loss, a rising concern for human health. Upon mycobacterial infection in mammals, the protective immune responses involving the innate and adaptive immune systems are highly complex and therefore not fully understood. This complexity results from the versatility and resilience of mycobacteria to hostile conditions as well as from the immune cell heterogeneity arising from the distinct developmental origins according with the concept of layered immunity. Similar to the differing responses of neonates versus adults during tuberculosis development, the pathogenesis and inflammatory responses are stage-specific in Xenopus laevis during infection by the NTM M. marinum. That is, both in human fetal and neonatal development and in tadpole development, responses are characterized by hypo-responsiveness and a lower capacity to contain mycobacterial infections. Similar to a mammalian fetus and neonates, T cells and myeloid cells in Xenopus tadpoles and axolotls are different from the adult immune cells. Fetal and amphibian larval T cells, which are characterized by a lower T cell receptor (TCR) repertoire diversity, are biased toward regulatory function, and they have distinct progenitor origins from those of the adult immune cells. Some early developing T cells and likely macrophage subpopulations are conserved in adult anurans and mammals, and therefore, they likely play an important role in the host-pathogen interactions from early stages of development to adulthood. Thus, we propose the use of developing amphibians, which have the advantage of being free-living early in their development, as an alternative and complementary model to study the role of immune cell heterogeneity in host-mycobacteria interactions.

Site- and Individual-Level Contaminations Affect Infection Prevalence of an Emerging Infectious Disease of Amphibians

Emerging infectious disease outbreaks are one of multiple stressors responsible for amphibian declines globally. In the northeastern United States, ranaviral diseases are prevalent in amphibians and other ectothermic species, but there is still uncertainty as to whether their presence is leading to population-level effects. Further, there is also uncertainty surrounding the potential interactions among disease infection prevalence in free-ranging animals and habitat degradation (co-occurrence of chemical stressors). The present study was designed to provide field-based estimates of the relationship between amphibian disease and chemical stressors. We visited 40 wetlands across three protected areas, estimated the prevalence of ranavirus among populations of larval wood frogs and spotted salamanders, and assessed chemical and biological stressors in wetland habitats and larval amphibians using a suite of selected bioassays, screening tools, and chemical analyses. Ranavirus was detected on larval amphibians from each protected area with an estimated occupancy ranging from 0.27 to 0.55. Considerable variation in ranavirus occupancy was also observed within and among each protected area. Of the stressors evaluated, ranavirus prevalence was strongly and positively related to concentrations of metalloestrogens (metals with the potential to bind to estrogen receptors) and total metals in wetland sediments and weakly and negatively related to total pesticide concentrations in larval amphibians. These results can be used by land managers to refine habitat assessments to include such environmental factors with the potential to influence disease susceptibility. Environ Toxicol Chem 2022;41:781-791. © 2022 SETAC. This article has been contributed to by US Government employees and their work is in the public domain in the USA.

Habitat fragmentation in the Brazilian Atlantic Forest is associated with erosion of frog immunogenetic diversity and increased fungal infections

Habitat fragmentation and infectious diseases threaten wildlife globally, but the interactions of these threats are poorly understood. For instance, while habitat fragmentation can impact genetic diversity at neutral loci, the impacts on disease-relevant loci are less well-studied. We examined the effects of habitat fragmentation in Brazil's Atlantic Forest on amphibian genetic diversity at an immune locus related to antigen presentation and detection (MHC IIB Exon 2). We used a custom high-throughput assay to sequence a fragment of MHC IIB and quantified Batrachochytrium dendrobatidis (Bd) infections in six frog species in two Atlantic Forest regions. Habitat fragmentation was associated with genetic erosion at MHC IIB Exon 2. This erosion was most severe in forest specialists. Significant Bd infections were detected only in one Atlantic Forest region, potentially due to relatively higher elevation. In this region, forest specialists showed an increase in both Bd prevalence and infection loads in fragmented habitats. Reduced population-level MHC IIB diversity was associated with increased Bd infection risk. On the individual level, MHC IIB heterozygotes exhibited a trend toward reduced Bd infection risk, although this was marginally non-significant. Our results suggest that habitat fragmentation increases Bd infection susceptibility in amphibians, mediated at least in part through erosion of immunogenetic diversity. Our findings have implications for management of fragmented populations in the face of emerging infectious diseases.

Balancing selection, genetic drift, and human-mediated introgression interplay to shape MHC (functional) diversity in Mediterranean brown trout

The extraordinary polymorphism of major histocompatibility complex (MHC) genes is considered a paradigm of pathogen-mediated balancing selection, although empirical evidence is still scarce. Furthermore, the relative contribution of balancing selection to shape MHC population structure and diversity, compared to that of neutral forces, as well as its interaction with other evolutionary processes such as hybridization, remains largely unclear. To investigate these issues, we analyzed adaptive (MHC-DAB gene) and neutral (11 microsatellite loci) variation in 156 brown trout (Salmo trutta complex) from six wild populations in central Italy exposed to introgression from domestic hatchery lineages (assessed with the LDH gene). MHC diversity and structuring correlated with those at microsatellites, indicating the substantial role of neutral forces. However, individuals carrying locally rare MHC alleles/supertypes were in better body condition (a proxy of individual fitness/parasite load) regardless of the zygosity status and degree of sequence dissimilarity of MHC, hence supporting balancing selection under rare allele advantage, but not heterozygote advantage or divergent allele advantage. The association between specific MHC supertypes and body condition confirmed in part this finding. Across populations, MHC allelic richness increased with increasing admixture between native and domestic lineages, indicating introgression as a source of MHC variation. Furthermore, introgression across populations appeared more pronounced for MHC than microsatellites, possibly because initially rare MHC variants are expected to introgress more readily under rare allele advantage. Providing evidence for the complex interplay among neutral evolutionary forces, balancing selection, and human-mediated introgression in shaping the pattern of MHC (functional) variation, our findings contribute to a deeper understanding of the evolution of MHC genes in wild populations exposed to anthropogenic disturbance.

Using PacBio SMRT data for identification of class I MHC alleles in a wildlife species, Zalophus californianus (California sea lion)

High allelic polymorphism and association with disease susceptibility has made the genes encoding major histocompatibility complex (MHC) antigen presentation molecules in humans, domesticated animals, and wildlife species of wide interest to ecologists, evolutionary biologists, and health specialists. The often multifaceted polygenism and extreme polymorphism of this immunogenetic system have made it especially difficult to characterize in non-model species. Here we compare and contrast the workflows of traditional Sanger sequencing of plasmid-cloned amplicons to Pacific Biosciences SMRT circular consensus sequencing (CCS) in their ability to capture alleles of MHC class I in a wildlife species where characterization of these genes was absent. We assessed two California sea lions (Zalophus californianus), a species suffering from a high prevalence of an aggressive cancer associated with a sexually transmitted gamma herpesvirus. In this pilot study, SMRT CCS proved superior in identifying more alleles from each animal than the more laborious plasmid cloning/Sanger workflow (12:7, 10:7), and no alleles were identified with the cloning/Sanger approach that were not identified by SMRT CCS. We discuss the advantages and disadvantages of each approach including cost, allele rarefaction, and sequence fidelity.

Depauperate major histocompatibility complex variation in the endangered reticulated flatwoods salamander (Ambystoma bishopi)

Reticulated flatwoods salamander (Ambystoma bishopi) populations began decreasing dramatically in the 1900s. Contemporary populations are small, isolated, and may be susceptible to inbreeding and reduced adaptive potential because of low genetic variation. Genetic variation at immune genes is especially important as it influences disease susceptibility and adaptation to emerging infectious pathogens, a central conservation concern for declining amphibians. We collected samples from across the extant range of this salamander to examine genetic variation at major histocompatibility complex (MHC) class Iα and IIβ exons as well as the mitochondrial control region. We screened tail or toe tissue for ranavirus, a pathogen associated with amphibian declines worldwide. Overall, we found low MHC variation when compared to other amphibian species and did not detect ranavirus at any site. MHC class Iα sequencing revealed only three alleles with a nucleotide diversity of 0.001, while MHC class IIβ had five alleles with a with nucleotide diversity of 0.004. However, unique variation still exists across this species' range with private alleles at three sites. Unlike MHC diversity, mitochondrial variation was comparable to levels estimated for other amphibians with nine haplotypes observed, including one haplotype shared across all sites. We hypothesize that a combination of a historic disease outbreak and a population bottleneck may have contributed to low MHC diversity while maintaining higher levels of mitochondrial DNA variation. Ultimately, MHC data indicated that the reticulated flatwoods salamander may be at an elevated risk from infectious diseases due to low levels of immunogenetic variation necessary to combat novel pathogens.