HIV-Infection, Atherosclerosis and the Inflammatory Pathway: Candidate Gene Study in a Spanish HIV-Infected Population

Background Higher prevalence of atherosclerosis and higher cardiovascular risk is observed in HIV-infected individuals. The biological mechanisms underlying these processes are unclear. Several studies have implicated genetic variants in the inflammatory genes in cardiovascular disease and in HIV natural course infection. Methods & Findings In this study we have tested the possible association between genetic variants in several inflammatory genes and asymptomatic cardiovascular disease measured by carotid intima media thickness (cIMT) and atherosclerotic plaque presence as dependent variables in 213 HIV-infected individuals. A total of 101 genetic variants in 25 candidate genes have been genotyped. Results were analyzed using Plink and SPSS statistical packages. We have found several polymorphisms in the genes ALOX5 (rs2115819 p = 0.009), ALOX5AP (rs9578196 p = 0.007; rs4769873 p = 0.004 and rs9315051 p = 0.0004), CX3CL1 (rs4151117 p = 0.040 and rs614230 p = 0.015) and CCL5 (rs3817655 p = 0.018 and rs2107538 p = 0.018) associated with atherosclerotic plaque. cIMT mean has been associated with CRP (1130864 p = 0.0003 and rs1800947 p = 0.008), IL1RN (rs380092 p = 0.002) and ALOX5AP (rs3885907 p = 0.02) genetic variants. Conclusions In this study we have found modest associations between genetic variants in several inflammatory genes and atherosclerotic plaque or cIMT. Nevertheless, our study adds evidence to the association between inflammatory pathway genetic variants and the atherosclerotic disease in HIV-infected individuals.


Background
Cardiovascular disease (CVD), has been identified as a major cause of death in HIV-infected people [1][2][3]. HIV-infected individuals have accelerated atherogenesis, which is associated to a high risk of suffering a cardiovascular event such as coronary artery disease (CAD), peripheral vascular disease and stroke. The prevalence of these events in HIV-infected subjects is higher than in the general population and has an earlier onset [4][5][6]. The biological mechanisms underlying such risk among HIV-infected people are unclear [7].
Several studies implicate inflammation processes in CVD [7][8][9][10]. Non-HIV infected cohort studies have demonstrated that markers of inflammation are strongly predictive of CVD events and mortality [8]. Inflammatory markers are elevated in HIVinfected patients in comparison with non-infected individuals [7,11]. It has been hypothesized that this increased inflammation may be the explanation for the elevated cardiovascular risk in HIV-infected individuals [9].
Mechanisms of immune activation and inflammation have been proposed as the cause of earlier CAD in HIV [7]. HIV is thought to play a crucial role in the pathogenesis of atherosclerosis in HIV-infection [6,11]. In addition to the high prevalence of traditional CVD risk factors in HIV-infected individuals, several factors such as immunosupression, inflammation, HIV ability to induce foam cell transformation, cumulative exposure to antiretroviral drugs, and mitochondrial and metabolic dysfunctions have been hypothesized to be involved in HIV associated atherosclerosis [11][12][13].
The immunological and inflammatory pathways have many shared genes that may interact in the pathogenesis of atherosclerosis in HIV-infected individuals. The aim of this study was to assess the implication of genetic variants in relevant inflammatory genes in the atherosclerotic disease of HIV infected subjects.

Study population
We performed a cross-sectional study with 213 Spanish Caucasians HIV-infected individuals attended in Hospital Universitari MútuaTerrassa (Terrassa, Catalonia, Spain). This project was approved by the local ethics committee (Comité ético de Investigación Clínica del Hospital Universitario MútuaTerrassa -Approval number: EO/0915). All Participants gave written informed consent for genetic testing. At the time of enrolment, demographic, clinical and biochemical variables were collected from each patient by interviews and from medical notes. Simultaneously ultrasonographic measures (cIMT and atherosclerotic plaque presence) were performed.

Carotid artery Ultrasound
Carotid Intima Media Thickness (cIMT) is the most widely used surrogate marker of atherosclerosis. It relies on the fact that the  presence of atherosclerosis in one vascular bed will correlate with the atherosclerosis in another vascular bed. High-resolution Bmode ultrasound of the carotid arteries has been well validated in the non infected population as a surrogate marker of cardiovascular risk [24]. cIMT measurements and atherosclerotic plaque presence were obtained for each patient using a B-mode ultrasound recording with a 7 to 14-MZ transductor. For the imaging studies, patients were placed in the supine position with their head in the midline position and tilted slightly upwards, and the heart in systole. For each participant, a total of eight measures were obtained: left and right carotid primitive and bulb region in both proximal and distal walls. The cIMT value was defined as the mean of all values excluding those corresponding to atherosclerotic plaque thickness. According to Manheim cIMT consensus, atherosclerotic plaque presence is defined as a focal structure encroaching into the arterial lumen by at least 0.5 mm or 50% of the surrounding IMT value [25]. It was recorded as a bimodal variable (yes/no). All measures were performed by the same operator with an experimental intra-variability of 4.7%.

