Doubling dolutegravir dosage reduces the viral reservoir in ART-treated people with HIV

  1. Céline Fombellida-Lopez
  2. Aurelija Valaitienė
  3. Lee Winchester
  4. Nathalie Maes
  5. Patricia Dellot
  6. Céline Vanwinge
  7. Aurélie Ladang
  8. Etienne Cavalier
  9. Fabrice Susin
  10. Dolores Vaira
  11. Marie-Pierre Hayette
  12. Catherine Reenaers
  13. Michel Moutschen
  14. Courtney V Fletcher
  15. Alexander O Pasternak  Is a corresponding author
  16. Gilles Darcis  Is a corresponding author
  1. Laboratory of Immunology and Infectious Diseases, GIGA-Institute, University of Liège, Belgium
  2. Laboratory of Experimental Virology, Department of Medical Microbiology, Amsterdam UMC, University of Amsterdam, Netherlands
  3. Antiviral Pharmacology Laboratory, University of Nebraska Medical Center, United States
  4. Department of Biostatistics and Medico-Economic Information, University Hospital of Liège, Belgium
  5. Department of General Internal Medicine and Infectious Diseases, University Hospital of Liège, Belgium
  6. GIGA Flow Cytometry Platform, University of Liège, Belgium
  7. Department of Clinical Chemistry, University Hospital of Liège, Belgium
  8. Laboratory of Clinical Microbiology, University Hospital of Liège, Belgium
  9. Department of Gastroenterology, University Hospital of Liège, Belgium
8 figures, 4 tables and 1 additional file

Figures

Study design.

(A) Inclusion criteria, study groups, and sampling time points. (B) Time points where markers were measured.

Figure 2 with 1 supplement
Antiretroviral drug concentrations in plasma and rectal tissue.

Plasma (A) and tissue (B) concentrations of DTG, as well as plasma (C) and tissue (D) concentrations of 3TC, were compared between days 0 (D0) and 84 (D84) in intensified and control groups. Wilcoxon tests were used to calculate statistical significance.

Figure 2—figure supplement 1
Plasma concentrations of DTG and 3TC in intensified (blue) and control (red) groups.

Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10.

Figure 3 with 3 supplements
Longitudinal dynamics of HIV reservoir markers.

Fold change of total HIV DNA in PBMCs (A), US HIV RNA in PBMCs (B), and US RNA/total DNA ratios (C) from baseline on days 1, 28, 56, and 84 of the study in the intensified (blue) and control (red) groups. Median values and interquartile ranges (IQRs) are shown. Linear mixed-effects modelling on log10-transformed values was used to calculate statistical significance. p values at the bottom of the graphs show the significance of between-group comparisons and those on top of the graphs show the significance of comparisons of the changes from baseline with zero in each group separately (intercept-only analysis). An upward or downward facing arrow next to a p value indicates a statistically significant increase or decrease from baseline, respectively. Participant numbers in both groups per time point are indicated below the graphs. Comparisons of intact HIV DNA in PBMCs (D) and of total HIV DNA in rectal tissue (E) between days 0 (D0) and 84 (D84) in the intensified and control groups. Open circles depict undetectable values, assigned the values corresponding to 50% of the assay detection limits. Wilcoxon tests were used to calculate statistical significance. All p values are marked red if significant (<0.05).

Figure 3—figure supplement 1
Longitudinal dynamics of total HIV DNA in PBMCs.

(A) Levels of total HIV DNA in intensified (blue) and control (red) groups. Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10. (B) Total HIV DNA was compared between days 0 (D0) and 84 (D84) in intensified and control groups. Open circles depict undetectable values, assigned the values corresponding to 50% of the assay detection limits. Wilcoxon tests were used to calculate statistical significance. (C) Time-weighted changes of total HIV DNA from baseline were compared between intensified and control groups and compared with 0 in each group separately. (D) Day 1/day 0, day 28/day 0, day 56/day 0, and day 84/day 0 ratios of total HIV DNA were compared between intensified and control groups and compared with 1 in each group separately. Panels C and D: Mann–Whitney tests and one-sample Wilcoxon tests were used to calculate between- and within-group statistical significance, respectively. If a within-group p value is <0.2, an upward or downward facing arrow is depicted next to this p value, indicating a trend to increase or decrease from baseline, respectively. p values are marked red if significant (<0.05). All statistical analysis was performed on log10-transformed values of total HIV DNA. (E) Individual trajectories of total HIV DNA (fold change from baseline) in intensified (blue) and control (red) groups.

