The advent and widespread utilization of combination antiretroviral therapy (cART) has dramatically reduced the incidence of opportunistic infectious diseases and AIDS-associated malignancies in PLWH, leading to a subsequent rise in the prevalence of serious non-AIDS complications as the primary driver of morbidity and mortality in this population [1,2]. Among PLWH, elevated levels of biomarkers associated with inflammation and coagulation including interleukin-6 (IL-6), D-dimer, and high-sensitivity C-reactive protein (hsCRP) have all been associated with an increased risk of cardiovascular disease (CVD), cancer, and all-cause mortality [3,4]. These biomarkers have also been associated with a higher risk of similar diseases and mortality in the general population [5–8]; an interesting discovery that has stimulated a field of inquiry into understanding the role of chronic inflammation as the primary driver of serious non-AIDS-related complications in PLWH [9–11].
Although the underlying mechanisms by which these biomarkers are involved in the pathogenesis of non-AIDS complications remains an active area of research, studies have shown that elevations in these biomarkers decline with cART . However, despite plasma HIV RNA levels less than 50 copies/ml over long periods of time, many patients continue to exhibit persistently elevated levels relative to the general population, suggesting that viral replication alone is not the sole contributor to chronic inflammation in this population . Prior studies have shown that in certain ethnic populations, elevations of these biomarkers may have a partial genetic basis [14,15], although such genetic associations have not been as clearly established or examined in the HIV-infected population.
Genome-wide association studies (GWAS) have been utilized in the HIV-infected population to assess host determinants of a variety of phenotypes including host susceptibility to infection, innate virologic control , and cART pharmacogenomics [17,18]. However, to our knowledge, no similar investigation has previously been conducted to investigate the potential genetic basis for inflammation in this population, although similar genomic/GWAS analyses have been attempted for a host of other diseases [19–21]. Furthermore, GWAS analyses have been used to identify loci associated with elevated levels of hsCRP, D-dimer, and IL-6 in the general population [22–26]. The CRP, HNF1A, and APOE loci have been commonly identified to be associated with levels of hsCRP [22,23,26] and blood CRP concentration is associated with colorectal cancer . HNF1A is a transcription factor, which can bind to the CRP promoter  and APOE can attenuate unsolvable inflammation by complexing with activated C1q . Interestingly, hsCRP has been used to diagnose maturity-onset diabetes of the young (MODY) as a result of HNF1A . In another study, three loci, F3, F5, and FGA, were identified to be associated with D-dimer levels in healthy adults . Consistently, F3, F5, and FGA belong to the Complement and Coagulation Cascades pathway. In this study, we aimed to both qualitatively and quantitatively investigate host genetic contributions to elevation of biomarkers noted to be strongly associated with serious non-AIDS-related complications using GWAS analyses of data combined from three large international HIV studies: Evaluation of Subcutaneous Proleukin in a Randomized International Trial (ESPRIT) study , the Strategies for Management of Antiretroviral Therapy (SMART)  study, and Strategic Timing of Antiretroviral Therapy (START) study .
Materials and methods
The participants enrolled in this GWAS were enrolled in one of three international HIV treatment trials: ESPRIT , SMART , and START . Written informed consent was obtained from every participant in each study. All informed consents were reviewed and approved by participant site ethics review committees.
Genotyping and quality control
The study and analytical design are illustrated in Fig. 1. A total of 7768 participants were successfully genotyped using a custom Affymetrix Axiom SNP array at Advanced BioMedical Laboratories. The array consisted of 770 558 probe sets enriched with markers related to immune functions.
Genotypes were called using Axiom Analysis Suite (version 126.96.36.199, Thermo Fisher Scientific). Quality controls were performed using PLINK (version 1.9) . Individuals with any of the following were excluded: sex mismatch, autosome SNP call rate less than 96%, duplicates, estimated by pairwise identity-by-descent (IBD) (pi-hat at least 0.90). These exclusions resulted in 7720 participants. For each study, individuals with any of the following were further excluded: an autosome heterozygosity rate outside three standard deviations of the mean of individuals in each assigned ancestral group: 93 participants were excluded; cryptic relatedness individuals (pi-hat >0.1875) with a lower quality: 435 participates were excluded (see Fig. 1).
ADMIXTURE (version 1.3.0)  was used to assign ancestry for each individual and principal component analysis (PCA) was applied to account for population stratification. The details of SNP quality control, ancestry assignment, and PCA analysis are described in detail within Supplemental Material and Methods, http://links.lww.com/QAD/B883.
Genome-wide association studies
Single variant analysis was performed for participants in each ancestral group of each study using multiple linear regression with an additive model (1 df) in PLINK, with the first 10 principal components, and the following other covariates measured at study entry: age, sex, natural log transformed CD4+ cell count, natural log transformed viral load, BMI, smoker at baseline (missing in ESPRIT), hepatitis B virus (HBV), hepatitis C virus (HCV), diabetes, and history of CVD. Most covariates were used in all biomarker association analyses except the following: HCV was only used in IL-6 and D-dimer analyses, history of CVD was only used in hsCRP and IL-6 analyses, and history of diabetes was only used in D-dimer analysis. These covariate selections were based on bidirectional step-wise variable selection using AIC on potential variables in R , previous literature [35–37] and biological plausibility. Since not all the potential covariates of interest were collected in each trial, models were fit on each trial separately and covariates were included if any of the three stepwise regression models selected the potential covariate of interest. The variance inflation factor was also investigated to determine collinearity among the covariates chosen via stepwise regression and collinearity did not appear to be a significant problem. The levels of hsCRP, IL6, and D-dimer were natural log transformed.
To increase statistical power, results from the linear regression analyses across the three studies were combined in a fixed-effects meta-analysis using METAL  with inverse variance weighting for each ancestry group: European (EUR), African (AFR), and Admixed American (AMR). The EUR, AFR, and AMR meta-analyses results were further combined in a second transethnic fixed-effects meta-analysis to identify the SNPs associated with hsCRP, D-dimer, and IL-6 levels. Genomic control adjustment was turned on in METAL to adjust for any inflation. Lambda gc (λgc)  and quantile--quantile (Q--Q) plots were used to assess the system bias in the GWAS results. Manhattan plots were generated to visualize GWAS results. Q--Q plots and Manhattan plots were generated with the CMPlot R package  (version 3.3.4). A genome-wide significance (GWS) threshold was set at P less than 5 × 10−8. A locus was defined by significant SNPs within 100 kb of each other and annotated by AVIA  (v3.0, http://avia-abcc.ncifcrf.gov/) with the locus name based on the closest gene to the SNP with the most significant P value (lead SNP). Locus SNPs were subjected to stepwise conditional analysis using GCTA-COJO to determine if there were any secondary association signals (SAS) .
Study participants and quality control
After exclusion of sex mismatches and duplicates, 7720 out of 7768 genotyped participants remained from ESPRIT, SMART, and START (third row in Fig. 1). Baseline demographic, clinical, and laboratory characteristics including biomarkers of interest (D-dimer, hsCRP, and IL-6) of the remaining 7720 participants were acquired from previous studies [12,35,36] and were summarized (Supplemental Table 1, http://links.lww.com/QAD/B883). As a limitation of the data, hsCRP, D-dimer, and IL-6 measurements were not available for all 7720 participants listed in this table. The distribution of the major biomarkers in these three datasets are shown in Fig. 2. As the biomarker levels from the participants in the START study were measured before ART was initiated, median and IQR viral load values are higher than the other two studies and the viral load and CD4+ cell count distributions in START are different from those of the SMART and ESPRIT studies. Viral load and CD4+ cell count were used as covariates in the GWAS analysis to account for these differences between the studies. The meta-analysis process (see below) used in this study will further minimize the effect of the different treatment strategies between studies on biomarker levels. A total of 528 participants were excluded from further analysis after performing quality control analysis for cryptic relatedness and outlying heterozygosity (Fig. 1).
Quality control and imputation of SNPs was described in Supplemental Materials and Methods and illustrated in Supplemental Figure 1, http://links.lww.com/QAD/B883).
Meta-analysis of three multiethnic studies
Participants were partitioned into the five ancestral groups based on individual ancestry assignment. AFR, AMR, and EUR groups remained for association analysis whereas South and East Asian groups were excluded because of low sample size: 140 participants (Fig. 1). The GWAS analysis of each ethnic group from the three studies were combined by meta-analysis to identify the SNPs associated with hsCRP, D-dimer, and IL-6 levels.
To assess any systematic biases that may be present in our GWAS results, Q--Q plots were generated and λgc calculated (Fig. 3). These plots show that the ethnic-specific and transethnic analysis have normal distributions of P values. All related λgc are close to 1.0, demonstrating there was no genomic inflation in the GWAS analysis results.
Single nucleotide polymorphisms associated with levels of high-sensitivity C-reactive protein
Three gene loci associated with levels of hsCRP were identified for the transethnic group: CRP, HNF1 Homeobox A (HNF1A), and apolipoprotein E (APOE) (Fig. 4a, Supplemental Data 1, http://links.lww.com/QAD/B884). There was no evidence of heterogeneity across ethnicities as determined by the Cochran's Q-test for heterogeneity (HetPVal) for any of the lead SNPs at these loci (Table 1). Due to the small sample size, not all loci could be identified in all ethnic groups. However, there were SNPs in the CRP and HNF1A loci, which reached GWS in the EUR group only and exhibited heterogeneity across ethnic groups (Supplemental Data 1, http://links.lww.com/QAD/B884, Supplemental Figure 2, http://links.lww.com/QAD/B883).
Table 1 -
Lead variant of loci or second associated signal identified to be associated with the level of high-sensitivity C-reactive protein, D-dimer, and interleukin-6
||Number of SNPs
||Build 37 position (bp)
||159 677 654
||159 678 816
||121 424 406
||45 411 941
||95 052 476
||169 519 049
||155 488 821
||79 156 992
||154 395 212
Effect allele is minor allele. P, β, EAF, and HetPval are P value, effect size, effect allele frequency, and Cochran's Q-test heterogeneity P value, respectively. There is no evidence of heterogeneity across ethnicities for the lead variants. Covariates included in the model: first 10 principal components, and the following measured at study entry: age, sex, natural log transformed CD4+ cell count, natural log transformed viral load, BMI, smoker baseline (missing in ESPRIT), hepatitis B virus (HBV), hepatitis C virus (HCV), diabetes and history of CVD. HCV was only used in IL-6 and D-dimer analyses, history of CVD was only used in hsCRP and IL-6 analyses, and history of diabetes was only used in D-dimer analysis. EAF, effect allele frequencies; hsCRP, high-sensitivity C-reactive protein; IL-6, interleukin-6.
The 77 identified SNPs in the CRP locus can be divided into two subsets with different effect sizes and ancestral frequencies (Supplemental Data 1, spreadsheet ‘CRP’, http://links.lww.com/QAD/B884). The most significant SNP in the first subset, rs6667499 (P value: P = 2.38 × 10−25, effect size: β = 0.569), is located downstream of the CRP gene. Effect allele frequencies (EAF) for rs6667499 in transethnic, AFR, AMR, and EUR groups were 4.4, 16.9, 1.2, and 0.3%, respectively. Results for this SNP from AMR and EUR groups were not available because of low EAF (Table 1). The most significant SNP in the second subset was rs2794520 (P = 5.43 × 10−24, β = −0.224), also located downstream of the CRP gene. EAF for rs2794520 in transethnic, AFR, AMR, and EUR groups were 30.9, 20, 35.6, and 34.3%, respectively (Table 1). As a result of the differences in effects and ethnic EAF, we treated the first subset of 32 SNPs as the first SAS of the CRP locus with rs6667499 as the lead SNP and the second subset of 45 SNPs as the second SAS with rs2794520 as the lead SNP. Within the second SAS, the variants can be further divided into two types: positive and negative β. The relationship between these variants were assessed by LDpair in LDlink  (Supplemental Figure 5A, http://links.lww.com/QAD/B883). In most cases (659 out of 664), the minor allele of rs3116635(G) (lead positive β group) commigrates with the major allele of rs2794520(C) (lead negative β group), explaining the opposite effect size associated with the minor alleles of the two variants (Supplemental Figure 5A, http://links.lww.com/QAD/B883). An additional variant, rs12093699, was identified as a EUR group-specific variant associated with hsCRP levels with P = 2.01 × 10−08 and β = 0.148 (Supplemental Data 1, Spreadsheet ‘CRP’, http://links.lww.com/QAD/B884). This variant is located downstream of the CRP gene and showed evidence of heterogeneity across ethnic groups (HetPval = 2.52 × 10−06) with AFR and AMR effect sizes equal to −0.118 and 0.037, respectively.
HNF1A intron variant, rs2393776 (P = 9.90 × 10−10, β = −0.129), is the lead variant of 34 SNPs at the HNF1A locus (Supplemental Data 1 Spreadsheet ‘HNF1A’, http://links.lww.com/QAD/B884). EAF for transethnic, AFR, AMR, and EUR groups were 38.6, 29.6, 45.4, and 40.9%, respectively (Table 1). Like the second SAS in the CRP locus, the HNF1A locus is made up of two sets of SNPs with opposite effects and in most cases (606 out of 619), the minor allele of rs11065394(G) (lead positive β) commigrates with the major allele of rs2393776(A) (lead negative β), explaining the difference in effect size associated with the minor alleles of the two variants (Supplemental Figure 5B, http://links.lww.com/QAD/B883). Six additional variants with a positive effect were identified within the HNF1A locus, reaching GWS in the EUR but not the transethnic analysis. The lead variant for this group, rs2264779 (P = 1.10 × 10−09, β = −0.155), exhibited heterogeneity between the ethnicities (HetPval = 1.41 × 10−3) with AFR and AMR effect sizes equal to −0.014 and −0.010, respectively (Supplemental Data 1 Spreadsheet ‘HNF1A’, http://links.lww.com/QAD/B884).
An APOE missense variant, rs429358 (P = 1.57 × 10−22, β = −0.282), was the lead variant in the APOE locus of 11 identified SNPs (Supplemental Data 1, Spreadsheet ‘APOE’, http://links.lww.com/QAD/B884). No evidence of heterogeneity was found across ethnicities. EAF for transethnic, AFR, AMR, and EUR datasets were 15.4, 22.9, 12.1, and 13.1%, respectively (Table 1).
Single nucleotide polymorphisms associated with levels of D-dimer
Twenty-one variants associated with D-dimer levels were identified at four loci for the transethnic group: Coagulation Factor III (F3), Coagulation factor V (F5), Fibrinogen Beta Chain (FGB), and Glucosaminyl (N-Acetyl) Transferase 1 (GCNT1) (Fig. 4b, Supplemental Data 2, http://links.lww.com/QAD/B885). There was no evidence of heterogeneity across ethnicities for any of the lead SNPs at these loci (Table 1). Due to the small sample size, not all loci could be identified in all ethnic groups (Supplemental Data 2, http://links.lww.com/QAD/B885, Supplemental Figure 3, http://links.lww.com/QAD/B883).
The variant, rs2022309 (P = 1.86 × 10−13, β = 0.092), upstream of the F3 gene, was the lead variant for the F3 locus of 17 SNPs. EAF for transethnic, AFR, AMR, and EUR groups were 23.3, 8.9, 26.8, and 28.1%, respectively (Table 1).
The only variant within the F5 loci, rs6025 (P = 3.57 × 10−11, β = 0.257), is a nonsynonymous SNP in the F5 gene (Supplemental Data 2, http://links.lww.com/QAD/B885). EAF for the transethnic and EUR groups were 1.8 and 2.5%, respectively (Table 1). EAF for the AFR and AMR datasets were less than 1%. Therefore, they were not analyzed further.
A synonymous variant within the FGB gene, rs6056 (P = 5.61 × 10−09, β = −0.082), is the lead of two SNPs for the FGB locus (Supplemental Data 2, http://links.lww.com/QAD/B885). EAF of this variant in the transethnic AFR, AMR, and EUR groups were 16.3, 8.9, 15.4, and 19.2%, respectively (Table 1).
The variant, rs4745559 (P = 1.89 × 10−08, β = 0.085), downstream of the GCNT1 gene was the only variant in the GCNT1 locus to reach GWS with an EAF in the transethnic, AFR, AMR, and EUR datasets of 14.6, 17.3, 20.4 and 12.8%, respectively (Supplemental Data 2, http://links.lww.com/QAD/B885, Table 1, http://links.lww.com/QAD/B883)
In addition to these variants identified in the transethnic analysis, rs78474816, intronic to the Cation Channel Sperm Associated Auxiliary Subunit Gamma (CATSPERG) gene reached GWS in the AFR group (P = 4.22 × 10−12, β = 1.379) and exhibited a high level of heterogeneity (HetPVal = 2.12 × 10−11) across ethnicities with AMR and EUR effect sizes equal to 0.185 and −0.021, respectively (Supplemental Data 2, http://links.lww.com/QAD/B885, Supplemental Figure 3A, http://links.lww.com/QAD/B883).
Single nucleotide polymorphisms associated with interleukin-6 levels
Forty variants associated with IL-6 levels were identified within the IL6 Receptor (IL6R) locus (Fig. 4c, Supplemental Data 3, http://links.lww.com/QAD/B886). The lead SNP, rs4133213 (P = 5.84 × 10−17, β = 0.100), is an intronic variant in IL6R. EAF for transethnic, AFR, AMR, and EUR groups were 39.4, 21, 58.2, and 43.6, respectively (Table 1). Additionally, a subset of the variants (14 out of the 40), in strong linkage disequilibrium (LD) with each other, exhibited an opposite effect for their respective effect alleles and were in lower LD across ethnicities as compared with the lead and other variants within the locus. However, conditional analysis placed all 40 variants in the same group. LDpair revealed that in most cases (884 out of 1012), the minor allele of rs4133213(A) (top positive β) commigrates with the major allele of rs6427658(C) (top negative β), explaining the opposite effect of the minor alleles for the two variants (Supplemental Figure 5C, http://links.lww.com/QAD/B883). Due to the small sample size, IL6R could not be identified in AFR and AMR groups (Supplemental Data 3, http://links.lww.com/QAD/B886, Supplemental Figure 4, http://links.lww.com/QAD/B883).
In addition, rs76497186, downstream of noncoding RNA (ncRNA), PROX1-AS1 at chr1 : 213816733, reached GWS in the AFR group (P = 1.19 × 10−08, β = 0.823) and exhibited a high level of heterogeneity (HetPVal = 1.46 × 10−07) across ethnicities with AMR and EUR effect sizes equal to 0.0311 and −0.0086, respectively. An intronic variant in TRAPPC9, rs28368302, reached GWS in the AMR group (P = 3.43 × 10−08, β = −0.232) and also exhibited a high level of heterogeneity (HetPVal = 4.32 × 10−07) across ethnicities with AFR and EUR effect sizes equal to 0.025 and 0.004, respectively (Supplemental Data 3, http://links.lww.com/QAD/B886, Supplemental Figure 4A and B, http://links.lww.com/QAD/B883).
In this study, multiple loci associated with hsCRP, D-dimer and IL-6 levels have been identified in PLWH. Data from three ethnically diverse HIV studies were combined to increase statistical power (Fig. 1). To mitigate the effect of confounding factors, we included age, sex, natural log transformed CD4+ cell count, natural log transformed viral load, BMI, smoker at baseline (missing in ESPRIT), hepatitis B virus (HBV), hepatitis C virus (HCV), diabetes and history of CVD as covariates in the GWAS analyses and many of these covariates have statistically significant effects on the biomarker levels of interest (Supplemental Figure 6, http://links.lww.com/QAD/B883).
A total of 183 variants within 8 loci associated with hsCRP, D-dimer, and IL-6 levels were identified. Not surprisingly, some lead variants in the transethnic group also reached GWS in EUR and/or AFR groups as EUR and AFR participants make up 66 and 25% of the total, respectively. These variants are more significant in the transethnic group than EUR/AFR groups in most cases except for the first SAS in the CRP locus associated with hsCRP level (Table 1). In addition, seven EUR group-specific variants associated with levels of hsCRP (six within HNF1A and one at the CRP locus), one AFR group-specific SNP (rs78474816) associated with D-dimer levels, as well as one AFR group-specific SNP (rs76497186) and one AMR-specific SNP (rs28368302) associated with IL-6 levels have been identified. All of these identified loci and the pathways and mechanism relevant to the regulation of biomarker levels were summarized in Supplemental Table 2, http://links.lww.com/QAD/B883.
Three loci (CRP, HNF1A, and APOE) associated with hsCRP levels in PLWH were identified. Previously, five loci (LEPR, IL6R, CRP, HNF1A, and APOE-CI-CII cluster) were identified from 17 967 participants  and 18 loci (metabolic related: APOC1, HNF1A, LEPR, GCKR, HNF4A, and PTPN2; immune system: CRP, IL6R, NLRP3, IL1F10, and IRF1; chronic inflammation: PPP1R3B, SALL1, PABPC4, ASCL1, RORA, and BCL7B) were identified from more than 80 000 participants  in the general population. Less loci have been identified in this HIV study, which could be the result of small sample size, leading to lower statistical power. HNF1A is a transcription factor, which can bind to the CRP promoter and the common coding variants of the HNF1A gene have been proven to be associated with multiple cardiovascular risk phenotypes [27,44,45]. APOE regulates CRP levels by complexing with activated C1q  and APOE gene variants have been shown to be related to coronary heart disease  (Supplemental Table 2, http://links.lww.com/QAD/B883).
In the current study, four loci (F3, F5, FGB, and GCNT1) associated with D-dimer levels were identified in PLWH. Previously, Smith et al. reported that three genes (F3, F5, and FGA) were associated with D-dimer levels in healthy adults. Lange et al. reported the FGG-10034C/T variant, was associated with higher plasma D-dimer levels in European-American adults. Our findings are consistent with these reports and extend them with the identification of FGB and GCNT1. FGB belongs to the same family as FGA and FGG and they are next to each other around 4q31.3. F3, F5, and FGB (FGA, FGG) belong to the complement and coagulation cascades pathway and regulate D-dimer and fibrinogen levels [24,47–50]. Another novel locus, GCNT1, is a member of the β-1,6-N-acetylglucosaminyltransferase gene family and GCNT1 knockout mice have been shown to have increased susceptibility to Mycobacterium tuberculosis infection and the complete deficiency of GCNT1 was associated with increased lung expression of the neutrophil chemoattractant CXCL2 . Moreover, Notch signaling regulates Gcnt1-mediated core-2 O-glycosylation in activated T cells and the core-2 O-glycoform of CD43 could be a sensitive indicator of Notch signaling . Furthermore, GCNT1 expression in prostate cancer positively correlates with cancer progression and prostate-specific antigen recurrence  (Supplemental Table 2, http://links.lww.com/QAD/B883). In addition, rs78474816 within CATSPERG was associated with D-dimer levels in the AFR but not the transethnic group. Although this variant was detected in a limited number of participants, it still reached GWS because of its strong effect. Although CATSPERG is associated with CATSPER1 channel protein, only expressed in testis , seems irrelevant to D-dimer levels, rs78474816 is also located upstream of PSMD8 (Proteasome 26S Subunit, non-ATPase 8), which encodes a non-ATPase subunit of the 19S regulator. As the proteasome participates in the protein degradation process  and HIV can interact with proteasome via Tat [56–58], vif [59–62] and integrase , this variant might play an important role in regulation of D-dimer in PLWH.
Forty variants, led by rs4133213, within the IL6R locus were associated with IL-6 levels in the transethnic group. The variant rs4133213 has previously been identified associated with plasma sIL-6R levels with rs2228145 as the lead variant . Interestingly, we found that variant rs2228145 was in high LD with the lead SNP, rs4133213 (Supplemental Figure 7, http://links.lww.com/QAD/B883) and also highly significant in our study. SNP rs2228145 is a missense variant (p.ASP358Ala) within IL6R. It has been reported to modulate IL-6 levels using a total of 1979 older Chinese individuals aged 50–92 years , associated with IL-6 levels and sIL-6R in both African Americans and European Americans  and in a study of 1273 participants from the InCHIANTI Italian cohort . The increased level of IL-6 could be the result of increased levels of sIL-6R as the binding of IL-6 by sIL-6R could protect IL-6 from degradation. When IL-6 binds to sIL-6R, it can trigger the trans-signaling pathway, which is responsible for the pro-inflammatory action of IL-6 . Blockade of the IL6R signaling seems to have a causal role in the development of coronary heart disease in the general population and inhibition of IL6R is proposed as a target for prevention of CHD [68,69]. In addition to the IL6R locus identified in the transethnic group, two novel ethnic-specific loci associated with IL-6 levels have been identified: one variant (rs76497186) downstream of a ncRNA, PROX1-AS1, in the AFR group and one intronic variant (rs28368302) within TRAPPC9 in the AMR group. PROX1_AS1 is located upstream of PROX1, which is a homeodomain transcription factor expressed in various tissues during mouse development including young neurons of the subventricular region of the CNS, developing eye lens, pancreas, liver, heart, and transiently in the skeletal muscles . Prox1 and vascular endothelial growth factor receptor-3 (VEGFR-3) are two primary mediators of lymphangiogenesis. Prox1 can be activated by NF-kB pathway, which is induced by inflammatory stimuli. Prox1 can then increase the expression of the VEGFR-3 in lymphatic endothelial cells, leading to the enhanced responsiveness of preexisting lymphatic endothelium to VEGFR-3-specific VEGF (VEGF-C and VEGF-D) . TRAPPC9 has been shown to enhance the NF-kB signaling pathway [72–74] and IL-6 gene expression can be activated through the NF-kB transcription factor . IL-6 was previously found to be correlated with blood pressure in male individuals with essential hypertension  and hypomethylation of IL-6 promoter was associated with prehypertension in young adults . Consistently, two intronic SNPs within TRAPPC9 (rs10088725 and rs6578061) were identified to be associated with blood pressure within 750 000 trans-ethnic subjects . It is possible that TRAPPC9 regulates blood pressure via regulation of IL-6 expression. In addition, the variant rs11166927 in the TRAPPC9 region has also been identified to be associated with nonalcoholic fatty liver disease in Hispanic boys . This could be caused by the effect of TRAPPC9 on IL-6 via NF-κB.
In conclusion, multiple SNPs were associated with levels of hsCRP, D-dimer, and IL-6 in PLWH from three ethnically diverse studies. Most loci have been identified in the general population. However, five novel loci have also been found to be associated with D-dimer levels (FGB and GCNT1 in the transethnic group, and CATSPERG in the AFR group) and IL-6 levels (PROX1-AS1 in the AFR group and TRAPPC9 in the AMR group). These findings support the hypothesis that host genetics may partially contribute to chronic inflammation in HIV+ individuals. Given the link between these biomarkers and serious non-AIDS complications, further exploration of the relationship between these 11 loci and the events may help to identify potential targets for intervention.
We thank Ms Ju Qiu and Dr Xiaoli Jiao in LHRI for scientific discussion, Dr Norman Gerry in Advanced BioMedical Laboratory for generation of genotyping data, and Ms Shweta Sharma Mistry, Ms Jacqueline Nordwall, Ms Deborah Wentworth and Mr Greg Thompson at the University of Minnesota for providing participant biomarker and other clinically relevant data.
This project was funded in whole or in part with federal funds from the National Cancer Institute, NIH under contract no. HHSN261200800001E. The project is also supported in part by other NIH grants UM1-AI068641, UM1-AI120197, 1U01-AI36780, U01-AI046362, U01-AI042170 and U01-AI46957.
The content of this publication does not necessarily reflect the views or policies of the US Department of Health and Human Services nor does mention of trade names, commercial products, or organizations imply endorsement by the US government. This research was supported in part by the National Institute of Allergy and Infectious Disease.
B.S. and X.H. analyzed data and wrote manuscript draft, K.S. intellectually involved in the project design and wrote manuscript draft, J.N. and J.L. intellectually involved in the project design, L.H. analyzed data and wrote the draft manuscript, A.W.R. measured biomarker levels, T.I. was intellectually involved in the project and reviewed the manuscript, W.C. oversaw data analysis and wrote the manuscript, H.C.L. designed, oversaw the project and wrote the manuscript. We are grateful to all study participants and all investigators of the ESPRIT, SMART and START Study groups. Please see N Engl J Med 2009; 361 : 1548-59 for the complete list of ESPRIT investigators; N Engl J Med 2006; 355 : 2283-96 for the complete list of SMART investigators; and N Engl J Med 2015; 373 : 795-807 for the complete list of START investigators.
Conflicts of interest
There are no conflicts of interest.
1. Antiretroviral Therapy Cohort Collaboration. Causes of death in HIV-1-infected patients treated with antiretroviral therapy, 1996–2006: collaborative analysis of 13 HIV cohort studies
. Clin Infect Dis
2. Lundgren JD, Babiker AG, Gordin F, Emery S, Grund B, et al. INSIGHT START Study Group. Initiation of antiretroviral therapy in early asymptomatic HIV infection
. N Engl J Med
3. Kuller LH, Tracy R, Belloso W, De Wit S, Drummond F, Lane HC, et al. INSIGHT SMART Study Group. Inflammatory and coagulation biomarkers and mortality in patients with HIV infection
. PLoS Med
4. Grund B, Baker JV, Deeks SG, Wolfson J, Wentworth D, Cozzi-Lepri A, et al. INSIGHT SMART/ESPRIT/SILCAAT Study Group. Relevance of interleukin-6 and d-dimer for serious non-AIDS morbidity and death among HIV-positive adults on suppressive antiretroviral therapy
. PLoS One
5. Danesh J, Kaptoge S, Mann AG, Sarwar N, Wood A, Angleman SB, et al. Long-term interleukin-6 levels and subsequent risk of coronary heart disease: two new prospective studies and a systematic review
. PLoS Med
6. Danesh J, Whincup P, Walker M, Lennon L, Thomson A, Appleby P, et al. Fibrin D-dimer and coronary heart disease: prospective study and meta-analysis
7. Kaptoge S, Di Angelantonio E, Lowe G, Pepys MB, Thompson SG, et al. The Emerging Risk Factors Collaboration. C-reactive protein concentration and risk of coronary heart disease, stroke, and mortality: an individual participant meta-analysis
8. Li Y, Zhong X, Cheng G, Zhao C, Zhang L, Hong Y, et al. Hs-CRP and all-cause, cardiovascular, and cancer mortality risk: a meta-analysis
9. Gilotra TS, Geraci SA. C-reactive protein as an independent cardiovascular risk predictor in HIV+ patients: a focused review of published studies
. J Clin Med Res
10. Savvoulidis P, Butler J, Kalogeropoulos A. Cardiomyopathy and heart failure in patients with HIV infection
. Can J Cardiol
11. Nasi M, De Biasi S, Gibellini L, Bianchini E, Pecorini S, Bacca V, et al. Ageing and inflammation in patients with HIV infection
. Clin Exp Immunol
12. Baker JV, Sharma S, Grund B, Rupert A, Metcalf JA, Schechter M, et al. Systemic inflammation, coagulation, and clinical risk in the START Trial
. Open Forum Infect Dis
13. Neuhaus J, Jacobs DR Jr, Baker JV, Calmy A, Duprez D, La Rosa A, et al. Markers of inflammation, coagulation, and renal function are elevated in adults with HIV infection
. J Infect Dis
14. Bennermo M, Held C, Ericsson CG, Silveira A, Hamsten A, Tornvall P. Genotype-specific increase in plasma concentrations of activated coagulation factor VII in response to experimental inflammation. A link between infection and acute myocardial infarction?
. Thromb Haemost
15. Lange LA, Reiner AP, Carty CL, Jenny NS, Cushman M, Lange EM. Common genetic variants associated with plasma fibrin D-dimer concentration in older European- and African-American adults
. J Thromb Haemost
16. Ekenberg C, Tang MH, Zucco AG, Murray DD, MacPherson CR, Hu X, et al. Association between single-nucleotide polymorphisms in HLA alleles and human immunodeficiency virus type 1 viral load in demographically diverse, antiretroviral therapy-naive participants from the Strategic Timing of AntiRetroviral Treatment Trial
. J Infect Dis
17. McLaren PJ, Coulonges C, Bartha I, Lenz TL, Deutsch AJ, Bashirova A, et al. Polymorphisms of large effect explain the majority of the host genetic contribution to variation of HIV-1 virus load
. Proc Natl Acad Sci U S A
18. van Manen D, van ’t Wout AB, Schuitemaker H. Genome-wide association studies on HIV susceptibility, pathogenesis and pharmacogenomics
19. Mayerle J, den Hoed CM, Schurmann C, Stolk L, Homuth G, Peters MJ, et al. Identification of genetic loci associated with Helicobacter pylori serologic status
20. Wojcik GL, Marie C, Abhyankar MM, Yoshida N, Watanabe K, Mentzer AJ, et al. Genome-wide association study reveals genetic link between diarrhea-associated entamoeba histolytica infection and inflammatory bowel disease
21. Nagao Y, Nishida N, Toyo-Oka L, Kawaguchi A, Amoroso A, Carrozzo M, et al. Genome-wide association study identifies risk variants for Lichen Planus in patients with hepatitis C virus infection
. Clin Gastroenterol Hepatol
22. Elliott P, Chambers JC, Zhang W, Clarke R, Hopewell JC, Peden JF, et al. Genetic Loci associated with C-reactive protein levels and risk of coronary heart disease
23. Dehghan A, Dupuis J, Barbalic M, Bis JC, Eiriksdottir G, Lu C, et al. Meta-analysis of genome-wide association studies in >80 000 subjects identifies multiple loci for C-reactive protein levels
24. Smith NL, Huffman JE, Strachan DP, Huang J, Dehghan A, Trompet S, et al. Genetic predictors of fibrin D-dimer levels in healthy adults
25. Jiang CQ, Lam TH, Liu B, Lin JM, Yue XJ, Jin YL, et al. Interleukin-6 receptor gene polymorphism modulates interleukin-6 levels and the metabolic syndrome: GBCS-CVD
. Obesity (Silver Spring)
26. Nimptsch K, Aleksandrova K, Boeing H, Janke J, Lee YA, Jenab M, et al. Association of CRP genetic variants with blood concentrations of C-reactive protein and colorectal cancer risk
. Int J Cancer
27. Odom DT, Zizlsperger N, Gordon DB, Bell GW, Rinaldi NJ, Murray HL, et al. Control of pancreas and liver gene expression by HNF transcription factors
28. Yin C, Ackermann S, Ma Z, Mohanta SK, Zhang C, Li Y, et al. ApoE attenuates unresolvable inflammation by complex formation with activated C1q
. Nat Med
29. Besser REJ, Jones J, McDonald TJ, Smith R, Shepherd MH, Hattersley AT. Using highly sensitive C-reactive protein measurement to diagnose MODY in a family with suspected type 2 diabetes
. BMJ Case Rep
30. INSIGHT-ESPRIT Study Group, SILCAAT Scientific Committee, Abrams D, Levy Y, Losso MH, Babiker A, et al. Interleukin-2 therapy in patients with HIV infection
. N Engl J Med
31. El-Sadr WM, Lundgren J, Neaton JD, Gordin F, Abrams D, et al. Strategies for Management of Antiretroviral Therapy (SMART) Study Group. CD4+ count-guided interruption of antiretroviral treatment
. N Engl J Med
32. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses
. Am J Hum Genet
33. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals
. Genome Res
34. Ihaka R, Gentleman R. R: a language for data analysis and graphics
. J Comput Graph Stat
35. Borges AH, O’Connor JL, Phillips AN, Baker JV, Vjecha MJ, Losso MH, et al. Factors associated with D-dimer levels in HIV-infected individuals
. PLoS One
36. Borges AH, O’Connor JL, Phillips AN, Ronsholt FF, Pett S, Vjecha MJ, et al. Factors associated with plasma IL-6 levels during HIV infection
. J Infect Dis
37. De Luca A, de Gaetano Donati K, Colafigli M, Cozzi-Lepri A, De Curtis A, Gori A, et al. The association of high-sensitivity c-reactive protein and other biomarkers with cardiovascular disease in patients treated for HIV: a nested case-control study
. BMC Infect Dis
38. Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans
39. Yang J, Weedon MN, Purcell S, Lettre G, Estrada K, Willer CJ, et al. GIANT Consortium. Genomic inflation factors under polygenic inheritance
. Eur J Hum Genet
40. Yin L, Kant DH. Available at: https://github.com/YinLiLin/R-CMplot
. [Accessed 26 April 2019]
41. Vuong H, Che A, Ravichandran S, Luke BT, Collins JR, Mudunuri US. AVIA v2.0: annotation, visualization and impact analysis of genomic variants and genes
42. Ning Z, Lee Y, Joshi PK, Wilson JF, Pawitan Y, Shen X. A selection operator for summary association statistics reveals allelic heterogeneity of complex traits
. Am J Hum Genet
43. Machiela MJ, Chanock SJ. LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants
44. Reiner AP, Gross MD, Carlson CS, Bielinski SJ, Lange LA, Fornage M, et al. Common coding variants of the HNF1A gene are associated with multiple cardiovascular risk phenotypes in community-based samples of younger and older European-American adults: the Coronary Artery Risk Development in Young Adults Study and The Cardiovascular Health Study
. Circ Cardiovasc Genet
45. Armendariz AD, Krauss RM. Hepatic nuclear factor 1-alpha: inflammation, genetics, and atherosclerosis
. Curr Opin Lipidol
46. Xu M, Zhao J, Zhang Y, Ma X, Dai Q, Zhi H, et al. Apolipoprotein E gene variants and risk of coronary heart disease: a meta-analysis
. Biomed Res Int
47. Raffield LM, Zakai NA, Duan Q, Laurie C, Smith JD, Irvin MR, et al. NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium, Hematology & Hemostasis TOPMed Working Group∗. D-Dimer in African Americans: whole genome sequence analysis and relationship to cardiovascular disease risk in the Jackson Heart Study
. Arterioscler Thromb Vasc Biol
48. Zaleski AL, Pescatello LS, Thompson PD, Taylor BA. Protective effect of compression socks in a marathon runner with a genetic predisposition to thrombophilia due to Factor V Leiden
. Phys Sportsmed
49. de Moerloose P, Casini A, Neerman-Arbez M. Congenital fibrinogen disorders: an update
. Semin Thromb Hemost
50. Danik JS, Pare G, Chasman DI, Zee RY, Kwiatkowski DJ, Parker A, et al. Novel loci, including those related to Crohn disease, psoriasis, and inflammation, identified in a genome-wide association study of fibrinogen in 17 686 women: the Women's Genome Health Study
. Circ Cardiovasc Genet
51. Fonseca KL, Maceiras AR, Matos R, Simoes-Costa L, Sousa J, Cá B, et al. Deficiency in the glycosyltransferase Gcnt1 increases susceptibility to tuberculosis through a mechanism involving neutrophils
. Mucosal Immunol
52. Perkey E, Maurice De Sousa D, Carrington L, Chung J, Dils A, Granadier D, et al. GCNT1-mediated O-glycosylation of the Sialomucin CD43 is a sensitive indicator of notch signaling in activated T cells
. J Immunol
53. Kojima Y, Yoneyama T, Hatakeyama S, Mikami J, Sato T, Mori K, et al. Detection of Core2 beta-1,6-N-acetylglucosaminyltransferase in post-digital rectal examination urine is a reliable indicator for extracapsular extension of prostate cancer
. PLoS One
54. Wang H, Liu J, Cho KH, Ren D. A novel, single, transmembrane protein CATSPERG is associated with CATSPER1 channel protein
. Biol Reprod
55. Driscoll J, Finley D. A controlled breakdown: antigen processing and the turnover of viral proteins
56. Seeger M, Ferrell K, Frank R, Dubiel W. HIV-1 tat inhibits the 20 S proteasome and its 11 S regulator-mediated activation
. J Biol Chem
57. Huang X, Seifert U, Salzmann U, Henklein P, Preissner R, Henke W, et al. The RTP site shared by the HIV-1 Tat protein and the 11S regulator subunit alpha is crucial for their effects on proteasome function including antigen processing
. J Mol Biol
58. Apcher GS, Heink S, Zantopf D, Kloetzel PM, Schmid HP, Mayer RJ, et al. Human immunodeficiency virus-1 Tat protein interacts with distinct proteasomal alpha and beta subunits
. FEBS Lett
59. Conticello SG, Harris RS, Neuberger MS. The Vif protein of HIV triggers degradation of the human antiretroviral DNA deaminase APOBEC3G
. Curr Biol
60. Lecossier D, Bouchonnet F, Clavel F, Hance AJ. Hypermutation of HIV-1 DNA in the absence of the Vif protein
61. Madani N, Kabat D. An endogenous inhibitor of human immunodeficiency virus in human lymphocytes is overcome by the viral Vif protein
. J Virol
62. Simon JH, Gaddis NC, Fouchier RA, Malim MH. Evidence for a newly discovered cellular anti-HIV-1 phenotype
. Nat Med
63. Mulder LC, Muesing MA. Degradation of HIV-1 integrase by the N-end rule pathway
. J Biol Chem
64. van Dongen J, Jansen R, Smit D, Hottenga JJ, Mbarek H, Willemsen G, et al. The contribution of the functional IL6R polymorphism rs2228145, eQTLs and other genome-wide SNPs to the heritability of plasma sIL-6R levels
. Behav Genet
65. Reich D, Patterson N, Ramesh V, De Jager PL, McDonald GJ, Tandon A, et al. Health, Aging and Body Composition (Health ABC) Study. Admixture mapping of an allele affecting interleukin 6 soluble receptor and interleukin 6 levels
. Am J Hum Genet
66. Rafiq S, Frayling TM, Murray A, Hurst A, Stevens K, Weedon MN, et al. A common variant of the interleukin 6 receptor (IL-6r) gene increases IL-6r and IL-6 levels, without other inflammatory effects
. Genes Immun
67. Schett G. Physiological effects of modulating the interleukin-6 axis
. Rheumatology (Oxford)
2018; 57: (Suppl_2): ii43–ii50.
68. Sarwar N, Butterworth AS, Freitag DF, Gregson J, Willeit P, et al. IL6R Genetics Consortium Emerging Risk Factors Collaboration. Interleukin-6 receptor pathways in coronary heart disease: a collaborative meta-analysis of 82 studies
69. Swerdlow DI, Holmes MV, Kuchenbaecker KB, Engmann JE, Shah T, Sofat R, et al. Interleukin-6
Receptor Mendelian Randomisation Analysis (IL6R MR) Consortium. The interleukin-6 receptor as a target for prevention of coronary heart disease: a mendelian randomisation analysis
70. Oliver G, Sosa-Pineda B, Geisendorf S, Spana EP, Doe CQ, Gruss P. Prox 1, a prospero-related homeobox gene expressed during mouse development
. Mech Dev
71. Flister MJ, Wilber A, Hall KL, Iwata C, Miyazono K, Nisato RE, et al. Inflammation induces lymphangiogenesis through up-regulation of VEGFR-3 mediated by NF-kappaB and Prox1
72. Zhang Y, Liu S, Wang H, Yang W, Li F, Yang F, et al. Elevated NIBP/TRAPPC9 mediates tumorigenesis of cancer cells through NF(B signaling
73. Qin M, Zhang J, Xu C, Peng P, Tan L, Liu S, et al. Knockdown of NIK and IKKβ0-binding protein (NIBP) reduces colorectal cancer metastasis through down-regulation of the canonical NF-κB signaling pathway and suppression of MAPK signaling mediated through ERK and JNK
. PLoS One
74. Hu WH, Pendergast JS, Mo XM, Brambilla R, Bracchi-Ricard V, Li F, et al. NIBP, a novel NIK and IKK(beta)-binding protein that enhances NF-(kappa)B activation
. J Biol Chem
75. Libermann TA, Baltimore D. Activation of interleukin-6 gene expression through the NF-kappa B transcription factor
. Mol Cell Biol
76. Gibas-Dorna M, Nowak D, Piatek J, Pupek-Musialik D, Krauss H, Kopczynski P. Plasma ghrelin and interleukin-6 levels correlate with body mass index and arterial blood pressure in males with essential hypertension
. J Physiol Pharmacol
77. Omar W, Abdullah A, Talib NA, Shah ASM, Rahman JA. Leucocytic DNA methylation of interleukin-6 promoter reduction in pre-hypertensive young adults
. Malays J Med Sci
78. Giri A, Hellwege JN, Keaton JM, Park J, Qiu C, Warren HR, et al. Trans-ethnic association study of blood pressure determinants in over 750,000 individuals
. Nat Genet
79. Wattacheril J, Lavine JE, Chalasani NP, Guo X, Kwon S, Schwimmer J, et al. Genome-wide associations related to hepatic histology in nonalcoholic fatty liver disease in Hispanic boys
. J Pediatr