Characterisation of major histocompatibility complex class I in the Australian cane toad, Rhinella marina

The Amphibian Major Histocompatibility Complex-A Review and Future Outlook

The major histocompatibility complex (MHC) is a cluster of functionally related genes encoding proteins which, among other functions, mediate immune system activation. While the MHC of many vertebrates has been extensively studied, less is known about the amphibian MHC. This represents an important knowledge gap because amphibians mark the evolutionary transition from an aquatic to a terrestrial lifestyle and often maintain a biphasic lifestyle. Hence, they tend to be exposed to both aquatic and terrestrial pathogen communities, providing opportunities to gain fundamental insights into how the immune system responds to different environmental challenges. Moreover, amphibians are globally threatened by invasive pathogens and the MHC may play a role in combating population decline. In this review, we summarize the current state of knowledge regarding the amphibian MHC and identify the major differences with other vertebrates. We also review how the number of MHC gene copies varies across amphibian groups and how MHC-based variation relates to amphibian ontogeny, behaviour, disease, and phylogeography. We conclude by identifying knowledge gaps and proposing priorities for future research.

Whole-mitogenome analysis unveils previously undescribed genetic diversity in cane toads across their invasion trajectory

Invasive species offer insights into rapid adaptation to novel environments. The iconic cane toad (Rhinella marina) is an excellent model for studying rapid adaptation during invasion. Previous research using the mitochondrial NADH dehydrogenase 3 (ND3) gene in Hawai'ian and Australian invasive populations found a single haplotype, indicating an extreme genetic bottleneck following introduction. Nuclear genetic diversity also exhibited reductions across the genome in these two populations. Here, we investigated the mitochondrial genomics of cane toads across this invasion trajectory. We created the first reference mitochondrial genome for this species using long-read sequence data. We combined whole-genome resequencing data of 15 toads with published transcriptomic data of 125 individuals to construct nearly complete mitochondrial genomes from the native (French Guiana) and introduced (Hawai'i and Australia) ranges for population genomic analyses. In agreement with previous investigations of these populations, we identified genetic bottlenecks in both Hawai'ian and Australian introduced populations, alongside evidence of population expansion in the invasive ranges. Although mitochondrial genetic diversity in introduced populations was reduced, our results revealed that it had been underestimated: we identified 45 mitochondrial haplotypes in Hawai'ian and Australian samples, none of which were found in the native range. Additionally, we identified two distinct groups of haplotypes from the native range, separated by a minimum of 110 base pairs (0.6%). These findings enhance our understanding of how invasion has shaped the genetic landscape of this species.

Climate match is key to predict range expansion of the world's worst invasive terrestrial vertebrates

Understanding the determinants of the range expansion of invasive alien species is crucial for developing effective prevention and control strategies. Nevertheless, we still lack a global picture of the potential factors influencing the invaded range expansion across taxonomic groups, especially for the world's worst invaders with high ecological and economic impacts. Here, by extensively collecting data on 363 distributional ranges of 19 of world's worst invasive terrestrial vertebrates across 135 invaded administrative jurisdictions, we observed remarkable variations in the range expansion across species and taxonomic groups. After controlling for taxonomic and geographic pseudoreplicates, model averaging analyses based on generalized additive mixed-effect models showed that species in invaded regions having climates more similar to those of their native ranges tended to undergo a larger range expansion. In addition, as proxies of propagule pressure and human-assisted transportation, the number of introduction events and the road network density were also important predictors facilitating the range expansion. Further variance partitioning analyses validated the predominant role of climate match in explaining the range expansion. Our study demonstrated that regions with similar climates to their native ranges could still be prioritized to prevent the spread of invasive species under the sustained global change.

Transcriptome annotation reveals minimal immunogenetic diversity among Wyoming toads, Anaxyrus baxteri

Briefly considered extinct in the wild, the future of the Wyoming toad (Anaxyrus baxteri) continues to rely on captive breeding to supplement the wild population. Given its small natural geographic range and history of rapid population decline at least partly due to fungal disease, investigation of the diversity of key receptor families involved in the host immune response represents an important conservation need. Population decline may have reduced immunogenetic diversity sufficiently to increase the vulnerability of the species to infectious diseases. Here we use comparative transcriptomics to examine the diversity of toll-like receptors and major histocompatibility complex (MHC) sequences across three individual Wyoming toads. We find reduced diversity at MHC genes compared to bufonid species with a similar history of bottleneck events. Our data provide a foundation for future studies that seek to evaluate the genetic diversity of Wyoming toads, identify biomarkers for infectious disease outcomes, and guide breeding strategies to increase genomic variability and wild release successes.

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.

