BRIEF RESEARCH REPORT article

Front. Cell. Infect. Microbiol., 14 March 2024

Sec. Clinical Infectious Diseases

Volume 14 - 2024 | https://doi.org/10.3389/fcimb.2024.1340610

Plasma metabolomic profile is near-normal in people with HIV on long-term suppressive antiretroviral therapy

  • 1. Unidad de Infección Viral e Inmunidad, Centro Nacional de Microbiología (CNM), Instituto de Salud Carlos III (ISCIII), Madrid, Spain

  • 2. Centro de Investigación Biomédica en Red en Enfermedades Infecciosas (CIBERINFEC), Instituto de Salud Carlos III (ISCIII), Madrid, Spain

  • 3. Unidad de Enfermedades Infecciosas/VIH, Hospital General Universitario “Gregorio Marañón”, Madrid, Spain

  • 4. Instituto de Investigación Sanitaria Gregorio Marañón (IiSGM), Madrid, Spain

  • 5. Servicio de Medicina Interna-Unidad de VIH, Hospital Universitario La Paz, Madrid, Spain

  • 6. Instituto de Investigación Sanitaria La Paz (IdiPAZ), Madrid, Spain

  • 7. Centre of Metabolomics and Bioanalysis (CEMBIO), Facultad de Farmacia, Universidad San Pablo-CEU, CEU Universities, Boadilla del Monte, Spain

Abstract

Background:

Combination antiretroviral therapy (ART) has transformed human immunodeficiency virus (HIV) infection in people with HIV (PWH). However, a chronic state of immune activation and inflammation is maintained despite achieving HIV suppression and satisfactory immunological recovery. We aimed to determine whether the plasma metabolomic profile of PWH on long-term suppressive ART and immunologically recovered approximates the normality by comparison with healthy controls with similar age and gender.

Methods:

We carried out a cross-sectional study in 17 PWH on long-term ART (HIV-RNA <50 copies/mL, CD4+ ≥500 cells/mm3, and CD4+/CD8+ ≥1) and 19 healthy controls with similar age and gender. Metabolomics analysis was performed by gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS). The statistical association analysis was performed by principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and Generalized Linear Models (GLM) with a gamma distribution (log-link). Significance levels (p-value) were corrected for multiple testing (q-value).

Results:

PCA and PLS-DA analyses found no relevant differences between groups. Adjusted GLM showed 14 significant features (q-value<0.20), of which only three could be identified: lysophosphatidylcholine (LysoPC) (22:6) (q-value=0.148), lysophosphatidylethanolamine (LysoPE) (22:6) (q-value=0.050) and hydroperoxy-octadecatrienoic acid (HpOTrE)/dihydroperoxy-octadecatrienoic acid (DiHOTrE)/epoxy-octadecadienoic acid (EpODE) (q-value=0.136). These significant identified metabolites were directly correlated to plasma inflammatory biomarkers in PWH and negatively correlated in healthy controls.

Conclusion:

PWH on long-term ART have a metabolomic profile that is almost normal compared to healthy controls. Nevertheless, residual metabolic alterations linked to inflammatory biomarkers persist, which could favor the development of age-related comorbidities among this population.

Introduction

People with human immunodeficiency virus (HIV) (PWH) on combination antiretroviral therapy (ART) usually achieve undetectable viral load and CD4+ ≥500 cells/mm3, decreasing the risk of developing acquired immune deficiency syndrome (AIDS)-related events, mortality, and transforming HIV infection into a chronic disease (). However, ART does not eradicate HIV from the body, promoting a chronic state of immune activation and inflammation that leads to the development of non-AIDS comorbidities despite viral suppression (). Besides, PWH on suppressive ART generally presents persistently elevated CD8+ T cell counts and a low CD4+/CD8+ ratio, linked to a higher immune dysfunction (immune activation, inflammation, immunosenescence, among others), viral reservoir size, aging, comorbidities, and mortality (; ; ). Therefore, besides CD4+ ≥500 cells/mm3, CD4+/CD8+ ratio ≥1 is another goal for immune reconstitution in PWH on suppressive ART (). In this context, previous studies have reported that only one-third of PWH on long-term ART achieve a restoration of the CD4+/CD8+ ratio, its normalization being slow (; ).

Metabolomics is based on studying and analyzing metabolites and metabolic pathways involved in a specific process. Serum and plasma derived from PWH revealed altered metabolites involved in lipid and mitochondrial pathways as well as fatty acids and organic acids (). The ability of ART to rectify HIV-induced metabolic dysregulation is unclear, and a robust characterization of the metabolic alterations experienced is needed to determine the effect of ART on these pathways. In this regards, Peltenburg et al. found increased lipid metabolites in PWH after 12 months of ART (). Regarding younger cohorts of PWH, significant changes in the levels of several metabolites were found between HIV untreated patients, HIV patients on ART, and healthy controls (). In PWH on long-term suppressive ART, reported plasma metabolomic abnormalities related to amino acids and energy metabolism, urea, and tricarboxylic acid cycle compared to healthy controls. Their findings also showed alterations in the lipid complex, which could be markers of inflammation, oxidative stress, and immune cell function. have also reported an independent association between HIV infection and hypertension, hypertriglyceridemia, and abdominal obesity. Nevertheless, they included PWH with a CD4+/CD8+ ratio <1, which could play a significant role in the differences. However, to our knowledge, no previous studies have evaluated metabolomic dysregulation in PWH on long-term ART with CD4+/CD8+ ≥1. In this population subgroup, these studies must assess the risk of inflammaging, immunosenescence, and age-related comorbidities.

This study aimed to determine whether the plasma metabolomic profile of PWH on long-term suppressive ART and immunologically recovered approximates the normality by comparison with healthy controls with similar age and gender.

Methods

Study subjects

We carried out a cross-sectional study in PWH (n=17) on long-term suppressive ART and significant immunological recovery in two Hospitals in Madrid (Hospital Universitario “La Paz” and Hospital General Universitario “Gregorio Marañón”). The selection criteria of PWH were: i) ART with HIV viral load <50 copies/ml and CD4+ T-cell counts ≥500 cells/mm3 during more than one year before blood extraction, and ii) CD4+/CD8+ ratio recovery to normal levels (≥1) at time of blood extraction. Patients with active hepatitis B virus (HBV) or hepatitis C virus (HCV) coinfections were excluded.

To evaluate normal plasma metabolite levels, we also selected a group of age- and gender-matched healthy controls (HC-group, n=19) that were negative for HIV, HBV, and HCV.

The study was approved by the Research Ethics Committee of the Institute of Health Carlos III (CEI PI 23_2011, CEI PI 41_2020-v2) and was carried out according to the Declaration of Helsinki. Before registration, all participants signed written consent.

Clinical data and samples

Participant characteristics were collected from medical records. Peripheral blood samples were collected in EDTA tubes, and plasma samples were separated by centrifugation and stored at -80°C in the Spanish HIV HGM Biobank until use.

Non-targeted metabolomics

The list of reagents and standards, metabolite extraction, and sample preparation are available in Appendix A. Metabolomic analysis was performed by two complementary analytical platforms: gas chromatography–mass spectrometry (GC-MS) system (Agilent Technologies 7890A) and liquid chromatography–mass spectrometry (LC-MS) (LC: 1290 infinity II Agilent, MS: Agilent 6550 iFunnel). Detailed methods can be found in Appendix A.

In GC-MS, the deconvolution and identification were performed using MassHunter Quantitative Unknowns Analysis (B.07.00, Agilent), alignment with MassProfiler Professional software (version 13.0, Agilent), and peak integration using MassHunter Quantitative Analysis (version B.07.00, Agilent). In LC-MS, the Molecular Feature Extraction and the Recursive Feature Extraction algorithms in the MassHunter Profinder software (B.08.00, Agilent) were used for the deconvolution and alignment of the raw data. After data reprocessing, the metabolic features were filtered (full description in Supplementary File - Appendix A).

Multiplex immunoassays and ELISA

ProcartaPlexTM multiplex immunoassay (Bender MedSystems GmbH, Vienna, Austria) was used to measure several plasma biomarkers according to the manufacturer’s specifications using a Luminex 200TM analyzer (Luminex Corporation, Austin, TX, United States). The plasma biomarkers measured by ELISA multiplex were anti-inflammatory/suppressor markers – interleukin 10 (IL-10), transforming growth factor-beta 1 (TGF-β1), IL-1 receptor antagonist (IL-1RA) and IL-4 –, pro-inflammatory chemokine markers – human interferon-inducible protein 10 (IP-10), monocyte chemoattractant protein-1 (MCP-1)] and IL-8 –, pro-inflammatory cytokine markers – IL-1β, IL- 18, IL-6, tumor necrosis factor-alpha (TNF-α), interferon-gamma (IFN-γ), IL-12p70, IL-2 and IL-17A –, endothelial dysfunction markers – soluble vascular cell adhesion molecule-1 (sVCAM-1), soluble intercellular adhesion molecule-1 (sICAM-1) and soluble tumor necrosis factor receptor-1 (sTNFR-1) –, and coagulopathy markers – D-Dimer and plasminogen activator inhibitor-1 (PAI-1)–.

A Commercial ELISA was used to measure bacterial translocation markers – sCD14 and fatty acid-binding protein 2 (FABP2) (Raybiotech, Georgia, USA) and lipopolysaccharide-binding protein (LBP) (R&D Systems, Minneapolis, USA) – and the anti-transforming growth factor beta 1 (TGF-β1; Bender MedSystems GmbH, Vienna, Austria) as the multiplex immunoassay was not available. The lipopolysaccharide (LPS; Hycult Biotech, Uden, The Netherlands) was evaluated by a Limulus amebocyte lysate (LAL) chromogenic endpoint ELISA.

Statistical analysis

For the group description, variables were expressed as median [25th; 75th percentile] for continuous and as absolute numbers [percentage] for categorical data. The Mann–Whitney U and Chi-square tests were used to analyze continuous and categorical variables, respectively.

For the metabolomics analysis, variables from GC-MS and LC-MS were log-transformed (log10) and auto-scaled to make individual features more comparable. Next, we performed an unsupervised analysis by principal component analysis (PCA) and a supervised analysis by partial least squares discriminant analysis (PLS-DA) for features detected in GC-MS and LC-MS [positive and negative electrospray ionization (ESI)]. The optimal number of PLS-DA components was determined with the leave-one-out cross-validation (LOOCV) method, using R2 and Q2 values as performance measures. Permutation was carried out by separation distance (B/W) with a permutation number of 1000 to confirm the model’s validity.

Generalized Linear Models (GLM) with gamma distribution (log-link) were used to independently analyze the differences between the study groups for each metabolite. This test provides the arithmetic mean ratio (AMR) and its significance level (p-value), which was corrected for multiple testing using the False Discovery Rate (FDR) with Benjamini and Hochberg procedure (q-value). Additionally, GLM models were adjusted by baseline characteristics (age, gender, and body mass index) previously selected by a stepwise method by the Akaike information criterion (AIC) (forward, p<0.05; q-value<0.20).

Correlation between significant metabolites and plasma biomarkers was performed using the Spearman correlation test. Those suitable correlations (r>0.5 or r<-0.5) and a significance value (p<0.05; q-value<0.20) were considered relevant.

