Oral Microbiota Is Associated With Immune Recovery in Human Immunodeficiency Virus-Infected Individuals

The role of the oral microbiota in HIV-infected individuals deserves attention as either HIV infection or antiretroviral therapy (ART) may have effect on the diversity and the composition of the oral microbiome. However, few studies have addressed the oral microbiota and its interplay with different immune responses to ART in HIV-infected individuals. Salivary microbiota and immune activation were studied in 30 HIV-infected immunological responders (IR) and 34 immunological non-responders (INR) (≥500 and < 200 CD4 + T-cell counts/μl after 2 years of HIV-1 viral suppression, respectively) with no comorbidities. Metagenome sequencing revealed that the IR and the INR group presented similar salivary bacterial richness and diversity. The INR group presented a significantly higher abundance of genus Selenomonas_4, while the IR group manifested higher abundances of Candidatus_Saccharimonas and norank_p_Saccharimonas. Candidatus_Saccharimonas and norank_p_Saccharimonas were positively correlated with the current CD4 + T-cells. Candidatus_Saccharimonas was positively correlated with the markers of adaptive immunity CD4 + CD57 + T-cells, while negative correlation was found between norank _p_Saccharimonas and the CD8 + CD38 + T-cells as well as the CD4/CD8 + HLADR + CD38 + T-cells. The conclusions are that the overall salivary microbiota structure was similar in the immunological responders and immunological non-responders, while there were some taxonomic differences in the salivary bacterial composition. Selenomona_4, Candidatus_Saccharimonas, and norank _p_Saccharimonas might act as important factors of the immune recovery in the immunodeficiency patients, and Candidatus_Saccharimonas could be considered in the future as screening biomarkers for the immune responses in the HIV-infected individuals.

Introduction

The widespread use of potent antiretroviral therapy (ART) which has ability to achieve viral suppression and immune reconstitution has made human immunodeficiency virus (HIV) infection become a chronic manageable disease. ART can improve the immune function of HIV infected subjects, but significant individual difference exists in the extent of immunological recovery. Despite the persistent manifestation of virological suppression after receiving ART, HIV infected individuals with low increases of CD4 + T-cells are considered as immunological non-responders (INR), showing contrast to immunological responders (IR) (Cenderello and De Maria, 2016). An unanimously agreed definition of INR has not been reached by far. Therefore, the acceptable range of the prevalence of INR is from 10 to 40% in the cohorts (Yang et al., 2020). The most restrictive definition of INR refers to patients whose absolute CD4 + T-cells fail to reach 200 cells/μl with an undetectable plasma viral load after 2 years of receiving ART (Florence et al., 2003; Kaufmann et al., 2003; Tsukamoto et al., 2009).

Thanks to high-throughput sequencing, thorough and comprehensive studies of complex human microbiology have been conducted during the past years. In recent years, studies mainly focused on characterizing the impact of HIV infection on host–microbe interactions in the gut (Dillon et al., 2016) and some studies found that the gut microbiota are associated with immune recovery in HIV-infected patients (Ji et al., 2018; Lu et al., 2018; Xie et al., 2021). Oral microbiomes were also reported to have strong associations with human immune system functions, therefore, they were correlated with human immune system diseases such as rheumatic diseases (RDs), acute lymphoblastic leukemia (ALL), and HIV infection (Gao et al., 2018). Alterations in the oral microbiome distinguished individuals with rheumatoid arthritis (RA) from the healthy controls; correlations were shown between these alterations and clinical measures including immune response, and they could be used to stratify individuals based on their responses to the therapy, especially with microbial triggers being implicated in RA (Bellando-Randone et al., 2021). Furthermore, characterization of the oral microbiome in ALL patients demonstrated a structural imbalance of the oral microbiota, indicating the importance of immune status in shaping the structure of the oral microbiota. Although valuable insights on immune status have been presented by these studies, it is surprising that studies which focus on the role of the oral microbiome and their relationship with immune response in HIV are relatively rare. Few research addresses the overall oral microbiota structure in HIV infection; to our knowledge, the oral microbiota dysbiosis and the interaction with different immune responses to ART is still poorly defined. This study aims to evaluate the microbial composition in the salivary samples collected from HIV-infected immunological non-responders and immunological responders. We hypothesize that the compositional changes of the salivary microbiota could be associated with different immune responses of HIV-infected individuals receiving ART. The study adopted 16S ribosomal RNA (rRNA) targeted sequencing and flow cytometry to explore the oral microbiome and their relationship with immune activation in patients who are immunodiscordant and who are immunoconcordant. The compositional changes of salivary microbiota and their association with different immune responses in HIV-infected patients with ART is characterized for the first time in this study.

Materials and Methods

Recruitment of Subjects

64 HIV-infected individuals who were diagnosed by the Disease Control and Prevention Center of Zhejiang Province (30 immunological responders and 34 immunological non-responders) were all recruited from the HIV clinic of the First Affiliated Hospital of Zhejiang University from November 2015 to October 2017. All subjects start ART during the chronic phase of HIV infection. In this study, IR and INR were defined as patients who with the average of the last two CD4 + T-cell counts/μl equal or is more than 500 or less than 200 and after 2 years of receiving complete viral suppression therapy, respectively. The selection excluded candidates with either one of the following conditions: below 18 years old; showing opportunistic infection symptoms; infected with hepatitis B or C; having history of using antibiotics, immunosuppressive regimen, probiotics, prebiotics, or symbiotics in the past 6 months; BMI higher than 30; showing oral active inflammation. None of the patients had obvious symptoms of oral mucosal diseases and periodontal disease (redness, swelling, and bleeding) when the clinical samples were collected. However, formal dental examinations were not performed to rule out the mild periodontal symptoms.

Ethics Statement

This study conforms to the ethical norms of the 1975 Helsinki Declaration. The research protocol was approved by the Institutional Review Committee of The First Affiliated Hospital of Zhejiang University on October 7, 2015. All participants provided written informed consents before participating in the study. All the data used for analysis were anonymized.

Salivary Samples Collection and DNA Extraction

Participants were required to refrain from eating, drinking, smoking before saliva collection. The amount of each salivary sample collected from participants before their clinic visit was 5 ml. The samples were stored in sterile containers of –80°C until DNA extraction by QiaAmp DNA Mini Kit (QIAGEN, Hilden, Germany) following instructions of the manufacturer. NanoDrop (Thermo Fisher Scientific) was used to determine the concentration and purity of DNA as well as 1.0% agarose gel electrophoresis for the integrity and size of DNA. After the procedures above, the DNA samples were frozen at –20°C for further analysis.

16S Ribosomal RNA Gene Sequencing

The bacterial 16S rRNA gene high-throughput sequencing was conducted by Shanghai Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). The bacterial 16S rRNA gene sequences spanning the variable regions V3–V4 were amplified using the primer 338F (5′- ACTCCTACGGGAGGCAGCAG-3′), and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) as previous recorded (Xie et al., 2021). The amplicons were extracted from 2% agarose gels, purified by the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) and quantified using QuantiFluor™-ST (Promega, United States) based on the guidelines of the manufacturer’s protocol. In equimolar amounts, purified amplicons were sent to paired-end sequencing (2 × 300) on an Illumina MiSeq platform (Illumina, San Diego, United States).

Bacterial Translocation, Viral Load, and Flow Cytometry

Sera samples of immunological responders and non-responders were collected to measure the bacterial translocation markers. Following the standard protocols, human Lipopolysaccharides (LPS) ELISA Kit (CUSABIO; Wuhan, China) and Human soluble CD14 (sCD14) ELISA Kit (MultiSciences, Hangzhou, China) were used to test plasma LPS and sCD14. Flow cytometry and Cobas Amplicor (Roche Molecular Systems Inc., Branchburg, New Jersey, United States) were used to quantify CD4 + /CD8 + T-cells and HIV-1 RNA, respectively. Fresh anticoagulated whole blood was used to quantify the expressing markers of immune activation (CD25 +, CD38 +, HLADR +, or CD38 + /HLA-DR +) of CD4 + and CD8 + T-cells and immune senescence (CD57 +) by BD FACS Canto II flow cytometer (BD Biosciences, California, United States). The antibodies needed during the experiment were purchased from Biolegend (San Diego, CA), including CD3-FITC, CD4- PerCP/Cy5.5, CD8-Brilliant Violet 510™, CD38-Brilliant Violet 421, CD25-PE, HLA-DR-APC/Fire™ 750, and CD57-allophycocyanin (APC).

Bioinformatics and Statistics

