Cecilia Rosel-Pech1,2†
Sandra Pinto-Cardoso3*†
Monserrat Chávez-Torres3
Nadia Montufar3
Iván Osuna-Padilla3
Santiago Ávila-Ríos3
Gustavo Reyes-Terán1
Charmina Aguirre-Alvarado2
Norma Angelica Matías Juan4Héctor Pérez-Lorenzana5
José Guillermo Vázquez-Rosales6
Vilma Carolina Bekker-Méndez2*- 1Posgrado en Ciencias Biológicas, Unidad de Posgrado, Universidad Nacional Autónoma de México (UNAM), Ciudad de México, Mexico
- 2Unidad de Investigación Médica en Inmunología e Infectología, Hospital de Infectología “Dr. Daniel Méndez Hernández”, Centro Médico Nacional “La Raza”, Instituto Mexicano del Seguro Social (IMSS), Ciudad de México, Mexico
- 3Centro de Investigación en Enfermedades Infecciosas, Instituto Nacional de Enfermedades Respiratorias Ismael Cosío Villegas, Ciudad de México, Mexico
- 4Hospital de Infectología “Dr. Daniel Méndez Hernández”, Centro Médico Nacional “La Raza”, Instituto Mexicano del Seguro Social (IMSS), Ciudad de México, Mexico
- 5UMAE Hospital General Dr. Gaudencio González Garza, Centro Médico Nacional La Raza, Instituto Mexicano del Seguro Social (IMSS), Ciudad de México, Mexico
- 6Hospital de Pediatría “Doctor Silvestre Frenk Freund”, Centro Médico Nacional Siglo XXI, Instituto Mexicano del Seguro Social, México, Mexico
Introduction: Our understanding of HIV-associated gut microbial dysbiosis in children perinatally-infected with HIV (CLWH) lags behind that of adults living with HIV. Childhood represents a critical window for the gut microbiota. Any disturbances, including prolonged exposure to HIV, antiretroviral drugs, and antibiotics are likely to have a significant impact on long-term health, resulting in a less resilient gut microbiome. The objective of our study was to characterize the gut microbiota in CLWH, and compare it with HIV-unexposed and -uninfected children.
Methods: We enrolled 31 children aged 3 to 15 years; 15 were CLWH and 16 were HUU. We assessed dietary patterns and quality; quantified soluble and cellular markers of HIV disease progression by flow cytometry, enzyme-linked immunosorbent and multiplex-bead assays, and profiled the gut microbiota by 16S rRNA sequencing. We explored relationships between the gut microbiota, antibiotic exposure, dietary habits, soluble and cellular markers and host metadata.
Results: Children had a Western-type diet, their median health eating index score was 67.06 (interquartile range 58.76-74.66). We found no discernable impact of HIV on the gut microbiota. Alpha diversity metrics did not differ between CLWH and HUU. Sex impacted the gut microbiota (R-squared= 0.052, PERMANOVA p=0.024). Male children had higher microbial richness compared with female children. Two taxa were found to discriminate female from male children independently from HIV status: Firmicutes for males, and Bacteroides for females. Markers of HIV disease progression were comparable between CLWH and HUU, except for the frequency of exhausted CD4+ T cells (PD-1+) which was increased in CLWH (p=0.0024 after adjusting for confounders). Both the frequency of exhausted CD4+ and activated CD4+ T cells (CD38+ HLADR+) correlated positively with the relative abundance of Proteobacteria (rho=0.568. false discovery rate (FDR)-adjusted p= 0.029, and rho=0.62, FDR-adjusted p=0.0126, respectively).
Conclusion: The gut microbiota of CLWH appears similar to that of HUU, and most markers of HIV disease progression are normalized with long-term ART, suggesting a beneficial effect of the latter on the gut microbial ecology. The relationship between exhausted and activated CD4+ T cells and Proteobacteria suggests a connection between the gut microbiome, and premature aging in CLWH.
1 Introduction
In July 2022, the World Health Organization (WHO) reported the estimated global human immunodeficiency virus (HIV) prevalence across all ages; there are 38.4 million (M) people living with HIV (PLWH) (1). Of those, 1.7M are children aged 15 and less. In 2021, there was a 32% and 52% global decrease in new infections and HIV-related deaths relative to 2010, respectively. Still, 650,000 PLWH died of HIV-related complications, of these 98,000 were children. This decline in HIV-related deaths over time is largely due to the increased availability and effectiveness of antiretroviral therapy (ART) and continuous efforts towards reaching the 95-95-95 targets of the HIV service cascade. Globally, 75% (28.7M) of PWLH received ART in 2021. In Mexico, the estimated number of PLWH is 227,378. The cumulative number of perinatally-infected children aged 0 to 14 years is 4,229 (2). In 2023, 14 children were perinatally-infected with HIV [CLWH, (2)]. The estimated ART coverage among adults and children is 61% and 49% respectively, 54% are virally-suppressed, and the mother-to child transmission is 0.3%. Although the gap in life expectancy is gradually decreasing, PLWH (adults and children alike) are still at a higher risk of noncommunicable diseases (NCD) compared to the general population (3). This is partly due to persistent immune dysfunction and chronic inflammation which fail to normalize despite effective suppressive ART (4, 5).
In 2006, a seminal paper published in Nature Medicine concluded that microbial translocation, the passage of microbial products and/or bacteria from the gut to the systemic circulation, is a key driver of chronic inflammation and HIV pathogenesis, linking the gut microbiome to HIV disease progression (6, 7) and NCD (5, 8–10). This study (and others) conceptualized the current paradigm that functional (massive depletion of gut T helper 17 cells) and structural (loss of tight junctions) events secondary to HIV infection cause lasting damage to the gut mucosa, providing favorable conditions for alterations in the composition and function of gut microbial communities to occur (5, 11–14). These alterations are commonly referred to as microbial dysbiosis. Collectively, studies in adults PLWH (ALWH) have shown a decrease in bacterial (12) and fungal diversity (15, 16), an expansion of the virome (9, 17), an increase in pathobionts and a decrease in anti-inflammatory commensals (14), compared with seronegative individuals. Microbial dysbiosis is partially restored by ART. Also, studies have shown that ART is linked to a distinct microbial dysbiosis independent from HIV (18–21). Less is known in CLWH, as only a handful of studies have been published (22–24). Contrary to ALWH, whose studies have been published in high income countries (for the most part), studies in CLWH have been conducted in Africa [Cameron (24), Zimbabwe (25), South Africa (26)], India (22), and Italy (27). Collectively, studies suggest similarities between ALWH and CLWH, despite differences in mode of HIV transmission (vertical versus horizontal), absence of some confounding variables (sexual preference) and the presence of confounders specific to CLWH [use of co-trimoxazole, a prophylactic antibiotic recommended by the WHO in countries with high burden of malaria and bacterial infections (28)].
Understanding of HIV-associated gut microbial dysbiosis in CLWH lags behind that of ALWH. Since childhood represents a critical window for the gut microbiota (development and maturation), and the immune system (education), disturbances are likely to have a significant impact on long-term health, quality of life and longevity (29–31). Also, prolonged exposure to HIV, ART, antibiotics (32) and maternal factors (33), may result in a less resilient gut microbiome and a compromised gut microbiome-immune system axis, increasing tolerance to pathobionts. This, in turn, may lead to increased prevalence of cardiometabolic alterations, and other NCD (and their premature clinical presentation compared to the general population). The objective of our study was to characterize the gut microbiota in CLWH in Mexico, and compare it with HIV-unexposed and -uninfected children (HUU). We used dietary questionnaires and healthy eating index (HEI) scores to assess dietary patterns and quality; flow cytometry to quantify CD4+ and CD8+ T cell immune activation, senescence, exhaustion and cell cycling; enzyme-linked immunosorbent assays to quantify markers of inflammation, bacterial and fungal translocation, enterocyte integrity, and adipose tissue inflammation, multiplex-bead assay to quantify soluble markers of immune activation, and profiled the gut microbiota using 16S rRNA sequencing using fecal samples. We also explored relationships between the gut microbiota, antibiotic exposure, dietary habits, soluble and cellular markers of HIV disease progression and host metadata.
2 Methods
2.1 Ethics statement
This study was conducted according to the Declaration of Helsinki and its later amendments. Ethical approval for this study was obtained from the Instituto Nacional de Enfermedades Respiratorias (INER, approval number C23-20) and the Instituto Mexicano del Seguro Social (IMSS, approval number R-2018-785-130). Written informed consent and age-appropriate assent were obtained from all participants and/or parents before study enrollment.
2.2 Study participants
We enrolled a total 31 children aged 3–15 years between 2019 and 2020; 15 were perinatally infected with HIV; here referred to as children living with HIV (CLWH) and 16 were HIV-unexposed and uninfected children (HUU). All children were enrolled from three pediatric services within IMSS. CLWH were enrolled from Servicios de Infectología Pediátrica, Centro Médico Nacional (CMN) La Raza and Servicios de Infectología Pediátrica, CMN, Siglo XXI. HUU were enrolled from Servicio de Cirugía Pediátrica, Hospital General, CMN “La Raza”.
2.3 Research samples PBMC, plasma and stools
Whole blood (8mL) was collected by a phlebotomist by venipuncture in Vacutainer® EDTA tubes (ethylenediaminetetraacetic acid, BD, New Jersey, United States). Blood samples were immediately processed to separate plasma by centrifugation (3,000 rpm for 10 min) and cryopreserved at −80°C. Next, peripheral blood mononuclear cells (PBMC) were isolated by Ficoll density gradient centrifugation (STEMCELL technologies, Vancouver, Canada) and stored at -196°C. Stool samples were collected by each participant and processed immediately. First, stools were classified according to the Bristol Stool Chart (34), then 220 mgs of feces were aliquoted in 2mL Eppendorf LoBind tubes, and stored at -80°C.
2.4 Nutritional assessment
Anthropometric and dietary assessment were performed by qualified dietitians. First, height and weight were measured using a Seca mechanical scale with a stadiometer (Seca, Connecticut, United States). Body mass index (BMI) and BMI-for-age (BMI/A) were calculated in accordance to the WHO´s reference values using the WHO AnthroPlus (35). Dietary assessment was conducted using semi-quantitative methods including a food frequency questionnaire (FFQ) and three 24-hour dietary recalls validated for the Mexican population as previously described (36). Estimated nutrient and energy intake were determined using the NutriSurvey software (37) and total energy consumed was expressed in kilocalories (kcal), macronutrients, and micronutrients in grams (g) and milligrams (mg), respectively. Diet quality was assessed using the 2010 Healthy Eating Index [HEI, (38)] and information derived from three 24-hour dietary recalls. HEI is a score system that recapitulates how well individuals adhere to dietary recommendations for Americans. This score system includes two types of components distributed across 12 groups (total fruit, whole fruit, total vegetables, greens and beans, whole grains, dairy, total protein foods, seafood, plant proteins, fatty acids, refined grains, sodium, and empty calories) and ranges from 0 to 100, 100 denoting excellent compliance to dietary recommendations. The HEI score was further categorized into healthy (HEI score above 80), needs changes (HEI score between 50-79) and unhealthy (HEI score below 50).
2.5 Flow cytometry
Canonical markers of immune activation (CD38 and HLADR), senescence (CD57), exhaustion (programmed cell death 1, PD-1) and cell cycling (Ki67) were quantified in both CD4+ and CD8+ T lymphocytes using flow cytometry as previously published (39) with minor modifications. Monoclonal antibodies used in this study are listed in Supplementary Table 1. The LIVE/Dead™ fixable blue dead stain was used to exclude dead cells (viability marker, ThermoFisher Scientific). A dump gate was added to exclude unwanted populations. Briefly, PBMCs were thawed and rested for at least 30 min. Cells were washed in 15mL of phosphate buffered solution (PBS), centrifuged at 500 x g for 5 min at room temperature (RT). After discarding the supernatant, cells were stained with a cocktail of monoclonal extracellular antibodies for 20 min in the dark. Cells were washed twice in PBS and permeabilized with Fix/Perm kit (ThermoFisher Scientific) for intracellular staining (Ki-67) and incubated for 60 min. Cells were washed twice in PBS, and fixed in 300µL of 1% paraformaldehyde (Sigma-Aldrich, St. Louis, Missouri, USA) and kept in the dark at 4°C until acquisition. Flow cytometry acquisition was performed on a BD LSRFortessa (BD Biosciences, San Jose, CA, USA) within 24 h of staining. A minimum of 3 million events were acquired per sample. Quality controls were performed using BD cytometer Setup and Tracking Beads and rainbow beads (BD biosciences). A compensation matrix was calculated and applied using BD Comp Beads (BD Biosciences). Data was analyzed using FlowJo™ v10.9 (BD, Ashland, Oregon, Unites States). Raw FCS (flow cytometry standard) files were first quality-controlled using FlowAI v2.3.1 (40) with default parameters, except for the changepoint penalty which was set to 500. A mean of 2,529,571 events were retained after quality-control. Populations were pre-gated for singlets and morphology by using forward versus side scatter, followed by the exclusion of dead cells (leaving only aqua dye negative events). Live cells were gated on CD3+ and then on CD4+ or CD8+ T-cells. A representative gating strategy is shown in Supplementary Figure 1. FMO (fluorescence minus one) were used to accurately set the gates for positive populations.
2.6 Quantification of markers by ELISA
Markers of microbial and fungal translocation (soluble CD14, sCD14, 1,3-β-D-Glucan), epithelial barrier integrity (intestinal fatty binding acid protein, I-FABP), inflammation (C reactive protein, CRP), and adipose tissue inflammation (sCD163, monocyte chemoattractant protein-1, MCP-1 and adiponectin) were measured in plasma using commercial enzyme-linked immunosorbent assays (ELISA) according to their respective package inserts. Human Quantikine Immunoassays were used for sCD14, sCD163, CRP, MCP-1, adiponectin (R&D systems, Minneapolis, Minnesota, USA), Hycult Biotech (Uden, The Netherlands) for I-FABP and MyBioSource (San Diego, CA, USA) for 1,3-β-D-Glucan. Plasma samples were diluted 250-fold for sCD14, 25-fold for sCD163, 3-fold for I-FABP and MCP-1, and 150-fold for CRP and adiponectin. For 1,3-β-D-Glucan, plasma was used neat. Assays were read at 450nm (Epoch, BioTek, Winooski, Vermont, USA) and analyzed using the Gen5 software, version 2.07. Wavelength correction was applied by resting the optical density at 570nm. Samples were run in duplicate for all assays.
2.7 Quantification of cytokines by multiplex-bead assay
Plasma samples were thawed, mixed and clarified by centrifugation (14,000 x g for 10 min at 4°C). The Human Ultrasensitive Cytokine Magnetic 10-Plex panel (ThermoFisher Scientific) was used to quantify granulocyte-macrophage colony-stimulating factor (GM-CSF), interferon [IFN]-γ, interleukin [IL]-1β, IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, and tumor necrosis factor [TNF]-α) in plasma. This 10-Plex cytokine assay was run according to the manufacturer´s instructions. All samples were run in duplicate and read with the Luminex™ 200™. Cytokine concentrations were determined from each standard curve using a 5-parameter logistic weighted algorithm (Xponent ™ v3.1 software, Luminex™, Austin, Texas, USA). Cytokines with at least 75% of values within the inferior and superior limit of detection set by each respective standard curve were used for subsequent analyses. These included IL-2, IL-4, IL-5, IL-6, IFN-γ, and TNF-α.
2.8 16S rRNA gene sequencing
DNA was extracted using the PowerSoil DNA Extraction Kit (MO BIO Laboratories, Carlsbad, CA, USA) and libraries were prepared as described in the MiSeq 16S rRNA gene Amplicon Sequencing protocol. Primers for the V3–V4 region primers (41) were: forward (341F): 5′-CCTACGGGNGGCWGCAG-3′; and reverse (805R): 5′- GACTACHVGGGTATCTAATCC-3′. Primers were purchased with overhang adapter sequences, the forward overhang sequence was: 5′- TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-followed by the 341F primer; and the reverse overhang sequence was 5′-GTCTCGTGGGCTCGGAGATGTGTATAAG AGACAG-followed by the 805R primer (Invitrogen, Waltham, MA, USA). Briefly, each DNA sample was amplified in triplicate by PCR in a final volume of 25μL. The PCR mix included the V3–V4 overhang adapter sequences (1μM), 5ng/μL of DNA, and a Platinum Taq DNA High Fidelity Polymerase (Invitrogen). PCR cycling was limited to 25 cycles. PCR amplification was verified on a 2% agarose gel (Sigma-Aldrich, MO, USA). Triplicate PCR reactions were pooled per sample and purified using AgenCourt AMPure XP beads (Beckman Dickson, Atlanta, GA, USA) as instructed by the manufacturer. Dual indices were attached by PCR to the purified PCR amplicon using the Nextera XT Index Kit (Illumina, San Diego, CA, USA). The index PCR clean-up step was repeated to guarantee complete purification of the indexed libraries using AgenCourt AMPure XP beads. Library quantification was performed using the Qubit dsDNA HS Buffer and Standards kit (Invitrogen) as instructed by the manufacturer. Next, each library was normalized to 4nM by dilution and pooled together. Libraries size and molarity were checked on a 2100 Bioanalyzer instrument (Agilent Technologies, Santa Clara, CA, USA) and sequenced on the Illumina MiSeq™ platform using a final library molarity of 14pM to obtain paired-end sequences (2×300 cycles). The internal PhiX control was used at 25% (also at 14pM). Negative controls were included in parallel to account for contamination from reagents during DNA extraction and library preparation. Control samples did not amplify by PCR nor did they yield a quantifiable library, verified on the 2100 Bioanalyzer. Library preparation, and sequencing was performed at the Centro de Investigación en Enfermedades Infecciosas, INER. Raw 16S sequences were deposited at the National Center for Biotechnology Information-Sequence Read Archive (NCBI-SRA) under project PRJNA868762 (https://www.ncbi.nlm.nih.gov/sra/PRJNA868762).
2.9 16S rRNA gene sequence analysis
Raw demultiplex FastQ files (read 1 and read 2) were processed using QIIME2 [Quantitative Insights into Microbial Ecology, version 2022.2, (42)]. First, sequences were imported into QIIME2, primers were removed using qiime cutadapt trim-paired, and checked for sequence quality and length. Next, DADA2 (43) was used to denoise, remove chimera and PhiX sequences, merge paired-end reads and construct an amplicon variant sequences (ASVs) table. Negative controls yielded less than 1000 raw sequences (read 1 and read2) and no paired-sequences were obtained after DADA2. One participant (CLWH) was excluded due to low number of sequences (<1000). After excluding this sample, a total of 4344355 paired-sequences were retained, with a median of 143,403.5 (min, 65,696-max, 238,800) per sample. Taxonomy was assigned using the SILVA database (44), release 138, at 99% identity, pre-trained for the V3-V4 16S rRNA region. A phylogenetic tree was constructed using qiime phylogeny align-to-tree-mafft-fasttree. Qiime2 artifacts were imported into R for further manipulation, and graph visualization (R version 4.3.0). Rarefaction was performed at a sampling depth of 59,126 sequences per sample. Alpha-diversity was estimated using Shannon index, observed ASVs (richness), and Faith phylogenetic distance, while beta-diversity was estimated using Bray-Curtis dissimilarity index (vegan v.2.6-4, R). We also explored other beta-diversity metrics (unweighted and weighted Unifrac) and found they compared well with Bray-Curtis, with no additional information gained from these two metrics (data not shown). Beta-diversity was visualized using principal coordinate analysis (PCoA). Permutational Multivariate Analysis of Variance (PERMANOVA, adonis2, vegan, with 10,000 permutations) was used to test differences between groups. Multi-factor models were also run with multiple explanatory variables (sex, HIV status and age) using the adonis2() function to test for the interaction between these 3 variables. Homogeneity of group dispersion was verified with the betadisper function (vegan, R). Linear discriminant analysis (LDA) effect size [LEfSe, (45)] was performed using default parameters and an LDA score equal to or above 4. To test for discriminant taxa by HIV status (and sex assigned at birth: female vs male, hereby referred to as sex only), LEfSe was performed on the unrarefied ASV table, filtered to remove ASV found in fewer than 10% of samples. This analysis was run on the Galaxy web application (https://huttenhower.sph.harvard.edu/galaxy/).
2.10 Statistical analysis
Data were expressed as number (percentage), median (interquartile range), mean (plus or minus standard deviation) or number (percentage) as appropriate. Groups were compared using the Wilcoxon Rank Sum Test for continuous variables and Fisher´s Exact test or Chi-square test (as required) for categorical values (GraphPad Prism v8, San Diego, California, United States). Correlation coefficients and p values were computed with the nonparametric Spearman test using “cor.test” function in R, p values were adjusted using the false discovery rate (FDR). Regression analyses were used to assess relationships between HIV status and numeric outcome variables using the “lm” function in R, and adjusted for sex, age, BMI/A (model 1) or sex, age, BMI/A and antibiotic use within the last 3 months (model 2). Given the impact of sex on the gut microbiota structure (as assessed by PERMANOVA), we also fitted regression analyses between numeric outcome variables of interest and sex using the “lm” function, adjusting for HIV status, age, BMI/A (model 1) and HIV status, age, BMI/A and antibiotic exposure within the last 3 months (model 2). Data were visualized using “ggplot2”, ggcorrplot2” and “ampvis2”.
3 Results
3.1 Cohort characteristics
A total of 31 children were enrolled; 15 were CLWH and 16 were HUU. Children were aged between 3 and 15 years old, with a mean of 9.53 years. Table 1 summarizes the cohort characteristics. This cohort was well-balanced in terms of number of participants per group, sex, age, BMI/A, and immune status (as assessed by CD4+ T cell count). All CLWH were on antiretroviral therapy (ART) for a median of 8.65 years (interquartile range 7.57-10.75) and were virally suppressed (plasma viral load <200 copies/mL). Ten CLWH were on protease inhibitors (PI, 66.7%), 4 on integrase strand transfer inhibitors (InSTI, 26.7%) and 1 on non-nucleoside reverse transcriptase inhibitors (NNRTI, 6.7%). No CLWH had CD4+ T cell counts below 200 cells/µL. Supplementary Table 2 summarizes the HIV clinical parameters of CLWH at HIV diagnosis and at study enrollment.
Table 1 Cohort characteristics.
3.2 Dietary patterns and dietary quality
Overall, children in this study (CLWH and HUU alike) had a Western-type diet, characterized by excessive intake of refined carbohydrates, animal protein, sugar, saturated fats and inadequate fiber consumption (Supplementary Tables 3, 4). Ten children had overweight and 2 had obesity (according to their BMI/A). The overall median HEI score was 67.06 (interquartile range 58.76-74.66). Two children (1 CLWH and 1 HUU) had a healthy diet. All other children had poor diets or needed changes (Table 1). After adjusting for multiple comparisons, no differences in dietary pattern were seen between CLWH and HUU (Supplementary Tables 3, 4).
3.3 Immune senescence remains elevated in CLWH despite antiretroviral therapy
Next, we assessed markers linked to HIV disease progression and pathogenesis, and overall health. We quantified soluble and cellular markers in plasma and peripheral blood mononuclear cells respectively, to determine levels of microbial translocation (bacterial and fungal), gut permeability, immune activation (plus senescence, exhaustion and cell cycling) and inflammation. We also included markers of adipose tissue inflammation (sCD163, adiponectin and MCP-1). The median and interquartile range of each marker is shown in Table 2. When comparing CLWH and HUU, no differences were found, except for the frequency of exhausted CD4+ T cells (PD-1+), which was elevated in CLWH compared with HUU (p=0.0092). The association between HIV status and the frequency of exhausted CD4+ T cells (PD-1+) was further evaluated using regression analysis and remained significant after adjusting for sex, age, BMI/A, and recent antibiotic use (p=0.0024, Supplementary Table 5).
Table 2 Soluble and cellular markers linked to HIV disease progression and pathogenesis.
3.4 No discernable impact of HIV on the gut microbiota
Next, we profiled the gut microbiota using fecal samples and 16S rRNA sequencing. Alpha diversity was higher in CLWH for all three indices (richness, shannon and phylogenetic distance), although no differences were found when comparing with HUU (Supplementary Figure 2A). Regression analyses confirmed no clear association between HIV and alpha diversity (Supplementary Table 6). Next, bacterial community structure was visualized using principal coordinate analysis (PCoA, Bray-Curtis dissimilarity index). We found no clear evidence of clustering by HIV status (CLWH compared with HUU, R-squared= 0.034, PERMANOVA p=0.46, Supplementary Figure 2B). The Bray-Curtis dissimilarity index was similar between CLWH and HUU, showing a high degree of interindividual variation within groups (Supplementary Figure 2C). Next, we explored the taxonomic composition at phylum level (Supplementary Figure 3). Firmicutes and Bacteroidota were predominant with an overall mean relative abundance of 50.67% and 44.26%, respectively. These two phyla accounted for over 94% of all phyla present. At genus level, Bacteroides was predominant (29.69%), followed by Faecalibacterium (10.24%), Prevotella (6.23%), Blautia (4.18%), Agathobacter (3.27%), Parabacteroides (2.60%), Alistipes (2.39%), UCG-002 (Oscillospiraceae, 1.94%), Ruminococcus (1.93%), and Roseburia (1.79%). Figure 1A summarizes the overall relative abundance of the top 40 genera; these account for 89.44% of all genera present. Heatmap showing the relative abundance of the top 40 genera in each child is shown in Supplementary Figure 4. To assess if there were discriminant ASVs between CLWH and HUU (i.e., by HIV status), we performed linear discriminant analysis (LDA) effect size (LEfSe, using default parameters, the threshold on the LDA was set to 4 or above). No discriminant taxa between CLWH and HUU were found using LEfSe.
Figure 1 Sex impacts the gut microbiota of children independently of HIV status. (A) Heatmap showing the mean overall relative abundance of the top 40 genera stratified by HIV status and faceted by sex. (B) PCoA showing the impact of sex on the gut microbiota structure based on the Bray-Curtis dissimilarity index (PERMANOVA p=0.024). Sex explains 5.2% of the variation in bacterial community structure. Points are colored by sex (females in red, males in blue), ellipses represent 95% confidence interval for each group. (C) Alpha diversity metric comparisons between female and male children using Wilcoxon rank sum test, showing the male children have increased richness compared to female children. (D) Linear discriminant analysis (LDA) effect size (LEfSe), showing the two taxa that discriminate female and male children: Firmicutes for male children, Bacteroides for female children. Graphs were created using ampvis2 version 2.8, and ggplot2 version 3.4.2. %, percentage; CLWH, children living with HIV; F, female; HUU, HIV-unexposed and uninfected children; M, male; PCoA, principal coordinate analysis; Faith's phylogenetic diversity (PD) PERMANOVA, Permutational Multivariate Analysis of Variance.
3.5 Sex impacts the gut microbiota of children
Next, we explored the influence of other variables known to impact the gut microbiota, including sex, age, BMI/A, previous antibiotic use, ART, place of residence, immune status, maternal factors, siblings, and stool consistency. Information regarding previous antibiotic usage and stool consistency is summarized in Supplementary Tables 6, 7, respectively. To assess the impact of these variables on bacterial community membership, we used PERMANOVA (10,000 permutations). Results on the impact of each variable in single-variable models are shown in Supplementary Table 8. Sex was the only variable that significantly impacted the gut microbiota (R-squared= 0.052, PERMANOVA p=0.024, Figure 1B), though sex accounted for a relatively small amount of the variation (5.2%) and the PERMANOVA p was not adjusted for multiple comparisons. While HIV status and age did not significantly impact the gut microbiota in single-variable models, we also explored the interaction of sex, HIV status and age using adonis2 function. This multi-factor analysis yielded similar results to the single-variable model (R2 = 0.048, PERMANOVA p=0.044), suggesting that sex impacts the gut microbiota independently of age and HIV status. We also found that male children had higher alpha diversity (richness only) compared with female children (Figure 1C). Two ASV were found to discriminate female from male children using LEfSe (LDA>4): Firmicutes for males, and Bacteroides for females (Figure 1D). To confirm the discriminatory value of these two ASVs and the impact of sex on microbial richness we used regression analyses. After adjusting for HIV status, age, BMI/A and recent antibiotic use, associations between sex and microbial richness, Firmicutes and Bacteroides were lost (Supplementary Table 9).
3.6 Proteobacteria are associated with CD4+ T cell activation and exhaustion
Finally, we related the gut microbiota to markers of gut health, HIV disease progression and pathogenesis, dietary patterns and host metadata. Relationships were explored by nonparametric correlation to provide an overview of microbiota-diet-host associations. We limited the number of correlations and included the 4 predominant phyla (Firmicutes, Bacteroidota, Proteobacteria, and Actinobacteria), the top 10 genera, alpha diversity, most markers linked to HIV disease progression and pathogenesis, diet-related markers (BMI/A, HEI score, daily caloric intake, 3 macronutrients, sugar) and host metadata (age, immune status, antibiotic usage). Correlations were performed using nonparametric Spearman, and p values were corrected [false discovery rate (FDR)] as shown in Figure 2. We found strong correlations (positive and negative) among Firmicutes, Bacteroidota, and their respective genera, and between Firmicutes, Bacteroidota (and their respective genera) and alpha diversity. In particular, the genus UCG-002 was positively correlated with all 3 indices of alpha diversity (rho>0.8, FDR-adjusted p<0.001 for all 3 correlations). Bacteroides was inversely correlated with Prevotella (rho=-0.822, FDR-adjusted p<0.001) and Roseburia (rho=-0.63, FDR-adjusted p=0.0084) and Shannon (rho=-0.587, FDR-adjusted p=0.025). Both the frequency of activated CD4+ T cells (as measured by CD38+ HLADR+) and exhausted CD4+ T cells (as measured by PD-1+) were positively correlated with Proteobacteria (rho=0.568. FDR-adjusted p= 0.029, and rho=0.62, FDR-adjusted p=0.0126, respectively).
Figure 2 Heatmap showing relationships between the gut microbiota, markers of HIV disease progression/pathogenesis, diet and host metrics. Correlations were performed using nonparametric Spearman, and p values were adjusted for multiple correlations using the false discovery rate (FDR). We included the 4 predominant phyla (Firmicutes, Bacteroidota, Proteobacteria, and Actinobacteria), the top 10 genera, alpha diversity metrics, markers linked to HIV disease progression and pathogenesis, diet and host metadata. Correlations with an unadjusted p>0.05 were removed for ease of visualization (empty cells). Significant FDR-adjusted p values are marked with a black box, and the level of significance (* FDR-adjusted p<0.05, ** FDR-adjusted p<0.01, *** FDR-adjusted p<0.001) is provided. Positive correlations are shown in blue, and negative correlations are shown in red. The heatmap was created using ggcorrplot2 version 0.1.2 in R. %, percentage; Abx, antibiotics; β, beta; BMI/A, body mass index for age; CD, cluster of differentiation; CRP, C-reactive protein; FDR, false discovery rate; I-FABP, intestinal fatty acid binding protein; HEI, healthy eating index; MCP-1, monocyte chemoattractant protein-1; PD, phylogenetic distance; PD-1, programmed cell death-1; s, soluble.
4 Discussion
The work presented here characterized the gut microbiota in CLWH, and assessed its relation to markers of HIV disease progression and pathogenesis (and others), dietary patterns and host metadata, comparing it with HUU. We found no clear evidence that HIV impacted the gut microbiota composition, structure or diversity, nor did we find discriminant taxa between CLWH and HUU. Our results are both in agreement and in contrast with previously published studies in CLWH. Kaur and colleagues (22) showed that alpha diversity was decreased in CLWH compared to uninfected children, in both CLWH naïve to ART and CLWH on ART. CLWH naïve to ART had the least alpha diversity. They also reported that CLWH had a distinct pattern of gut microbiota compared to uninfected children, characterized by increased Bacteroidetes, and decreased Firmicutes, and Actinobacteria at phylum level, and increased Prevotella, and decreased Bacteroides, and Ruminococcus at genus level. Abange and colleagues (24) also reported that CLWH had lower alpha diversity and altered beta diversity. They found that CLWH had lower relative fractions of butyrate producers and obligate anaerobes, that Prevotella was decreased, together with gammaproteobacteria (Haemophilus influenzae) and Enterobacteriaceae, while Akkermansia and Faecalibacterium were increased. In contrast, Flygel and colleagues (25) suggested that prolonged ART (> 10 years) might restore microbial richness in CLWH. Flygel and colleagues also reported increased beta-diversity and a rectal microbiota enriched in 3 genera from the phylum Actinobacteria (Corynebacterium, Lawsonella, Collinsella) and 4 genera from the phylum Firmicutes (Finegoldia, Anaerococcus, Erysipelotrichaceae, and Lachnoclostridium). Sessa and colleagues (46) reported higher alpha diversity in CLWH, and two distinct clusters characterized by Akkermansia muciniphila in cluster 1 and Bacteroides, Faecalibacterium (and others) in cluster 2. Goosen et al. (26), reported that alpha diversity did not differ between CLWH and uninfected controls, while beta diversity was significantly higher in CLWH. Finally, Tincati and colleagues (27) found similar bacterial microbiota composition and alpha diversity in fecal samples. They also interrogated the plasma microbiome and found no difference in alpha and beta diversity. However, they did find a predominance of Gammaproteobacteria and Pseudomonas spp. in HIV-exposed infected compared to HIV-exposed uninfected children.
Discrepancies between studies are expected. These could be due to differences in the number of participants included (small versus large cohorts), age of children, inclusion of CLWH naïve to ART, or only CLWH on ART, specific ART regimens and length of ART duration, whether CLWH on ART were virally suppressed or viremic, use of antibiotics, and the choice of HIV-uninfected controls [whether HIV-unexposed and uninfected or HIV-exposed and uninfected (27)]. Furthermore, other variables including geographic location and even ethnic groups within the same country (24) are also likely to impact results, precluding their generalization. Kaur and colleagues included both CLHW naïve to ART and on suppressive ART and compared them to age- and sex-matched uninfected children. In this cohort, children were excluded if antibiotics were used one month prior to sample collection. Abange and colleagues included CLWH on ART, predominantly NNRTI-based regimen, for a median of 5 years and non-exposed HIV-negative children. Some CLWH were viremic (n=13) and most had CD4>500 cells/mm3. Flygel and colleagues investigated the gut microbiota in CLWH on long-term ART (most CLWH had chronic lung disease) and HIV-uninfected children, and used rectal swabs instead of fecal samples. Sessa and colleagues included aviremic CLWH for a median of 9 years (46). Goosen and colleagues (26) included aviremic CLWH aged 8 to 13 years, and excluded children with recent antibiotic use, and those with BMI/A <-3 or >2. Finally, Tincati and colleagues (27) included HIV-exposed infected children all on ART but not all aviremic, and HIV-exposed uninfected children (born from HIV-infected mothers), and excluded children with previous antibiotic exposure (2 weeks), gastrointestinal disorders and co-infections. Compared to other studies, our cohort included a small number of participants, CLWH were all aviremic, on PI-based ART for the most part. Interestingly, Abange and colleagues (24) reported that CLWH on PI-based regimen had lower alpha diversity and altered composition when compared to CLWH on NNRTI, as well as their HIV non-infected counterparts. In particular, they found that Blautia, Prevotella, Oscillospira, and Faecalibacterium were more abundant in NNRTI compared to PI/r, with CLWH on PI/r had a lower abundance of microorganisms predicted to be butyrate producers and obligate anaerobes. Given the predominant use of PI-based regimen in our cohort and the small number of CLWH, subgroup analysis looking at the potential impact of ART on the gut microbiota was not possible. Interestingly, we reported differences in the gut microbiota of Mexican ALWH on PI/r-based and EFV-based regimen (36). Although alpha diversity was not significantly different between ART regimens, we did find higher microbial translocation (assessed by plasma sCD14) and enterocyte damage (assessed by plasma levels of I-FABP) in those ALWH on PI/r-based regimen. Flygel and colleagues found no association between beta diversity and type of ART, nor years on ART (25). Both studies were however in agreement with the fact that duration of ART was not associated with alpha and beta diversity, suggesting that in pediatric HIV and viremia were not predictive of microbiota diversity (12), in contrast with ALWH (47, 48).
Some studies, including ours, suggest a beneficial impact of long-term ART on the gut microbiota. Some genera associated with intestinal health were also found in CLWH [Akkermansia, Faecalibacterium, (24)]. In agreement with other studies (25–27), we found that long-term ART might restore all metrics of alpha diversity. The median time on ART was similar between these studies and ours. We found no differences in alpha and beta diversity, and a high abundance of Faecalibacterium. Akkermansia was not among the top 20 genera in our cohort, but we were able to detect this genus, the mean overall abundance was 0.2% in CLWH and 0.7% in HUU. Faecalibacterium (Firmicutes) and Akkermansia (Verrucomicrobia) are two genera thought to be beneficial, due to their protective role in metabolism, immune system and gut barrier (49). Kaur and colleagues found decreased abundance of Faecalibacterium in CLWH on ART compared to HIV-uninfected controls (22). Interestingly, Sessa and colleagues found that the gut microbial group represented predominantly by Akkermansia (A. muciniphila) was linked to higher markers of inflammation, vascular endothelial activation, microbial translocation, but not T cell activation (46). Their results suggest that overabundance of a beneficial commensal could have negative implications. Indeed, Akkermansia has been linked to deleterious effects (50).
We report that markers linked to HIV disease progression and pathogenesis were similar to those of HUU, which suggests that prolonged ART might normalize their levels. Our results are in agreement with Fitzgerald and colleagues (51) who found that immune activation decreased in pediatric HIV-infected children after ART initiation. Weinberg and colleagues also found that several markers including IFN-, IL-6, TNF-, sCD163, CRP, sCD14 and MCP-1 decreased after virologic control (52). Others found that ALWH on stable ART showed evidence of altered gut permeability and fungal translocation (53) and increased monocyte and T cell activation and higher sCD14 (54). In our study, no differences in inflammation, immune activation, bacterial and fungal translocation, and enterocyte damage were observed between CLWH and HUU. However, our results also differ from Tincati and colleagues (27). They found elevated levels of sCD14, IL-6 in HIV-exposed infected children compared with HIV-exposed uninfected children, however levels of I-FABP were similar. They raise an important point in their discussion regarding the choice of the control uninfected group used to compare with HIV-exposed infected children, as we alluded to earlier. Most studies enroll HIV-unexposed and uninfected children. As shown in ALWH, careful consideration is needed to enroll the most appropriate control group (55). Studies in perinatal HIV should ideally include both HIV-exposed uninfected and HIV-unexposed uninfected children, when possible. We did find increased frequency of exhausted T cell, a marker of HIV disease progression in CLWH. This is in agreement with other studies in CLWH (56, 57). Exhausted T cells are also a marker of disease progression in ALWH (58). In CLWH, PD-1+ T cells are heterogenous population of central and effector memory T cells that have impaired proliferative capacity and secrete proinflammatory cytokines. Interestingly, we also found that Proteobacteria, a phylum known to harbor pathobionts, was associated with the frequency of activated CD4+ T cells and exhausted CD4+ T cells, linking the gut microbiome to persistent immune exhaustion and premature aging. In our cohort, the mean overall abundance of Proteobacteria was 2.40%, when stratifying by HIV, we found an increase in CLWH (2.8%) compared with HUU (2.1%). At genus level, we found an unclassified Alphaproteobacteria, Sutterella (order of Burkholderiales) and Escherichia-Shigella (family Enterobacteriaceae), known to translocate and elicit pro-inflammatory responses (59, 60). Interestingly, Tincati and colleagues (27) found that the plasma microbiome was enriched in Gammaproteobacteria and Pseudomonas spp. in HIV-exposed infected children and Alphaproteobacteria in HIV-exposed uninfected children. They also reported a weak trend towards an inverse correlation between plasma Pseudomonas spp. and CD4+ T cell counts. Collectively, our results suggest that although some markers and gut microbiota metrics might normalize or differences may not be apparent or quantifiable (in particular if studies are underpowered like ours), CLWH are at increased risk of premature aging. Also, the CD4/CD8 ratio did not normalize, which is concordant with the fact that most CLWH in our study started ART later. Our findings could also be attributable, at least in part to the inclusion of CLWH and HUU who had recently taken antibiotics. Interestingly, studies have shown that co-trimoxazole increased alpha diversity (23), decreased inflammation (61) and suppressed the gut-resident viridans group streptococcal species (VGS) associated with intestinal inflammation (62). At the same time, antibiotics are a confounder in gut microbiota studies (63) and have long-term adverse effects of the gut microbiota (32). Most studies report that antibiotic use causes a reduction in microbial diversity and richness, including a decrease in beneficial bacterial taxa, resulting in metabolic alterations, increased susceptibility to colonization by pathogens and an increase in bacterial resistance to antibiotics. Children in this study (CLWH and HUU alike) reported a high number of antibiotics courses, even considering the recall bias that parents might have had at the time of the interview. Antibiotics taken by children in this study were primarily β-lactams, broad-spectrum antibiotics which interrupt bacterial cell-wall formation (64). When assessing the impact of antibiotics on the gut microbiota, we found none. Anecdotally, stool consistency was not indicative of diarrhea associated with antibiotics (type 6-7), most children had normal stool consistency (type 3-4). These results were unexpected as the cumulative effects on the microbiome as a result of repeated administration of multiple courses of antibiotics have been demonstrated29,30. Consequences have also been reported with children facing the risk of experiencing temporary or chronic alterations in their microbiome, which could predispose them to develop chronic disorders 27,28. Our findings could be attributed, at least in part, to the fact that amoxicillin-clavulanate was the most frequently prescribed antibiotic, and it has been reported to cause no or only mild to moderate changes in the composition of the intestinal microbiome (65, 66).
Diet is one of the key modulators of the gut microbiota (63). In this study, children had a diet that is representative of the Mexican urban population (67). Both diet and its effect on the gut microbiota have clearly been linked to obesity, metabolic syndrome, and type 2 diabetes (68, 69). We found no associations between dietary intake amounts or diet quality (HEI score) and gut microbiota metrics in the context of HIV infection. We also found that markers of adipose tissue inflammation were comparable between CLWH and HUU. Generoso and colleagues also found that sCD163 normalized in Kenyan CLWH on ART (70). According to the Encuestas Nacionales de Salud y Nutrición (Ensanut, National Health and Nutrition Surveys), Instituto Nacional de Salud Pública (INSP, National Institute of Health) in 2018-2019, the prevalence of overweight and obesity in children aged 5 to 11 years was 18.1 and 17.5% (respectively) and in adolescents aged 12 to 19 years was 23.8 and 14.6%, respectively (71). In our study, no CLWH had obesity, compared to 2 HUU who had obesity (no significant differences were observed between groups). Childhood obesity is a growing health problem worldwide including Mexico (72). Childhood obesity in the context HIV infection is further alarming given the increased prevalence of cardiometabolic diseases and other NCD. Improving the health of CLWH (and HUU) is very important. One possible avenue is dietary interventions, which can be very successful short-term, but are not largely sustainable over time (73, 74). Food and nutrition literacy should only be explored, as empowering children with knowledge and skills related to nutrition and dietary habits, are likely to promote healthier eating behaviors from an early age (75). Interventions targeting the gut microbiota using different approaches have had relatively little success so far (13, 76). A recent study reported an attenuation of bacterial dysbiosis in vertically HIV-infected children (n=24), randomized to receive a daily nutritional supplementation or placebo for 4 weeks (77). Authors postulated that the gut microbiota is less resilient in children and better amiable to change and therefore would respond better to interventions. Interventions to promote a healthier gut microbiome, together with healthier lifestyle to increase life expectancy, quality in CLWH is necessary as a pediatric HIV cure is still lacking (78). CLWH will face many challenges as they transition to adulthood, and efforts should be made to minimize NCD and other comorbidities, which will negatively impact their quality of life. Interestingly, no correlations were found between bacterial genera, alpha diversity and dietary variables. We did find that alpha diversity metrics were positively correlated with UCG-002, a genus of the Oscillospiraceae family. UCG-002 has been negatively associated with obesity and insulin resistance, suggesting that UCG-002 may be involved in mechanisms that improve insulin sensitivity (79, 80). Another study reported an enrichment of UCG-002 in pubertal participants, as well as an increased abundance of Bacteroides and Prevotella in individuals at different stages of puberty (81). The authors suggest that these bacterial genera may represent a more mature stage of the gut microbiota. It is therefore plausible, though highly speculative, that the positive association between UCG-002 and alpha diversity may indicate, at least in part, the start of the transitional change from childhood to adulthood, impacting the gut microbiota. The mean of our cohort was 9.53 years.
We found that sex impacted the gut microbiota with female children being Bacteroides predominant and having lesser microbial richness compared to their male counterparts, although after adjusting for HIV status, age, BMI/A and recent antibiotic use, these associations were lost. Recent studies have shown that sex (and sex hormones) contribute to differences in the gut microbiome (82), as early as birth and throughout life. Interestingly, differences in inflammation profiles by sex have also been reported (52, 83). Because our primary objective was to understand the impact of HIV infection, associations between markers of HIV disease progression and sex or age were not explored. Differential effects of age have also been reported (84). We acknowledge that this is a limitation. In HIV infection, Bacteroides is a common enterotype, together with Prevotella, this latter is linked to sexual preference in ALWH (55). Bacteroides abundance has also been linked to consumption of protein- and animal-fat rich diets (85) and an increased Prevotella/Bacteroides has been linked to obesity (86). Interestingly, a study in Mexican children also reported the Bacteroides was highly abundant in children with type 1 diabetes (87). Additionally, we found a negative correlation between Prevotella and Bacteroides, which has been previously reported (86, 88). Our understanding of which variables are responsible for these taxonomic shifts remains incomplete. Indeed, Sessa and colleagues (46) failed to identify clinical variables that would explain the presence of two bacterial clusters in CLWH (A. muciniphila, less diverse, and Bacteroides and other genera, more diverse). Factors that impact the gut microbiota are likely to have synergistic and cumulative effects, making teasing out the impact of each factor difficult. Disentangling these effects are very important in the immediate future to improve our understanding of HIV-associated microbial dysbiosis.
This study has several limitations. First, our cohort included a small number of participants, and is underpowered to detect differences. Also, our cohort was not matched on age. The inability to balance groups on confounders might have introduced bias in our analyses. Second, CWH and HUU were recently exposed to antibiotics, which might have masked differences in bacterial community structure. Also, although our results suggest microbial signatures linked to CD4+ T cell activation and exhaustion (Proteobacteria) and sex (Bacteroides), these should be interpreted carefully given the high degree of interindividual variation in gut microbiota, which are not accurately represented when visualizing data (mean relative abundance) using heatmaps or barplots. Also, we were limited to reporting taxonomic classification to genus level, a limitation inherent to the16S targeted sequencing. Subgroup analyses to explore the impact of different variables (ART, diet) were not possible due to low number of participants in each subgroup. Finally, our study is cross-sectional and does not allow to establish causality. At the same time, our study also has strengths, as we included a comprehensive use of biologics: dietary patterns and quality; several markers of inflammation, immune activation (soluble and cellular), microbial and fungal translocation, enterocyte integrity, and adipose tissue inflammation, host metadata and explored relationships with the gut microbiota.
5 Conclusion
Given the prolonged exposure to HIV and ART, and prophylactic exposure to antimicrobials, the gut microbiota of CLWH may be less resilient, increasing tolerance to pathobionts. This, in turn, may lead to increased prevalence of cardiometabolic alterations, and NCD, or their premature symptomatic presentation compared to the general population. Here we report that the gut microbiota of CLWH appears similar to that of HUU, and that most markers of HIV disease progression normalize on ART. Our findings suggest that long-term ART administration (8 years) in CLWH, might have a beneficial impact on the gut microbial ecology and immune activation/inflammation. Sex is also likely to impact the gut microbiota. CLWH do exhibit signs of premature aging (frequency of exhausted CD4+ T cells), which puts them at risk for increased comorbidities in adulthood. Interventions to reverse HIV-associated premature aging, in synergy with interventions targeting overweight and obesity are likely to positively impact the quality of life and longevity of CLWH.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: PRJNA868762 (SRA).
Ethics statement
The studies involving humans were approved by Ethics and Scientific committees of INER and IMSS. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from the recruitment of participants to a research protocol. Written informed consent and age-appropriate assent were obtained from all participants and/or parents before study enrollment.
Author contributions
VB-M, SP-C, CR-P: Designed the study. CR-P, NAMJ, HP-L, JV-R: Recruited participants. CR-P: performed interviews and questionnaires. CR-P, IO-P: Performed dietary assessments. CA-A: Performed laboratory determination of lymphocyte populations. MC-T, NM, SP-C: Performed laboratory experiments. SP-C, CR-P: Analyzed the data. VB-M, GR-T, SA-R: Provided critical intellectual contributions. GRT, VB-M, SP-C: Funding acquisition. CR-P, SP-C: Wrote the manuscript with contributions from all co-authors. All authors approved the final version of this manuscript.
Funding
This work was supported by grants from Apoyo financiero para el desarrollo de protocolos de investigación y desarrollo tecnológico sobre temas específicos de los temas prioritarios de salud del Instituto Mexicano del Seguro Social (Fondo de Investigación en Salud FIS/IMSS/PROT/PRIO/19/097 R-2018-785-130), Programa Presupuestario F003 “Programas Nacionales Estratégicos de Ciencia,. Tecnología y Vinculación con los Sectores Social, Público y Privado” Dirección de Ciencia de Frontera, Consejo Nacional de Humanidades Ciencias y Tecnologías (CONAHCYT, Convocatoria Ciencia de Frontera 2019 formely FORDECYT-PRONACES/140641/2020) and the Programa Presupuestal P016, Anexo 13 del Decreto del Presupuesto de Egresos de la Federación.
Acknowledgments
We thank all children, parents and guardians who decided to participate in this study. Authors would also like to thank Perla M. Del Rio Estrada, PhD (CIENI, INER), Dafné Diaz-Rivera (Biol., CIENI, INER), for helpful discussion on flow cytometry data analysis, and Olivia Briceño, PhD (CIENI, INER), for technical assistance with the flow cytometer. Cecilia Rosel-Pech is a doctoral student from the Posgrado en Ciencias Biológicas, Universidad Nacional Autónoma de México. CRP thanks the Consejo Nacional de Humanidades Ciencias y Tecnologías (CONAHCYT) for her doctoral scholarship (CVU 545509). The authors would also like to thank Dr. María Teresa Ramos Cervantes for her support in the logistics for clinical data collection, and to laboratory technician Nayeli Cabrera Velázquez for her contribution in blood sample collection, both affiliated with the National Medical Center “La Raza,” IMSS. This paper is part of the requirements for obtaining a Doctoral Degree at the Posgrado en Ciencias Biológicas, UNAM.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1244473/full#supplementary-material
References
1. HIV. Estimated number of people (all ages) living with HIV (2023). Available at: https://www.who.int/data/gho/data/indicators/indicator-details/GHO/estimated-number-of-people–living-with-hiv.
2. Centro Nacional para la Prevención y Control del VIH y el sida | Gobierno | gob.mx (2023). Available at: https://www.gob.mx/censida.
3. Marcus JL, Leyden WA, Alexeeff SE, Anderson AN, Hechter RC, Hu H, et al. Comparison of overall and comorbidity-free life expectancy between insured adults with and without HIV infection, 2000-2016. JAMA Netw Open 15 junio (2020) 3(6):e207954. doi: 10.1001/jamanetworkopen.2020.7954
4. Peters BA, Burk RD, Kaplan RC, Qi Q. The gut microbiome, microbial metabolites, and cardiovascular disease in people living with HIV. Curr HIV/AIDS Rep abril (2023) 20(2):86–99. doi: 10.1007/s11904-023-00648-y
5. MacCann R, Landay AL, Mallon PWG. HIV and comorbidities - the importance of gut inflammation and the kynurenine pathway. Curr Opin HIV AIDS (2023) 18(2). doi: 10.1097/COH.0000000000000782
6. Brenchley JM, Price DA, Schacker TW, Asher TE, Silvestri G, Rao S, et al. Microbial translocation is a cause of systemic immune activation in chronic HIV infection. Nat Med diciembre (2006) 12(12):1365–71. doi: 10.1038/nm1511
7. Steele AK, Lee EJ, Vestal B, Hecht D, Dong Z, Rapaport E, et al. Contribution of intestinal barrier damage, microbial translocation and HIV-1 infection status to an inflammaging signature. PloS One (2014) 9(5):e97171. doi: 10.1371/journal.pone.0097171
8. 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(193):193ra91. doi: 10.1126/scitranslmed.3006438
9. 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(3):311–22. doi: 10.1016/j.chom.2016.02.011
10. Gelpi M, Vestad B, Hansen SH, Holm K, Drivsholm N, Goetz A, et al. Impact of human immunodeficiency virus-related gut microbiota alterations on metabolic comorbid conditions. Clin Infect Dis Off Publ Infect Dis Soc Am (2020) 71(8):e359–67. doi: 10.1093/cid/ciz1235
11. Zevin AS, McKinnon L, Burgener A, Klatt NR. Microbial translocation and microbiome dysbiosis in HIV-associated immune activation. Curr Opin HIV AIDS (2016) 11(2):182–90. doi: 10.1097/COH.0000000000000234
12. Tuddenham S, Koay WL, Sears CHIV. Sexual orientation, and gut microbiome interactions. Dig Dis Sci (2020) 65(3):800–17. doi: 10.1007/s10620-020-06110-y
13. Rosel-Pech C, Chávez-Torres M, Bekker-Méndez VC, Pinto-Cardoso S. Therapeutic avenues for restoring the gut microbiome in HIV infection. Curr Opin Pharmacol octubre (2020) 54:188–201. doi: 10.1016/j.coph.2020.09.010
14. Vemuri R, Shankar EM, Chieppa M, Eri R, Kavanagh K. Beyond just bacteria: functional biomes in the gut ecosystem including virome, mycobiome, archaeome and helminths. Microorganisms (2020) 8(4):483. doi: 10.3390/microorganisms8040483
15. Hager CL, Ghannoum MA. The mycobiome in HIV. Curr Opin HIV AIDS (2018) 13(1):69–72. doi: 10.1097/COH.0000000000000432
16. Li S, Yang X, Moog C, Wu H, Su B, Zhang T. Neglected mycobiome in HIV infection: Alterations, common fungal diseases and antifungal immunity. Front Immunol (2022) 13:1015775. doi: 10.3389/fimmu.2022.1015775
17. Stern J, Miller G, Li X, Saxena D. Virome and bacteriome: two sides of the same coin. Curr Opin Virol (2019) 37:37–43. doi: 10.1016/j.coviro.2019.05.007
18. 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(4):562–70. doi: 10.4161/gmic.32132
19. Pinto-Cardoso S, Klatt NR, Reyes-Terán G. Impact of antiretroviral drugs on the microbiome: unknown answers to important questions. Curr Opin HIV AIDS (2018) 13(1):53–60. doi: 10.1097/COH.0000000000000428
20. Villanueva-Millán MJ, Pérez-Matute P, Recio-Fernández E, Lezana Rosales JM, Oteo JA. Differential effects of antiretrovirals on microbial translocation and gut microbiota composition of HIV-infected patients. J Int AIDS Soc (2017) 20(1):21526. doi: 10.7448/IAS.20.1.21526
21. Ishizaka A, Koga M, Mizutani T, Parbie PK, Prawisuda D, Yusa N, et al. Unique gut microbiome in HIV patients on antiretroviral therapy (ART) suggests association with chronic inflammation. Microbiol Spectr (2021) 9(1):e0070821. doi: 10.1128/Spectrum.00708-21
22. Kaur US, Shet A, Rajnala N, Gopalan BP, Moar P, D H, et al. High Abundance of genus Prevotella in the gut of perinatally HIV-infected children is associated with IP-10 levels despite therapy. Sci Rep (2018) 8(1):17679. doi: 10.1038/s41598-018-35877-4
23. D’Souza AW, Moodley-Govender E, Berla B, Kelkar T, Wang B, Sun X, et al. Cotrimoxazole prophylaxis increases resistance gene prevalence and α-diversity but decreases β-diversity in the gut microbiome of human immunodeficiency virus-exposed, uninfected infants. Clin Infect Dis Off Publ Infect Dis Soc Am (2020) 71(11):2858–68. doi: 10.1093/cid/ciz1186
24. Abange WB, Martin C, Nanfack AJ, Yatchou LG, Nusbacher N, Nguedia CA, et al. Alteration of the gut fecal microbiome in children living with HIV on antiretroviral therapy in Yaounde, Cameroon. Sci Rep (2021) 11(1):7666. doi: 10.1038/s41598-021-87368-8
25. Flygel TT, Sovershaeva E, Claassen-Weitz S, Hjerde E, Mwaikono KS, Odland JØ, et al. Composition of gut microbiota of children and adolescents with perinatal human immunodeficiency virus infection taking antiretroviral therapy in Zimbabwe. J Infect Dis (2020) 221(3):483–92. doi: 10.1093/infdis/jiz473
26. Goosen C, Proost S, Baumgartner J, Mallick K, Tito RY, Barnabas SL, et al. Associations of HIV and iron status with gut microbiota composition, gut inflammation and gut integrity in South African school-age children: a two-way factorial case-control study. J Hum Nutr Diet Off J Br Diet Assoc (2023) 36(3):819–32. doi: 10.1111/jhn.13171
27. Tincati C, Ficara M, Ferrari F, Augello M, Dotta L, Tagliabue C, et al. Gut-dependent inflammation and alterations of the intestinal microbiota in individuals with perinatal HIV exposure and different HIV serostatus. AIDS Lond Engl 15 noviembre (2022) 36(14):1917–25. doi: 10.1097/QAD.0000000000003324
28. World Health Organization. Guidelines on post-exposure prophylaxis for HIV and the use of co-trimoxazole prophylaxis for HIV-related infections among adults, adolescents and children: recommendations for a public health approach. In: December 2014 supplement to the 2013 consolidated guidelines on the use of antiretroviral drugs for treating and preventing HIV infection, vol. 42. World Health Organization, Geneva (2014).
29. Derrien M, Alvarez AS, de Vos WM. The gut microbiota in the first decade of life. Trends Microbiol (2019) 27(12):997–1010. doi: 10.1016/j.tim.2019.08.001
30. Chiappini E, Bianconi M, Dalzini A, Petrara MR, Galli L, Giaquinto C, et al. Accelerated aging in perinatally HIV-infected children: clinical manifestations and pathogenetic mechanisms. Aging (2018) 10(11):3610–25. doi: 10.18632/aging.101622
31. Dalzini A, Petrara MR, Ballin G, Zanchetta M, Giaquinto C, De Rossi A. Biological aging and immune senescence in children with perinatally acquired HIV. J Immunol Res (2020) 2020:8041616. doi: 10.1155/2020/8041616
32. Morreale C, Giaroni C, Baj A, Folgori L, Barcellini L, Dhami A, et al. Effects of perinatal antibiotic exposure and neonatal gut microbiota. Antibiot Basel Switz (2023) 12(2):258. doi: 10.3390/antibiotics12020258
33. Machiavelli A, Duarte RTD, Pires MM de S, Zárate-Bladés CR, Pinto AR. The impact of in utero HIV exposure on gut microbiota, inflammation, and microbial translocation. Gut Microbes (2019) 10(5):599–614. doi: 10.1080/19490976.2018.1560768
34. Lewis SJ, Heaton KW. Stool form scale as a useful guide to intestinal transit time. Scand J Gastroenterol (1997) 32(9):920–4. doi: 10.3109/00365529709011203
35. World health organization (WHO) (2023). Available at: https://www.who.int/es.
36. Pinto-Cardoso S, Lozupone C, Briceño O, Alva-Hernández S, Téllez N, Adriana A, et al. Fecal Bacterial Communities in treated HIV infected individuals on two antiretroviral regimens. Sci Rep 6 marzo (2017) 7:43741. doi: 10.1038/srep43741
37. Nutrition surveys and assessment (2023). Available at: https://www.nutrisurvey.de/.
38. McGuire S. U.S. Department of Agriculture and U.S. Department of Health and Human Services, Dietary Guidelines for Americans, 2010. 7th Edition, (Washington, DC: U.S. Government Printing Office). (2011) 2(3):293–4. doi: 10.3945/an.111.000430
39. Briceño O, Chávez-Torres M, Peralta-Prado A, Garrido-Rodríguez D, Romero-Mora K, Pinto-Cardoso S, et al. Associations between recent thymic emigrants and CD4+ T-cell recovery after short-term antiretroviral therapy initiation. AIDS Lond Engl (2020) 34(4):501–11. doi: 10.1097/QAD.0000000000002458
40. Monaco G, Chen H, Poidinger M, Chen J, de Magalhães JP, Larbi A. flowAI: automatic and interactive anomaly discerning tools for flow cytometry data. Bioinforma Oxf Engl (2016) 32(16):2473–80. doi: 10.1093/bioinformatics/btw191
41. Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res (2013) 41(1):e1. doi: 10.1093/nar/gks808
42. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol (2019) 37(8):852–7. doi: 10.1038/s41587-019-0209-9
43. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods (2016) 13(7):581–3. doi: 10.1038/nmeth.3869
44. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res (2013) 41(Database issue):D590–6. doi: 10.1093/nar/gks1219
45. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol (2011) 12(6):R60. doi: 10.1186/gb-2011-12-6-r60
46. Sessa L, Reddel S, Manno E, Quagliariello A, Cotugno N, Del Chierico F, et al. Distinct gut microbiota profile in antiretroviral therapy-treated perinatally HIV-infected patients associated with cardiac and inflammatory biomarkers. AIDS Lond Engl (2019) 33(6):1001–11. doi: 10.1097/QAD.0000000000002131
47. 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 Lond Engl (2015) 29(18):2409–18. doi: 10.1097/QAD.0000000000000869
48. Borgognone A, Noguera-Julian M, Oriol B, Noël-Romas L, Ruiz-Riol M, Guillén Y, et al. Gut microbiome signatures linked to HIV-1 reservoir size and viremia control. Microbiome (2022) 10(1):59. doi: 10.1186/s40168-022-01247-6
49. Effendi RMRA, Anshory M, Kalim H, Dwiyana RF, Suwarsa O, Pardo LM, et al. Akkermansia muciniphila and Faecalibacterium prausnitzii in Immune-Related Diseases. Microorganisms (2022) 10(12):2382. doi: 10.3390/microorganisms10122382
50. Si J, Kang H, You HJ, Ko G. Revisiting the role of Akkermansia muciniphila as a therapeutic bacterium. Gut Microbes 14(1):2078619. doi: 10.1080/19490976.2022.2078619
51. Fitzgerald FC, Lhomme E, Harris K, Kenny J, Doyle R, Kityo C, et al. Microbial translocation does not drive immune activation in Ugandan children infected with HIV. J Infect Dis (2019) 219(1):89–100. doi: 10.1093/infdis/jiy495
52. Weinberg A, Giganti MJ, Sirois PA, Montepiedra G, Canniff J, Agwu A, et al. Coordination of inflammatory responses in children with perinatally acquired HIV infection. AIDS Lond Engl (2022) 36(8):1117–27. doi: 10.1097/QAD.0000000000003229
53. Dirajlal-Fargo S, El-Kamari V, Weiner L, Shan L, Sattar A, Kulkarni M, et al. Altered intestinal permeability and fungal translocation in Ugandan children with human immunodeficiency virus. Clin Infect Dis Off Publ Infect Dis Soc Am (2020) 70(11):2413–22. doi: 10.1093/cid/ciz561
54. Dirajlal-Fargo S, Albar Z, Bowman E, Labbato D, Sattar A, Karungi C, et al. Increased monocyte and T-cell activation in treated HIV+ Ugandan children: associations with gut alteration and HIV factors. AIDS Lond Engl (2020) 34(7):1009–18. doi: 10.1097/QAD.0000000000002505
55. Noguera-Julian M, Rocafort M, Guillén Y, Rivera J, Casadellà M, Nowak P, et al. Gut microbiota linked to sexual preference and HIV infection. EBioMedicine (2016) 5:135–46. doi: 10.1016/j.ebiom.2016.01.032
56. Sperk M, van Domselaar R, Neogi U. Immune checkpoints as the immune system regulators and potential biomarkers in HIV-1 infection. Int J Mol Sci (2018) 19(7):2000. doi: 10.3390/ijms19072000
57. Foldi J, Kozhaya L, McCarty B, Mwamzuka M, Marshed F, Ilmet T, et al. HIV-infected children have elevated levels of PD-1+ Memory CD4 T cells with low proliferative capacity and high inflammatory cytokine effector functions. J Infect Dis (2017) 216(6):641–50. doi: 10.1093/infdis/jix341
58. Day CL, Kaufmann DE, Kiepiela P, Brown JA, Moodley ES, Reddy S, et al. PD-1 expression on HIV-specific T cells is associated with T-cell exhaustion and disease progression. Nature (2006) 443(7109):350–4. doi: 10.1038/nature05115
59. Klase Z, Ortiz A, Deleage C, Mudd JC, Quiñones M, Schwartzman E, et al. Dysbiotic bacteria translocate in progressive SIV infection. Mucosal Immunol (2015) 8(5):1009–20. doi: 10.1038/mi.2014.128
60. Ortiz AM, Brenchley J. Microbial translocation: translating SIV to HIV. Curr Opin HIV AIDS (2018) 13(1):15–21. doi: 10.1097/COH.0000000000000424
61. Bourke CD, Gough EK, Pimundu G, Shonhai A, Berejena C, Terry L, et al. Cotrimoxazole reduces systemic inflammation in HIV infection by altering the gut microbiome and immune activation. Sci Transl Med (2019) 11(486):eaav0537. doi: 10.1126/scitranslmed.aav0537
62. Gough EK, Bourke CD, Berejena C, Shonhai A, Bwakura-Dangarembizi M, Prendergast AJ, et al. Strain-level analysis of gut-resident pro-inflammatory viridans group Streptococci suppressed by long-term cotrimoxazole prophylaxis among HIV-positive children in Zimbabwe. Gut Microbes 11(4):1104–15. doi: 10.1080/19490976.2020.1717299
63. Bokulich NA, Chung J, Battaglia T, Henderson N, Jay M, Li H, et al. Antibiotics, birth mode, and diet shape microbiome maturation during early life. Sci Transl Med (2016) 8(343):343ra82. doi: 10.1126/scitranslmed.aad7121
64. Bush K, Bradford PA. β-lactams and β-lactamase inhibitors: an overview. Cold spring harb perspect med. (2016) 6(8):a025247. doi: 10.1101/cshperspect.a025247
65. Elvers KT, Wilson VJ, Hammond A, Duncan L, Huntley AL, Hay AD, et al. Antibiotic-induced changes in the human gut microbiota for the most commonly prescribed antibiotics in primary care in the UK: a systematic review. BMJ Open (2020) 10(9):e035677. doi: 10.1136/bmjopen-2019-035677
66. McDonnell L, Gilkes A, Ashworth M, Rowland V, Harries TH, Armstrong D, et al. Association between antibiotics and gut microbiome dysbiosis in children: systematic review and meta-analysis. Gut Microbes (2021) 13(1):1–18. doi: 10.1080/19490976.2020.1870402
67. Pérez-Tepayo S, Rodríguez-Ramírez S, Unar-Munguía M, Shamah-Levy T. Trends in the dietary patterns of Mexican adults by sociodemographic characteristics. Nutr J (2020) 19(1):51. doi: 10.1186/s12937-020-00568-2
68. Gentile CL, Weir TL. The gut microbiota at the intersection of diet and human health. Science (2018) 362(6416):776–80. doi: 10.1126/science.aau5812
69. Moszak M, Szulińska M, Bogdański P. You are what you eat-the relationship between diet, microbiota, and metabolic disorders-A review. Nutrients (2020) 12(4):1096. doi: 10.3390/nu12041096
70. Generoso M, Álvarez P, Kravietz A, Mwamzuka M, Marshed F, Ahmed A, et al. High soluble CD163 levels correlate with disease progression and inflammation in Kenyan children with perinatal HIV-infection. AIDS Lond Engl (2020) 34(1):33–8. doi: 10.1097/QAD.0000000000002378
71. Shamah-Levy T, Rivera-Dommarco J, Bertozzi S. [Mexico’s National Health and Nutrition Survey 2018-19: analysis of its main results]. Salud Publica Mex. (2020) 62(6):614–7. doi: 10.21149/12280
72. The Lancet Diabetes Endocrinology. Childhood obesity: a growing pandemic. Lancet Diabetes Endocrinol (2022) 10(1):1. doi: 10.1016/S2213-8587(21)00314-4
73. Pereira AR, Oliveira A. Dietary interventions to prevent childhood obesity: A literature review. Nutrients (2021) 13(10):3447. doi: 10.3390/nu13103447
74. Agyemang P, Kwofie EM, Baum JI. Transitioning to sustainable healthy diets: A model-based and conceptual system thinking approach to optimized sustainable diet concepts in the United States. Front Nutr (2022). doi: 10.3389/fnut.2022.874721
75. Asakura K, Mori S, Sasaki S, Nishiwaki Y. A school-based nutrition education program involving children and their guardians in Japan: facilitation of guardian-child communication and reduction of nutrition knowledge disparity. Nutr J (2021) 20(1):92. doi: 10.1186/s12937-021-00751-z
76. Serrano-Villar S, Talavera-Rodríguez A, Gosalbes MJ, Madrid N, Pérez-Molina JA, Elliott RJ, et al. Fecal microbiota transplantation in HIV: A pilot placebo-controlled study. Nat Commun (2021) 12(1):1139. doi: 10.1038/s41467-021-21472-1
77. Sainz T, Diaz L, Rojo D, Clemente MI, Barbas C, Gosalbes MJ, et al. Targeting the gut microbiota of vertically HIV-infected children to decrease inflammation and immunoactivation: A pilot clinical trial. Nutrients (2022) 14(5):992. doi: 10.3390/nu14050992
78. Goulder PJ, Lewin SR, Leitman EM. Paediatric HIV infection: the potential for cure. Nat Rev Immunol (2016) 16(4):259–71. doi: 10.1038/nri.2016.19
79. Maya-Lucas O, Murugesan S, Nirmalkar K, Alcaraz LD, Hoyo-Vadillo C, Pizano-Zárate ML, et al. The gut microbiome of Mexican children affected by obesity. Anaerobe (2019) 55:11–23. doi: 10.1016/j.anaerobe.2018.10.009
80. Atzeni A, Bastiaanssen TFS, Cryan JF, Tinahones FJ, Vioque J, Corella D, et al. Taxonomic and functional fecal microbiota signatures associated with insulin resistance in non-diabetic subjects with overweight/obesity within the frame of the PREDIMED-plus study. Front Endocrinol (2022) 13:804455. doi: 10.3389/fendo.2022.804455
81. Ou Y, Belzer C, Smidt H, de Weerth C. Development of the gut microbiota in the first 14 years of life and its relations to internalizing and externalizing difficulties and social anxiety during puberty. Eur Child Adolesc Psychiatry (2023). doi: 10.1007/s00787-023-02205-9
82. Valeri F, Endres K. How biological sex of the host shapes its gut microbiota. Front Neuroendocrinol (2021) 61:100912. doi: 10.1016/j.yfrne.2021.100912
83. Santinelli L, Ceccarelli G, Borrazzo C, Innocenti GP, Frasca F, Cavallari EN, et al. Sex-related differences in markers of immune activation in virologically suppressed HIV-infected patients. Biol Sex Differ (2020) 11(1):23. doi: 10.1186/s13293-020-00302-x
84. Krebs SJ, Slike BM, Sithinamsuwan P, Allen IE, Chalermchai T, Tipsuk S, et al. Sex differences in soluble markers vary before and after the initiation of antiretroviral therapy in chronically HIV-infected individuals. AIDS Lond Engl (2016) 30(10):1533–42. doi: 10.1097/QAD.0000000000001096
85. Wu GD, Chen J, Hoffmann C, Bittinger K, Chen YY, Keilbaugh SA, et al. Linking long-term dietary patterns with gut microbial enterotypes. Science (2011) 334(6052):105–8. doi: 10.1126/science.1208344
86. Dong TS, Guan M, Mayer EA, Stains J, Liu C, Vora P, et al. Obesity is associated with a distinct brain-gut microbiome signature that connects Prevotella and Bacteroides to the brain’s reward center. Gut Microbes (2022) 14(1):2051999. doi: 10.1080/19490976.2022.2051999
87. Mejía-León ME, Petrosino JF, Ajami NJ, Domínguez-Bello MG, de la Barca AMC. Fecal microbiota imbalance in Mexican children with type 1 diabetes. Sci Rep (2014) 4:3814. doi: 10.1038/srep03814
Keywords: pediatric HIV infection, mother-to-child-transmission, gut microbiome, immune activation, inflammation, antibiotics, healthy eating index
Citation: Rosel-Pech C, Pinto-Cardoso S, Chávez-Torres M, Montufar N, Osuna-Padilla I, Ávila-Ríos S, Reyes-Terán G, Aguirre-Alvarado C, Matías Juan NA, Pérez-Lorenzana H, Vázquez-Rosales JG and Bekker-Méndez VC (2023) Distinct fecal microbial signatures are linked to sex and chronic immune activation in pediatric HIV infection. Front. Immunol. 14:1244473. doi: 10.3389/fimmu.2023.1244473
Received: 22 June 2023; Accepted: 09 August 2023;
Published: 29 August 2023.
Edited by:
Nicholas Funderburg, The Ohio State University, United StatesReviewed by:
Tiffany Hensley-McBain, McLaughlin Research Institutes, United StatesModupe Coker, The State University of New Jersey, United States
Copyright © 2023 Rosel-Pech, Pinto-Cardoso, Chávez-Torres, Montufar, Osuna-Padilla, Ávila-Ríos, Reyes-Terán, Aguirre-Alvarado, Matías Juan, Pérez-Lorenzana, Vázquez-Rosales and Bekker-Méndez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sandra Pinto-Cardoso, sandra.pintocardoso.cieni@gmail.com; Vilma Carolina Bekker-Méndez, bekkermendez@yahoo.com
†These authors have contributed equally to this work and share first authorship