Pulmonary microbiome and transcriptome signatures reveal distinct pathobiologic states associated with mortality in two cohorts of pediatric stem cell transplant patients
Authors: Zinter Matt S. , Dvorak Christopher C. , Mayday Madeline Y. , Reyes Gustavo , Simon Miriam R. , Pearce Emma M. , Kim Hanna , Shaw Peter J. , Rowan Courtney M. , Auletta Jeffrey J. , Martin Paul L. , Godder Kamar , Duncan Christine N. , Lalefar Nahal R. , Kreml Erin M. , Hume Janet R. , Abdel-Azim Hisham , Hurley Caitlin , Cuvelier Geoffrey D.E. , Keating Amy K. , Qayed Muna , Killinger James S. , Fitzgerald Julie C. , Hanna Rabi , Mahadeo Kris M. , Quigg Troy C. , Satwani Prakash , Castillo Paul , Gertz Shira J. , Moore Theodore B. , Hanisch Benjamin , Abdel-Mageed Aly , Phelan Rachel , Davis Dereck B. , Hudspeth Michelle P. , Yanik Greg A. , Pulsipher Michael A. , Sulaiman Imran , Segal Leopoldo N. , Versluys Birgitta A. , Lindemans Caroline A. , Boelens Jaap J. , DeRisi Joseph L. 
Published: November 2023

Lung injury is a major determinant of survival after pediatric hematopoietic cell transplantation (HCT). A deeper understanding of the relationship between pulmonary microbes, immunity, and the lung epithelium is needed to improve outcomes. In this multicenter study, we collected 278 bronchoalveolar lavage (BAL) samples from 229 patients treated at 32 children’s hospitals between 2014-2022. Using paired metatranscriptomes and human gene expression data, we identified 4 patient clusters with varying BAL composition. Among those requiring respiratory support prior to sampling, in-hospital mortality varied from 22-60% depending on the cluster (p=0.007). The most common patient subtype, Cluster 1, showed a moderate quantity and high diversity of commensal microbes with robust metabolic activity, low rates of infection, gene expression indicating alveolar macrophage predominance, and low mortality. The second most common cluster showed a very high burden of airway microbes, gene expression enriched for neutrophil signaling, frequent bacterial infections, and moderate mortality. Cluster 3 showed significant depletion of commensal microbes, a loss of biodiversity, gene expression indicative of fibroproliferative pathways, increased viral and fungal pathogens, and high mortality. Finally, Cluster 4 showed profound microbiome depletion with enrichment of Staphylococci and viruses, gene expression driven by lymphocyte activation and cellular injury, and the highest mortality. BAL clusters were modeled with a random forest classifier and reproduced in a geographically distinct validation cohort of 57 patients from The Netherlands, recapitulating similar cluster-based mortality differences (p=0.022). Degree of antibiotic exposure was strongly associated with depletion of BAL microbes and enrichment of fungi. Potential pathogens were parsed from all detected microbes by analyzing each BAL microbe relative to the overall microbiome composition, which yielded increased sensitivity for numerous previously occult pathogens. These findings support personalized interpretation of the pulmonary microenvironment in pediatric HCT, which may facilitate biology-targeted interventions to improve outcomes.