The 16S rRNA high-throughput sequencing raw fastq files were demultiplexed and quality-filtered by QIIME (version 1.9.1).¹ Operational taxonomic units (OTUs) were clustered with a 97% similarity cut-off using UPARSE.² The taxonomy of each 16S rRNA gene sequence was analyzed using the RDP Classifier (version 2.2)³ and compared with the SILVA rRNA database⁴ with the confidence threshold being 70%. After eliminating the interference sequence, alpha diversity estimator calculations were performed using Mothur v.1.30.2. Phylogenetic beta diversity measures, such as Bray-Curtis distance metrics analysis; the representative sequences of OTUs were used for each sample, respectively. Principal Coordinates analysis (PCoA) was performed to visualize the microbial communities following the distance matrices, as well as Linear discriminant analysis effect size (LEfSe) on Galaxy to calculate bacteria taxa with significantly different abundances between groups (Segata et al., 2011). In this study, alpha values for the factorial Kruskal Wallis at 0.05 and linear discriminant analysis (LDA) effect size threshold of 2.0 for discriminative features were applied for all bacteria that were discussed. The calculation of correlations between the variables were conducted with Spearman’s rank-correlation analysis. Spearman correlation matrix with p-adjust lower than 0.05 and ρ-value above 0.2 were used to filter strong correlations. The discriminatory function of the biomarkers was evaluated through calculating the area under the receiver operating characteristic (ROC) curve (AUC) using pROC of R package. The comparisons between groups were conducted through the Chi-square test, Independent-Samples T-test, Wilcoxon rank sum test and Mann-Whitney U-test in the R package and SPSS 21.0 software (SPSS Inc., Chicago, IL, United States). Differences were considered significant when P < 0.05.

Results

General Clinical Features of the Patients

The cross-sectional study subjects included 30 immunological responders and 34 immunological non-responders with HIV infection. The characteristics such as gender, age, body mass index (BMI), and smoking status are relatively matched between the two groups (Table 1). The MSM transmission route rate is 53.3% vs. 44.1% (p = 0.461).

TABLE 1

Table 1. Clinical characteristics data summary.

ART medications were composed by two Nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs) and a non-nucleoside reverse transcriptase inhibitor (NNRTI), or a ritonavir-boosted protease inhibitor (PI). No differences are observed in the duration of ART and the type of ART drugs between the IR and the INR groups (p = 0.193 and p = 0.821). The HIV RNA viral load levels were considered as undetectable (<20 copies/ml) in all samples of the patients. In the IR group, the number of Nadir CD4 + T cell is greatly higher than that in the INR group (298 vs. 42.5, p < 0.001) (Table 1).

Salivary Microbiota Analysis in the Immunological Responders and Immunological Non-responders Groups

To characterize the salivary microbiota composition, 3,884,775 high-quality16S rRNA sequences were obtained from all 64 participants. An average length of 445 bp and an average of 60,700 sequences per sample were adopted for further analysis. Rarefaction was conducted on the OTU table to 27,986 reads per sample to avoid any methodological artifacts. Specifically, 504 OTUs in the IR group and 503 OTUs in the INR group were defined at a similarity level of 97%. The diversity (Sobs index, Simpson’s index of diversity and Shannon index) and richness (Chao1, ACE estimator and Good’s coverage) of the salivary microbiota in each group at the OTU level were analyzed, and the IR group present similar Alpha and Beta diversity compared to the INR group (Table 2). It can be seen from the principal coordinate (PCoA) analysis by Bray-Curtis matrices that there is no significant difference between the two groups (PERMANOVA, pseudo-F: 1.64608, R² = 0.02586, p = 0.109, Figure 1). Visualization of the relative abundances of dominant taxa at the genus level in the oral microbiome is presented in Figure 2. In order to identify the key phylotypes responsible for the difference in distinguishing the saliva microbiota between the two groups, linear discriminant analysis (LDA) effect size (LEfSe) was performed and the threshold was 2. The INR group present a significantly higher abundance of genus Selenomonas_4, while the IR group has higher abundances of genus Candidatus_Saccharimonas, norank _p_Saccharimonas, and Desulfobulbus (Figure 3). The Wilcoxon Rank Sum test was also used to detect the taxa with significant differences in the relative abundances between the two groups (using confidence interval method). The abundance of Saccharimonas is more abundant in the IR group than the INR group at the phylum level (Figure 4A). Compared with the INR group, the abundances of genus Candidatus_Saccharimonas and norank_p_Saccharimonas are dramatically increased in the IR group, while the abundance of genus Selenomonas_4 is dramatically decreased in the IR group (Figure 4B).

TABLE 2

Table 2. Salivary microbiota 16S rRNA gene high-throughput sequencing data summary.

FIGURE 1

Figure 1. Principal coordinates analysis (PCoA) of salivary microbiota in the IR and INR groups. No significant difference of bacterial communities between the immunological responders (IR) and immunological non-responders (INR) groups. IR, immunological responders; INR, immunological non-responders.

FIGURE 2

Figure 2. The salivary microbiota composition at genus level in the IR and INR groups. IR, immunological responders; INR, immunological non-responders.

FIGURE 3

Figure 3. Linear discriminative analysis (LDA) effect size (LEfSe) analysis in the IR and INR groups. LDA scores for the significant taxa in the immunological non-responders (INR) group are represented on the positive scale (red), and LDA-negative scores represent enriched taxa in the immunological responders (IR) group (blue). IR, immunological responders; INR, immunological non-responders.

FIGURE 4

