Host immune responses to enzootic and invasive pathogen lineages vary in magnitude, timing, and efficacy.

Infectious diseases of wildlife continue to pose a threat to biodiversity worldwide, yet pathogens are far from uniform in virulence or host disease outcome. Within the same pathogen species, virulence can vary considerably depending on strain or lineage, in turn eliciting variable host responses. One pathogen that has caused extensive biodiversity loss is the amphibian-killing fungus, Batrachochytrium dendrobatidis (Bd), which is comprised of a globally widespread hypervirulent lineage (Bd-GPL), and multiple geographically restricted, enzootic lineages. Whereas host immunogenomic responses to Bd-GPL have been characterized in a number of amphibian species, immunogenomic responses to geographically restricted, enzootic Bd lineages are less clear. To examine lineage-specific host immune responses to Bd, we exposed a species of pumpkin toadlet, Brachycephalus pitanga, which is endemic to Brazil's Southern Atlantic Forest, to either the Bd-GPL or the enzootic Bd-Asia-2/Brazil (hereafter Bd-Brazil) lineage. Using temporal samples from early, mid, and late infection stages, we quantified functional immunogenomic responses over the course of infection using differential gene expression tests and coexpression network analyses. Host immune responses varied significantly with Bd lineage. Relative to controls, toadlet responses to Bd-Brazil were weak at early infection (25 genes significantly differentially expressed), peaked by mid-stage infection (414 genes), and were nearly fully resolved by late-stage infection (nine genes). In contrast, responses to Bd-GPL were magnified and delayed; toadlets significantly differentially expressed 111 genes early, 87 genes at mid-stage infection, and 726 genes by late-stage infection relative to controls. Given that infection intensity did not vary between mid- and late-stage disease in either Bd-Brazil or Bd-GPL treatments, this suggests that pumpkin toadlets may be at least partially tolerant to the enzootic Bd-Brazil lineage. In contrast, late-stage immune activation against Bd-GPL was consistent with immune dysregulation previously observed in other species. Our results demonstrate that both the timing of immune response and the particular immune pathways activated are specific to Bd lineage. Within regions where multiple Bd lineages co-occur, and given continued global Bd movement, these differential host responses may influence not only individual disease outcome, but transmission dynamics at the population and community levels.

Cryptic lineages and standing genetic variation across independent cane toad introductions.

Widespread introduced species can be leveraged to investigate the genetic, ecological and adaptive processes underlying rapid evolution and range expansion, particularly the contributions of genetic diversity to adaptation. Rhinella marina, the cane toad, has been a focus of invasion biology for decades in Australia. However, their introduction history in North America is less clear. Here, we investigated the roles of introduction history and genetic diversity in establishment success of cane toads across their introduced range. We used reduced representation sequencing (ddRAD) to obtain 34,000 SNPs from 247 toads in native (French Guiana, Guyana, Ecuador, Panama, Texas) and introduced (Bermuda, southern Florida, northern Florida, Hawai'i, Puerto Rico) populations. Unlike all other cane toad introductions, we found that Florida populations were more closely related to native Central American lineages (R. horribilis), than to native Southern American lineages (R. marina). Furthermore, we found high levels of diversity and population structure in the native range, corroborating suggestions that R. marina is a species complex. We also found that introduced populations exhibit only slightly lower genetic diversity than native populations. Together with demographic analyses, this indicates founding populations of toads in Florida were larger than previously reported. Lastly, within R. marina, only one of 245 putatively adaptive SNPs showed fixed differences between native and introduced ranges, suggesting that putative selection in these introduced populations is based upon existing genetic variation. Our findings highlight the importance of genetic sequencing in understanding biological introductions and hint at the role of standing genetic variation in range expansion.

Temporal and spatial diversification along the Amazonia-Cerrado transition in Neotropical treefrogs of the Boana albopunctata species group.

Despite extensive research on biodiversity in Neotropical forests, biodiversity in seasonally dry, open biomes in South America has been underestimated until recently. We leverage a widespread group, Boana albopunctata, to uncover cryptic lineages and investigate the timing of diversification in Neotropical anurans with a focus on dry diagonal biomes (Cerrado, Caatinga and Chaco) and the ecotone between Amazonia and the Cerrado. We inferred a multilocus phylogeny of the B. albopunctata species group that includes 15 of 18 described species, recovered two cryptic species, and reconstructed the timing of diversification among species distributed across multiple South American biomes. One new potential species (B. aff. steinbachi), sampled in the Amazonian state of Acre, clustered within the B. calcara-fasciata species complex and is close to B. steinbachi. A second putative new species (B. aff. multifasciata), sampled in the Amazonia-Cerrado ecotone, is closely related to B. multifasciata. Lastly, we place a recently identified Cerrado lineage (B. aff. albopuncata) into the B. albopunctata species group phylogeny for the first time. Our ancestral range reconstruction showed that species in the B. albopuctata group likely dispersed from Amazonia-Cerrado into the dry-diagonal and Atlantic Forest. Intraspecies demography showed, for both B. raniceps and B. albopunctata, signs of rapid expansion across the dry diagonal. Similarly, for one clade of B. multifasciata, our analyses support an invasion of the Cerrado from Amazonia, followed by a rapid expansion across the open diagonal biomes. Thus, our study recovers several recent divergences along the Amazonia-Cerrado ecotone in northern Brazil. Tectonic uplift and erosion in the late Miocene and climate oscillations in the Pleistocene corresponded with estimated divergence times in the dry diagonal and Amazonia-Cerrado ecotone. Our study highlights the importance of these threatened open formations in the generation of biodiversity in the Neotropics.

Temperature dependence of metabolic rate in tropical and temperate aquatic insects: Support for the Climate Variability Hypothesis in mayflies but not stoneflies.

A fundamental gap in climate change vulnerability research is an understanding of the relative thermal sensitivity of ectotherms. Aquatic insects are vital to stream ecosystem function and biodiversity but insufficiently studied with respect to their thermal physiology. With global temperatures rising at an unprecedented rate, it is imperative that we know how aquatic insects respond to increasing temperature and whether these responses vary among taxa, latitudes, and elevations. We evaluated the thermal sensitivity of standard metabolic rate in stream-dwelling baetid mayflies and perlid stoneflies across a ~2,000 m elevation gradient in the temperate Rocky Mountains in Colorado, USA, and the tropical Andes in Napo, Ecuador. We used temperature-controlled water baths and microrespirometry to estimate changes in oxygen consumption. Tropical mayflies generally exhibited greater thermal sensitivity in metabolism compared to temperate mayflies; tropical mayfly metabolic rates increased more rapidly with temperature and the insects more frequently exhibited behavioral signs of thermal stress. By contrast, temperate and tropical stoneflies did not clearly differ. Varied responses to temperature among baetid mayflies and perlid stoneflies may reflect differences in evolutionary history or ecological roles as herbivores and predators, respectively. Our results show that there is physiological variation across elevations and species and that low-elevation tropical mayflies may be especially imperiled by climate warming. Given such variation among species, broad generalizations about the vulnerability of tropical ectotherms should be made more cautiously.

Narrow thermal tolerance and low dispersal drive higher speciation in tropical mountains.

Species richness is greatest in the tropics, and much of this diversity is concentrated in mountains. Janzen proposed that reduced seasonal temperature variation selects for narrower thermal tolerances and limited dispersal along tropical elevation gradients [Janzen DH (1967) Am Nat 101:233-249]. These locally adapted traits should, in turn, promote reproductive isolation and higher speciation rates in tropical mountains compared with temperate ones. Here, we show that tropical and temperate montane stream insects have diverged in thermal tolerance and dispersal capacity, two key traits that are drivers of isolation in montane populations. Tropical species in each of three insect clades have markedly narrower thermal tolerances and lower dispersal than temperate species, resulting in significantly greater population divergence, higher cryptic diversity, higher tropical speciation rates, and greater accumulation of species over time. Our study also indicates that tropical montane species, with narrower thermal tolerance and reduced dispersal ability, will be especially vulnerable to rapid climate change.

Incapacitating effects of fungal coinfection in a novel pathogen system.

As globalization lowers geographic barriers to movement, coinfection with novel and enzootic pathogens is increasingly likely. Novel and enzootic pathogens can interact synergistically or antagonistically, leading to increased or decreased disease severity. Here we examine host immune responses to coinfection with two closely related fungal pathogens: Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal). Both pathogens have had detrimental effects on amphibian populations, with Bd now largely enzootic, while Bsal is currently spreading and causing epizootics. Recent experimental work revealed that newts coinfected with Bd and Bsal had significantly higher mortality than those infected with either pathogen alone. Here we characterize host immunogenomic responses to chytrid coinfection relative to single infection. Across several classes of immune genes including pattern recognition receptors, cytokines, and MHC, coinfected host gene expression was weakly upregulated or comparable to that seen in single Bd infection, but significantly decreased when compared to Bsal infection. Combined with strong complement pathway downregulation and keratin upregulation, these results indicate that coinfection with Bd and Bsal compromises immune responses active against Bsal alone. As Bsal continues to invade naïve habitats where Bd is enzootic, coinfection will be increasingly common. If other Bd-susceptible species in the region have similar responses, interactions between the two pathogens could cause severe population and community-level declines.

Close-kin mating, but not inbred parents, reduces hatching rates and offspring quality in a threatened tortoise.

Inbreeding depression, the reduction in fitness due to mating of related individuals, is of particular conservation concern in species with small, isolated populations. Although inbreeding depression is widespread in natural populations, long-lived species may be buffered from its effects during population declines due to long generation times and thus are less likely to have evolved mechanisms of inbreeding avoidance than species with shorter generation times. However, empirical evidence of the consequences of inbreeding in threatened, long-lived species is limited. In this study, we leverage a well-studied population of gopher tortoises, Gopherus polyphemus, to examine the role of inbreeding depression and the potential for behavioural inbreeding avoidance in a natural population of a long-lived species. We tested the hypothesis that increased parental inbreeding leads to reduced hatching rates and offspring quality. Additionally, we tested for evidence of inbreeding avoidance. We found that high parental relatedness results in offspring with lower quality and that high parental relatedness is correlated with reduced hatching success. However, we found that hatching success and offspring quality increase with maternal inbreeding, likely due to highly inbred females mating with more distantly related males. We did not find evidence for inbreeding avoidance in males and outbred females, suggesting sex-specific evolutionary trade-offs may have driven the evolution of mating behaviour. Our results demonstrate inbreeding depression in a long-lived species and that the evolution of inbreeding avoidance is shaped by multiple selective forces.

Diverse genotypes of the amphibian-killing fungus produce distinct phenotypes through plastic responses to temperature.

Phenotypes are the target of selection and affect the ability of organisms to persist in variable environments. Phenotypes can be influenced directly by genes and/or by phenotypic plasticity. The amphibian-killing fungus Batrachochytrium dendrobatidis (Bd) has a global distribution, unusually broad host range, and high genetic diversity. Phenotypic plasticity may be an important process that allows this pathogen to infect hundreds of species in diverse environments. We quantified phenotypic variation of nine Bd genotypes from two Bd lineages (Global Pandemic Lineage [GPL] and Brazil) and a hybrid (GPL-Brazil) grown at three temperatures (12, 18 and 24°C). We measured five functional traits including two morphological traits (zoospore and zoosporangium sizes) and three life history traits (carrying capacity, time to fastest growth and exponential growth rate) in a phylogenetic framework. Temperature caused highly plastic responses within each genotype, with all Bd genotypes showing phenotypic plasticity in at least three traits. Among genotypes, Bd generally showed the same direction of plastic response to temperature: larger zoosporangia, higher carrying capacity, longer time to fastest growth and slower exponential growth at lower temperatures. The exception was zoospore size, which was highly variable. Our findings indicate that Bd genotypes have evolved novel phenotypes through plastic responses to temperature over very short timescales. High phenotypic variability likely extends to other traits and may facilitate the large host range and rapid spread of Bd.

