Beyond the effects of HIV infection and integrase inhibitors

Aug 24, 2023

Scientific Reports volume 13, Article number: 14327 (2023) Cite this article

Metrics details

Oral microbiome is the second largest microbial community in humans after gut. Human immunodeficiency virus (HIV) infection triggers an impairment of the immune system which could favour the growth and the colonization of pathogens in the oral cavity, and this dysbiosis has been associated with oral manifestations that worsen the quality of life of these patients. Antiretroviral therapy (ART) could also drive changes in specific oral bacterial taxa associated with such periodontal diseases. Integrase strand transfer inhibitors (INSTIs), therapy of choice in the treatment of naive HIV-patients, are able to reverse the impact of HIV infection on systemic inflammation, gut permeability, and gut bacterial diversity/richness. The objective of this study was to analyse the effects of HIV infection per se and INSTIs on salivary bacteriome composition, taking into consideration other factors such as smoking, that could also have a significant impact on oral microbiome. To accomplish this objective, 26 non-HIV-infected volunteers and 30 HIV-infected patients (15 naive and 15 under INSTIs-regimen) were recruited. Salivary samples were collected to measure lysozyme levels. Oral bacteriome composition was analysed using 16S rRNA gene sequencing. Naive HIV-infected patients showed statistically higher levels of lysozyme compared to controls (p < 0.001) and INSTIs-treated patients (p < 0.05). Our study was unable to detect differences in α nor β-diversity among the three groups analysed, although significant differences in the abundance of some bacterial taxonomical orders were detected (higher abundance in the phylum Pseudomonadota, in the order Acholeplasmatales, and in the genera Ezakiella and Acholeplasma in the naive group compared to controls; and higher abundance in the phylum Mycoplasmatota, in the order Acholeplasmatales, and in the genera Acholeplasma and uncultured Eubacteriaceae bacterium in the INTIs-treated HIV-infected patients compared to controls). These differences seem to be partially independent of smoking habit. HIV infection and INSTIs effects on oral microbiota seem not to be very potent, probably due to the modulation of other factors such as smoking and the greatest outward exposure of the oral cavity.

HIV infection triggers an impairment of the immune system which could favour the growth and the colonization of pathogens in the oral cavity. In fact, this dysbiosis has been associated with oral manifestations such as oropharyngeal candidiasis1,2,3. Besides, antiretroviral treatment (ART) itself can increase oral microbial translocation due to the inhibition of epithelial cell repair and proliferation4, 5. HIV and ART drive changes in specific bacterial taxa in oral microbiome associated with periodontal disease6. In fact, it has been recently reported that shifts in oral microbiome after ART initiation are complex, and may play an important role in immune function and inflammatory disease6. The analysis of oral microbiota through supra and subgingival plaque samples is often time-consuming7, 8. Because of that, saliva-based analysis has gained considerable attention since it is simple, non-invasive, and inexpensive9,10,11.

Previous works from our group demonstrated that ARTs based on integrase strand transfer inhibitors (INSTIs) were associated with levels of systemic inflammation, bacterial translocation, and microbial diversity similar to those observed in uninfected controls12. Besides, these studies showed a clear impact of HIV-infection and INSTIs-based treatments on gut bacteriome13 and virome14. However, how long-term ART, and specifically INSTIs, modulate oral microbiota in HIV-infected patients and the implications of these effects on health deserve further research. Thus, the objective of this work was to analyse the effects of HIV infection and INSTIs-based therapies in first line of treatment on salivary bacteriome composition. Besides, seeing that some habits such as smoking could alter oral microbiota15, we also investigated the impact of smoking in oral microbiota composition of HIV-infected patients with and without INSTIs-based treatments.

HIV-infected patients (naive and under ART) as well as the non-infected volunteers used as “healthy” controls were recruited from the Infectious Diseases Department at Hospital Universitario San Pedro (HUSP) and from the Health Care Center “Siete Infantes de Lara” (Logroño, Spain) from March 2019 to February 2021, as previously described13, 14. Supplementary Fig. 1 shows a flowchart of patient recruitment, and the characteristics of these groups can be consulted in references13, 14.