Figure 4. The difference of saliva microbial taxa between the IR and INR groups at the phylum and genus levels. The Wilcoxon Rank Sum test was performed to detect taxa with significant differences in relative abundances at the phylum (A) and the genus levels (B) between the two groups (using confidence interval method). IR, immunological responders; INR, immunological non-responders. *P < 0.05 and ^**P < 0.01.

Comparison of Adaptive Immunity and Bacterial Translocation Markers in the Immunological Responders and Immunological Non-responders Groups

As expected, the amount of the current CD4 + T-cell counts and the CD4/CD8 ratio in the INR group is lower than those in the IR group (p < 0.001). The proportion of the CD8 + CD38 + T-cell and the CD8 + CD57 + T-cell is significantly higher in the INR group than those in the IR group (p = 0.032 and p = 0.001). The proportion level of the CD4 + immune activation (CD4 + T-cell by the expression of CD25 +, HLA-DR +, CD38 +, or HLA-DR + /CD38 +) shows similarity in the INR and IR groups. Lipopolysaccharide (LPS), which is commonly used as the major antigens driving chronic immune activation, is significantly higher in the INR group compared with that in the IR group (p = 0.027) (Table 3).

TABLE 3

Table 3. Comparison of adaptive immunity and bacterial translocation markers.

Associations Between Oral Microbiome and Adaptive Immunity

The Spearman’s correlation test was used to investigate the correlation between the relative abundance of the different genera and adaptive immunity markers. The p-value was corrected by FDR, p-adjust lower than 0.05 with ρ-value above 0.2 were considered relevant and showed in Figure 5. As the abundance of genus Selenomonas_4 is low in the two groups, the multiple correlation analyses were not performed. Candidatus_Saccharimonas and norank _p_Saccharimonas are positively correlated with the current CD4 + T cells (p-adjust = 0.033 and 0.025, respectively). Candidatus_Saccharimonas is positively correlated with the markers of the adaptive immunity CD4 + CD57 + T-cells (p-adjust = 0.049), while norank _p_Saccharimonas is negatively correlated with the CD8 + CD38 + T-cells and the CD4/CD8 + HLADR + CD38 + T-cells, respectively (p-adjust = 0.001 and 0.032, respectively) (Figure 5). To explore the potential function of the saliva microbiome for discriminating the IR and INR status, a random forest model was created based on the microbiome and the top 5 genus was shown in Figure 6A. The ROC analysis shows that Candidatus_Saccharimonas could be used to discriminate the IR from the INR group [ROC-plot AUC value of 0.7 (95% CI, 0.56–0.83), Figure 6B].

FIGURE 5

Figure 5. Two-way correlation network analysis between the saliva microbiota and the adaptive immunity markers. Some adaptive immunity markers are correlated with the specific genera of saliva microbiota. Positive correlation is shown in the red line, negative correlation is shown in the green line. Spearman’s correlation was used, and associations with p-adjust lower than 0.05 and ρ-value above 0.2 were considered relevant.

FIGURE 6

Figure 6. Candidatus_Saccharimonas could be used to discriminate the IR from the INR group. Barplot of variable importance on the genus level (A) and the receiving operational curve analysis were performed (B) [area under curve (AUC) = 0.7]. Diagonal lines represent random classifications (AUC = 0.5).

Discussion

It is well-known that there is a dynamic interaction between the host and the microbiota, which is a major factor of people’s health (Dethlefsen et al., 2007). Current available studies have confirmed that discernible alterations of the composition of the salivary microbiota are inherent to a range of systemic disorders. Studies on oral microbiomes of human immune system diseases such as RA, ALL, and HIV have also indicated that the oral microbiota is strongly related to the immune responses in immunodeficiency patients (Acharya et al., 2017; Corrêa et al., 2017; Annavajhala et al., 2020; Bellando-Randone et al., 2021). Furthermore, the hypothesis that oral microbiota is associated with different immune responses to ART is supported by the high frequency of opportunistic oral infections in HIV-infected patients and its association with CD4 + T cells levels (Berberi and Noujeim, 2015). As stated, the role of the oral microbiota played in HIV-infected patients deserves much attention because the diversity and composition of the oral microbiome can be changed by HIV infection or by ART (Heron and Elahi, 2017); however, the results are highly variable. Previous studies have observed the generally similar structure of salivary microbiota as well as the differences in the relative abundances of several bacterial taxa between HIV-infected subjects and uninfected controls (Li et al., 2014; Kistler et al., 2015; Mukherjee et al., 2018; Lewy et al., 2019). However, significant differences were also found by several studies in the saliva bacterial communities between HIV-infected and uninfected individuals (Hegde et al., 2014; Mukherjee et al., 2014; Beck et al., 2015; Goldberg et al., 2015; Kistler et al., 2015; Li et al., 2021). What’s more, some studies reported significant distinctions in the prevalence and the distribution of the saliva bacterial communities among HIV-infected individuals before and after the antiretroviral therapy (Li et al., 2014; Kistler et al., 2015; Mukherjee et al., 2018). Few studies have addressed the interaction of oral microbiota with different immune responses to ART received by the HIV-infected individuals. In this paper, we demonstrate the oral microbiota structure and its relationship with immune response in HIV infection.

