Secondary Logo

Journal Logo

BASIC SCIENCE

Effects of HIV viremia on the gastrointestinal microbiome of young MSM

Cook, Ryan R.a; Fulcher, Jennifer A.b,c; Tobin, Nicole H.d; Li, Fand; Lee, Davidd; Javanbakht, Marjana; Brookmeyer, Rone; Shoptaw, Stevef,g; Bolan, Roberth,i; Aldrovandi, Grace M.d,*; Gorbach, Pamina M.a,*

Author Information
doi: 10.1097/QAD.0000000000002132

Abstract

Introduction

The trillions of bacteria, viruses and fungi inhabiting the human gastrointestinal tract have a profound impact on our health and the development of disease. Disruption in the homeostasis of these microbes, a state of ‘dysbiosis,’ has been associated with a broad range of illnesses, including localized gastrointestinal conditions, neurocognitive disorders, cancer, autoimmune disorders and cardiovascular disease [1]. There is tremendous variability in the diversity and composition of the microbiome, even between healthy individuals [2], and the effects of different exposures, behaviours and personal characteristics on the composition and function of the microbiome are incompletely understood.

Chronic inflammation is a hallmark of HIV infection and continues despite suppressive antiretroviral therapy (ART). The gastrointestinal tract is a primary site for HIV replication resulting in a significant loss of CD4+ T cells vital to a healthy mucosal immune system. Depletion of regulatory immune cells and pathways leads to decreased epithelial barrier function allowing translocation of microbes and microbial products, which contributes to the chronic inflammatory response [3,4]. HIV replication may also result in a state of dysbiosis [5–7], which has been correlated with increases in markers of disease progression, microbial translocation and immune activation [5,7–9]. It has been hypothesized that HIV-associated immune dysfunction induces this dysbiosis, and dysbiosis causes further dysfunction [5], thereby driving persistent systemic inflammation in HIV-infected individuals [10].

Many studies have found that overall microbial diversity is reduced in HIV-infected individuals [11–13], and HIV has been associated with a shift from commensals such as Bacteroides to pro-inflammatory taxa such as Prevotella[11,14,15]. However, results have been inconsistent, with some studies showing the opposite or no effects of HIV on microbial diversity [16–18] and a little effect on composition [19,20]. Certainly, differences in sample collection, sequencing, postprocessing and analytic methodology may be responsible for much of the heterogeneity in results [21]. Although some research has focused on the role of sampling variability and type-I error on irreproducibility of findings [22,23], comparatively little attention has been paid to sources of systematic error such as incomparable study cohorts or confounders such as sexual behaviour, substance use, diet, race/ethnicity and age. Such confounders are highly prevalent in observational human studies of the microbiome and may have a larger effect than HIV itself [19]. Due to these limitations, the specific effects of HIV on the microbiome warrant further investigation.

Increased relative abundance of pro-inflammatory bacteria has been correlated with increased viremia, whereas the opposite has been found for potentially beneficial bacteria [9,13,14]. Given these findings and the effects of viral replication on mucosal immunity, it stands to reason that the level of viremia may be a significant determinant of HIV-associated dysbiosis. However, the effect of viremia has not been thoroughly explored. Numerous studies comparing cohorts of HIV-infected individuals that are either on ART or ART-naive have shown that ART does not result in full ‘reconstitution’ of the microbiome, even if the virus is suppressed [13,24,25]. Additional studies focused on elite controllers showed that the microbial composition among individuals with controlled viremia is more similar to HIV-uninfected than viremic individuals [26,27]. However, there are likely to be important biological differences between elite controllers and other HIV-infected individuals that may limit the generalizability of these findings. In order to accurately characterize the effects of HIV on the microbiome, a more detailed examination of the effects of viremia is needed.

To this end, we compared intestinal microbial composition between HIV-uninfected, HIV-infected with undetectable viremia (HIV RNA < 20 copies/ml), HIV-infected with suppressed viremia (HIV RNA ≥20–200 copies/ml) and HIV-infected viremic individuals (HIV RNA > 200 copies/ml). We utilized data from a cohort composed entirely of MSM and employed inverse probability of treatment weighting (IPTW) to control for a robust set of clinical and behavioural confounders. We hypothesized that alterations to the microbiome would be present in all HIV-infected subgroups as compared to HIV-uninfected controls, and the severity of dysbiosis would increase with increasing viremia.

Materials and methods

Study population

Specimens for this study were obtained from an ongoing prospective cohort (The mSTUDY, NIDA U01 DA036267). The mSTUDY was approved by a UCLA Institutional Review Board (IRB) and all individuals provided written informed consent at study entry. Participants are recruited from community clinics in Los Angeles and complete biannual assessments, including a comprehensive physical examination and medical history, urine toxicology panel, clinical laboratory tests including plasma HIV RNA, specimen collection and detailed behavioural questionnaire. Data presented in this manuscript were collected from baseline study visits completed between August 2014 and July 2017. Additional details on sample selection and HIV RNA quantification are provided in the supplemental content, http://links.lww.com/QAD/B423.

