Microbiome function predicts amphibian chytridiomycosis disease dynamics

[Background] The fungal pathogenBatrachochytrium 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. [MediaObject not available: see fulltext.]K.A.B. was funded by a CASE studentship from NERC, NERC Biomolecular Analysis Facility grant (NBAF939) and an E.P. Abraham Junior Research Fellowship from St Hilda’s College, University of Oxford. M.C.F and T.W.J.G. were funded by NERC award NE/E006701/1 and the Biodiversa project RACE: Risk Assessment of Chytridiomycosis to European Amphibian Biodiversity. T.W.J.G was also funded by Research England and NERC NE/S000062/1. D.S.S. and A.L. received funding through the project People, Pollution, and Pathogens financed through the call “Mountains as Sentinels of Change” by the Belmont-Forum (ANR-15-MASC-0001 - P3, DFG-SCHM3059/6-1, NERC-1633948, NSFC-41661144004). D.S.S. holds the AXA Chair for Functional Mountain Ecology funded by the AXA Research Fund through the project GloMEc and M.C.F. is a fellow in the CIFAR ‘Fungal Kingdoms’ Program