Fresh saliva samples were received at CIBIR, aliquoted in tubes (approximately 500 μl) and stored at − 80 °C for further analysis. After defrosting, lysozyme concentration in saliva was measured through enzyme-linked immunosorbent assay (ELISA) from ABCAM (Ab108880) according to the manufacturer’s instructions. The intra-assay and the inter-assay coefficient of variation (CV) expressed as percentage are 4.7% and 9.5%, respectively.

After defrosting, DNA was extracted using a modified version of the DNeasy® Blood and Tissue kit (Qiagen, Venlo, The Netherlands). In this protocol, 250 µl of saliva is transferred to an Eppendorf tube. Then, 1 ml of PBS is added and the mixture is centrifuged at 1800 g for 5 min. The supernatant is decanted and the pellet is resuspended in 180 µl of PBS and transferred to another Eppendorf tube. Next, 25 µl of Proteinase K and 200 µl of Buffer AL are added, followed by vortexing for 20 s and centrifugation. The mixture is incubated at 56 °C for 10 min. After that, 200 µl of ethanol is added, vortexed for 20 s, and centrifuged. The resulting solution is transferred to a column and centrifuged, at 6000 g for 1 min and the collected tube with the filtrate is discarded. The column is then placed in a new collection tube, 500 µl of Buffer AW1 is added and centrifuged, and the collected tube with the filtrate is discarded. Then, the column is again placed in a new collection tube, 500 µl of Buffer AW2 is added and centrifuged at 20,000 g for 3 min, and the collected tube with the filtrate is discarded. Finally, the column is placed in a new Eppendorf tube, 50 µl of Buffer AE is added, and incubated at room temperature for 1 min. The mixture is then centrifuged at 6000 g for 1 min and the elution process is repeated. After this process, purity, concentration, and quality were determined by a Qubit 3.0 fluorometer (kit dsDNA HS, Thermo Fisher Scientific, MA, USA) and a Fragment Analyzer (HS Genomic DNA 50 Kb kit, Agilent, USA). Then, samples were amplified for the 16S rRNA hypervariable regions V3-V416 and sequencing was performed using an Illumina sequencer (MiSeq, 2 × 300 pb, paired end) at the Genomics & Bioinformatics Core Facility at CIBIR.

The first step in computational analysis was to check the quality of reads by the quality control tool FastQC program17. Then, the Qiime2 pipeline18 was used along the bioinformatic analysis. Firstly, the raw sequences already demultiplexed (mapping the barcodes to the samples they belong) by the Illumina sequencer were denoised using the DADA2 software19. Specifically, it was performed the trimming of sequencing adapters and primer regions, the filtering of noisy reads, the dereplicate of our sequences to reduce repetition, the joint of paired reads, the identification of amplicon sequence variants (ASVs), and the elimination of chimeras. Secondly, we used the SILVA database20 trained with the V3-V4 amplification primers using during the wet-lab process to do the taxonomic assignation. The alpha and beta diversity were analysed: α-diversity is a measure of sample-level species richness, whereas β-diversity describes inter-subject similarity of microbial composition and facilitates the identification of broad differences between samples. The measure of α-diversity was analysed using Observed Features, Chao1 index, Fisher’s alpha, Pielou’s index, Shannon index, and Simpson index and differences between groups were assessed through Kruskal–Wallis test. Observed features, Chao1 index and Fisher’s alpha are based in richness, Pielou’s index is based in evenness and Shannon index and Simpson index are based in diversity (richness + evenness). The measure of β-diversity was analysed using PERMANOVA (999 permutations) on Bray Curtis metric and visualized using Principal Coordinate Analysis (PCoA) by R software (version 4.0.5) and R Studio (version 1.4.1105). Finally, the analysis of the differential composition of bacteria was carried out with the ANCOM methodology at phylum, order, and genus taxonomic levels. This methodology accounts for the underlying structure in the data and is widely used for comparing the composition of microbiomes in two or more populations, with no assumptions of population distribution21. ANCOM just runs a bunch of pairwise tests for the following sub-hypothesis:

where \({x}_{i}\) denotates the ith ASV abundance from sample x, \({x}_{j}\) denotates de jth ASV abundance from sample x, \({y}_{i}\) denotates the ith ASV abundance from sample y and \({y}_{j}\) denotates de jth ASV abundance from sample y. The W value is just a count of the number of times that \({H}_{0(ij)}\) is rejected for the ith ASV.

Characteristics of the population and lysozyme levels are presented as mean ± standard error of the man (SEM). Categorical variables were analysed using the Chi-square or Fisher’s exact test. Normal distribution of quantitative variables was checked using the Shapiro–Wilk test. Comparison between two groups were performed using unpaired t test or Mann–Whitney U test depending on the normality of the data. Comparison between three or more groups were analysed using ANOVA followed by Tukey post-hoc regardless the normality of the data. p values < 0.05 and false discovery rates (FDRs) < 0.05 were considered as statistically significant. Statistical analysis was performed using R software (version 4.0.5) and R Studio (version 1.4.1105).

This study was performed following the Helsinki Declaration and was approved by the Committee for Ethics in Drug Research in La Rioja (CEImLAR) (28 February 2019, reference number 349). All participants provided their written informed consent.

The main characteristics of the recruited population (Supplementary Table 1 and Supplementary Table 2) were previously described13, 14. Of note, smoking habits were higher in the naive group and INSTIs-treated group compared to controls (11.54 vs. 46.67% and 66.67%, respectively). On the other hand, none of the individuals recruited reported periodontal diseases. Finally, none of the patients reported differential dietary patterns.

A statistically significant increase in salivary lysozyme concentrations was observed in naive HIV-infected patients compared to the control population (p < 0.001) and also compared to INSTIs-treated patients (p < 0.05) (Fig. 1). No statistically significant difference was detected between the control and the INSTIs-treated group (Fig. 1).

Levels of lysozyme in the studied population comparing control group, naive group, and INSTIs-treated group. ***p < 0.001 versus control and #p < 0.05 versus INSTIs. INSTIs integrase strand transfer inhibitors-based treatment.

After DNA extraction, all samples had a DNA concentration between 11.40 and 336.00 ng/μl with a mean of 86.77 ng/μl and a median of 73.70 ng/μl. A total of 12.5 ng of each one was used for the amplification PCR, the first step at the library preparation protocol. Besides, attending to the DNA quality, measured as DNA integrity through the GQN parameter in the FA system, all samples were ranked between 0.6 and 3.6 values. In samples where a specific fragment of 450–500 bp was amplified, the assessment of genomic integrity may not be as significant since the focus is on the amplification of a specific target region rather than the overall integrity of the genome. As an intermediate step, after the amplicon PCR, but before the indexing step, the presence of the amplicon was verified using Fragment Analyzer (Amplicon Library kit, 35–5000 bp, AATI), resulting in fragments ranging from 500 to 540 bp in size. This step confirmed the successful amplification of the target region and ensured the presence of the desired amplicon for further analysis or sequencing. After library preparation, all the final libraries were quantified once again using Qubit fluorometry and qualified using Fragment Analyzer (Amplicon Library kit, 35–5000 bp, AATI) to determine the average amplicon size. This information (concentration and medium library size) allowed us to subsequently normalize the libraries to the same concentration (nM) in order to prepare an equimolar pool at 10 pM. In this case, the libraries had a concentration ranging from 31 to 210 ng/uL, with a mean value of 94.80 ng/uL, as determined by fluorometry. Additionally, the average amplicon size ranged from 600 to 630 bp, taking into account that this size includes 120 bp corresponding to the sequencing adapters included in the library preparation for the V3-V4 amplicon by PCR. The individual average amplicon size with the library concentration (ng/uL) was used to calculate the molarity of each sample, normalized at 10 nM and a final equimolar pool at 10 pM, that was introduced into the MiSeq sequencer.

The raw sequencing reads obtained from the saliva samples in this study ranged from a minimum of 134,247 reads to a maximum of 396,314 reads. After adapter trimming, removal of low-quality regions, and processing using the DADA2 protocol, a total of between 66,634 and 187,188 reads (corresponding to a percentage range relative to the initial 38.33–50.34% of the initial reads) were maintained.

No statistical differences were observed on α-diversity of oral bacteriome (Observed features, Chao1 index, Fisher’s alpha, Pielou’s evenness, Simpson index, and Shannon index) (Supplementary Fig. 2).

The analysis of β-diversity did not either reveal a different clustering pattern between the three groups analysed (Supplementary Fig. 3).

Concerning oral microbiota composition and/or abundance, a total of 20 phyla and 68 orders were detected. The three most abundant phyla in saliva were Bacteroidota (12.61–64.16%), Bacillota (7.37–66.14%), and Pseudomonadota (0.12–54.40%). In the order level, Bacteroidales (10.05–63.13%), Lactobacillales (2.65–56.11%), and Betaproteobacteriales (0.04–42.54%) were the most abundant.

When comparing controls to naive patients, a significant higher abundance in the phylum Pseudomonadota, in the order Acholeplasmatales (phylum Bacillota), and in the genera Ezakiella (order Tissierellia, phylum Bacillota) and Acholeplasma (order Acholeplasmatales, phylum Bacillota) were observed in the naive group compared to controls (Table 1). On the other hand, when controls were compared against INSTIs-treated HIV-infected patients, a higher abundance in the phylum Mycoplasmatota, in the order Acholeplasmatales (phylum Bacillota), and in the genera Acholeplasma (order Acholeplasmatales, phylum Bacillota) and uncultured Eubacteriaceae bacterium (order Eubacteriales, phylum Bacillota) were detected in the INSTIs-treated HIV-infected patients (Table 1). No differences were observed between both HIV-infected groups. No lower abundances were observed in HIV-infected patients.

As smoking has an impact on oral microbiota22, we assessed if smoking habits could modulate those actions exerted by HIV infection per se and/or INSTIs. We revealed that oral microbiota coming from smoking individuals (irrespectively of their HIV status) showed a statistically significant increase in all α-diversity indexes analysed (p < 0.01 for Observed features, Chao1 index, and Fisher’s alpha and p < 0.05 Pielou’s evenness, Shannon index, and Simpson index) compared to non-smokers (Supplementary Fig. 4). Besides, the analysis of β-diversity revealed a different clustering between the two groups (p < 0.001) (Supplementary Fig. 5) suggesting a clear effect of tobacco on oral microbiota composition/community.

Similarly, when only HIV-infected patients were split into two groups according to their smoking habits (and regardless ART), those samples from smoking individuals presented a statistically significant increase in the richness indexes analysed (p < 0.05 for Observed features, Chao1 index, and Fisher’s alpha) (Fig. 2) compared to non-smoking HIV-infected population. Moreover, the analysis of β-diversity revealed a different clustering between the two groups (p < 0.05) (Fig. 3).

Different indexes of α-diversity from bacteria in salivary samples of HIV-infected patients regarding smoking habit. *p < 0.05 versus non-smokers.

PCoAs from bacteria in salivary samples from HIV-infected patients regarding smoking habit (accounting for 30% of the total variation [Component 1 = 18% and Component 2 = 12%]). Results are plotted according to the first two principal components. Each circle represents a sample: red circles represent the non-smoking individuals and blue circles represent the smoking individuals. The clustering of samples is represented by their respective 95% confidence interval ellipse. *p < 0.05 smokers versus non-smokers.

Considering the clear impact of smoking on oral microbiota, we next analysed the controls, naive HIV-infected patients, and INSTIs-treated patients taking into account their smoking status, comparing firstly the non-smokers and, secondly, the smokers.

No differences were detected in the α-diversity of the non-smoking population (Supplementary Fig. 6). Besides, the analysis of β-diversity did not either reveal a different clustering between the three groups (Supplementary Fig. 7). However, the analysis of differential abundance revealed some taxonomical orders differentially represented between the groups, although the magnitudes of such differences were not very potent. Specifically, the INSTIs-treated group showed a higher abundance of phylum Mycoplasmatota (W = 1) and a lower abundance of phyla Epsilonbacteraeota (W = 4), Bacteroidota (W = 1), Spirochaetota (W = 1), Cemmatimonadetes (W = 1), Chlamydiota (W = 1), Acidobacteriota (W = 1), and Campylobacterota (W = 1) and of order Campylobacterales (W = 10) when compared to controls. On the other hand, a higher abundance of phyla Epsilonbacteraeota (W = 2) and Bacteroidota (W = 1) were observed in the INSTIs-treated group compared to naive group.

On the other hand, no significant differences were observed in either α-diversity or β-diversity or differential abundance when only smokers were analysed (Supplementary Fig. 8 and 9).

Contrary to the results observed in our previous works focused on gut microbiota13, 14, the effects of HIV infection and INSTIs on oral microbiota seem to be not very significant, probably due to the modulation of other factors such as smoking, and the greatest outward exposure of the oral cavity.

Lysozyme is an antimicrobial protein, expressed by various cells including neutrophils, macrophages, and epithelial cells. It is abundant in saliva and plays an important role in the host constitutive defense system23. Previous studies have revealed an increase in salivary lysozyme levels in immunocompromised patients24, 25, which could serve as a compensatory mechanism to alleviate the decrease in the number of CD4 T lymphocytes. Besides, it was revealed that lysozyme from chicken egg white, human milk, and human neutrophils possesses antiviral activity against HIV-infection in vitro26. Considering that our study was focused on oral microbiota, we believed interesting to quantify this marker as an indirect measure of the HIV-infection and more specifically, the mechanisms that are set in motion to try to counteract the infection action and the potential association with oral microbiota. In fact, we observed a statistically significant increase of salivary lysozyme concentration in naive HIV-infected patients compared to controls and INSTIs in parallel with the decrease observed in CD4 T lymphocytes in these patients, suggesting a possible relationship between lysozyme, infection, and immunological response. Interestingly, INSTI-treated patients showed a significant decrease in the salivary levels of this enzyme along with the concomitant increase of CD4 T lymphocytes. To our knowledge, this is the first study that explored the effects of INSTIs on salivary lysozyme concentrations. Our results could highlight the potential usage of this enzyme as a non-invasive biomarker of the immunological recovery observed after INSTIs, but more studies are needed in this regard.

Our study was unable to detect differences in α nor in β-diversity associated with HIV infection and/or INSTIs-based treatment. Similarly, Presti et al.6, were unable to detect differences in α nor in β-diversity between the naive group and the ART-treated group (treated with EFV/FTC/tenofovir(TDF)). Imahashi et al.27, did not either show differences in species richness between HIV negative people and ART-treated HIV-infected subjects. However, Li et al.28, revealed a statistically significant decrease in Chao1 index and Shannon index between the naive group and the control group and also a different clustering between both groups. The differences could be explained by the fact that the populations among studies were different. In fact, in the study of Li et al.28, only HIV-infected MSM patients from Beijing were included, whereas in our study the percentage of MSM were 60.00% in the naive group and only 46.67% in the ART-treated group. In addition, only 15 HIV-infected patients in our study were not Caucasians (none of them from Asia). Thus, the potential influence of the origin of subjects (in terms of “dietary habits”) and/or the sexual orientation and mode of transmission of HIV infection in the composition of oral microbiota needs to be further investigated. In addition, the hypervariable regions of the 16S rRNA gene analysed were also different between both studies (V4-V5 vs. V3-V4 in our study).

