- Open Access
Correlation of the lung microbiota with metabolic profiles in bronchoalveolar lavage fluid in HIV infection
Microbiomevolume 4, Article number: 3 (2016)
While 16S ribosomal RNA (rRNA) sequencing has been used to characterize the lung’s bacterial microbiota in human immunodeficiency virus (HIV)-infected individuals, taxonomic studies provide limited information on bacterial function and impact on the host. Metabolic profiles can provide functional information on host-microbe interactions in the lungs. We investigated the relationship between the respiratory microbiota and metabolic profiles in the bronchoalveolar lavage fluid of HIV-infected and HIV-uninfected outpatients.
Targeted sequencing of the 16S rRNA gene was used to analyze the bacterial community structure and liquid chromatography-high-resolution mass spectrometry was used to detect features in bronchoalveolar lavage fluid. Global integration of all metabolic features with microbial species was done using sparse partial least squares regression. Thirty-nine HIV-infected subjects and 20 HIV-uninfected controls without acute respiratory symptoms were enrolled. Twelve mass-to-charge ratio (m/z) features from C18 analysis were significantly different between HIV-infected individuals and controls (false discovery rate (FDR) = 0.2); another 79 features were identified by network analysis. Further metabolite analysis demonstrated that four features were significantly overrepresented in the bronchoalveolar lavage (BAL) fluid of HIV-infected individuals compared to HIV-uninfected, including cystine, two complex carbohydrates, and 3,5-dibromo-l-tyrosine. There were 231 m/z features significantly associated with peripheral blood CD4 cell counts identified using sparse partial least squares regression (sPLS) at a variable importance on projection (VIP) threshold of 2. Twenty-five percent of these 91 m/z features were associated with various microbial species. Bacteria from families Caulobacteraceae, Staphylococcaceae, Nocardioidaceae, and genus Streptococcus were associated with the greatest number of features. Glycerophospholipid and lineolate pathways correlated with these bacteria.
In bronchoalveolar lavage fluid, specific metabolic profiles correlated with bacterial organisms known to play a role in the pathogenesis of pneumonia in HIV-infected individuals. These findings suggest that microbial communities and their interactions with the host may have functional metabolic impact in the lung.
Combination anti-retroviral therapy (ART) has reduced the incidence of opportunistic infections in individuals infected with human immunodeficiency virus (HIV), but certain pneumonias and chronic lung diseases remain important clinical problems. Subtle immune deficits may alter the microbial composition in the lung, and these microbes or their products may impact the host.
Culture-independent techniques, such as 16S ribosomal RNA (rRNA) gene sequencing, have been instrumental in identifying bacterial species present in the lung microbiome . Alterations in complex microbial communities have been identified in many diseases compared to healthy individuals ; however, most studies of the lung microbiome have primarily described bacterial communities. These taxonomic analyses of microbial communities provide useful information, but have limited insights into the functional impact of the bacterial microbiome. Metabolomics is the scientific study of small-molecule by-products of host or microbial processes that are part of an integrated biosystems approach; they provide a unique molecular fingerprint for an individual or tissue. The metabolome, which represents a collection of small molecule features found within an organism, is dynamic and can reflect various physiological pathways including bacterial metabolism and the host inflammatory response. New high-throughput technologies enable metabolome-wide association studies (MWAS) that can support integrated models for personalized medicine by correlating specific metabolites with a specific disease.
There are limited data available on the metabolome of HIV-infected individuals, but metabolomic analyses can help elucidate interactions between inflammation and cellular metabolism . A few studies have characterized the profile of metabolites present in oral wash samples  and cerebrospinal fluid . Most recently, we used liquid chromatography-high-resolution mass spectrometry (LC-FTMS) to perform metabolomic analyses of bronchoalveolar lavage (BAL) fluid of otherwise healthy HIV-infected and HIV-uninfected individuals . We demonstrated that overall metabolic profiles in BAL fluid were different between these two groups and showed that pyochelin, a siderophore produced by Pseudomonas aeruginosa, was elevated in the lungs of HIV-infected individuals compared to HIV-uninfected individuals . While this study did not investigate microbial communities in the lung, many of the significant metabolites identified were not associated with any known human metabolites and therefore likely produced by microbial species. Given the potential metabolic products from bacteria to stimulate host metabolism, investigation of the relationship of metabolomics and the microbiota has the potential to provide novel insights into pathogenic processes.
Thus, we investigated the respiratory microbiota and global metabolic profiles in the BAL of HIV-infected and HIV-uninfected individuals to determine if specific metabolites and bacteria genera were associated with HIV infection. This study is the first to evaluate the global functional metabolic interactions between the microbial species in the lung and the metabolome in HIV-infected individuals.
Demographics of HIV-infected and HIV-uninfected participants
Fifty-nine participants were enrolled (39 HIV-infected participants and 20 HIV-uninfected); 47 were used for the training set (32 HIV-infected and15 HIV-uninfected subjects) while 12 were used for the validation (test) set (7 HIV-infected and 5 HIV-uninfected). There were no differences in race, sex, or smoking status between the groups (Table 1). The majority of HIV-infected subjects (87.2 %) reported receiving ART at the time of study, and their median CD4 count was 600 cells/μl (interquartile range (IQR) 404–853). There was only one HIV-infected individual with a CD4 count below 200 cells/μl Additional file 1.
Lung microbiota did not distinguish HIV-infected from HIV-uninfected individuals
Taxonomic assignment of the operational taxonomic units (OTUs) revealed a total of 153 bacterial taxa across all samples. The microbial lung communities in these participants fell into one of two “pulmotypes,” dominated by either the order Actinomycetales or the order Bacteroidales; in the latter, we see a combination of the family Prevotellaceae and the family Streptococcaceae. These types did not sort by HIV status. When using principal coordinate analysis (PCoA) of the weighted UniFrac distance between all OTUs identified, the samples also did not cluster or separate by HIV status (Fig. 1). Unlike principal component analysis (PCA), PCoA can use the weighted UniFrac distance between samples, which incorporates both phylogenetic and compositional differences in communities. Comparison of intra-group beta diversity showed a significant difference in UniFrac distances between the HIV-infected and HIV-uninfected subjects, suggesting a higher degree of heterogeneity among the microbial composition of the lower airways in HIV-infected subjects. Furthermore, there was no difference in biodiversity, as measured by Simpson’s Diversity Index (Table 1). There was also no difference in community composition based on institution where the samples were collected. These findings indicate that there were no significant differences between the lung microbiota of our HIV-infected and HIV-uninfected groups.
Metabolomics of bronchoalveolar lavage fluid distinguish HIV-infected from HIV-uninfected individuals
We compared the metabolomics features found in the BALs of the HIV-infected and HIV-uninfected groups. Of the 5930 mass-to-charge ratio (m/z) features identified by C18 chromatography, 12 were significantly different between the HIV-infected and uninfected individuals at a false discovery rate (FDR) = 0.2 (Appendix 1). Unadjusted raw data for these 12 features were also evaluated to determine statistical significance. A logistic regression model using the top three features (p < 0.0002; FDR 0.05) as predictors gave an area under the curve (AUC) of 0.89 and 0.80, respectively, for the 47 training samples and 12 independent test samples (7 HIV-infected vs. 5 HIV-uninfected) used only for validation (see Additional file 2). The 10-fold cross-validation classification accuracy for the 47 training samples was 90.8 %, and the overall classification accuracy for the 12 test samples was 83.3 % using the logistic regression model. Further, PCoA of 12 m/z features from C18 chromatography shows not only that the lung metabolome distinguishes HIV status but also that there is decreased variability in the metabolome among HIV-infected subjects compared to HIV-uninfected subjects (Fig. 2). Using the PCoA analysis allowed for a direct comparison to the microbiota data. Network analysis of the 12 differentiating m/z features based on Pearson correlation showed an additional 79 m/z features that correlated with these features with a |r| > 0.40 at FDR = 0.2.
Further metabolite analysis demonstrated that four features were significantly overrepresented in the BAL fluid of HIV-infected individuals compared to HIV-uninfected, including cystine, two complex carbohydrates and 3,5-dibromo-l-tyrosine (see Additional file 2). Cystine was measured in the BAL fluid using LC-FTMS, as done previously  and showed increased concentrations among HIV-infected subjects compared to HIV-uninfected subjects; this increase was however not statistically significant (7.3 (IQR 3.2–16.9) μM vs. 4.3 (IQR 3.6–9.7) μM, p = 0.53).
There were 231 m/z features significantly associated with peripheral blood CD4 cell counts using sparse partial least squares regression (sPLS) at a variable importance on projection (VIP) threshold of 2 (Fig. 3). One-way hierarchical cluster analysis (HCA) showed that these features were grouped into 58 clusters based on the CD4 cell count. Metabolite annotation and pathway enrichment analysis using Mummichog mapped these 231 m/z features to a number of inflammatory pathways, including fatty acid activation (Table 2).
Integration of microbiome and metabolome
Global integration using sPLS regression of all 5930 m/z features with the 153 microbial taxa identified showed that 109 microbial genera were associated with all features. Of the 91 m/z features associated with HIV, 23 were also found to be associated with 29 microbial genera at a correlation of 0.30 (Fig. 4a). The bacteria associated at the highest correlation with the greatest number of these 23 features included the families Caulobacteraceae (11 features), Staphylococcaceae (12 features), Nocardioidaceae (12 features), and the genus Streptococcus (7 features). Streptococcus and its associated m/z features were in a separate network compared to the other three bacterial families, as illustrated in the network analysis (Fig. 4b). Using Mummichog, metabolite annotation and pathway enrichment analysis was done using the m/z features that were associated with these three bacteria families—Staphylococcaceae, Caulobacteraceae, and Nocardioidaceae. Staphylococcaceae is shown in Fig. 4b as the representative family since pathways were similar for all three. This analysis shows that the metabolic pathways affected most significantly included lineolate, glycerophospholipid, and fatty acid metabolism. Pathway analysis could not be done for Streptococcus given the limited number of features.
We investigated the correlation between the lung microbiota and global metabolic pathways in the BAL of HIV-infected individuals to determine if specific features were altered in relation to HIV infection and to the microbial composition in the lung. The interactions of microorganisms with the host and the environment could lead to the development of a variety of diseases and contribute to their severity. Metabolomic analyses allow us to measure both metabolic changes in the host and in bacteria. A recent study in cystic fibrosis patients identified significant relationships among P. aeruginosa-isolate growth behavior, forced expiratory volume in one second, and differing metabolomics clusters , suggesting that bacterial species in the lungs are affecting metabolic composition and potentially resulting in adverse clinical outcomes. We previously found that metabolomic profiles in BAL fluid differentiated HIV-infected from HIV-uninfected individuals , but it was unclear if those patterns related to bacterial communities. The current study has combined these two modalities and identified a number of clinically relevant features that are altered with respect to specific lung microbial communities in HIV-infected individuals.
We found that specific metabolic pathways were not only associated with HIV status but specific metabolites were also increased in the BAL fluid of HIV-infected individuals compared to HIV-uninfected individuals, one of the most intriguing being cystine. Cysteine, an essential amino acid, is a rate-limiting component of glutathione synthesis and plays an important role in maintaining the detoxification of free radicals and reactive oxygen species [8, 9]. It is also an integral part of fatty acid synthesis . Oxidative stress can be generated through several mechanisms including the oxidation of extracellular thiols, in particular cysteine to cystine , and HIV infection is known to cause oxidative stress [11–13]. Cysteine supply is limiting for important lymphocyte functions and impaired in HIV [14, 15]. In fact, treatment with N-acetylcysteine has been studied in ART-naïve HIV-infected patients [16, 17] and shown relative increases in CD4+ T cell numbers.
In our study, we observed that oxidized cysteine, or cystine, concentrations were increased in HIV-infected individuals and that CD4 counts were also associated with fatty acid metabolism. Further, fatty acid metabolism in addition to glycerophospholipid and linoleate metabolic pathways also correlated with the lung microbial species in these HIV-infected individuals; these have previously been shown to be altered systemically in HIV-infected individuals . High dietary linoleic acid can lead to excessive prostaglandin-E2 production, which can depress cellular immunity specifically associated with the Th1 subset of T cells, which is also reduced in HIV infection [19, 20]. Alterations in surfactant glycerophospholipid have been seen in AIDS patients with acute Pneumocystis pneumonia . Most recently, using combined metabolomic and qRT-PCR analyses, HIV-1 Tat protein was noted to cause significant alterations in glycerophospholipid metabolism . These data suggest that HIV infection alters these metabolic pathways systemically. We now show that these perturbations, which are associated with pathogenic bacteria, are also present in the lung.
The bacterial organisms associated with the highest number of metabolic features—Caulobacteraceae, Staphylococcaceae, and Nocardioidaceae—are known causative organisms in the pathogenesis of pneumonia in HIV-infected patients. Staphylococcaceae are gram-positive organisms that include the genus Staphylococcus. Together with Streptococcus, these are two of the most frequent organisms seen in the era of ART and are often responsible for the greatest incidence of community-acquired pneumonia in HIV-infected adults and children . Similarly, Nocardioidaceae is a family of bacteria with several genera including Norcardia and Rhodococcus. Norcardia are gram-positive organisms with partial acid-fast branching filaments that can be isolated in the soil worldwide, while Rhodococcus equi is a gram-positive coccobacillus that can also cause pneumonia in HIV-infected individuals . Caulobacteraceae, an alpha-Proteobacteria, is common in aquatic and soil environments. It was recently shown that HIV-infected individuals who use cocaine display intestinal dysbiosis resulting in an abundance of Proteobacteria in their feces compared to HIV-uninfected controls . It is possible that this organism, although mainly found in the gut, may be important for HIV pathogenesis in the lung as well. The Caulobacteraceae were not found in the reagent controls or in the relevant BAL controls.
Based on our results, we postulate that HIV infection alters inflammatory and oxidant metabolic pathways within the lung, which then shift the functional properties of the lung microbiome and/or the host response to the lung microbiome (Fig. 5). The resulting increase in inflammation and oxidative stress may further impair lung host immunity leading to an increased risk of pneumonia or to chronic lung disease associated with HIV. In addition, alterations in redox balance may further influence the function of the microbiome and metabolome leading to a dynamic imbalance within the lung.
This study has several limitations. With respect to the microbiota data, we did not explicitly account for contamination (from the environment or upper airways). Although some investigators have noted that BAL fluid samples are not affected by upper airway contamination , others have not corroborated this finding . However, by limiting our analyses to the 153 genera with >5 % relative abundance in all BAL samples, we were able to separate the BAL samples from the reagent controls. Another limitation is that this study only focuses on bacterial communities while other organisms, such as fungi, may also be important . With respect to metabolomics, this study does not distinguish microbial from host metabolic products; however, given the prominent associations between the lung microbiota and metabolic profiles, it demonstrates potentially important relationships for future study. In addition, authenticated standards are not available for many of the features we identified, and so the 91 m/z features of interest have not all been validated. Metabolome data was also not adjusted for BAL or serum urea, which could have introduced error. Further, although 12 m/z features were significantly different between HIV-infected individuals and controls, these results were marginally significant at a FDR of 20 % and should be interpreted with some caution.
In summary, we have performed the first study to investigate the respiratory microbiome-metabolome interaction in HIV-infected individuals. We found 12 metabolic features in BAL that differentiated between HIV-infected and HIV-uninfected participants, and another 79 features were identified in network analysis. Further metabolite analysis demonstrated that cystine was one of four features that were significantly overrepresented in the BAL fluid of HIV-infected individuals compared to HIV-uninfected and cystine concentrations were increased among HIV-infected subjects compared to HIV-uninfected subjects, albeit not statistically significant. Metabolome associations with CD4 counts also showed alterations in inflammatory pathways. Several HIV-differentiating metabolites, including those involved in inflammatory pathways, correlated with lung microbes known to cause lung infections. Further investigation into other aspects of innate lung immunity and the interaction between these organisms and the host immune response will be critical as this field moves forward. These findings may lead directly to a better understanding of the interactions of microbial communities, host physiology, and lung disease in HIV.
Fifty-nine participants from the Pittsburgh Lung HIV cohort were recruited to undergo bronchoscopy at the University of Pittsburgh, the University of California San Francisco, and the University of California Los Angeles [28, 29]. The lung HIV cohort enrolled a total of 396 HIV-infected and HIV-uninfected individuals from 2008 through 2011 from the Multicenter AIDS Cohort Study (MACS) and the Women’s Interagency HIV Study (WIHS) cohorts. Inclusion criteria for the bronchoscopy study included the absence of acute respiratory symptoms or fever in the previous four weeks and use of systemic antimicrobials (excluding antiretroviral treatments) or immunosuppressive medications within 6 months. We obtained clinical and demographic data by standardized subject interview and medical record review. Written informed consent was obtained from all participants after approval of human subjects’ protection protocols from review boards of the University of Pittsburgh, the University of California San Francisco, the University of California Los Angeles, and Emory University. The study was registered at clinicaltrials.gov (NCT00870857). Microbiota data from a subset of these individuals has been previously reported as part of the NHLBI’s Lung HIV Microbiome Project collaboration .
Subjects were asked to fast and to refrain from smoking for at least 12 h before sample collection. Bronchoscopy with BAL was performed according to standardized procedures designed to minimize oral contamination . BAL was performed by sequentially instilling and then withdrawing 50 ml aliquots of sterile normal saline up to a total of 200 ml into the right middle lobe or lingula.
Sample processing and sequencing
DNA was extracted from samples using the Power Soil Extraction kit (Mo Bio, USA). Sequencing was performed using the Roche 454 FLX Titanium platform with primers for the amplification of hypervariable regions 1 through 3 (V1–3) of the 16S rRNA gene, as previously described . Sequence reads were trimmed and analyzed using the QIIME software package . These were then clustered into OTUs and used to characterize the diversity of the lung microbiota in the BAL. Specifically, to estimate the richness and evenness of the microbial community for each BAL sample, we calculated Simpson’s Diversity Index. To determine the phylogenetic and compositional differences between the microbial communities across the samples, we calculated weighted UNIFRAC distances  between pairs of samples. Intra-group beta diversities were also compared. Taxonomic assignment of the OTUs was performed down to the genus level, and genera with less than 5 % relative abundance in all samples were removed, leaving a total of 153 microbial genera across samples.
Samples were analyzed by LC-FTMS, as previously described . m/z features were collected by a Thermo Q Exactive Hybrid Quadrupole-Orbitrap Mass Spectrometer (Thermo Fisher, San Diego, CA) from m/z 85 to 1275 over 10 min. Peak extraction, noise removal, and quantification of ion intensities were performed by an adaptive processing software package (apLCMS) with xMSanalyzer designed for use with LC-FTMS data . For clarification, data extracted on m/z features, which are not definitively identified metabolites, will be referred to as “features” and not metabolites. See Additional file 2.
Biostatistics and bioinformatics
The major predictor variable was HIV infection status. Two-sample Wilcoxon rank-sum tests were performed to compare continuous demographic characteristics between HIV-infected and HIV-uninfected individuals. Pearson chi-square tests were used to compare categorical characteristics. Statistical analyses were conducted in NCSS statistical software, except as indicated below. All reported p values were two-sided; and p-values of less than 0.05, corrected for multiple hypotheses testing as necessary, were considered statistically significant.
High-resolution mass spectral data were filtered to include only m/z features where at least one group (HIV-infected or HIV-uninfected) had values for 50 % of the samples. A total of 5930 m/z features were selected and used for all statistical analyses. LIMMA package  in R (linear models for microarray data) in Bioconductor was used to identify differentially expressed features at a significance threshold of 0.2 after Benjamini-Hochberg FDR adjustment. Two-way HCA was performed to identify clusters of individuals associated with discriminating clusters of features. Principal component analysis using Pirouette version 4.0 (InfoMetrix) was performed, as previously described . Metabolite annotation and pathway enrichment analysis was done using Mummichog software implementation of an approach that uses the collective power of combining metabolite identification and metabolic pathway/network analysis into one site. The sPLS regression method  implemented in the R package mixOmics was used to perform integration and visualization of microbiome × metabolome associations. See Additional file 2.
Availability of data
The sequence data supporting the results of this study is available in NCBI sequence read archive (SRA) under accession SRP065274.
area under the curve
hierarchical cluster analysis
human immunodeficiency virus
liquid chromatography-high-resolution mass spectrometry
- m/z :
Multicenter AIDS Cohort Study
metabolome-wide association studies
operational taxonomic units
principal component analysis
principal coordinate analysis
sparse partial least squares regression
variable importance on projection
Women’s Interagency HIV Study
Huse SM, Ye Y, Zhou Y, Fodor AA. A core human microbiome as viewed through 16S rRNA sequence clusters. PLoS ONE. 2012;7:e34242. doi:10.1371/journal.pone.0034242.
Morris A, Beck JM, Schloss PD, Campbell TB, Crothers K, Curtis JL, et al. Comparison of the respiratory microbiome in healthy nonsmokers and smokers. Am J Respir Crit Care Med. 2013;187:1067–75. doi:10.1164/rccm.201210-1913OC.
Johnson JM, Strobel FH, Reed M, Pohl J, Jones DP. A rapid LC-FTMS method for the analysis of cysteine, cystine and cysteine/cystine steady-state redox potential in human plasma. Clin Chim Acta. 2008;396:43–8. doi:10.1016/j.cca.2008.06.020.
Ghannoum MA, Mukherjee PK, Jurevic RJ, Retuerto M, Brown RE, Sikaroodi M, et al. Metabolomics reveals differential levels of oral metabolites in HIV-infected patients: toward novel diagnostic targets. OMICS. 2011. doi: 10.1089/omi.2011.0035
Wikoff WR, Pendyala G, Siuzdak G, Fox HS. Metabolomic analysis of the cerebrospinal fluid reveals changes in phospholipase expression in the CNS of SIV-infected macaques. J Clin Invest. 2008;118:2661–9. doi:10.1172/JCI34138.
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 Retroviruses. 2014;30:579–85. doi:10.1089/AID.2013.0198.
Kozlowska J, Rivett DW, Vermeer LS, Carroll MP, Bruce KD, Mason AJ, Rogers GB (2013) A relationship between growth behaviour and cystic fibrosis patient lung function identified in a metabolomic investigation. Metabolomics 9. doi: 10.1007/s11306-013-0538-5.
Jones DP. Redefining oxidative stress. Antioxid Redox Signal. 2006;8:1865–79.
Jones DP. Radical-free biology of oxidative stress. Am J Physiol Cell Physiol. 2008;295:C849–68. doi:10.1152/ajpcell.00283.2008.
Berg JM, Tymoczko JL, Stryer L. Fatty acids are synthesized and degraded by different pathways. In: Biochemistry 5th Edition. New York: W H Freeman; 2002. Available from: http://www.ncbi.nlm.nih.gov/books/NBK22554.
Coaccioli S, Crapa G, Fantera M, Del Giorno R, Lavagna A, Standoli ML, et al. Oxidant/antioxidant status in patients with chronic HIV infection. Clin Ter. 2010;161:55–8.
Cribbs SK, Guidot DM, Martin GS, Lennox J, Brown LA. Anti-retroviral therapy is associated with decreased alveolar glutathione levels even in healthy HIV-infected individuals. PLoS ONE. 2014;9:e88630. doi:10.1371/journal.pone.0088630.
Lassiter C, Fan X, Joshi PC, Jacob BA, Sutliff RL, Jones DP, et al. HIV-1 transgene expression in rats causes oxidant stress and alveolar epithelial barrier dysfunction. AIDS Res Ther. 2009;6:1.
Droge W, Eck HP, Gmunder H, Mihm S. Dysregulation of plasma amino acid levels in HIV-infection and cancer and its relevance for the immune system. Amino Acids. 1991;1:193–8. doi:10.1007/BF00806916.
Droge W, Eck HP, Gmunder H, Mihm S. Modulation of lymphocyte functions and immune responses by cysteine and cysteine derivatives. Am J Med. 1991;91:140S–4.
Droge W, Eck HP, Mihm S. HIV-induced cysteine deficiency and T-cell dysfunction—a rationale for treatment with N-acetylcysteine. Immunol Today. 1992;13:211–4. doi:10.1016/0167-5699(92)90156-2.
Kinscherf R, Fischbach T, Mihm S, Roth S, Hohenhaus-Sievert E, Weiss C, et al. Effect of glutathione depletion and oral N-acetyl-cysteine treatment on CD4+ and CD8+ cells. FASEB J. 1994;8:448–51.
Amara A, Chaugier C, Ragnaud JM, Geffard M. Circulating autoantibodies directed against conjugated fatty acids in sera of HIV-1-infected patients. Clin Exp Immunol. 1994;96:379–83.
Katamura K, Shintaku N, Yamauchi Y, Fukui T, Ohshima Y, Mayumi M, et al. Prostaglandin E2 at priming of naive CD4+ T cells inhibits acquisition of ability to produce IFN-gamma and IL-2, but not IL-4 and IL-5. J Immunol. 1995;155:4604–12.
Khoo SH, Wilkins EG, Fraser IS, Hamour AA, Stanford JL. Lack of skin test reactivity to common mycobacterial antigens in human immunodeficiency virus infected individuals with high CD4 counts. Thorax. 1996;51:932–5.
Hoffman AG, Lawrence MG, Ognibene FP, Suffredini AF, Lipschik GY, Kovacs JA, et al. Reduction of pulmonary surfactant in patients with human immunodeficiency virus infection and Pneumocystis carinii pneumonia. Chest. 1992;102:1730–6.
Liao W, Tan G, Zhu Z, Chen Q, Lou Z, Dong X, et al. Combined metabonomic and quantitative real-time PCR analyses reveal systems metabolic changes in Jurkat T-cells treated with HIV-1 Tat protein. J Proteome Res. 2012;11:5109–23. doi:10.1021/pr300173c.
Cilloniz C, Torres A, Polverino E, Gabarrus A, Amaro R, Moreno E, et al. Community-acquired lung respiratory infections in HIV-infected patients: microbial aetiology and outcome. Eur Respir J. 2014;43:1698–708. doi:10.1183/09031936.00155813.
Torres-Tortosa M, Arrizabalaga J, Villanueva JL, Galvez J, Leyes M, Valencia ME, et al. Prognosis and clinical evaluation of infection caused by Rhodococcus equi in HIV-infected patients: a multicenter study of 67 cases. Chest. 2003;123:1970–6.
Volpe GE, Ward H, Mwamburi M, Dinh D, Bhalchandra S, Wanke C, et al. Associations of cocaine use and HIV infection with the intestinal microbiota, microbial translocation, and inflammation. J Stud Alcohol Drugs. 2014;75:347–57.
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:821–30. doi:10.1513/AnnalsATS.201501-029OC.
Cui L, Lucht L, Tipton L, Rogers MB, Fitch A, Kessinger C et al. Topographical diversity of the respiratory tract mycobiome and alteration in HIV and lung disease. Am J Respir Crit Care Med. 2015. doi:10.1164/rccm.201409-1583OC
Fitzpatrick ME, Gingo MR, Kessinger C, Lucht L, Kleerup E, Greenblatt RM, et al. HIV infection is associated with diffusing capacity impairment in women. J Acquir Immune Defic Syndr. 2013;64:284–8. doi:10.1097/QAI.0b013e3182a9213a.
Crothers K, McGinnis K, Kleerup E, Wongtrakool C, Hoo GS, Kim J, et al. HIV infection is associated with reduced pulmonary diffusing capacity. J Acquir Immune Defic Syndr. 2013;64:271–8. doi:10.1097/QAI.0b013e3182a9215a.
Lozupone C, Cota-Gomez A, Palmer BE, Linderman DJ, Charlson ES, Sodergren E, et al. Widespread colonization of the lung by Tropheryma whipplei in HIV infection. Am J Respir Crit Care Med. 2013;187:1110–7. doi:10.1164/rccm.201211-2145OC.
Iwai S, Fei M, Huang D, Fong S, Subramanian A, Grieco K, et al. Oral and airway microbiota in HIV-infected pneumonia patients. J Clin Microbiol. 2012;50:2995–3002. doi:10.1128/JCM.00278-12.
Schloss PD, Gevers D, Westcott SL. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies. PLoS ONE. 2011;6:e27310. doi:10.1371/journal.pone.0027310.
Lozupone C, Knight R. UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol. 2005;71:8228–35. doi:10.1128/AEM.71.12.8228-8235.2005.
Uppal K, Soltow QA, Strobel FH, Pittard WS, Gernert KM, Yu T, et al. xMSanalyzer: automated pipeline for improved feature detection and downstream analysis of large-scale, non-targeted metabolomics data. BMC Bioinformatics. 2013;14:15. doi:10.1186/1471-2105-14-15.
Smyth G. Limma: linear models for microarray data. In: Gentleman VCR, Dudoit S, Irizarry R, Huber W, editors. Bioinformatics and Computational Biology Solutions using R and Bioconductor. 2005. p. 397–420.
Soltow Q, Strobel F, Mansfield K, Wachtman L, Park Y, Jones DP. High-performance metabolic profiling with dual chromatography-Fourier-transform mass spectrometry (DC-FTMS) for study of the exposome. Metabolomics. 2011. In press.
Le Cao KA, Rossouw D, Robert-Granie C, Besse P. A sparse PLS for variable selection when integrating omics data. Stat Appl Genet Mol Biol. 2008;7:Article. doi:10.2202/1544-6115.1390.
The authors would like to thank Vi Linh (Emory), Stephen Stone and Serena Fong, Alisa Malki, Stephen Stone, and Sofya Tokman (UCSF), Cathy Kessinger, Nicholas Leo, Danielle Camp, M. Patricia George, Lorrie Lucht, Matthew Gingo, Robert Hoffman, Meghan Fitzpatrick, John Ries, and Andrew Clarke (Pittsburgh), and John Dermand and Eric Kleerup (UCLA).
Data in this manuscript were collected by the MACS with centers (Principal Investigators) at the University of California, Los Angeles (Roger Detels), U01-AI35040; University of Pittsburgh (Charles Rinaldo), U01-AI35041; and the Center for Analysis and Management of MACS, Johns Hopkins University Bloomberg School of Public Health (Lisa Jacobson), UM1-AI35043. The MACS is funded primarily by the National Institute of Allergy and Infectious Diseases (NIAID), with additional co-funding from the National Cancer Institute (NCI). Targeted supplemental funding for specific projects was also provided by the National Heart, Lung, and Blood Institute (NHLBI) and the National Institute on Deafness and Communication Disorders (NIDCD). MACS data collection is also supported by UL1-TR000424 (JHU CTSA). Website located at http://www.statepi.jhsph.edu/macs/macs.html. The contents of this publication are solely the responsibility of the authors and do not represent the official views of the National Institutes of Health (NIH). Data in this manuscript were collected by the Women’s Interagency HIV Study (WIHS). The contents of this publication are solely the responsibility of the authors and do not represent the official views of the National Institutes of Health (NIH). WIHS (principal investigators): U01-AI-103408; Connie Wofsy Women’s HIV Study, Northern California (Ruth Greenblatt, Bradley Aouizerat, and Phyllis Tien). The WIHS is funded primarily by the National Institute of Allergy and Infectious Diseases (NIAID), with additional co-funding from the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), the National Cancer Institute (NCI), the National Institute on Drug Abuse (NIDA), and the National Institute on Mental Health (NIMH). WIHS data collection was also supported by UL1-TR000004 (UCSF CTSA).
Sources of funding
NIH HL113451 (D.J.); NIH R01 HL090339, K24 HL123342 (A.M.); U01 HL098962 (A.M., E.G.); K24HL123342; HL087713 (LH); NIAID and NCI UO1-AI-35042, 5-MO1-RR-00722 (GCRC), UL1TR000124 (UCLA CTRC), UO1-AI-35043, UO1-AI-37984, UO1-AI-35039, UO1-AI-35040, UO1-AI-37613, UO1-AI-35041 (Multicenter AIDS Cohort); NIAID and NICHHD UO1-AI-35004, UO1-AI-31834, UO1-AI-34994, UO1-AI-34989, UO1-AI-34993, UO1-AI-42590, and UO1-HD-32632 (WIHS). The corresponding author had full access to all the data and is responsible for the decision to submit this manuscript for publication.
Data in this manuscript were collected by the MACS with centers at Los Angeles (U01-AI35040): University of California Los Angeles, Schools of Public Health and Medicine: Roger Detels (PI), Otoniel Martínez-Maza (Co-PI), Aaron Aronow, Robert Bolan, Elizabeth Breen, Anthony Butch, Beth Jamieson, Eric N. Miller, John Oishi, Harry Vinters, Dorothy Wiley, Mallory Witt, Otto Yang, Stephen Young, Zuo Feng Zhang; Pittsburgh (U01-AI35041): University of Pittsburgh, Graduate School of Public Health: Charles R. Rinaldo (PI), Lawrence A. Kingsley (Co-PI), James T. Becker, Ross D. Cranston, Jeremy J. Martinson, John W. Mellors, Anthony J. Silvestre, Ronald D. Stall; and the Data Coordinating Center (UM1-AI35043): The Johns Hopkins University Bloomberg School of Public Health: Lisa P. Jacobson (PI), Alvaro Munoz (Co-PI), Alison, Abraham, Keri Althoff, Christopher Cox, Jennifer Deal, Gypsyamber D’Souza, Priya Duggal, Janet Schollenberger, Eric C. Seaberg, Sol Su, Pamela Surkan. The MACS is funded primarily by the National Institute of Allergy and Infectious Diseases (NIAID), with additional co-funding from the National Cancer Institute (NCI). MACS data collection is also supported by UL1-TR000424 (JHU CTSA). Website located at http://www.statepi.jhsph.edu/macs/macs.html.
Data in this manuscript were collected by the Women’s Interagency HIV Study (WIHS) Collaborative Study Group with centers (principal investigators) at The Connie Wofsy Study Consortium of Northern California (Ruth Greenblatt). The WIHS is funded by the National Institute of Allergy and Infectious Diseases (UO1-AI-35004, UO1-AI-31834, UO1-AI-34994, UO1-AI-34989, UO1-AI-34993, and UO1-AI-42590) and by the Eunice Kennedy Shriver National Institute of Child Health and Human Development (UO1-HD-32632). The study is co-funded by the National Cancer Institute, the National Institute on Drug Abuse, and the National Institute on Deafness and Other Communication Disorders. Funding is also provided by the National Center for Research Resources (UCSF-CTSI Grant Number UL1 RR024131).
The authors declare that they have no competing interests.
The following authors contributed to the study design, literature search, data collection, data analysis, data interpretation, and manuscript preparation: SKC, KU, SL, DPJ, EG, and AM. The following authors contributed to data collection, data analysis, data interpretation, and manuscript preparation: LH, LT, AF, RMG, LK, and DMG. All authors read and approved the final manuscript.
Elodie Ghedin and Alison Morris share senior authorship.
Bal public dataset. (XLSX 12 kb)
Figure S1. Receiver operating characteristic (ROC) curve using the top discriminatory features between HIV-infected (n = 32) and HIV-uninfected (n = 15) individuals. The top three features (p < 0.0002; FDR 0.05) were used as predictors in a logistic regression model that gave an area under the curve (AUC) = 0.89 and a 10-fold cross-validation classification accuracy of 90.8 % for the 47 training samples and AUC of 0.8 and classification accuracy of 83.3 % for the 12 test set samples. Dotted line training set (n = 47). Solid line validation (test) set (n = 12). Figure S2. Metabolomics of bronchoalveolar (BAL) fluid of HIV-infected individuals compared to HIV-uninfected. Further metabolite analysis demonstrates that four features were significantly over-represented in the bronchoalveolar lavage (BAL) fluid of HIV-infected individuals compared to HIV-uninfected, including cystine, two complex carbohydrates and 3,5-dibromo-l-tyrosine. (PDF 106 kb)