Figure 3—figure supplement 2
Longitudinal dynamics of US HIV RNA in PBMCs.

(A) Levels of US HIV RNA in intensified (blue) and control (red) groups. Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10. (B) US HIV RNA was compared between days 0 (D0) and 84 (D84) in intensified and control groups. Open circles depict undetectable values, assigned the values corresponding to 50% of the assay detection limits. Wilcoxon tests were used to calculate statistical significance. (C) Time-weighted changes of US HIV RNA from baseline were compared between intensified and control groups and compared with 0 in each group separately. (D) Day 1/day 0, day 28/day 0, day 56/day 0, and day 84/day 0 ratios of US HIV RNA were compared between intensified and control groups and compared with 1 in each group separately. Panels C and D: Mann–Whitney tests and one-sample Wilcoxon tests were used to calculate between- and within-group statistical significance, respectively. If a within-group p value is <0.2, an upward or downward facing arrow is depicted next to this p value, indicating a trend to increase or decrease from baseline, respectively. p values are marked red if significant (<0.05). All statistical analysis was performed on log10-transformed values of US HIV RNA. (E) Individual trajectories of US HIV RNA (fold change from baseline) in intensified (blue) and control (red) groups.

Figure 3—figure supplement 3
Longitudinal dynamics of US RNA/total DNA ratio in PBMCs.

(A) Levels of US RNA/total DNA ratio in intensified (blue) and control (red) groups. Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10. (B) US RNA/total DNA ratios were compared between days 0 (D0) and 84 (D84) in intensified and control groups. Open circles depict undetectable values. Wilcoxon tests were used to calculate statistical significance. (C) Time-weighted changes of US RNA/total DNA ratio from baseline were compared between intensified and control groups and compared with 0 in each group separately. (D) Day 1/day 0, day 28/day 0, day 56/day 0, and day 84/day 0 ratios of US RNA/total DNA ratio were compared between intensified and control groups and compared with 1 in each group separately. Panels C and D: Mann–Whitney tests and one-sample Wilcoxon tests were used to calculate between- and within-group statistical significance, respectively. If a within-group p value is <0.2, an upward or downward facing arrow is depicted next to this p value, indicating a trend to increase or decrease from baseline, respectively. p values are marked red if significant (<0.05). All statistical analysis was performed on log10-transformed values of US RNA/total DNA ratios. (E) Individual trajectories of US RNA/total DNA ratio (fold change from baseline) in intensified (blue) and control (red) groups.

Figure 4 with 1 supplement
Longitudinal dynamics of cellular markers of immune activation and exhaustion.

Graphs show changes of (A) CD4+ cell markers and (B) CD8+ cell markers from baseline on days 1, 28, 56, and 84 of the study in the intensified (blue) and control (red) groups. Median values and interquartile ranges (IQRs) are shown. Participant numbers: intensified group, n=10; control group, n=10. Linear mixed-effects modelling was used to calculate statistical significance. p values at the bottom of the graphs show the significance of between-group comparisons and those on top of the graphs show the significance of comparisons of the changes from baseline with zero in each group separately (intercept-only analysis). An upward or downward facing arrow next to a p value indicates a statistically significant increase or decrease from baseline, respectively. p values are marked red if significant (<0.05).

Figure 4—figure supplement 1
Levels of cellular markers of immune activation and exhaustion.

Levels of cellular markers of immune activation and exhaustion on CD4+ (A) and CD8+ (B) T cells in intensified (blue) and control (red) groups. Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10.