In our study, the three most abundant phyla were Bacteroidetes (now known as Bacteroidota), Firmicutes (now known as Bacillota), and Pseudomonadota, which correlates with the results from Li et al.28 and Presti et al.6. We observed some taxa increased in the naive group and in the INSTIs-treated group compared to controls. It is worth mentioning the increase observed in the genus Acholeplasma in HIV-infected patients irrespectively of ART. It has been demonstrated that Acholeplasma laidlawii could facilitate the entry of HIV infection into cells through the bind to its glycoglycerolipids29 and treatment is not able to reduce its presence. The clinical impact (systemic and locally) of these higher abundances should be deeply evaluated.

We also observed that smoking influences salivary bacteriome diversity and composition, as previously revealed by Gopinath et al.15, in the general population. However, and despite this fact, we also observed that the impact of HIV infection and INSTIs on oral microbiota were partially independent of smoking since no effects were observed in either α-diversity or β-diversity in both smokers and non-smokers separately analysed. Of note, some differences (not very statistically potent according to the W value) were detected in the abundance of some bacterial taxonomical orders in the non-smoking population whereas no modifications were observed in the smokers suggesting that smoking could modulate or even “buffer/mask” the effects of HIV and INSTIs on oral microbiota. In addition, the less potent effects observed on oral microbiota could be explained by the fact that the oral cavity presents a higher interindividual variability as it is more exposed to environment/diet. Besides, ART pills do not dissolve in saliva, so they have more a systemic effect than a local effect in oral cavity.

This study has some limitations. Some of the patients were recruited during COVID-19 pandemic. However, none of them reported symptoms related to COVID-19 before sample collection and they were not vaccinated in the previous month. Besides, the number of patients in each group is small (what increases the probability of type-II error), but large enough to detect differences among groups. Moreover, some differences between the controls and HIV-infected groups were detected, such as gender, age, and smoking habits, all of them factors that could have an impact on microbiota composition. However, the two HIV-infected groups were well-balanced in terms of these factors. Besides, saliva sampling, although being practical and non-invasive, may not capture the complete diversity of anaerobic bacteria in the mouth, so it cannot be ruled out that other methods will give different results. Finally, patients fulfilled a survey focused on dietary habits/patterns to detect those habits that could have an impact on microbiota, such as veganism or excessive consumptions of prebiotics and/or probiotics. None of the patients reported such differential habits.

Lysozyme concentration is elevated in naive patients compared to control population suggesting a compensatory mechanism to alleviate the decrease in the number of CD4 T lymphocytes. In fact, INSTI-treatment increased CD4 T cells in the patients and, at the same time, diminished the levels of lysozyme corroborating the association among lysozyme-inmune response and infection. On the other hand, HIV infection and INSTIs impact on oral microbiota was not very potent, probably due to the modulation of other factors such as smoking and the greatest outward exposure of the oral cavity. In fact, smoking increases salivary bacteriome diversity and it also influences bacteriome composition. Therefore, smoking could modulate or even “buffer/mask” the effects of HIV and INSTIs on oral microbiota and this could be the reason for the mild effects observed on salivary microbiome.

The datasets generated during and/or analysed during the current study are available in the NCBI SRA repository, http://www.ncbi.nlm.nih.gov/bioproject/819232, and also available from corresponding author on reasonable request.

Pérez-Matute, P., Íñiguez, M., Villanueva-Millán, M. J. & Oteo, J. A. Chapter 32—The oral, genital and gut microbiome in HIV infection. Microbiome metabolome diagnosis 307–323 (Academic Press, 2019).

Google Scholar

Li, Y. et al. HIV infection and microbial diversity in saliva. J. Clin. Microbiol. Am. Soc. Microbiol. 52, 1400–1411 (2014).

Article Google Scholar

Nasi, M., Pinti, M., Mussini, C. & Cossarizza, A. Persistent inflammation in HIV infection: Established concepts, new perspectives. Immunol. Lett. 161, 184–188 (2014).

Article CAS PubMed Google Scholar

Mitchell, D. et al. HIV nucleoside reverse transcriptase inhibitors efavirenz and tenofovir change the growth and differentiation of primary gingival epithelium. HIV Med. 15, 196–202 (2014).

Article CAS PubMed Google Scholar

Ghosh, S. K. et al. Comparison of epigenetic profiles of human oral epithelial cells from HIV-positive (on HAART) and HIV-negative subjects. Epigenetics 8, 703–709 (2013).

Article CAS PubMed PubMed Central Google Scholar

Presti, R. M. et al. Alterations in the oral microbiome in HIV-infected participants after antiretroviral therapy administration are influenced by immune status. AIDS 32, 1279–1287 (2018).

Article CAS PubMed Google Scholar

Schafer, C. A. et al. Saliva diagnostics: Utilizing oral fluids to determine health status. Monogr. Oral Sci. 24, 88–98 (2014).

Article PubMed Google Scholar

Yakob, M., Fuentes, L., Wang, M. B., Abemayor, E. & Wong, D. T. W. Salivary biomarkers for detection of oral squamous cell carcinoma—Current state and recent advances. Curr. Oral Heal. Rep. 1, 133–141 (2014).

Article Google Scholar

Yoshizawa, J. M. et al. Salivary biomarkers: Toward future clinical and diagnostic utilities. Clin. Microbiol. Rev. 26, 781–791 (2013).

Article CAS PubMed PubMed Central Google Scholar

Giannobile, W. V., McDevitt, J. T., Niedbala, R. S. & Malamud, D. Translational and clinical applications of salivary diagnostics. Adv. Dent. Res. 23, 375–380 (2011).

Article CAS PubMed PubMed Central Google Scholar

Miller, C. S. et al. Current developments in salivary diagnostics. Biomark. Med. 4, 171–189 (2010).

Article PubMed Google Scholar

Villanueva-Millán, M. J., Pérez-Matute, P., Recio-Fernández, E., Rosales, J. M. L. & Oteo, J. A. Differential effects of antiretrovirals on microbial translocation and gut microbiota composition of HIV-infected patients. J. Int. AIDS Soc. 20, 21526 (2017).

Article PubMed PubMed Central Google Scholar

Villoslada-Blanco, P. et al. Integrase inhibitors partially restore bacterial translocation, inflammation and gut permeability induced by HIV infection: Impact on gut microbiota. Infect. Dis. Ther. 11, 1541–1557 (2022).

Article PubMed PubMed Central Google Scholar

Villoslada-Blanco, P. et al. Impact of HIV infection and integrase strand transfer inhibitors-based treatment on the gut virome. Sci. Rep. 12, 21658 (2022).

Article ADS PubMed PubMed Central Google Scholar

Gopinath, D. et al. Compositional profile of mucosal bacteriome of smokers and smokeless tobacco users. Clin. Oral. Investig. 26, 1647–1656 (2022).

Article PubMed Google Scholar

16S Metagenomic Sequencing Library Preparation. https://support.illumina.com/documents/documentation/chemistry_documentation/16s/16s-metagenomic-library-prep-guide-15044223-b.pdf.

FastQC. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

Qiime2. Qiime2. https://qiime2.org/.

Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. https://doi.org/10.1038/nmeth.3869 (2016).

Article CAS PubMed PubMed Central Google Scholar

SILVA. https://www.arb-silva.de/.

Mandal, S. et al. Analysis of composition of microbiomes: A novel method for studying microbial composition. Microb. Ecol. Health Dis. 26, 27663 (2015).

PubMed Google Scholar