A previous study indicated that after receiving 24 weeks of ART, the salivary microbiome in the three HIV-infected participants with persistently low CD4 + T-cell count had significantly higher bacterial richness and Shannon diversity; when compared to those with CD4 counts that remained or recovered to greater than 200 cells/μl. Several taxa with different abundancies, such as Porphyromonas species, discriminated between the baseline and the posttreatment samples; this suggested that the salivary microbiome can be an important factor in the CD4 + T-cell count recovery after ART in the study (Presti et al., 2018). However, the major limitation of this study is the insufficient number of studied subjects and the short time period of patients receiving ART. Thus, more studies on larger cohorts are necessary for a better understanding of the potential roles of different immune responses to ART on oral microbiome. Making the most use of the non-invasive and unsuspicious functions of saliva sampling, we collected the salivary samples from the HIV-infected immunological non-responders and immunological responders and studied the salivary microbiome using high-throughput sequencing technology. To our knowledge, this is the first study that has fully utilized the next-generation sequencing technology to characterize and compare the community composition of the salivary microbiota in the higher number of HIV-infected immunological responders and non-responders. These findings are particularly important given that the adaptive immunity markers assessments were performed on all participants, offering an opportunity to evaluate the relationship between the salivary microbiota and the immunologic markers in the immunological responders and non-responders.

The IR group presented similar salivary bacterial richness and diversity when compared with the INR group in this study. A similar Alpha diversity of salivary bacterial community was also described in the HIV-infected IR (n = 18) and INR (n = 9) individuals in a previous study (Jiménez-Hernández et al., 2019). While another study reported that the three HIV-infected participants with persistently low CD4 + T-cell counts had significantly higher salivary bacterial richness and Shannon diversity after their 24 weeks of ART, when compared with the participants whose CD4 + T-cell counts higher than 200 cells/μl (Presti et al., 2018). It claimed that the HIV infection and highly active antiretroviral therapy (HAART) had significant effects on salivary microbial colonization and composition, and Selenomonas noticeably increased after HAART (Li et al., 2014), which were reported to be depleted in the oral microbiome of the HIV-associated periodontitis (Noguera-Julian et al., 2017). For patients with systemic lupus erythematosus (SLE), increased numbers of Selenomonas were directly correlated with the elevated levels of inflammatory cytokines IL-6, IL-17, and IL-33 (Corrêa et al., 2017). In this study, when we visualized the relative abundances of dominant taxa at the genus level in the oral microbiome, we found that the INR group presented a significantly higher abundance of genus Selenomonas_4. Gut Candidatus_Saccharimonas was reported decreased in rats acute necrotizing pancreatitis (Chen et al., 2017), which indicated that Candidatus_Saccharimonas did play an important role in maintaining normal intestinal functions. A previously published study stated that the relative abundance of gut Candidatus_Saccharimonas was negatively correlated with the expression levels of cadherin-11, IL-17α, and TLR2 in the adjuvant-induced arthritis rat model (Huang et al., 2019). In this study, the IR group had higher abundances of genus Candidatus_Saccharimonas and norank _p_Saccharimonas when compared with the INR group. In addition, Candidatus_Saccharimonas had positive correlation with the CD4 + T-cells and CD4 + CD57 + T-cells, while norank _p_Saccharimonas was positively correlated with the CD4 + T-cells but negatively correlated with the CD8 + CD38 + T-cells and CD4/CD8 + HLADR + CD38 + T-cells. While host biomarkers were subjected to the individual biological variations, oral microbiome was relatively conserved among unrelated individuals. Recent advances of saliva analysis have played a key role in the definitions of biomarkers for the diagnosis, prognosis, and the treatment of human immune system diseases (Bellando-Randone et al., 2021). Expansion of specific microbial consortia in the saliva may act as imprints of the underlying immuno-inflammatory processes, especially in HIV. This study indicated that Selenomona, Candidatus_Saccharimonas, and norank _p_Saccharimonas might all played important roles in the immune recovery of the immunodeficiency patients, and Candidatus_Saccharimonas could be considered in the future as screening biomarkers for immune responses in HIV-infected individuals, leading to the future design of effective individualized treatment strategies such as probiotics for the immunological non-responders. However, a previous study on the effect of the prebiotic modulation of the salivary microbiota in HIV-infected patients with diverse immunopathogenesis stated that: Streptococcus anginosus was in correlation with the CD4 + T cells, Veillonella parvula with the CD4 + CD25 + T cells, and Prevotella pallens, Prevotella copri, and Prevotella nigrencens with the markers of adaptive immunity such as CD4 + CD25 + T cells, CD4 + CD57 + T cells, and CD4 + HLADR + CD38 + T cells, respectively (Jiménez-Hernández et al., 2019). The inconsistent results may be caused by the different number of subjects in the two studies and their different definitions of INR. More studies need to be carried out to determine the precise cause of this inconsistency.

