The Global Panzootic Lineage of the fungus Batrachochytrium dendrobatidis (Bd-GPL) is threatening amphibians worldwide. In contrast, four lineages (Bd-Brazil, Bd-CH, Bd-Cape, and Bd-Korea) that diverged early in the history of Bd have not yet been directly linked to amphibian declines. Bd likely evolves in response to strong selective pressure imposed by hosts and the environment, leading to differences among pathogen phenotypes and genotypes that may directly affect virulence. Here, we report on variation in phenotype, genotype, and virulence of Bd-Brazil and Bd-GPL. Specifically, we (i) used a controlled infection experiment to compare virulence between one Bd-Brazil and three Bd-GPL isolates on a North American amphibian host (Lithobates sylvaticus), (ii) tested for relative phenotypic and genotypic differentiation among Bd isolates from Brazil, and (iii) tested for possible correlations between environmental variables and Bd phenotypes. We found substantial variation in virulence among Bd-GPL isolates and found that our Bd-Brazil isolate showed virulence comparable to an average North American Bd-GPL. North American hosts infected with a Bd-GPL isolate from Panama did not show significant mortality. Bd phenotypes varied significantly across sampling locations; these phenotypes were neither spatially clustered nor correlated with any environmental variables. Additionally, we found a surprising lack of correlation between genotypic divergence and zoospore and zoosporangium sizes in our sample. Although Bd-Brazil was less virulent, infecting L. sylvaticus than one Bd-GPL isolate, this endemic lineage still caused ~50% mortality in our experimental North American hosts. This indicates that Bd-Brazil has the potential to kill amphibians if introduced to naïve wild populations. Our findings underscore that characterizing virulence of multiple Bd isolates and lineages is important for understanding the evolutionary history and diversity of Bd.

Predators can alter nutrient cycles simply by inducing stress in prey. This stress accelerates prey's protein catabolism, nitrogen waste production, and nitrogen cycling. Yet predators also reduce the feeding rates of their prey, inducing food deprivation that is expected to slow protein catabolism and nitrogen cycling. The physiology of prey under predation risk thus balances the influences of predation risk and food deprivation, and this balance is central to understanding the role of predators in nutrient cycles. We explored the separate and combined effects of predation risk and food deprivation on prey physiology and nutrient cycling by exposing guppies (Poecilia reticulata) to predation risk and food deprivation in a 2 × 2 design. We simulated predation risk using chemical cues from a natural predator of guppies, and we created food deprivation by rationing food availability. We measured guppy response as food consumption, growth, tissue energy density, tissue carbon:nitrogen, and nitrogen (N) excretion and assimilation. We found that N-linked physiological processes (N consumption, assimilation, excretion) were strongly affected by predation risk, independent of food consumption. Guppies excreted substantially less under predation risk than they did under food deprivation or control conditions. These results suggest that predation risk, per se, triggers physiological changes in guppies that increase N retention and decrease N excretion. We suggest that slower N metabolism under predation risk is an adaptive response that minimizes protein loss in the face of predictable, predator-induced food restriction. Notably, N metabolism shares common hormonal control with food seeking behavior, and we speculate that increased N retention is a direct and immediate result of reduced food seeking under predation risk. Contrary to predation-stress-based hypotheses for how predators affect nutrient cycling by prey, our result indicates that even short-term exposure to predators may decelerate, rather than accelerate, the speed of N cycling by suppressing N turnover by prey.

Immune gene evolution can be critical to species survival in the face of infectious disease. In particular, polymorphism in the genes of the major histocompatibility complex (MHC) helps vertebrates combat novel and diverse pathogens by increasing the number of pathogen-derived proteins that can initiate the host's acquired immune response. In this study, we used a combination of presumably adaptive and neutral markers to investigate MHC evolution in populations of five salamander species within the Ambystoma velasci complex, a group consisting of 15 recently diverged species, several of which are endangered. We isolated 31 unique MHC class II β alleles from 75 total individuals from five species in this complex. MHC heterozygosity was significantly lower than expected for all five species, and we found no clear relationship between number of MHC alleles and species range, life history, or level of heterozygosity. We inferred a phylogeny representing the evolutionary history of Ambystoma MHC, with which we found signatures of positive selection on the overall gene, putative peptide-binding residues, and allelic lineages. We identified several instances of trans-species polymorphism, a hallmark of balancing selection observed in other groups of closely related species. In contrast, we did not detect comparable allelic diversity or signatures of selection on neutral loci. Additionally, we identified 17 supertypes among the 44 unique Ambystoma alleles, indicating that these sequences may encode functionally distinct MHC variants. We therefore have strong evidence that positive selection is a major evolutionary force driving patterns of MHC polymorphism in this recently radiated species complex.

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.