Abstract
The inflammatory context of HIV infection has been posited to contribute to the higher comorbidity risk noted in HIV-infected populations. One possible pathway may involve 1,25-dihydroxyvitamin D [1,25(OH)2D], which plays a wide biologic role in many tissues. We sought to investigate whether inflammation was associated with vitamin D metabolites in a cohort of HIV-infected (HIV+) men receiving treatment and HIV-uninfected (HIV−) men. Vitamin D metabolites, including 25-hydroxyvitamin D [25(OH)D] and 1,25(OH)2D, were measured along with 24 inflammatory markers among Multicenter AIDS Cohort Study participants. Exploratory factor analysis reduced inflammatory marker data to a smaller set of inflammatory processes (IPs). Multivariate linear regression was used to evaluate associations between vitamin D metabolites and IPs. There were 466 HIV+ and 100 HIV− men, who contributed 658 stored samples from 1998 to 2008. We found three IPs with IP 1 characterized by sTNF-R2, sIL-2Rα, sCD27, BAFF, sgp130, sCD14, CXCL10 (IP-10), and sIL-6R. While none of the three IPs was associated with 25(OH)D levels in either HIV+ or HIV−, higher levels of IP 1 were significantly associated with the reduced levels of 1,25(OH)2D in HIV+, and a similar although nonsignificant trend was seen in HIV−. The association between 1,25(OH)2D and inflammation found among HIV-infected men suggests a possible mechanism whereby inflammation leads to the increased comorbidity risk noted among HIV-infected individuals.
Introduction
V
Given evidence that both poor vitamin D status and residual inflammation may contribute to age-related morbidity risk, the question arises as to whether inflammation and vitamin D metabolism act through a single pathway to increase risk or whether they represent separate mechanisms. Vitamin D has two main forms as follows: 25-hydroxyvitamin D [25(OH)D], the most widely assessed but inactive form, and 1,25-dihydroxyvitamin D [1,25(OH)2D], the active form that regulates calcium-phosphate homeostasis. 1 In vitro studies have demonstrated immunomodulating effects of 1,25(OH)2D, as the active form acts on a variety of immune cells. 9 For instance, 1,25(OH)2D3 can promote monocyte-to-macrophage differentiation and downregulate the production of pro-inflammatory cytokines and chemokines. 9 In some contexts, inflammation also appears to interfere with activation and degradation of 1,25(OH)2D, 10 thus impeding the regulatory effects of vitamin D.
Existing data on the relationship between vitamin D and inflammation have come primarily from general elderly populations or from trials of vitamin D supplementation, but the nature of the relationship remains ambiguous. 11 HIV-infected populations represent a group with enduring inflammation and, thus, may be an ideal population for looking at this relationship. One recent prospective study among HIV treated individuals reported no change in inflammatory biomarkers following supplementation of vitamin D and restoration of 25(OH)D to normal levels. 12 Examining change in inflammatory biomarker levels and corresponding change in both inactive and active vitamin D metabolites within individuals would more definitively establish the potential links between inflammation and vitamin D metabolism.
In this study, we sought to evaluate associations of markers of inflammation and immune activation with vitamin D metabolites in HIV virologically suppressed men following treatment and HIV-uninfected counterparts, to elucidate relationships between general and HIV-specific patterns of inflammation and vitamin D metabolites.
Methods
Study population and design
A detailed description of the MACS has been provided elsewhere.
13
Briefly, it is a longstanding prospective cohort study that has followed 7,350 men who have sex with men from four metropolitan locations in the United States with a goal of understanding the clinical course of HIV-1 infection. Data on demographics, behavioral risks, clinical outcomes, and laboratory measurements are collected semiannually. Details on MACS can be found at:
The main analysis of present study was restricted to post-HAART samples of HAART users who had experienced viral suppression (defined as <400 copies/mL) following treatment and samples from comparable HIV-uninfected men. Secondary analysis included samples from a subset of above HAART users with both pre- and post-HAART measurements to evaluate relationships of change in inflammatory biomarkers and change in vitamin D metabolites around HAART initiation.
Vitamin D metabolites
Stored serum samples collected between 1999 and 2008 were tested for vitamin D metabolites, including 25(OH)D2, 25(OH)D3, 1,25(OH)2D2, and 1,25(OH)2D3. Testing was completed at University of Washington between October 2014 and February 2015, using immunoaffinity-liquid chromatography-tandem mass spectrometry with yielding limit of quantitation (LOQ) <1 ng/mL and coefficient of variation (CV) <5% for 25(OH)D and LOQ <4 pg/mL and CV <12% for 1,25(OH)2D. 14 Total levels of 25(OH)D and 1,25(OH)2D were obtained by summing their corresponding subforms. Consistent with the Institute of Medicine, 15 we defined vitamin D deficiency as 25(OH)D <20 ng/mL. To remove influences of race and season on vitamin D metabolites, 16 levels of 25(OH)D and 1,25(OH)2D were standardized by race and season by adding residuals of regressing vitamin D on race and season to intercepts.
Inflammatory biomarkers
A total of 24 serologic inflammatory biomarkers were quantified in the MACS inflammation substudy using two multiplex assays and a separate assay for C-reactive protein (CRP). 17 Measurement of 15 cytokines, including IL-1β, IL-2, IL-6, CXCL8 (IL-8), IL10, IL-12p70, GM-CSF, INF-γ, TNF-α, CXCL10 (IP-10), CCL11 (eotaxin), CCL2 (MCP-1), CCL13 (MCP-4), CCL4 (MIP-1β), and CCL17 (TARC), was completed using the pro-inflammatory 9-plex and Chemokine 7-plex (Meso-Scale Diagnostics, LLC, Rockville, MD). Another six soluble receptors (sCD14, sgp130, sIL-2Rα, sIL-6R, sTNF-R2, and sCD27), a chemokine CXCL13 (BLA-BCA1), and a cytokine (BAFF) were determined using the fluorescent bead-based multiplex Luminex xMAP system (Fluorokine MAP; R&D Systems, Minneapolis, MN) and analyzed using a Bio-Plex 200 Luminex instrument and Bio-Plex software (Bio-Rad, Hercules, CA). CRP was measured using a high-sensitivity nephelometric assay (Quest Diagnostics, Dade Behring, Inc., Newark, DE). To reduce variations, a single assay lot was chosen to implement testing, and samples for a single person were analyzed on the sample plate. Biomarker values were natural log transformed to achieve approximate normality distribution.
For the present study, 15 of 24 measured biomarkers were included in the main analysis. Excluded biomarkers were those with >10% of undetectable values [IL-1β (46%), GM-CSF (44%), IFN-γ (38%), IL-2 (23%), and IL-12p70 (11%)] and those for which there was no evidence from the literature or MACS data of associations with vitamin D (CXCL13, IL10, CCL17, and CRP).
Remaining missing values due to above or below the limits of detection were imputed using truncated log-normal distributions specific to the distribution of the biomarker.
Covariates
Covariates for inclusion in the analysis were drawn from prior work in the MACS evaluating correlates of vitamin D metabolite levels. 16 Covariate values were assessed at the time point of the vitamin D metabolite measurements. Enzyme-linked immunosorbent assays along with confirmatory Western blot determined HIV status (HIV infected vs. HIV uninfected). Race (white, black, and other race) and age (≥50 vs <50 years) were self-reported. Season was determined based on the timing of serum sample collection and categorized as spring (March, April, and May), summer (June, July, and August), fall (September, October, and November), and winter (December, January, and February). Body mass index (BMI) was grouped into three categories: normal (<25.0 kg/m 2 ), overweight (25.0–29.9 kg/m 2 ), and obese (≥30.0 kg/m 2 ). A few individuals with BMI <18.5 kg/m 2 were grouped into the normal BMI category. Hepatitis C virus (HCV) infection was considered present if a plasma RNA test was positive. Plasma HIV RNA was measured using the Roche assays (27% Roche second generation, LOQ <50 copies/mL; 68% Roche COBAS TaqMan, LOQ <20 copies/mL), and viral suppression was defined as <400 copies/mL. CD4+ T lymphocytes were measured using standardized flow cytometry. 18 The DHHS/Kaiser Panel guidelines 19 were used to define HAART use.
Statistical analyses
Characteristics of the study population were presented by HIV status using frequency and percent and compared using chi-square tests. The geometric mean difference in each biomarker and levels of 25(OH)D and 1,25(OH)2D by HIV status was statistically compared using nonparametric Wilcoxon tests.
We used exploratory factor analysis (EFA) to reduce the number of inflammatory exposure indicators and summarize their contribution to inflammatory processes (IPs). Implementation of EFA methods for inflammatory biomarkers in the MACS has been previously described. 20 Briefly, inflammatory biomarker levels were standardized [mean = 0, standard deviation (SD) = 1]. The strength and direction of correlations between biomarkers were then used to infer the presence of underlying latent processes or factors. A weighted linear combination of biomarkers was used to estimate individual levels of each latent process (factor scores). The derived weights (factor loadings) represent the strength and direction of the correlation of each biomarker with the underlying latent IP.
Next, standardized continuous 25(OH)D and 1,25(OH)2D levels were separately regressed on EFA-identified IPs. Generalized estimating equation was used to account for multiple visits contributed by some individuals. We selected covariates into the final model based on literature 21 and previous study results. 16 Final models were adjusted for age >50 years, BMI, and HCV infection and were run stratified by HIV status.
To evaluate individual biomarker contributions, relationships of individual biomarkers with vitamin D metabolite levels were also explored, accounting for multiple comparisons using the Benjamini–Hochberg method. 22 A sensitivity analysis was conducted, treating vitamin D metabolite levels as categorical, with 25(OH)D <20 ng/mL defined as vitamin D deficient and 1,25(OH)2D < 33 pg/mL (lowest 10th percentile) defined as low 1,25(OH)2D. Also a secondary analysis of change in selected individual biomarkers and change in vitamin D metabolites around HAART initiation was conducted. We considered p-value <.05 as indicative of statistical significance. All statistical analyses were completed using SAS 9.4 (SAS Institute, Cary, NC).
Results
Description of study population
A total of 566 men contributed 658 stored serum samples to the analysis: 466 post-HAART samples from HIV virologically suppressed men following treatment and 192 samples from HIV-uninfected men (Table 1). The age distribution was similar between the two groups. There were significant differences in the distribution of race and BMI by HIV status. A higher prevalence of HCV infection was seen among HIV-uninfected men as a result of oversampling in the MACS substudy of inflammation and immune activation. 7 Among HIV-infected men, the median time from HAART initiation to the biomarker measurement was 2.2 years [interquartile range (IQR): 1.7–2.3], and the median CD4+ T cell count at the time of the biomarker measurements was 531 cells/mm 3 (IQR: 367–717). Among the HIV-infected men, 205 had both pre- and post-HAART samples from which to evaluate the relationship between change in levels of inflammatory biomarkers and change in vitamin D metabolites.
Refers to post-HAART samples from HIV-infected men with suppressed viral load (≤400 copies/mL).
Refers to samples from HIV-uninfected men.
BMI, body mass index; HAART, highly active antiretroviral therapy; IQR, interquartile range; 25(OH)D, 25-hydroxyvitamin D; 1,25(OH)2D, 1,25-dihydroxyvitamin D.
Difference in vitamin D metabolites and biomarkers
The median standardized levels of 25(OH)D and 1,25(OH)2D did not statistically differ by HIV status (Table 1). The median change in levels of 25(OH)D and 1,25(OH)2D among 205 HIV virologically suppressed men with two measurements was 0.44 ng/mL (IQR: −5.75, 5.37) and 2.12 pg/mL (IQR: −8.98, 11.99), respectively. As reported previously, 7 despite experiencing HAART, more than a half of individual biomarkers measured were higher among HIV-infected men than HIV-uninfected men (Table 2).
SD, standard deviation.
Exploratory factor analysis
Among 658 samples included, 31 measurements from 6 biomarkers yielded values below or above the limits of detection and were imputed. The correlation between 15 biomarkers and factor loadings indicating the contribution of each biomarker to underlying IPs is consistent with previous reports in the MACS. 20 Three IPs were identified. IP 1 was primarily characterized by sTNF-R2, sCD27, sIL-2Rα, BAFF, sgp130, sCD14, CXCL10 (IP-10), and sIL-6R; IP 2 was primarily characterized by TNF-α, IL-6, IL-8, and MIP-1; and IP 3 was mainly characterized by MCP-1, eotaxin, and MCP-4.
Relationship between vitamin D metabolites and IPs
In the multivariate analysis, none of the IPs was associated with levels of 25(OH)D in either HIV group. However, levels of IP 1 were significantly associated with lower levels of 1,25(OH)2D among HIV-infected men, such that every SD increase in IP 1 was associated with a 1.82 pg/mL lower level of 1,25(OH)2D [95% confidence interval (CI): −3.44 to −0.21]. A similar magnitude of effect was seen among HIV-uninfected men, but the point estimate did not reach statistical significance in the smaller HIV-uninfected samples (β = −2.18 pg/mL, 95% CI: −4.86 to 0.51). IP 2 and 3 were not significantly associated with 1,25(OH)2D levels in either HIV group (Fig. 1).

Multivariate regression of inflammatory processes on levels of vitamin D metabolites by HIV status. Each inflammatory process was standardized with mean 0 and standard deviation 1.
Sensitivity analyses that categorized 25(OH)2D (<20 pg/mL) and 1,25(OH)2D (<33 pg/mL, the lowest 10th percentile) to evaluate the assumption of linearity of effects provided analogous results: none of IPs was significantly related to deficiency of 25(OH)D in either HIV group, while IP 1 levels were borderline significantly associated with the lowest 10th percentile of 1,25(OH)2D (odds ratio = 1.36, 95% CI: 0.96–1.96) among HIV-infected men.
Relationship between vitamin D metabolites and individual biomarkers from IPs
In the multivariate analysis, we found only a weak association of CXCL10 (IP-10) with levels of 25(OH)D among HIV-uninfected men and that of eotaxin with 25(OH)D levels among HIV-infected men. However, these associations were not statistically significant after adjusting the significance threshold to account for multiple comparisons.
In contrast, two markers, including sTNF-R2 and sCD27, that contributed to IP 1 were significantly associated with lower levels of 1,25(OH)2D levels among HIV-infected men. Among HIV-uninfected men, sTNF-R2 from IP 1 was associated with decreased levels of 1,25(OH)2D, while TNF-α from IP 2 was associated with lower levels of 1,25(OH)2D, and MIP-1 from IP 2 had an opposite effect. After accounting for multiple comparisons, the only significant relationships were with sTNF-R2 among HIV-infected men and TNF-α among HIV-uninfected men. There was no differential effect of biomarkers on either 25(OH)D or 1,25(OH)2D by HIV status (Table 3).
The concentration of each biomarker had been natural log-transformed and then standardized with mean = 0 and standard deviation = 1.
The model was adjusted for age, BMI, and HCV infection. An example of interpretation is that one SD increase of Eotaxin was significantly associated with 0.82 pg/mL higher levels of 25(OH)D in HIV+.
Bold font refers to P-value <0.05 from the multivariate regression model.
CI, confidence interval; HCV, hepatitis C virus.
Relationship between change in vitamin D metabolite and change in sTNF-R2 or sCD27
Neither change in sTNF-R2 or sCD27 was associated with change in levels of 25(OH)D. However, change in sCD27—a member of the tumor necrosis factor receptor superfamily that is expressed primarily on T cells—was significantly associated with change in 1,25(OH)2D (β = −2.98 pg/mL, 95% CI: −5.38 to −0.58) (Table 4).
Change in sTNF-R2 or sCD27 had been standardized with mean = 0 and standard deviation = 1.
The model was adjusted for age, BMI, and HCV infection. The interpretation is similar as above.
Analysis was restricted to 205 HIV+ who had pre- and post-HAART samples and available sTNF-R2 and sCD27.
Bold font refers to P-value <0.05 from the multivariate regression model.
Discussion
In our study, we found an association of markers important for immune responses toward infection and inflammation with levels of 1,25(OH)2D—the active vitamin D metabolite—among HIV-infected men following treatment, but not with levels of 25(OH)D. This relationship appeared to be driven by sTNF-R2 and sCD27, as reflected in the individual biomarker analysis and change analysis around HAART initiation.
In the MACS cohort, HIV infection did not appear to increase the risk of poor 25(OH)D and 1,25(OH)2D status associated with inflammation, at least among those with well-treated HIV. A similar point estimate for the association of IP 1 with lower levels of 1,25(OH)2D was observed among HIV-uninfected to that seen among HIV-infected men; both groups also had comparable standardized levels of 1,25(OH)2D. Thus the inflammatory relationship noted in our study does not appear to be specific to HIV infection, although HIV infection is associated with higher levels of many inflammatory markers.
Some evidence suggested that intracellular bacteria can disrupt vitamin D metabolism leading to increased extrarenal production of 1,25(OH)2D, consequently depleting 25(OH)D while disabling kidney's control of 1,25(OH)2D production. 23 Intracellular bacteria may block vitamin D receptors and cause an increase in inflammatory cytokines (e.g., IFN-γ). These cytokines are able to upregulate CYP27B1, a key enzyme converting 25(OH)D to 1,25(OH)2D in the extrarenal tissues. As our data show, the median levels of 25(OH)D, 18.5 ng/mL among HIV-infected men and 18.0 ng/mL among HIV-uninfected men, were relatively lower than the adjusted mean levels of 25(OH)D, 25.2 ng/mL, in the general men population of United States, 24 while the median levels of 1,25(OH)2D, 48 pg/mL, in our study samples were at the higher end of the reported normal range, 20–60 pg/mL. 25
However, our results were not fully consistent with such a bacterial infection model. Several antibacterial cytokines, including IL-6 and IL-8 (all indicators of IP 2), were not associated with 25(OH)D or 1,25(OH)2D. Of the important chemokine identified in the pathogenesis of bacterial infection, 26 only MIP-1α was positively associated with 1,25(OH)2D, but not with 25(OH)D. Likely, the relationship between cytokines and vitamin D metabolites is complex, involving multiple pathways. For example, TNF-α can either upregulate or downregulate the production of 1,25(OH)2D under different conditions. 27 Furthermore, 1,25(OH)2D as an immunomodulatory agent can inhibit the production of TNF-α. 28 Thus relationships noted in our analysis between an IP related to immune activation (i.e., IP 1) and 1,25(OH)2D but not 25(OH)D levels may represent an alternative mechanism disrupting production of 1,25(OH)2D. Alternatively, the mild-to-moderate inflammatory state of this well-treated and suppressed HIV-infected samples may not reach a threshold sufficient for impacting 25(OH)D levels, as seen in other highly inflammatory conditions (e.g., multiple sclerosis, acute infantile congestive heart failure, inflammatory bowel disease, and cystic fibrosis). 29
Examining individual inflammatory markers, it appeared that sTNF-R2 and sCD27 drove the association of IP 1 with levels of 1,25(OH)2D. Looking at changes in marker levels, the association between change in sCD27 and change in levels of 1,25(OH)2D was consistent with a relationship between immune markers and vitamin D metabolites. Both sCD27 and sTNF-R2 belong to members of the TNF-R superfamily, 30 and their elevation represents a state of immune activation and chronic inflammation. 31 The sTNF-R2, in certain conditions, may act to prolong the half-life of circulated TNF-α or to facilitate the interaction of TNF-α with receptors on the cell surface. 32
The strengths of present study include measurements of both vitamin D metabolites and a large number of inflammatory and immune activation biomarkers in a well-characterized HIV cohort. Furthermore, the inclusion of an HIV-uninfected group provides a comparison from which to assess differences due to the context of HIV infection and associated inflammation. However, we did not capture other biomarkers that are important for a full characterization of vitamin D metabolism, including serum calcium and phosphorus, PTH and FGF23, the latter of which play a critical role in regulation of renal production of 1,25(OH)2D. 33
Of note, the three IPs were characterized in a data-driven approach instead of a hypothesis-driven approach. The clustering of each IP purely depended on their intercorrelation. Thus identifying biological pathways for these identified processes may be difficult.
In conclusion, we found that an IP mainly driven by sTNF-R2 and sCD27 was associated with lower levels of the active form of vitamin D. These results indicate a potential for immune activation and inflammation to influence vitamin D metabolism and suggest a vitamin D endocrine pathway through which the inflammation could result in the excess risk of comorbidities. Longitudinal studies with repeated measures of inflammatory markers and vitamin D metabolites, as well as comorbidity ascertainment, could establish whether HIV-associated inflammation contributes to comorbidity risk through a pathway involving vitamin D metabolism.
Footnotes
Acknowledgments
The authors appreciate efforts from all collaborators, staffs, and participants of the Multicenter AIDS Cohort Study (MACS). Data in present study were collected by the MACS centers (Principal Investigators, grant number): institution, comprising Baltimore (J.B.M., U01-AI35042): Johns Hopkins Bloomberg School of Public Health; Chicago (Steven M. Wolinsky, U01-AI35039): Northwestern University; Los Angeles (Roger Detels and Otoniel Martinez-Maza, U01-AI35040): University of California, Los Angeles; and Pittsburgh (Charles R. Rinaldo and L.A.K., U01-AI35041): University of Pittsburgh. Data were managed by Data Coordinating Center at Johns Hopkins Bloomberg School of Public Health (L.P.J. and Gypsyamber D’ Souza, UM1-AI35043). The MACS is primarily funded by the National Institute of Allergy and Infectious Diseases, with additional co-funding from the National Cancer Institute, the National Institute on Drug Abuse, and the National Institute of Mental Health. Targeted supplemental funding for specific projects was also provided by the National Heart, Lung, and Blood Institute and the National Institute on Deafness and Communication Disorders. MACS data collection is also supported by UL1-TR001079 (JHU ICTR) from the National Center for Advancing Translational Sciences a component of the National Institutes of Health (NIH) and NIH Roadmap for Medical Research. The funding for this substudy was supported by the National Institute of Allergy and Infectious Diseases (R21-AI-109817). The contents of this article are solely the responsibility of the authors and do not necessarily represent the official views of the NIH. The MACS website is located at
Authors' Contributions
L.Z. and A.G.A. contributed to the analysis and composition of the article. A.G.A., L.P.J., T.T.B., and J.B.M. were responsible for the design and conduct of the study. A.N.H. carried out the measurement of serum samples. A.G.A., L.P.J., T.T.B, M.D.W., F.J.P., L.A.K., and A.T. provided critical intellectual comments on revisions of the article.
Author Disclosure Statement
No competing financial interests exist.
