Abstract:Abstract Background The fungal pathogen Batrachochytrium dendrobatidis ( Bd ) threatens amphibian biodiversity and ecosystem stability worldwide. Amphibian skin microbial community structure has been linked to the clinical outcome of Bd infections, yet its overall functional importance is poorly understood. Methods Microbiome taxonomic and functional profiles were assessed using high-throughput bacterial 16S rRNA and fungal ITS2 gene sequencing, bacterial shotgun metagenomics and skin mucosal metabolomics. We sampled 56 wild midwife toads ( Alytes obstetricans ) from montane populations exhibiting Bd epizootic or enzootic disease dynamics. In addition, to assess whether disease-specific microbiome profiles were linked to microbe-mediated protection or Bd -induced perturbation, we performed a laboratory Bd challenge experiment whereby 40 young adult A. obstetricans were exposed to Bd or a control sham infection. We measured temporal changes in the microbiome as well as functional profiles of Bd -exposed and control animals at peak infection. Results Microbiome community structure and function differed in wild populations based on infection history and in experimental control versus Bd- exposed animals. Bd exposure in the laboratory resulted in dynamic changes in microbiome community structure and functional differences, with infection clearance in all but one infected animal. Sphingobacterium , Stenotrophomonas and an unclassified Commamonadaceae were associated with wild epizootic dynamics and also had reduced abundance in laboratory Bd -exposed animals that cleared infection, indicating a negative association with Bd resistance. This was further supported by microbe-metabolite integration which identified functionally relevant taxa driving disease outcome, of which Sphingobacterium and Bd were most influential in wild epizootic dynamics. The strong correlation between microbial taxonomic community composition and skin metabolome in the laboratory and field is inconsistent with microbial functional redundancy, indicating that differences in microbial taxonomy drive functional variation. Shotgun metagenomic analyses support these findings, with similar disease-associated patterns in beta diversity. Analysis of differentially abundant bacterial genes and pathways indicated that bacterial environmental sensing and Bd resource competition are likely to be important in driving infection outcomes. Conclusions Bd infection drives altered microbiome taxonomic and functional profiles across laboratory and field environments. Our application of multi-omics analyses in experimental and field settings robustly predicts Bd disease dynamics and identifies novel candidate biomarkers of infection.

The major histocompatibility complex: a double-edged sword in fungal disease susceptibility

Evolutionary analysis of major histocompatibility complex variants in chytrid-resistant and susceptible amphibians

Limited Polymorphism of the Functional MHC Class II B Gene in the Black-Spotted Frog (pelophylax Nigromaculatus) Identified by Locus-Specific Genotyping

Lost but not forgotten: MHC genotypes predict overwinter survival despite depauperate MHC diversity in a declining frog

Pathogenic selection promotes adaptive immune variations against serious bottlenecks in early invasions of bullfrogs

Gene expression differs in susceptible and resistant amphibians exposed to Batrachochytrium dendrobatidis

Major Histocompatibility Complex Genomic Investigation of Endangered Chinese Alligator Provides Insights into the Evolution of Tetrapod Major Histocompatibility Complex and Survival of Critically Bottlenecked Species

Microbiome function predicts amphibian chytridiomycosis disease dynamics

Habitat fragmentation in the Brazilian Atlantic Forest is associated with erosion of frog immunogenetic diversity and increased fungal infections

Meta‐analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Coevolution of a generalist pathogen with many hosts: the case of the amphibian chytrid Batrachochytrium dendrobatidis

Contribution of Multiple Inter-Kingdom Horizontal Gene Transfers to Evolution and Adaptation of Amphibian-Killing Chytrid, Batrachochytrium Dendrobatidis.

What evolutionary processes maintain MHC IIꞵ diversity within and among populations of stickleback?

The MHC Class Ia Genes in Chenfu’s Treefrog (zhangixalus Chenfui) Evolved Via Gene Duplication, Recombination, and Selection

Genomic Transition to Pathogenicity in Chytrid Fungi

Priority effects among amphibian microbes lead to varying levels of Bd infection

Complex history of the amphibian-killing chytrid fungus revealed with genome resequencing data

The commensal skin microbiome of amphibian mountain populations and its association with the pathogen Batrachochytrium dendrobatidis

Evolution by Gene Duplication, Recombination and Selection in MHC Class I Genes of Odorrana Margaretae

Molecular Characterization of the Major Histocompatibility Complex Class Ia Gene in the Black-Spotted Frog, Pelophylax nigromaculata

Presence of low virulence chytrid fungi could protect European amphibians from more deadly strains