- Open Access
Differences in airway microbiome and metabolome of single lung transplant recipients
Respiratory Research volume 21, Article number: 104 (2020)
Recent studies suggest that alterations in lung microbiome are associated with occurrence of chronic lung diseases and transplant rejection. To investigate the host-microbiome interactions, we characterized the airway microbiome and metabolome of the allograft (transplanted lung) and native lung of single lung transplant recipients.
BAL was collected from the allograft and native lungs of SLTs and healthy controls. 16S rRNA microbiome analysis was performed on BAL bacterial pellets and supernatant used for metabolome, cytokines and acetylated proline-glycine-proline (Ac-PGP) measurement by liquid chromatography-high-resolution mass spectrometry.
In our cohort, the allograft airway microbiome was distinct with a significantly higher bacterial burden and relative abundance of genera Acinetobacter & Pseudomonas. Likewise, the expression of the pro-inflammatory cytokine VEGF and the neutrophil chemoattractant matrikine Ac-PGP in the allograft was significantly higher. Airway metabolome distinguished the native lung from the allografts and an increased concentration of sphingosine-like metabolites that negatively correlated with abundance of bacteria from phyla Proteobacteria.
Allograft lungs have a distinct microbiome signature, a higher bacterial biomass and an increased Ac-PGP compared to the native lungs in SLTs compared to the native lungs in SLTs. Airway metabolome distinguishes the allografts from native lungs and is associated with distinct microbial communities, suggesting a functional relationship between the local microbiome and metabolome.
Recent studies have implicated the lung microbiome in the occurrence of chronic lung diseases such as idiopathic pulmonary fibrosis (IPF), chronic obstructive pulmonary disease, cystic fibrosis and chronic lung allograft dysfunction (CLAD) in lung transplant recipients [1,2,3,4]. There is emerging evidence that early alterations in lung microbiome and/or dysbiosis modulates inflammatory mediators leading to pathogenesis and/or progression of chronic lung diseases including CLAD [5,6,7]. Our group recently reported that a shift to a Proteobacteria dominant allograft microbiome was associated with CLAD in lung transplant recipients . Likewise, microbial adaptations and changes in bacterial diversity have been implicated in progression of fibrosis in IPF subjects . However, the mechanisms involved in the microbiome-host interactions leading to chronic lung inflammation are not well understood.
Delineation of the microbiome signatures and taxonomic profiles of bacterial communities in various disease states is an important first step but does not directly provide an insight into bacteriome-allograft-host interaction. Bacterial colonization and/or infection leading to pathology often results in physiological changes in the host, including alterations in metabolic profile [9, 10]. A better understanding of these metabolic shifts associated with a specific infection can improve our understanding of disease pathophysiology through the identification of by-products of host and microbial metabolism, while also providing vital information about the unique metabolites produced with these ever-changing interactions . In particular, the delineation of bacterial metabolism and the host/allograft inflammatory response is critical to better define factors modulating the local microbiome. Metabolites or metabolomic profiles identified in a defined cohort can also serve as potential biomarkers for disease characterization and/or novel therapeutic targets . In the context of a potential clinical application, unique metabolome signatures in urine have been found to distinguish Streptococcus pneumoniae from Staphylococcus aureus lung infection . Similarly, lung metabolome analysis from HIV subjects have shown that pyochelin, a siderophore produced by Pseudomonas aeruginosa, is elevated in HIV-infected individuals compared to HIV-uninfected individuals . Likewise, specific metabolome pathways have been identified in lung transplant recipients with CLAD , though further investigations are needed to delineate its relevance to CLAD pathobiology. Identification and correlation of novel metabolome profiles associated with specific pathological microbiome signatures may help guide personalized treatment of host disease states .
Dysregulation of muco-ciliary clearance and subsequent increases in bacterial burden are well described in advanced lung disease [17, 18]. Moreover, colonization of pathobionts in diseased native lungs of single lung transplant recipients (SLTs) may alter and/or contribute to the microbiome of the allografts [19, 20]. However, to date, the interaction of the metabolome and microbiome in the allografts and native lungs has not been evaluated in patients after lung transplantation. In this study, we utilized the airway microbiome and metabolome signatures of the native lungs and transplanted lungs (allografts) of SLTs as a model system to answer fundamental questions regarding the inherent lung metabolome and its influence on the lung microbiome. We hypothesized that the airway metabolome would correlate with the abundance of distinct microbiome signatures in the native and allografts lungs of SLTs.
Material and methods
Subjects for the study were recruited from the adult lung transplant program at the University of Alabama at Birmingham between September 2014 to July 2016. Six consecutive adult (> 18 years) single lung transplant recipients undergoing bronchoscopy due to a decline in pulmonary function were recruited for this study. Written consent was obtained for sample collection under an institutional review board–approved protocol (IRB No. X120606006, University of Alabama at Birmingham). Subject details are elucidated in Table 1. Bronchoalveolar lavage (BAL) samples from non-transplant volunteers attending the UAB lung health center clinic and University of South Florida were collected. Demographics of the non-transplant volunteers are detailed in supplementary text.
Bronchoalveolar lavage (BAL) fluid from the allograft (A) (i.e. the transplanted lung) and the native lung (N) of each SLT subject were collected using two separate bronchoscopes. Bronchoscopic control samples (C) were collected from each bronchoscope used prior to the procedure (25 ml of sterile saline flushed via the bronchoscope and collected). For standardization, all bronchoscopies were conducted via the oral route and separate scopes used for sampling the A and N sides of each subject. A total of 120 ml of saline (4 aliquots of 30 cc) were instilled for the BAL on each side (right middle lobe and lingula) and remnant BAL fluid was collected from the last aliquot fraction from the right middle lobe and lingula of allograft and native lung. Similarly, non-transplant normal lung volunteers BAL (H) were collected and processed.
Processing of samples
Two aliquots (5 ml each) of N, A and H BAL fluid samples were centrifuged at 1000 rpm for 5 min to separate the eukaryotic cellular fraction. The supernatants were centrifuged again at 15,000 rpm for 10 min to pellet the bacterial component. Similarly, 2 aliquots (5 ml each) of control (C) samples from each bronchoscope used were collected and centrifuged at 15,000 rpm for 10 min to pellet the bacterial component. All bacterial pellets were then stored at − 80 °C. Supernatants from each samples after collection of bacterial pellet were also stored at − 80 °C. DNA was extracted from the bacterial pellets. The two processed aliquots from each BAL sample served as technical replicates. Microbiome sequencing as below was performed from each BAL sample in duplicates. DNA from pellets was extracted using Zymo DNA extraction (CA) kit. Polymerase chain reaction (PCR) amplification of the V4 region of the 16S rRNA gene was performed and products used for microbial DNA sequencing using illumina miseq platform. Data was analyzed used QIIME pipeline and CLC genomics workbench platforms . 16S quantitative PCR to estimate bacterial count was performed using published methodologies [21,22,23]. Supernatant from the above samples were used to measure Ac-PGP levels, multiplex cytokine assays and untargeted metabolome analysis using liquid chromatography-high resolution mass spectrometry using established techniques. Details of our experimental protocols are provided in the supplementary materials.
To compare β-diversities between individual patient allograft and native lung samples, weighted UniFrac distances were calculated between all pairs of samples and then each sample type was plotted separately in 3-dimensional (3D) space by principal coordinate analysis. The plots were then transformed by Procrustes analysis to achieve maximum alignment. Within the 3D plots, blue color represents 1 sample and the red color represents the other sample, and the 2 points from individual subjects are connected by a bar. If both plots are similar, then the relative distance will be small. The overall similarity is summarized by the M2 value, and statistical goodness-of-fit is measured by a Monte Carlo label permutation approach. To identify individual OTUs at the phylum and/or genus levels that were distinctive between the 2 airway compartments in native and allograft, wilcoxon test (non-parametric). Differences between the native, allograft and normal was calculated using non-parametric ANOVA (Kruskal Wallis) with subgroup analysis conducted by Dunn multiple comparison test. Significant differences in community membership identified via constrained ordination were confirmed by using PERMANOVA (permutational multivariate analysis of variance) and plotted in R package vegan via RDA function.
Features in the LC-MS metabolomics data were aligned and peak areas determined using XCMSonline (https://xcmsonline.scripps.edu/) . Statistical analysis of peak areas was carried using the univariate (Volcano plot) and multivariate (sparse Partial Least Squares-Discriminant Analysis, sPLS-DA) programs of MetaboAnalyst (http://www.metaboanalyst.ca) . Metlin (https://metlin.scripps.edu/)  was used to identify the individual m/z features. Mummichog (http://www.mummichog.org, version 1.0.9)  was used for pathway analysis. Fischer’s Exact test was used for comparison and generation of enrichment P values. P-values for all pathways were then modeled as a Gamma distribution and then adjusted for the permutations. The correlations and interaction plot between proportions of a metabolome and a bacterial phylum were generated using the online software XMWAS (https://kuppal.shinyapps.io/xmwas/) . GRAPHPAD PRISM version 7.0 (GraphPad, Inc., San Diego, CA software) was used for statistical analysis and generation of figures, with statically significance achieved at p value = 0.05.
Differences between the allograft and native lung microbiome
We first determined the bacterial count in Allograft (A), Native (N), normal (H) groups using 16 S quantitative PCR. Overall, bacterial 16S DNA levels reported as log 16S copies/ml of BAL were significantly different among the groups (P = 0.02, A 4.15 × 107 vs N 2.0 × 107, P = 0.1, A 4.15 × 107 vs H1.8 × 105, P = 0.04, N vs H, P=NS) and higher in the allograft compared to normal controls (Fig. 1a). Next, to determine whether A, N and H differed in diversity of bacterial communities between samples, we used the Shannon diversity index , a measure of α-diversity within a sample that represents both species richness and evenness (Fig. 1b). A Shannon diversity index of ≥3.5 indicates a highly diverse bacterial community . Although shannon diversity index of the allograft was lower (2.4) compared to was higher in H (5.38) and N (3.8) samples, these were not statistically different (P = 0.06).
We evaluated the differences in the microbiome compositions between A and N of each individual lung transplant recipient. Using Procrustes analysis, we mapped individual samples on a principal component analyses (PCoA) plot using weighted UniFrac distances. The PCoA showed wide separation between A and N microbiomes of each individual subject (Fig. 1c). Monte Carlo label permutation was used to calculate a M2 value. The M2 value varies between 0 to 1 and a higher number suggests greater variation between the samples groups. The M2 value was 0.715 suggesting a wide variation between the A and N samples derived from individual patients. In addition, the pooled analysis using redundancy analysis (RDA) showed a cluster pattern for A versus the N microbiomes. Significant differences in community membership between the constrained ordination plot were confirmed using PERMANOVA, P = 0.05 (Fig. 1d). When comparing the A and N microbiome to normal lung microbiome (H), again a distinct clustering pattern was noted for each group (PERMANOVA H vs A, P = 0.0009, H Vs N, P = 0.003, Fig. 1e and f). Significant microbiome signature differences between the A and N groups are detailed in Table ST1 in the supplemental materials.
To confirm the validity of our findings and rule out contamination from the bronchoscope, we assessed the microbial signatures present in the bronchoscope washes (C) collected prior to the collection of BAL fluid samples. PCoA analysis showed significant weighted UniFrac distances between the C microbiome of the bronchoscopes when compared to their respective N and A microbiome of individual subjects (Figure S1A and B and Tables ST2 and ST3 in Supplementary Material).
Allograft lung has a higher abundance of genera Acinetobacter, Pseudomonas and increased ac-PGP
Emergence of microbiome dominated by a single and/or a group of bacteria, gives rise to microbial dysbiosis in the gut and lung; this is associated with immune dysregulation and has been linked to several chronic inflammatory disease states [31, 32]. The healthy lung microbiome is largely composed of the phyla Firmicutes, Bacteroidetes and Proteobacteria . To understand the local microbiome of the native lung and allograft, we studied the overall community composition structure in samples from the allograft of SLTs and compared them to the corresponding native lung microbiome. In addition, we analyzed their differences with the normal lung microbiome. Overall, the relative abundance of Phyla Proteobacteria (A 77% vs N 52% vs H 39%, p = 0.02) and Firmicutes (A 14% vs N 29% vs H 41%, p = 0.009) were significantly different between the allograft, native and normal lung controls (Fig. 2a, b, c). Previously, an increased Firmicutes/Bacteroidetes (F/B) ratio has been suggested as a marker for dysbiosis in the gut . We calculated the F/B ratio in our cohorts. Additionally, we calculated the Proteobacteria/Firmicutes (P/F) ratio given that several members of the phyla Proteobacteria have been implicated in the pathogenesis of respiratory diseases. In our sample both F/B (A (3.25) vs N (2.8) vs H (2.9)) ratio and P/F (A (5.5) vs N (1.7) vs H (1.4)) were higher in the allografts compared to both the native and healthy lung controls. At the genera level, Acinetobacter and Pseudomonas were significantly higher in the allografts compared to both the native and normal lung controls (Fig. 2d and e). The hospital-based BAL cultures were negative for Pseudomonas, Acinetobacter or any other significant pathogenic bacteria (Table 1).
Next, we measured the levels of Ac-PGP, a collagen breakdown matrikine peptide that is known to upregulate neutrophil chemotaxis and increase lung vascular permeability [35, 36]. We found that Ac-PGP was elevated (p = 0.03) in the allograft as compared to native lungs (Fig. 2f). Additionally, we measured other neutrophilic chemoattractants, IL-8 and LTB4 in the native and allograft samples and these were not found to be statistically different between the two groups (supplementary Figure 2A and B). Thereafter, we measured the levels of pro-inflammatory cytokines (VEGF, TNF -α, IFN-γ) in the BAL. Overall, the levels of all three cytokines were statistically different between A, N and H. Subgroup analysis showed that VEGF was significantly increased in the A compared to H. In contrast, the levels of IFN-γ were reduced in N compared to H (Fig. 2g).
Airway metabolome distinguishes native lungs and allografts
Airway metabolome can provide an insight in to the inflammatory state of the lung . Using an untargeted metabolomics approach and analysis previously described [38, 39], we identified a total of 4886 mass/charge (m/z) features in LC–MS data collected in both the negative ion and positive ion modes. Of these, 676 m/z features in the negative ion mode and 586 m/z features that were significantly different (p < 0.05) between the allografts and native lungs. After reorganizing the downloaded LC–MS data from XCMSonline to remove early (non-bound) metabolite ions and ions collected in the solvent wash period, a volcano plot (Fig. 3a, b) revealed several negatively charged ions and positively charged ions with fold changes > 1.5 and p-values < 0.05. PCoA and PLS-DA analyses were performed to assess chemometric separation among the allograft and native lung BAL fluid samples. A sparse PLS-DA plot showed complete separation of the allograft and native lung BAL samples in both the negative and positive ion mode (Fig. 3c and d).
The molecular identity of the top variable importance in projection (VIP) score metabolites (Figure S3A, B) observed by nanoLC–MS that were significantly changed in the negative and positive ion mode were determined based on their measured accurate mass using the Metlin database (https://metlin.scripps.edu/) or by fragmentation pattern of the ion. The m/z features 332.278, 316.283, 296.257, 314.268, 298.273 were identified as sphingosine like metabolites and found to be significantly increased in the native lungs compared to the allograft (Figure S3A and Table 2).
Metabolite annotation and pathway enrichment analysis were performed using mummichog analysis to understand the role of the identified ions in different metabolic pathways. Mummichog identified 21 pathways associated with the negative ion metabolites that were found in the native and allograft BAL fluid samples (Table 3). Metabolites identified were known to be involved in fatty acid activation/metabolism, beta-oxidation, Krebs cycle and amino acid and nucleotide metabolism.
Sphingosine metabolites negatively correlate with abundance of Proteobacteria
To understand the correlation between the microbiome and metabolome, we performed integration analysis using sparse PLS-DA regression of all m/z features (4885) and the 16 Phyla found in the native lung and allograft samples. We found that overall 864 m/z features were associated with the various bacterial phyla  at a correlation of > 0.4. The phyla Proteobacteria, Firmicutes and Bacteroidetes had positive and negative correlations with several m/z features both in the native and allograft lungs (Table ST4, ST5). Next, we performed integration analysis of the top 25 m/z features by VIP score each in the positive and negative mode and the significantly different microbial genera  between the two groups. The VIP score is a measure of a variable’s importance in the PLS-DA model. The VIP score of a variable is calculated as a weighted sum of the squared correlations between the PLS-DA components and the original variable. It summarizes the contribution a variable makes to the model. Of the top 50 metabolites by VIP score in the positive and negative mode, 35 m/z features correlated with 18 bacterial genera in the allografts and 26 m/z features correlated with 20 bacterial genera in the native lung. In the native lungs, the genera associated with most features at a correlation ≥0.4 were Acinetobacter (22 features) and Pseudomonas (17 features) (Fig. 4a). Whereas in the allografts, the genera associated with the greatest number of features at a correlation ≥0.4 were Acinetobacter (18 features), Pseudomonas (17 features), Bacteroidetes (18 features) and Streptococcus (17 features). The genera Burkholderia, Enterobacter and some of their associated m/z features were in a separate network compared to the others as illustrated in the network plot 4B. In the native lungs, several sphingosine-like metabolites (e.g., m/z 296.257, 298.273, 314.268) had a negative correlation with bacterial genera (Pseudomonas, Acinetobacter) from the phyla Proteobacteria (Table ST6, Fig. 4a). Similarly, amino acid metabolites m/z features (201.1, 186.115) had positive correlation with genus Acinetobacter and Pseudomonas in the allograft lung.
Emerging evidence suggests a role for the lung microbiome in the pathogenesis of chronic inflammatory lung diseases including CLAD [5, 8, 40]. In this study, we utilized the airway microbiome and metabolome signatures of the native lungs and transplanted lungs (allografts) of SLTs as a model system to answer fundamental questions regarding the inherent lung metabolome and its influence on the lung microbiome. We characterized the associations between the airway microbiome and metabolome of the allograft and the native lungs of SLTs. In our cohort, we found that the airway microbiome of the native and allograft lungs were distinct with a significantly higher abundance of genus Pseudomonas and Acinetobacter (phyla Proteobacteria) and elevated levels of VEGF and ac-PGP in the allograft. Furthermore, the native lung metabolome differed from the allograft with a higher abundance of sphingosine and sphingosine-like metabolites and its presence negatively correlated with the abundance of bacteria Pseudomonas and Acinetobacter. These hypothesis generating observations lay the foundation for future studies to evaluate the cause-effect relationship between the airway metabolome and microbiome.
To our knowledge, this is the first report characterizing the microbial heterogeneity of the native and transplanted lung in the same individual. The results from our cohort suggests that the allograft’s airway microbiome is distinct with a significantly increased relative abundance of genera Acinetobacter & Pseudomonas (phyla Proteobacteria) (Figs. 1a, 2d and e). In comparison, the microbiome in native lungs of SLTs and healthy individuals have a greater proportion of the phyla Firmicutes and Bacteroidetes. Alterations in the Firmicutes/Bacteroidetes (F/B) ratio have previously been proposed as a marker for intestinal dysbiosis in various disease states . In the lung, where a higher proportion of Proteobacteria is present than the gut , the Proteobacteria/Firmicutes (P/F) ratio may be a better marker for dysbiosis, given the predominance of Firmicutes in the healthy lung . In our cohort, the allografts had an elevated F/B and P/F ratios compared to the native lungs and normal lung controls suggesting greater dysbiosis. Donor lungs at the time of implantation harbor a different microbiome than the recipient, indicative of greater microbial variation between the allograft and native lung at the time of transplantation. Following transplantation, common thinking is that over time the microbiomes of the allograft and native lung would become similar. However, in our cohort, despite the subjects being more than a year post-lung transplantation, the microbiome remained varied between the allograft and native lung. Likewise, we noted an increased bacterial biomass burden and dysbiosis in the allografts compared to the native and healthy lung controls. These results are concordant with other studies linking elevated bacterial biomass to microbial dysbiosis . Increased bacterial biomass and colonization with pathogenic bacteria such as Acinetobacter and Pseudomonas has been associated with allograft dysfunction in lung transplant recipients .
Cytokine VEGF and tripeptide Ac-PGP were found to be significantly elevated in allograft BAL samples. We and others have shown that Ac-PGP, a matrikine tripeptide, mediates inflammation in acute and chronic lung diseases including CLAD and bacterial infections [43, 44]. Interestingly, elevated VEGF levels has been linked to lung infections with Pseudomonas  as well as allograft rejection including post lung transplant primary graft dysfunction and bronchiolitis obliterans syndrome, a form of chronic lung allograft dysfunction [46,47,48]. Likewise, dysbiosis in the lung and gut has been implicated in heightening inflammatory states in several chronic diseases [31, 32]. The increase in the ratio of the pro-inflammatory bacteria such as Pseudomonas and low stimulatory bacteria such as Prevotella and Streptococcus are known to upregulate the inflammatory gene expression profile . Although we see an association between elevated Proteobacteria the inflammatory markers Ac-PGP and VEGF, additional investigations need to be conducted to establish a causal link between a Proteobacteria-dominant microbiome and pro-inflammation.
To further understand the functional impact of the microbiome on the host, we conducted metabolome analyses. The airway metabolome differentiated the allograft and the native lung of the participants of our study (Fig. 3b and c). Several of these metabolites had significant negative and positive correlations with various bacterial genera (Fig. 4a, b). In the positive ion mode, the top m/z features that differentiated the allograft from the native lung (high VIP scores) were dominated by sphingosine-like molecules that were found to be increased in the native lungs (Table 2 and Figure S3A &B). These sphingosine metabolites had a negative correlation with bacterial genera (Pseudomonas, Acinetobacter) from the phyla Proteobacteria in the native lungs. Sphingolipids are bioactive lipids known to be part of the plasma membrane lipid bilayer in eukaryotic cells  and are cleaved by sphingosine kinases, with several of these molecules having key roles in the regulation of oxidative stress and immune function . Sphingosine improves the host response to Pseudomonas infections and augments neutrophil killing of reactive oxygen species resistant Pseudomonas [51,52,53]. An increased presence of sphingosine-like molecules in the metabolome of native lungs may account for the lower Proteobacteria signature found in them as compared to the allografts. Likewise, some metabolite features had positive correlation with the microbial genera (Pseudomonas and Acinetobacter) in the allograft lungs. Metabolite m/z 201.1, a alanine metabolite and m/z 186.115, a threonine metabolite are involved in maintenance of bacterial cell wall structure and cellular stiffness, promote bacterial proliferation [54, 55] and have regulatory role in T cell activation . Mummichog analysis of the metabolome suggested the presence of metabolites that were increased in the allograft compared to the native lung (Table 3). Pathways in the amino acid metabolism, fatty acid activation and fatty acid beta-oxidation were increased in the allograft compared to the native lung. Methionine and cysteine metabolism pathways are known to regulate oxidative stress through the methionine/glutathinone trans-sulfuration pathway. Disruption of this pathway can lead to increased oxidative injury. Likewise, increased fatty acid activation can result in mobilization of cell membrane derived lipid signaling molecules such as sphingospine-1-kinase and arachidonic acid derived eicosainoides [57, 58]. More studies are needed to further understand if the variations in metabolome are in part influenced by the difference in airway microbiome or vice-versa.
As the primary antigen presenting cells of the lung, resident alveolar macrophages originate in the embryogenesis period and are sparsely replenished by the bone marrow during adult life . Studies in lung transplant recipients have shown that most of the alveolar macrophages even 3.5 years after lung transplant are donor derived in the allografts . An alternate explanation to the difference in the lung microbiome in the native and allograft lungs could be due to the variable antigen presenting/regulatory nature of donor-derived and recipient alveolar macrophage phenotypes in the allografts and native lungs.
Our study has limitations, including sample size and design. We did not find any differences in the microbiome or inflammatory signatures in those diagnosed with CLAD at sampling or at follow-up. However, with a cross-sectional study design and inadequate power, we cannot comment on the dynamic changes in the microbiome and metabolome of native lungs and allografts and their association with chronic lung allograft dysfunction. However, this study is hypothesis generating and provides insightful data to further advance our understanding of microbiome-metabolome interactions in the lung. One of the unique challenges in the field of lung microbiome is presence of low biomass in BAL samples compared to samples from the GI tract [61, 62]. Interpretation of results obtained from a single low biomass sample without appropriate controls can be misleading. To circumvent this issue, we collected control samples from the separate bronchoscopes used during each A and N sample collection. The bronchoscope control microbiome signatures were found to be dissimilar to the A and N airway microbiome (Figure S2A, B). Furthermore, the bacterial 16S DNA levels in the bronchoscope control samples were below the lower limit of PCR quantification suggesting absence of bronchoscope contamination (data not shown). BAL samples from normal lung control subjects, collected and processed using a similar methodology were also analyzed and compared to the allograft/native lung microbiome. We processed technical replicates of each biological sample to evaluate variability in both the relative abundance and bacterial 16S DNA quantification. Although, antimicrobial prophylaxis and immunosuppression regimen can also influence the microbiome , in our subjects these were given systemically, and hence would impact both the native and the allograft similarly. While we found that lung allografts harbored a more Proteobacteria dominant microbiome and higher PGP levels, longitudinal studies are needed to investigate the causal association, role of varied prophylactic antibiotics, immunosuppression on microbial composition in the lung and whether this inflammation translates into allograft rejection. Our healthy controls did not have Ac-PGP measurements. Nevertheless, our previous studies have shown absent or extremely low levels of Ac-PGP in healthy lung airway fluid . It has been reported that BAL samples obtained from separate geographic regions within the same lung can demonstrate highly dissimilar microbial communities [33, 65]. Although, it is possible that the variations in microbiome in the native lung and allograft may be related to expected variation in sampling from different regions of the lung, these spatial variations are known to be less significant compared to variations across individuals . Finally, all metabolites that correlated with bacterial genera were not able to be identified due to the current limitations of untargeted mass spectrometry and the adjusted P values had limitations due to the smaller samples size . Although, we found correlations between sphingosine like molecules with bacterial genera in the native lungs, due to the very low concentration of these sphingosine metabolites in the allograft, these bacterial-sphingosine correlations could not be accurately predicted in the allograft and targeted LC-MS analysis was not performed to validate the metabolite. Nevertheless, these observations are hypothesis generating and warrant further validation in future mechanistic studies.
Allograft lungs have a distinct microbiome signature, and an increased pro-inflammatory milieu compared to the native lungs in SLTs. Airway metabolome distinguishes the allografts from native lungs and is associated with distinct microbial communities, suggesting a functional relationship between the local microbiome and metabolome. Our studies have characterized, for the first time, the native lung and allograft microbiome of SLTs and provide fundamental insight into host-microbiome and metabolome interactions, a potentially important early feature for long-term allograft viability and host-pathogen related lung injury. These hypothesis generating results pave way for future well powered longitudinal studies to elucidate dynamic changes in the lung microbiome and metabolome, determine their associated interactions, and measure their subsequent impact on lung allograft health and function.
Availability of data and materials
The data that support the findings of this study are openly available in UAB CCTS repository at https://data.genome.uab.edu/bmi/microbiome/result2017/M146_analysis/Sharma2nd_analysis/ANALYSIS/report_files/initial_overview.html
Chronic Lung Allograft Dysfunction
Single lung transplant recipients
Annual Data Report of the US Organ Procurement and Transplantation Network. Preface. Am J Transplant. 2014;14(Suppl 1):5–7.
Burton CM, Carlsen J, Mortensen J, Andersen CB, Milman N, Iversen M. Long-term survival after lung transplantation depends on development and severity of bronchiolitis obliterans syndrome. J Heart Lung Transplant. 2007;26(7):681–6.
Huang YJ, Erb-Downward JR, Dickson RP, Curtis JL, Huffnagle GB, Han MK. Understanding the role of the microbiome in chronic obstructive pulmonary disease: principles, challenges, and future directions. Transl Res. 2017;179:71–83.
O'Dwyer DN, Ashley SL, Gurczynski SJ, Xia M, Wilke C, Falkowski NR, et al. Lung microbiota contribute to pulmonary inflammation and disease progression in pulmonary fibrosis. Am J Respir Crit Care Med. 2019;199(9):1127–38.
Bernasconi E, Pattaroni C, Koutsokera A, Pison C, Kessler R, Benden C, et al. Airway microbiota determines innate cell inflammatory or tissue remodeling profiles in lung transplantation. Am J Respir Crit Care Med. 2016;194(10):1252–63.
Charlson ES, Diamond JM, Bittinger K, Fitzgerald AS, Yadav A, Haas AR, et al. Lung-enriched organisms and aberrant bacterial and fungal respiratory microbiota after lung transplant. Am J Respir Crit Care Med. 2012;186(6):536–45.
Litvak Y, Byndloss MX, Tsolis RM, Baumler AJ. Dysbiotic Proteobacteria expansion: a microbial signature of epithelial dysfunction. Curr Opin Microbiol. 2017;39:1–6.
Sharma NS, Wille KM, Athira S, Zhi D, Hough KP, Diaz-Guzman E, et al. Distal airway microbiome is associated with immunoregulatory myeloid cell responses in lung transplant recipients. J Heart Lung Transplant. 2017;S1053-2498(17):31898–3.
Ceccarani C, Foschi C, Parolin C, D'Antuono A, Gaspari V, Consolandi C, et al. Diversity of vaginal microbiome and metabolome during genital infections. Sci Rep. 2019;9(1):14095.
Ding S, Ma Y, Liu G, Yan W, Jiang H, Fang J. Lactobacillus brevis alleviates DSS-induced colitis by reprograming intestinal microbiota and influencing serum Metabolome in murine model. Front Physiol. 2019;10:1152.
Schoeman JC, Harms AC, van Weeghel M, Berger R, Vreeken RJ, Hankemeier T. Development and application of a UHPLC-MS/MS metabolomics based comprehensive systemic and tissue-specific screening method for inflammatory, oxidative and nitrosative stress. Anal Bioanal Chem. 2018;410(10):2551–68.
Stringer KA, McKay RT, Karnovsky A, Quemerais B, Lacy P. Metabolomics and its application to acute lung diseases. Front Immunol. 2016;7:44.
Slupsky CM, Cheypesh A, Chao DV, Fu H, Rankin KN, Marrie TJ, et al. Streptococcus pneumoniae and Staphylococcus aureus pneumonia induce distinct metabolic responses. J Proteome Res. 2009;8(6):3029–36.
Cribbs SK, Park Y, Guidot DM, Martin GS, Brown LA, Lennox J, et al. Metabolomics of bronchoalveolar lavage differentiate healthy HIV-1-infected subjects from controls. AIDS Res Hum Retrovir. 2014;30(6):579–85.
Ciaramelli C, Fumagalli M, Viglio S, Bardoni AM, Piloni D, Meloni F, et al. (1) H NMR to evaluate the Metabolome of Bronchoalveolar lavage fluid (BALf) in bronchiolitis Obliterans syndrome (BOS): toward the development of a new approach for biomarker identification. J Proteome Res. 2017;16(4):1669–82.
Shaffer M, Armstrong AJS, Phelan VV, Reisdorph N, Lozupone CA. Microbiome and metabolome data integration provides insight into health and disease. Transl Res. 2017;189:51–64.
Evans CM, Fingerlin TE, Schwarz MI, Lynch D, Kurche J, Warg L, et al. Idiopathic pulmonary fibrosis: a genetic disease that involves Mucociliary dysfunction of the peripheral airways. Physiol Rev. 2016;96(4):1567–91.
Dransfield MT, Wilhelm AM, Flanagan B, Courville C, Tidwell SL, Raju SV, et al. Acquired cystic fibrosis transmembrane conductance regulator dysfunction in the lower airways in COPD. Chest. 2013;144(2):498–506.
Westney GE, Kesten S, De Hoyos A, Chapparro C, Winton T, Maurer JR. Aspergillus infection in single and double lung transplant recipients. Transplantation. 1996;61(6):915–9.
Frost AE, Keller CA, Noon GP, Short HD, Cagle PT. Outcome of the native lung after single lung transplant. Multiorgan Transplant Group. Chest. 1995;107(4):981–4.
Dickson RP, Erb-Downward JR, Freeman CM, Walker N, Scales BS, Beck JM, et al. Changes in the lung microbiome following lung transplantation include the emergence of two distinct Pseudomonas species with distinct clinical associations. PLoS One. 2014;9(5):e97214.
Simon TD, Van Yserloo B, Nelson K, Gillespie D, Jensen R, McAllister JP 2nd, et al. Use of quantitative 16S rRNA PCR to determine bacterial load does not augment conventional cerebrospinal fluid (CSF) cultures among children undergoing treatment for CSF shunt infection. Diagn Microbiol Infect Dis. 2014;78(2):188–95.
Grice EA, Kong HH, Renaud G, Young AC, Program NCS, Bouffard GG, et al. A diversity profile of the human skin microbiota. Genome Res. 2008;18(7):1043–50.
Forsberg EM, Huan T, Rinehart D, Benton HP, Warth B, Hilmers B, et al. Data processing, multi-omic pathway mapping, and metabolite activity analysis using XCMS online. Nat Protoc. 2018;13(4):633–51.
Chong J, Soufan O, Li C, Caraus I, Li S, Bourque G, et al. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46(W1):W486–W94.
Guijas C, Montenegro-Burke JR, Domingo-Almenara X, Palermo A, Warth B, Hermann G, et al. METLIN: a technology platform for identifying Knowns and unknowns. Anal Chem. 2018;90(5):3156–64.
Li S, Park Y, Duraisingham S, Strobel FH, Khan N, Soltow QA, et al. Predicting network activity from high throughput metabolomics. PLoS Comput Biol. 2013;9(7):e1003123.
Uppal K, Ma C, Go YM, Jones DP, Wren J. xMWAS: a data-driven integration and differential network analysis tool. Bioinformatics. 2018;34(4):701–2.
Kim BR, Shin J, Guevarra R, Lee JH, Kim DW, Seol KH, et al. Deciphering diversity indices for a better understanding of microbial communities. J Microbiol Biotechnol. 2017;27(12):2089–93.
Li K, Bihan M, Yooseph S, Methe BA. Analyses of the microbial diversity across the human microbiome. PLoS One. 2012;7(6):e32118.
O'Dwyer DN, Zhou X, Wilke CA, Xia M, Falkowski NR, Norman KC, et al. Lung Dysbiosis, Inflammation, and Injury in Hematopoietic Cell Transplantation. Am J Respir Crit Care Med. 2018;198(10):1312–21.
Dopkins N, Nagarkatti PS, Nagarkatti M. The role of gut microbiome and associated metabolome in the regulation of neuroinflammation in multiple sclerosis and its implications in attenuating chronic inflammation in other inflammatory and autoimmune disorders. Immunology. 2018;154(2):178–85.
Charlson ES, Bittinger K, Haas AR, Fitzgerald AS, Frank I, Yadav A, et al. Topographical continuity of bacterial populations in the healthy human respiratory tract. Am J Respir Crit Care Med. 2011;184(8):957–63.
Ley RE, Turnbaugh PJ, Klein S, Gordon JI. Microbial ecology: human gut microbes associated with obesity. Nature. 2006;444(7122):1022–3.
Sharma NS, Lal CV, Li JD, Lou XY, Viera L, Abdallah T, et al. The neutrophil Chemoattractant peptide Proline-glycine-Proline is associated with acute respiratory distress syndrome (ARDS). Am J Physiol Lung Cell Mol Physiol. 2018;315(5):L653–61.
Hahn CS, Scott DW, Xu X, Roda MA, Payne GA, Wells JM, et al. The matrikine N-alpha-PGP couples extracellular matrix fragmentation to endothelial permeability. Sci Adv. 2015;1(3):e1500175.
Evans CR, Karnovsky A, Kovach MA, Standiford TJ, Burant CF, Stringer KA. Untargeted LC-MS metabolomics of bronchoalveolar lavage fluid differentiates acute respiratory distress syndrome from health. J Proteome Res. 2014;13(2):640–9.
Prasain JK, Wilson LS, Arabshahi A, Grubbs C, Barnes S. Mass spectrometric evidence for the modification of small molecules in a cobalt-60-irradiated rodent diet. J Mass Spectrom. 2017;52(8):507–16.
Lal CV, Kandasamy J, Dolma K, Ramani M, Kumar R, Wilson L, et al. Early airway microbial metagenomic and metabolomic signatures are associated with development of severe bronchopulmonary dysplasia. Am J Physiol Lung Cell Mol Physiol. 2018;315(5):L810–L5.
Gottlieb J, Mattner F, Weissbrodt H, Dierich M, Fuehner T, Strueber M, et al. Impact of graft colonization with gram-negative bacteria after lung transplantation on the development of bronchiolitis obliterans syndrome in recipients with cystic fibrosis. Respir Med. 2009;103(5):743–9.
Petersen C, Round JL. Defining dysbiosis and its influence on host immunity and disease. Cell Microbiol. 2014;16(7):1024–33.
Vos R, Vanaudenaerde BM, Geudens N, Dupont LJ, Van Raemdonck DE, Verleden GM. Pseudomonal airway colonisation: risk factor for bronchiolitis obliterans syndrome after lung transplantation? Eur Respir J. 2008;31(5):1037–45.
Lal CV, Xu X, Jackson P, Atkinson TP, Faye-Petersen OM, Kandasamy J, et al. Ureaplasma infection-mediated release of matrix metalloproteinase-9 and PGP: a novel mechanism of preterm rupture of membranes and chorioamnionitis. Pediatr Res. 2017;81(1–1):75–9.
Hardison MT, Galin FS, Calderon CE, Djekic UV, Parker SB, Wille KM, et al. The presence of a matrix-derived neutrophil chemoattractant in bronchiolitis obliterans syndrome after lung transplantation. J Immunol. 2009;182(7):4423–31.
Martin C, Thevenot G, Danel S, Chapron J, Tazi A, Macey J, et al. Pseudomonas aeruginosa induces vascular endothelial growth factor synthesis in airway epithelium in vitro and in vivo. Eur Respir J. 2011;38(4):939–46.
Krenn K, Klepetko W, Taghavi S, Lang G, Schneider B, Aharinejad S. Recipient vascular endothelial growth factor serum levels predict primary lung graft dysfunction. Am J Transplant Off J Am Soc Transplant Am Soc Transplant Surg. 2007;7(3):700–6.
Krebs R, Tikkanen JM, Nykanen AI, Wood J, Jeltsch M, Yla-Herttuala S, et al. Dual role of vascular endothelial growth factor in experimental obliterative bronchiolitis. Am J Respir Crit Care Med. 2005;171(12):1421–9.
Xu H, Abuduwufuer A, Lv W, Zhou Z, Yang Y, Zhang C, et al. The role of HIF-1alpha-VEGF pathway in bronchiolitis obliterans after lung transplantation. J Cardiothorac Surg. 2019;14(1):27.
Hannun YA, Obeid LM. Principles of bioactive lipid signalling: lessons from sphingolipids. Nat Rev Mol Cell Biol. 2008;9(2):139–50.
Pchejetski D, Kunduzova O, Dayon A, Calise D, Seguelas MH, Leducq N, et al. Oxidative stress-dependent sphingosine kinase-1 inhibition mediates monoamine oxidase A-associated cardiac cell apoptosis. Circ Res. 2007;100(1):41–9.
LaBauve AE, Wargo MJ. Detection of host-derived sphingosine by Pseudomonas aeruginosa is important for survival in the murine lung. PLoS Pathog. 2014;10(1):e1003889.
Pewzner-Jung Y, Tavakoli Tabazavareh S, Grassme H, Becker KA, Japtok L, Steinmann J, et al. Sphingoid long chain bases prevent lung infection by Pseudomonas aeruginosa. EMBO Mol Med. 2014;6(9):1205–14.
Becker KA, Li X, Seitz A, Steinmann J, Koch A, Schuchman E, et al. Neutrophils kill reactive oxygen species-resistant Pseudomonas aeruginosa by Sphingosine. Cell Physiol Biochem. 2017;43(4):1603–16.
Odermatt PD, Arjes HA, Chang F, Huang KC. Who’s Your DadA? d-Alanine Levels Regulate Bacterial Stiffness. mBio. 2018;9(5):e02127–18.
Shockman GD. Bacterial cell wall synthesis: the effect of threonine depletion. J Biol Chem. 1959;234:2340–2.
Ron-Harel N, Ghergurovich JM, Notarangelo G, LaFleur MW, Tsubosaka Y, Sharpe AH, et al. T cell activation depends on extracellular alanine. Cell Rep. 2019;28(12):3011–21 e4.
Trompette A, Gollwitzer ES, Yadava K, Sichelstiel AK, Sprenger N, Ngom-Bru C, et al. Gut microbiota metabolism of dietary fiber influences allergic airway disease and hematopoiesis. Nat Med. 2014;20(2):159–66.
Koizumi C, Yamada M, Ishizaki K, Ueda T, Sakurai K. Anti-infective control in human bronchiolar epithelial cells by mucin phenotypic changes following uptake of N-acetyl-L-cysteine. Free Radic Res. 2015;49(12):1449–58.
Guilliams M, De Kleer I, Henri S, Post S, Vanhoutte L, De Prijck S, et al. Alveolar macrophages develop from fetal monocytes that differentiate into long-lived cells in the first week of life via GM-CSF. J Exp Med. 2013;210(10):1977–92.
Nayak DK, Zhou F, Xu M, Huang J, Tsuji M, Hachem R, et al. Long-term persistence of donor alveolar macrophages in human lung transplant recipients that influences donor-specific immune responses. Am J Transplant Off J Am Soc Transplant Am Soc Transplant Surg. 2016;16(8):2300–11.
Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87.
Eisenhofer R, Minich JJ, Marotz C, Cooper A, Knight R, Weyrich LS. Contamination in low microbial biomass microbiome studies: issues and recommendations. Trends Microbiol. 2019;27(2):105–17.
Zaza G, Dalla Gassa A, Felis G, Granata S, Torriani S, Lupo A. Impact of maintenance immunosuppressive therapy on the fecal microbiome of renal transplant recipients: comparison between an everolimus- and a standard tacrolimus-based regimen. PLoS One. 2017;12(5):e0178228.
Sharma NS, Lal CV, Li JD, Lou XY, Viera L, Abdallah T, et al. The neutrophil chemoattractant peptide proline-glycine-proline is associated with acute respiratory distress syndrome. Am J Physiol Lung Cell Mol Physiol. 2018;315(5):L653–L61.
Dickson RP, Erb-Downward JR, Freeman CM, McCloskey L, Beck JM, Huffnagle GB, et al. Spatial variation in the healthy human lung microbiome and the adapted island model of lung biogeography. Ann Am Thorac Soc. 2015;12(6):821–30.
Caesar LK, Kellogg JJ, Kvalheim OM, Cech NB. Opportunities and limitations for untargeted mass spectrometry metabolomics to identify biologically active constituents in complex natural product mixtures. J Nat Prod. 2019;82(3):469–84.
We acknowledge Taylor Berryhill for the analyses of acetyl-PGP.
USF New Investigator Award (N.S.S), Angus Cooper Award UAB (N.S.S) for design, data collection, analyses and interpretation. NIH Shared Instrumentation Grant (S10 027822, S. B.) and UAB Health Services Foundation General Endowment Fund (S. B.) for data analyses.
Ethics approval and consent to participate
Ethics approval and consent approval was obtained from the IRBs at University of Alabama at Birmingham and University of South Florida.
Consent for publication
Not applicable. No patient identifier is revealed in the manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Additional file 2.
Additional file 3.
Additional file 4.
Additional file 5.
Additional file 6.
Additional file 7.
Additional file 8.
Additional file 9.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Sharma, N.S., Vestal, G., Wille, K. et al. Differences in airway microbiome and metabolome of single lung transplant recipients. Respir Res 21, 104 (2020). https://doi.org/10.1186/s12931-020-01367-3