This study provides a better understanding of the oral microbial profiles and their relationship with the immune responses in the immunocompromised patients. Although the correlations between the salivary microbiota and the biomarkers of adaptive immunity were observed, we could not establish a causal influence between the oral microbiome and the immune system in the HIV-infected patients. Thus, further research on the exact mechanisms involved in the interaction between the immune system and oral microbiota is still required. We find the low number of participants, absence of control group, and cross-sectional analysis caused limitations to the explanation of our findings. Intra-person variability will need larger longitudinal studies with control group, which should involve plaque collection and the observation of the clinical and environmental changes that account for intra-person variability over time. Apart from taxonomic characterization, there should be more studies on identifying bacterial pathways and the resulting metabolites that promote disease and immunity.

In summary, this study focuses on the overall oral microbiota structure and its interactions with different immune responses to ART. While there were some taxonomic differences in the salivary bacterial composition, the results suggested that the overall structure of the salivary microbiota in the immunological responders was similar with those that were in the immunological non-responders. Selenomona_4, Candidatus_Saccharimonas, and norank_p_Saccharimonas might act as important factors of the immune recovery in the immunodeficiency patients, and Candidatus_Saccharimonas could be considered in the future as screening biomarkers for immune responses in the HIV-infected individuals.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: The raw reads were deposited into the NCBI Sequence Read Archive (SRA) database and is accessible with the following link: https://dataview.ncbi.nlm.nih.gov/object/PRJNA744583?reviewer=3gmrc9a67n88rhe2gr9ehub4ok.

Ethics Statement

The research protocol was approved by the Institutional Review Committee of The First Affiliated Hospital of Zhejiang University on October 7, 2015. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

YX participated in the designing of the study and wrote the manuscript. JS performed the statistical analysis and revised the manuscript. JS and CH collected the biopsy samples and carried out the experiment. BR and BZ participated in the designing and reviewing of the manuscript. All authors read through and approved the final manuscript.

Funding

This study was funded by the National Natural Science Foundation of China (Young Scientist Fund, 81500491) and the National Natural Science Foundation of China (Major Research Plan, 2018ZX10715-014).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

This study gratefully acknowledges the patients who participated in the research.

Abbreviations

HIV, human immunodeficiency virus; ART, antiretroviral therapy; HAART, highly active antiretroviral therapy; INR, immunological non-responders; IR, immunological responders; SLE, systemic lupus erythematosus; RDs, rheumatic diseases; ALL, acute lymphoblastic leukemia; RA, rheumatoid arthritis; BMI, body mass index; MSM, men who have sex with men; LPS, Lipopolysaccharide; PCoA, Principal coordinate analysis; LDA, Linear discriminant analysis; LEfSe, linear discriminant analysis effect size; OUT, operational taxonomic unit; 16S rRNA, 16S ribosomal RNA; ROC, receiver operating characteristic; NRTIs, nucleoside reverse transcriptase inhibitors; NNRTIs, non-nucleoside reverse transcriptase inhibitors; AZT, Zidovudine; TDF, Tenofovir Disoproxil Fumarate; 3TC, Lamivudine; EFV, Efavirenz; LPV/r, Lopinavir/ritonavir.

Footnotes

1. ^ http://qiime.org/install/index.html
2. ^ http://drive5.com/uparse/
3. ^ http://rdp.cme.msu.edu/
4. ^ http://www.arb-silva.de

References

Acharya, A., Chan, Y., Kheur, S., Jin, L. J., Watt, R. M., and Mattheos, N. (2017). Salivary microbiome in non-oral disease: a summary of evidence and commentary. Arch. Oral Biol. 83, 169–173. doi: 10.1016/j.archoralbio.2017.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Annavajhala, M. K., Khan, S. D., Sullivan, S. B., Shah, J., Pass, L., Kister, K., et al. (2020). Oral and gut microbial diversity and immune regulation in patients with HIV on antiretroviral therapy. mSphere 5, e798–e719. doi: 10.1128/msphere.00798-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Beck, J. M., Schloss, P. D., Venkataraman, A., Twigg, H., Jablonski, K. A., Bushman, F. D., et al. (2015). Multicenter comparison of lung and oral microbiomes of HIV-infected and HIV-uninfected individuals. Am. J. Respir. Crit. Care Med. 192, 1335–1344. doi: 10.1164/rccm.201501-0128OC

