|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Molecular Diagnostics and Genetics |
1 Centre for Cardiovascular Genetics, Department of Medicine, University College London, London, UK; 2 Periodontology Unit, UCL Eastman Dental Institute, London, UK; 3 Centre for Paediatric and Adolescent Rheumatology, Windeyer Institute, London, UK; 4 Department of Surgery, University College London, The Heart Hospital, London, UK; 5 Department of Cardiology, The Whittington Hospital NHS Trust, London, UK.
aAddress correspondence to this author at: Centre for Cardiovascular Genetics, Department of Medicine, University College London, 5 University St., London WC1E 6JF, UK. Fax 44-(0)20-7679-6212; e-mail Andrew.J.P.Smith{at}ucl.ac.uk.
| Abstract |
|---|
|
|
|---|
Methods: We sequenced a region 6-kb upstream of the IL6 [interleukin 6 (interferon, beta 2)] transcription start site in a search for functional variants and detected 3 common variants: –6331T>C, –6101A>T, and –5617/–5616C/A>T/G. IL6 –6331T>C (C allele frequency, 0.20; 95% confidence interval, 0.16–0.24) showed strong negative linkage disequilibrium with –174G>C (D' = –0.97) and was studied further in 309 individuals who underwent coronary artery bypass grafting.
Results: Patients with the TT genotype had higher IL-6 concentrations 6 h after surgery than those with the CC genotype (mean, 199.4 ng/L vs 114.9 ng/L; P = 0.02). A similar association was seen in a cohort of 173 patients who underwent intensive periodontal therapy: Individuals with the CC genotype had significantly lower IL-6 concentrations 24 h after therapy than TT patients (mean, 0.78 ng/L vs 5.00 ng/L; P < 0.0001). A similar trend was observed in 203 healthy individuals from northern Europe (1.29 ng/L for the TT genotype vs 0.89 ng/L for the CC genotype; P = 0.07). Reporter assays that used a sequence flanking the –6331 single-nucleotide polymorphism spliced upstream to the IL-6 minimal promoter driving luciferase gene expression demonstrated a 1.3-fold increase in promoter activity (P < 0.01) for constructs containing –6331T. Electrophoretic mobility shift assays revealed enhanced binding of transcription factor Oct-1 to the T allele.
Conclusions: IL6 –6331T is associated with increased IL-6 concentrations in an acute inflammatory state via a mechanism involving binding of the Oct-1 transcription factor. This finding may help resolve conflicting studies based on the IL6 –174G>C variant.
| Introduction |
|---|
|
|
|---|
, but the fact that IL-6 remains in the plasma substantially longer makes this molecule a good marker of inflammation(2). Increased plasma IL-6 concentrations act as prognostic markers in several clinical conditions, including intra-abdominal sepsis(3), metastatic breast cancer(4), and diffuse large cell lymphoma(5). The development and severity of coronary artery disease are also associated with increased IL-6 concentrations(6), as is the transition to plaque instability(7), possibly because of the action of IL-6 in vascular smooth muscle cell proliferation, leukocyte recruitment, and endothelial activation.
IL-6 binds to cell surface complexes consisting of the IL-6 receptor and gp130, thereby activating the JAK/STAT and MAPK cascades. Located on chromosome 7p21, the IL62 gene [interleukin 6 (interferon, beta 2)] contains 5 exons, and the IL6 coding sequence contains only a few very rare polymorphisms (http://genome.ucsc.edu/, 2006 assembly). This fact has led to the hypothesis that the observed differences in IL-6 concentration among individuals are influenced by variation in the IL6 promoter in combination with variation in genes upstream of IL-6 synthesis and environmental stimuli, such as a smoking habit.
IL6 promoter polymorphisms have been well characterized in a number of healthy and disease populations, and 2 single-nucleotide polymorphisms (SNPs) have been studied in detail—the common –174G>C variant(8) and the less frequent –572G>C allele(9). The –174C allele has been associated with higher IL-6 serum concentrations in several studies(10)(11)(12), with no effect in 1 study(13), and with lower concentrations in another(14). The –572G>C SNP has also been associated with the plasma IL-6 concentration, with the rare C allele associated with higher IL-6 concentrations(15). A recent study also evaluated 2 other promoter variants, –1363G>T and –1480CTdel, and although these variants were not shown to be functional, they did form a haplotype with the other promoter polymorphisms that showed an increased association with systemic-onset juvenile arthritis(16). This finding suggests the presence of functional polymorphisms at this locus that influence IL-6 concentration but that have yet to be identified.
Our recent work has identified a region between –5202 and –5307 bp upstream of the IL6 transcription start site that plays an important role in IL6 gene regulation. In conjunction with the proximal promoter, this region increases gene expression in vitro by up to 20-fold compared with the proximal promoter alone. The aim of the present study was to examine whether any functional genetic variation is present near this region that affects IL-6 concentration.
| Materials and Methods |
|---|
|
|
|---|
| Results |
|---|
|
|
|---|
|
genotyping of il6 snpS and linkage disequilibrium analysis
We genotyped these 4 SNPs for a UK cohort of patients who underwent coronary artery bypass grafting surgery (CABG)(10), a cohort of patients who received intensive periodontal therapy (IPT)(17)(18)(19)(20), and the control population from a European case-control study of myocardial infarction (HIFMECH)(21). The baseline characteristics of the patients (CABG, n = 321; IPT, n = 173; HIFMECH, n = 421) for whom IL-6 concentrations and at least 1 genotype were measured are shown in Table 1
.
|
We successfully genotyped 309 of the UK CABG patients for all IL6 SNPs investigated in this study: IL6 –6331T>C, –6101A>T, –5617C>T, –5616A>G, –572G>C, and –174G>C. All genotype distributions were in Hardy–Weinberg equilibrium at the following allele frequencies: –6331T>C, 0.20; –6101A>T, 0.20; –5617C>T, 0.36; –5616A>G, 0.36; –572G>C, 0.08; and –174G>C, 0.36. These frequencies were similar to those reported in HapMap (http://www.hapmap.org/; data not shown). Genotyping of the contiguous –5617C>T and –5616A>G alleles revealed the presence of 2 naturally occurring diplotypes, C/A and T/G, which were subsequently analyzed as a diallelic SNP. There was complete linkage disequilibrium (LD) between –6331T>C and –6101A>T, and subsequent analysis focused only on –6331T>C (C allele frequency, 0.20; 95% confidence interval, 0.16–0.24), because an in silico analysis revealed that the T>C change eliminates an Oct-1 binding site.
effect of genotype on plasma il-6 concentration in cabg patients
IL-6 concentrations and the IL6 genotype at baseline, 6 h, and 24 h before CABG are shown in Table 1A in the online Data Supplement. A previous study examined 127 DNA samples obtained from the same cohort used in the present study and identified an association of –572G>C and –174G>C with the plasma IL-6 concentration at 6 h(10). With current data on IL-6 concentrations and at least one IL6 genotype measured for 321 patients, we reexamined these 2 SNPs along with the new SNPs described above for an association with IL-6 concentration (variation in sample size reflects missing IL-6 concentration and/or genotype data, because IL-6 concentration data were not available for all patients at every time point). These data show a reduction in the significance of the association of previous –174/–572 genotypic data with IL-6 concentration. The –5617C>T/–5616A>G alleles showed no significant association with IL-6 concentration at any time (all P > 0.5; see Table 1A in the online Data Supplement), but a comparison of homozygote genotypes revealed an association with –6331T>C (data not shown in this table). We found a nonsignificant (P = 0.13) trend of association between IL-6 concentration and genotype at baseline, with a 34% higher plasma IL-6 concentration for the TT genotype than for the CC genotype [mean (SE), 4.64 (0.23) ng/L vs 3.45 (0.40) ng/L]. The TT genotype was associated with a 73% higher IL-6 plasma concentration at 6 h compared with the CC genotype [199.4 (12.3) ng/L vs 114.9 (36.9) ng/L; P = 0.02]. The IL-6 concentrations of the TT and CC genotypes were not significantly different at 24 h after surgery.
effect of genotype on plasma il-6 concentration in ipt patients
Serum IL-6 concentrations for the IPT cohort showed no association with the IL6 genotype at the baseline visit (Table 1B in the online Data Supplement). As previously reported, periodontal therapy produced a moderate inflammatory response that lasted 1 week. We have previously shown no effect of –174G>C on IL-6 concentration for a smaller sample of individuals(18). In the present analysis with a larger sample, we confirmed the lack of an association between the –174/–572 genotypic data and acute increases in IL-6 concentration; however, we did observe a significant association between –6331T>C genotype and serum IL-6 concentration at 24 h after therapy, with CC patients exhibiting significantly lower IL-6 concentrations than TT patients [0.78 (0.55) ng/L vs 5.00 (0.57) ng/L; P < 0.0001]. Patients with the CC genotype also had significantly lower IL-6 concentrations at 7 days after therapy than those with the TT genotype [0.40 (0.56) ng/L vs 1.53 (0.24) ng/L; P = 0.02].
effect of genotype on plasma il-6 concentration in the hifmech study
To establish whether the findings in this UK cohort could be replicated, we genotyped individuals from another study (HIFMECH) and selected only the control groups to avoid the confounding effects of the medication (421 of the 575 individuals had measurements of IL-6 concentration and at least one genotype). There was no statistically significant association between genotype and plasma IL-6 concentration (see Table 1C
in the online Data Supplement); however, the TT and CC genotypes at position –6331 were a strong indicator of IL-6 concentration for the northern European cohort (P = 0.07), a result similar to that seen in the baseline data for the CABG cohort. The data for the southern European cohort did not replicate this trend, however; the lowest IL-6 concentrations occurred in heterozygous individuals.
analysis of haplotype and il-6 concentration in cabg patients
For the haplotype analysis, we included only SNPs for which evidence existed for a functional role [–6331T>C (described here), –572G>C, and –174G>C (Fig. 1
)]. There was significant LD between these SNPs, and 4 common haplotypes were observed. The –6331C allele was always observed with –572G and –174G, and –572C was always observed with –174G. After CABG, the patients who carried the haplotype containing –6331C (CGG, 21% of the sample) had IL-6 concentrations that were 34% lower at 6 h than individuals carrying the common haplotype (TGG, 38% of the sample; P = 0.01) and had IL-6 concentrations at 24 h that were 22% lower (P = 0.03; Table 2A
, Fig. 2A
).
|
|
analysis of haplotype and il-6 concentration in ipt patients
The haplotype analysis of the IPT cohort confirmed the findings with the CABG cohort, with a lower acute release of IL-6 being associated with the haplotype containing –6331C. Indeed, the CGG haplotype (21% of the sample) was associated with a 51% reduction in IL-6 concentration compared with the common TGG haplotype (39% of the sample) at 24 h after therapy (P = 0.002; Table 2B
, Fig. 2B
). Because the data for the CABG and IPT cohorts are time series, we repeated the analysis by calculating the area under the IL-6 curve and confirmed that IL-6 concentrations for the CGG haplotype were significantly lower than for the common TGG haplotype in both the CABG cohort (P = 0.03) and the IPT cohort (P = 0.02) (see Table 2 in the online Data Supplement). No effects of any other haplotypes were observed at any of the other times.
analysis of haplotype and il-6 concentration in the hifmech cohort
The results of our haplotype analysis of the HIFMECH cohort repeated the trend found in the CABG study (Table 2C
). IL-6 concentrations were 39% lower in individuals with the CGG haplotype than with the common TGG haplotype, although this difference did not reach statistical significance (P = 0.07). We observed a continuation of this trend, however, in both the northern and southern European cohorts, with 45% and 40% reductions, respectively, in IL-6 concentration relative to the reference haplotype, suggesting a genuine haplotypic effect.
luciferase reporter assay
To evaluate whether the IL6 –6331T>C allele was itself a functional polymorphism or simply a marker for another site elsewhere at this locus, we performed luciferase reporter assays with constructs containing the IL6 minimal promoter and a 196-bp fragment containing the –6331T>C variant at the 5' end. All 4 combinations containing –174G/C and –6331T/C alleles were created. The constructs were transfected into HeLa and Huh7 cell lines, both of which express the IL6 gene (data not shown). Cells were cotransfected with an Oct-1 expression plasmid and stimulated with IL-1β to promote IL-6 production. Fig. 3
shows luciferase expression for each plasmid construct in both stimulated and unstimulated HeLa and Huh7 cell lines, relative to that of an unstimulated wild-type vector (–6331T/–174G) without Oct-1 cotransfection. All constructs that were cotransfected with the Oct-1 expression vector produced significantly higher luciferase activity than those cotransfected with the empty vector (data not shown). There was no significant difference in luciferase activity between any of the constructs in either unstimulated cell line. Compared with the –6331C/–174G construct, IL-1β stimulation significantly increased the expression by the –6331T/–174G construct by 1.2-fold (P = 0.013) and 1.1-fold (P = 0.01) in the HeLa and Huh7 cell lines, respectively (Fig. 3
). Similarly, the –6331T/–174C construct had 1.3-fold (P = 0.0014) and 1.3-fold (P = 0.001) higher expression in HeLa and Huh7 cells, respectively, than the –6331C/–174C construct. Performing the luciferase reporter assay with IL-1β stimulation for 3 h or 6 h did not alter the observed relative differences in luciferase expression (data not shown).
|
electrophoretic mobility shift assay
Our in silico analysis of the wild-type and variant sequences with MatInspector (Genomatix Software) indicated that transcription factor Oct-1 binds to the T allele but not to the C allele at position –6331. To investigate the effect of this SNP on transcription factor binding, we performed an EMSA analysis on a 22-bp sequence surrounding the SNP. We used a known Oct-1 consensus sequence as a positive control and carried out supershift EMSA assays with an Oct-1–specific monoclonal antibody to confirm binding specificity. As shown in Fig. 4
, both the Oct-1 consensus sequence and the –6331T allele probe were shifted after incubation with nuclear extract from IL-1β–stimulated HeLa cells, and we observed a supershift after adding Oct-1–specific monoclonal antibody. The –6331C allele probe was not shifted after incubation with HeLa nuclear extract.
|
| Discussion |
|---|
|
|
|---|
We examined a UK CABG cohort, hypothesizing that the increase in plasma IL-6 concentration following this proinflammatory insult would increase our power to detect an association with SNPs that may play a role in producing the documented variation in IL-6 concentration. This study identified the –6331T>C SNP to be an important predictor of IL-6 concentration at 6 h after surgery and demonstrated the TT genotype to be associated with a 74% higher IL-6 concentration than the CC genotype. The genotype ranking of IL-6 concentrations at baseline was similar to that at 6 h after surgery, but the differences did not reach statistical significance. It is notable that the patients with the CC genotype eventually reached the same IL-6 concentrations by 24 h after surgery as the TT and CT genotypes had reached by 6 h. Therefore, the –6331 genotype may be not only a predictor of IL-6 concentration after acute stimulation but also a predictor of the rate at which the IL-6 concentration increases after stimulation. To confirm these findings, we repeated this analysis with a cohort of individuals who underwent IPT. Although IPT patients do not constitute such a severe model of acute inflammation, our group had previously characterized the IPT model to be associated with a significant inflammatory response at 1 week. We previously had shown in a smaller IPT cohort a positive association between common SNPs and the acute release of other inflammatory mediators(29). In the IPT model, we confirmed a strong association of the –6331C allele and haplotypes containing this allele with the release of substantially lower serum IL-6 concentrations after a moderate inflammatory stimulus. As in the CABG cohort, we observed no association between any of the SNPs and baseline IL-6 values.
We carried out the same analysis with the control group of a case-control study of myocardial infarction, in which IL-6 was measured for 6 months following myocardial infarction. We confirmed the –6331T>C SNP as a predictor of IL-6 concentration in the northern Europe cohort but not in southern Europeans, possibly because of diet, lifestyle differences, and genetic interactions. Southern Europeans have significantly lower systolic blood pressures and cholesterol and triglyceride concentrations (Table 1A
) and are genetically distinct. As with the CABG study at baseline, the genotype/phenotype correlation did not reach statistical significance, indicating that significant concentrations of phenotypic variation are dependent on acute inflammatory stimuli. The chronic inflammatory response to smoking had little effect on –6331T>C and IL-6 concentrations in all of the cohorts, and we also noted no association between any IL6 polymorphism studied and the concentration of C-reactive protein (data not shown). Because the numbers in some groups were small, we assessed differences in IL-6 concentration with the Kruskal–Wallis test. We chose a significance level of P < 0.05. We made no adjustments for multiple comparisons, because such adjustments have been suggested to lead to more errors in interpretation(30). The results should therefore be interpreted in the context of the number of tests conducted.
We carried out a haplotype analysis with the 2 IL6 SNPs that have previously been associated with IL-6 concentration [–174G>C(10)(11)(12)(14) and –572G>C(15)] and the functional –6331T>C SNP, because if the genotypes are truly functional, they should show a greater degree of genetic association with the phenotype. In all 3 of the cohorts examined, haplotype was not strongly associated with IL-6 concentration, although it was evident that the haplotype containing –6331C was consistently associated with lower IL-6 concentrations. This result indicates that the –174G>C and –572G>C variants are of less importance than previously thought(14)(15). A recent study by Malarstig et al. identified an association between –572 and IL-6 concentration (P = 0.07), the significance of which increased when acute coronary syndrome patients with higher IL-6 concentrations were selected (P = 0.01)(15). The –572G>C variant does not alter a transcription factor–binding site in vitro (unpublished data), and it is likely that the SNP is simply in strong LD with a functional SNP. Indeed, there is high LD between –572 and –6331, with –572C always being associated with –6331T.
A recent study by Müller-Steinhardt et al. examined IL6 haplotypes and IL-6 concentrations in healthy blood donors after stimulating their monocytes and T lymphocytes with lipopolysaccharide for 4 h(31). In an analysis of IL6 –572/–174 haplotypes, the investigators observed an association of significantly lower concentrations of IL-6 secreted with the G/G haplotype than with all other haplotypes combined; however, the SNPs used by Müller-Steinhardt et al. do not distinguish the low IL-6–secreting haplotype containing –6331C that we have described.
To assess the in vitro ability of the –6331T>C SNP to confer alterations in IL6 expression patterns, we designed a reporter assay with constructs consisting of a region 100 bp on either side of –6331T>C spliced to the minimal IL6 promoter. The results indicated no difference between T and C constructs in unstimulated cells, but IL-1β stimulation produced a noticeable difference for the 2 constructs containing –6331T, compared with the 2 constructs containing the –6331C variant. Interestingly, there was no significant difference in IL6 expression between the –174 G and C alleles when spliced to this upstream element, suggesting that the previously reported effects of this SNP are diluted when larger constructs are used. The reporter assay has limitations as a vehicle for testing the effects of SNPs on the concentrations of products of expressed genes because of the restriction on insert size, the large numbers of possible SNP combinations, and the role that long-range chromatin structure may play in this situation (data not shown); however, the effects of IL-1β stimulation on the reporter constructs seem to reflect a pattern similar to what is seen in patients with acute inflammation in that the –6331T allele confers higher IL6 expression than the C allele.
The EMSA analysis revealed that transcription factor Oct-1 bound to –6331T to a much greater degree than to –6331C. Oct-1 plays a major role in controlling housekeeping genes and tissue-specific genes, including genes encoding several other cytokines: TNF[tumor necrosis factor (TNF superfamily, member 2)](32), IL2 (interleukin 2)(33), IL3 [interleukin 3 (colony-stimulating factor, multiple)], IL5 [interleukin 5 (colony-stimulating factor, eosinophil)], CSF2 [colony stimulating factor 2 (granulocyte-macrophage)](34), and IL8 (interleukin 8)(35). Oct-1 now appears to play a role in IL6 expression. Binding to the –6331T allele was increased when we used nuclear extract from IL-1β–stimulated HeLa cells compared with the use of extract from unstimulated cells (data not shown), and this finding correlates with both in vitro expression data and data on the effect of acute inflammation on plasma IL-6 concentrations in individuals with the –6331T allele. This result may be explained by several mechanisms. Oct-1 may be up-regulated by inflammatory stimuli, DNA/Oct-1 complexes at the –6331T/C position may require the presence of inflammation-induced cofactors such as AP-1 or nuclear factor
B for increased DNA-binding affinity, or Oct-1 binding requires these or other cofactors to activate the IL6 transcriptional machinery.
Variation in inflammatory responses between individuals can have considerable effects on disease outcomes for a number of conditions, and such variation may play a role in the outcomes of cardiac surgery(36)(37), in which IL-6 concentrations remain high in patients with postoperative complications. Indeed, preoperative IL-6 concentrations can act as a prognostic marker for postoperative complications (data not shown). Lowering the IL-6 concentration with medication may be a useful way to reduce adverse outcomes in such patients, and knowledge of the genetic contribution to IL-6 variation may facilitate this process.
The data we have presented show that the effect of the –6331T>C polymorphism on IL-6 concentration is at the transcriptional level following acute stimulation. Numerous studies have examined the –174G>C variant as an indicator of IL-6 concentration and its association with disease. The rare alleles –6331C and –174C never occur on the same chromosome, and studies that associate –174G>C with a particular phenotype may in fact be using this SNP as a marker for a SNP with greater functionality, such as –6331T>C. Alternatively, both the –174G>C and –6331T>C variants may act in concert to affect gene expression. Fully addressing which SNPs are functional requires a haplotypic analysis of the entire IL6 gene and its regulatory region, and this investigation is currently under way.
| Acknowledgments |
|---|
Financial Disclosures: None declared.
Acknowledgments: We thank Winship Herr for his kind gift of the Oct-1 expression plasmid. The key investigators in the HIFMECH study are Anders Hamsten, Steve E. Humphries, Irène Juhan-Vague, Maurizio Margaglione, Giovanni di Minno, John Yudkin, and Elena Tremoli. We acknowledge all of the staff of the Periodontology Unit at the UCL Eastman Dental Institute for their help and support.
| Footnotes |
|---|
2 Human genes: IL6, interleukin 6 (interferon, beta 2); TNF, tumor necrosis factor (TNF superfamily, member 2); IL2, interleukin 2; IL3, interleukin 3 (colony-stimulating factor, multiple); IL5, interleukin 5 (colony-stimulating factor, eosinophil); CSF2, colony stimulating factor 2 (granulocyte-macrophage); IL8, interleukin 8. ![]()
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
B. P. Fairfax, F. O. Vannberg, J. Radhakrishnan, H. Hakonarson, B. J. Keating, A. V.S. Hill, and J. C. Knight An integrated expression phenotype mapping approach defines common variants in LEP, ALOX15 and CAPNS1 associated with induction of IL-6 Hum. Mol. Genet., February 15, 2010; 19(4): 720 - 730. [Abstract] [Full Text] [PDF] |
||||
![]() |
J-Q He, M G Foreman, K Shumansky, X Zhang, L Akhabir, D D Sin, S F P Man, D L DeMeo, A A Litonjua, E K Silverman, et al. Associations of IL6 polymorphisms with lung function decline and COPD Thorax, August 1, 2009; 64(8): 698 - 704. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |