Introduction

Approximately 40 million people around the world are living with HIV, with more than half being women of reproductive age3.

In 2018, a surveillance study in Botswana reported a higher incidence rate of neural tube defects (NTD) in infants born to mothers living with HIV taking DTG from conception, compared to other ART regimens, and women without HIV5,13.

Metabolic complications, including weight gain and hyperglycemia, have been reported with DTG use. Initiation of DTG-based treatment in both ART-experienced and treatment-naïve individuals has been linked to significant weight gain, with female sex, black race, increasing age, and concomitant use of tenofovir alafenamide, being risk factors17,28.

The mechanisms contributing to the reported DTG-associated metabolic alterations are not fully understood. Treatment with DTG has been shown to cause changes to adipose tissue composition, function, and signaling33 Whether the reported DTG-associated cell- and tissue-specific changes would have an impact on whole-body physiology leading to hyperglycemia, diabetes, or metabolic syndrome is not yet known.

With the goal of identifying potential mechanisms through which DTG could increase the risk for NTDs and to better understand the impact of DTG on female metabolic health, the aim of this study was to evaluate the effects of therapeutic and supratherapeutic doses of DTG-based ART on weight gain and glucose homeostasis over time in female mice.

Methods

Ethics statement and reporting

All animal work was approved by the University Health Network Animal Use Committee (#2575.29) and all methods were performed in accordance with Canadian Council on Animal Care guidelines. ARRIVE guidelines were followed in the reporting of this study.

Animals

Animal experiments were performed with in-house bred female C57BL/6J mice, original colony obtained from Jackson Laboratory RRID:IMSR JAX:000,664. Mice were individually housed starting at 6 weeks of age with ad libitum access to food (Teklad LM-485 Mouse/Rat Sterilizable Diet #7012) and water in a 12h light cycle-controlled facility. Mice were gavage trained with water 2 weeks prior to experiment initiation to minimize any potential stress associated with handling.

Treatment and sample size

Prescription DTG and tenofovir disoproxil fumarate (TDF)/ emtricitabine (FTC) pills were crushed and suspended in distilled water for daily oral gavage. Drug dosing was calculated to yield human plasma levels as previously described34. The therapeutic dose group (1xDTG) received 2.5 mg/kg DTG + 50/33.3 mg/kg TDF/FTC once daily (yielding DTG peak plasma concentrations of ~ 3,000 ng/ml and trough concentrations of ~ 150 ng/ml), and the supratherapeutic dose group (5xDTG) received 12.5 mg/kg + 50/33.3 mg/kg TDF/FTC once daily (yielding DTG peak plasma concentrations of ~ 12.000 ng/ml and trough concentrations of ~ 250 ng/ml). The control group received water. At 8 weeks of age, mice were randomly assigned to control (N = 15), 1xDTG (N = 13), or 5xDTG (N = 15) and were gavaged daily with 100µL of appropriate treatment for 9 weeks.

Assuming a mean of 1 and standard deviation of 0.2, we estimated that a sample size of 10 per group would give us an 80% power to detect a 25% change from the control group. A sample size of 15 per group would give us a power of 93% at an alpha of 0.05.

Metabolic assessment

Mice were assessed weekly for food intake, body weight, and fasted glucose. Food intake was assessed by weighing the food at the start and end of the week using a digital scale. Body weight was recorded using a digital scale. Blood glucose levels were assessed in the morning following a 4-6h fast using a glucometer (OneTouch UltraMini).

Biweekly, mice were overnight fasted for 13-15h and istered an oral glucose tolerance test (OGTT). Glucose was dosed at 0.5g/kg, which we determined to be well tolerated by lean mice during optimization experiments. Blood glucose levels were assessed using a glucometer at time point 0, 10, 20, 30, 60, 90, and 120 min post dosing. Area under the curve (AUC) above baseline glucose of the OGTT was calculated by first subtracting the zero time point glucose levels for each animal35. Hyperglycemia was defined as blood glucose measurements greater than mean + 2 standard deviations (SD) of glucose values of the control group over the course of the study.

Tissue collection

Mice were overnight fasted (13-15h) and euthanized at 9 weeks of treatment by CO2 inhalation followed by cervical dislocation. Liver, muscle, interscapular brown fat, gonadal white fat, and hypothalamus were collected and flash-frozen in liquid nitrogen. A portion of each lobe of the liver was frozen in OCT (Tissue-Tek 4583) for histological analysis. Mouse and organ weights were measured by digital scale.

Quantification of gene expression by qPCR

Collected tissues were ground using mortar and pestle on dry ice. Ground tissue was suspended in Trizol (Invitrogen Cat#15,596,018) and homogenized. Liver and muscle RNA was extracted by chloroform and purified by MinElute kit (Qiagen Cat#74,204). Tissues with high lipid content: hypothalamus, white adipose tissue, and brown adipose tissue were extracted using the Zymogen kit (Cat#R2050). RNA was quantified and quality checked by nanodrop, DNAse treated, and reverse transcribed into cDNA (BioRad Cat#1,725,035). Gene expression was quantified by RT-qPCR using SybrGreen mix (BioRad Cat#1,725,271) on Biorad CFX Opus 384. Full list of primers and sequences are described in Supplemental Table S1. Geometric mean of three housekeeping genes (Actb, Hprt, Ywhaz) was used to calculate relative expression using the ΔΔCt method36.

Quantification of protein expression by western blotting

Ground tissue was suspended in lysis buffer (80mM Tris HCl, 2% SDS, 10% glycerol) with protease inhibitors (Roche Cat#11,836,170,001) and homogenized. The samples were purified by several cycles of centrifugation followed by a freeze–thaw cycle. Protein concentration was quantified by Pierce BCA assay (Thermo Cat#23,227). 30µg of protein in 14µl of NuPage sample buffer (ThermoFisher Cat#NP0007), along with PageRuler (Thermo Cat#26,616) protein ladder, was loaded in acrylamide gels (BioRad Cat#4,561,096) and run at 100V. Proteins were blotted onto PVDF membranes (BioRad Cat#1,704,157) and stained with Ponceau (Sigma Cat#SLCL7892) for secondary loading control validation. The following antibodies were used: anti-FAS (Cell signaling Cat#3189S), anti-ACC1 (Cell signaling Cat#4190S), anti-DGAT1 (Proteintech Cat#11,561–1-AP). Blocking and secondary antibody incubation was done for 1 h at room temperature. Primary antibody incubation was done overnight at 4°C. Washes were done with TBST. The blots were developed with ECL (BioRad Cat#170–5060) and imaged with BioRad ChemiDoc. Each antibody was tested on a gradient of protein concentrations to that the blots were within the combined linear range. β-actin (anti-β-actin:HRP (Abcam Cat#20,272)) was used as the loading control.

ELISA

Leptin (R&D Cat#MOB00B), corticosterone (Enzo Cat#ADI-900–097), and insulin (Alpco Cat#80-INSMSU-E01) were quantified in mouse plasma by ELISA according to the manufacturer’s instructions.

Histological analysis

A section from each of the liver lobes was collected at sacrifice, embedded in OCT (Tissue-Tek 4583), and frozen. Livers were sectioned (8µm slices) using a Leica cryostat. The sections were air-dried for at least 30 min prior to Oil-Red-O (ORO, Sigma Cat# 01,391-250ml) staining according to published protocols37. Four to five lobes per mouse were included in the analysis. Each slide was digitally scanned and each lobe was annotated to remove edges, folds, or tears to ensure consistency of digital analysis. Number and size of each lipid droplet and ORO stain density was quantified for each lobe using HALO software (Indiga Labs) and an average across all lobes was reported.

Statistical analysis

Shapiro–Wilk test was used to assess the normality of data distribution. For normally distributed data, one-way ANOVA with a pairwise post-hoc with a Bonferroni correction was used to test for statistical significance between the three treatment arms. For non-parametric distributions, Kruskal–Wallis test with Dunn’s post-hoc with Bonferroni correction was used. Two-way ANOVA with Tukey’s post-test was used to compare between treatment arms over time. For correlation analyses, Pearson’s correlation was used. Generalized linear models were used to calculate difference in means with 95% confidence intervals compared to controls. For mRNA expression analyses, data were log-transformed prior to statistical analyses. Heat maps of the log-transformed means were generated. Comparison between groups was performed using multiple t-tests without assumption of consistent SD between groups. We corrected for multiple comparisons using the false discovery rate (FDR) approach using the method of Benjamini, Krieger, and Yekutieli, and with an FDR of 10%. Volcano plots, where the mean difference was plotted on the x-axis and the -log (q value) on the y-axis for gene expression, were employed to visualize genes that differed significantly between DTG-treated mice and controls. All statistical analyses were performed in R-4.0.5, Prism v8.2, and Stata v13.

Results

DTG-based ART in female mice is associated with a transient dysregulation of glycemic control, but not with excess weight gain or food intake

To assess metabolic alterations due to DTG-based ART, single-housed female C57Bl/6J mice were treated with either control, 1xDTG, or 5xDTG regimens once daily starting at 8 weeks of age. Their weekly food intake, body weight, and morning fasted blood glucose were recorded over 8 weeks. Female mice were selected for this study as female sex was one of the risk factors identified for DTG-related weight gain and we were interested in looking at glucose dysregulation as a risk factor for NTD.

Weekly food intake did not differ between treatments (Fig. 1C), although none reached significance.

Fig. 1
figure 1

DTG-based ART does not affect body weight, food intake, or 4-6 h fasted glucose, but causes a transient increase in overnight fasted glucose. Female C57BL/6J mice were treated with either control (grey circles), 1xDTG (red squares), or 5xDTG (blue triangles) for 9 weeks and their body weight (A), food intake (B), and 4-6 h fasted blood glucose (C) were measured weekly. Overnight (13-15 h) fasted glucose was measured every 2 weeks (E), and at time of sacrifice at week 9 (F). For (A-D) data are shown as means with SEM. Statistical analyses with 2-way ANOVA with Tukey’s post-test. For (E–F), data presented as dot plots with median shown. The dotted line represents the mean + 2SD of glucose for the control group (a value of 10.4 mmol/L). Mice with glucose over 10.4 mmol/L were deemed to be hyperglycemic. Statistical analyses by Kruskal–Wallis test with Dunn’s post-test. * p < 0.05. N = 15 for control, N = 13 for 1xDTG, and N = 15 for 5xDTG.

Full size image

Overnight (13-15h) fasted blood glucose (measured prior to istration of glucose in the OGTT) rose over the course of the study more prominently in the 1xDTG-treated group, but also mildly in the 5xDTG group (Fig. 1F).

Overnight fasted mice were also istered an oral glucose tolerance test (OGTT) at 2, 4, 6, and 8 weeks of treatment. There were no significant differences in the OGTT (Fig. 2B), between treatment arms at any time point in the study.

Fig. 2
figure 2

DTG treatment does not affect glucose tolerance. Female C57BL/6J mice treated with either control (grey circles), 1xDTG (red squares), or 5xDTG (blue triangles) were fasted overnight (13-15 h) and then istered an oral glucose tolerance test (OGTT; 0.5 g/kg) on weeks 2, 4, 6, and 8 following treatment initiation. Blood glucose response over 120 min is shown in (A). The baseline adjusted area under the curve (AUC) for each treatment arm over time is shown in (B). Data are shown as mean with SEM. N = 15 for control, N = 13 for 1xDTG, and N = 15 for 5xDTG. Statistical comparisons by 2-way ANOVA with Tukey’s post-test. * p < 0.05, ** p < 0.01.

Full size image

DTG treatment is associated with a progressive reduction in plasma leptin and an increase in corticosterone levels

As we observed a transient increase in overnight fasted glucose in the 1xDTG group and in some mice in the 5xDTG group, we next investigated changes to plasma hormones associated with regulation of glucose. We first measured insulin, the main glucoregulatory hormone, in overnight fasted plasma over 9 weeks of treatment. Insulin levels were similar between treatment groups, with a small non-significant increase seen at week 6 in the 5xDTG group and week 9 in the 1xDTG group (Sup. Fig. S1A, B).

Leptin is an adipocyte produced hormone with glucoregulatory function. Lower leptin secretion has been reported in DTG-treated adipocytes3D).

Fig. 3
figure 3

Lower leptin and higher corticosterone levels with DTG treatment. Female C57BL/6J mice treated with either control (grey circles), 1xDTG (red squares), or 5xDTG (blue triangles) for a period of 9 weeks and leptin and corticosterone were assessed in overnight fasted plasma. (A) Difference in mean in fasted plasma leptin levels from control (dotted line) with 95% confidence interval (CI) for 1xDTG (red) and 5xDTG (blue) treated mice at week 2, 4, 6, and 8 of treatment. (B) Fasted plasma leptin levels at time of sacrifice at week 9 of treatment. (C) Leptin mRNA levels (log transformed) measured in white adipose tissue collected at time of sacrifice at week 9 of treatment. (D) Correlation between plasma leptin collected at sacrifice and white adipose tissue leptin mRNA. (E) Difference in mean in fasted plasma corticosterone levels from control (dotted line) with 95% CI for 1xDTG (red) and 5xDTG (blue) treated mice at week 2, 4, 6, and 8 of treatment. (F) Fasted plasma corticosterone levels at time of sacrifice at week 9 of treatment. (G) Correlation between white adipose tissue leptin mRNA and week 9 plasma corticosterone levels. For (A) and (E) mean difference with 95%CI calculated using generalized linear models. Not crossing the dotted line indicates significant difference from control. For (B), (C), and (F) box plots show median (line), interquartile range (box), and range (whiskers), with individual data points shown. Statistical comparisons for (B), (C), and (F) by ANOVA with Tukey’s post-test. Pearson r with p-value shown for (D) and (G). N = 14 for control, N = 12 for 1xDTG, N = 14 for 5xDTG. * p < 0.05, ** p < 0.01.

Full size image

Due to the glucoregulatory role of leptin through cortisol (corticosterone in mice),3G), as well as plasma leptin levels (r = -0.45, p = 0.0035).

DTG treatment increases leptin receptor expression in muscle and liver.

Since we observed persistently low leptin levels with DTG treatment, even-though euglycemia was restored by week 8 and 9, we next examined if DTG-treated mice could have adapted to the lower leptin levels by upregulating leptin receptor expression. We observed significantly higher leptin receptor expression in the muscle for both DTG groups compared to controls (Fig. 40.

Fig. 4
figure 4

DTG treatment is associated with increase in leptin receptor mRNA levels in muscle and liver. Leptin receptor mRNA levels (log-transformed) in muscle (A) and liver (C) of mice treated with either control (grey), 1xDTG (red), or 5xDTG (blue) following 9 weeks of treatment. Box plots show median (line), interquartile range (box), and range (whiskers), with individual data points shown. Statistical comparison by ANOVA with Tukey’s post-test. Correlations between muscle (B) and liver (D) leptin receptor mRNA levels and plasma leptin at week 9, assessed using Pearson r. N = 14 for control, N = 12 for 1xDTG, N = 15 for 5xDTG. ** p < 0.01, *** p < 0.001.

Full size image

DTG-based ART exposure has a mild effect on hepatic and muscle gene expression.

Since we observed a transient hyperglycemia in the DTG-treated mice, we speculated that compensatory mechanisms may have contributed to restoring or maintaining euglycemia. The liver is the key regulator of endogenous glucose production in a fasted state and integrates peripheral and central endocrine and neural signals for maintenance of glucose homeostasis. We measured gene expression of proteins involved in gluconeogenesis, fatty acid catabolism, and lipogenesis in livers collected at sacrifice following 9 weeks of treatment. We hypothesized that a reduction in endogenous glucose production and/or fatty acid catabolism could potentially offset the mechanisms contributing to hyperglycemia and either restore euglycemia or protect mice from developing hyperglycemia. Means with 95% CI of mRNA levels for all genes assessed are shown in Table 2 for muscle, and Supplemental Table S2 for brown adipose tissue, Table S3 for hypothalamus, and Table S4 for white adipose tissue.

Table 1 Hepatic gene expression after 9 weeks of treatment.
Full size table
Table 2 Muscle gene expression after 9 weeks of treatment.
Full size table

In the liver, we did not observe changes to the gluconeogenic phosphoenolpyruvate carboxykinase (Pepck) or glucose-6-phosphatase (G6p) between treatment arms, however both G6p and Pepck showed greater spread in the data in the DTG groups (Table 1).

Fig. 5
figure 5

Hepatic and muscle gene expression in control, 1xDTG, and 5xDTG mice. Quantification of mRNA expression levels in liver (A–D) and muscle (E–G) tissue collected at time of sacrifice at 9 weeks of treatment in overnight fasted animals. (A) Fatty acid synthase (Fas), (B) diacylglycerol O-acyltransferase 1 (Dgat1), (C) acetyl-CoA carboxylase (Acc1), (D) insulin receptor substrate 2 (Irs2), (E) Cd36, (F) Irs1, (G) peroxisome proliferator-activated receptor gamma (Ppparγ). Data are log-transformed arbitrary units (AU) and are shown as box plots with the line indicating median, the box interquartile range, and the whiskers range, with individual data points shown. Statistical analyses using ANOVA with Tukey’s post-test. N = 14–15 for control, N = 12–13 for 1xDTG, N = 15 for 5xDTG. * p < 0.05, ** p < 0.01.

Full size image

To examine if changes in liver mRNA levels translated to changes in protein levels, we quantified FASN, ACC1, and DGAT1 protein levels in liver lysates using Western blot. ACC1 and FASN were significantly lower in the 5xDTG group only (Fig. S2A and B). DGAT1 levels did not differ between treatment groups (Fig. S2C).

In muscle, we observed higher levels of Cd36 that reached significance in the 1xDTG group, higher levels of Irs1 in both DTG groups, and significantly lower levels of Pparγ in the 1xDTG group only (Fig. 2).

Other potential pathways that could stabilize glycemia include increased brown adipose lipid uptake and thermogenesis, and changes to hypothalamic energy balance. We did not observe changes to transcripts of uncoupling protein-1 (U1), responsible for thermogenesis, nor of brown adipose tissue markers cell death-inducing DFFA-like effector A (Cidea) and type II iodothyronine deiodinase (Dio2) in brown adipose tissue (Table S2). In hypothalamic tissue, we observed a non-significant upregulation of orexigenic agouti-related peptide (Agrp) mRNA expression in both DTG-treated groups compared to control, suggesting a negative energy balance in some animals. No differences were observed in the levels of Pomc, Bdnf, or Trh (Table S3). As noted above, leptin (Lep) transcripts were downregulated in gonadal white adipose tissue in DTG-treated mice, although this did not quite reach significance (Fig. 3C). No differences were observed between groups in adiponectin, Pparγ, Fas, or Srebp1c levels in white adipose tissue (Table S4).

Differential regulation of metabolic genes in the liver of 5xDTG mice that maintained euglycemia.

The greater spread in the expression data observed in the DTG groups, prompted us to examine if gene expression patterns may differ by glycemia status (euglycemic vs. hyperglycemic). Expression data for all the genes analyzed in the liver were log transformed and the means displayed in a heat map stratified by treatment arm and glycemia state (Fig. 6B top left ). Tdo was the only gene significantly upregulated, with the upregulation seen in the hyperglycemic–5xDTG.

Fig. 6
figure 6

Downregulation of gluconeogenic and lipogenic genes in the liver of DTG-treated mice that maintained euglycemia. (A) Heat map showing means of log-transformed mRNA expression levels in the liver stratified by treatment arm (control, 1xDTG, 5xDTG) and glycemia state (euglycemic (e) or hyperglycemic (h)). All control mice remained euglycemic. (B) Volcano plots showing mean difference for each group compared to control on the x-axis, and -log(q-value) on the y-axis for all genes. A positive mean difference indicates upregulation compared to control, a negative mean difference indicates downregulation compared to control. Anything above the horizontal dotted lines indicates significant difference. Statistical comparisons using false discovery rate of 10%. N = 6 for 1xDTG-euglycemic, N = 10 for 5xDTG-euglycemic, N = 6 for 1xDTG-hyperglycemic, N = 5 for 5xDTG-hyperglycemic.

Full size image

A similar analysis in muscle (Fig. 7) showed an upregulation of Cd36 and Lepr in both the euglycemic–5xDTG and the hyperglycemic–1xDTG groups. However, the hyperglycemic–1xDTG group also showed a downregulation of Pparγ and Glut4, suggesting perhaps that Cd36 and LepR upregulation may be protective, but only in the absence of Pparγ and/or Glut4 downregulation.

Fig. 7
figure 7

Gene expression profiles in muscle stratified by treatment and glycemia status. (A) Heat map showing means of log-transformed mRNA expression levels in the muscle stratified by treatment arm (control, 1xDTG, 5xDTG) and glycemia state (euglycemic (e) or hyperglycemic (h)). All control mice remained euglycemic. (B) Volcano plots showing mean difference for each group compared to control on the x-axis, and -log(q-value) on the y-axis for all genes. A positive mean difference indicates upregulation compared to control, a negative mean difference indicates downregulation compared to control. Anything above the horizontal dotted lines indicates significant difference. Statistical comparisons using false discovery rate of 10%. N = 6 for 1xDTG-euglycemic, N = 10 for 5xDTG-euglycemic, N = 6 for 1xDTG-hyperglycemic, N = 5 for 5xDTG-hyperglycemic.

Full size image

No hepatic lipid accumulation in DTG-treated mice

As we observed increased corticosterone and lower leptin combined with mild downregulation in gene expression of factors associated with lipogenesis (Fas, Acc1) with DTG treatment, we stained liver sections with oil-red-O (ORO) to assess triacylglycerol accumulation. We did not observe any differences in average lipid droplet size or ORO stain density (Fig. S3).

Discussion

In this study we explored the impact of DTG-based ART on weight gain and glucose homeostasis, in female mice over a 9-week treatment window. We istered DTG at a therapeutic (1xDTG) or a supratherapeutic (5xDTG) dose in combination with TDF/FTC. DTG treatment, at either dose, was not associated with excess weight gain, but was associated with a transient fasted hyperglycemia that was more prominent in the 1xDTG group. DTG treatment was also associated with a gradual decline in plasma leptin in both DTG groups that correlated with lower leptin mRNA expression in white adipose tissue, and a gradual increase in plasma corticosterone in 1xDTG group that remained significantly elevated at week 9 of treatment. We observed evidence of induction of compensatory mechanisms that were more prominent in the mice that maintained euglycemia in the 5xDTG group, including downregulation of genes involved in gluconeogenesis and lipogenesis in the liver. We also observed increased expression of leptin receptor in muscle and liver with DTG treatment, that may have served to improve leptin sensitivity in the context of low leptin levels.

Several clinical studies have reported greater weight gain with initiation of DTG-based regimens in both ART-experienced and naive patients, especially in women and older individuals45. It has been proposed that weight gain in the context of DTG use may be related to the body’s compensatory response to a period of weight suppression. The absence of HIV infection in our model prevents us from investigating this possibility.

In clinical studies, hyperglycemia, metabolic syndrome, and diabetes following initiation of DTG have been reported and occurred in patients with normal BMI, experiencing weight gain or loss, and without prior history of insulin resistance; it is therefore unclear whether DTG-associated weight gain increases the risk for diabetes or if DTG is interfering with glucose uptake and secretion through another mechanism17,47. Our data suggest that DTG impacts glucose homeostasis independent of weight gain, and that lower leptin and higher corticosterone may contribute this dysregulation. This is consistent with some case reports of patients experiencing hyperglycemia while having normal insulin and having had lost weight since DTG-treatment initiation49. Our findings longer term follow-up of glycemic status in people initiating DTG-based ART, as incident hyperglycemia may be transient.

We observed a gradual decrease in plasma leptin with DTG treatment that persisted throughout the experiment. Declines in leptin levels have been observed in patients switching from protease inhibitor-based regimens to DTG or raltegravir-based regimens51. In our study DTG-induced declines in leptin plasma levels coincided with a rise in glucose levels, suggesting that the lower leptin levels may have affected glucose clearance leading to hyperglycemia. However, while leptin levels remained low in both DTG groups, all mice returned to euglycemia by week 9. Increased leptin receptor expression in muscle and liver may have contributed to this, causing a compensatory increase in leptin sensitivity.

DTG was also associated with gradual increases in corticosterone levels that were more prominent in the 1xDTG group. Corticosterone/cortisol also plays a role in glucose homeostasis. Cortisol acts on the liver to increase gluconeogenesis and decrease glycogen synthesis, in muscle to decrease glucose uptake and consumption, and in adipose tissue to increase lipolysis52. The higher levels of corticosterone observed in the 1xDTG group may have contributed to the slightly elevated glucose levels observed in the 1xDTG group at week 9. However, corticosterone levels did not correlate with glucose levels at any other point during treatment.

A portion of our DTG-treated mice, particularly those in the 5xDTG group, did not develop hyperglycemia. We observed downregulation of key genes and proteins involved in de novo lipogenesis including ACC1 and FASN in both DTG groups but more prominently in the 5xDTG group. These changes suggest a possible downregulation of de novo lipogenesis and increase in fatty acid oxidation, that may have protected DTG-treated mice from becoming hyperglycemic, and may have contributed to the transiency of the glucose dysregulation55. The downregulation of Fas we observed here suggests that compensatory mechanisms are engaged following DTG treatment that ameliorate at least some of the drug-associated metabolic alterations.

One of our aims for this study was to identify potential mechanisms through which DTG could have caused the initial NTD signal5,7. This is ed by our findings here of a stronger and more varied downregulation of lipogenic and gluconeogenic genes in the 5xDTG group that maintained euglycemia. It is possible that a transient hyperglycemia, and the induction of a compensatory response with longer term DTG use, could explain in part the NTD signal that was observed at the initial roll out of DTG but waned over time.

Our study has some limitations. This study was carried out in healthy mice and thus does not capture the effects of HIV or the interaction of HIV and antiretroviral therapy. Although this may be a limitation when generalizing findings to people living with HIV, experiments in healthy animals provide better characterization of the direct effects of DTG-based ART. A limitation of this study is the absence of body composition and indirect calorimetry data, which would provide greater insight in interpreting our results. In addition, while we looked at the expression profile of a large number of metabolic related genes in a variety of relevant tissue, the activity of these enzymes is transcriptionally and post-transcriptionally regulated. Further studies focusing on DTG effects on lipogenesis and fatty acid oxidation are merited.

In summary, DTG treatment of female mice was associated with transient hyperglycemia, and persistent decreases in leptin and elevations in corticosterone. Prospective monitoring of glucose, as well as leptin and cortisol, in patients initiating DTG-based treatment would be of interest. Further investigations of the effects on DTG on the hypothalamic–pituitary–adrenal axis should also be considered, given our observed effects on corticosterone/cortisol. Taken together, DTG effects on metabolic health appear mild and affecting only a portion of the treated mice, with evidence of transient effects and induction of compensatory mechanisms.