PubMed Abstract | CrossRef Full Text | Google Scholar

Bellando-Randone, S., Russo, E., Venerito, V., Matucci-Cerinic, M., Iannone, F., Tangaro, S., et al. (2021). Exploring the oral microbiome in rheumatic diseases, state of art and future prospective in personalized medicine with an ai approach. J. Pers. Med. 11:625. doi: 10.3390/jpm11070625

PubMed Abstract | CrossRef Full Text | Google Scholar

Berberi, A., and Noujeim, Z. (2015). Epidemiology and Relationships between CD4+ Counts and Oral Lesions among 50 patients infected with human immunodeficiency virus. J. Int. Oral Heal. JIOH 7, 18–21.

Google Scholar

Cenderello, G., and De Maria, A. (2016). Discordant responses to cART in HIV-1 patients in the era of high potency antiretroviral drugs: clinical evaluation, classification, management prospects. Expert Rev. Anti. Infect. Ther. 14, 29–40. doi: 10.1586/14787210.2016.1106937

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Huang, C., Wang, J., Zhou, H., Lu, Y., Lou, L., et al. (2017). Dysbiosis of intestinal microbiota and decrease in paneth cell antimicrobial peptide level during acute necrotizing pancreatitis in rats. PLoS One 12:e0176583. doi: 10.1371/journal.pone.0176583

PubMed Abstract | CrossRef Full Text | Google Scholar

Corrêa, J. D., Calderaro, D. C., Ferreira, G. A., Mendonça, S. M. S., Fernandes, G. R., Xiao, E., et al. (2017). Subgingival microbiota dysbiosis in systemic lupus erythematosus: association with periodontal status. Microbiome 5:34. doi: 10.1186/s40168-017-0252-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Dethlefsen, L., McFall-Ngai, M., and Relman, D. A. (2007). An ecological and evolutionary perspective on humang-microbe mutualism and disease. Nature 449, 811–818. doi: 10.1038/nature06245

PubMed Abstract | CrossRef Full Text | Google Scholar

Dillon, S. M., Frank, D. N., and Wilson, C. C. (2016). The gut microbiome and HIV-1 pathogenesis: a two-way street. AIDS 30, 2737–2751. doi: 10.1097/QAD.0000000000001289

PubMed Abstract | CrossRef Full Text | Google Scholar

Florence, E., Lundgren, J., Dreezen, C., Fisher, M., Kirk, O., Blaxhult, A., et al. (2003). Factors associated with a reduced CD4 lymphocyte count response to HAART despite full viral suppression in the EuroSIDA study. HIV Med. 4, 255–262. doi: 10.1046/j.1468-1293.2003.00156.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, L., Xu, T., Huang, G., Jiang, S., Gu, Y., and Chen, F. (2018). Oral microbiomes: more and more importance in oral cavity and whole body. Protein Cell 9, 488–500. doi: 10.1007/s13238-018-0548-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, B. E., Mongodin, E. F., Jones, C. E., Chung, M., Fraser, C. M., Tate, A., et al. (2015). The oral bacterial communities of children with well-controlled HIV infection and without HIV infection. PLoS One 10:e0131615. doi: 10.1371/journal.pone.0131615

PubMed Abstract | CrossRef Full Text | Google Scholar

Hegde, M. C., Kumar, A., Bhat, G., and Sreedharan, S. (2014). Oral microflora: a comparative study in HIV and normal patients. Indian J. Otolaryngol. Head Neck Surg. 66, 126–132. doi: 10.1007/s12070-011-0370-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Heron, S. E., and Elahi, S. (2017). HIV infection and compromised mucosal immunity: oral manifestations and systemic inflammation. Front. Immunol. 8:241. doi: 10.3389/fimmu.2017.00241

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Y., Li, M., Zhou, L., Xu, D., Qian, F., Zhang, J., et al. (2019). Effects of qingluo tongbi decoction on gut flora of rats with adjuvant-induced arthritis and the underlying mechanism. Evidence-Based Complement. Altern. Med. 2019:6308021. doi: 10.1155/2019/6308021

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, Y., Zhang, F., Zhang, R., Shen, Y., Liu, L., Wang, J., et al. (2018). Changes in intestinal microbiota in HIV-1-infected subjects following cART initiation: influence of CD4+ T cell count article. Emerg. Microbes Infect. 7:113. doi: 10.1038/s41426-018-0117-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiménez-Hernández, N., Serrano-Villar, S., Domingo, A., Pons, X., Artacho, A., Estrada, V., et al. (2019). Modulation of saliva microbiota through prebiotic intervention in HIV-Infected individuals. Nutrients 11:1346. doi: 10.3390/nu11061346

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaufmann, G. R., Perrin, L., Pantaleo, G., Opravil, M., Furrer, H., Telenti, A., et al. (2003). CD4 T-lymphocyte recovery in individuals with advanced HIV-1 infection receiving potent antiretroviral therapy for 4 years: the Swiss HIV cohort study. Arch. Intern. Med. 163, 2187–2195. doi: 10.1001/archinte.163.18.2187

PubMed Abstract | CrossRef Full Text | Google Scholar

Kistler, J. O., Arirachakaran, P., Poovorawan, Y., Dahlén, G., and Wade, W. G. (2015). The oral microbiome in human immunodeficiency virus (HIV)-positive individuals. J. Med. Microbiol. 64, 1094–1101. doi: 10.1099/jmm.0.000128

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewy, T., Hong, B. Y., Weiser, B., Burger, H., Tremain, A., Weinstock, G., et al. (2019). Oral microbiome in HIV-infected women: shifts in the abundance of pathogenic and beneficial bacteria are associated with aging, HIV load, CD4 count, and antiretroviral therapy. AIDS Res. Hum. Retroviruses 35, 276–286. doi: 10.1089/aid.2017.0200

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Zhu, J., Su, B., Wei, H., Chen, F., Liu, H., et al. (2021). Alteration in oral microbiome among men who have sex with men with acute and chronic HIV infection on antiretroviral therapy. Front. Cell. Infect. Microbiol. 11:695515. doi: 10.3389/fcimb.2021.695515

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Saxena, D., Chen, Z., Liu, G., Abrams, W. R., Phelan, J. A., et al. (2014). HIV infection and microbial diversity in saliva. J. Clin. Microbiol. 52, 1400–1411. doi: 10.1128/JCM.02954-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, W., Feng, Y., Jing, F., Han, Y., Lyu, N., Liu, F., et al. (2018). Association between gut microbiota and CD4 recovery in HIV-1 infected patients. Front. Microbiol. 9:1451. doi: 10.3389/fmicb.2018.01451

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukherjee, P. K., Chandra, J., Retuerto, M., Sikaroodi, M., Brown, R. E., Jurevic, R., et al. (2014). Oral mycobiome analysis of HIV-infected patients: identification of pichia as an antagonist of opportunistic fungi. PLoS Pathog. 10:e1003996. doi: 10.1371/journal.ppat.1003996

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukherjee, P. K., Chandra, J., Retuerto, M., Tatsuoka, C., Ghannoum, M. A., and McComsey, G. A. (2018). Dysbiosis in the oral bacterial and fungal microbiome of HIV-infected subjects is associated with clinical and immunologic variables of HIV infection. PLoS One 13:e0200285. doi: 10.1371/journal.pone.0200285

PubMed Abstract | CrossRef Full Text | Google Scholar

Noguera-Julian, M., Guillén, Y., Peterson, J., Reznik, D., Harris, E. V., Joseph, S. J., et al. (2017). Oral microbiome in HIV-associated periodontitis. Medicine (United States) 96:e5821. doi: 10.1097/MD.0000000000005821

PubMed Abstract | CrossRef Full Text | Google Scholar

Presti, R. M., Handley, S. A., Droit, L., Ghannoum, M., Jacobson, M., Shiboski, C. H., et al. (2018). Alterations in the oral microbiome in HIV-infected participants after antiretroviral therapy administration are influenced by immune status. AIDS 32, 1279–1287. doi: 10.1097/QAD.0000000000001811

PubMed Abstract | CrossRef Full Text | Google Scholar

Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsukamoto, H., Clise-Dwyer, K., Huston, G. E., Duso, D. K., Buck, A. L., Johnson, L. L., et al. (2009). Age-associated increase in lifespan of naïve CD4 T cells contributes to T-cell homeostasis but facilitates development of functional defects. Proc. Natl. Acad. Sci. U.S.A. 106, 18333–18338. doi: 10.1073/pnas.0910139106

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, Y., Sun, J., Wei, L., Jiang, H., Hu, C., Yang, J., et al. (2021). Altered gut microbiota correlate with different immune responses to HAART in HIV-infected individuals. BMC Microbiol. 21:11. doi: 10.1186/s12866-020-02074-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Su, B., Zhang, X., Liu, Y., Wu, H., and Zhang, T. (2020). Incomplete immune reconstitution in HIV/AIDS patients on antiretroviral therapy: challenges of immunological non-responders. J. Leukoc. Biol. 107, 597–612. doi: 10.1002/JLB.4MR1019-189R

PubMed Abstract | CrossRef Full Text | Google Scholar