Abstract
People with HIV (PWH) on combination antiretroviral therapy (cART) are living longer lives due to modern cART advances and increased routine medical care. The full landscape of aging with HIV is unclear; given that HIV emerged relatively recently in human history and initially had a high mortality rate, there has not been a substantially aged population to evaluate. In this study, we set out to perform high-throughput plasma analyte profiling by multiplex analysis, focusing on various T helper (Th)-related cytokines, chemokines, and proinflammatory and anti-inflammatory cytokines. The primary goals being to provide reference ranges of these analytes for aging PWH cohorts, as well as testing the utility of high-throughput multiplex plasma assays. The cohort used in this study comprised age-matched healthy donors (32.6–73.5 years of age), PWH on cART (26.7–60.2 years of age), and viremic PWH (27.5–59.4 years of age). The patients in each group were then stratified across the age span to examine age-related impacts of these plasma biomarkers. Our results largely indicate feasibility of plasma analyte monitoring by multiplex and demonstrate a high degree of person-to-person variability regardless of age and HIV status. Nonetheless, we find multiple associations with age, duration of known infection, and viral load, all of which appear to be driven by either prolonged HIV disease progression or long-term use of cART.
Introduction
Combination antiretroviral therapy (cART) has dramatically improved the prognosis of people with HIV (PWH) due to effective viral suppression and CD4+ T cell recovery. 1 Increased survival rate, along with increased incidences of HIV acquisition among people older than 50 years, is now shifting the population of PWH to an older demographic. 2 While aging alone causes an accrual of alterations throughout the body over time, which are associated with a variety of chronic diseases and potential changing susceptibility to infectious diseases, 3,4 the compounded effect of long-term HIV infection and associated chronic inflammation in PWH as they age is not fully understood. 5,6
Chronic immune dysfunction associated with HIV infection may exacerbate age-related comorbidities and contribute to early onset of aging-related chronic conditions, 7 –9 and emerging data support contribution that cART itself may influence the aging process. 7 Pre-exposure prophylaxis has been proven to be the most effective HIV prevention method for high-risk individuals such as the aging population, although this population is also at a higher risk of developing side effects from exposure to the active ingredient tenofovir. 10
With access to adequate treatment and care, the average life expectancy of PWH is now comparable to those without HIV infection, although management of comorbidities remains a complex issue. 11 –13 Biological changes in aging are often reflected in the plasma proteome, where plasma sampled at different stages of life reveals a variance in biological pathways and associations with age-related diseases and phenotypic traits. 4 Indeed, plasma proteomics is an established source for early biomarkers of systemic disease, 14 and has previously revealed higher inflammatory protein concentrations within PWH on long-term cART compared to healthy donors (HD), including mucosal defense chemokines. 15
Both viremic and virally suppressed PWH have shown age-related monocyte activation markers comparable to older HD. 16 To determine whether immune-related proteins are impacted by HIV across the age spectrum, we conducted a high-throughput multiplex immunoassay 17 to measure cytokines and chemokines in the plasma from HIV-uninfected donors, PWH with viral suppression by cART, and viremic PWH over a wide age distribution.
Methods
Study participant characteristics
Cryopreserved human plasma samples were obtained from the Hawai'i Aging with HIV-1 Cohort (HAHC) at University of Hawai'i. 18 Details of the HAHC study enrollment and clinical characterization were previously published 19 and approved by the University of Hawai'i Institutional Review Board #10675. All participants signed institutional review board-approved informed consent forms before participation. CD4+ T lymphocyte counts were obtained in real time by standard technique from a local CAP-certified reference laboratory.
Peripheral blood samples were collected in EDTA-coated collection tubes and samples were processed according to AIDS Clinical Trial Group SOPs. 20 Demographic and clinical details of study participants are provided in Supplementary Table S1. Overall patient demographics such as BMI, smoking status, and diabetes status were consistent among the three groups (Supplementary Table S1). Samples were grouped into three groups: HD (n = 16), PWH on cART and virally suppressed (cART; n = 20), and PWH off cART with viremia (Viremic; n = 14).
Cytomegalovirus antibody titer quantification
Matched cryopreserved plasma samples were used for the quantification of human cytomegalovirus (CMV) IgG and IgM antibody titers (Quest Diagnostics; order code 6732; Marlborough, MA) as previously published. 18
Multiplex assay determination of plasma cytokine/chemokine concentrations
As a microsphere bead capture assay, the ProcartaPlex multiplex assays require only 25 μL of plasma and takes only 4 h to yield analyzed results. 21,22 With Luminex multianalyte profiling (xMAP) technology, simultaneous analysis of up to 50 analytes per 20–50 μL sample can be accurately performed, 23 while an enzyme-linked immunosorbent assay (ELISA) generally requires a minimum sample volume of 100 μL. 24,25 As a precise method of plasma cytokine profiling, Luminex is more cost- and time-effective than ELISA and other immunoassays, making it ideally suitable for high-throughput biomarker identification, 26,27 especially in situations where sample availability is restricted.
For this assay, plasma samples were taken from −80°C storage and thawed at room temperature (RT) before being deactivated for viral contaminates by diluting 50/50 with 2% triton for 1 h at RT. For generation of standard curves, manufacturer-provided lyophilized standards were reconstituted and prepared by fourfold serial dilutions for a total of 8 standards, as described in the kit insert. These standards were analyzed in duplicate on a 96-well optical plate alongside the prepared plasma samples, which were also run in duplicate. Two additional wells were run without standards or samples to provide background measurements. The concentrations of plasma cytokines and chemokines were measured by Luminex xMAP technology, utilizing the Procartaplex Human Cytokine/Chemokine Convenience Panel 1A 34-plex (Thermo Fisher EPXR340-12167-901), in accordance with the manufacturer's instructions.
Briefly, 50 μL of vortexed magnetic capture beads was added to each well and placed on a handheld magnet to collect beads at the bottom of the well. The plate was then washed with provided wash buffer (WB) twice. Fifty microliters of samples was then added in duplicate to assigned wells.
For standards, 25 μL of each standard plus 25 μL of Universal Assay Buffer (UAB) was added to the appropriate wells. Duplicate background wells had 50 μL of UAB added. The plate was then sealed, placed on a shaker, and incubated at 4°C for 20 h. The plate was then washed twice, with WB and 25 μL of detection antibody added to each well, and shaken for 30 min at RT. The plate was then washed twice and 50 μL of Streptavidin-PE was added to each well and incubated at RT on the shaker for 30 min.
Finally, the plate was washed twice, 120 μL of Reading Buffer added to each well, and then the plate was incubated on shaker at RT for 5 min. Analysis of plate wells was performed using the Luminex 200 (Luminex Corporation), which was calibrated, and performance was validated as per instrument validation protocol. Measurements were reported with the xPONENT 4.2 software (Luminex Corporation). To evaluate the robustness of the assay, proportion of samples under the lower limit of detection (LLOD) is reported in Supplementary Table S3.
Bioinformatics analyses
Final concentrations of analytes were loaded into R v4.1.2. 28 Fold change values were calculated by first taking the mean expression of each analyte for each group and age group. Mean expression values were then divided by the matched age group of the HD cohort and the resulting fold changes were transformed with the log2 function. p Values were determined using two-tailed t-tests. Term enrichment analysis was performed with the pathfindR package. 29
Statistical tests
Statistical tests and violin plots were performed in GraphPad Prism v9.4.0. p Values were calculated using ordinary one-way analysis of variance followed by adjustment with Tukey's multiple comparison test. Spearman correlation was performed in GraphPad Prism for correlation analyses between analytes and clinical markers such as age, duration of known infection, CD4+ T cell counts, and viral loads (VLs).
Results
Cohort demographics
The cohort used for this study comprised three major age-mapped groups: HIV-uninfected donors (HD; n = 16; 32.6–73.5 years of age), PWH on cART (cART; n = 20; 26.7–60.2 years of age), and PWH off cART and with viremia (Viremic; n = 14; 27.5–59.4 years of age). Viremic patients were off cART and had viremia at the time of sample acquisition. The median age of each group was comparable (Supplementary Table S1). Other demographic information such as percent male, CMV IgG titers, VL measures, CD4 counts, and duration of known infection is available in Supplementary Table S1.
Analyte concentrations
Most analyte concentration measurements were not significantly different between groups, but clear observable trends were informative to provide expected reference ranges for each of the analyte measured (Figs. 1 − 4; Supplementary Table S2). For the purposes of data analysis, analytes were separated into groupings based on similarity of action in vivo 21 : CD4+ T helper (Th) lymphocyte type 1 and 2 (Th1/Th2) associated cytokines (Fig. 1), Th9/Th17/Th22 and regulatory T cell (Treg) associated cytokines (Fig. 2), inflammatory cytokines (Fig. 3), and chemokines (Fig. 4). Concentrations of interleukin-2 (IL-2), IL-6, IL-8, tumor necrosis factor-alpha (TNF-α) (Fig. 1), and IL-15 (Fig. 3) were consistently higher among individuals older than 45 years compared to those younger than 45 years for all groups, although they did not reach statistical significance.

