Sequence capture identifies fastidious chytrid fungi directly from host tissue.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches.

Reduction of genetic diversity in 'Alala (Hawaiian crow; Corvus hawaiiensis) between the late 1800s and the late 1900s.

Genetic and genomic data are increasingly used to aid conservation management of endangered species by providing insights into evolutionary histories, factors associated with extinction risks, and potential for future adaptation. For the 'Alala, or Hawaiian crow (Corvus hawaiiensis), genetic concerns include negative correlations between inbreeding and hatching success. However, it is unclear if low genetic diversity and inbreeding depression are consequences of a historical population bottleneck, or if 'Alala had historically low genetic diversity that predated human influence, perhaps as a result of earlier declines or founding events. In this study, we applied a hybridization-based sequence capture to generate a genome-wide single nucleotide polymorphism (SNP) dataset for comparing historical specimens collected in the 1890s, when 'Alala were more numerous, to samples taken between 1973 and 1998, when 'Alala population densities were near the lowest documented levels in the wild, prior to all individuals being collected for captive rearing. We found low genome-wide diversity in both sample groups, however, the modern sample group (1973 to 1998 cohort) exhibited relatively fewer polymorphic alleles, a lower proportion of polymorphic loci, and lower observed heterozygosity, consistent with a population decline and potential bottleneck effects. These results combined with a current low population size highlight the importance of continued efforts by conservation managers to mitigate inbreeding and maintain founder representation to preserve what genetic diversity remains.

Microbiomes associated with avian malaria survival differ between susceptible Hawaiian honeycreepers and sympatric malaria-resistant introduced birds.

Of the estimated 55 Hawaiian honeycreepers (subfamily Carduelinae) only 17 species remain, nine of which the International Union for Conservation of Nature considers endangered. Among the most pressing threats to honeycreeper survival is avian malaria, caused by the introduced blood parasite Plasmodium relictum, which is increasing in distribution in Hawai'i as a result of climate change. Preventing further honeycreeper decline will require innovative conservation strategies that confront malaria from multiple angles. Research on mammals has revealed strong connections between gut microbiome composition and malaria susceptibility, illuminating a potential novel approach to malaria control through the manipulation of gut microbiota. One honeycreeper species, Hawai'i 'amakihi (Chlorodrepanis virens), persists in areas of high malaria prevalence, indicating they have acquired some level of immunity. To investigate if avian host-specific microbes may be associated with malaria survival, we characterized cloacal microbiomes and malaria infection for 174 'amakihi and 172 malaria-resistant warbling white-eyes (Zosterops japonicus) from Hawai'i Island using 16S rRNA gene metabarcoding and quantitative polymerase chain reaction. Neither microbial alpha nor beta diversity covaried with infection, but 149 microbes showed positive associations with malaria survivors. Among these were Escherichia and Lactobacillus spp., which appear to mitigate malaria severity in mammalian hosts, revealing promising candidates for future probiotic research for augmenting malaria immunity in sensitive endangered species.

Gene expression reveals immune response strategies of naïve Hawaiian honeycreepers experimentally infected with introduced avian malaria.

The unprecedented rise in the number of new and emerging infectious diseases in the last quarter century poses direct threats to human and wildlife health. The introduction to the Hawaiian archipelago of Plasmodium relictum and the mosquito vector that transmits the parasite has led to dramatic losses in endemic Hawaiian forest bird species. Understanding how mechanisms of disease immunity to avian malaria may evolve is critical as climate change facilitates increased disease transmission to high elevation habitats where malaria transmission has historically been low and the majority of the remaining extant Hawaiian forest bird species now reside. Here, we compare the transcriptomic profiles of highly susceptible Hawai'i 'amakihi (Chlorodrepanis virens) experimentally infected with P. relictum to those of uninfected control birds from a naïve high elevation population. We examined changes in gene expression profiles at different stages of infection to provide an in-depth characterization of the molecular pathways contributing to survival or mortality in these birds. We show that the timing and magnitude of the innate and adaptive immune response differed substantially between individuals that survived and those that succumbed to infection, and likely contributed to the observed variation in survival. These results lay the foundation for developing gene-based conservation strategies for Hawaiian honeycreepers by identifying candidate genes and cellular pathways involved in the pathogen response that correlate with a bird's ability to recover from malaria infection.

Conservation genomics and systematics of a near-extinct island radiation.

Conservation benefits from incorporating genomics to explore the impacts of population declines, inbreeding, loss of genetic variation and hybridization. Here we use the near-extinct Mariana Islands reedwarbler radiation to showcase how ancient DNA approaches can allow insights into the population dynamics of extinct species and threatened populations for which historical museum specimens or material with low DNA yield (e.g., scats, feathers) are the only sources for DNA. Despite their having paraphyletic mitochondrial DNA (mtDNA), nuclear single nucleotide polymorphisms (SNPs) support the distinctiveness of critically endangered Acrocephalus hiwae and the other three species in the radiation that went extinct between the 1960s and 1990s. Two extinct species, A. yamashinae and A. luscinius, were deeply divergent from each other and from a third less differentiated lineage containing A. hiwae and extinct A. nijoi. Both mtDNA and SNPs suggest that the two isolated populations of A. hiwae from Saipan and Alamagan Islands are sufficiently distinct to warrant subspecies recognition and separate conservation management. We detected no significant differences in genetic diversity or inbreeding between Saipan and Alamagan, nor strong signatures of geographical structuring within either island. However, the implications of possible signatures of inbreeding in both Saipan and Alamagan, and long-term population declines in A. hiwae that pre-date modern anthropogenic threats require further study with denser population sampling. Our study highlights the value that conservation genomics studies of island radiations have as windows onto the possible future for the world's biota as climate change and habitat destruction increasingly fragment their ranges and contribute to rapid declines in population abundances.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra.

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

Transcriptional response of individual Hawaiian Culex quinquefasciatus mosquitoes to the avian malaria parasite Plasmodium relictum.

BACKGROUND: Plasmodium parasites that cause bird malaria occur in all continents except Antarctica and are primarily transmitted by mosquitoes in the genus Culex. Culex quinquefasciatus, the mosquito vector of avian malaria in Hawai'i, became established in the islands in the 1820s. While the deadly effects of malaria on endemic bird species have been documented for many decades, vector-parasite interactions in avian malaria systems are relatively understudied. METHODS: To evaluate the gene expression response of mosquitoes exposed to a Plasmodium infection intensity known to occur naturally in Hawai'i, offspring of wild-collected Hawaiian Cx. quinquefasciatus were fed on a domestic canary infected with a fresh isolate of Plasmodium relictum GRW4 from a wild-caught Hawaiian honeycreeper. Control mosquitoes were fed on an uninfected canary. Transcriptomes of five infected and three uninfected individual mosquitoes were sequenced at each of three stages of the parasite life cycle: 24 h post feeding (hpf) during ookinete invasion; 5 days post feeding (dpf) when oocysts are developing; 10 dpf when sporozoites are released and invade the salivary glands. RESULTS: Differential gene expression analyses showed that during ookinete invasion (24 hpf), genes related to oxidoreductase activity and galactose catabolism had lower expression levels in infected mosquitoes compared to controls. Oocyst development (5 dpf) was associated with reduced expression of a gene with a predicted innate immune function. At 10 dpf, infected mosquitoes had reduced expression levels of a serine protease inhibitor, and further studies should assess its role as a Plasmodium agonist in C. quinquefasciatus. Overall, the differential gene expression response of Hawaiian Culex exposed to a Plasmodium infection intensity known to occur naturally in Hawai'i was low, but more pronounced during ookinete invasion. CONCLUSIONS: This is the first analysis of the transcriptional responses of vectors to malaria parasites in non-mammalian systems. Interestingly, few similarities were found between the response of Culex infected with a bird Plasmodium and those reported in Anopheles infected with human Plasmodium. The relatively small transcriptional changes observed in mosquito genes related to immune response and nutrient metabolism support conclusions of low fitness costs often documented in experimental challenges of Culex with avian Plasmodium.

A genome-wide investigation of adaptive signatures in protein-coding genes related to tool behaviour in New Caledonian and Hawaiian crows.

Very few animals habitually manufacture and use tools. It has been suggested that advanced tool behaviour co-evolves with a suite of behavioural, morphological and life history traits. In fact, there are indications for such an adaptive complex in tool-using crows (genus Corvus species). Here, we sequenced the genomes of two habitually tool-using and ten non-tool-using crow species to search for genomic signatures associated with a tool-using lifestyle. Using comparative genomic and population genetic approaches, we screened for signals of selection in protein-coding genes in the tool-using New Caledonian and Hawaiian crows. While we detected signals of recent selection in New Caledonian crows near genes associated with bill morphology, our data indicate that genetic changes in these two lineages are surprisingly subtle, with little evidence at present for convergence. We explore the biological explanations for these findings, such as the relative roles of gene regulation and protein-coding changes, as well as the possibility that statistical power to detect selection in recently diverged lineages may have been insufficient. Our study contributes to a growing body of literature aiming to decipher the genetic basis of recently evolved complex behaviour.

Sustained immune activation is associated with susceptibility to the amphibian chytrid fungus.

The disease chytridiomycosis caused by the fungus Bd has devastated amphibian populations worldwide. Functional genomic contributions to host susceptibility remain enigmatic and vary between species and populations. We conducted experimental Bd infections in Rana yavapaiensis, a species with intraspecific variation in chytridiomycosis susceptibility, to assess the skin and spleen transcriptomic response to infection over time. We predicted that increased immune gene expression would be associated with a positive disease outcome, but we instead found that surviving frogs had significantly reduced immune gene expression compared to susceptible frogs and to uninfected controls. MHC class IIß gene expression was also significantly higher in susceptible frogs compared to surviving frogs. Furthermore, susceptible frogs expressed a significantly larger number of distinct class IIß alleles, demonstrating a negative correlation between class IIß expression, functional diversity, and survival. Expression of the MHC class IIß locus previously associated with Bd disease outcomes was a significant predictor of Bd infection intensity at early infection stages but not at late infection stages, suggesting initial MHC-linked immune processes are important for ultimate disease outcomes. We infer through disease association and phylogenetic analysis that certain MHC variants are linked to the immune expression that was negatively associated with survival, and we hypothesize that frogs that did not express these alleles could better survive infections. Our study finds that MHC expression at early and late infection stages predicts Bd infection intensity, and suggests that generating a sustained immune response against Bd may be counterproductive for surviving chytridiomycosis in this partially susceptible species.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks.

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

Adaptive Radiation Genomics of Two Ecologically Divergent Hawai'ian Honeycreepers: The 'akiapola'au and the Hawai'i 'amakihi.

The Hawai'ian honeycreepers (drepanids) are a classic example of adaptive radiation: they adapted to a variety of novel dietary niches, evolving a wide range of bill morphologies. Here we investigated genomic diversity, demographic history, and genes involved in bill morphology phenotypes in 2 honeycreepers: the 'akiapola'au (Hemignathus wilsoni) and the Hawai'i 'amakihi (Chlorodrepanis virens). The 'akiapola'au is an endangered island endemic, filling the "woodpecker" niche by using a unique bill morphology, while the Hawai'i 'amakihi is a dietary generalist common on the islands of Hawai'i and Maui. We de novo sequenced the 'akiapola'au genome and compared it to the previously sequenced 'amakihi genome. The 'akiapola'au is far less heterozygous and has a smaller effective population size than the 'amakihi, which matches expectations due to its smaller census population and restricted ecological niche. Our investigation revealed genomic islands of divergence, which may be involved in the honeycreeper radiation. Within these islands of divergence, we identified candidate genes (including DLK1, FOXB1, KIF6, MAML3, PHF20, RBP1, and TIMM17A) that may play a role in honeycreeper adaptations. The gene DLK1, previously shown to influence Darwin's finch bill size, may be related to honeycreeper bill morphology evolution, while the functions of the other candidates remain unknown.

Comparing Adaptive Radiations Across Space, Time, and Taxa.

Adaptive radiation plays a fundamental role in our understanding of the evolutionary process. However, the concept has provoked strong and differing opinions concerning its definition and nature among researchers studying a wide diversity of systems. Here, we take a broad view of what constitutes an adaptive radiation, and seek to find commonalities among disparate examples, ranging from plants to invertebrate and vertebrate animals, and remote islands to lakes and continents, to better understand processes shared across adaptive radiations. We surveyed many groups to evaluate factors considered important in a large variety of species radiations. In each of these studies, ecological opportunity of some form is identified as a prerequisite for adaptive radiation. However, evolvability, which can be enhanced by hybridization between distantly related species, may play a role in seeding entire radiations. Within radiations, the processes that lead to speciation depend largely on (1) whether the primary drivers of ecological shifts are (a) external to the membership of the radiation itself (mostly divergent or disruptive ecological selection) or (b) due to competition within the radiation membership (interactions among members) subsequent to reproductive isolation in similar environments, and (2) the extent and timing of admixture. These differences translate into different patterns of species accumulation and subsequent patterns of diversity across an adaptive radiation. Adaptive radiations occur in an extraordinary diversity of different ways, and continue to provide rich data for a better understanding of the diversification of life.

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

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIß locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIß alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIß functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

Fungal disease and temperature alter skin microbiome structure in an experimental salamander system.

Pathogens compete with host microbiomes for space and resources. Their shared environment impacts pathogen-microbiome-host interactions, which can lead to variation in disease outcome. The skin microbiome of red-backed salamanders (Plethodon cinereus) can reduce infection by the pathogen Batrachochytrium dendrobatidis (Bd) at moderate infection loads, with high species richness and high abundance of competitors as putative mechanisms. However, it is unclear if the skin microbiome can reduce epizootic Bd loads across temperatures. We conducted a laboratory experiment to quantify skin microbiome and host responses (P. cinereus: n = 87) to Bd at mimicked epizootic loads across temperatures (13, 17 and 21°C). We quantified skin microbiomes using 16S rRNA gene metabarcoding and identified operational taxonomic units (OTUs) taxonomically similar to culturable bacteria known to kill Bd (anti-Bd OTUs). Prior to pathogen exposure, temperature changed the microbiome (OTU richness decreased by 12% and the abundance of anti-Bd OTUs increased by 18% per degree increase in temperature), but these changes were not predictive of disease outcome. After exposure, Bd changed the microbiome (OTU richness decreased by 0.1% and the abundance of anti-Bd OTUs increased by 0.2% per 1% increase in Bd load) and caused high host mortality across temperatures (35/45: 78%). Temperature indirectly impacted microbiome change and mortality through its direct effect on pathogen load. We did not find support for the microbiome impacting Bd load or host survival. Our research reveals complex host, pathogen, microbiome and environmental interactions to demonstrate that during epizootic events the microbiome will be unlikely to reduce pathogen invasion, even for putatively Bd-resistant species.

Parallel evolution of gene classes, but not genes: Evidence from Hawai'ian honeycreeper populations exposed to avian malaria.

Adaptation in nature is ubiquitous, yet characterizing its genomic basis is difficult because population demographics cause correlations with nonadaptive loci. Introduction events provide opportunities to observe adaptation over known spatial and temporal scales, facilitating the identification of genes involved in adaptation. The pathogen causing avian malaria, Plasmodium relictum, was introduced to Hawai'i in the 1930s and elicited extinctions and precipitous population declines in native honeycreepers. After a sharp initial population decline, the Hawai'i 'amakihi (Chlorodrepanis virens) has evolved tolerance to the parasite at low elevations where P. relictum exists, and can sustain infection without major fitness consequences. High-elevation, unexposed populations of 'amakihi display little to no tolerance. To explore the genomic basis of adaptation to P. relictum in low-elevation 'amakihi, we genotyped 125 'amakihi from the island of Hawai'i via hybridization capture to 40,000 oligonucleotide baits containing SNPs and used the reference 'amakihi genome to identify genes potentially under selection from malaria. We tested for outlier loci between low- and high-elevation population pairs and identified loci with signatures of selection within low-elevation populations. In some cases, genes commonly involved in the immune response (e.g., major histocompatibility complex) were associated with malaria presence in the population. We also detected several novel candidate loci that may be implicated in surviving malaria infection (e.g., beta-defensin, glycoproteins and interleukin-related genes). Our results suggest that rapid adaptation to pathogens may occur through changes in different immune genes, but in the same classes of genes, across populations.

Diversity and temporal dynamics of primate milk microbiomes.

Milk is inhabited by a community of bacteria and is one of the first postnatal sources of microbial exposure for mammalian young. Bacteria in breast milk may enhance immune development, improve intestinal health, and stimulate the gut-brain axis for infants. Variation in milk microbiome structure (e.g., operational taxonomic unit [OTU] diversity, community composition) may lead to different infant developmental outcomes. Milk microbiome structure may depend on evolutionary processes acting at the host species level and ecological processes occurring over lactation time, among others. We quantified milk microbiomes using 16S rRNA high-throughput sequencing for nine primate species and for six primate mothers sampled over lactation. Our data set included humans (Homo sapiens, Philippines and USA) and eight nonhuman primate species living in captivity (bonobo [Pan paniscus], chimpanzee [Pan troglodytes], western lowland gorilla [Gorilla gorilla gorilla], Bornean orangutan [Pongo pygmaeus], Sumatran orangutan [Pongo abelii], rhesus macaque [Macaca mulatta], owl monkey [Aotus nancymaae]) and in the wild (mantled howler monkey [Alouatta palliata]). For a subset of the data, we paired microbiome data with nutrient and hormone assay results to quantify the effect of milk chemistry on milk microbiomes. We detected a core primate milk microbiome of seven bacterial OTUs indicating a robust relationship between these bacteria and primate species. Milk microbiomes differed among primate species with rhesus macaques, humans and mantled howler monkeys having notably distinct milk microbiomes. Gross energy in milk from protein and fat explained some of the variations in microbiome composition among species. Microbiome composition changed in a predictable manner for three primate mothers over lactation time, suggesting that different bacterial communities may be selected for as the infant ages. Our results contribute to understanding ecological and evolutionary relationships between bacteria and primate hosts, which can have applied benefits for humans and endangered primates in our care.

Effects of host species and environment on the skin microbiome of Plethodontid salamanders.

The amphibian skin microbiome is recognized for its role in defence against pathogens, including the deadly fungal pathogen Batrachochytrium dendrobatidis (Bd). Yet, we have little understanding of evolutionary and ecological processes that structure these communities, especially for salamanders and closely related species. We investigated patterns in the distribution of bacterial communities on Plethodon salamander skin across host species and environments. Quantifying salamander skin microbiome structure contributes to our understanding of how host-associated bacteria are distributed across the landscape, among host species, and their putative relationship with disease. We characterized skin microbiome structure (alpha-diversity, beta-diversity and bacterial operational taxonomic unit [OTU] abundances) using 16S rRNA gene sequencing for co-occurring Plethodon salamander species (35 Plethodon cinereus, 17 Plethodon glutinosus, 10 Plethodon cylindraceus) at three localities to differentiate the effects of host species from environmental factors on the microbiome. We sampled the microbiome of P. cinereus along an elevational gradient (n = 50, 700-1,000 m a.s.l.) at one locality to determine whether elevation predicts microbiome structure. Finally, we quantified prevalence and abundance of putatively anti-Bd bacteria to determine if Bd-inhibitory bacteria are dominant microbiome members. Co-occurring salamanders had similar microbiome structure, but among sites salamanders had dissimilar microbiome structure for beta-diversity and abundance of 28 bacterial OTUs. We found that alpha-diversity increased with elevation, beta-diversity and the abundance of 17 bacterial OTUs changed with elevation (16 OTUs decreasing, 1 OTU increasing). We detected 11 putatively anti-Bd bacterial OTUs that were present on 90% of salamanders and made up an average relative abundance of 83% (SD ± 8.5) per salamander. All salamanders tested negative for Bd. We conclude that environment is more influential in shaping skin microbiome structure than host differences in these congeneric species, and suggest that environmental characteristics that covary with elevation influence microbiome structure. High prevalence and abundance of anti-Bd bacteria may contribute to low Bd levels in these populations of Plethodon salamanders.

White-Nose Syndrome Fungus in a 1918 Bat Specimen from France.

White-nose syndrome, first diagnosed in North America in 2006, causes mass deaths among bats in North America. We found the causative fungus, Pseudogymnoascus destructans, in a 1918 sample collected in Europe, where bats have now adapted to the fungus. These results are consistent with a Eurasian origin of the pathogen.

Interacting effects of wildlife loss and climate on ticks and tick-borne disease.

Both large-wildlife loss and climatic changes can independently influence the prevalence and distribution of zoonotic disease. Given growing evidence that wildlife loss often has stronger community-level effects in low-productivity areas, we hypothesized that these perturbations would have interactive effects on disease risk. We experimentally tested this hypothesis by measuring tick abundance and the prevalence of tick-borne pathogens (Coxiella burnetii and Rickettsia spp.) within long-term, size-selective, large-herbivore exclosures replicated across a precipitation gradient in East Africa. Total wildlife exclusion increased total tick abundance by 130% (mesic sites) to 225% (dry, low-productivity sites), demonstrating a significant interaction of defaunation and aridity on tick abundance. When differing degrees of exclusion were tested for a subset of months, total tick abundance increased from 170% (only mega-herbivores excluded) to 360% (all large wildlife excluded). Wildlife exclusion differentially affected the abundance of the three dominant tick species, and this effect varied strongly over time, likely due to differences among species in their host associations, seasonality, and other ecological characteristics. Pathogen prevalence did not differ across wildlife exclusion treatments, rainfall levels, or tick species, suggesting that exposure risk will respond to defaunation and climate change in proportion to total tick abundance. These findings demonstrate interacting effects of defaunation and aridity that increase disease risk, and they highlight the need to incorporate ecological context when predicting effects of wildlife loss on zoonotic disease dynamics.

Antifungal Bacteria on Woodland Salamander Skin Exhibit High Taxonomic Diversity and Geographic Variability.

Diverse bacteria inhabit amphibian skin; some of those bacteria inhibit growth of the fungal pathogen Batrachochytrium dendrobatidis Yet there has been no systematic survey of anti-B. dendrobatidis bacteria across localities, species, and elevations. This is important given geographic and taxonomic variations in amphibian susceptibility to B. dendrobatidis Our collection sites were at locations within the Appalachian Mountains where previous sampling had indicated low B. dendrobatidis prevalence. We determined the numbers and identities of anti-B. dendrobatidis bacteria on 61 Plethodon salamanders (37 P. cinereus, 15 P. glutinosus, 9 P. cylindraceus) via culturing methods and 16S rRNA gene sequencing. We sampled co-occurring species at three localities and sampled P. cinereus along an elevational gradient (700 to 1,000 meters above sea level [masl]) at one locality. We identified 50 anti-B. dendrobatidis bacterial operational taxonomic units (OTUs) and found that the degree of B. dendrobatidis inhibition was not correlated with relatedness. Five anti-B. dendrobatidis bacterial strains occurred on multiple amphibian species at multiple localities, but none were shared among all species and localities. The prevalence of anti-B. dendrobatidis bacteria was higher at Shenandoah National Park (NP), VA, with 96% (25/26) of salamanders hosting at least one anti-B. dendrobatidis bacterial species compared to 50% (7/14) at Catoctin Mountain Park (MP), MD, and 38% (8/21) at Mt. Rogers National Recreation Area (NRA), VA. At the individual level, salamanders at Shenandoah NP had more anti-B. dendrobatidis bacteria per individual (µ = 3.3) than those at Catoctin MP (µ = 0.8) and at Mt. Rogers NRA (µ = 0.4). All salamanders tested negative for B. dendrobatidis Anti-B. dendrobatidis bacterial species are diverse in central Appalachian Plethodon salamanders, and their distribution varied geographically. The antifungal bacterial species that we identified may play a protective role for these salamanders.IMPORTANCE Amphibians harbor skin bacteria that can kill an amphibian fungal pathogen, Batrachochytrium dendrobatidis Some amphibians die from B. dendrobatidis infection, whereas others do not. The bacteria that can kill B. dendrobatidis, called anti-B. dendrobatidis bacteria, are thought to influence the B. dendrobatidis infection outcome for the amphibian. Yet how anti-B. dendrobatidis bacterial species vary among amphibian species and populations is unknown. We determined the distribution of anti-B. dendrobatidis bacterial species among three salamander species (n = 61) sampled at three localities. We identified 50 unique anti-B. dendrobatidis bacterial species and found that all of the tested salamanders were negative for B. dendrobatidis Five anti-B. dendrobatidis bacterial species were commonly detected, suggesting a stable, functional association with these salamanders. The number of anti-B. dendrobatidis bacteria per individual varied among localities but not among co-occurring salamander species, demonstrating that environment is more influential than host factors in structuring the anti-B. dendrobatidis bacterial community. These anti-B. dendrobatidis bacteria may serve a protective function for their salamander hosts.