Identification of major histocompatibility complex genotypes associated with resistance to an amphibian emerging infectious disease.

Amphibian chytridiomycosis, caused by Batrachochytrium dendrobatidis (Bd), emerged from Asia and spread globally. By comparing functional MHC IIß1 alleles from an Asian Bd-resistant anuran species (Bufo gargarizans) with those of an Australasian Bd-susceptible species (Litoria caerulea), we identified MHC genotypes associated with Bd resistance. These alleles encode a glycine deletion (G90ß1) and adjacent motifs in the deepest pathogen-derived peptide-binding groove. Every Bd-resistant individual, but no susceptible individuals, possessed at least one allele encoding the variant. We detected trans-species polymorphism at the end of the MHC IIß1 sequences. The G90ß1 deletion was encoded by different alleles in the two species, suggesting it may have evolved independently in each species rather than having been derived from a common ancestor. These results are consistent with a scenario by which MHC adaptations that confer resistance to the pathogen have evolved by convergent evolution. Immunogenetic studies such as this are critical to ongoing conservation efforts.

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 species, Bufo gargarizans, from South Korea with those of a Bd-susceptible Australasian 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 distance of those alleles known for potential 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 the full MHC IIß1 enables us to better understand how host MHC may adapt to emerging infectious diseases.

Effects of land-use types and the exotic species, Hypochaeris radicata, on plant diversity in human-transformed landscapes of the biosphere reserve, Jeju Island, Korea.

Land-use and plant invasion influence biodiversity. Understanding the effects of land-use types and invasive plants on the ecosystem is crucial for better management and the development of strategic plans for increasing biodiversity in Jeju Island, Korea, a designated Biosphere Reserve by the United Nations Education, Scientific, and Cultural Organization. The effect of the most dominant invasive exotic species, Hypochaeris radicata, on the four land-use types of Jeju Island was investigated. Plant composition, soil characteristics, and plant diversity among four land-use types (cropland, green space, neglected land, and residential) were compared. Among the land-use types, croplands had the most diverse plant composition and the highest richness in exotic and native plant species. Croplands, such as tangerine orchards, which are widely distributed throughout Jeju Island, showed the highest plant diversity because of medium intensity disturbance caused by weed removal. The relative cover of H. radicata did not differ between land-use types. However, H. radicata invasion was negatively related with plant species richness, making this invasive species a threat to the biodiversity of native herbs present in land-use areas. H. radicata adapts to areas with a broad range of soil properties and a variety of land-use types. Therefore, it is crucial to monitor land-use types and patterns of plant invasion to guide the implementation of consistent management and conservation strategies for maintaining ecosystem integrity of the transformed habitat in Jeju Island.

Novel chytrid pathogen variants and the global amphibian pet trade.

Global wildlife trade spreads emerging infectious diseases that threaten biodiversity. The amphibian chytrid pathogen Batrachochytrium dendrobatidis (Bd) has caused population declines and species extinctions worldwide except in Asia. Fire-bellied toads (Bombina orientalis), exported in large numbers from Asia, are tolerant of Bd and carry hypervirulent ancestral chytrid BdAsia-1 variants. We assayed the virulence of a new isolate of BdAsia-1 on the model Australasian frog host Litoria caerulea. Infected individuals (n = 15) all showed rapid disease progression culminating in death, whereas sham-inoculated individuals (n = 10) presented no clinical signs of disease and all survived (log rank test, χ2 = 15.6, df = 1, p < 0.0001). The virulence of the new isolate of BdAsia-1 is comparable to the one we assayed previously (χ2 = 0.0, df = 1, p = 0.91). Internationally traded wildlife, even when they appear healthy, can carry hypervirulent variants of pathogens. Once new pathogen variants escape into the environment, native species that have had no opportunity to evolve resistance to them may perish. Our study suggests that hypervirulent pathogens are being spread by the international pet trade. Notifiable wildlife diseases attributable to locally endemic pathogens often fail to generate conservation concern so are rarely subject to border surveillance or import controls. Because of the danger novel variants pose, national border control agencies need to implement disease screening and quarantine protocols to ensure the safety of their endemic fauna.

Variantes Patógenas Nuevas de Quitridios y el Mercado Mundial de Anfibios Mascota Resumen El mercado mundial de fauna dispersa enfermedades infecciosas emergentes que amenazan a la biodiversidad. El quitridio patógeno de anfibios Batrachochytrium dendrobatidis (Bd) ha causado declinaciones poblacionales y la extinción de especies en todo el mundo excepto Asia. El sapo Bombina orientalis, exportado en grandes cantidades desde Asia, es tolerante al Bd y carga genéticamente las variantes ancestrales hipervirulentas de quitridio BdAsia-1. Analizamos la virulencia de una nueva cepa de BdAsia-1 con el modelo de la rana australo-asiática hospedera Litoria caerulea. Todos los individuos infectados (n = 15) mostraron una progresión acelerada de la enfermedad que culminaba con la muerte, mientras que los individuos con inoculación simulada (n = 10) no presentaron señales clínicas de la enfermedad y todos sobrevivieron (prueba log de rango, χ2 = 15.6, df = 1, p < 0.0001). La virulencia de la nueva cepa de BdAsia-1 es comparable a la que analizamos previamente (χ2 = 0.0, df = 1, p = 0.91). La fauna comercializada internacionalmente, incluso cuando parece estar saludable, puede portar variantes hipervirulentas de los patógenos. Una vez que un patógeno nuevo se introduce al ambiente, pueden perecer las especies nativas que no han tenido la oportunidad de evolucionar la resistencia a estos patógenos. Nuestro estudio sugiere que los patógenos hipervirulentos se están dispersando mediante el mercado internacional de mascotas. Con frecuencia las enfermedades silvestres notificables que pueden atribuirse a los patógenos endémicos no generan interés para la conservación, así que rara vez están sujetas a la vigilancia fronteriza o el control de importación. Debido al riesgo que representan las variantes nuevas, las agencias nacionales de control fronterizo necesitan implementar evaluaciones patológicas y protocolos de cuarentena para asegurar la seguridad de su fauna endémica.

Early-diverging fungal phyla: taxonomy, species concept, ecology, distribution, anthropogenic impact, and novel phylogenetic proposals.

The increasing number of new fungal species described from all over the world along with the use of genetics to define taxa, has dramatically changed the classification system of early-diverging fungi over the past several decades. The number of phyla established for non-Dikarya fungi has increased from 2 to 17. However, to date, both the classification and phylogeny of the basal fungi are still unresolved. In this article, we review the recent taxonomy of the basal fungi and re-evaluate the relationships among early-diverging lineages of fungal phyla. We also provide information on the ecology and distribution in Mucoromycota and highlight the impact of chytrids on amphibian populations. Species concepts in Chytridiomycota, Aphelidiomycota, Rozellomycota, Neocallimastigomycota are discussed in this paper. To preserve the current application of the genus Nephridiophaga (Chytridiomycota: Nephridiophagales), a new type species, Nephridiophaga blattellae, is proposed.

Coevolution between MHC Class I and Antigen-Processing Genes in Salamanders.

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

Phylogeographic study of the Bufo gargarizans species complex, with emphasis on Northeast Asia.

We conduct a phylogeographic and population genetic study of the Asiatic toad (Bufo gargarizans) to understand its evolutionary history, and the influence of geology and climate. A total of 292 individuals from 94 locations were genotyped for two mitochondrial loci (cytb, ND2) and five nuclear introns (Sox9-2, Rho-3, CCNB2-3, UCH-2, and DBI-2), and we performed a suite of phylogenetic, population genetic, and divergence dating analyses. The phylogenetic trees constructed using mitochondrial loci inferred B. gargarizans being divided into two major groups: China mainland and Northeast Asia (Northeast China, Russia, and Korean Peninsula). As with previous studies of this species, we recover population genetic structure not tied to geographic region. Additionally, we discover a new genetic clade restricted to Northeast Asia that points towards the Korean Peninsula being a glacial refugium during the Pleistocene. The weak phylogeographic pattern of B. gargarizans is likely the result of multiple biological, anthropogenic, and historical factors - robust dispersal abilities as a consequence of physiological adaptations, human translocation, geologic activity, and glacial cycles of the Pleistocene. We highlight the complex geologic and climatic history of Northeast Asia and encourage further research to understand its impact on the biodiversity in the region.

Structural implications of traditional agricultural landscapes on the functional diversity of birds near the Korean Demilitarized Zone.

Bird assemblages are sensitive to changes in landscape composition and the environment, such as those that result from drought. In this study, the relationship between landscape composition and avian functional diversity in traditional agricultural ecosystems in the Civilian Control Zone (CCZ) of Korea was examined. In addition, the resilience of biodiversity to changes in landscape elements resulting from drought conditions was investigated. The traditional agricultural landscape (TAL) of the sites studied was divided into three types: TAL 1 had a high proportion of rice paddies, TAL 2 included large forest areas, and TAL 3 represented areas with drylands. Of these, TAL 1 showed the highest species richness and functional richness, but these measures were most vulnerable to drought. Meanwhile, TAL 2 showed that the bird communities were more tolerant under drought event. This study shows that to conserve and enhance the diversity of birds in traditional agricultural landscapes of Northeast Asia, active management of forest areas is needed to protect bird populations. In addition, commercial pressures to develop this area will require urgent biodiversity conservation plans to protect the unique biodiversity of the Korean CCZ. This study thus provides landscape management guidance for conservation planning.

Ancestral chytrid pathogen remains hypervirulent following its long coevolution with amphibian hosts.

Many amphibian species around the world, except in Asia, suffer morbidity and mortality when infected by the emerging infectious pathogen Batrachochytrium dendrobatidis (Bd). A lineage of the amphibian chytrid fungus isolated from South Korean amphibians (BdAsia-1) is evolutionarily basal to recombinant global pandemic lineages (BdGPL) associated with worldwide amphibian population declines. In Asia, the Bd pathogen and its amphibian hosts have coevolved over 100 years or more. Thus, resilience of Asian amphibian populations to infection might result from attenuated virulence of endemic Bd lineages, evolved immunity to the pathogen or both. We compared susceptibilities of an Australasian amphibian, Litoria caerulea, known to lack resistance to BdGPL, with those of three Korean species, Bufo gargarizans, Bombina orientalis and Hyla japonica, after inoculation with BdAsia-1, BdGPL or a blank solution. Subjects became infected in all experimental treatments but Korean species rapidly cleared themselves of infection, regardless of Bd lineage. They survived with no apparent secondary effects. By contrast, L. caerulea, after infection by either BdAsia-1 or BdGPL, suffered deteriorating body condition and carried progressively higher Bd loads over time. Subsequently, most subjects died. Comparing their effects on L. caerulea, BdAsia-1 induced more rapid disease progression than BdGPL. The results suggest that genomic recombination with other lineages was not necessary for the ancestral Bd lineage to evolve hypervirulence over its long period of coevolution with amphibian hosts. The pathogen's virulence may have driven strong selection for immune responses in endemic Asian amphibian host species.

Changes in soil taxonomic and functional diversity resulting from gamma irradiation.

Little is known of the effects of ionizing radiation exposure on soil biota. We exposed soil microcosms to weekly bursts of 60Co gamma radiation over six weeks, at three levels of exposure (0.1 kGy/hr/wk [low], 1 kGy/hr/wk [medium] and 3 kGy/hr/wk [high]). Soil DNA was extracted, and shotgun metagenomes were sequenced and characterised using MG-RAST. We hypothesized that with increasing radiation exposure there would be a decrease in both taxonomic and functional diversity. While bacterial diversity decreased, diversity of fungi and algae unexpectedly increased, perhaps because of release from competition. Despite the decrease in diversity of bacteria and of biota overall, functional gene diversity of algae, bacteria, fungi and total biota increased. Cycles of radiation exposure may increase the range of gene functional strategies viable in soil, a novel ecological example of the effects of stressors or disturbance events promoting some aspects of diversity. Moreover, repeated density-independent population crashes followed by population expansion may allow lottery effects, promoting coexistence. Radiation exposure produced large overall changes in community composition. Our study suggests several potential novel radiation-tolerant groups: in addition to Deinococcus-Thermus, which reached up to 20% relative abundance in the metagenome, the phyla Chloroflexi (bacteria), Chytridiomycota (fungi) and Nanoarcheota (archaea) may be considered as radiation-tolerant.

Community Ecology of Deinococcus in Irradiated Soil.

Deinococcus is a genus of soil bacteria known for radiation resistance. However, the effects of radiation exposure on its community structure are unknown. We exposed soil to three levels of gamma radiation, 0.1 kGy/h (low), 1 kGy/h (medium), and 3 kGy/h (high), once a week for 6 weeks and then extracted soil DNA for 16S rRNA amplicon sequencing. We found the following: (1) Increasing radiation dose produced a major increase in relative abundance of Deinococcus, reaching ~ 80% of reads at the highest doses. Differing abundances of the various Deinococcus species in relation to exposure levels indicate distinct "radiation niches." At 3 kGy/h, a single OTU identified as D. ficus overwhelmingly dominated the mesocosms. (2) Corresponding published genome data show that the dominant species at 3 kGy/h, D. ficus, has a larger and more complex genome than other Deinococcus species with a greater proportion of genes related to DNA and nucleotide metabolism, cell wall, membrane, and envelope biogenesis as well as more cell cycle control, cell division, and chromosome partitioning-related genes. Deinococcus ficus also has a higher guanine-cytosine ratio than most other Deinococcus. These features may be linked to genome stability and may explain its greater abundance in this apparently competitive system, under high-radiation exposures. (3) Genomic analysis suggests that Deinococcus, including D. ficus, are capable of utilizing diverse carbon sources derived from both microbial cells killed by the radiation (including C5-C12-containing compounds, like arabinose, lactose, N-acetyl-D-glucosamine) and plant-derived organic matter in the soil (e.g., cellulose and hemicellulose). (4) Overall, based on its metagenome, even the most highly irradiated (3 kGy/h) soil possesses a wide range of the activities necessary for a functional soil system. Future studies may consider the resilience and sustainability of such soils in a high-radiation environment.

Fungal Elevational Rapoport pattern from a High Mountain in Japan.

Little is known of how fungal distribution ranges vary with elevation. We studied fungal diversity and community composition from 740 to 2940 m above sea level on Mt. Norikura, Japan, sequencing the ITS2 region. There was a clear trend, repeated across each of the fungal phyla (Basidiomycota, Ascomycota, Zygomycota, Chytridomycota and Glomeromycota), and across the whole fungal community combined, towards an increased elevational range of higher elevation OTUs, conforming to the elevational Rapoport pattern. It appears that fungi from higher elevation environments are more generalized ecologically, at least in terms of climate-related gradients. These findings add to the picture from latitudinal studies of fungal ranges, which also suggest that the classic Rapoport Rule (broader ranges at higher latitudes) applies on a geographical scale. However, there was no mid-elevation maximum in diversity in any of the phyla studied, and different diversity trends for the different phyla, when different diversity indices were used. In terms of functional guilds, on Norikura there were trends towards increased saprotrophism (Zygomycota), symbiotrophism (Basidiomycota), symbiotrophism and saprotrophism (Ascomycota) and pathotrophism (Chytridiomycota) with elevation. The causes of each of these trends require further investigation from an ecological and evolutionary viewpoint.

Community richness of amphibian skin bacteria correlates with bioclimate at the global scale.

Animal-associated microbiomes are integral to host health, yet key biotic and abiotic factors that shape host-associated microbial communities at the global scale remain poorly understood. We investigated global patterns in amphibian skin bacterial communities, incorporating samples from 2,349 individuals representing 205 amphibian species across a broad biogeographic range. We analysed how biotic and abiotic factors correlate with skin microbial communities using multiple statistical approaches. Global amphibian skin bacterial richness was consistently correlated with temperature-associated factors. We found more diverse skin microbiomes in environments with colder winters and less stable thermal conditions compared with environments with warm winters and less annual temperature variation. We used bioinformatically predicted bacterial growth rates, dormancy genes and antibiotic synthesis genes, as well as inferred bacterial thermal growth optima to propose mechanistic hypotheses that may explain the observed patterns. We conclude that temporal and spatial characteristics of the host's macro-environment mediate microbial diversity.

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi.

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

Recent Asian origin of chytrid fungi causing global amphibian declines.

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

Skin Bacterial Community Reorganization Following Metamorphosis of the Fire-Bellied Toad (Bombina orientalis).

In organisms with complex life histories, dramatic changes in microbial community structure may occur with host development and immune system maturation. Amphibian host susceptibility to diseases such as chytridiomycosis may be affected by the reorganization of skin microbial community structure that occurs during metamorphosis. We tracked changes in the bacterial communities inhabiting skin of Korean fire-bellied toads (Bombina orientalis) that we infected as tadpoles with different strains of Batrachochytrium dendrobatidis (Bd), the pathogenic fungus that causes chytridiomycosis. We found that B. orientalis undergoes a major change in skin bacterial community composition between 5 and 15 days following metamorphosis. Richness indices and phylogenetic diversity measures began to diverge earlier, between aquatic and terrestrial stages. Our results further reveal differences in skin bacterial community composition among infection groups, suggesting that the effect of Bd infection on skin microbiome composition may differ by Bd strain. Additional studies are needed to further investigate the structural and temporal dynamics of microbiome shifts during metamorphosis in wild and captive amphibian populations. Analyses of the ontogeny of microbiome shifts may contribute to an understanding of why amphibians vary in their susceptibility to chytridiomycosis.

Major histocompatibility complex variation and the evolution of resistance to amphibian chytridiomycosis.

Chytridiomycosis, caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd), has been implicated in population declines and species extinctions of amphibians around the world. Susceptibility to the disease varies both within and among species, most likely attributable to heritable immunogenetic variation. Analyses of transcriptional expression in hosts following their infection by Bd reveal complex responses. Species resistant to Bd generally show evidence of stronger innate and adaptive immune system responses. Major histocompatibility complex (MHC) class I and class II genes of some susceptible species are up-regulated following host infection by Bd, but resistant species show no comparable changes in transcriptional expression. Bd-resistant species share similar pocket conformations within the MHC-II antigen-binding groove. Among susceptible species, survivors of epizootics bear alleles encoding these conformations. Individuals with homozygous resistance alleles appear to benefit by enhanced resistance, especially in environmental conditions that promote pathogen virulence. Subjects that are repeatedly infected and subsequently cleared of Bd can develop an acquired immune response to the pathogen. Strong directional selection for MHC alleles that encode resistance to Bd may deplete genetic variation necessary to respond to other pathogens. Resistance to chytridiomycosis incurs life-history costs that require further study.

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.