Genotyping
Genomic DNA was extracted from whole blood using the QIAamp DNA Blood Mini kit (Qiagen, Izasa, Barcelona, Spain) following the manufacturer's instructions. Genotyping was performed using Competitive Allele PCR methodology (KBioscience, Hoddesdon, United Kingdom). Genes from the inflammatory pathway were selected. A total of 101 genetic variants or single nucleotide polymorphisms (SNPs) in 25 candidate genes were tested (Table 1). TagSNPs were selected from HapMap for maximum coverage of the genes (Selection criteria: R 2 .0.8 and minor allele frequency (maf).0.05). We also included some relevant SNPs from other studies in HIV and non-HIV infected populations.

Statistical analysis and quality control
Statistical analyses were performed using the G*Power Calculator [26], SPSS 20 (SPSS, Chicago, IL) and Plink 1.07 [27] statistical packages. The study had an 80% power (considering a = 0.05 and b = 0.95, two-sided) to capture the effect of SNPs with maf$0.1 and Odds ratio (OR)$2. For rarer alleles (maf, 0.1), the study had a 75% power to capture genetic effects with OR$3. For continuous variables the study had a 95% power (considering a = 0.05 and b = 0.95, two-sided) to capture the effect of SNPs with maf$0.1.
Genotype call rate, Hardy-Weinberg equilibrium (HWE) and maf was assessed for each SNP using Plink 1.07. Those SNPs with call rates,95% or maf,0.01 and individuals with more than 5% missing genotypes were excluded from the analyses.
Chi square analyses or fisher exact tests (for those alleles with less than 5 counts in any group) were used to investigate the possible association of atherosclerotic plaque with genetic variants. cIMT variable was recorded as a continuous variable and it was normally distributed. Linear regression analyses were used to compare cIMT values with gene alleles.
Age, abdominal obesity, metabolic syndrome, hypertension, diabetes, total cholesterol, non-HDL cholesterol, LDL cholesterol and triglycerides were significantly associated with atherosclerotic plaque presence (p,0.05 in all cases). Age, hypertension, diabetes, dyslipidemia, total cholesterol, non-HDL cholesterol and LDL cholesterol were significantly associated with cIMT (p,0.05 in all cases).
Lipid values were highly correlated (correlation coefficient. 0.75 in all cases). Only dyslipidemia was included as a clinical adjusting variable in the regression models. Metabolic syndrome was collinear with the variables diabetes, dyslipidemia, abdominal obesity and hypertension. Metabolic syndrome was not included as a clinical variable in the regression model for atherosclerotic plaque presence.
For atherosclerotic plaque presence, age, abdominal obesity, hypertension, diabetes and dyslipidemia were included in logistic regression analyses as adjusting clinical variables for all the findings. For cIMT, age, hypertension, diabetes, and dyslipidemia were included in the regression analyses as adjusting clinical variables for all the findings. For the importance of HAART in the vascular pathogenesis, it was included in both adjustments.

Study Population
Demographic, lipid profile and HIV baseline characteristics of the study population are shown in Table 2. Mean age was 45.34 (SD68.20) years. Males constituted 77.9% of the samples. 82.2% of the individuals were receiving Highly Active Antiretoviral Therapy (HAART) at the time of the study. Atherosclerotic plaque presence (indicator of atherosclerotic lesion [28]) and cIMT means (indicator of arterial remodelling [29]) were the dependent variables investigated as markers of CVD. The study population was classified in two groups according to the bimodal variable atherosclerotic plaque presence. Baseline characteristics are shown in Table 1. Patients with atherosclerotic plaque presence conformed 39% of the cohort (n = 83), with 73.5% males and a mean age of 49.67 (SD68.38) years. Patients without atherosclerotic plaque had a mean age of 42.58 (SD66.79) years and 80.8% were males.
SNP Quality control. Three of the genotyped SNPs were excluded from the analysis, because of a call rate,95% (IL1RN rs3087263, CX3CR1 rs2669845 and IL6 rs2069833) and two because of a maf,0.01 (CXCL10 rs11548618 and TNF rs4248160). Two SNPs (IL1B rs55778004 and IL1RN rs4252019) were monomorphic in our population. The remaining 94 SNPs were included in the analysis.
Single marker analyses by logistic regression using age, abdominal obesity, hypertension, diabetes, dyslipidemia and HAART as adjusting clinical variables showed 8 positive associations.

Discussion
In this study we have tested the possible association between genetic variants in several inflammatory genes and CVD measured by cIMT and atherosclerotic plaque presence in 213 HIV-infected individuals. We have found several polymorphisms in the ALOX5, ALOX5AP, CX3CL1 and CCL5 genes significantly associated with atherosclerotic plaque, whereas cIMT mean has been associated with CRP, IL1RN and ALOX5AP genetic variants.