Figure 5 with 1 supplement
Longitudinal dynamics of inflammation markers in plasma and tissue.

(A) Plasma inflammation markers. Graphs show changes of markers from baseline on days 1, 28, 56, and 84 of the study in the intensified (blue) and control (red) groups. Median values and interquartile ranges (IQRs) are shown. Participant numbers: intensified group, n=10; control group, n=10. Linear mixed-effects modelling was used to calculate statistical significance. p values at the bottom of the graphs show the significance of between-group comparisons and those on top of the graphs show the significance of comparisons of the changes from baseline with zero in each group separately (intercept-only analysis). An upward or downward facing arrow next to a p value indicates a statistically significant increase or decrease from baseline, respectively. p values are marked red if significant (<0.05). (B) Tissue inflammation markers. Markers were compared between days 0 (D0) and 84 (D84) in the intensified and control groups. Open circles depict undetectable values, assigned the values corresponding to 50% of the assay detection limits. Wilcoxon tests were used to calculate statistical significance.

Figure 5—figure supplement 1
Levels of inflammation markers in plasma in intensified (blue) and control (red) groups.

Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10.

Figure 6 with 1 supplement
Longitudinal dynamics of clinical markers.

Graphs show changes of markers from baseline on days 1, 28, 56, and 84 of the study in the intensified (blue) and control (red) groups. Median values and interquartile ranges (IQRs) are shown. Participant numbers: intensified group, n=10; control group, n=10. Linear mixed-effects modelling was used to calculate statistical significance. p values at the bottom of the graphs show the significance of between-group comparisons and those on top of the graphs show the significance of comparisons of the changes from baseline with zero in each group separately (intercept-only analysis). An upward or downward facing arrow next to a p value indicates a statistically significant increase or decrease from baseline, respectively. p values are marked red if significant (<0.05).

Figure 6—figure supplement 1
Levels of clinical markers in intensified (blue) and control (red) groups.

Median values and interquartile ranges are shown. Participant numbers: intensified group, n=10; control group, n=10.

Figure 7 with 2 supplements
Spearman correlogram of baseline parameters.

A heat map is used to indicate the strengths of associations between parameters. Red indicates a negative correlation, and blue indicates a positive correlation.

Figure 7—figure supplement 1
Spearman correlogram of time-weighted changes from baseline.

A heat map is used to indicate the strengths of associations between parameters. Red indicates a negative correlation, and blue indicates a positive correlation.

Figure 7—figure supplement 2
Spearman correlogram of changes between days 0 and 84.

A heat map is used to indicate the strengths of associations between parameters. Red indicates a negative correlation, and blue indicates a positive correlation.

Appendix 3—figure 1
Gating strategy of the flow cytometry measurements.

Tables

Table 1
Clinical characteristics of participants at baseline.
All(n = 20)Control group (n = 10)DTG group(n = 10)p*
Sex, male19 (95.0)9 (90.0)10 (100.0)1.00
Age, years52 (43–60; 25–73)52 (40–62; 25–73)51 (46–59; 36–66)0.94
BMI, kg/m224 (22–26; 18–36)25 (22–27; 19–30)23 (22–26; 18–36)0.62
Ethnicity
 Caucasian17 (85.0)9 (90.0)8 (80.0)1.00
 African2 (10.0)1 (10.0)1 (10.0)
 Maghrebi1 (5.0)0 (0.0)1 (10.0)