Specimen collection and DNA preparation

The majority (76%) of rectal swabs (FLOQSwabs; Copan Diagnostics, Murrieta, California, USA) were collected via anoscopy under direct mucosal visualization and without preparatory enema at approximately 8 cm from the anal verge. Due to an mSTUDY protocol change, others (24%) were participant self-collected at approximately 4–5 cm from the anal verge. Collection method was taken into account in the analysis (see Tables 1 and S1, http://links.lww.com/QAD/B423). Swabs were immediately frozen neat at -80°C until processing in bulk. For DNA processing, the samples were transferred to Lysing Matrix E tubes (MP Biomedicals, Burlingame, California, USA) containing RLT lysis buffer (Qiagen, Hilden, Germany) and bead-beated on a TissueLyser (Qiagen). DNA was then extracted using the AllPrep DNA/RNA/Protein kit (Qiagen) per manufacturer's protocol.

Table 1
Table 1:
Participant characteristics, N = 383 MSM in Los Angeles, California.

16S rRNA gene sequencing and data processing

Microbiome profiling was performed by sequencing of the V4 region of the 16S rRNA gene as previously described [28,29]. Briefly, the V4 region was amplified in triplicate reactions using Golay-barcoded primers 515F/806R. Negative controls from the DNA extraction and PCR steps, as well as independent aliquots of a bacterial mock community [30] were processed alongside the samples to identify contaminant sequences and ensure data reproducibility. PCR products were then pooled and sequenced on the Illumina MiSeq platform using 2 x 150 bp v2 chemistry. The sequences were demultiplexed with Golay error correction using QIIME v1.9.1 [31], and Divisive Amplicon Denoising Algorithm (DADA2) version 1.8 was used for error correction, exact sequence inference, read merging, and chimera removal [32]. Following contaminant removal (see supplemental content, http://links.lww.com/QAD/B423), the amplicon sequence variant (ASV) table comprised 19 955 039 total merged read pairs (mean per sample = 52 375; range 10 906–124 889). Taxonomic assignment was performed using RDP trainset 16 [33]. Rarefaction was performed at a depth of 10 906 reads for alpha diversity analyses. For all other analyses, estimates of relative library sizes ('size factors’) were obtained by calculating geometric means of pairwise read count ratios [34].

Behavioural and clinical covariates

Demographic and behavioural covariates included in the analyses were age, race/ethnicity, employment status, country of origin, a dichotomous variable for homelessness in past six months, number of receptive anal intercourse (RAI) acts in past month, a dichotomous variable for RAI within the past 7 days, frequency of methamphetamine, marijuana and cocaine use in the past 6 months, tobacco smoking and binge drinking. All demographic and behavioural data were self-reported by participants using a computer-aided self-interview (CASI); measures are described in the supplemental content, http://links.lww.com/QAD/B423. Dichotomous variables for obesity (defined as BMI > 30 or waist circumference > 40 inches), and antibiotic use in the past month were also included in the analyses; these data were collected by clinical staff.

Statistical analyses

The primary analyses were unadjusted and inverse probability of treatment-weighted comparisons of microbiome diversity and composition between HIV-negative, HIV-positive undetectable, HIV-positive suppressed and HIV-positive viremic participants. The R package ‘phyloseq’ (R Foundation, Vienna, Austria) was used to calculate alpha diversity statistics, distance matrices and for ordination. Differences in alpha diversity between groups were examined with Kruskal–Wallis tests followed by comparisons of median values using quantile regression (R package ‘quantreg’). Permutational multivariate analysis of variance (PERMANOVA) was used to test for overall differences in microbial composition between HIV groups (R package ‘vegan’). Zero-inflated negative binomial (ZINB) models were fit in order to test for differential abundance in bacterial genera between groups with multinomial least absolute shrinkage and selection operator (LASSO) models employed as a confirmatory analysis (R packages ‘pscl’ and ‘glmnet’). ZINB and LASSO model selection and analytic procedures are described in the supplemental content, http://links.lww.com/QAD/B423.

IPTW [35] is a method of confounder control wherein the study sample is reweighted in order to create a ‘pseudo-population’ in which treatment/exposure, here referring to the four HIV viremia groups, is independent of confounding variables (see supplemental content, http://links.lww.com/QAD/B423). We used IPTW to control for all variables described in the Behavioral and Clinical Covariates section as well as sample collection method. Weights were estimated using generalized boosted models (R package ‘twang’), and balance between groups was assessed by computing standardized mean differences for each covariate in the weighted sample (R package ‘tableone’). Table 1 and Table S1, http://links.lww.com/QAD/B423 provide information on covariate balance before and after weighting. Robust variance estimates for inference tests in weighted ZINB analyses were obtained via the sandwich estimator (R package ‘sandwich’). Additional detail about the IPTW estimation and modelling procedures is provided in the supplemental content, http://links.lww.com/QAD/B423. Missing covariate data were imputed using the Chained Equations method [36] (R package ‘mice’); the proportion of missing data for each covariate is shown in Table S1, http://links.lww.com/QAD/B423. In order to account for multiple testing, alpha diversity and ZINB p values were corrected with the Benjamini–Hochberg false discovery rate (FDR) method [37]; FDR-adjusted P values are labelled as q values. We utilized a threshold of two-sided P or q less than 0.1 for significance testing; accordingly, we also display false coverage rate (FCR)-adjusted 90% confidence intervals [38] where relevant. All statistical analyses were completed using R v.3.4.3.

Results

Sample characteristics

N = 383 participants were included in this study; 200 were HIV-negative, 66 were HIV-positive undetectable (HIV RNA < 20 copies/ml), 72 were HIV-positive suppressed (HIV RNA ≥20–200 copies/ml) and 45 were HIV-positive viremic (HIV RNA > 200 copies/ml). All participants were MSM with an average age 31 (standard deviation = 7). Most were Hispanic (49%) or non-Hispanic black (39%). Table 1 provides further detail on participant characteristics. Generally, HIV-infected participants, especially those with higher levels of viremia, were more likely to be older, unemployed, recently homeless and to report methamphetamine use and frequent binge drinking than their HIV-uninfected peers. Among the HIV-infected participants, the mean number of years since diagnosis was 7.5 (sd = 5.7), mean log10 viral load was 2.0 (sd = 1.2) and mean CD4+ cell count was 625 cells/μl (sd = 287). No participants were ART-naive and 90% of participants reported current ART.

Effects of HIV viremia on overall microbial composition

The relative composition of each individual's microbiome is displayed in Figures 1a and S1, http://links.lww.com/QAD/B423, and average composition within each HIV viremia category in Fig. 1b. Prevotella is the most highly represented bacterial genus among most participants, with increasing relative amounts of Bacteroides, Bifidobacteria and Fusobacteria in those towards the right side of the axis (Fig. 1a). Higher levels of Alloprevotella and Porphyromonas are apparent in the HIV-positive viremic group, and lower levels of Bacteroides are apparent in all HIV-positive groups relative to HIV-negative controls (Fig. 1b).

Fig. 1
Fig. 1:
Rectal microbial composition of study participants, N = 383.(a) Columns represent the relative composition of each individual's microbiome at the genus level. Individuals are ordered by the first principal coordinate of a Bray–Curtis pairwise distance matrix. Genera representing less than 1% of the composition on average across samples were combined into ‘Other’. (b) Average microbial composition within each HIV viremia category. Bacterial genera representing less than 1% of the overall relative composition or present in less than 10% of the samples were grouped into ‘Other’.

To quantitatively examine the influence of HIV viremia on differences in microbial composition between-individuals, we used PERMANOVA with Bray–Curtis distance. HIV viremia explained a significant amount of variability in microbial composition in both unadjusted (R2 = 0.014, P = 0.001) and covariate-adjusted analyses (R2 = 0.011, P = 0.017) (Table S2, http://links.lww.com/QAD/B423). Figure 2a displays ordination of the samples by principal coordinates analysis (Bray–Curtis distance), where HIV-negative and HIV-positive viremic groups are distinct, while HIV-positive undetectable and HIV-positive suppressed are more similar.

Fig. 2
Fig. 2:
Associations between HIV viremia and overall microbial composition.(a) Ordination of the Bray–Curtis distance between samples using principal coordinates analysis. PCoA, principal coordinate axis. Ellipses are 95% confidence regions for each group assuming points follow a multivariate t distribution. (b) Boxplots of richness metrics. Boxes represent the lower, median and upper quartile of the data and whiskers are 1.5*interquartile range.

Comparisons of alpha diversity suggest a tendency for HIV-positive individuals to have higher diversity in metrics that do not account for evenness (observed count and Chao1 statistic) (Fig. 2b). Kruskal–Wallis analyses revealed significant differences in observed and Chao1 richness by HIV group (Table S3, http://links.lww.com/QAD/B423). Quantile regression was further used to investigate these differences and revealed higher median observed and Chao1 values for HIV-positive suppressed versus HIV-negative individuals (q = 0.022 in IPTW-adjusted analyses). No other significant differences were found in any group. Shannon and Simpson indices did not vary greatly between groups.

Differences in specific bacterial taxa associated with HIV viremia

Zero-inflated binomial (ZINB) models were utilized to identify bacterial genera that were differentially abundant among the HIV viremia groups. HIV-positive undetectable participants showed significant enrichment in Finegoldia and Streptococcus and depletion in Bacteroides, Brachyspira, Campylobacter, Helicobacter, Parasutterella and Turicibacter when compared with HIV-uninfected. After IPTW adjustment for behavioural and clinical confounders, depletions in Campylobacter, Parasutterella and Brachyspira remained significant (Fig. 3a; Table S4, http://links.lww.com/QAD/B423). HIV-positive suppressed participants (HIV RNA ≥ 20–200 copies/ml ) had increased Pseudoflavonifractor and decreased Bacteroides, Barnsiella, Brachyspira, Campylobacter, Escherichia/Shigella, Flaonifractor, Helicobacter, Oxalobacter, Parabacteroides, Turicibacter and Victivallis relative to HIV-negative individuals. Following IPTW adjustment, depletions in Barnesiella, Helicobacter and Brachyspira remained significant (Fig. 3b; Table S4, http://links.lww.com/QAD/B423).

Fig. 3
Fig. 3:
Comparisons of individual bacterial genera between HIV viremia categories.Forest plots of results of zero-inflated negative binomial models comparing genus-level bacterial counts between HIV-negative and (a) HIV-positive undetectable (HIV RNA < 20 copies/ml), (b) HIV-positive suppressed (HIV RNA ≥20–200 copies/ml) and (c) HIV-positive viremic (HIV RNA > 200 copies/ml) participants. Inverse probability of treatment-weighted effect sizes and 90% false coverage rate-adjusted confidence intervals (truncated at -4, 4) are plotted, with statistical significance (q < 0.1) indicated in colour. Effect sizes are log ratios of normalized genera counts.

HIV-positive viremic men (HIV RNA > 200 copies/ml) had the most distinct microbial signature relative to HIV-negative, showing significant enrichment of Corynebacterium, Dietzia, Finegolda, Murdochiella, Negativicoccus, Peptoniphilus, Porphyromonas and Prevotella as well as depletion of Arcanobacterium, Brachyspira, Bacteroides, Campylobacter, Faecalibacterium, Helicobacter and Succinivibrio. With IPTW adjustment, enrichment of Murdochiella, Peptoniphilus, Porphyromonas and Prevotella and depletion of Arcanobacterium, Bacteroides, Brachyspira, Faecalibacterium and Helicobacter were significant (Figure 3C; Table S4, http://links.lww.com/QAD/B423). For some bacteria, including Faecalibacterium, Peptoniphilus, Porphyromonas, Prevotella and Streptococcus, effect size (i.e. degree of enrichment or depletion) increased with increasing viremia (Figure S2, http://links.lww.com/QAD/B423).

An IPTW-adjusted multinomial least absolute shrinkage and selection operator (LASSO) model was also used as an additional method of feature selection to compare with ZINB findings. Among HIV-positive participants with undetectable or suppressed viremia, no genera were significant in both the adjusted ZINB and LASSO models. However, among the HIV-positive viremic group, differences in Bacteroides, Peptoniphilus, Porphyromonas and Prevotella were consistent across analytic strategies (Fig. 4).

Fig. 4
Fig. 4:
Summary of zero-inflated negative binomial (ZINB) and least absolute shrinkage and selection operator (LASSO) model results.Enriched taxa are those with positive effect sizes (relative to HIV-negative), depleted are those with negative effect sizes. Genera with no effect in either analysis are not shown. UNW, unadjusted; WT, IPTW adjusted.

It was also of interest to determine whether HIV-positive participants with low levels of viremia (HIV RNA ≥ 20–200 copies/ml) had distinct microbial signatures from those who were HIV-positive but undetectable (HIV RNA < 20 copies/ml). One genus, Sneathea, was significantly different between these groups in adjusted ZINB analyses (q < 0.1). The LASSO model identified depletions in Gemmiger in HIV-positive suppressed as compared to undetectable participants (Figure S3, http://links.lww.com/QAD/B423).

Discussion

In this study examining the effects of HIV on the rectal microbiome in a cohort of 383 young, mostly minority MSM, we found important differences in microbial composition between HIV-uninfected and HIV-infected men, which varied depending on level of viremia. HIV viremia category accounted for about 1% of the variability in microbiome composition, an effect size that is consistent with previous studies [19]. As hypothesized, microbiome perturbations were most evident among HIV-positive viremic men, and least evident in HIV-positive men with undetectable viremia. Importantly, we utilized IPTW to account for multiple confounding factors in our analyses, which decreased the likelihood that the results we report are attributable to clinical or behavioural covariates affecting the microbiome such as sexual behaviour, substance use or obesity.

High diversity is generally associated with a healthy rectal microbiome [2], and reduced richness and diversity has been reported in studies comparing HIV-infected and uninfected persons [11,12,25,39–41]. Still, other studies report no differences in diversity associated with HIV-infection [16,17], while others have suggested that differences in diversity may be related to sampling location [14], or HIV-treatment status [24,40]. We found few significant differences in diversity metrics in our study, and findings did not follow a clear dose–response pattern with level of viremia. As we were able to adjust for multiple confounders in our analyses, we can be reasonably confident that previously reported determinants of diversity such as sexual behaviour [19] and substance use [42] had limited influence on our findings. Our results suggest that once these confounding factors are taken into account, bacterial diversity and richness may not be substantially impacted by HIV infection itself.

One of the more consistent findings across studies of HIV and the microbiome has been enrichment in Prevotella and depletion in Bacteroides among both untreated and treated HIV-infected individuals [9,11,14]. Prevotella species are considered pro-inflammatory [43,44], while Bacteroides species have been shown to induce regulatory T-cell differentiation and IL-10 production [45,46]. Previous work has suggested that observed alterations to the Prevotella/Bacteroides ratio may have been due to sexual behaviour rather than HIV [19]; however, others have shown decreased Bacteroides among HIV-infected MSM who were matched with MSM controls [9]. Our study examined exclusively MSM and controlled for recency and frequency of RAI in our analyses; therefore, our study provides additional evidence that HIV may directly alter the Prevotella/Bacteroides ratio independent of sexual behaviour. Although we found decreased Bacteroides in all HIV-infected individuals, the effect was similar between undetectable and suppressed participants and only statistically significant after adjustment for confounding in the viremic group. In addition, we found increasing relative amounts of Prevotella with increasing levels of viremia, which were only significant in the viremic group. Our findings are consistent with previous research showing that Prevotella may normalize with ART [13], whereas depletions in Bacteroides persist even with therapy [24].

Of the 78 genera tested, ZINB and LASSO models identified Porphyromonas as the genus with the largest difference between HIV-positive viremic and HIV-negative individuals. Porphyromonas is a well known modifier of inflammatory cytokines [47]. In fact, Porphyromonas gingivalis has been identified as a potential cause of systemic inflammation and metabolic disorders associated with periodontal disease [48] and implicated in inflammatory processes leading to the development of atherosclerosis [49]. Furthermore, administration of P. gingivalis to mice was shown to induce gastrointestinal dysbiosis and contribute to intestinal permeability [50]. The association between increasing levels of HIV viremia and Porphyromonas may therefore represent an important mechanism behind HIV-associated chronic inflammation deserving of further study.

Of particular interest in our study is the examination of low-level viremia individuals who are not undetectable (HIV RNA ≥20–200 copies/ml). It is notable that this group, although distinct from HIV-uninfected, was very similar to the undetectable viremia (HIV RNA < 20 copies/ml) group in ordination analysis. Only a single genus, Sneatha, was statistically significantly different between these two groups when directly compared. Clinically, persistent low-level viremia may increase risk of subsequent virologic failure [51], but a recent large study showed no difference in progression to AIDS or incidence of non-AIDS events in persons with low-level viremia compared with undetectable [52]. Our analysis suggests that microbial composition is similar between those with low-level and undetectable viremia, but remains distinguishable from HIV-uninfected individuals. Although those with low-level viremia can still have microbial translocation and inflammatory biomarkers [53], the overall decreased dysbiosis in low-level viremia may correspond to reduced chronic inflammation, which lessens clinical progression.

Our results should be interpreted with consideration of the following limitations. First, we did not have diet information on our cohort. We included race/ethnicity, country of origin, employment status, obesity and homelessness (all of which may impact dietary intake) in our covariate adjustment set to mitigate this limitation to the best of our ability. Second, the IPTW procedure will only achieve perfect balance between exposure groups in nonparametric settings with large sample sizes relative to the number of relevant confounders, and there is potential for residual confounding even by variables we controlled for in our study. However, we note that many of the most significant confounders (e.g. age, antibiotic use, sexual behaviour, alcohol use, obesity) were well balanced after weighting. The ability of 16S sequencing to provide species-level resolution is limited; thus, we conducted analyses at the genus level. We acknowledge that differentially abundant genera do not necessarily indicate differences in functionally important species. Finally, we did not have time since ART initiation for our cohort, and thus cannot determine whether participants were viremic because they recently started treatment, because of nonadherence, or because treatment was failing.

Our study also has numerous strengths. Primarily, we utilized data from a large cohort of regionally, socioeconomically and behaviourally similar individuals, increasing the internal validity of our findings. We employed a novel technique, IPTW, to incorporate a large amount of clinical and behavioural data into our analyses. With IPTW, analyses are ‘marginal structural models’ instead of conditional on covariates, as in multiple regression. Modelling microbiome data marginally offers several advantages including the ability to control for many confounding factors without inducing overfitting bias [54] or losing efficiency due to overstratification [55]. Addressing sources of systematic error using IPTW may improve reproducibility in future studies of HIV and the microbiome. We also stratified our HIV-infected participants by level of viremia, allowing us to examine differences between HIV-uninfected and HIV-positive undetectable, suppressed and viremic individuals. This stratification leads to better understanding of the relationship between active viral replication and dysbiosis, namely, that dysbiosis increases with increasing viremia. Finally, we were able to replicate our major findings using two distinct analytic strategies. Genera identified as differentially abundant in both analyses may be more likely to be true discoveries.

This study contributes to a growing body of literature describing the effects of HIV on microbial dysbiosis. We show that, even when taking into account the influence of multiple confounding factors, HIV is associated with intestinal dysbiosis in a dose-dependent manner. Although great strides have been made in the management of chronic HIV infection, the life expectancy among HIV-infected individuals remains reduced relative to their HIV-negative peers [56]. This reduction has largely been attributed to increased rates of inflammation-related comorbidities observed among people living with HIV [57], and the microbiome likely plays a key role in modulating interactions between HIV and the immune system. Therefore, understanding the ways in which HIV and the microbiome interact may be a crucial step towards developing intervention strategies to reduce the burden of HIV-associated morbidity and mortality.

Acknowledgements

Concept and design: R.C., J.F., N.H.T., F.L., M.J., R. Brookmeyer, G.M.A., P.M.G. Data acquisition: S.S., R. Bolan, P.M.G., G.M.A. Statistical analysis: R.C., F.L., D.L., R. Brookmeyer. Drafting of manuscript: R.C., J.E., N.H.T., F.L., D.L., M.J., R. Brookmeyer, S.S., R. Bolan, G.M.A., P.M.G. Critical revision of the manuscript: R.C., J.F., N.H.T., F.L., D.L., M.J., R. Brookmeyer, S.S., R. Bolan, G.M.A., P.M.C.

This work was supported by the National Institute on Drug Abuse (1R36 DA046310 and 2U01 DA036267) and the UCLA Center for HIV Identification, Prevention and Treatment Services (CHIPTS; National Institute of Mental Health P30 MH58107). Additional support provided by the UCLA AIDS Institute and UCLA CFAR Microbiome and Mucosal Immunology Core (P30 AI028697). J.A.F. was supported in part by National Institute of Allergy and Infectious Diseases (NIAID) (K08 AI124979). Overall support for the International Maternal Pediatric Adolescent AIDS Clinical Trials Group (IMPAACT) was provided by the National Institute of Allergy and Infectious Diseases (NIAID) of the National Institutes of Health (NIH) under Award Numbers UM1AI068632 (IMPAACT LOC), UM1AI068616 (IMPAACT SDMC) and UM1AI106716 (IMPAACT LC), with cofunding from the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) and the National Institute of Mental Health (NIMH). The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Conflicts of interest

There are no conflicts of interest.

References

1. Lynch SV, Pedersen O. The human intestinal microbiome in health and disease. N Engl J Med 2016; 375:2369–2379.
2. Human Microbiome Project Consortium. Structure, function and diversity of the healthy human microbiome. Nature 2012; 486:207–214.
3. Brenchley JM, Douek DC. Microbial translocation across the GI tract. Annu Rev Immunol 2012; 30:149–173.
4. Marchetti G, Tincati C, Silvestri G. Microbial translocation in the pathogenesis of HIV infection and AIDS. Clin Microbiol Rev 2013; 26:2–18.
5. Dillon SM, Frank DN, Wilson CC. The gut microbiome and HIV-1 pathogenesis: a two-way street. AIDS 2016; 30:2737–2751.
6. Gootenberg DB, Paer JM, Luevano JM, Kwon DS. HIV-associated changes in the enteric microbial community: potential role in loss of homeostasis and development of systemic inflammation. Curr Opin Infect Dis 2017; 30:31–43.
7. Zevin AS, McKinnon L, Burgener A, Klatt NR. Microbial translocation and microbiome dysbiosis in HIV-associated immune activation. Curr Opin HIV AIDS 2016; 11:182–190.
8. d’Ettorre G, Ceccarelli G, Giustini N, Serafino S, Calantone N, De Girolamo G, et al. Probiotics reduce inflammation in antiretroviral treated, HIV-infected individuals: results of the ‘Probio-HIV’ clinical trial. PLoS One 2015; 10:e0137200.
9. Vujkovic-Cvijin I, Dunham RM, Iwai S, Maher MC, Albright RG, Broadhurst MJ, et al. Dysbiosis of the gut microbiota is associated with HIV disease progression and tryptophan catabolism. Sci Transl Med 2013; 5:193ra191.
10. Deeks SG. HIV infection, inflammation, immunosenescence, and aging. Annu Rev Med 2011; 62:141–155.
11. Mutlu EA, Keshavarzian A, Losurdo J, Swanson G, Siewe B, Forsyth C, et al. A compositional look at the human gastrointestinal microbiome and immune activation parameters in HIV infected subjects. PLoS Pathog 2014; 10:e1003829.
12. Yu G, Fadrosh D, Ma B, Ravel J, Goedert JJ. Anal microbiota profiles in HIV-positive and HIV-negative MSM. AIDS 2014; 28:753–760.
13. Nowak P, Troseid M, Avershina E, Barqasho B, Neogi U, Holm K, et al. Gut microbiota diversity predicts immune status in HIV-1 infection. AIDS 2015; 29:2409–2418.
14. Dillon SM, Lee EJ, Kotter CV, Austin GL, Dong Z, Hecht DK, et al. An altered intestinal mucosal microbiome in HIV-1 infection is associated with mucosal and systemic immune activation and endotoxemia. Mucosal Immunol 2014; 7:983–994.
15. Vazquez-Castellanos JF, Serrano-Villar S, Latorre A, Artacho A, Ferrus ML, Madrid N, et al. Altered metabolism of gut microbiota contributes to chronic immune activation in HIV-infected individuals. Mucosal Immunol 2015; 8:760–772.
16. Dinh DM, Volpe GE, Duffalo C, Bhalchandra S, Tai AK, Kane AV, et al. Intestinal microbiota, microbial translocation, and systemic inflammation in chronic HIV infection. J Infect Dis 2015; 211:19–27.
17. 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–357.
18. Lozupone CA, Li M, Campbell TB, Flores SC, Linderman D, Gebert MJ, et al. Alterations in the gut microbiota associated with HIV-1 infection. Cell Host Microbe 2013; 14:329–339.
19. Noguera-Julian M, Rocafort M, Guillen Y, Rivera J, Casadella M, Nowak P, et al. Gut microbiota linked to sexual preference and HIV infection. EBioMedicine 2016; 5:135–146.
20. Hoel H, Hove-Skovsgaard M, Hov JR, Gaardbo JC, Holm K, Kummen M, et al. Impact of HIV and type 2 diabetes on gut microbiota diversity, tryptophan catabolism and endothelial dysfunction. Sci Rep 2018; 8:6725.
21. Pollock J, Glendinning L, Wisedchanwet T, Watson M. The madness of microbiome: attempting to find consensus ‘best practice’ for 16s microbiome studies. Appl Environ Microbiol 2018; 84:pii: e02627-17.
22. Thorsen J, Brejnrod A, Mortensen M, Rasmussen MA, Stokholm J, Al-Soud WA, et al. Large-scale benchmarking reveals false discoveries and count transformation sensitivity in 16S rRNA gene amplicon data analysis methods used in microbiome studies. Microbiome 2016; 4:62.
23. Hawinkel S, Mattiello F, Bijnens L, Thas O. A broken promise: microbiome differential abundance methods do not control the false discovery rate. Brief Bioinform 2017; [Epub ahead of print].
24. Lozupone CA, Rhodes ME, Neff CP, Fontenot AP, Campbell TB, Palmer BE. HIV-induced alteration in gut microbiota: driving factors, consequences, and effects of antiretroviral therapy. Gut Microbes 2014; 5:562–570.
25. McHardy IH, Li X, Tong M, Ruegger P, Jacobs J, Borneman J, et al. HIV Infection is associated with compositional and functional shifts in the rectal mucosal microbiota. Microbiome 2013; 1:26.
26. Manner IW, Baekken M, Kvale D, Oektedalen O, Pedersen M, Nielsen SD, et al. Markers of microbial translocation predict hypertension in HIV-infected individuals. HIV Med 2013; 14:354–361.
27. Vesterbacka J, Rivera J, Noyan K, Parera M, Neogi U, Calle M, et al. Richer gut microbiota with distinct metabolic profile in HIV infected elite controllers. Sci Rep 2017; 7:6269.
28. Bender JM, Li F, Martelly S, Byrt E, Rouzier V, Leo M, et al. Maternal HIV infection influences the microbiome of HIV-uninfected infants. Sci Transl Med 2016; 8:349ra100.
29. Pannaraj PS, Ly M, Cerini C, Saavedra M, Aldrovandi GM, Saboory AA, et al. Shared and distinct features of human milk and infant stool viromes. Front Microbiol 2018; 9:1162.
30. Bender JM, Li F, Adisetiyo H, Lee D, Zabih S, Hung L, et al. Quantification of variation and the impact of biomass in targeted 16S rRNA gene sequencing studies. Microbiome 2018; 6:155.
31. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 2010; 7:335–336.
32. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods 2016; 13:581–583.
33. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 2007; 73:5261–5267.
34. Chen L, Reeve J, Zhang L, Huang S, Wang X, Chen J. GMPR: a robust normalization method for zero-inflated count data with application to microbiome sequencing data. PeerJ 2018; 6:e4600.
35. Robins JM, Hernan MA, Brumback B. Marginal structural models and causal inference in epidemiology. Epidemiology 2000; 11:550–560.
36. Van Buuren S, Brand JPL, Groothuis-Oudshoorn CGM, Rubin DB. Fully conditional specification in multivariate imputation. J Stat Comput Simulat 2006; 76:1049–1064.
37. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc B 1995; 57:289–300.
38. Benjamini Y, Yekutieli D. False discovery rate–adjusted multiple confidence intervals for selected parameters. J Am Stat Assoc 2005; 100:71–81.
39. Sun Y, Ma Y, Lin P, Tang YW, Yang L, Shen Y, et al. Fecal bacterial microbiome diversity in chronic HIV-infected patients in China. Emerg Microbes Infect 2016; 5:e31.
40. Monaco CL, Gootenberg DB, Zhao G, Handley SA, Ghebremichael MS, Lim ES, et al. Altered virome and bacterial microbiome in human immunodeficiency virus-associated acquired immunodeficiency syndrome. Cell Host Microbe 2016; 19:311–322.
41. Dubourg G, Lagier JC, Hue S, Surenaud M, Bachar D, Robert C, et al. Gut microbiota associated with HIV infection is significantly enriched in bacteria tolerant to oxygen. BMJ Open Gastroenterol 2016; 3:e000080.
42. Xu Y, Xie Z, Wang H, Shen Z, Guo Y, Gao Y, et al. Bacterial diversity of intestinal microbiota in patients with substance use disorders revealed by 16S rRNA gene deep sequencing. Sci Rep 2017; 7:3628.
43. Dillon SM, Lee EJ, Kotter CV, Austin GL, Gianella S, Siewe B, et al. Gut dendritic cell activation links an altered colonic microbiome to mucosal and systemic T-cell activation in untreated HIV-1 infection. Mucosal Immunol 2016; 9:24–37.
44. Scher JU, Sczesnak A, Longman RS, Segata N, Ubeda C, Bielski C, et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. Elife 2013; 2:e01202.
45. Troy EB, Kasper DL. Beneficial effects of Bacteroides fragilis polysaccharides on the immune system. Front Biosci (Landmark Ed) 2010; 15:25–34.
46. Mazmanian SK, Round JL, Kasper DL. A microbial symbiosis factor prevents intestinal inflammatory disease. Nature 2008; 453:620–625.
47. Wang PL, Ohura K. Porphyromonas gingivalis lipopolysaccharide signaling in gingival fibroblasts-CD14 and Toll-like receptors. Crit Rev Oral Biol Med 2002; 13:132–142.
48. Arimatsu K, Yamada H, Miyazawa H, Minagawa T, Nakajima M, Ryder MI, et al. Oral pathobiont induces systemic inflammation and metabolic changes associated with alteration of gut microbiota. Sci Rep 2014; 4:4828.
49. Choi J, Chung SW, Kim SJ, Kim SJ. Establishment of Porphyromonas gingivalis-specific T-cell lines from atherosclerosis patients. Oral Microbiol Immunol 2001; 16:316–318.
50. Nakajima M, Arimatsu K, Kato T, Matsuda Y, Minagawa T, Takahashi N, et al. Oral administration of P. gingivalis induces dysbiosis of gut microbiota and impaired barrier function leading to dissemination of enterobacteria to the liver. PLoS One 2015; 10:e0134234.
51. Laprise C, de Pokomandy A, Baril JG, Dufresne S, Trottier H. Virologic failure following persistent low-level viremia in a cohort of HIV-positive patients: results from 12 years of observation. Clin Infect Dis 2013; 57:1489–1496.
52. Bernal E, Gomez JM, Jarrin I, Cano A, Munoz A, Alcaraz A, et al. Low level viremia is associated with clinical progression in HIV-infected patients receiving antiretroviral treatment. J Acquir Immune Defic Syndr 2018; 78:329–337.
53. Reus S, Portilla J, Sanchez-Paya J, Giner L, Frances R, Such J, et al. Low-level HIV viremia is associated with microbial translocation and inflammation. J Acquir Immune Defic Syndr 2013; 62:129–134.
54. Babyak MA. What you see may not be what you get: a brief, nontechnical introduction to overfitting in regression-type models. Psychosom Med 2004; 66:411–421.
55. De Stavola BL, Cox DR. On the consequences of overstratification. Biometrika 2008; 95:992–996.
56. Harrison KM, Song R, Zhang X. Life expectancy after HIV diagnosis based on national HIV surveillance data from 25 states, United States. J Acquir Immune Defic Syndr 2010; 53:124–130.
57. Guaraldi G, Orlando G, Zona S, Menozzi M, Carli F, Garlassi E, et al. Premature age-related comorbidities among HIV-infected persons compared with the general population. Clin Infect Dis 2011; 53:1120–1126.

* Grace M. Aldrovandi and Pamina M. Gorbach contributed equally to the writing of this article.

Keywords:

causal inference; dysbiosis; HIV; inverse probability of treatment weighting; MSM; microbiome; propensity scores

Supplemental Digital Content

Copyright © 2019 Wolters Kluwer Health, Inc.