Abstract
The antiretroviral (ARV) drug combination of abacavir sulfate, dolutegravir, and lamivudine [ABC/DTG/3TC; Tri combination Anti-retroviral therapy (TC-ART)] has revolutionized HIV treatment by effectively targeting different stages of viral replication. Despite its therapeutic efficiency for maintaining low viremia in the mother during pregnancy, there are concerns for long-term liabilities in offspring that are indirectly exposed during vulnerable periods of development. The commensal microbiota plays a crucial role in maintaining overall gut health, and disruption of the microbiome is often linked to various extraintestinal effects such as immune dysregulation and inflammation. We recently reported the effects of this drug combination in altering fecal microbiome composition of aged rats perinatally exposed to ABC/DTG/3TC-ART. The fecal microbiome can provide only a snapshot of the composition of microbial community at the end of the digestive tract, which may not reflect the microbial population interacting with ileal mucosa. Thus, the current work reports the effects of this drug combination in the gut mucosa-associated microbiome of the same animals, which showed significant microbial diversity and species richness in high dose exposed female adult offspring, along with dose-dependent changes in Firmicutes/Bacteroidetes ratio. The high dose exposure also showed an increase in opportunistic bacterial species in male animals. Overall, we found that, similar to the fecal microbiome, perinatal exposure to TC-ART led to sex- and dose-dependent alterations in the gut mucosa-associated microbial population in aged rats, suggesting that early life exposure to these drugs may influence gut mucosa-associated immune responses and intestinal permeability.
Impact statement
This study shows that prenatal and lactational exposure to the Anti-retroviral therapy (TC-ART) drug combination (ABC/3TC/DTG) is associated with sex-dependent alterations in the ileal mucosa-associated microbiome composition when offspring are evaluated at 12Ā months of age. These findings emphasize the importance of including microbiome analyses in studies aiming at evaluating the long-term health effects of TC-ART exposure during early development.
Introduction
Human immunodeficiency virus (HIV) is a retrovirus that primarily targets the immune system, making it difficult for the body to fight against infections and diseases [1]. Without treatment, HIV can progress to AIDS (acquired immunodeficiency syndrome), which can be fatal [2]. Anti-retroviral (ARV) drugs have revolutionized HIV treatment by suppressing the replication of HIV, preventing it from progressing to AIDS. They have turned HIV from a fatal disease into a manageable chronic condition [2]. One such combination of drugs that has been used is abacavir sulphate, dolutegravir, and lamivudine [3]. Abacavir and lamivudine inhibit HIV reverse transcriptase and cause premature termination of viral DNA synthesis by incorporating non-functional nucleotide analogs into growing DNA chain [3, 4], while also having off target effects on the mitochondria [5, 6]. Millions of individuals globally are on Tri combination Anti-retroviral therapy (TC-ART), including pregnant women, to reduce the risk of mother to child transmission [7ā9].
The gut microbiome is important due to its beneficial functions and its involvement in maintaining gastrointestinal homeostasis, health, and disease progression [10]. ARV drugs particularly nucleoside reverse transcriptase inhibitors (NRTIs) (lamivudine), are known to affect host mitochondria by inhibiting mitochondrial DNA polymerase γ, leading to mitochondrial DNA depletion and dysfunction [11, 12]. Considering the evolutionary origin of mitochondria from symbiotic bacteria and similar biochemical pathways between mitochondria and prokaryotes, it is plausible that gut microbiome can be inadvertently affected by ARV drugs. However, most of the anti-retroviral drugs are not investigated for changes in gut bacterial community and overall gut health. Recent studies have shown that exposure to ARVās can have an impact on the hostās gastrointestinal health including altered microbiome and cytokine profiles [13, 14]. Microbial dysbiosis can have several downstream effects, including changes in immune function, intestinal mucosal permeability, and imbalance in other functional organs [15ā17]. For instance, dysbiosis in the ratio of Firmicutes and Bacteroidetes is associated with cardiac diseases [18], while increased Bacteroides levels are observed in liver ailments, with the degree of alteration correlating to the severity of the condition [19]. Similarly, decreases in butyrate producing bacteria have been associated with increased immune activation and microbial translocation [20], suggesting a relationship between microbial dysbiosis, and changes in cytokine profiles, and intestinal permeability.
The intestinal barrier is a complex multilayer system with an external physical āmucosal barrier,ā epithelial cells, and an internal immunological functional barrier in the form of Payerās patches [21]. The mucus layer on the top of the gut mucosa hosts several commensal bacteria and serves as a defense line along with the tight epithelial layer. Disruption in the barrier can compromise the integrity of the gut ecosystem, which is associated with several pathological conditions including lung, cardiac, neuronal, and autoimmune diseases [22, 23]. It is widely accepted that newborns are sterile at their birth, and they develop the microbiota through maternal exposure, surrounding environment, and through feeding practices [22]. Interestingly, the gut microbiome also plays a major role in infant immune system development, highlighting the importance of maternal gut health in shaping the offspring health [24, 25]. It is becoming increasingly evident that changes in the gut microbiome influence the development of various diseases in later life [26]. While most of the research is primarily focused on the direct effect of the drug on the host, there is limited information on the effect of the drug on the gut health of the progeny. We recently showed that perinatal exposure to this ARV combination altered the fecal microbiome of aged rats (12Ā months old) in a sex- and dose-dependent manner [27]. This study expanded those findings by exploring the changes in gut mucosal microbiota of the same animals, providing knowledge on the effects of drug exposure in the intestinal mucosa associated microbiome rather than focusing on the fecal -microbiome at the distal end of the intestinal tract.
Materials and methods
Animal housing, care, treatment, euthanasia, and sample collection
Time-mated Sprague Dawley rats (Hsd:SD) were obtained from Envigo (Indianapolis, IN). All animals (pregnant female rats and their male and female offspring) were housed in the animal facility at AmplifyBio, West Jefferson, OH, an independent, scientific contract research organization. Rats were approximately 11ā14Ā weeks of age upon receipt. The facilityās Institutional Animal Care and Use Committee (IACUC) reviewed the protocol and approved it. The IACUC number for this protocol is T06055. Animals were housed in polycarbonate cages with irradiated hardwood bedding chips (Sani ChipsĀ®; Envigo, Madison, WI). During gestation and lactation, rats were provided natural crinkled kraft paper for enrichment (Crink-lānestā¢, The Andersons; Maumee, Ohio). Offspring remained with their respective dams until postnatal day (PND) 21. After weaning, F1 offspring were provided polycarbonate rectangular shelters (Rat Retreatsā¢, Bio-Serve; Flemington, NJ) as enrichment and were group housed by sex, up to 5 per cage. During gestation and lactation, animals were fed irradiated NIH-07 pellets or wafers (Zeigler Bros., Gardners, PA). After weaning, animals were fed NTP-2000 (Zeigler Bros., Gardners, PA). Rats were provided municipal water ad libitum from an automatic watering system; both water and feed were analyzed for known contaminants that could interfere with or affect the outcome of the study, and none were found. Animals used in this investigation of the mucosa-associated microbiome were from a larger toxicology study that will be reported elsewhere (manuscript under preparation). The experimental design is outlined in Figure 1.
FIGURE 1

This figure illustrates the experimental flow for this study, which includes three treatment groups: control, low dose, and high dose of tricombo drug treatment. Pregnant rats were treated until the rat pups were weaned on postnatal day 21. DNA, was collected from the ileal tissue of 12-month-old rat offspring for microbiome analysis.
HIV dose selection and exposure to animals
There were 3 groups in this study (n = 5 per group). Group 1 (control): 0Ā mg/kg/day dosed with vehicle; Group 2 (low dose): 150/12.5/75Ā mg/kg ABC/DTG/3TC; Group 3 (high dose): 300/25/150Ā mg/kg ABC/DTG/3TC. Vehicle was 0.2% methylcellulose/0.1% Tween 80. The treatment regimen and administered doses were based on our previously published study, which used the same animal cohort and provided detailed justification for dose selection [27]. In summary, the doses were chosen based on human therapeutic relevance and previously established toxicological data.
TC-ART was administered via oral gavage to time-mated female Sprague Dawley rats (Hsd:SD) starting on gestational day 6 (GD6) and continued through gestation, and lactation (PND21) (NTP study T06055). F1 offspring were never directly dosed in this study. Dams were weighed approximately every 3Ā days, and dosing volumes (5Ā mL/kg) were based on the most recent body weight. One male and one female offspring from each dam were used for the study.
Collection of ileal tissue for DNA extraction
After the pups were weaned from the dams, they were housed up to five per sex per cage and provided food and water ad libitum. Once the offspring reached 12Ā months of age, which is comparable to approximately 26Ā years in humans [28], they were euthanized by exposure to carbon dioxide, strictly adhering to the guidelines set by the IACUC. In a sterile environment, the ileum was exposed and flushed with sterile PBS to remove feces. The samples were immediately placed in pre-labelled tubes, frozen in liquid nitrogen and transported overnight on dry ice to the National Center for Toxicological Research (NCTR), Jefferson, AR, for analysis.
Ileal mucosa-associated microbial DNA extraction for population analysis
Maintaining the sterile conditions, a small section of the ileal tissue was placed in CTAB extraction solution and finely minced using a scalpel. The minced tissue was transferred into bead-beating tubes (lysing matrix B) and processed using a Fast prep Machine at 6 m/s for 45 s. Thus, obtained lysate was incubated with proteinase k (at 65°C for 20 min), followed by RNase-A (at 37°C for 15 min) as described earlier [29]. To this an equal volume of phenol-chloroform-isopropanol (PCI) solution was added and centrifuged at 12,000xg. The supernatant aqueous solution was separated, and the DNA was precipitated using sodium acetate, isopropanol, and a polyacryl carrier. The DNA pellet was washed with 70% ethanol and then dissolved in nuclease- free water. The DNA was quantified with a Qubit fluorometer (Thermo Fisher Scientific, Waltham, Massachusetts) and subsequently used for sequencing.
Rat mucosa-associated microbial population amplicon sequencing
The samples were processed and standardized using a modified methodology based on Gohl et al. 2016 [30]. The KAPA HiFi Polymerase from KAPA Biosystems was used to amplify the variable region 4 (515f/806r) of the 16S rRNA gene. The libraries were subjected to sequencing on a MiSeq (Illumina) instrument, with paired-end 2 Ć 250 reads, using the MiSeq Reagent Kit V3 (Illumina, 600 cycle kit), aiming to achieve a target depth of 50,000 reads. The sequencing reads underwent adapter and primer removal using Cutadapt, and any reads with a Q score below 30 were discarded. The amplicon sequence variations (ASVs) were generated using dada2 v1.16.0, using quality-controlled reads. As part of the dada2 workflow, the paired FASTQ reads underwent trimming and filtering to eliminate reads that had ambiguous nucleotides (represented by āNās) or more than 2 sequencing mistakes per read. The dada2 learn error rate model was utilized to estimate the error profile prior to employing the core dada2 algorithm for inferring the sample composition. The forward and reverse readings were combined, and any chimeras were eliminated prior to taxonomy classification. The ASV taxonomy was determined up to the level of species using the SILVA v.138 database, following the technique developed earlier [31]. A minimum bootstrapping support of 50% was required for taxonomic assignment. Species-level taxonomy was awarded to ASVs based only on 100% identity and unambiguous reference matching.
Statistical/data analysis
Microbiome data was visualized using either microbiome analyst [32] or metaboanalyst websites [33]. For the microbiome analysis, the software used an ANOVA, followed by post hoc pairwise comparisons using Mann-Whitney and Kruskal-Wallis statistical tests. For the comparison between two groups either EdgeR or student T-tests were conducted. The samples for the microbiome analyst were normalized to remove 19 low abundance and 11 low variance values based on mean and standard deviation respectively. While for the metaboanalyst the samples normalization was performed using the normalization by sum. Following that the data scaling was applied using mean centering and log transformation options to ensure approximate normal distribution.
Results
Comparative analysis of the treatment groups (n = 5 per group) was conducted to assess the impact of gestational and lactational exposure of TC-ART on microbial profiles, focusing on sex-specific impacts on the offspring. During the analysis one control female and one high dose female animal showed deviation in the microbial profile compared to the other animals in respective experimental groups. However, the outlier test showed the two animals are not statistically classified outliers. We have excluded the animals for the alpha diversity but to maintain the integrity, we have included the animals in our dataset for the rest of the analysis. The inclusion of these animals for later analysis did not alter our conclusions, rather it highlights the natural variability within the controlled experimental groups.
Sex-specific microbiota differences and differential response to exposure are observed in male and female animals
Sex-specific differences were observed between the microbiota of the control animals. A total of 37 significant features were differentiated between the female and male animals (Supplementary Table S1). At the phyla and species level 2 significant features were identified, 3 at family level, and 32 features varied at the genus level. The cladogram analysis revealed notable differences between the male and the female animals at the genus level (Figure 2). Notably the Bacteroidota and Bifidobacteriales were significantly increased in the female animals compared to the male animals while the Clostridia, Bacilli, and Eubacterium coprostanoligenes group were decreased in the female animals. The male and the female animals have responded differently to TC-ART exposure. While the exposure in male animals resulted in increase of opportunistic bacteria the female animals showed decrease in the abundance of opportunistic bacterial species. To have a more granular view on the differences we analyzed the changes in the microbiome upon exposure at different taxonomic levels.
FIGURE 2