The statistical analysis was done with MetaboAnalyst 4.0 software (http://www.metaboanalyst.ca/) and R statistical package version v3.5.1 (R Foundation for Statistical Computing, Vienna, Austria).

Metabolite identification

The significant metabolites (q-value<0.2) were identified. In GC-MS, the identification was made based on FiehnLib () and NIST 14 libraries. In LC-MS, the list of accurate masses was searched using the CEU Mass Mediator search tool (http://ceumass.eps.uspceu.es/; error ± 5 ppm) to obtain tentative identifications. Each of them were manually curated based on their MS adducts (; ). In the cases that it was applicable, the elution order was also considered to discard spurious identifications. Eventually, the biological role of each compound was evaluated, and unrelated identifications such as pesticides, drugs, or not possible chemical structures were excluded. The metabolites are reported in agreement with the criteria of the Metabolomics Standards Initiative (10.1007/s11306-007-0070-6) with a confidence level grade 2 (putatively annotated compounds), which certitude is increased after manual curation of the final list.

Results

Patient characteristics

The epidemiological and clinical data of participants are shown in Table 1. In brief, the median age of PWH was 57 years, 64.7% were males, 47.1% had prior AIDS diagnosis, and the median time on ART was 10.7 years. Although all PWH had a CD4+/CD8+ ratio ≥1, healthy controls had significantly higher values (p-value<0.001).

Table 1

CharacteristicHealthy controlPWHp-value
No.1917
Age (years)56 (51.5; 58.5)57 (55.0; 58.0)0.656
Gender (male)10 (52.6%)11 (64.7%)0.693
BMI (kg/m2)24.9 (23.8; 27.1)25.6 (22.8; 27.9)0.634
Comorbidities
Obesity (BMI>30)2 (10.5%)3 (17.6%)0.889
Arterial hypertension2 (10.5%)1 (5.9%)0.999
Diabetes0 (0.0%)1 (5.9%)0.999
Risk group of HIV infection
Heterosexual8 (47%)
Homosexual9 (53%)
Time of HIV infection (years)15.1 (10.9; 21.4)
Antiretroviral therapy
PI-based5 (29.4%)
2NRTI+II-based1 (5.9%)
2NRTI+NNRTI-based9 (52.9%)
Others2 (11.8%)
Antidepressant treatment2 (11.8%)
Antihypertensive treatment3 (17.6%)
Time on cART (years)10.7 (6.9; 16.3)
Lymphocytes counts
CD4+ T-cells (%)45.6 (42.9; 48.3)37.7 (34.9; 43.3)0.003
CD8+ T-cells (%)15.7 (12.5; 21.1)23.5 (20.4; 28.6)0.002
CD4+/CD8+2.7 (2.2; 3.6)1.6 (1.2; 2.1)<0.001
CD4+/CD8+ ≥119 (100%)17 (100%)0.999
CD4+ T-cells/mm3977 (804; 1062)
CD4+ ≥500 cells/mm317 (100%)
HIV markers
Prior AIDS diagnosis8 (47.1%)
Nadir CD4+ T-cells/mm3309 (64; 402)
Nadir CD4+ ≤200 T-cells/mm37 (41.2%)
uVL (HIV-RNA <50 copies/mL)17 (100%)

Clinical and epidemiological characteristics of PLW on long-term suppressive ART and healthy controls.

Statistics: Values are expressed as absolute number (percentage) and median (interquartile range). P-values were calculated by Chi-square tests, Fisher's Exact Test or Mann-Whitney tests. Statistically significant differences are shown in bold.

Abbreviations: BMI, body mass index; HIV, human immunodeficiency virus; HIV-RNA, HIV plasma viral load; AIDS, acquired immune deficiency syndrome; NNRTI, non-nucleoside analogue HIV reverse transcriptase inhibitor; NRTI, nucleoside analogue HIV reverse transcriptase inhibitor; PI, protease inhibitor; II, integrase inhibitor; PWH, people with HIV; uVL, undetectable HIV viral load (<50 copies/mL).

Reliability analysis

PCA indicated that quality control (QC) samples were tightly clustered together in the center of the plot, thus validating the signal stability and technical reproducibility (Supplementary Figure 1).

Metabolite association analysis

PCA showed similarity between the sample groups for all the platforms used (Figure 1A). PLS-DA was performed for features detected in GC-MS (R2 = 0.377 and Q2 = 0.002; one component), LC-MS ESI+ (R2 = 0.658 and Q2=-0.076; three components), and LC-MS ESI- (R2 = 0.884 and Q2 = 0.232; five components) (Figure 1B). However, PLS-DA could not be validated by permutation for any of the platforms: GC-MS (p=0.323), LC-MS ESI+ (p=0.221), and LC-MS ESI- (p=0.548) (Supplementary Figure 2). Therefore, PCA and PLS-DA showed no significant differences between the two study groups. Similarly, differences between types of ART (PI-based, 2NRTI+NNRTI-based, and others) (Supplementary Figure 3) and nadir CD4+ levels (<200 cells/mm3 and ≥200 cells/mm3) (Supplementary Figure 4) in PWH group were not found.

Figure 1

GLM analysis adjusted by the most relevant covariates showed 63 significant features (p<0.05), of which 14 had a q<0.20 after correcting by FDR (Table 2). Of these, identification data were only obtained for three metabolites. Briefly, while lysophosphatidylcholine [LysoPC (22:6)] and lysophosphatidylethanolamine [LysoPE (22:6)] showed reduced levels [aAMR=0.59 (p=0.005; q=0.148) and aAMR=0.68 (p=0.001; q=0.050), respectively], an oxidized lipid not fully identified had increased levels [aAMR=1.11 (p=0.004; q=0.136)] among PWH compared to healthy controls (Supplementary Figure 5). The possible tentative identifications for this oxidized lipid were: 12-hydroperoxy-octadecatrienoic acid (HpOTrE), 13-HpOTrE, 13S-HpOTrE, 15,16-epoxy-octadecadienoic acid (EpODE), 16-HpOTrE, 7,8-dihydroperoxy-octadecatrienoic acid (DiHOTrE), 9H-12(13)-EpODE, 9-HpOTrE, or 9S-HpOTrE (HpOTrE/DiHOTrE/EpODE).

Table 2

FeatureTechnologyMassRT (min)aAMRIC2.5IC97.5pqIdentification
Oleic acidGC-MS33920.430.650.450.930.0240.310Oleic acid
Palmitic acidGC-MS31318.840.710.530.960.0320.310Palmitic acid
Palmitoleic acidGC-MS31118.650.600.380.950.0380.310Palmitoleic acid
p-CresolGC-MS1658.190.640.430.960.0390.310p-Cresol
Threonic acidGC-MS29213.550.670.500.890.0100.271Threonic acid
ThreonineGC-MS21811.350.820.680.990.0460.316Threonine
Unknown_24.7GC-MS41524.700.500.300.830.0120.271
Unknown_7.04GC-MS897.040.500.260.950.0400.310
Unknown_7.18GC-MS897.180.490.250.940.0390.310
Unknown_8.108GC-MS1178.110.620.450.850.0060.271
414.2041/0.22399998LC-MS ESI+414.20410.221.181.031.350.0240.270
352.2021/0.25400043LC-MS ESI+352.20210.252.391.344.270.0060.148Unknown
416.2075/0.22399998LC-MS ESI+416.20750.221.271.071.510.0100.194Unknown
211.1932/6.314987LC-MS ESI+211.19326.311.121.011.240.0460.307
519.3325/5.4370084LC-MS ESI+519.33255.440.730.580.940.0190.270
267.2559/6.318981LC-MS ESI+267.25596.321.141.021.280.0330.289
297.268/0.2760002LC-MS ESI+297.2680.281.141.031.260.0150.240
103.0999/5.572011LC-MS ESI+103.09995.570.760.600.950.0230.270
519.333/5.573987LC-MS ESI+519.3335.570.730.550.970.0400.295
442.1838/0.26800057LC-MS ESI+442.18380.271.231.061.440.0120.220
307.2872/6.5200105LC-MS ESI+307.28726.521.191.031.370.0250.270
567.3322/5.5489993LC-MS ESI+567.33225.550.590.420.830.0050.148LysoPC(22:6)
326.2033/0.25500023LC-MS ESI+326.20330.261.131.041.240.0080.188Unknown
269.2718/6.687009LC-MS ESI+269.27186.691.171.021.350.0300.285
295.287/6.8239846LC-MS ESI+295.2876.821.161.021.330.0280.281
796.155/9.532984LC-MS ESI+796.1559.531.431.022.020.0450.307
729.2363/8.525017LC-MS ESI+729.23638.531.341.031.740.0380.295
798.1522/9.532984LC-MS ESI+798.15229.531.491.052.130.0340.289
806.2534/9.532984LC-MS ESI+806.25349.531.401.041.880.0320.289
117.0785/11.895997LC-MS ESI+117.078511.901.111.051.16<0.0010.027Unknown
354.0627/10.401996LC-MS ESI+354.062710.401.331.021.720.0420.300
814.2073/10.39901LC-MS ESI+814.207310.401.431.061.930.0240.270
309.3033/7.241983LC-MS ESI+309.30337.241.181.011.360.0400.295
320.2461/0.27700037LC-MS ESI+320.24610.281.341.071.680.0150.240
283.2871/7.157991LC-MS ESI+283.28717.161.181.021.360.0320.289
946.1904/11.094979LC-MS ESI+946.190411.091.571.042.380.0400.295
944.1871/11.094977LC-MS ESI+944.187111.091.371.051.770.0240.270
281.2718/7.0589924LC-MS ESI+281.27187.061.211.011.460.0450.307
360.2254/0.26599964LC-MS ESI+360.22540.271.081.021.140.0100.194Unknown
103.0996/5.4420066LC-MS ESI+103.09965.440.730.610.890.0030.135Unknown
140.1059/11.895014LC-MS ESI+140.105911.901.141.051.240.0060.148Unknown
662.4444/11.645997LC-MS ESI+662.444411.651.131.011.270.0490.307
404.2528/0.28799918LC-MS ESI+404.25280.291.251.041.500.0210.270
877.2731/10.406988LC-MS ESI+877.273110.411.501.241.82<0.0010.027Unknown
525.2854/5.403004LC-MS ESI+525.28545.400.680.550.840.0010.050LysoPE(22:6)
879.2761/10.407983LC-MS ESI+879.276110.411.601.262.02<0.0010.027Unknown
879.2739/10.407982LC-MS ESI+879.273910.411.551.291.86<0.001<0.001Unknown
870.1751/10.402992LC-MS ESI+870.175110.41.361.051.770.0270.280
340.2399/7.2280173LC-MS ESI-340.23997.230.230.060.920.0450.708
177.0801/0.25799963LC-MS ESI-177.08010.260.410.200.850.0220.636
188.0143/0.24000052LC-MS ESI-188.01430.240.500.300.820.0090.397
108.0572/0.24200036LC-MS ESI-108.05720.240.550.360.830.0080.397
886.5581/11.944002LC-MS ESI-886.558111.942.831.176.860.0270.636
251.1548/6.224987LC-MS ESI-251.15486.220.230.060.840.0330.636
404.2717/7.6960144LC-MS ESI-404.27177.700.230.090.610.0060.397
889.5749/11.930998LC-MS ESI-889.574911.932.221.044.740.0470.708
499.9411/0.257LC-MS ESI-499.94110.260.500.360.69<0.0010.045Unknown
884.5422/11.938984LC-MS ESI-884.542211.942.551.115.810.0330.636
506.3392/8.925018LC-MS ESI-506.33928.930.270.090.880.0360.636
240.0729/7.214998LC-MS ESI-240.07297.210.430.230.820.0160.545
300.2091/6.3210077LC-MS ESI-300.20916.320.180.060.550.0050.397
713.4483/6.718995LC-MS ESI-713.44836.721.611.052.470.0360.636
310.2146/0.25799963LC-MS ESI-310.21460.261.111.041.180.0040.136HpOTrE/DiHOTrE/
EpODE *

Association of individual metabolites with HIV infection among PWH on long-term ART and immunologically recovered compared to healthy controls.

Statistics: Generalized Linear Models (GLM) with a gamma distribution (log-link) (dependent variable: plasma metabolites; independent variable: HIV-infection), adjusted by epidemiological characteristics (age, gender and body mass index). P-values were adjusted by FDR correction for multiple comparisons (Benjamini and Hochberg). Statistically significant differences are shown in bold.

Abbreviations: RT, retention time; aAMR, adjusted ratio of the arithmetic means; CI, confidence interval; p-value, level of significance; q-value, adjusted p-value by FDR correction; GC-MS, gas chromatography–mass spectrometry; LC-MS ESI+, liquid chromatography–mass spectrometry, positive electrospray ionization; LC-MS ESI+, liquid chromatography–mass spectrometry, negative electrospray ionization; LysoPC, lysophosphocoline; LysoPE, lysophosphatidylethanolamine; HpOTrE, hydroperoxy-octadecatrienoic acid; DiHOTrE, dihydroperoxy-octadecatrienoic acid; EpODE, epoxy-octadecadienoic acid.

* 12-HpOTrE / 13-HpOTrE / 13S-HpOTrE / 15,16-EpODE / 16-HpOTrE / 7,8-DiHOTrE / 9H-12(13)-EpODE / 9-HpOTrE / 9S-HpOTrE.

Correlation analysis between metabolites and plasma biomarkers

The plasma biomarkers concentrations in both HC and PWH groups are shown in Supplementary Table 1. The correlations between significant identified metabolites and plasma biomarkers are shown in Figure 2 (full description in Supplementary Tables 2, 3). Several correlations were found significant, even after FDR correction (r>0.5 or r<-0.5; p<0.05; q-value<0.20). LysoPC (22:6) was negatively correlated with MCP-1 in PWH and the HC-group (p=0.037 and p=0.002, respectively). Besides, while no significant correlations were found for LysoPE (22:6) in PWH, negative correlations were found between LysoPE (22:6) and IL-12p70 (p=0.012), IL-17A (p=0.040), and sCD14 (p=0.002) in the HC-group. HpOTrE/DiHOTrE/EpODE was also positively correlated with sICAM-1 (p=0.004) and sTNFR-I (p=0.009) in PWH. However, we found negative correlations between HpOTrE/DiHOTrE/EpODE and sICAM-1 (p=0.001), sVCAM-1 (p=0.016), IL-12p70 (p=0.018), IL-1β (p=0.011), IL-8 (p=0.027), IL-4 (p=0.021) in the HC-group.

Figure 2

Discussion

We found little differences in the metabolic profile between PWH with immunological recovery after long-term suppressive ART and healthy controls.

PCA and PLS-DA multivariate analysis showed no relevant results for any platforms used. Previous metabolomic studies have shown a clear separation between PWH and HC groups (; ), even for PWH on long-term successful antiretroviral therapy (), in which alterations in amino-acid levels, energetics, and lipids have been found. However, PWH did not achieve a CD4/CD8 ratio >1 in most studies, which could contribute to the differences observed in these studies in contrast to ours. Besides, several articles have described a metabolomic signature associated with immunological CD4+ T-cell recovery after long-term of antiretroviral therapy. However, it has been also studied in PWH whose immunological recovery did not exceed the CD4 T-cell count of 500 (; ; ), which limits comparisons with our study. In this sense, different overlapping metabolic profiles between PWH and HC has been found probably due to different disease stage of individuals (), which indicates that the characteristics of patients and their level of immunological recovery are crucial to interpret the findings. Thus, further studies including PWH with longer periods of successful ART and improved immune reconstitution (CD4+/CD8+ ≥1) would be needed to corroborate the near-normal metabolomic profile found in our study.

However, although, to our knowledge, no previous metabolomic studies have been performed in PWH on long-term ART with CD4+/CD8+ ≥1compared to healthy controls, the finding of an almost normalization of the metabolic profile in this subgroup of PWH is concordant with previous studies, in which a normalization of different immune-related molecular markers has been described. In this same cohort, Brochado-Kith et al. showed that peripheral blood mononuclear cells gene expression and peripheral blood biomarkers in PWH, with normalized CD4+/CD8+ ratio, had a similar profile compared to healthy controls (). Serrano-Villar et al. found that PWH on ART with a normalized CD4+/CD8+ ratio demonstrated traits of a nearly healthy immune system (). Sperk et al. showed that some pro-inflammatory cytokines and chemokines return to healthy levels in PWH with nearly twenty years of ART ().

Additionally, the GLM analysis of each metabolite showed scarce differences between groups. Decreased levels of lysoPCs (22:6) and lysoPE (22:6) were found in PWH on long-term ART compared to healthy controls. Although these specific metabolites have not been described in previous studies, an altered level of phosphatidylcholine and phosphatidylethanolamine has been found comparing PWH on long-term successful antiretroviral therapy and HC groups (). Likewise, Lu et al. described that glycerophospholipid metabolism was one of the pathways with highest impact (), which is in line with our findings. Additionally, lower level of these lysophospholipids has been also associated to other human diseases, such as metabolic, cardiovascular, and neurodegenerative disorders, all non-AIDS-defining events (NADEs) among PWH (). Likewise, lower lysoPCs (22:6) and lysoPE (22:6) levels have also been associated with more advanced cirrhosis stages among HIV/HCV-coinfected patients (). In addition, decreased concentrations of different LysoPCs species have been associated with the risk of obesity (), linked to inflammaging, insulin resistance, metabolic syndrome, and non-alcoholic fatty liver disease, among others. Regarding ART, several diabetes-associated lipid species are perturbed in ART-treated PWH, as ART disrupts lipid metabolism ().

The metabolites whose levels differed between groups have been associated with inflammaging and immune activation (; ; ), which can lead to premature aging and NADEs (). In this setting, we analyzed the correlation of these metabolites with several plasma biomarkers and found a significant negative correlation between LysoPC (22:6) and the pro-inflammatory chemokine MCP-1 in PWH, supporting the inflammaging and immune activation state in these patients. No significant correlations were found for LysoPE (22:6) in PWH. Regarding the oxidized lipid, significant positive correlations were found with sICAM-1 and sTNFR-1 in PWH, both being endothelial dysfunction markers and associated with increased risk of cardiovascular disease, cancers, and atherosclerosis, among others (; ; ). Interestingly, phospholipid metabolism has been previously found to be altered in studies including younger PWH on long-term ART compared HC (mean age of 45 years in both groups) (; ), which could indicate that similar molecular mechanisms may also occur in younger cohorts, although further studies are needed.

Although the introduction of ART has increased the life expectancy of PWH (), our study suggests it does not restore health ad-integrum. We observed that even patients with a normalized CD4+/CD8+ ratio presented alteration of specific metabolites that could be involved in the pathogenesis of different age-related comorbidities due to chronic immune activation, immunosenescence, and inflammaging.

Therefore, it is essential to carry out further studies that corroborate the role of metabolic changes during ART among PWH. Likewise, metabolomics can offer an alternative view of the inflammatory state of patients.

Some limitations should be considered for a correct interpretation of the data. Firstly, the sample size was limited, which could have restricted the statistical power to detect metabolic differences between groups. In addition, the modest sample size may also increase the false positive rate, but our positive findings were FDR-corrected, lending robustness to our results. Second, this study has a cross-sectional design, which may introduce some bias and limit the interpretation of our findings. Thirdly, clinical data related to the individual cognitive status were unavailable, which would be interesting as neurocognitive disorders have been shown to impact the lipidome. Finally, more studies would be needed in other patient cohorts, such as among PWH with different age ranges, to confirm whether similar results are found.

Conclusions

In conclusion, our data suggests that PWH on long-term ART, with CD4+/CD8+ ratio ≥1, have a metabolomic profile that is almost normal compared to healthy controls. Nevertheless, residual metabolic alterations linked to inflammatory biomarkers persist, which could favor the development of age-related comorbidities among this population.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by The study was approved by the Research Ethics Committee of the Institute of Health Carlos III (CEI PI 23_2011, CEI PI 41_2020-v2) and was carried out according to the Declaration of Helsinki. Before registration, all participants signed written consent. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

AV-B: Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft. RM-E: Data curation, Formal analysis, Investigation, Visualization, Writing – original draft. JB: Conceptualization, Investigation, Writing – review & editing. JG-G: Investigation, Writing – review & editing. OB-K: Formal analysis, Investigation, Writing – review & editing. DR: Investigation, Writing – review & editing. AF-R: Investigation, Writing – review & editing. LP-L: Investigation, Writing – review & editing. VH: Investigation, Writing – review & editing. CB: Investigation, Writing – review & editing. SR: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing. MAJ-S: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Instituto de Salud Carlos III (ISCIII; grant numbers CP17CIII/00007, PI18CIII/00028 and PI21CIII/00033 to MAJS, PI17/00657 and PI20/00474 to JB, PI17/00903 and PI20/00507 to JG-G, PI18CIII/00020 to AF-R, and PI17CIII/00003 and PI20CIII/00004 to SR) and Ministerio de Ciencia e Innovación (AEI, PID2021-126781OB-I00 to AF-R). CB and DR acknowledge funding from the Ministerio de Ciencia e Innovación (RTI2018-095166-B-I00). The study was also funded by the CIBER -Consorcio Centro de Investigación Biomédica en Red- (CB 2021), Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación and Unión Europea – NextGenerationEU (CB21/13/00044). MAJ-S is Miguel Servet researcher supported and funded by ISCIII (grant numbers CP17CIII/00007). RME is César Nombela researcher supported and funded by Comunidad de Madrid (CAM) (grant number 2023-T1/SAL-GL28980).

Acknowledgments

This study would not have been possible without the collaboration of all the patients, medical and nursery staff, and data managers who have taken part in the project. We want to acknowledge the patients in this study for their participation and the HIV BioBank integrated into the Spanish AIDS Research Network and collaborating Centres (http://hivhgmbiobank.com/donor-area/hospitals-and-centres-transferring-samples/?lang=en) for the generous gifts of clinical samples used in this work. The HIV BioBank, integrated into the Spanish AIDS Research Network, is partially funded by the RD16/0025/0019 project as part of the Plan Nacional R + D + I and cofinanced by ISCIII- Subdirección General de Evaluación and Fondo Europeo de Desarrollo Regional (FEDER).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2024.1340610/full#supplementary-material

Abbreviations

AIDS, Acquired immune deficiency syndrome; AMR, Arithmetic mean ratio; ART, Combination antiretroviral therapy; DiHOTrE, Dihydroperoxy-octadecatrienoic acid; FDR, False discovery rate; GC-MS, Gas chromatography–mass spectrometry; GLM, Generalized Linear Models; EpODE, Epoxy-octadecadienoic acid; HC-group, Healthy control group; HIV, Human immunodeficiency virus; HpOTrE, Hydroperoxy-octadecatrienoic acid; IL, Interleukin; LOOCV, Leave-one-out cross-validation; LC-MS, Liquid chromatography–mass spectrometry; LysoPC, Lysophosphatidylcholine; LysoPE, Lysophosphatidyletanolamine; MCP-1, Monocyte chemoattractant protein-1; PLS-DA, Partial least squares discriminant analysis; PWH, People with HIV; PCA, Principal component analysis; sCD14, Soluble cluster of differentiation 14; sICAM-1, Soluble intercellular adhesion molecule-1;sVCAM-1, Soluble vascular cell adhesion molecule-1.

References

  • 1

    BabuH.SperkM.AmbikanA. T.RachelG.ViswanathanV. K.TripathyS. P.et al. (2019). Plasma metabolic signature and abnormalities in HIV-infected individuals on long-term successful antiretroviral therapy. Metabolites9 (10), 210. doi: 10.3390/metabo9100210

  • 2

    BarberM. N.RisisS.YangC.MeikleP. J.StaplesM.FebbraioM. A.et al. (2012). Plasma lysophosphatidylcholine levels are reduced in obesity and type 2 diabetes. PloS One7 (7), e41456. doi: 10.1371/journal.pone.0041456

  • 3

    BroChado-KithO.MartinezI.BerenguerJ.MedranoL. M.Gonzalez-GarciaJ.Garcia-BroncanoP.et al. (2020). Near normalization of peripheral blood markers in HIV-infected patients on long-term suppressive antiretroviral therapy: a case-control study. AIDS34 (13), 18911897. doi: 10.1097/QAD.0000000000002645

  • 4

    BrownT. T.GlesbyM. J. (2011). Management of the metabolic effects of HIV and HIV drugs. Nat. Rev. Endocrinol.8 (1), 1121. doi: 10.1038/nrendo.2011.151

  • 5

    BuiT. M.WiesolekH. L.SumaginR. (2020). ICAM-1: A master regulator of cellular responses in inflammation, injury resolution, and tumorigenesis. J. Leukoc. Biol.108 (3), 787799. doi: 10.1002/JLB.2MR0220-549R

  • 6

    CabyF. (2017). CD4+/CD8+ ratio restoration in long-term treated HIV-1-infected individuals. AIDS31 (12), 16851695. doi: 10.1097/QAD.0000000000001533

  • 7

    CassolE.MisraV.HolmanA.KamatA.MorgelloS.GabuzdaD. (2013). Plasma metabolomics identifies lipid abnormalities linked to markers of inflammation, microbial translocation, and hepatic function in HIV patients receiving protease inhibitors. BMC Infect. Dis.13, 203. doi: 10.1186/1471-2334-13-203

  • 8

    ChenJ.LiuX.QinS.RuanG.LuA.ZhangJ.et al. (2022). A novel prognostic score including the CD4/CD8 for AIDS-related lymphoma. Front. Cell Infect. Microbiol.12, 919446. doi: 10.3389/fcimb.2022.919446

  • 9

    GelpiM.AfzalS.LundgrenJ.RonitA.RoenA.MocroftA.et al. (2018). Higher Risk of Abdominal Obesity, Elevated Low-Density Lipoprotein Cholesterol, and Hypertriglyceridemia, but not of Hypertension, in People Living With Human Immunodeficiency Virus (HIV): Results From the Copenhagen Comorbidity in HIV Infection Study. Clin. Infect. Dis.67 (4), 579586. doi: 10.1093/cid/ciy146

  • 10

    GodzienJ.CiborowskiM.ArmitageE. G.JorgeI.CamafeitaE.BurilloE.et al. (2016). A single in-vial dual extraction strategy for the simultaneous lipidomics and proteomics analysis of HDL and LDL fractions. J. Proteome Res.15 (6), 17621775. doi: 10.1021/acs.jproteome.5b00898

  • 11

    GrossM. D.BielinskiS. J.Suarez-LopezJ. R.ReinerA. P.BaileyK.ThyagarajanB.et al. (2012). Circulating soluble intercellular adhesion molecule 1 and subclinical atherosclerosis: the Coronary Artery Risk Development in Young Adults Study. Clin. Chem.58 (2), 411420. doi: 10.1373/clinchem.2011.168559

  • 12

    HanW. M.ApornpongT.KerrS. J.HiransuthikulA.GatechompolS.DoT.et al. (2018). CD4/CD8 ratio normalization rates and low ratio as prognostic marker for non-AIDS defining events among long-term virologically suppressed people living with HIV. AIDS Res. Ther.15 (1), 13. doi: 10.1186/s12981-018-0200-4

  • 13

    HerbertC.LuiesL.LootsD. T.WilliamsA. A. (2023). The metabolic consequences of HIV/TB co-infection. BMC Infect. Dis.23 (1), 536. doi: 10.1186/s12879-023-08505-4

  • 14

    HewerR.VorsterJ.SteffensF. E.MeyerD. (2006). Applying biofluid 1H NMR-based metabonomic techniques to distinguish between HIV-1 positive/AIDS patients on antiretroviral treatment and HIV-1 negative individuals. J. Pharm. BioMed. Anal.41 (4), 14421446. doi: 10.1016/j.jpba.2006.03.006

  • 15

    KindT.WohlgemuthG.LeeD. Y.LuY.PalazogluM.ShahbazS.et al. (2009). FiehnLib: mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal. Chem.81 (4), 1003810048. doi: 10.1021/ac9019522

  • 16

    LawS. H.ChanM. L.MaratheG. K.ParveenF.ChenC. H.KeL. Y. (2019). An updated review of lysophosphatidylcholine metabolism in human diseases. Int. J. Mol. Sci.20 (5), 1149. doi: 10.3390/ijms20051149

  • 17

    LeeJ. M.LeeH.KangS.ParkW. J. (2016). Fatty acid desaturases, polyunsaturated fatty acid regulation, and biotechnological advances. Nutrients8 (1), 23. doi: 10.3390/nu8010023

  • 18

    LuL.YangY.YangZ.WuY.LiuX.LiX.et al. (2023). Altered plasma metabolites and inflammatory networks in HIV-1 infected patients with different immunological responses after long-term antiretroviral therapy. Front. Immunol.14, 1254155. doi: 10.3389/fimmu.2023.1254155

  • 19

    LuW.MehrajV.VybohK.CaoW.LiT.RoutyJ. P. (2015). CD4:CD8 ratio as a frontier marker for clinical outcome, immune dysfunction and viral reservoir size in virologically suppressed HIV-positive patients. J. Int. AIDS Soc.18 (1), 20052. doi: 10.7448/IAS.18.1.2005220052

  • 20

    MunshiS. U.RewariB. B.BhaveshN. S.JameelS. (2013). Nuclear magnetic resonance based profiling of biofluids reveals metabolic dysregulation in HIV-infected persons and those on anti-retroviral therapy. PloS One8 (5), e64298. doi: 10.1371/journal.pone.0064298

  • 21

    NystromS.GovenderM.YapS. H.KamarulzamanA.RajasuriarR.LarssonM. (2021). HIV-infected individuals on ART with impaired immune recovery have altered plasma metabolite profiles. Open Forum Infect. Dis.8 (7), ofab288. doi: 10.1093/ofid/ofab288

  • 22

    PeltenburgN. C.SchoemanJ. C.HouJ.MoraF.HarmsA. C.LoweS. H.et al. (2018). Persistent metabolic changes in HIV-infected patients during the first year of combination antiretroviral therapy. Sci. Rep.8 (1), 16947. doi: 10.1038/s41598-018-35271-0

  • 23

    QianS.ChenX.WuT.SunY.LiX.FuY.et al. (2021). The accumulation of plasma acylcarnitines are associated with poor immune recovery in HIV-infected individuals. BMC Infect. Dis.21 (1), 808. doi: 10.1186/s12879-021-06525-6

  • 24

    Rodriguez-GallegoE.GomezJ.PachecoY. M.PeraireJ.ViladesC.Beltran-DebonR.et al. (2018). A baseline metabolomic signature is associated with immunological CD4+ T-cell recovery after 36 months of antiretroviral therapy in HIV-infected patients. AIDS32 (5), 565573. doi: 10.1097/QAD.0000000000001730

  • 25

    SalekR. M.SteinbeckC.ViantM. R.GoodacreR.DunnW. B. (2013). The role of reporting standards for metabolite annotation and identification in metabolomic studies. Gigascience2 (1), 13. doi: 10.1186/2047-217X-2-13

  • 26

    SalgueroS.RojoD.BerenguerJ.Gonzalez-GarciaJ.Fernandez-RodriguezA.BroChado-KithO.et al. (2020). Plasma metabolomic fingerprint of advanced cirrhosis stages among HIV/HCV-coinfected and HCV-monoinfected patients. Liver Int.40 (9), 22152227. doi: 10.1111/liv.14580

  • 27

    Serrano-VillarS.SainzT.LeeS. A.HuntP. W.SinclairE.ShacklettB. L.et al. (2014). HIV-infected individuals with low CD4/CD8 ratio despite effective antiretroviral therapy exhibit altered T cell subsets, heightened CD8+ T cell activation, and increased risk of non-AIDS morbidity and mortality. PloS Pathog.10 (5), e1004078. doi: 10.1371/journal.ppat.1004078

  • 28

    SperkM.ZhangW.NowakP.NeogiU. (2018). Plasma soluble factor following two decades prolonged suppressive antiretroviral therapy in HIV-1-positive males: A cross-sectional study. Med. (Baltimore)97 (5), e9759. doi: 10.1097/MD.0000000000009759

  • 29

    ToledoE.WangD. D.Ruiz-CanelaM.ClishC. B.RazquinC.ZhengY.et al. (2017). Plasma lipidomic profiles and cardiovascular events in a randomized intervention trial with the Mediterranean diet. Am. J. Clin. Nutr.106 (4), 973983. doi: 10.3945/ajcn.116.151159

  • 30

    WilliamsA.KoekemoerG.LindequeZ.ReineckeC.MeyerD. (2012). Qualitative serum organic acid profiles of HIV-infected individuals not on antiretroviral treatment. Metabolomics8, 804818. doi: 10.1007/s11306-011-0376-2

  • 31

    ZicariS.SessaL.CotugnoN.RuggieroA.MorrocchiE.ConcatoC.et al. (2019). Immune activation, inflammation, and non-AIDS co-morbidities in HIV-infected patients under long-term ART. Viruses11 (3), 200. doi: 10.3390/v11030200

Summary

Keywords

antiretroviral therapy, CD4+/CD8+ ratio, HIV, inflammation, metabolomics

Citation

Virseda-Berdices A, Martín-Escolano R, Berenguer J, González-García J, Brochado-Kith O, Rojo D, Fernández-Rodríguez A, Pérez-Latorre L, Hontañón V, Barbas C, Resino S and Jiménez-Sousa MÁ (2024) Plasma metabolomic profile is near-normal in people with HIV on long-term suppressive antiretroviral therapy. Front. Cell. Infect. Microbiol. 14:1340610. doi: 10.3389/fcimb.2024.1340610

Received

18 November 2023

Accepted

23 February 2024

Published

14 March 2024

Volume

14 - 2024

Edited by

Alejandro Vallejo, Ramón y Cajal Institute for Health Research, Spain

Reviewed by

Anna Rull, Pere Virgili Health Research Institute (IISPV), Spain

Maria Abad Fernandez, University of North Carolina at Chapel Hill, United States

Updates

Copyright

*Correspondence: María Ángeles Jiménez-Sousa, ; Salvador Resino,

†These authors have contributed equally to this work

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics