Network Analysis to Identify Multi-Omic Correlations in the Lower Airways of Children With Cystic Fibrosis

The leading cause of morbidity and mortality in cystic fibrosis (CF) is progressive lung disease secondary to chronic airway infection and inflammation; however, what drives CF airway infection and inflammation is not well understood. By providing a physiological snapshot of the airway, metabolomics can provide insight into these processes. Linking metabolomic data with microbiome data and phenotypic measures can reveal complex relationships between metabolites, lower airway bacterial communities, and disease outcomes. In this study, we characterize the airway metabolome in bronchoalveolar lavage fluid (BALF) samples from persons with CF (PWCF) and disease control (DC) subjects and use multi-omic network analysis to identify correlations with the airway microbiome. The Biocrates targeted liquid chromatography mass spectrometry (LC-MS) platform was used to measure 409 metabolomic features in BALF obtained during clinically indicated bronchoscopy. Total bacterial load (TBL) was measured using quantitative polymerase chain reaction (qPCR). The Qiagen EZ1 Advanced automated extraction platform was used to extract DNA, and bacterial profiling was performed using 16S sequencing. Differences in metabolomic features across disease groups were assessed univariately using Wilcoxon rank sum tests, and Random forest (RF) was used to identify features that discriminated across the groups. Features were compared to TBL and markers of inflammation, including white blood cell count (WBC) and percent neutrophils. Sparse supervised canonical correlation network analysis (SsCCNet) was used to assess multi-omic correlations. The CF metabolome was characterized by increased amino acids and decreased acylcarnitines. Amino acids and acylcarnitines were also among the features most strongly correlated with inflammation and bacterial burden. RF identified strong metabolomic predictors of CF status, including L-methionine-S-oxide. SsCCNet identified correlations between the metabolome and the microbiome, including correlations between a traditional CF pathogen, Staphylococcus, a group of nontraditional taxa, including Prevotella, and a subnetwork of specific metabolomic markers. In conclusion, our work identified metabolomic characteristics unique to the CF airway and uncovered multi-omic correlations that merit additional study.