Functional Variants in the Lymphotoxin- Gene Predict Cardiovascular Disease in Dialysis Patients
Yongmei Liu*,
Yvette Berthier-Schaad,,
Laura Plantinga,
Nancy E. Fink,
Russell P. Tracy,
Wen Hong Kao,
Michael J. Klag,
Michael W. Smith,|| and
Josef Coresh
* Wake Forest University School of Medicine, Winston-Salem, North Carolina; Johns Hopkins Medical Institutions, Baltimore, Maryland; Laboratory of Genomic Diversity; || Basic Research Program, SAIC-Frederick, National Cancer Institute, Frederick, Maryland; and University of Vermont, Burlington, Vermont
Address correspondence to: Dr. Josef Coresh, Johns Hopkins University, 2024 East Monument Street, Suite 2-600, Baltimore, MD 21205. Phone: 410-955-0495; Fax: 410-955-0476; E-mail: coresh{at}jhu.edu
Received for publication April 5, 2006.
Accepted for publication August 7, 2006.
TNF- that is encoded by lymphotoxin- gene (LTA) regulates adhesionmolecules and IL-6. Previously, a genome-wide case-control studyshowed that LTA gene variants predisposed to cardiovasculardisease (CVD). In a prospective study of 775 dialysis patients,LTA and IL-6 gene variants were tested as independent predictorsof CVD risk. Four polymorphisms in the LTA gene and one in theIL-6 gene were genotyped. CVD events were ascertained from medicalrecords. During a mean follow-up of 2.6 yr, 294 first-incidentCVD events occurred. The LTA 26Asn variant predicted higheradjusted CVD risk (hazard ratio HR 1.33 for each additionalcopy of Asn allele; 95% confidence interval 1.14 to 1.55; P= 0.0003). Two other nonsynonymous polymorphisms in the LTA,13Agr and 51Pro, were associated with lower inflammatory activityand CVD risk. LTA haplotypes (based on all four single-nucleotidepolymorphisms) were associated with inflammatory markers andpredicted CVD risk (P = 0.005) after adjustment. These LTA genotypeassociations were independent of the IL-6 174G/C genotypeassociation that was reported recently. LTA and IL-6 gene variantsindependently predicted risk for CVD among dialysis patients,suggesting that susceptibility in multiple inflammatory pathwayscontribute to the development of CVD.
Serum measures of the inflammatory process have been associatedwith increased cardiovascular disease (CVD) risk in many studies,but it is unclear whether the systemic inflammatory responseis causative or just a marker of the local vascular inflammatoryresponse that is seen in all phases of atherosclerosis. Proinflammatorycytokines, such as TNF and IL-6, have been implicated as importantin atherogenesis and thrombogenesis. We showed previously thata common IL-6 gene polymorphism at position 174 (G/C)in the promoter region (on chromosome 7) predicts incident CVDin incident dialysis patients (1). TNF- is a primary proinflammatorycytokine that stimulates the expression of the "messenger" cytokine,IL-6 (2). Lymphotoxin- (LTA) gene, encoding TNF-, is locatedon chromosome 6p21. The joint effect of LTA and IL-6 gene variationon CVD risk, if any, is unclear.
A recent experimental study reported that LTA gene knockoutmice showed a 62% reduction in fatty streak lesion size (3).A genome-wide screen of 92,788 markers in a Japanese case-controlstudy recently identified the LTA gene as a candidate locusassociated with susceptibility to myocardial infarction (MI)(4). In following up on significant findings, 10 common single-nucleotidepolymorphisms (SNP) of the LTA gene were studied, and 252A/Gin intron 1 was found to be almost completely linked to a missensemutant (Thr26Asn) in exon 3 and associated with MI; subsequentfunctional assays showed a higher transcriptional activity associatedwith the 252G allele and higher mRNA expression of vascularcell-adhesion molecule 1 and E-selectin associated with the26Asn allele. The association of the 252G allele or 26Asn allelewith coronary artery disease or MI was confirmed in anotherJapanese case-control study (5) and a trio family study of whiteindividuals (6) but not in three other case-control studiesin white individuals (79). Because most of those studieswere retrospective, the selection of case patients and controlsubjects could bias the association, and prospective studiesare warranted. Moreover, two nonsynonymous polymorphisms inthe coding region, Cys13Arg and His51Pro, have not yet beenexamined in an association study. Examination of LTA in dialysispatients is helpful in assessing genetic effects in a proinflammatorysetting.
In a prospective study, we investigated the association of LTAvariants with levels of inflammatory markers at baseline andwith incident CVD and cause-specific mortality during 2.6 yrof follow-up in a cohort of multiethnic dialysis patients. Wefurther examined whether the LTA variants predicted risk forCVD independent of IL-6 174G/C variant. Confirmationof the predisposition to CVD associated with the LTA gene ina prospective study will provide a better understanding of therole of TNF- as a key mediator in atherogenesis and/or thrombogenesis.
Study Design and Population
The Choices for Healthy Outcomes in Caring for ESRD (CHOICE)study is a national, prospective cohort study of 1041 incidentdialysis patients who were aged 19 to 95 yr and recruited from81 dialysis clinics associated with Dialysis Clinic, Inc. (DCI,Nashville, TN), New Haven CAPD (New Haven, CT), or Saint Raphaelshospital (New Haven, CT). The study design and enrollment criteriahave been described elsewhere (10). Participants were enrolledfrom October 1995 to June 1998, a median of 45 d after initiationof dialysis (95% within 3.5 mo), and were followed up throughNovember 2000. Only those in the DCI facilities (n = 923) wereable to participate in the specimen bank (n = 898). Genotypeinformation on IL-6 and LTA was available on 775 of 898 participantsfor whom blood was drawn before a dialysis session at a medianof 5.0 mo from the initiation of dialysis (95% within 8.7 mo).Their baseline demographic and clinical data were obtained fromquestionnaires as well as hospital and clinic records. The levelof comorbidity was assessed by a trained nurse on the basisof medical records and clinic staff reports using the Indexof Co-Existent Disease, a standardized and validated four-levelscale that has been tested in multiple studies (11). The studywas approved by the Institutional Review Board at the JohnsHopkins Medical Institutions and the dialysis clinics, and participantsprovided written informed consent.
DNA Analysis
TaqMan assays were used for genotyping the four SNP252A/Gin intron 1 (rs909253), Cys13Arg in exon 2 (rs2229094), His51Proin exon 3 (rs2229092), and Thr26Asn (rs1041981) in exon 3 (atamino acid position 60)of the LTA gene. A length-modifiedsingle base extension protocol (12) was used for genotypingthe 174G/C polymorphism in the IL-6 gene. The statistics,based on 45 blindly split samples from the CHOICE cohort, were100% for the LTA SNP and 96% for the IL-6 174G/C polymorphism.The call rates for the five SNP were 99% or above.
Biochemical Measurements
IL-6 was measured in serum by an ultrasensitive ELISA method(R&D Systems, Minneapolis, MN) with a coefficient of variation(CV; interassay CV) of 7%. Serum high-sensitivity C-reactiveprotein (hs-CRP) was measured using a colorimetric competitiveELISA (interassay CV 8.9%). Plasma fibrinogen was measured usingan automated clot-rate assay (interassay CV 2.9%). Serum P-selectinwas measured using an ultrasensitive, solid-phase sandwich ELISA(interassay CV 8.9%).
Outcome Ascertainment
Death events were ascertained through the dialysis clinics andCenter for Medicare & Medicaid Services (CMS) data. CVDevents included any of the following events that occurred duringthe follow-up period: MI, cerebrovascular accident, coronaryartery bypass graft, percutaneous coronary angioplasty, peripheralartery bypass, amputation, abdominal aortic aneurysm repair,carotid endarterectomy, and sudden coronary death. Medical recordsfrom hospitalizations were requested and, when available, reviewedand adjudicated by two members of the studys outcomescommittee using uniformly applied criteria that were modifiedfrom the Cardiovascular Health Study (13) and HEMO study (11).The statistic for the event adjudication was 95%.
Of 303 deaths, 141 (69% of 204 in-hospital deaths) terminalhospitalization medical records were available. For patientswhose terminal medical records were unavailable, the underlyingcause of death was determined using the Center for Medicare& Medicaid Services death notification form (Form 2746),completed at the time of death by the dialysis unit staff. Theunderlying cause of death was considered to be CVD when thedeath was adjudicated as being due to coronary heart disease,cerebrovascular accident, or peripheral vascular disease.
Statistical Analyses
All regression analyses were performed using STATA 7.0 statisticalsoftware (StataCorp, College Station, TX). P 0.05 was markedas statistically significant and an analysis of the false-positivereporting probability (14) is discussed. The distributions ofhs-CRP, IL-6, and P-selectin levels were highly skewed to theright, and logarithmic transformations were applied for normalization.Regression coefficients from a linear regression model of thelog-transformed levels were converted to percentage change inthe geometric mean by calculating 1 exp(coefficient).P value for trend was estimated by assigning genotypes withthe numeric values 1, 2, and 3.
Survival analysis was conducted from initiation of dialysisto the first CVD event or death. Individuals were censored asa result of renal transplant (n = 129), loss to follow-up (n= 3), or death attributed to causes other than CVD (n = 113).Kaplan-Meier analysis and the log-rank test were used to analyzethe associations adjusted for race. Several groups of covariateswere selected a priori for multivariate Cox proportional hazardsregression models. The first group included the following variables:Age, race, gender, dialysis modality (hemodialysis versus peritonealdialysis), cigarette smoking, dialysis clinics, diabetes, andcomorbidity score. The second group added seven other risk factors:The 174G/C polymorphism of the IL-6 gene, prevalent CVD,congestive heart failure, body mass index (BMI), systolic BP(SBP), serum cholesterol, and albumin. The third group addedthe four inflammatory makers: hs-CRP, IL-6, fibrinogen, andP-selectin. For all survival analyses, the proportionality assumptionof the Cox model was confirmed by inspection of log (log[survivalfunction]) curves and Schoenfeld residuals.
Haplotypes were constructed separately for black and white patientsusing the PHASE program, version 2.0.2 (15). Individuals wereassigned the most likely pair of haplotypes (when the probabilityof assignment was >90%). Haplotype dose was modeled 0, 1,and 2 for LTA. In the haplotype analyses, to use the commonhaplotype as a reference and adjust for other haplotypes, allof the haplotypes except for the common haplotype were includedin the regression models. To minimize the type 1 error rateas a result of multiple comparisons procedure, we tested overallsignificance at the gene using an omnibus likelihood ratio test.
To assess for potential population stratification, we genotypeda panel of 87 ancestry-informative SNP to measure admixture.The degree of individual genetic white-to-black admixture wasestimated using Bayesian methods implemented in the STRUCTUREprogram (version 2) (16). We also tested the null hypothesisthat the LTA gene effect is similar across race by stratifyingthe analyses according to race and evaluating the gene-by-raceinteraction terms. In addition, we explored the potential geneenvironmentand genegene interactions as part of a full presentationof the data.
Table 1 shows the characteristics of the 775 individuals byrace. Black patients tended to be younger and were more likelyto be female, current smokers, and on hemodialysis. Black patientsalso less frequently presented with history of CVD and congestiveheart failure and had higher BMI and SBP and lower levels ofIL-6, CRP, fibrinogen, and P-selectin.
Table 1. Patient characteristics by race (n = 775)a
The genotype frequencies of the 174G/C polymorphismsand the four LTA SNP were consistent with Hardy-Weinberg expectations.Allele frequencies in both subgroups were similar to those reportedpreviously (17). The variants 174C in the IL-6 gene and51Pro in the LTA gene were less frequent in black than in whitepatients, but the two variants 252G and 26Asn in LTA were morecommon in black patients (Table 1). Because the two polymorphisms252A/G and Thr26Asn were in almost total linkage disequilibrium(LD; r2 = 1), data are not presented further for 252A/G. Thepair-wise LD, r2, for the remaining three SNP ranged from 0.02to 0.17 in black patients and from 0.04 to 0.20 in white patients.A four-SNP haplotype had the following common variants: Hap1as carriers of all conserved alleles, Hap2 as carriers of 13Argbut not 51Pro, Hap3 as carriers of both 13Arg and 51Pro, andHap4 as carriers of 26Asn. Hap1 and Hap3 were less common andHap4 was more common in black than in white patients.
LTA Polymorphisms and Levels of Inflammatory Markers
The two novel nonsynonymous polymorphisms, 13Arg and 51Pro,were associated with lower IL-6, CRP, and fibrinogen levels(Figure 1). The 26Asn allele was marginally associated withhigher P-selectin levels in an additive mode of inheritance,and this association was limited to white patients (P = 0.047for trend).
Figure 1. Adjusted levels of inflammatory markers in 775 dialysis patients according to lymphotoxin- gene (LTA) genotypes (after adjustment for age, gender, race, dialysis modality, cigarette smoking, dialysis clinics, diabetes, and comorbidity score).
LTA Polymorphisms and the Risk for CVD and CVD Mortality
During a mean of 2.6 yr of follow-up in 775 individuals, 294individuals developed a new event, which comprised 155 individualswith coronary heart disease events, 52 with cerebral vasculardisease events, and 87 with peripheral vascular disease events.During this period of follow-up, 303 deaths occurred, 152 ofwhich were attributed to CVD. The Cys13Arg genotype was notstatistically associated with risk for CVD (Figure 2A), butthose who carried 51Pro allele (Figure 2B) had lower risk forCVD. Compared with the genotype 26Thr/Thr, 26Thr/Asn and 26Asn/Asnwere associated with higher risk for CVD (Figure 2C). Afteradjustment for age, race, gender, dialysis modality (hemodialysisversus peritoneal dialysis), cigarette smoking, dialysis clinics,diabetes, and comorbidity score, these associations remainedsignificant (Figure 2D). Compared with 26Thr homozygotes, thehazard ratios (HR) for CVD were 1.42 (95% confidence interval[CI] 1.09 to 1.85) for 26Thr/Asn heterozygotes and 1.75 (95%CI 1.27 to 2.40) for 26Asn/Asn homozygotes. Under an additivemode of inheritance, the HR associated with the 26Asn allelewas 1.33 (95% CI 1.14 to 1.55; P = 0.0003) for risk for incidentCVD overall. The similar pattern was seen for CVD mortality(HR 1.31; P = 0.006). This association was not apparent fornon-CVD mortality and was weaker for all-cause mortality (datanot shown). Table 2 shows the results stratified by race; the26Asn variant association with CVD/CVD mortality remained significantin white but not in black patients. This may be due to lowerpower in black patients; interaction with race was not statisticallysignificant.
Figure 2. Cumulative incidence of cardiovascular disease (CVD; adjusted to white patients) according to LTA genotypes Cys13Arg (A), His51Pro (B), Thr26Asn (C) and adjusted hazard ratio (HR; adjusted for age, gender, race, dialysis modality, cigarette smoking, dialysis clinics, diabetes, and comorbidity score) of CVD risk by LTA genotypes (D) in 775 dialysis patients.
Table 2. Adjusted HR of incident CVD and CVD mortality associated with the Thr26Asn in the LTA genea
LTA Thr26Asn, IL-6 174G/C, and the Risk for CVD and CVD Mortality
The previously seen associations of the IL-6 174C allelewith higher CVD risk remained after adjustment for the Thr26Asnpolymorphism of the LTA gene (HR 1.36; P = 0.01 for 174Ccarriers; HR 1.31; P = 0.0002 per one-copy increment of 26Asn).The combined effect of the two at-risk SNP, Thr26Asn and 174G/C,is displayed in Figure 3. The effects of LTA 26Asn allele onboth incident CVD and CVD mortality are consistent across theindividuals with the genotype 174G/G and the 174Ccarriers. The greatest joint effect, compared with those whowere homozygotes for both LTA 26Thr/Thr and IL-6 G/G, was seenfor LTA 26Asn/Asn and IL-6 174C carriers on incidentCVD (HR 2.17; 95% CI 1.41 to 3.89; P = 0.001) and CVD mortality(HR 2.75; 95% CI 1.38 to 6.06; P = 0.005). Carriers of eitherthe LTA 26Asn or IL-6 174C allele were at higher riskfor CVD (HR 1.75 [95% CI 1.16 to 2.63; P = 0.008] overall; HR2.15 [95% CI 1.33 to 3.47; P = 0.002] for white patients; HR1.09 [95% CI 0.59 to 1.99; P = 0.79] for black patients) andhigher risk for CVD mortality (HR 2.62 [95% CI 1.49 to 4.62;P = 0.001] overall; HR 2.30 [95% CI 1.22 to 4.31; P = 0.01]for white patients; HR 6.69 [95% CI 1.01 to 44.43; P = 0.049]for black patients), compared with noncarriers at both loci.There was no statistical interaction between the two gene variants.
Figure 3. HR of CVD and CVD mortality in 757 dialysis patients according to the LTA (Thr25Asn) and IL-6 (174G/C) genotypes (adjusted for age, gender, race, cigarette smoking, dialysis modality, dialysis clinics, history of diabetes, and comorbidity score). Genotypes combined multiplicatively with no interaction, P values are for the combined effect of both genotypes: *P < 0.05; **P < 0.01.
LTA Haplotypes, Levels of Inflammatory Markers, Incident CVD, and CVD Mortality
Haplotype analyses revealed patterns that were similar to thesingle-locus analyses (Table 3). Haplotypes other than the fourpreviously mentioned haplotypes were omitted from the analysesbecause of their small size, and adjusting for them as one groupdid not change the results. Compared with Hap1 homozygotes,each additional copy of Hap2 was associated with 14% lower CRPlevels (P = 0.037) and 19 mg/dl lower fibrinogen (P = 0.029)levels; Hap3 was associated with 24% lower IL-6 (P = 0.029),24% lower CRP (P = 0.07), and 45 mg/dl lower fibrinogen (P =0.005) levels; and Hap4 was not associated with inflammatorymarkers, except for a marginal association with higher serumP-selectin levels (P = 0.08). The four haplotypes explained1 to 2% of the variation in IL-6, CRP, fibrinogen, and P-selectinlevels overall. In white patients, all of the observed associationsin overall analyses remained significant; in addition, Hap4was significantly associated with higher P-selectin levels (P= 0.023). In black patients, only the relationship between Hap3and lower IL-6 levels was significant (P = 0.008).
Table 3. Association of LTA Haplotypes with levels of inflammatory markers, incident CVD, and CVD mortality in 756 dialysis patientsa
Compared with the homozygous Hap1 haplotype, Hap3 was marginallyassociated with lower (P = 0.081), whereas Hap4 was associatedwith higher (P = 0.029), risk for CVD (Table 3). A similar patternwas seen for CVD mortality (Table 3). These estimates remainedafter additional adjustment for the 174G/C polymorphismof the IL-6 gene, prevalent CVD, congestive heart failure, BMI,SBP, serum cholesterol, and albumin. Risk estimates were slightlyattenuated and were marginally significant after additionallycontrolling for any one or all four of the inflammatory markers(data not shown).
Secondary Analyses
We also checked for potential geneenvironment interactionby running stratified analyses. Similar results were seen inthe various subgroups stratified by age (<60 versus60 yr),gender, smoking status, and presence of diabetes or CVD. Inaddition, all regression models were rerun with adjustment ofgenetic admixture. No significant changes were observed withthis additional adjustment (data not shown).
On the basis of DNA variations in the LTA gene (23 polymorphisms)that were retrieved from the resequencing data of 24 AfricanAmerican and 23 European individuals, (17) a single-haplotypeblock can be defined using the CI method (18) implemented inHaploview 3.2. Five more polymorphisms were selected to capturecommon variation in the LTA gene, on the basis of two primarycriteria: (1) The minor allele frequency of an SNP >10% and(2) the minimum level of association between assayed and unassayedSNP, measured by the LD statistic r2 of 0.65 (19). These polymorphismsincluded polymorphisms at sites 306 (rs2009658), 2490 (rs746868),2756 (rs3053942), 3842 (rs3093544), and 4371 (rs1799964). Althoughthese polymorphisms added information on haplotype heterozygosity,the additional haplotypes did not increase the predictive valueof LTA variants for CVD risk.
Ours is the first prospective study to replicate the effectof the haplotype 252G in intron 1 and 26Asn (Hap4) on increasedCVD risk reported by the large-scale, genome-wide Japanese study(4). In our study, the predictive value of the variants forCVD risk was independent of the IL-6 174G/C polymorphismand traditional cardiovascular risk factors. The risk was explainedpartially by systemic inflammatory markers, such as IL-6, CRP,fibrinogen, and P-selectin. Our results also are consistentwith another Japanese case-control study (5) and a trio familystudy of 460 families recruited from the United Kingdom, Sweden,Germany, and Italy (6) but not with three case-control studiesin white individual, in which no association was found (79).This difference may be related to the smaller sample size (148to 199 cases) and exclusion of fatal cases in the case-controlstudies.
Several studies have suggested that the haplotype 252G in intron1 and 26Asn (Hap4) influence gene function. The 252G alleleis associated with higher TNF- levels in human, and greaterproduction of TNF- and TNF- in in vitro studies, compared withthe 252A allele, presumably through higher transcriptional regulatoryactivity (4,2024). Ozaki et al. (4) further showed thatbinding of an unknown nuclear factor to the 252G allele is tighterthan to the 252A allele, pointing to a novel regulatory sequenceelement in intron 1. The 26Asn allele is almost in total LDwith the 252G allele, which makes it almost impossible to studyits individual independent effect in an epidemiologic study.In an in vitro setting, however, biologic activity of a threonine-to-asparaginesubstitution at codon 26 was suggested; the 26Asn allele wasassociated with a higher expression of mRNA coding for vascularcell-adhesion molecule 1 and E-selectin (4). These proteinshave a major role in leukocyte rolling and attachment in atherogenesis(25). Expression of P-selectin, an adhesion receptor that isexpressed on activated endothelial cells and a member of theadhesion protein family, also is regulated by TNF. The higherlevel of P-selectin associated with 26Asn allele that was observedin our study suggests that the high-risk haplotype may be atherogenicpartially through the adhesion protein pathway.
We also report associations of the two nonsynonymous polymorphismsof LTA, Cys13Arg and His51Pro, with circulating inflammatorymarkers and CVD risk. 13Arg and 51Pro were associated with lowerinflammatory states (as assessed by lower levels of IL-6, CRP,and fibrinogen) and lower CVD risk and CVD mortality, whichwas particularly prominent for carriers of both 13Arg and 51Pro(defining Hap3). Given that these variants are in coding regions,our findings suggest that they may alter TNF-s structuralstability and decrease inflammation. However, this first reportof the possible risk implications of these polymorphisms needsto be replicated.
This study has several limitations. We measured products ofactivated TNF-, IL-6, CRP, fibrinogen, and P-selectin levelson only one occasion; the surrogate nature of the measurementand unmeasured variability within individuals may underestimatethe relationship between the LTA variants and inflammatory activities.Given the long duration of the processes that lead up to ESRD,people with variant alleles may have died of CVD before developingESRD. Such selective mortality can lead to an underestimationof the risk that is associated with variant alleles. Becauseterminal medical records for approximately half of the patientswere unavailable for detailed review of their underlying causeof death, nondifferential misclassification bias might occur,which could have underestimated the risk for CVD death or non-CVDdeath associated with variant alleles. In addition, this genotypingstudy did not include all of the study participants as a resultof missing genotype information. The baseline characteristicsof patients who were not included in the genotyping study weresimilar to those who were included in terms of age, gender,race, smoking, BMI, BP, serum total cholesterol, history ofdiabetes, and CVD event rate (all comparisons, P > 0.38).Furthermore, we have limited power to examine racial difference,so the ability to generalize these findings to black individualsis uncertain. Sensitivity analysis that adjusted for percentageof admixture among black patients and was estimated using 87markers did not change the main results. Finally, our studyhas a limited power to detect a small effect for rare polymorphisms.However, our overall sample would have at least 80% power todetect an HR of 1.29 in CVD risk for the variants with a minorallele frequency >25% at a two-tailed of 0.05.
The observed association could be due to a highly linked diseasevariant in or near LTA. Some resequencing efforts have concludedthat three polymorphisms, 10G/A, 252A/G, and Thr26Asn, occuron the same haplotype in general and that the 10G/A did notexpress a biologic activity (4,17,26). Furthermore, Ozaki etal. (4) also screened the flanking region of the LTA gene andidentified two more variants, 63T/A in the promoter regionof NFKBIL1 and 23G/C in the promoter region of BAT1,that were in tight association with the LTA 252A/G variant withinone extended block of intense LD (also known as a haplotypeblock). The two variants in the NFKBIL1 (encoding nuclear factorof k light polypeptide gene enhancer in B cells, inhibitor-like1) and BAT1 gene (encoding HLA-B associated transcript 1), however,did not show biologic activity that would explain their associationwith higher risk for MI. Furthermore, TNF- and LTA loci arelocated in the MHC. A recent sequence analysis of 69 samples(27) found no association between TNF- or LTA alleles and specifichuman leukocyte antigen (HLA) genes (HLA-B and HLA-DR) alleles.Although TNF- alleles had some linkage with LTA alleles, itwas not genotyped because no observed polymorphism in the formerlocus was significantly associated with MI in the Japanese genome-wideassociation study.
Unlike the measurements of serum levels of IL-6 and TNF thatare strongly correlated, the variants in the LTA gene (252A/Gand Thr26Asn) and the variant in the IL-6 gene (174G/C)are uncorrelated because they reside on different chromosomesand thereby provide an option to examine their independent effects.The multiplicative effect of the two gene variants suggeststhat TNF- has additional proatherogenic effects (perhaps mediatedthrough adhesion molecules) in addition to its effect on IL-6.Their independent predictive values for clinically evident CVDand CVD mortality in our study, together with the reported biologicfunctions of these variants in vitro and in vivo, support theimportance of cytokines in determining CVD risk. Therefore,combined evaluation of both genes ought to be more useful inidentifying CVD high-risk groups than screening for either alone.
Interpretation of genetic association studies in the face ofmultiple comparisons can be aided by an analysis of the false-positivereport probability suggested by Wacholder et al. (14). If weassume a most likely relative hazard (substituted for odds ratio)worthy of clinical reporting to be 1.5 for CVD in dialysis patients,then we can estimate the false-positive report probability asa function of the previous probability that a given SNP is associatedwith CVD. For the LTA 26Asn variant (P = 0.0003), the associationreported here is likely to have a false-positive report probabilitiesare 0.002 if a high (0.1) previous probability of associationwith CVD is assumed. Assigning a high previous probability issupported by previous human, animal, and expression studiesof LTA 26Asn (3,4,2024). The false-positive report probabilitiesare 0.001, 0.03, and 0.22 if previous probabilities of 0.25(very high), 0.01 (medium), and 0.001 (low) are assumed. Forthe omnibus test of LTA haplotypes with CVD (P = 0.0053), thefalse-positive report probability is 0.02, 0.05, and 0.35 forprevious probabilities of 0.25 (very high), 0.1 (high), and0.01 (medium). The associations with intermediate phenotypesare less statistically significant and may reflect a higherfalse-positive report rate and therefore should be interpretedmore cautiously. We explored these pathways as the most likelyintermediates and therefore think that finding associationsin the expected directions is supportive of inflammationsplaying a role in the risk that is associated with LTA, butconfirmation of the specific associations of these and otherintermediate markers is needed.
This study was limited to dialysis patients, who experienceboth a high level of inflammation and high risk for CVD andmortality and may have different pathophysiology of CVD fromthe general population. The findings suggest that genetic susceptibilityplays an important role in the response of different patientsto the proinflammatory and proatherogenic state that is inducedby kidney failure and dialysis treatment. More prospective studiesand clinical trials in both the high-inflammation populationsand other populations are needed to understand fully the interplaybetween genetic susceptibility at LTA and IL-6 and CVD risk.Confirmation of the independent effect of this genetic susceptibilitywill strengthen the understanding of the multiple causal pathwaysthat lead to CVD risk and could prompt strategies for targetingeffectors of inflammation directly in prevention and therapyof CVD.
Acknowledgments
CHOICE is supported by RO1-HL-62985 (National Heart, Lung, andBlood Institute), RO1-DK-59616 (National Institute of Diabetesand Digestive and Kidney Diseases), R01-HS-08365 (Agency forHealthcare Research and Quality), and a grant from Baxter HealthcareCorp. J.C. is supported in part as an American Heart Associationestablished investigator (01-4019-7N). M.J.K. is supported byK24-DK-02856 (National Institute of Diabetes and Digestive andKidney Diseases. R.P.T. is supported by HL 46696 and HL 58329.This research was supported by the Intramural Research Programof the National Institutes of Health, National Cancer Institute,Center for Cancer Research. This publication has been fundedin part with federal funds from the National Cancer Institute,National Institutes of Health under contract no. NO1-CO-12400.
We thank the patients, staff, laboratory, and physicians whoparticipated in the CHOICE Study at Dialysis Clinic, Inc., andJohns Hopkins University and the Cardiovascular Endpoint Committee:Current members: Bernard G. Jaar, MD, MPH; Yongmei Liu, MD;Joseph A. Eustace, MD, MHS; Richard M. Ugarte, MD; Melanie H.Katzman, MD, MHS; and J. Craig Longenecker, MD, PhD. Formermembers of the Committee include Michael Klag, MD, MPH; NeilR. Powe, MD, MPH, MBA; Michael J. Choi, MD; Renuka Sothinathan,MD, MHS; and Caroline Fox, MD, MPH. Cardiovascular events adjudicatorsare Nancy E. Fink, MPH, and Laura C. Plantinga, ScM.
Footnotes
Published online ahead of print. Publication date availableat www.jasn.org.
Liu Y, Berthier-Schaad Y, Fallin MD, Fink NE, Tracy RP, Klag MJ, Smith MW, Coresh J: IL-6 haplotypes, inflammation, and risk for cardiovascular disease in a multiethnic dialysis cohort.
J Am Soc Nephrol 17
: 863
870, 2006[Abstract/Free Full Text]
Libby P, Ridker PM: Novel inflammatory markers of coronary risk: Theory versus practice.
Circulation 100
: 1148
1150, 1999[Free Full Text]
Schreyer SA, Vick CM, LeBoeuf RC: Loss of lymphotoxin-alpha but not tumor necrosis factor-alpha reduces atherosclerosis in mice.
J Biol Chem 277
: 12364
12368, 2002[Abstract/Free Full Text]
Ozaki K, Ohnishi Y, Iida A, Sekine A, Yamada R, Tsunoda T, Sato H, Sato H, Hori M, Nakamura Y, Tanaka T: Functional SNPs in the lymphotoxin-alpha gene that are associated with susceptibility to myocardial infarction.
Nat Genet 32
: 650
654, 2002[CrossRef][Medline]
Iwanaga Y, Ono K, Takagi S, Terashima M, Tsutsumi Y, Mannami T, Yasui N, Goto Y, Nonogi H, Iwai N: Association analysis between polymorphisms of the lymphotoxin-alpha gene and myocardial infarction in a Japanese population.
Atherosclerosis 172
: 197
198, 2004[CrossRef][Medline]
A trio family study showing association of the lymphotoxin-alpha N26 (804A) allele with coronary artery disease.
Eur J Hum Genet 12
: 770
774, 2004[CrossRef][Medline]
Padovani JC, Pazin-Filho A, Simoes MV, Marin-Neto JA, Zago MA, Franco RF: Gene polymorphisms in the TNF locus and the risk of myocardial infarction.
Thromb Res 100
: 263
269, 2000[CrossRef][Medline]
Koch W, Kastrati A, Bottiger C, Mehilli J, von Beckerath N, Schomig A: Interleukin-10 and tumor necrosis factor gene polymorphisms and risk of coronary artery disease and myocardial infarction.
Atherosclerosis 159
: 137
144, 2001[CrossRef][Medline]
Longenecker JC, Coresh J, Powe NR, Levey AS, Fink NE, Martin A, Klag MJ: Traditional cardiovascular disease risk factors in dialysis patients compared with the general population: The CHOICE Study.
J Am Soc Nephrol 13
: 1918
1927, 2002[Abstract/Free Full Text]
Rocco MV, Yan G, Gassman J, Lewis JB, Ornt D, Weiss B, Levey AS: Comparison of causes of death using HEMO Study and HCFA end-stage renal disease death notification classification systems. The National Institutes of Health-funded Hemodialysis. Health Care Financing Administration.
Am J Kidney Dis 39
: 146
153, 2002[Medline]
Lindblad-Toh K, Winchester E, Daly MJ, Wang DG, Hirschhorn JN, Laviolette JP, Ardlie K, Reich DE, Robinson E, Sklar P, Shah N, Thomas D, Fan JB, Gingeras T, Warrington J, Patil N, Hudson TJ, Lander ES: Large-scale discovery and genotyping of single-nucleotide polymorphisms in the mouse.
Nat Genet 24
: 381
386, 2000[CrossRef][Medline]
Tracy RP, Lemaitre RN, Psaty BM, Ives DG, Evans RW, Cushman M, Meilahn EN, Kuller LH: Relationship of C-reactive protein to risk of cardiovascular disease in the elderly. Results from the Cardiovascular Health Study and the Rural Health Promotion Project.
Arterioscler Thromb Vasc Biol 17
: 1121
1127, 1997[Abstract/Free Full Text]
Wacholder S, Chanock S, Garcia-Closas M, El GL, Rothman N: Assessing the probability that a positive report is false: An approach for molecular epidemiology studies.
J Natl Cancer Inst 96
: 434
442, 2004[Abstract/Free Full Text]
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data.
Am J Hum Genet 68
: 978
989, 2001[CrossRef][Medline]
Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data.
Genetics 155
: 945
959, 2000[Abstract/Free Full Text]
Carlson CS, Eberle MA, Rieder MJ, Smith JD, Kruglyak L, Nickerson DA: Additional SNPs and linkage-disequilibrium analyses are necessary for whole-genome association studies in humans.
Nat Genet 33
: 518
521, 2003[CrossRef][Medline]
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome.
Science 296
: 2225
2229, 2002[Abstract/Free Full Text]
Carlson CS, Eberle MA, Rieder MJ, Yi Q, Kruglyak L, Nickerson DA: Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium.
Am J Hum Genet 74
: 106
120, 2004[CrossRef][Medline]
Knight JC, Keating BJ, Rockett KA, Kwiatkowski DP: In vivo characterization of regulatory polymorphisms by allele-specific quantification of RNA polymerase loading.
Nat Genet 33
: 469
475, 2003[CrossRef][Medline]
Messer G, Spengler U, Jung MC, Honold G, Blomer K, Pape GR, Riethmuller G, Weiss EH: Polymorphic structure of the tumor necrosis factor (TNF) locus: An NcoI polymorphism in the first intron of the human TNF-beta gene correlates with a variant amino acid in position 26 and a reduced level of TNF-beta production.
J Exp Med 173
: 209
219, 1991[Abstract/Free Full Text]
Lee SH, Park SH, Min JK, Kim SI, Yoo WH, Hong YS, Park JH, Cho CS, Kim TG, Han H, Kim HY: Decreased tumour necrosis factor-beta production in TNFB*2 homozygote: An important predisposing factor of lupus nephritis in Koreans.
Lupus 6
: 603
609, 1997[Abstract/Free Full Text]
Warzocha K, Ribeiro P, Bienvenu J, Roy P, Charlot C, Rigal D, Coiffier B, Salles G: Genetic polymorphisms in the tumor necrosis factor locus influence non-Hodgkins lymphoma outcome.
Blood 91
: 3574
3581, 1998[Abstract/Free Full Text]
Temple SE, Cheong KY, Almeida CM, Price P, Waterer GW: Polymorphisms in lymphotoxin alpha and CD14 genes influence TNFalpha production induced by Gram-positive and Gram-negative bacteria.
Genes Immun 4
: 283
288, 2003[CrossRef][Medline]
Abraham LJ, Du DC, Zahedi K, Dawkins RL, Whitehead AS: Haplotypic polymorphisms of the TNFB gene.
Immunogenetics 33
: 50
53, 1991[Medline]
Posch PE, Cruz I, Bradshaw D, Medhekar BA: Novel polymorphisms and the definition of promoter 'alleles of the tumor necrosis factor and lymphotoxin alpha loci: inclusion in HLA haplotypes.
Genes Immun 4
: 547
558, 2003[CrossRef][Medline]
Related Article
This Months Highlights
J. Am. Soc. Nephrol. 2006 17: 2947-2948.
[Full Text][PDF]