The cladogram represents the phylogenetic relationship between the male and the female animals between different taxa. The blue or highlighted branches represent significant decreased taxa in females while the red represent the significantly increased taxa. The branches are colored based on median log2 ratios and the intensity of color represents the abundance. Significant differences were observed in Bacteroidales and Bacili.
Mucosa associated microbiome diversity changes with 16s amplicon sequencing data
The ileal mucosa-associated bacterial population in the adult offspring perinatally exposed to TC-ART was assessed by 16s amplicon sequencing.
The changes in bacterial diversity, richness, and evenness in each group were evaluated by using different alpha diversity indices including total number of observed operational taxonomic units (OTUs), Chao1, Shannon, and Simpson as shown in Figure 3. The female high dose offspring exhibited significantly higher species richness compared to all other groups, with an observed OTU of 132 (Figure 3a); results that were further confirmed when evaluating the Shannon and Chao1 diversity indices (Figures 3b, c). Interestingly, female control animals showed higher species richness compared to male animals, indicating base line sex differences. The microbial species richness in the high dose exposed female offspring were significantly higher compared to other groups (Figures 3aāc; Supplementary Table S2), suggesting that maternal administration of the drugs at the high dose resulted in the enhanced species richness and rare species (low abundant species) in these animals. Conversely, no significant differences between groups were observed using the Simpson Index, which measures richness and evenness between bacterial communities. This suggests a similar level of species dominance within bacterial community among all tested animals (Figure 3d). The rarefaction curves also suggest higher species richness at higher sequence depths in both male and female animals in the high dose-exposed groups (Supplementary Figure S1).
FIGURE 3

Alpha diversity (a) Observed OTU abundance over sequence number. The figure shows the observed operational taxonomic unit (OTU) abundance as a function of sequence number for male and female animals across three experimental groups: control, low dose, and high dose treatment. Each line represents the mean OTU abundance for each group. The high dose exposed female animals display a higher OTU abundance compared to the other groups. (b) Shannon Diversity Index. (c) Chao1 Index (d) Simpson Diversity Index. Illustrates the diversity of microbial communities in male and female animals across the groups. The high dose treated female animals show significant differentiation from the other groups, indicating a notable shift in microbial diversity and richness associated with higher treatment levels. The data represents mean with standard deviation (n = 5; *p < 0.05).
Perinatal exposure to TC-ART leads to taxonomic level changes at phyla, genus, and species levels
The principal component analysis (PCA) at the phylum level showed a distinct shift in the clustering of the high dose exposed female animals (Figure 4a), correlating with significant alterations in Firmicutes/Bacteroidetes ratio (Figure 4b). The low dose showed a substantial increase of Bacteroidetes relative microbial abundance compared to the control animals, while the high dose showed an even greater increase. The Bacteroidetes change was consistent across all the female high dose exposed animals (Figure 4b). The relative abundance of Firmicutes was reduced with the treatment. Although the males also showed a similar trend, the variations were not as prominent as in female animals (Figure 4b; Supplementary Figure S2). Treated male animals had increased Proteobacteria phyla compared to the control animals, while females had a decrease (Figure 4b). Similarly, an increase in the phyla Deferribacterota was observed in the high dose exposed female rats relative to the control animals. Of note, males from the high dose-only and females from both dosed-groups have more unidentified phyla, indicating a diversified bacterial population associated with drug exposure. This, in turn, could potentially lead to metabolic disturbances.
FIGURE 4

Principal Component Analysis (PCA) of microbiome profiles and changes in gut microbiota composition in male and female rat pups exposed to Tricombo drug. (a) PCA of the microbial profiles. The PCA score plot illustrates the clustering of microbiome profiles from each group. The plot shows ellipses indicating the 95% confidence intervals with the high dose female (blue) distinctly separated from the other groups. (b) Relative abundance of gut microbiota at the phyla level: The bar diagram represents the relative abundance of different bacterial phyla in male and female rat pups exposed to low and high doses of the tricombo drug compared to control animals. The bacterial phyla include Firmicutes, Proteobacteria, Actinobacteriota, Bacteroidota, Cyanobacteria, Patescibacteria, Verrucomicrobiota, Desulfobacterota, and Deferribacterota. The High dose exposed female animals has shown significant variations in Firmicutes/Bacteroidota ratio compared to other groups. Data is represented as the percentage of total bacterial abundance for each group Control Male, Low dose Male, High dose Male, Control Female, Low dose Female, High dose Female.
The identified changes in bacterial diversity described above might differentially influence downstream effects. To determine if this was the case, we next analyzed bacteria composition at the genus and species level, which can provide information on their pathogenic or commensal origins. We observed notable sex- and treatment- dependent variabilities. The major genera showed comparable differences between the male and female control animals (Figures 5a, b). Beneficial genera such as Lactobacillus, Lachnospiraceae, Lactococcus, Romboutsia, and Bifidobacterium, that promote gut health showed dose-dependent alterations (Figures 5a, b). Interestingly, Lactobacillus genera has shown a dose dependent decrease in all groups while Romboutsia from Firmicutes phyla showed a decrease in male animals (Figure 5a). Conversely, the genera associated with inflammation or pathogenicity, including Escherichia (Escherichia schigella), Pseudomonas, Angelakisella, and Parabacteroides, were increased in the groups exposed to TC-ART. Notably, Pseudomonas, which is linked to gut dysbiosis and immune activation, showed a dose-dependent increase in both male and female subjects. Clostridium sensu stricto1 was decreased in females in both exposed groups but not changed in male groups (Figures 5a, b; Supplementary Figure S3a).
FIGURE 5

Relative abundance of bacterial genus in different TC-ART exposure groups. (a) The stacked bar plot illustrates the composition of gut microbiota in male and female subjects across groups: Control Male, Low dose Male, High dose Male, Control Female, Low dose Female, and High dose Female. Each color in the bars represents a different bacterial genus. Lactobacillus (yellow), Romboutsia (green), and Turicibacter (reddish pink), Lactococcus (orange), Clostridium sensu stricto 1 (grey) show significant alterations. High dose treatments decreased Lactobacillus and Lactococcus with more pronounced changes in females. These variations highlight the impact of exposure on gut microbiota composition. (b) Heatmap of bacterial genera distribution across different exposure doses and in control animals. The rows correspond to bacterial genera, and the columns to the samples within each group. The color scale indicates the abundance of each genus, ranging from blue (low abundance) to red (high abundance). A distinct microbial profile observed in the high dose exposed female animals compared to other groups, indicating a significant shift in bacterial genera due to the exposure. (c) The cladogram represents comparison between the control female and the high dose female animals. The blue highlighted branches represent significant taxa and the taxa with higher median log2 ratios in high dose female group are more abundant compared to the control animals indicated by color gradient and intensity. (d) The pie chart of top 20 abundance genera with the opportunistic bacterial genera highlighted between the groups (dark green). The abundance of the opportunistic bacteria was added and represented in the pie diagram. The exposed male animals have shown increased abundance of opportunistic bacterial genera while the females have shown low abundance with the exposure to TC-ART.
Genera such as Cupriavidus, Odoribacter, and Tyzzerella were uniquely present in both male and female high dose exposed groups. The abundance of some genera like Bacteroides, Butyricoccus, Turicibacter, and Eubacterium (Eubacterium ventriosum, Eubacterium siraeum) showed transient changes in response to treatment indicating sensitivity of the bacterial community following prenatal exposure to the ART (Figure 5b).
When single-factor statistical comparisons were performed using EdgeR, a total of 19 significant features were altered at the genus level between the control female and high dose female animals using PCA analysis (Supplementary Figure S3b). Additionally, 15 microbial features differed, including 9 variations at the genus level, were identified between the control male and the female animals, consistent with the notion of base line sex-specific differences described above. The cladogram (Figure 5c) depicts significant alterations in the relative abundance of several taxa such as Oscillospiraceae, Intestinimonas, Acettatifactor, Rikenella, Anaeroplasma, Acholeplasmatales, and Clostridia VadinBB60 group, are significantly enriched in high dose female group. The shift in Oscillospiraceae and Intestinimonas families in the high dose-exposed animals was notable. Also, females showed higher abundance of Bifidobacterium while males showed increased levels of Clostridia. Changes were observed in genera within the Bacteroidetes and Bacilli with more opportunistic genera in the control female animals compared to the male animals, which were further changed in the exposed groups (Figure 5d).
To have a more granular view at the taxonomy, we investigated abundance of bacterial species in each experimental group, which revealed dose- and sex-dependent changes (Figure 6a). Protective bacteria like Lactobacillus intestinalis and Akkermansia muciniphila decreased while other commensal bacteria Bacteroides uniformis, Bacteroides caccae, and Mucispirillum schaedleri were increased in treatment groups (Figures 6a, b). Interestingly, most of the commensal bacteria showed lower abundance in the low dose exposed animals compared to the high dose exposed animals.
FIGURE 6

Comparative analysis of mucosa associated microbial species across different groups. Bubble plots representing the relative abundance of (a) opportunistic bacterial species in mucosa tissue (b) commensal bacterial species across the treatment groups in mucosa tissue.
Discussion
The ability of TC-ART to cross the placental barrier and its presence in milk invite a better understanding of the risks associated with their use in pregnancy and lactation [34]. Numerous studies have explored the effects of TC-ART on the mother and the offspring, focusing on aspects like cardiovascular diseases, neuronal development, immune health [35]. These studies focus on the effects of the drugs on these tissues, but often neglect to include analysis on the effects on the gut health of the offspring. Given the importance of gut health and the microbiome in maintaining individual homeostasis, and their association with several pathological conditions when disrupted, it is important to consider how they may be affected by early drug exposure. Perinatal TC-ART exposure may alter the gut microbiome that could impact metabolism, intestinal permeability, and immune system contributing to long-term health outcomes. In this study, by investigating mucosa-associated microbial changes in rats that were exposed to TC-ART solely through placental and lactational transfer, we found significant changes in the composition of intestinal microbiome in adulthood.
Our data aligns with the understanding that the fecal microbiome is more abundant and diverse than the mucosa-associated microbiome. Although both fecal and mucosal microbiomes showed differences, the changes were more pronounced in the mucosal microbiome. This may suggest that the mucosal environment is more sensitive to high-dose exposure and exhibits differential responses, possibly due to the drug being primarily absorbed in the small intestine [27, 36]. Furthermore, several studies have shown that the mucosal environment is crucial for microbiome-immune interactions due to its proximity with the intestinal epithelium and pattern recognition receptors (PRRs) [37, 38]. The combined effect of mucosal environment along with the presence of TC-ART could have significantly contributed to observed shifts in mucosal microbiome. Interestingly, the mucosal data showing a higher proportion of opportunistic bacteria and decreased commensal bacteria, along with the differences in the fecal and mucosal microbiome of the control animals, highlights the importance of studying both to understand the dynamics of drug exposure.
Overall, our data showed sex-dependent and dose-dependent drug-induced changes in the gut microbiome, which were more pronounced in females than males. Notably, sex-dependent microbial variations between male and female control animals were identified, which may provide insights for understanding the sex-specific microbial functions or outcomes. Although female animals exhibited lower microbial abundance, their higher observed OTU counts compared to male samples suggests greater microbial diversity in females. Additionally, treatment in both male and females led to reduced microbial abundance indicating a non-favorable effect of the exposure on the gut microbiome. The opportunistic genera in the male animals increased with the exposure while it decreased in the female animals, suggesting sex-based differences in the effects of the drugs. Sex hormones such as estrogen and testosterone are known to influence gut microbial composition, diversity, and functions [39]. In healthy women, higher estrogen levels are associated with increased Bacteroidetes and microbial diversity, while in healthy men, elevated testosterone levels correlate with higher levels of Ruminococcus and Acinetobacter. However, in women, higher testosterone levels are associated with Escherichia and Shigella species, while Ruminococcus species; is negatively associated with elevated testosterone levels [40]. The increased Bacteroidetes and reduced Firmicutes ratio in high dose females observed in our study is often associated with several pathologies including IBS and immune profile alterations [41]. Generally, decreases in Firmicutes and increases in Bacteriodata are considered ideal but high changes in the composition might also disrupt the digestive function and other metabolic pathways. Some studies have revealed that similar microbial changes observed in female animals are associated with disease outcomes in patients with Alzheimerās, refractory epilepsy, and some neurological diseases [42ā44]. Alterations in Actinobacteria and Proteobacteria in the treated groups suggest microbial dysbiosis, which has been linked to inflammatory diseases [41], coronary heart failure, and other cardiac abnormalities [45ā47]. For example, Romboutsia ilealis species is often linked with better gut health, including improved metabolism, lowering dietary inflammation, and controlling blood sugar and fat [48, 49]. The enrichment of Romboutsia in the female animals could contribute to decreased proinflammatory cytokines production. Similarly, the changes in abundance of the Clostridium sensu stricto1, Escherichia, shigella, and Enterococcus observed in the treated groups have been associated with metabolic dysfunction and fatty acid liver disease [50, 51], as well as heart failure and coronary artery diseases [47]. [52]. Notably, these animals do show cardiac dysfunction later in life (Taube et al., in preparation), but whether the changes in the microbiome reported here play a role in that phenotype remain unclear.
The reduced abundance of Akkermansia mucinophila is often associated with vascular inflammation, endotoxemia, and barrier dysfunction [45, 53]. Akkermansia mucinophila is a mucin degrading bacterial community that modulates expression of tight junction genes like occludins [53]; thus can influence the intestinal permeability and induce translocation of pathogenic bacteria into the systemic circulation that could contribute the metabolic diseases. Another species, Faecalibaculum rodentium, known for its cholesterol degrading activity and protective effects against intestinal tumor growth, was slightly increased with the treatment in both male and female animals. This species is also known to induce oxidative stress and inflammation [54, 55]. The differential abundance of Aslercreutzia eqyolifaciens species is related to several disease conditions such as metabolic liver diseases [56]. Given that most of the bacteria are initially acquired from the mother and the surrounding environment, the coprophagy during the co-housing of dam and pups during early development could be another possibility of the transfer of altered microbial population (due to impact of drug on the intestinal microbiota) from dams to pups [57, 58]. These microbial changes may persist into adulthood influencing the microbiome composition of the offspring.
In conclusion, our study shows the effects of TC-ART on the mucosal microbiome of the adult offspring of animals exposed through gestation and lactation. The differences observed highlight the potential sensitivity of the mucosal environment to the drug exposure. The extent to which these changes might affect other phenotypes on the animals, including impacts on their cardiovascular health, is unclear but deserves further investigation. Since our data also suggest sex-specific differences, further studies on the influence of sex hormones to the phenotypes here described are also warranted.
Statements
Data Availability Statement
The datasets presented in this article are not readily available because the work is conducted at the government institute. Requests to access the datasets should be directed to the corresponding author.
Ethics Statement
All animals (pregnant female rats and their male and female offspring) were housed in the animal facility at AmplifyBio, West Jefferson, OH, an independent, scientific contract research organization. Rats were approximately 11-14 weeks of age upon receipt. The facilityās Institutional Animal Care and Use Committee (IACUC) reviewed the protocol and approved it. The IACUC number for this protocol is T06055. The study was conducted in accordance with the local legislation and institutional requirements.
Author Contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. YY, CM, and KK were supported by an appointment to the Postgraduate Research Program at NCTR administered by the Oak Ridge Institute for Science and Education (ORISE) through an interagency agreement with the U.S. Department of Energy and the U.S. Food and Drug Administration. This study was supported by NIEHS under an interagency agreement between FDA and NIEHS (FDA IAA#224-17-0502 and NIEHS IAA AES23002001/GT&C A2310-075-075-052689). It was also supported by contracts HHSN273201400015C and HHSN273201400022C. VS and HC were supported [in part] by the Intramural Research Program of the National Institutes of Health, NIEHS, under ZIA ES103379.
Acknowledgments
The authors would like to thank Dr. Radwa Hanafy (NCTR), Dr. Jyotshnabala Kanungo (NCTR), and Dr. Kelly Shipkowski (NIEHS) for reviewing the manuscript and providing valuable comments and suggestions.
Conflict of Interest
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Author Disclaimer
This article reflects the views of the authors and does not necessarily reflect those of the U.S. Food and Drug Administration. Any mention of a commercial product is for clarification only, it is not an endorsement for the use of it.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.ebm-journal.org/articles/10.3389/ebm.2025.10564/full#supplementary-material
References
1
Douek DC Roederer M Koup RA . Emerging concepts in the immunopathogenesis of AIDS. Annu Rev Med (2009) 60(1):471ā84. 10.1146/annurev.med.60.041807.123549
2
Deeks SG Lewin SR Havlir DV . The end of AIDS: HIV infection as a chronic disease. The lancet (2013) 382(9903):1525ā33. 10.1016/s0140-6736(13)61809-7
3
Arena CT . Triumeq (dolutegravir/abacavir/lamivudine) for the Treatment of HIV (2015). Available online at: https://www.clinicaltrialsarena.com/projects/triumeq-dolutegravirabacavirlamivudine-for-the-treatment-of-hiv/?cf-view. (Accessed April 9, 2024)
4
Amblard F Patel D Michailidis E Coats SJ Kasthuri M Biteau N et al HIV nucleoside reverse transcriptase inhibitors. Eur J Med Chem (2022) 240:114554. 10.1016/j.ejmech.2022.114554
5
Sun Y Xu M Duan Q Bryant JL Xu X . The role of autophagy in the progression of HIV infected cardiomyopathy. Front Cell Dev Biol (2024) 12:1372573. 10.3389/fcell.2024.1372573
6
Balogun K Serghides L . Comparison of the effects of three dual-nucleos(t)ide reverse transcriptase inhibitor backbones on placenta mitochondria toxicity and oxidative stress using a mouse pregnancy model. Pharmaceutics (2022) 14(5):1063. 10.3390/pharmaceutics14051063
7
Patel K Huo Y Jao J Powis KM Williams PL Kacanek D et al Dolutegravir in pregnancy as compared with current HIV regimens in the United States. New Engl J Med (2022) 387(9):799ā809. 10.1056/nejmoa2200600
8
Drugs.com. Abacavir/dolutegravir/lamivudine pregnancy and breastfeeding warnings. (2024).
9
Siegfried N van der Merwe L Brocklehurst P Sint TT . Antiretrovirals for reducing the risk of motherātoāchild transmission of HIV infection. The Cochrane database Syst Rev (2011)(7) CD003510. 10.1002/14651858.cd003510.pub3
10
Guinane CM Cotter PD . Role of the gut microbiota in health and chronic gastrointestinal disease: understanding a hidden metabolic organ. Therap Adv Gastroenterol (2013) 6(4):295ā308. 10.1177/1756283x13482996
11
Cossarizza A Moyle G . Antiretroviral nucleoside and nucleotide analogues and mitochondria. Aids (2004) 18(2):137ā51. 10.1097/00002030-200401230-00002
12
CĆ“tĆ© HC . Mechanisms of antiretroviral therapy-induced mitochondrial dysfunction. Curr Opin HIV AIDS (2007) 2(4):253ā60. 10.1097/coh.0b013e3281df3410
13
Ray S Narayanan A Giske CG Neogi U Sƶnnerborg A Nowak P . Altered gut microbiome under antiretroviral therapy: impact of efavirenz and zidovudine. ACS Infect Dis (2020) 7(5):1104ā15. 10.1021/acsinfecdis.0c00536
14
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. 10.1097/COH.0000000000000428
15
Brown K DeCoffe D Molcan E Gibson DL . Diet-induced dysbiosis of the intestinal microbiota and the effects on immunity and disease. Nutrients (2012) 4(8):1095ā119. 10.3390/nu4081095
16
StƤrkel P Leclercq S de Timary P Schnabl B . Intestinal dysbiosis and permeability: the yin and yang in alcohol dependence and alcoholic liver disease. Clin Sci (2018) 132(2):199ā212. 10.1042/cs20171055
17
Gong X Ma Y Deng X Li A Li X Kong X et al Intestinal dysbiosis exacerbates susceptibility to the anti-NMDA receptor encephalitis-like phenotype by changing blood brain barrier permeability and immune homeostasis. Brain Behav Immun (2024) 116:34ā51. 10.1016/j.bbi.2023.11.030
18
Yin J Liao SX He Y Wang S Xia GH Liu FT et al Dysbiosis of gut microbiota with reduced trimethylamine-N-oxide level in patients with large-artery atherosclerotic stroke or transient ischemic attack. J Am Heart Assoc (2015) 4(11):e002699. 10.1161/jaha.115.002699
19
Boursier J Mueller O Barret M Machado M Fizanne L AraujoāPerez F et al The severity of nonalcoholic fatty liver disease is associated with gut dysbiosis and shift in the metabolic function of the gut microbiota. Hepatology (2016) 63(3):764ā75. 10.1002/hep.28356
20
Kumari R Ahuja V Paul J . Fluctuations in butyrate-producing bacteria in ulcerative colitis patients of North India. World J Gastroenterol (2013) 19(22):3404ā14. 10.3748/wjg.v19.i22.3404
21
Seo K Seo J Yeun J Choi H Kim Y-I Chang S-Y . The role of mucosal barriers in human gut health. Arch Pharmacal Res (2021) 44(4):325ā41. 10.1007/s12272-021-01327-5
22
Leech B Schloss J Steel A . Association between increased intestinal permeability and disease: a systematic review. Adv Integr Med (2019) 6(1):23ā34. 10.1016/j.aimed.2018.08.003
23
Bischoff SC Barbara G Buurman W Ockhuizen T Schulzke JD Serino M et al Intestinal permeability--a new target for disease prevention and therapy. BMC Gastroenterol (2014) 14:189. 10.1186/s12876-014-0189-7
24
Alwazzan A . The influence of gut microbiome derived neurotransmitters on neonatal immune response. Pakistan BioMedical J (2024) 1. 10.54393/pbmj.v7i04.1075
25
Willers M Ulas T Vƶllger L Vogl T Heinemann AS Pirr S et al S100A8 and S100A9 are important for postnatal development of gut microbiota and immune system in mice and infants. Gastroenterology (2020) 159(6):2130ā45.e5. 10.1053/j.gastro.2020.08.019
26
Arrieta M-C Stiemsma LT Amenyogbe N Brown EM Finlay B . The intestinal microbiome in early life: health and disease. Front Immunol (2014) 5:427. 10.3389/fimmu.2014.00427
27
Muthumula CMR Yanamadala Y Gokulan K Karn K Cunny H Sutherland V et al Effect of in utero and lactational exposure to antiretroviral therapy on the gut microbial composition and metabolic function in aged rat offspring. Exp Biol Med (Maywood) (2025) 250:10468. 10.3389/ebm.2025.10468
28
Sengupta P . The laboratory rat: relating its age with human's. Int J Prev Med (2013) 4(6):624ā30.
29
Khare S Ficht TA Santos RL Romano J Ficht AR Zhang S et al Rapid and sensitive detection of Mycobacterium avium subsp. paratuberculosis in bovine milk and feces by a combination of immunomagnetic bead separation-conventional PCR and real-time PCR. J Clin Microbiol (2004) 42(3):1075ā81. 10.1128/jcm.42.3.1075-1081.2004
30
Gohl DM Vangay P Garbe J MacLean A Hauge A Becker A et al Systematic improvement of amplicon marker gene methods for increased accuracy in microbiome studies. Nat Biotechnol (2016) 34(9):942ā9. 10.1038/nbt.3601
31
Wang Q Garrity GM Tiedje JM Cole JR . Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol (2007) 73(16):5261ā7. 10.1128/aem.00062-07
32
MicrobiomeAnalyst. MicrobiomeAnalyst. (2024).
33
MetaboAnalyst. MetaboAnalyst. (2024).
34
Kobbe R Schalkwijk S Dunay G Eberhard JM Schulze-Sturm U Hollwitz B et al Dolutegravir in breast milk and maternal and infant plasma during breastfeeding. AIDS (2016) 30(17):2731ā3. 10.1097/qad.0000000000001259
35
Vannappagari V Thorne C . Pregnancy and neonatal outcomes following prenatal exposure to dolutegravir. JAIDS J Acquired Immune Deficiency Syndromes (2019) 81(4):371ā8. 10.1097/qai.0000000000002035
36
Healthcare V . Dolutegravir (DTG) + abacavir (ABC) + lamivudine (3TC) HIV/HCV medications guide: HIV/HCV medications guide (2022). Available online at: https://www.hivmedicationguide.com/medicament/dolutegravir-dtg-abacavir-abc-lamivudine-3tc/. (Accessed April 9, 2024)
37
Tian K Jing D Lan J Lv M Wang T . Commensal microbiome and gastrointestinal mucosal immunity: harmony and conflict with our closest neighbor. Immun Inflamm Dis (2024) 12(7):e1316. 10.1002/iid3.1316
38
Fukata M Arditi M . The role of pattern recognition receptors in intestinal inflammation. Mucosal Immunol (2013) 6(3):451ā63. 10.1038/mi.2013.13
39
Lahiani M Gokulan K Sutherland V Cunny HC Cerniglia CE Khare S . Early developmental exposure to triclosan impacts fecal microbial populations, IgA and functional activities of the rat microbiome. J Xenobiotics (2024) 14(1):193ā213. 10.3390/jox14010012
40
dāAfflitto M Upadhyaya A Green A Peiris M . Association between sex hormone levels and gut microbiota composition and diversity-A systematic review. J Clin Gastroenterol (2022) 56(5):384ā92. 10.1097/mcg.0000000000001676
41
Pittayanon R Lau JT Yuan Y Leontiadis GI Tse F Surette M et al Gut microbiota in patients with irritable bowel syndromeāa systematic review. Gastroenterology (2019) 157(1):97ā108. 10.1053/j.gastro.2019.03.049
42
Zhou K Jia L Mao Z Si P Sun C Qu Z et al Integrated macrogenomics and metabolomics explore alterations and correlation between gut microbiota and serum metabolites in adult epileptic patients: a pilot study. Microorganisms (2023) 11(11):2628. 10.3390/microorganisms11112628
43
Li Q Gu Y Liang J Yang Z Qin J . A long journey to treat epilepsy with the gut microbiota. Front Cell Neurosci (2024) 18:1386205. 10.3389/fncel.2024.1386205
44
Hung CC Chang CC Huang CW Nouchi R Cheng CH . Gut microbiota in patients with Alzheimer's disease spectrum: a systematic review and meta-analysis. Aging (Albany NY) (2022) 14(1):477ā96. 10.18632/aging.203826
45
Liu L He X Feng Y . Coronary heart disease and intestinal microbiota. Coron Artery Dis (2019) 30(5):384ā9. 10.1097/mca.0000000000000758
46
Khan I Khan I Jianye Z Xiaohua Z Khan M Hilal MG et al Exploring blood microbial communities and their influence on human cardiovascular disease. J Clin Lab Anal (2022) 36(4):e24354. 10.1002/jcla.24354
47
Lu Y-k Chen T-y Zhuang X-m Li J Jiang X-y Xu S-s et al Mechanism study of sijunzi decoction on rats with heart failure through regulating the imbalance of intestinal flora. (2021).
48
origin L . Romboutsia ilealis: how is this bacteria linked to healthy fats? (2022). Available online at: https://layerorigin.com/blogs/blog-layer-origin-nutrition/romboutsia-ilealis-how-is-the-bacteria-is-linked-to-healthy-fats. (Accessed April 9, 2024)
49
Gerritsen J Hornung B Renckens B van Hijum SA Martins dos Santos VAP Rijkers GT et al Genomic and functional analysis of Romboutsia ilealis CRIBT reveals adaptation to the small intestine. PeerJ (2017) 5:e3698. 10.7717/peerj.3698
50
Lanthier N Rodriguez J Nachit MY Hiel S Trefois P Neyrinck AM et al Lower abundance of Clostridium sensu stricto is associated with liver steatosis and fibrosis severity in a prospective cohort of obese patients with metabolic dysfunction-associated fatty liver disease. Acta Gastro-Enterologica Belgica (2021) 84:A20.
51
Kwan SY Jiao J Joon A Wei P Petty LE Below JE et al Gut microbiome features associated with liver fibrosis in Hispanics, a population at high risk for fatty liver disease. Hepatology (2022) 75(4):955ā67. 10.1002/hep.32197
52
Zhang B Wang X Xia R Li C . Gut microbiota in coronary artery disease: a friend or foe?Biosci Rep (2020) 40(5):BSR20200454. 10.1042/bsr20200454
53
Li J Lin S Vanhoutte PM Woo CW Xu A . Akkermansia muciniphila protects against atherosclerosis by preventing metabolic endotoxemia-induced inflammation in Apoeā/ā mice. Circulation (2016) 133(24):2434ā46. 10.1161/circulationaha.115.019645
54
Zagato E Pozzi C Bertocchi A Schioppa T Saccheri F Guglietta S et al Endogenous murine microbiota member Faecalibaculum rodentium and its human homologue protect from intestinal tumour growth. Nat Microbiol (2020) 5(3):511ā24. 10.1038/s41564-019-0649-5
55
Chai X Wen L Song Y He X Yue J Wu J et al DEHP exposure elevated cardiovascular risk in obese mice by disturbing the arachidonic acid metabolism of gut microbiota. Sci The Total Environ (2023) 875:162615. 10.1016/j.scitotenv.2023.162615
56
OƱate FP Chamignon C Burz SD Lapaque N Monnoye M Philippe C et al Adlercreutzia equolifaciens is an anti-inflammatory commensal bacterium with decreased abundance in gut microbiota of patients with metabolic liver disease. Int J Mol Sci (2023) 24(15):12232. 10.3390/ijms241512232
57
Wang M Zhang Y Miller D Rehman NO Cheng X Yeo J-Y et al Microbial reconstitution reverses early female puberty induced by maternal high-fat diet during lactation. Endocrinology (2020) 161(2):bqz041. 10.1210/endocr/bqz041
58
Bo T-B Zhang X-Y Kohl KD Wen J Tian S-J Wang D-H . Coprophagy prevention alters microbiome, metabolism, neurochemistry, and cognitive behavior in a small mammal. The ISME J (2020) 14(10):2625ā45. 10.1038/s41396-020-0711-6
Summary
Keywords
antiretroviral therapy (ART), ARV, abacavir, lamivudine, dolutegravir, microbiome, gut mucosa-associated microbes
Citation
Yanamadala Y, Muthumula CMR, Gokulan K, Karn K, Sutherland V, Cunny H, Santos JH and Khare S (2025) Gestational exposure to HIV drugs alters intestinal mucosa-associated microbial diversity in adult rat offspring. Exp. Biol. Med. 250:10564. doi: 10.3389/ebm.2025.10564
Received
03 March 2025
Accepted
07 July 2025
Published
13 August 2025
Volume
250 - 2025
Updates
Copyright
Ā© 2025 Yanamadala, Muthumula, Gokulan, Karn, Sutherland, Cunny, Santos and Khare.
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: Sangeeta Khare, sangeeta.khare@fda.hhs.gov
ā”These authors have contributed equally to this work
ā Present address: Vicki Sutherland, J&J Innovative Medicine, Spring House, PA, United States
Disclaimer
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.