Phenotypes in phylogeography: Species' traits, environmental variation, and vertebrate diversification.

Almost 30 y ago, the field of intraspecific phylogeography laid the foundation for spatially explicit and genealogically informed studies of population divergence. With new methods and markers, the focus in phylogeography shifted to previously unrecognized geographic genetic variation, thus reducing the attention paid to phenotypic variation in those same diverging lineages. Although phenotypic differences among lineages once provided the main data for studies of evolutionary change, the mechanisms shaping phenotypic differentiation and their integration with intraspecific genetic structure have been underexplored in phylogeographic studies. However, phenotypes are targets of selection and play important roles in species performance, recognition, and diversification. Here, we focus on three questions. First, how can phenotypes elucidate mechanisms underlying concordant or idiosyncratic responses of vertebrate species evolving in shared landscapes? Second, what mechanisms underlie the concordance or discordance of phenotypic and phylogeographic differentiation? Third, how can phylogeography contribute to our understanding of functional phenotypic evolution? We demonstrate that the integration of phenotypic data extends the reach of phylogeography to explain the origin and maintenance of biodiversity. Finally, we stress the importance of natural history collections as sources of high-quality phenotypic data that span temporal and spatial axes.

Extreme streams: species persistence and genomic change in montane insect populations across a flooding gradient.

The ecological and evolutionary consequences of extreme events are poorly understood. Here, we tested predictions about species persistence and population genomic change in aquatic insects in 14 Colorado mountain streams across a hydrological disturbance gradient caused by a one in 500-year rainfall event. Taxa persistence ranged from 39 to 77% across sites and declined with increasing disturbance in relation to species' resistance and resilience traits. For taxa with mobile larvae and terrestrial adult stages present at the time of the flood, average persistence was 84% compared to 25% for immobile taxa that lacked terrestrial adults. For two of six species analysed, genomic diversity (allelic richness) declined after the event. For one species it greatly expanded, suggesting resilience via re-colonisation from upstream populations. Thus, while resistance and resilience traits can explain species persistence to extreme disturbance, population genomic change varies among species, challenging generalisations about evolutionary responses to extreme events at landscape scales.

Globally invasive genotypes of the amphibian chytrid outcompete an enzootic lineage in coinfections.

Competition between genotypes is likely to be a key driver of pathogen evolution, particularly following a geographical invasion by distant strains. Theory predicts that competition between disease strains will result in the most virulent strain persisting. Despite its evolutionary implications, the role of strain competition in shaping populations remains untested for most pathogens. We experimentally investigated the in vivo competitive differences between two divergent lineages of the amphibian-killing chytrid fungus ( Batrachochytrium dendrobatidis, Bd). These Bd lineages are hypothesized to have diverged in allopatry but been recently brought back into secondary contact by human introduction. Prior studies indicate that a panzootically-distributed, global lineage of Bd was recently introduced into southern Brazil, and is competitively excluding enzootic lineages in the southern Atlantic Forest. To test for differences in competitive ability between invasive and enzootic Brazilian Bd isolates, we coinfected a model host frog system which we developed for this study ( Hymenochirus curtipes). We tracked isolate-specific zoospore production over the course of the coinfection experiment with chip-based digital PCR (dPCR). The globally invasive panzootic lineage had a competitive advantage in spore production especially during the first one to four weeks of infection, and on frogs that eventually succumbed to Bd infection. Our study provides new evidence that competitive pressure resulting from the human movement of pathogen strains can rapidly alter the genetics, community dynamics and spatial epidemiology of pathogens in the wild.

Molecular phylogeny of Neotropical rock frogs reveals a long history of vicariant diversification in the Atlantic forest.

The Brazilian Atlantic coastal forest is one of the most heterogeneous morphoclimatic domains on earth and is thus an excellent region in which to examine the role that habitat heterogeneity plays in shaping diversification of lineages and species. Here we present a molecular phylogeny of the rock frogs of the genus Thoropa Cope, 1865, native to the Atlantic forest and extending to adjacent campo rupestre of Brazil. The goal of this study is to reconstruct the evolutionary history of the genus using multilocus molecular phylogenetic analyses. Our topology reveals 12 highly supported lineages among the four nominal species included in the study. Species T. saxatilis and T. megatympanum are monophyletic. Thoropa taophora is also monophyletic, but nested within T. miliaris. Populations of T. miliaris cluster in five geographically distinct lineages, with low support for relationships among them. Although all 12 lineages are geographically structured, some T. miliaris lineages have syntopic distributions with others, likely reflecting a secondary contact zone between divergent lineages. We discuss a biogeographic scenario that best explains the order of divergence and the distribution of species in Atlantic forest and adjacent areas, and outline the implications of our findings for the taxonomy of Thoropa.

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.

Eco-evolutionary rescue promotes host-pathogen coexistence.

Emerging infectious pathogens are responsible for some of the most severe host mass mortality events in wild populations. Yet, effective pathogen control strategies are notoriously difficult to identify, in part because quantifying and forecasting pathogen spread and disease dynamics is challenging. Following an outbreak, hosts must cope with the presence of the pathogen, leading to host-pathogen coexistence or extirpation. Despite decades of research, little is known about host-pathogen coexistence post-outbreak when low host abundances and cryptic species make these interactions difficult to study. Using a novel disease-structured N-mixture model, we evaluate empirical support for three host-pathogen coexistence hypotheses (source-sink, eco-evolutionary rescue, and spatial variation in pathogen transmission) in a Neotropical amphibian community decimated by Batrachochytrium dendrobatidis (Bd) in 2004. During 2010-2014, we surveyed amphibians in Parque Nacional G. D. Omar Torríjos Herrera, Coclé Province, El Copé, Panama. We found that the primary driver of host-pathogen coexistence was eco-evolutionary rescue, as evidenced by similar amphibian survival and recruitment rates between infected and uninfected hosts. Average apparent monthly survival rates of uninfected and infected hosts were both close to 96%, and the expected number of uninfected and infected hosts recruited (via immigration/reproduction) was less than one host per disease state per 20-m site. The secondary driver of host-pathogen coexistence was spatial variation in pathogen transmission as we found that transmission was highest in areas of low abundance but there was no support for the source-sink hypothesis. Our results indicate that changes in the host community (i.e., through genetic or species composition) can reduce the impacts of emerging infectious disease post-outbreak. Our disease-structured N-mixture model represents a valuable advancement for conservation managers trying to understand underlying host-pathogen interactions and provides new opportunities to study disease dynamics in remnant host populations decimated by virulent pathogens.

Male Body Size Predicts Reproductive Success But Not Within-Clutch Paternity Patterns in Gopher Tortoises (Gopherus polyphemus).

In many vertebrates, body size is an important driver of variation in male reproductive success. Larger, more fit individuals are more likely to dominate mating opportunities, skewing siring success and resulting in lower effective population sizes and genetic diversity. The mating system of the gopher tortoise (Gopherus polyphemus) has been characterized as both female-defense and scramble-competition polygyny. Mating systems are typically not fixed and can be influenced by factors such as population density, demographic structure, and environmental conditions; however, most populations will have a predominant strategy that results from local conditions. We assessed how male body size influences patterns of paternity and reproductive success in a natural population of gopher tortoises in Florida, United States. Using microsatellites, we assigned parentage of 220 hatchlings from 31 nests collected during 2 reproductive seasons. Larger males were significantly more likely to sire offspring and sired more offspring than smaller males; however, the likelihood of a clutch being multiply sired was unrelated to male body size. We also found evidence of mate fidelity across years. Although paternity patterns in this high-density population are more consistent with defense polygyny, female monopoly by males was incomplete, with both large and small males contributing to multiply sired clutches. Additional behavioral data are needed to clarify the role of female mate selection in paternity outcomes. The context-dependence of mating systems underscores the need to compare parentage patterns across populations and to recognize the potential for more than 1 strategy to be employed within a single population.

Temperature variation, bacterial diversity and fungal infection dynamics in the amphibian skin.

Host-associated bacterial communities on the skin act as the first line of defence against invading pathogens. Yet, for most natural systems, we lack a clear understanding of how temperature variability affects structure and composition of skin bacterial communities and, in turn, promotes or limits the colonization of opportunistic pathogens. Here, we examine how natural temperature fluctuations might be related to changes in skin bacterial diversity over time in three amphibian populations infected by the pathogenic fungus Batrachochytrium dendrobatidis (Bd). Our focal host species (Eleutherodactylus coqui) is a direct-developing frog that has suffered declines at some populations in the last 20 years, while others have not experienced any changes. We quantified skin bacterial alpha- and beta-diversity at four sampling time points, a period encompassing two seasons and ample variation in natural infections and environmental conditions. Despite the different patterns of infection across populations, we detected an overall increase in bacterial diversity through time, characterized by the replacement of bacterial operational taxonomic units (OTUs). Increased frog body temperatures possibly allowed the colonization of bacteria as well as the recruitment of a subset of indicator OTUs, which could have promoted the observed changes in diversity patterns. Our results suggest that natural environmental fluctuations might be involved in creating opportunities for bacterial replacement, potentially attenuating pathogen transmission and thus contributing to host persistence in E. coqui populations.

Idiosyncratic responses to climate-driven forest fragmentation and marine incursions in reed frogs from Central Africa and the Gulf of Guinea Islands.

Organismal traits interact with environmental variation to mediate how species respond to shared landscapes. Thus, differences in traits related to dispersal ability or physiological tolerance may result in phylogeographic discordance among co-distributed taxa, even when they are responding to common barriers. We quantified climatic suitability and stability, and phylogeographic divergence within three reed frog species complexes across the Guineo-Congolian forests and Gulf of Guinea archipelago of Central Africa to investigate how they responded to a shared climatic and geological history. Our species-specific estimates of climatic suitability through time are consistent with temporal and spatial heterogeneity in diversification among the species complexes, indicating that differences in ecological breadth may partly explain these idiosyncratic patterns. Likewise, we demonstrated that fluctuating sea levels periodically exposed a land bridge connecting Bioko Island with the mainland Guineo-Congolian forest and that habitats across the exposed land bridge likely enabled dispersal in some species, but not in others. We did not find evidence that rivers are biogeographic barriers across any of the species complexes. Despite marked differences in the geographic extent of stable climates and temporal estimates of divergence among the species complexes, we recovered a shared pattern of intermittent climatic suitability with recent population connectivity and demographic expansion across the Congo Basin. This pattern supports the hypothesis that genetic exchange across the Congo Basin during humid periods, followed by vicariance during arid periods, has shaped regional diversity. Finally, we identified many distinct lineages among our focal taxa, some of which may reflect incipient or unrecognized species.

Polyandry, Predation, and the Evolution of Frog Reproductive Modes.

Frog reproductive modes are complex phenotypes that include egg/clutch characteristics, oviposition site, larval development, and sometimes, parental care. Two evident patterns in the evolution of these traits are the higher diversity of reproductive modes in the tropics and the apparent progression from aquatic to terrestrial reproduction, often attributed to higher fitness resulting from decreased predation on terrestrial eggs and tadpoles. Here, we propose that sexual selection-and not only natural selection due to predation-favors terrestrial breeding by reducing the loss of fitness due to polyandry. To examine this novel selective mechanism, we reconstructed the evolution of reproductive diversity in two frog families (Hylidae and Leptodactylidae) and tested for concerted evolution of egg and tadpole development sites with specific mating behaviors. We found that oviposition and tadpole development sites are evolving independently, do not show the same diversity and/or directionality in terms of terrestriality, and thus may be diversifying due to different selective mechanisms. In both families, terrestrial egg deposition is correlated with amplexus that is hidden from competing males, and in hylids, testes mass was significantly larger and more variable in males with exposed amplexus that are vulnerable to polyandry. Our results indicate that intrasexual selection has been an underappreciated mechanism promoting diversification of frog reproductive modes.

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.