HLA typing B57/01, negative (n = 19)19 (100.0)9 (100.0)10 (100.0)-
Smoking, smoker, or ex-smoker10 (50.0)4 (40.0)6 (60.0)0.66
Time since first positive HIV serology, years12.6 (6.8–18.2; 3.5–25.4)9.8 (7.3–13.3;
3.5–25.4)
15.6 (7.6–18.5; 3.5–21.5)0.58
Cumulative time of untreated
HIV infection, years
0.9 (0.1–3.5;
0.0–17.2)
1.8 (0.1–5.0;
0.0–17.2)
0.6 (0.1–2.1;
0.0–6.1)
0.47
Cumulative time of viral suppression, years7.4 (4.8–9.2;
2.9–16.5)
5.6 (4.7–7.7;
3.3–16.2)
7.9 (5.8–12.4;
2.9–16.5)
0.32
Continuous time of viral suppression before study, years3.9 (3.3–5.8;
0.4–10.2)
4.4 (3.7–6.9;
3.3–10.2)
3.4 (3.2–4.9;
0.4–7.9)
0.12
Time on DTG-containing ART regimen before study, years3.5 (3.4–3.7;
2.2–3.9)
3.5 (3.4–3.8;
2.2–3.9)
3.5 (3.3–3.5;
2.7–3.8)
0.41
Nadir CD4+ count, cells/mm3277 (147–410;
24–747)
277 (205–340;
24–521)
292 (136–508;
30–747)
0.65
First measured plasma viral load, log10 HIV RNA copies/ml4.52 (3.91–5.35; 2.48–6.60)4.52 (4.36–4.96; 3.00–6.60)4.55 (3.69–5.38; 2.48–5.70)0.85
CD4+ count, cells/mm3775 (659–1144; 225–1383)732 (672–970; 225–1267)886 (668–1217; 492–1383)0.39
CD8+ count, cells/mm3745 (637–1288; 377–2017)723 (535–931; 377–2017)939 (690–1310; 501–1853)0.39
CD4/CD8 ratio0.94 (0.66–1.47; 0.32–2.24)1.00 (0.64–1.36; 0.32–1.82)0.93 (0.69–1.56; 0.35–2.24)0.82
Total HIV DNA in PBMCs, copies/106 cells175.5 (70.0–467.3; 15.1–1090)99.4 (37.4–231.5; 15.1–333.0)421.0 (139.5–764.3; 61.8–1090)0.015
Intact HIV DNA in PBMCs, copies/106 cells55.2 (19.2–87.0; 5.17–624.2)29.5 (16.2–97.1; 10.3–624.2)58.6 (44.7–144.7; 5.17–292.6)0.36
US HIV RNA in PBMCs, copies/μg total RNA92.1 (13.3–290.3; 2.53–1580)72.2 (12.2–253.8; 6.18–486)124.6 (17.0–591.0; 2.53–1580)0.58
US RNA/total DNA ratio in PBMCs0.40 (0.07–1.47; 0.03–3.96)1.21 (0.11–2.81; 0.04–3.96)0.33 (0.06–0.70; 0.03–2.73)0.21
Total HIV DNA in rectal tissue, copies/106 cells477.0 (271.3–971.5; 7.27–1720)547.0 (291.0–1215; 53.4–1720)403.0 (194.6–952.5; 7.27–1510)0.44
Plasma DTG concentration, ng/mL3287 (2602–5087; 237–6593)3560 (2953–4455; 1969–5536)3266 (2245–5983; 237–6593)0.82
Tissue DTG concentration, ng/g634 (533–830; 303–1810)737 (534–852; 482–1810)581 (535–714; 303–1037)0.44
Plasma 3TC concentration, ng/ml316 (148–731;
50–1616)
431 (214–941; 102–1616)246 (117–439;
50–1156)
0.26
Tissue 3TC concentration, ng/g2114 (1417–2345; 90–4495)2193 (1590–2784; 1215–4495)2114 (1233–2317; 900–3187)0.50
CDC classification-
 A14 (20.0)1 (10.0)3 (30.0)
 A29 (45.0)6 (60.0)3 (30.0)
 A33 (15.0)1 (10.0)2 (30.0)
 B31 (5.0)0 (0.0)1 (10.0)
 C33 (15.0)2 (20.0)1 (10.0)
HBV status-
 Immune12 (60.0)8 (80.0)4 (40.0)
 Non-immune, not infected5 (25.0)2 (20.0)3 (30.0)
 Isolated HBc Ab2 (10.0)0 (0.0)2 (20.0)
 Cured hepatitis B1 (5.0)0 (0.0)1 (10.0)