Plasma concentrations of Th1/Th2 analytes. Violin plots of Th1/Th2 analyte concentrations among HD, cART, and viremic individuals younger than vs older than 45 years.

Plasma concentrations of Th9/Th17/Th22/Treg analytes. Violin plots of Th9/Th17/Th22/Treg analyte concentrations among HD, persons on cART, and viremic PWH younger than vs older than 45 years.

Plasma concentrations of inflammatory cytokines. Violin plots of inflammatory cytokine concentrations among HD, persons on cART, and viremic PWH younger than vs older than 45 years.

Plasma concentrations of chemokines. Violin plots of chemokines concentrations among HD, person on cART, and viremic PWH younger than vs older than 45 years.
Analyte correlations
HIV infection has been shown to have wide-ranging system changes throughout infection. 30 –32 To determine possible clinical correlations to the analytes in our study, we performed Spearman correlations between each analyte and CD4+ T cell counts, duration of known infection, and VLs in viremic patients (Table 1). Of note, we observed several analytes that had significant or near-significant correlations with the duration of known infection in the cART group, including CCL11 (R = 0.38, p = .0942), granulocyte macrophage colony-stimulating factor (GM-CSF) (R = −0.49, p = .0271), IL-1α (R = −0.49, p = .0286), and IL-4 (R = −0.44, p = .0542). These findings are not entirely unexpected due to previous reports often showing pathologically lower levels of some cytokines in treated PWH compared to untreated PWH. 33
Analyte Correlations Analysis
CCL11 also had a significant correlation with age among viremic patients (R = 0.75, p = .0030) (Table 1). Macrophage inflammatory protein-1 alpha (MIP-1α/CCL3) also had a significant correlation with age, in this case, among PWH on cART (R = −0.52, p = .0192) (Table 1). Surprisingly, few correlations in viremic PWH are significant or show moderate trends when examining VLs.
The chemokine C-X-C motif ligand 1 (GROα/CXCL1) had a significant correlation with HIV VL (R = 0.64, p = .0164), and was also significantly correlated with age within patients on cART (R = −0.63, p = .0030), as well as with CD4+ T cell counts within viremic patients (R = −0.75, p = .0031) (Table 1). Otherwise, IL-8 (R = −0.41, p = .0729) and IL-10 (R = −0.38, p = .0984) each had moderate correlations with CD4+ T cell counts among PWH on cART (Table 1). Taken together, our findings suggest that analyte correlations in PWH on cART were modulated differently compared to viremic PWH.
Analyte pathway analysis
Heatmap visualization of analyte concentrations (Fig. 5) highlights the variance by samples for all analytes. While the degree of variability was high, as shown above in correlative analyses, some associations do appear for analytes with higher expression in individuals older than 45 years for all groups. To get a better understanding of underlying expression patterns, analyte pathway analysis was subsequently performed (Supplementary Fig. S1). In PWH on cART, older than 45 years, there is an upregulation compared to age-matched HD in the “Viral protein interaction with cytokine and cytokine receptor,” and more specifically upregulations in the chemokines MCP-1/CCL2, CCL11, and RANTES/CCL5. These indicate a potential increase in immune trafficking, but whether the trafficking is to sites of HIV reservoir or other opportunistic infections related to age is currently unknown.

Luminex output. Heatmap plotted as individual samples showing analyte expressions per group: HD, PWH on cART, and viremic PWH. Each of these three groups was then divided into age groups, younger than 45 and older than 45 years. Colorimetric scale displays z-score values for that sample. Mean analyte expression was calculated per row and z-scores were calculated per sample. Color images are available online.
In addition, several Th1/Th2 related cytokines—interferon-gamma (IFN-γ), IL-4, IL-5, IL-12A, and IL-13—are all downregulated compared to HD. Similar analysis was performed for viremic PWH, although we did not observe similar results (data not shown), which may be a result of lower statistical power due to sample distributions. These findings suggest that PWH on cART exhibit increased immune trafficking, but lower inflammation, consistent with previous studies. 34,35
Principal component analysis (PCA) was performed using analyte concentrations without any input label for samples (sample ID, groups, and age). Label-free classification allows the samples to cluster based entirely on expression patterns instead of using classification labels as guides. After PCA reduction, group and age group labels were attached (Fig. 6). PCA analysis reveals that across age groups and HIV status, samples do not clearly cluster by any combination of age and HIV status (Fig. 6).

PCA output. PCA reduction using no group labels and all analytes.
However, there are two major clusters indicating potential for two distinct analyte profiles among subjects. PCA loadings comprised eigenvectors (direction) and eigenvalues (magnitude) of influence from each of the input variables. The direction and magnitude for each variable can be plotted to see their relative impact on the output PC dimensions. Upon examination of PCA loadings, this separation is driven more by analyte expression, high vs low, in the samples and the cluster with higher expression appears to favor samples from the cART group.
These results are consistent with prior studies, 35 which show that even though cART use reduces the chronic inflammation of HIV, it does not return it to levels seen in healthy controls. Notably, the analytes IL-2, TNF-α, IL-15, and several other interleukins have major impacts on the first principal component, while IL-1RA, RANTES/CCL5, and SDF-α are major driving factors for the second principal component. Based on these loadings, the first PC is largely driven by interleukin expression, while the second PC is largely driven by chemokine expression, highlighting the fact that inflammatory cytokine expression is still a major driver in PWH on cART separation from the other groups, regardless of age.
Discussion
Parallels between natural aging and the premature aging that correspond with HIV have been established, 36 and biomarkers with inflammation and immune activation may be associated with the accelerated aging observed during HIV infection. 36 In addition, HIV elite controllers have been observed to have higher levels of soluble inflammatory markers compared to HIV-suppressed or uninfected individuals, 37 and elevated levels of soluble biomarkers have been associated with occurrences of non-AIDS-defining events in PWH on cART. 38
There are many different types of cytokine profiles, each with critical roles in the immune response. 39 In this study, we examined profiles encapsulating Th1/Th2 helper T cells, Th9/Th17/Th22/Treg helper T cells, inflammatory cytokines, and chemokines. Th1/Th2 helper T cells have polarizing cytokine production profiles with an array of different immune responses. 40 The Th9/Th17/Th22/Treg class of helper T cells has activity that has been associated with immune balance and function in the context of HIV infection and its progression to AIDS. 41
Inflammatory cytokines have long been recognized for their vital role in protection against pathogens, although inappropriate activation can lead to acute and chronic disease, 21,42 and chronic inflammation has been associated with increased risk of disease in aging people. 43 Chemokines play vital roles as regulatory molecules in the function and trafficking of leukocytes, lymphocytes, and lymphoid tissues. 44 –46 HIV infection manipulates these functions to facilitate its own transmission and pathogenesis. 47 In combination, cytokines and chemokines regulate a healthy immune response throughout the body when functioning properly. 48,49
However, the dysregulation of cytokines and chemokines, which occurs during HIV infection, can increase pathogenesis by contributing to the depletion of CD4+ T cells and the elevation of VL. 50 Our focus in this study is to examine analyte expression profiles across three groups: HD, PWH on cART, and PWH with viremia. Subsequently we examined the impact of age for each of these groups by examining expression profiles comparing those younger than 45 years and older than 45 years. In addition, we also examined correlation of analyte expression with several clinical data points: age, CD4+ T cell count, duration of known infection, and VLs.
Of the Th1/Th2 analytes, IL-2 and IL-8 were both observed to have a significant negative correlation with CD4+ T cell counts among PWH on cART. IL-2 is a secretory cytokine that can stimulate T cells, B cells, and NK cells, and has in its recombinant form been utilized as a therapeutic for PWH concurrent with cART due to its ability to slow the apoptosis of HIV-infected cells. 51,52
Our findings with IL-2 indicate that PWH on cART may experience an increase in immune activation and proliferation, or the compensatory function of IL-2 as CD4+ T cell counts decrease. IL-8 is a proinflammatory cytokine, which has been identified as a marker of chronic inflammation within PHW on cART, 53 and blocking IL-8 in tandem with GROα/CXCL1 has been demonstrated to inhibit HIV replication within T cells and macrophages. 54 In addition, IL-8 and IFN-γ and each trended negatively with duration of known infection among viremic individuals.
IFN-γ is a cytokine secreted by activated T cells and NK cells with a wide range of immune response functions, such as macrophage activation, antiviral and antibacterial immunity mediation, and several other important immunoregulatory functions. 55 IFN-γ has been observed in higher concentrations within PWH compared to HD, 56 and PWH with high levels of circulating IFN-γ exhibited a lower CD4+ T cell count recovery after 4 years of being on cART. 57 In addition, IFN-γ trended positively with age among HD. This supports findings by past researchers that the production of IFN-γ is directly correlated with age. 43,58
The Th1/Th2 analytes GM-CSF and IL-4 each had a significant negative correlation with the duration of known infection among PWH on cART. GM-CSF promotes myeloid cell development and maturation, dendritic cell differentiation, and survival in vitro, and has a major role in the regulation of intestinal immune and inflammatory responses. 59 GM-CSF signaling in alveolar macrophages has been shown to be disrupted in HIV infection, 60 which may explain our observation of its decline over longer durations of infection. An imbalance in GM-CSF production or signaling may contribute to harmful inflammatory conditions, potentially giving it a role in the concept of inflammaging and/or chronic inflammation seen in long-term HIV infection. Meanwhile the cytokine IL-4 has been documented to regulate B and T cell growth, immunoglobulin secretion, eosinophil recruitment, and macrophage activation. 61,62
We observed a negative correlation with both duration of known infection and CD4 counts in PWH on cART, which supports a previous finding where serum levels of IL-4 were significantly lower in PWH on sustained highly active antiretroviral therapy (HAART) compared to before commencement of treatment. 63 Furthermore, IL-12p70 was found to be negatively trended with the duration of known infection among PWH on cART. IL-12p70 is the bioactive form of IL-12, comprising p35 and p40, and plays a role in the induction of IFN-γ production. 64,65 IL-12p70 is only transiently increased in viremic individuals compared with HD in correlation with VL, and its response to HIV infection may be overshadowed by other cytokines such as IL-23. 65
Among the Th9/Th17/Th22/Treg analytes, we observed several interesting trends. IL-9 trended negatively with the duration of known infection among viremic individuals and IL-10 trended negatively with CD4+ T cell counts among PWH on cART. Like IL-4, IL-10 is a proinflammatory cytokine that facilitates the differentiation of naive CD4+ T cells into Th2 subtypes. 66 IL-10 has been observed to significantly and gradually decrease within PWH on HAART.
IL-22 was found to be negatively trended with age among viremic individuals. This can be explained by the fact that HIV infection depletes mucosal Th22, Th17, and innate lymphoid cells, which are responsible for the production of inflammatory cytokines such as IL-22, although this deficiency can be rapidly restored by the administration of cART. 67 IL-23 was found to be negatively trended with the duration of known infection among PWH on cART. The inhibition of IL-23 signaling and expression has been reported within PHW on long-term HAART. 68,69
Regarding inflammatory cytokines, IL-1α had a significant negative correlation with the duration of known infection among PWH on cART. IL-1α induces proinflammatory effects and is central to the pathogenesis of numerous conditions characterized by body-wide organ or tissue inflammation. 70 Our findings in this study could be partially explained by the ability of cART to reduce systemic inflammation during acute HIV infection. 71 IL-31 trended negatively with the duration of known infection among viremic individuals. IL-31 is involved in maintaining hematopoietic stem cells and influencing proinflammatory activities of macrophages, monocytes, and T cells. 72,73 Serum concentrations of IL-31 have been reported to be higher in viremic patients compared to HD. 73
However, as IL-31 is primarily produced by CD4+ T cells, 74,75 and CD4+ T cells are depleted during HIV infection, 76,77 additional research is necessary to fully understand the potential relationship between IL-31 function and HIV infection. Tumor necrosis factor-beta (TNF-β) was found to be positively trended with VL among viremic individuals. TNF-β has been reported to stimulate HIV replication within T cells and monocyte-derived macrophages. 78
The last analyte grouping we analyzed was chemokines, among which CCL11 was found to be significantly positively correlated with age among viremic individuals, although it trended negatively with age among HD, indicating that it is impacted by some factor or combination of factors resulting from HIV infection. CCL11 plays a role in a skewed immune response toward a Th2 profile 79 ; Th2-type cytokines are associated with anti-inflammatory responses to mediate the proinflammatory responses of Th1-type cytokines. 80 In addition to its role in the immune response, CCL11 has been demonstrated to have associations with aging, neurogenesis and neurodegeneration, and microglia. 81 CCL11 also has a significant positive correlation with the duration of known infection among PWH on cART. High plasma levels of CCL11 have been related to a decrease in CD4+ T cells in HIV elite controllers with sustained viral suppression. 82
GROα/CXCL1 was found to be significantly negatively correlated with both age among PWH on cART, and with CD4+ T cell counts among viremic individuals. GROα/CXCL1 has a crucial role in inflammation, 83 and has been shown to have potential as a marker of disease activity. 84 GROα/CXCL1 is among a family of IL-8-related chemokines, which have been observed to be upregulated when exposed to HIV-infected T cells. 85 Perhaps the upregulation of GROα/CXCL1 initiated by HIV-infected T cells continues even as the T cells are depleted by HIV disease progression. 77 GROα/CXCL1 was also observed to have a significant positive correlation with VL among viremic individuals, further demonstrating its ability to be used as an indicator of disease activity. 84
Meanwhile, IP-10/CXCL10 and MCP-1/CCL2 were both found to be positively trended with VL among viremic individuals. Like GROα/CXCL1, IP-10/CXCL10 has also been identified as a robust soluble biomarker of monocyte activation with potential as a prognostic marker for HIV, as levels have been observed to increase along with VL. 86,87 MCP-1/CCL2 is a proinflammatory chemokine that supports the replication of HIV 88 and has been observed within aging PWH at significantly higher levels in comparison to age-stratified HD. 89
MIP-1α/CCL3 was found to be significantly negatively correlated with age among PWH on cART. MIP-1α/CCL3 is a macrophage-derived chemotactic chemokine with proinflammatory activities and a role in the effector immune response, which has been recognized as a potential marker for the detection of various diseases. 90,91 IP-10/CXCL10 trended negatively with CD4+ T cell counts among PWH on cART. IP-10/CXCL10 has been identified by other researchers at higher levels within PWH, who have confirmed lowered CD4+ T cell counts. 87 MIP-1β/CCL4 was found to be positively trended with the duration of known infection among viremic individuals.
MIP-1β/CCL4 is proinflammatory cytokine, which has been observed to have a significant positive correlation with plasma levels of the antiviral drug emtricitabine among individuals taking PrEP. 92 RANTES/CCL5 was found to be positively trended with the duration of known infection among PWH on cART. RANTES/CCL5 is expressed by inflammatory cells such as T cells and monocytes, and its interaction with the CCR5 receptor has been targeted as a strategy for long-term control of HIV infection. 93,94 These findings provide supportive evidence that a longer duration of known infection, and, indirectly, a longer duration of cART administration, corresponds with a reduction in inflammation markers.
Of all the analytes measured, concentrations of IL-2, IL-6, IL-8, IL-15, and TNF-α were higher among individuals older than 45 years compared to those younger than 45 years across groups. Raised levels of IL-6 have been associated with age-related disease, 95 –104 and IL-6 has also been considered a potential biomarker of clinical outcomes from HIV-related inflammation. 105 In addition, PWH receiving treatment have higher IL-6 levels than those who do not, and elevated IL-6 levels in PWH on cART are associated with advanced age and elevated body mass index. 101,105
The overactivation of TNF-α signaling has been associated with chronic inflammation and the eventual development of pathological complications. 106 TNF-α has been previously recorded within PWH taking long-term suppressive cART at levels comparable to those of elderly PWH who did not survive infection, and this may be due to persistent alteration of TNF-α levels in PWH causing tissue damage and influencing the development of non-AIDS-defining illnesses, while the replication of HIV itself is being controlled by cART. 107 IL-15 has established roles in the development and maintenance of NK and CD8+ T cells and has been observed to be significantly elevated within viremic individuals compared to HD. 108,109 Perhaps our observed increase in IL-15 could be the result of a natural attempt to overcome the reduced functionality of NK cells in aging adults. 110
There were some limitations that became evident during this study. Importantly, the majority of the subjects in this study were male. Considering that research has shown sex-based differences in several cytokine and immune profiles, 111,112 this discrepancy potentially influenced our results. While we preferred to have a balanced representation within our cohort, our samples were collected over 20 years before their use in this study and had since been depleted to the point that the male-to-female ratio had become skewed.
In addition, due to the cohort utilized, we did not have complete or extensive clinical information regarding cART regimens or duration, and so we are unable to determine whether the duration of cART impacted the levels of measured analytes. While we could access information about rough drug classes and cumulative duration of use, we did not have information regarding when and how long they were combined. Furthermore, without detailed treatment method information, we were unable to make comparisons of the older and younger members of each study group (HR, ART, and viremic). Our cohort was also lacking information regarding body mass index, clinical lipodystrophy, smoking, diabetes, or renal disease, all of which are factors that could affect cytokine/chemokine levels. 113 –122
Finally, the median time of infection differed between the cART and viremic PWH groups (5.39 and 11.1 years, respectively). Duration of infection may play a role in circulating cytokine levels, 33,50,123 although due to the original cohort aims, PWH not on cART were quickly enrolled in treatment regimens. As a result, the time of infection before treatment was much lower in the viremic PWH group. Overall, we made best use of what accompanying information we had available, although our findings would admittedly be more robust if there were substantive demographic and clinical data for our cohort.
Overall, we observed several significant correlations as well as nonsignificant trends between analytes and the conditions of age, CD4+ T cell counts, duration of known infection, and VL. As these data were acquired from plasma, it is more representative of systemic events, but can still be indicative that trafficking modulation plays a large role during HIV infection. Based on our findings, there may be moderate impacts on the Th1/Th2 cells, specifically inflammation from GM-CSF and cell survival from IL-2 and IL-4. While the sample size for this study was moderate, it may be underpowered to tease out subtle differences that may be masked by individual variability.
However, these data serve as an important reference for future studies seeking to examine changes in the immune landscape during aging with or without HIV. Our data further enlighten on the dynamics of the aging process with HIV, particularly immunological changes, which could help guide future strategies to address the complex clinical care that comes with this shifting age demographics of the PWH population.
Footnotes
Authors' Contributions
R.K.R. and L.C.N. designed the study. M.J.C., T.A.P., S.B., and C.M.S. facilitated sample evaluation and metadata. K.T. carried out Luminex assays. G.W. and K.W.K. analyzed data and prepared the article.
Data Availability
The data generated for this study are available upon request to the corresponding author.
Author Disclosure Statement
No competing financial interests exist.
Funding Information
This research was supported by National Institutes of Health (NIH) grants: R01AI120828, R01AI143457, and R01AI161010 (to R.K.R.).
Supplementary Material
Supplementary Figure S1
Supplementary Table S1
Supplementary Table S2
Supplementary Table S3
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