A comparison of nonlethal sampling methods for amphibian gut microbiome analyses

Noninvasive sampling methods for studying intestinal microbiomes are widely applied in studies of endangered species and in those conducting temporal monitoring during manipulative experiments. Although existing studies show that noninvasive sampling methods among different taxa vary in their accuracy, no studies have yet been published comparing nonlethal sampling methods in adult amphibians. In this study, we compare microbiomes from two noninvasive sample types (faeces and cloacal swabs) to that of the large intestine in adult cane toads, Rhinella marina. We use 16S rRNA gene sequencing to investigate how microbial communities change along the digestive tract and which nonlethal sampling method better represents large intestinal microbiota. We found that cane toads' intestinal microbiota was dominated by Bacteroidetes, Proteobacteria and Firmicutes and, interestingly, we also saw a high proportion of Fusobacteria, which has previously been associated with marine species and changes in frog immunity. The large and small intestine of cane toads had a similar microbial composition, but the large intestine showed higher diversity. Our results indicate that cloacal swabs were more similar to large intestine samples than were faecal samples, and small intestine samples were significantly different from both nonlethal sample types. Our study provides valuable information for future investigations of the cane toad gut microbiome and validates the use of cloacal swabs as a nonlethal method to study changes in the large intestine microbiome. These data provide insights for future studies requiring nonlethal sampling of amphibian gut microbiota.

A Glimpse of the Peptide Profile Presentation by Xenopus laevis MHC Class I: Crystal Structure of pXela-UAA Reveals a Distinct Peptide-Binding Groove

The African clawed frog, Xenopus laevis, is a model species for amphibians. Before metamorphosis, tadpoles do not efficiently express the single classical MHC class I (MHC-I) molecule Xela-UAA, but after metamorphosis, adults express this molecule in abundance. To elucidate the Ag-presenting mechanism of Xela-UAA, in this study, the Xela-UAA structure complex (pXela-UAAg) bound with a peptide from a synthetic random peptide library was determined. The amino acid homology between the Xela-UAA and MHC-I sequences of different species is <45%, and these differences are fully reflected in the three-dimensional structure of pXela-UAAg. Because of polymorphisms and interspecific differences in amino acid sequences, pXela-UAAg forms a distinct peptide-binding groove and presents a unique peptide profile. The most important feature of pXela-UAAg is the two-amino acid insertion in the α2-helical region, which forms a protrusion of ∼3.8 Å that is involved in TCR docking. Comparison of peptide-MHC-I complex (pMHC-I) structures showed that only four amino acids in β2-microglobulin that were bound to MHC-I are conserved in almost all jawed vertebrates, and the most unique feature in nonmammalian pMHC-I molecules is that the AB loop bound β2-microglobulin. Additionally, the binding distance between pMHC-I and CD8 molecules in nonmammals is different from that in mammals. These unique features of pXela-UAAg provide enhanced knowledge of T cell immunity and bridge the knowledge gap regarding the coevolutionary progression of the MHC-I complex from aquatic to terrestrial species.

Increased Adaptive Variation Despite Reduced Overall Genetic Diversity in a Rapidly Adapting Invader

Invasive species often evolve rapidly following introduction despite genetic bottlenecks that may result from small numbers of founders; however, some invasions may not fit this “genetic paradox”. The invasive cane toad (Rhinella marina) displays high phenotypic variation across its introduced Australian range. Here, we used three genome-wide datasets to characterize their population structure and genetic diversity. We found that toads form three genetic clusters: 1) native range toads, 2) toads from the source population in Hawaii and long-established areas near introduction sites in Australia, and 3) toads from more recently established northern Australian sites. Although we find an overall reduction in genetic diversity following introduction, we do not see this reduction in loci putatively under selection, suggesting that genetic diversity may have been maintained at ecologically relevant traits, or that mutation rates were high enough to maintain adaptive potential. Nonetheless, toads encounter novel environmental challenges in Australia, and the transition between genetic clusters occurs at a point along the invasion transect where temperature rises and rainfall decreases. We identify environmentally associated loci known to be involved in resistance to heat and dehydration. This study highlights that natural selection occurs rapidly and plays a vital role in shaping the structure of invasive populations.

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n = 20) and A. lemur (n = 20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.

ORIGINS AND GENETIC DIVERSITY OF INTRODUCED POPULATIONS OF THE PUERTO RICAN RED-EYED COQUÍ, ELEUTHERODACTYLUS ANTILLENSIS, IN SAINT CROIX (U.S. VIRGIN ISLANDS) AND PANAMÁ

The Red-eyed Coquí, Eleutherodactylus antillensis, is a terrestrial frog endemic to the Puerto Rican Bank (Puerto Rico and numerous islands and cays off its eastern coast), in the eastern Caribbean Sea. The species was likely introduced in Saint Croix, an island c. 100 km southeast of Puerto Rico, in the late 1930s, and in Panamá City, Panamá, in the late 1950s or early 1960s, but the source(s) of these introductions are unknown. We analyzed sequence data from one mtDNA locus and four nuDNA introns to infer the origin(s) of the Saint Croix and Panamá City populations and quantify their genetic diversity. Saint Croix and Panamanian populations do not share any haplotypes, and they cluster with different native populations, suggesting that they are derived from separate sources in the Puerto Rican Bank. Patterns of population structure trace the probable sources of E. antillensis in Saint Croix to islands off Puerto Rico's eastern coast, which include Vieques, Culebra, Saint Thomas, Saint John, Tortola, and Virgin Gorda, and possibly to eastern Puerto Rico as well. In contrast, Panamá City E. antillensis probably originated from either western or eastern Puerto Rico. Genetic diversity in the introduced populations is similar to or lower than in populations in the species' native range, indicating that genetic diversity has not increased in the alien frogs. Our findings may facilitate the development of preventive measures to minimize introductions of non-native amphibians in the Caribbean and Central America.

Multiple major histocompatibility complex class I genes in Asian anurans: Ontogeny and phylogeny

Amphibians, as the first terrestrial vertebrates, offer a window into early major histocompatibility complex (MHC) evolution. We characterized the MHC class I of two Korean amphibians, the Asiatic toad (Bufo gargarizans) and the Japanese tree frog (Hyla japonica). We found at least four transcribed MHC class I (MHC I) loci, the highest number confirmed in any anuran to date. Furthermore, we identified MHC I transcripts in terrestrial adults, and possibly in aquatic larvae, of both species. We conducted a phylogenetic analysis based on MHC I sequence data and found that B. gargarizans and H. japonica cluster together in the superfamily Nobleobatrachia. We further identified three supertypes shared by the two species. Our results reveal substantial variation in the number of MHC I loci in anurans and suggest that certain supertypes have particular physiochemical properties that may confer pathogen resistance.

Variation in Major Histocompatibility Complex diversity in invasive cane toad populations

Context The cane toad (Rhinella marina), a native species of central and southern America, was introduced to Australia in 1935 as a biocontrol agent after a complex history of prior introductions. The population rapidly expanded and has since spread through much of the Australian landmass, with severe impacts on the endemic wildlife, primarily via toxicity to predators. The invasion process has taken its toll on the cane toad, with changes in the immunological capacity across the Australian invasive population. Aims To investigate the immunogenetic underpinnings of these changes, we studied the diversity of the Major Histocompatiblity Complex (MHC) genes in introduced cane toad populations. Methods We studied the diversity of two MHC genes (the classical class I UA locus and a class II DAB locus) and compared these with neutral microsatellite markers in toads from the Australian site of introduction and the Australian invasion front. We also included toads from Hawai’i, the original source of the Australian toads, to infer founder effect. Key results Diversity across all markers was low across Australian and Hawai’ian samples, consistent with a reduction in genetic diversity through multiple founder effects during the course of the successive translocations. In Australia, allelic diversity at the microsatellite markers and the UA locus was reduced at the invasion front, whereas all three alleles at the DAB locus were maintained in the invasion-front toads. Conclusions Loss of allelic diversity observed at the microsatellite markers and the UA locus could be the result of drift and bottlenecking along the invasion process, however, the persistence of DAB diversity warrants further investigation to disentangle the evolutionary forces influencing this locus. Implications Through the use of different molecular markers, we provide a preliminary description of the adaptive genetic processes occurring in this invasive population. The extremely limited MHC diversity may represent low immunogenetic competence across the Australian population, which could be exploited for invasive species management.

The things they carried: The pathogenic effects of old and new parasites following the intercontinental invasion of the Australian cane toad (Rhinella marina)

Brought to Australia in 1935 to control agricultural pests (from French Guiana, via Martinique, Barbados, Jamaica, Puerto Rico and Hawai'i), repeated stepwise translocations of small numbers of founders enabled the cane toad (Rhinella marina) to escape many parasites and pathogens from its native range. However, the infective organisms that survived the journey continue to affect the dynamics of the toad in its new environment. In Australia, the native-range lungworm Rhabdias pseudosphaerocephala decreases its host's cardiac capacity, as well as growth and survival, but not rate of dispersal. The lungworm is most prevalent in long-colonised areas within the toads' Australian range, and absent from the invasion front. Several parasites and pathogens of Australian taxa have host-shifted to cane toads in Australia; for example, invasion-front toads are susceptible to spinal arthritis caused by the soil bacterium, Ochrobactrum anthropi. The pentastome Raillietiella frenata has host-shifted to toads and may thereby expand its Australian range due to the continued range expansion of the invasive toads. Spill-over and spill-back of parasites may be detrimental to other host species; however, toads may also reduce parasite loads in native taxa by acting as terminal hosts. We review the impact of the toad's parasites and pathogens on the invasive anuran's biology in Australia, as well as collateral effects of toad-borne parasites and pathogens on other host species in Australia. Both novel and co-evolved pathogens and parasites may have played significant roles in shaping the rapid evolution of immune system responses in cane toads within their invaded range.

•

Invasive cane toads have lost many parasites due to serial translocations.

•

One native lungworm (Rhabdias pseudosphaerocephala) has been retained.

•

Toads have also acquired novel parasites and pathogens from Australian hosts.

•

Toads either amplify parasite numbers or act as a parasite sink.

•

Differences in immune function exist between toad populations within Australia.

Characterisation of major histocompatibility complex class I genes in Japanese Ranidae frogs

The major histocompatibility complex (MHC) is a key component of adaptive immunity in all jawed vertebrates, and understanding the evolutionary mechanisms that have shaped these genes in amphibians, one of the earliest terrestrial tetrapods, is important. We characterised MHC class I variation in three common Japanese Rana species (Rana japonica, Rana ornativentris and Rana tagoi tagoi) and identified a total of 60 variants from 21 individuals. We also found evolutionary signatures of gene duplication, recombination and balancing selection (including trans-species polymorphism), all of which drive increased MHC diversity. A unique feature of MHC class I from these three Ranidae species includes low synonymous differences per site (d_S) within species, which we attribute to a more recent diversification of these sequences or recent gene duplication. The resulting higher d_N/d_S ratio relative to other anurans studied could be related to stronger selection pressure at peptide binding sites. This is one of the first studies to investigate MHC in Japanese amphibians and permits further exploration of the polygenetic factors associated with resistance to infectious diseases.

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.

A genetic perspective on rapid evolution in cane toads (Rhinella marina)

The process of biological invasion exposes a species to novel pressures, in terms of both the environments it encounters and the evolutionary consequences of range expansion. Several invaders have been shown to exhibit rapid evolutionary changes in response to those pressures, thus providing robust opportunities to clarify the processes at work during rapid phenotypic transitions. The accelerating pace of invasion of cane toads (Rhinella marina) in tropical Australia during its 80-year history has been well characterized at the phenotypic level, including common-garden experiments that demonstrate heritability of several dispersal-relevant traits. Individuals from the invasion front (and their progeny) show distinctive changes in morphology, physiology and behaviour that, in combination, result in far more rapid dispersal than is true of conspecifics from long-colonized areas. The extensive body of work on cane toad ecology enables us to place into context studies of the genetic basis of these traits. Our analyses of differential gene expression from toads from both ends of this invasion-history transect reveal substantial upregulation of many genes, notably those involved in metabolism and cellular repair. Clearly, then, the dramatically rapid phenotypic evolution of cane toads in Australia has been accompanied by substantial shifts in gene expression, suggesting that this system is well suited to investigating the genetic underpinnings of invasiveness.