HCV status1.00
 Not infected19 (95.0)10 (10.0)9 (90.0)
 Recovered1 (5.0)0 (0.0)1 (10.0)
  1. *

    Mann–Whitney tests were used for continuous variables and Fisher’s exact tests were used for categorical variables.

  2. Data are medians (interquartile ranges, followed by ranges) for continuous variables and numbers (percentages) for discrete variables.

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
AntibodyAPC/Fire 810 anti-human CD3 (Mouse monoclonal)Sony BiotechnologyCat# 23242901:25
AntibodyPE/Fire 700 anti-human CD4 (Mouse monoclonal)Sony BiotechnologyCat# 23233301:25
AntibodyPerCP anti-human CD8a (Mouse monoclonal)Sony BiotechnologyCat# 21046101:25
AntibodyPE Mouse anti-human CD279 (PD-1) (Mouse monoclonal)BD BiosciencesCat# 560795; RRID:AB_20339891:25
AntibodyBV421 Mouse anti-human TIGIT (Mouse monoclonal)BD BiosciencesCat# 747844; RRID:AB_28723071:25
AntibodyBB515 Mouse anti-human HLA-DR (Mouse monoclonal)BD BiosciencesCat# 564516; RRID:AB_27328461:25
AntibodyBV711 Mouse anti-human CD38 (Mouse monoclonal)BD BiosciencesCat# 563965; RRID:AB_27385161:25
Commercial assay or kitDNA-free DNA Removal KitThermo Fisher ScientificCat# AM1906
Commercial assay or kitTaqMan β-Actin Detection ReagentsThermo Fisher ScientificCat# 401846
Commercial assay or kitTaqMan Ribosomal RNA Control ReagentsThermo Fisher ScientificCat# 4308329
Commercial assay or kitPuregene Cell KitQIAGENCat# 158043
Commercial assay or kitReliaPrep gDNA Tissue Miniprep SystemPromegaCat# A2051
Commercial assay or kitProcartaPlex Human Inflammation Panel, 20plexThermo Fisher ScientificCat# EPX200-12185-901
Commercial assay or kitHuman CD14 ELISA Kit – QuantikineR&D SystemsCat# DC140
Chemical compound, drugPlatinum Quantitative PCR SuperMix-UDGThermo Fisher ScientificCat# 11730-025
Chemical compound, drugSuperScript III reverse transcriptaseThermo Fisher ScientificCat# 18080-085
Chemical compound, drugRandom primersThermo Fisher ScientificCat# 48190-011
Chemical compound, drugRNaseOUT Recombinant Ribonuclease InhibitorThermo Fisher ScientificCat# 10777-019
Chemical compound, drugddPCR Supermix for Probes (No dUTP)Bio-RadCat# 1863024
Chemical compound, drugBglI restriction enzymeThermo Fisher ScientificCat# ER0071
Chemical compound, drugTaqMan Gene Expression Assay, IL-1βThermo Fisher ScientificCat# 4331182Hs01555410_m1
Chemical compound, drugTaqMan Gene Expression Assay, IL-6Thermo Fisher ScientificCat# 4331182Hs00174131_m1
Chemical compound, drugTaqMan Gene Expression Assay, IFN-γThermo Fisher ScientificCat# 4331182Hs00989291_m1
Chemical compound, drugTaqMan Gene Expression Assay, IL-17αThermo Fisher ScientificCat# 4331182Hs00174383_m1
Chemical compound, drugTaqMan Gene Expression Assay, TNF-αThermo Fisher ScientificCat# 4331182Hs00174128_m1
Chemical compound, drugTaqMan Gene Expression Assay, GAPDHThermo Fisher ScientificCat# 4331182Hs02758991_g1
Software, algorithmPrism 10.2.0GraphPad Softwarehttps://www.graphpad.com/
RRID:SCR_002798
Statistics
Software, algorithmIBM SPSS Statistics 28.0.1.0IBMhttps://www.ibm.com/products/spss-statistics
RRID:SCR_016479
Statistics
Software, algorithmQuantaSoft 1.7.4Bio-RadRRID:SCR_025696ddPCR data analysis
Software, algorithmRotor-Gene 2.3.5QIAGENRRID:SCR_015740qPCR data analysis
Software, algorithmFlowJo 10.8.1Becton Dickinsonhttps://www.flowjo.com/
RRID:SCR_008520
Flow cytometry data analysis
Appendix 2—table 1
Comparison of models for the virological markers (between-group analysis).
Dependent variableExplanatory variablesTreatment only(repeated-measures model)Treatment + time point(as discrete variable)Treatment + time(as continuous variable)
Estimate (95% CI)p*Estimate (95% CI)pEstimate (95% CI)p
Total DNA
(log10 fold change from baseline)
Treatment, intensified vs. control group−0.38 (−0.55 to −0.21)2.3 × 10–5−0.38 (−0.55 to −0.20)5.8 × 10–5−0.38 (−0.56 to −0.20)6.5 × 10–5
Time, per day---0.150.001 (−0.002 to 0.004)0.41
US RNA
(log10 fold change from baseline)
Treatment, intensified vs. control group−0.56 (−0.98 to −0.14)0.010−0.49 (−0.94 to −0.04)0.035−0.49 (−0.93 to −0.05)0.028
Time, per day---0.99−0.001 (−0.008 to 0.006)0.82
US RNA/total DNA ratio
(log10 fold change from baseline)
Treatment, intensified vs. control group−0.36 (−0.77 to 0.06)0.090−0.29 (−0.73 to 0.16)0.20−0.28 (−0.71 to 0.16)0.21
Time, per day---0.91−0.002 (−0.009 to 0.005)0.62
  1. *

    p values were calculated by type III tests of fixed effects.

  2. Control group is assigned zero value.

