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.

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

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

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

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

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

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

Genomic Transition to Pathogenicity in Chytrid Fungi

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

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

Genome variation in the Batrachochytrium pathogens of amphibians

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

Positive Selection of TLR2 and MyD88 Genes Provides Insights Into the Molecular Basis of Immunological Adaptation in Amphibians

High chytrid prevalence and infection intensities in tadpoles of Mixophyes fleayi

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

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

Enzootic frog pathogen Batrachochytrium dendrobatidis in Asian tropics reveals high ITS haplotype diversity and low prevalence

Effects of amphibian phylogeny, climate and human impact on the occurrence of the amphibian-killing chytrid fungus

Environment predicts Batrachochytrium dendrobatidis lineage distribution and zones of recombination in South Africa

Reservoir-host amplification of disease impact in an endangered amphibian

Microbiome function predicts amphibian chytridiomycosis disease dynamics