Abstract
Cytokines may contribute to the severity of CD4 cell depletion with human immunodeficiency virus (HIV) infection, but quantitative relationships are not well defined. Serum and plasma from 181 HIV-infected individuals were tested with Millipore 30-plex Luminex cytokine assays. Within-individual correlations among cytokines were summarized by two-dimensional hierarchical cluster analysis. Associations with age, sex, race, CD4 count, and HIV viral load were determined with linear regression models. Tests for statistical significance were corrected for multiple comparisons, using a false discovery rate of 0.1. African-Americans had significantly higher levels than whites of six cytokines (IL-2, IL-5, IL-7, IL-15, fractalkine, and IFN-γ), and lower levels of MCP-1. Females had higher fractalkine levels than males. Age was not associated with levels of any cytokine. Six cytokines, including the T-helper (Th) type 1 cytokine IL-15, the Th2 cytokines IL-1ra and IL-10, the chemokines fractalkine and MCP-1, and the growth factor G-CSF were each inversely associated with CD4 count; no cytokine was directly associated with CD4 count. Fractalkine was directly associated with HIV viral load, adjusted for CD4 count. Cytokines clustered by primary function (e.g., Th1, Th2, proinflammatory, chemokines, or growth factors) whereas individuals clustered according to cytokine levels (generally high, intermediate, or low) had significantly different CD4 counts [medians (interquartile range) of 60 (17–162), 131 (62–321), and 155 (44–467), respectively; p<0.0001]. CD4 deficiency is associated with generalized increases in cytokines of various functions. Racial differences in cytokine response to HIV infection could contribute to disparities in disease progression.
Introduction
HIV
Materials and Methods
Subjects and samples
We analyzed stored serum (n=23) and plasma (n=158) samples obtained between 1983 and 2003 from 181 HIV-infected individuals. These individuals were originally selected as CD4-matched lymphoma-free controls for a study of AIDS non-Hodgkin lymphoma (NHL) in three National Cancer Institute HIV-infected cohorts. 9 Nineteen subjects were from a cohort of 133 HIV-infected male homosexuals from Manhattan and Washington, 10 103 were among 2126 HIV-infected hemophilia patients, 11 and 59 were in a study of 2803 patients at AIDS treatment and clinical trial sites. 12 The study participants were mainly white (83.4%), male (93.4%), and had a median age of 37.0 (IQR 31.4–42.9). Blood samples were transferred at ambient temperature to a central processing laboratory for separation of components within 24 h of phlebotomy. Serum and plasma were then immediately stored at −80°C until testing.
Cytokine assays
Cytokines were measured by Luminex fluorescent bead human cytokine immunoassays (MILLIPLEX MAP, Millipore Corp., Billerica, MA). Thirteen cytokines were part of a high-sensitivity multiplex panel [GM-CSF, IFN-γ, IL-1β, IL-2, IL-4, IL-5, IL-6, IL-7, IL-8, IL-10, IL-12 (p70), IL-13, and TNF-α], 16 were in a standard-sensitivity multiplex panel [EGF, eotaxin, fractalkine, G-CSF, IL-1α, IL-1ra, IL-12 (p40, free form), IL-15, IL-17, IP-10, MCP-1, MIP-1α, MIP-1β, soluble CD40L, TGF-α, and VEGF], and one was measured as a single analyte (RANTES). Measurements on blinded replicate samples included with test samples had coefficients of variation of 10–46%.
Statistical methods
All analyses were conducted using SAS 9.1 (SAS Institute, Cary, NC), BRB Array Tools 3.8.0 beta 3 (developed by the Biometric Research Branch of the U.S. National Cancer Institute
Using cytokine centering, Euclidean distance metrics, and complete linkage methods, we conducted hierarchical cluster analysis to construct classifications on two dimensions simultaneously (cytokines and subjects). 13 We tested whether the data exhibited more clustering than what is expected by chance from a multivariate normal distribution. As described in more detail by McShane et al., the multivariate normal distribution is unimodal and represents a single cluster. 14 To assess the association of individual cytokines with race, sex, age, and CD4 count, we employed class comparison analysis for categorical variables adjusted for material type (i.e., serum or plasma). 13 In this analysis, we identified cytokines that were significantly different among the classes by using a multivariate permutation test. 13 –15 The significance level was set to provide 80% confidence of no more than a 10% false discovery rate, defined as the proportion of cytokines incorrectly identified as significantly differing among the classes (i.e., the number of false positives divided by total positives). Although this technique assumed random variance t-statistics for each cytokine, the permutation test itself is nonparametric and does not require the assumption of Gaussian distribution. 16
We used a distance-based regression model to test whether the global pattern of cytokine measurements is influenced by CD4 count. The test for global association is based on a pseudo F statistic with its significance level evaluated through a permutation procedure with 10,000 iterations. 17
Antiretroviral therapy could be associated with altered cytokine levels. However, receipt of HAART could indicate severity of disease (i.e., treatment bias) in addition to direct effects on cytokine levels. We therefore examined variation of cytokine associations by year of blood draw as a proxy for treatment availability, rather than by actual treatment received. Accordingly, significant associations from the bivariate associations with CD4 count analysis, class comparison analysis, and quantitative trait analysis were adjusted for individual year of blood draw using linear regression models.
Results
Blood cytokine levels varied widely among the cytokines, with geometric means ranging from 1.0 (for IL-5) to 24,500 (for RANTES) pg/ml. However, interindividual distributions were more similar, with a geometric standard deviation between 2.3 and 5.3 (Table 1). CD4 count significantly declined with year of blood draw (β=−19, p<0001), with medians [interquartile range (IQR)] of 333 (226–459) for subjects enrolled during 1984–1990 (n=36), 65 (29–185) for 1991–1995 (n=70), and 72 (19–140) for 1996–2003 (n=75) (Supplementary Fig. S1; Supplementary Data are available online at
Excluding nondetectables.
The 30 cytokines were mutually correlated with one another, with the greatest pairwise correlations of 0.80 for sCD40L and RANTES, 0.78 for IFN-γ and interleukin (IL)-2, and 0.75 for IL-13 and IL-4. The cytokine with the lowest maximum correlation was G-CSF, based on its correlation of 0.30 with IL-17. All significant pairwise correlations were positive except for IL-2 with MCP-1, IL-2 with RANTES, IL-2 with Scd40l, and RANTES with IL12p70 (Table 2).
In the hierarchical cluster analysis, the global test of clustering was significant (p=0.001). We identified two main clusters of cytokines. The first cluster included four cytokines (GM-CSF, G-CSF, MIP-1b, and IL-4) that are known for their specific effects on growth, activation, differentiation, or chemotaxis of both neutrophils and B lymphocytes. Within the second cluster, cytokine levels were arranged in groups of similar function, including a predominantly Th1 group (IL-1β, IL-7, IL-12p70, IL-2, INF-γ, IL-15, and TGF-α), a Th2 group (IL-5 and IL-10), a proinflammatory group (IL-6, IL-8, RANTES, EOTAXIN, MCP-1, IP-10, and TNF-α), and a T cell chemokine group (fractalkine and IL-12p40) (Fig. 1, right-hand dendrogram).

Two-dimensional cluster analysis of 30 cytokine levels in 181 HIV-infected subjects. The right-hand dendrogram represents similarities among measured cytokines and the top dendrogram represents similarities among cytokine profiles of different subjects. Centered log-2 cytokine concentrations ranging from −17.2 to 16.6 are indicated in gradients of color, where red indicates higher values and green indicates lower values. Each row represents the values for a cytokine, and each column corresponds to a subject. Cluster 1 (green bar) indicates subjects with the highest cytokine concentrations, cluster 2 (orange bar) indicates subjects with the lowest concentrations, and cluster 3 (purple bar) indicates subjects with intermediate concentrations.
In the distance-based regression model, there was a significant global association between CD4 count and the measured cytokines (p<0.0001).
In the class comparison analysis adjusted for sample type, African-Americans had significantly higher (p range 0.0003 to 0.006) levels than whites of six cytokines (IL-2, IL-5, IL-7, IL-15, fractalkine, and IFN-γ) but significantly lower levels of MCP-1 (p=0.004) (Table 3). None of the cytokines significantly differed by age. When we analyzed our results by sex, only fractalkine was significantly higher in females compared to males [medians (IQR) 220 (114–296) versus 24 (0.1–231), respectively, p=0.004].
Cytokines were selected to provide 80% confidence that the false discovery rate was less than 10%.
Adjusted for sample type, CD4 count was inversely associated with significantly higher levels of six cytokines: the Th1 cytokine IL-15, the Th2 cytokines IL-1ra and IL-10, the chemokines fractalkine and MCP-1, and the growth factor G-CSF. The strongest association was for IL-1ra (9-fold difference between CD4< 200 and CD4 ≥200, p=0.0001) while the lowest association was for MCP-1 (1.5-fold difference, p=0.009) (Table 4). In the bivariate analysis, 17 of the remaining cytokines also varied inversely with CD4 count and only seven varied directly, but none of these 24 associations was significant with correction for multiple comparisons. Adjusted for blood draw date (as a surrogate for antiretroviral treatment availability), the associations of CD4 with G-CSF, IL-10, IL1-ra, and MCP-1 remained significant although attenuated, but the associations with fractalkine and IL-15 became nonsignificant. Results were qualitatively similar in a sensitivity analysis restricted to samples obtained in the pre-HAART era prior to 1996. Similarly, results did not change in a sensitivity analysis adjusted for draw year and study cohort.
Cytokines were selected to provide 80% confidence that the false discovery rate was less than 10%.
Of the 181 subjects, 179 fell into three clusters leaving two subjects that did not fit into any cluster (Fig. 1, top dendrogram). Compared to the second and third clusters, the first cluster (n=86) was characterized by higher levels for 28 of the 30 cytokines. The second cluster (n=15) had lower levels than either the first or the third cluster. Finally, the third cluster (n=78) was intermediate, with cytokine levels that were mostly lower than the first cluster and higher than the second cluster (Fig. 1). CD4 counts significantly differed among the three clusters of individuals, with median 60 (IQR, 17–162) in the first cluster with high levels of cytokines, median 155 (IQR, 44–467) in the second cluster with low levels, and median 131 (IQR, 62–321) in the third cluster with intermediate levels (p<0.0001).
HIV viral load did not affect the associations of cytokines with CD4 count, as only fractalkine was correlated (positively) with HIV viral load (r=0.2, p=0.007). Fractalkine remained significantly associated with both HIV viral load and CD4 count (β=0.50, p=0.05, and β=−0.01, p<0.0001, respectively) after mutual adjustment for both variables.
None of the subjects had cancer at the time of blood collection and only four developed cancer thereafter (latency intervals, 0.68 to 6.24 years). Exclusion of those four individuals did not substantially change the results.
Discussion
Cytokines, particularly certain chemokines, affect the risk of HIV acquisition and disease progression. 2 However, application of this biology to clinical practice has been hampered by their low levels in circulation, transient effects, and complex relationships. Our current study used new, high sensitivity, multiplex technology to quantify and characterize the relationships of circulating cytokine profiles to CD4 count and HIV viral load in well-characterized HIV-infected individuals.
We found an overall inverse relationship between the levels for 23 of 30 cytokines and CD4 count. These correlations were seen across cytokine functions, including significant associations with Th1, Th2, and proinflammatory cytokines and growth factors. Surprisingly, HIV viral load was associated only with fractalkine.
CD4 count declined significantly with the year of blood draw, reflecting the CD4 distribution of the matched lymphoma cases. Our findings of inverse associations of CD4 count with levels of both Th1 and Th2 cytokines are in contrast with some previous studies, and suggest that the theory of a Th1/Th2 switch is an oversimplification that does not explain the pathogenesis of HIV disease progression. 4 –6,18,19
African-Americans had higher levels of six cytokines (IL-5, fractalkine, IL-7, IL-2, IL-15, and IFN-γ) and lower levels of MCP-1. Associations with race may be due to genetic variation. Previous studies reported a race-specific distribution of allelic variants in cytokine genes. 20 –23 Alternatively, variations in cytokine levels by race could reflect differences in socioeconomic status (SES), which might indirectly affect disease severity. Issues associated with lower SES, such as being uninsured, limited access to healthcare, not receiving antiretroviral drugs, and lack of transportation to visit doctors could be related to HIV disease severity and thereby might affect cytokine levels. 24
Previous studies of people with HIV have not examined multiple circulating cytokines simultaneously. We evaluated a panel of 30 cytokines, finding many associations that are consistent with previous studies as well as some novel associations, such as lower CD4 with elevated IL-15. This range of findings underscores a strength of our study, namely our ability to begin untangling the cytokine network and its relationship to disease. In accord with our results, previous in vitro studies have demonstrated disease progression with excess production of proinflammatory cytokines such as TNF-α 25 and Th2 cytokines such as IL-10. 4
We found that fractalkine (CX3CL1), which is the only member of the CX3C chemokine family, was significantly associated not only with lower CD4 count but also with elevated HIV viral load. It had been suggested that fractalkine trafficking of HIV-infected lymphocytes propagates the dissemination of HIV in vivo. 26 Of the 30 cytokines that we evaluated, only fractalkine was associated with viral load, and this association remained significant after adjusting for CD4 count. Similarly, Widney et al. found that CXCL13—a B cell stimulatory chemokines—is elevated in HIV infection and correlated with CD4 count and HIV viral load. 27 In addition, CXCL13 levels were reported to be decreased with HAART. 27,28 However, some studies have found levels of beta-chemokines associated with CD4 changes, although not with HIV viral load, 29,30 and studies of other chemokines have yielded inconsistent results. 31–36
Our study has the strength of being the first to examine a comprehensive panel of circulating cytokines in HIV-infected individuals, enabling us to examine the individual and joint effects of multiple cytokines in vivo. Limitations of the study should be noted. Serum cytokine levels might be affected by antiretroviral therapy, although cytokine association with CD4 count did not vary by treatment availability, based on year of blood draw as a proxy variable. This adjustment suggests that cytokine associations with CD4 count are not strongly influenced by antiretroviral treatment, but a formal investigation into this assumption would be needed. Another limitation is the varying duration of sample storage prior to testing. However, Hosnijeh et al. previously reported that plasma and serum cytokine levels are stable and reproducible.37
In conclusion, we found HIV-related CD4 deficiency associated with generalized increases in Th1, Th2, and proinflammatory cytokines and growth factors. Racial differences in cytokine response to HIV infection could contribute to disparities in opportunistic infection and malignancy.
Footnotes
Acknowledgments
We would like to thank Dr. Eric Engels for his helpful suggestions on a previous version of this manuscript. This study was funded by the Intramural Research Program of the National Cancer Institute.
Author Disclosure Statement
No competing financial interests exist.
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.