Appendix 2—table 2
Comparison of models for the virological markers (within-group analysis).
Dependent variableGroupExplanatory variablesIntercept only (repeated-measures design)Intercept + time point
Estimate (95% CI)p*Estimate (95% CI)p
Total DNA
(log10 fold change from baseline)
Intensified groupIntercept−0.21 (−0.33 to −0.08)0.0022−0.17 (−0.44 to 0.11)0.0016
Time point---0.40
Control groupIntercept0.16 (0.05 to 0.27)0.00530.24 (0.01 to 0.48)0.019
Time point---0.42
US RNA
(log10 fold change from baseline)
Intensified groupIntercept−0.54 (−0.75 to −0.33)6.0 × 10–5−0.62 (−1.19 to −0.06)0.0069
Time point---0.60
Control groupIntercept0.05 (−0.26 to 0.37)0.730.24 (−0.44 to 0.92)0.71
Time point---0.88
US RNA/total DNA ratio
(log10 fold change from baseline)
Intensified groupIntercept−0.47 (−0.67 to −0.27)2.2 × 10–4−0.53 (−1.09 to 0.03)0.029
Time point---0.59
Control groupIntercept0.01 (−0.29 to 0.32)0.94−0.07 (−0.76 to 0.62)0.81
Time point---0.86
  1. *

    p values were calculated by type III tests of fixed effects.

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Céline Fombellida-Lopez
  2. Aurelija Valaitienė
  3. Lee Winchester
  4. Nathalie Maes
  5. Patricia Dellot
  6. Céline Vanwinge
  7. Aurélie Ladang
  8. Etienne Cavalier
  9. Fabrice Susin
  10. Dolores Vaira
  11. Marie-Pierre Hayette
  12. Catherine Reenaers
  13. Michel Moutschen
  14. Courtney V Fletcher
  15. Alexander O Pasternak
  16. Gilles Darcis
(2026)
Doubling dolutegravir dosage reduces the viral reservoir in ART-treated people with HIV
eLife 14:RP106931.
https://doi.org/10.7554/eLife.106931.3