Atherosclerotic Plaque
The ALOX5 rs2115819-C allele was associated with atherosclerotic plaque presence and the ALOX5AP rs9578196-T, rs4769873-T and rs9315051-G alleles with atherosclerotic plaque absence. Genetic variants in ALOX5 and ALOX5AP have been previously linked to CVD in non-HIV infected populations [21]. These two genes encode for two important proteins of the 5-LO pathway which has been previously linked to atherosclerosis [21]. Knock-out mice of ALOX5 showed less atherosclerotic plaque formation [21]. Studies in HIV-infected cultured cells have shown a diminished function of 5-LO pathway due to gp120 HIV-protein presence [30,31]. Our results, showing that ALOX5 and Table 4. Summary of single marker analyses in relation to cIMT. ALOX5AP may influence CVD in HIV patients, contribute to the growing evidence on the relevance of these genes on CVD. The CX3CL1 rs614230-C allele was found associated with absence of atherosclerotic plaque. This allele has been previously associated with smaller cIMT in a German non-HIV infected subjects (CAPS cohort) although this finding was not replicated in a French cohort (3C cohort) [18]. To our knowledge, the association found between CX3CL1 rs4151117-G allele and atherosclerotic plaque absence is a novel finding. CX3CL1 has been implicated in atherosclerotic plaque formation initial process [18,19]. CX3CR1 (CX3CL1 receptor) knockout mice and animal treated with CX3CR1 inhibitors show reduced atherosclerotic plaque formation [32,33]. The functional effects of the CX3CL1 rs614230 or rs4151117 are not known. However, these SNPs are located near a frame shift mutation that may diminish protein functionality. Our results seem to agree with the reported atheroprotective effects of CX3CL1 reduced activity.
The CCL5 rs3817655-A and rs2107538-T alleles were found associated with atherosclerotic plaque presence. The CCL5 rs2107538 polymorphism has been previously associated with higher plasma concentrations of CCL5 and increased risk of MI in Korean CAD patients [34] and Han Chinese MI patients [15]. Our results seem to agree with these findings as the T allele associated with plaque presence in our study is also associated with MI risk. The CCL5 rs3817655 finding is a novel association that has not been previously described. This SNP is near to two mutations that encode for new stop codons and truncated proteins. Our findings may reflect the functional effect of these mutations but further studies are required to confirm it.

Carotid Intima Media Thickness
The CRP rs1130864-T allele was found marginally associated with smaller cIMT and the rs1800947-C allele with greater cIMT. These results are in the opposite direction that in previous studies. The CRP rs1130864-T allele was found associated with increased CRP plasma levels [14], which are considered a marker of cardiovascular risk [35,36], and have been linked to a faster cIMT progression in HIV-infected individuals [37]. The CRP rs1800947-C allele have been previously linked to decreased CRP plasma levels [15,38], which are known to be atheroprotective. Although the functional consequences of these SNPs are not known, both of them are in regulatory regions. The presence of HIV could interact with the functionality of these genetic variants explaining the discrepancy between our findings and those published previously. These findings need to be replicated in independent cohorts to confirm the direction of the associations.
The IL1RN rs380092-T allele was found associated with smaller cIMT, although in previous studies has been associated with increased risk of MI [16,17,39]. The IL1RN rs380092 variant is located approximately 200 bp up-stream of another variant, rs121913161, which encodes for a new stop codon that significantly reduces IL1RN plasma levels. IL1RN has been suggested to be a proatherogenic cytokine [40] and diminished levels of IL1RN might be atheroprotective [41]. This would agree with our findings. However, this result needs to be repeated as this SNP was not found in HWE.
In studies performed by Shrestha et al. on American HIVinfected patients, an association was found between the RYR3 rs2229116 polymorphism and cIMT [22]. However, we were not able to replicate this finding. Aside from a false negative result, the discrepancy could be explained by the methodological differences in cIMT measures and/or different allele frequencies between American and Spanish populations.
This study has several limitations. Our cohort has a limited sample size, although it is very well characterized. High prevalence of cardiovascular risk traditional factors and the bias observed in HIV-infected populations should also be considered as a limitation and a possible source of error in our findings. The lack of an HIVnegative control group prevented drawing any conclusions about the contribution of HIV infection to early atherosclerosis. However, studies have suggested that the biological processes involved in HIV-related atherosclerosis are different from the ones observed in non-infected individuals. Moreover, it seems that the atherosclerotic plaque structure and characteristics are also different [42]. Finally, carotid ultrasonography measures are conditioned by its inter and intra-variability.
In summary, in this study we have found modest associations between genetic variants in several inflammatory genes and cardiovascular risk markers in an HIV-infected population. They need to be confirmed in a larger population and functional studies are needed to elucidate the consequences of these genetic changes. Nevertheless, our study adds evidence to the association between inflammatory pathway genetic variants and the atherosclerotic disease in HIV-infected individuals.