Ambrose, J. A. & Barua, R. S. The pathophysiology of cigarette smoking and cardiovascular disease: An update. J. Am. Coll. Cardiol. 43, 1731–1737 (2004).

Article CAS PubMed Google Scholar

Wiesner, J. & Vilcinskas, A. Antimicrobial peptides: The ancient arm of the human immune system. Virulence 1, 440–464 (2010).

Article PubMed Google Scholar

Arnold, R. R. et al. Secretory immunity and immunodeficiency. Adv. Exp. Med. Biol. 107, 401–410 (1978).

Article CAS PubMed Google Scholar

Schachtsiek, M., Hammes, W. P. & Hertel, C. Characterization of Lactobacillus coryniformis DSM 20001T surface protein Cpf mediating coaggregation with and aggregation among pathogens. Appl. Environ. Microbiol. 70, 7078–7085 (2004).

Article ADS CAS PubMed PubMed Central Google Scholar

Lee-Huang, S. et al. Lysozyme and RNases as anti-HIV components in beta-core preparations of human chorionic gonadotropin. Proc. Natl. Acad. Sci. U. S. A. 96, 2678–81 (1999).

Article ADS CAS PubMed PubMed Central Google Scholar

Imahashi, M. et al. Impact of long-term antiretroviral therapy on gut and oral microbiotas in HIV-1-infected patients. Sci. Rep. 11, 1–10 (2021).

Article Google Scholar

Li, S. et al. 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 (2021).

Article CAS PubMed PubMed Central Google Scholar

Shimizu, T. et al. Glycoglycerolipid from the membranes of Acholeplasma laidlawii binds to human immunodeficiency virus-1 (HIV-1) and accelerates its entry into cells. Curr. Microbiol. 48, 182–188 (2004).

Article CAS PubMed Google Scholar

Download references

We would like to thank all participants of this study and physicians involved in patient recruitment.

This work was supported by Fundación Rioja Salud and PVB was granted a predoctoral grant from Consejería de Desarrollo Económico e Innovación (Government of La Rioja).

Infectious Diseases, Microbiota and Metabolism Unit, Infectious Diseases Department, Center for Biomedical Research of La Rioja (CIBIR), C/Piqueras 98, CIBIR Building, Third Floor, 26006, Logroño, La Rioja, Spain

Pablo Villoslada-Blanco, Patricia Pérez-Matute, Emma Recio-Fernández, María Íñiguez & José A. Oteo

Centro de Salud Siete Infantes de Lara, Logroño, La Rioja, Spain

Pilar Blanco-Navarrete

Infectious Diseases Department, Hospital Universitario San Pedro, Logroño, La Rioja, Spain

Luis Metola, Valvanera Ibarra, Jorge Alba & José A. Oteo

Genomics and Bioinformatics Platform, Center for Biomedical Research of La Rioja (CIBIR), Logroño, La Rioja, Spain

María de Toro

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

P.P.M. and J.A.O. conceived the study concept. P.B.N., L.M., V.I., J.A. and J.A.O. recruited the patients. P.V.B., M.I., and E.R. conducted the prosecution of the samples and the DNA extraction from saliva. E.R. did the ELISAs to measure salivary lysozyme levels. P.V.B. and M.T. carried out the 16S rRNA gene sequencing and the bioinformatic analysis. P.V.B. and P.P.M. did the statistical analysis. P.V.B., P.P.M. and J.A.O. wrote the manuscript. All authors read an approved the final manuscript.

Correspondence to Patricia Pérez-Matute.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Villoslada-Blanco, P., Pérez-Matute, P., Recio-Fernández, E. et al. Beyond the effects of HIV infection and integrase inhibitors-based therapies on oral bacteriome. Sci Rep 13, 14327 (2023). https://doi.org/10.1038/s41598-023-41434-5

Download citation

Received: 20 March 2023

Accepted: 26 August 2023

Published: 31 August 2023

DOI: https://doi.org/10.1038/s41598-023-41434-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.