| 2009 JASN IMPACT FACTOR 7.689 | HOME AUTHOR INFO EDITORIAL BOARD SUBSCRIBE FEEDBACK ALERTS HELP | |||
| CURRENT ISSUE | ARCHIVES | JASN Express | ONLINE SUBMISSION | |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
CLINICAL EPIDEMIOLOGY |







,

,

* Department of Epidemiology and Welch Center for Prevention, Epidemiology & Clinical Research, Johns Hopkins University, Baltimore, Maryland;
Center for Population Studies, National Heart, Lung, and Blood Institute, Bethesda, Maryland, and National Heart, Lung, and Blood Institute's Framingham Heart Study, Framingham, Massachusetts;
Department of Medicine, University of Maryland School of Medicine, Baltimore, Maryland;
Department of Epidemiology and Carolina Center for Genome Sciences, University of North Carolina, Chapel Hill, North Carolina; || Division of Epidemiology and Community Health, University of Minnesota, Minneapolis, Minnesota; ¶ Department of Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts; ** Diabetes Unit and Center for Human Genetic Research, Massachusetts General Hospital, Boston, and Program in Medical and Population Genetics, Broad Institute of Harvard and MIT, Cambridge, Massachusetts; 
Human Genetics Center and Institute of Molecular Medicine, University of Texas Health Science Center, Houston, Texas; 
Geriatrics Research and Education Clinical Center, Veterans Administration Medical Center, Baltimore, Maryland; and 
Division of Endocrinology, Hypertension, and Metabolism, Brigham and Women's Hospital and Harvard Medical School, Boston, Massachusetts
Correspondence: Dr. Caroline S. Fox, National Heart, Lung, and Blood Institute's Framingham Heart Study, 73 Mount Wayte Avenue, Suite 2, Framingham, MA 01702; Phone: 508-935-3447; Fax: 508-626-1262; E-mail: foxca@nhlbi.nih.gov; or Dr. Linda Kao, 615 N. Wolfe Street, Room W6513, Baltimore, MD 21205; Phone: 410-614-0945; Fax: 410-955-0863; E-mail: wkao{at}jhsph.edu
Received for publication December 6, 2007. Accepted for publication April 23, 2008.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Variants in the transcription factor 7–like 2 (TCF7L2) gene, specifically single-nucleotide polymorphism (SNP) rs12255372 and rs7901346, have been identified as conferring type 2 diabetes risk and impaired glucose tolerance across many diverse populations.8–18 The TCF7L2 gene product TCF4 is a transcription factor that acts as a downstream effector in the canonical Wnt signaling pathway and may be important in the development of both type 2 diabetes and renal development and disease.19,20 It is therefore plausible that mutations in this gene could manifest themselves in reduced kidney function or kidney disease via their effects on hyperglycemia as well as independent of this mechanism.
The objective of our study was to examine whether variants in TCF7L2 are associated with CKD progression and measures of kidney function in population-based studies of individuals with and without diabetes. We specifically studied the intronic SNPs rs12255372, rs7903146, rs7901695, and rs11196205 originally reported to be associated with type 2 diabetes.8 Initial analyses were conducted in self-identified white and black participants of the population-based prospective Atherosclerosis Risk in Communities (ARIC) Study.21 Confirmation of results was sought in the Framingham Heart Study Offspring Cohort (FHS)22 and in the cross-sectional Heredity and Phenotype Intervention (HAPI) Heart Study.23
| RESULTS |
|---|
|
|
|---|
|
|
|
Similarly, in black ARIC participants, both rs7903146 and rs7901695 were significantly associated with CKD progression in multivariable adjusted analyses overall. The association remained significant for rs7901695 among individuals without diabetes at baseline. Among those with baseline diabetes, the risk alleles at rs7903146 and rs7901695 were significantly associated with CKD progression. Tests of statistical interaction between SNP and diabetes status were NS (smallest interaction P value 0.14 for rs7903146). Further adjustment for mean percentage of European ancestry did not alter the results observed among black ARIC participants.
The most significant association was observed for rs7901695 (r2 = 0.87 with rs7903146 in ARIC white participants and r2 = 0.49 in black participants; Figure 1) in both white (fully adjusted hazard ratio [HR] 1.19 per C allele; 95% confidence interval [CI] 1.06 to 1.34; P = 0.002) and black ARIC participants (fully adjusted HR 1.27 per C allele; 95% CI 1.09 to 1.48; P = 0.002; Table 2).
Role of Hyperglycemia in the Association of TCF7L2 Variants and CKD Progression in the ARIC Study
For further addressing whether TCF7L2 variants were associated with CKD progression among ARIC participants without baseline diabetes via hyperglycemia, several analyses for rs7903146 and rs7901695 were conducted among ARIC participants without baseline diabetes. First, diabetes and fasting glucose levels during follow-up visits were incorporated as time-dependent variables (Table 3, model 1). For rs7901695, the HR for CKD progression associated with each copy of the C allele remained significant at 1.19 in white participants (P = 0.010) and 1.22 in black participants (P = 0.042), whereas rs7901346 remained significant in white participants only.
|
Last, a very conservative analysis was performed by additionally excluding individuals with impaired fasting glucose (IFG; fasting glucose
100 mg/dl) at baseline, which led to an additional reduction of events by approximately 45%. The HR of CKD progression for both rs7901695 and rs7903146 decreased by approximately 40% in white participants while remaining of similar magnitude in black participants but were no longer statistically significant (Table 3, model 3). Further incorporation of time-dependent variables for hypertension and use of angiotensin-converting enzyme inhibitors did not substantially change these results.
Relation of SNPs in TCF7L2 and eGFR/Albuminuria in the ARIC Study
Supplemental Table 2 shows the association between the most CKD progression-associated TCF7L2 variants rs7903146 and rs7901695 and both continuous eGFR and the dichotomized outcome of eGFR <60 ml/min per 1.73 m2 at ARIC study visit 4, at which both urinary albumin and creatinine were measured. The association between the variants and either phenotype was not statistically significant. Results of further exploration of the association of rs7901695 and eGFR using a nonparametric lowess graph are shown in Supplemental Figure 2. None of the TCF7L2 variant risk alleles was associated with increased urinary albumin-to-creatinine ratio (ACR) as a marker of kidney damage at ARIC visit 4 (Supplemental Table 3).
Association between TCF7L2 SNPs and Kidney Function Is Confirmed in the FHS and HAPI Studies
FHS Offspring Cohort
In FHS, SNPs rs12255372 and rs7901346 were genotyped, which were also available in ARIC. Confirmation was therefore sought of significant associations between rs12255372 and rs7901346 with eGFR and serum cystatin C (cross-sectionally at examination cycle 7), as well as ACR (examination cycle 6). The same risk alleles as observed in ARIC at rs12255372 (T) and rs7903146 (T) were significantly associated with lower eGFR and higher cystatin C, both indicating lower kidney function, in multivariable adjusted analyses (Table 4). Limiting the analyses to participants without diabetes, both SNPs remained significantly associated with cystatin C (P = 0.02 for rs12255372 and 0.01 for rs7903146), but only rs7903146 remained marginally associated with eGFR (P = 0.066). In sensitivity analyses, the additional exclusion of individuals with hemoglobin A1c
6.0% and IFG from these analyses reduced the sample size by approximately 20 and 40%, respectively (Table 4). Whereas the association of both SNPs and cystatin C among those with hemoglobin A1c <6.0% remained significant (P < 0.007), regression coefficients among those with IFG were attenuated by approximately 30% and became NS for both variants (P < 0.18).
|
HAPI Heart Study
We further sought to confirm our findings using data from 861 HAPI participants who had clean genotype data available on 347,144 markers across the genome (Affymetrix 500K assay). Two of the SNPs genotyped in ARIC, rs11196205 and rs7901695, were included on the Affymetrix 500K chip. Moreover, 64 other SNP in the 800-kb region around the TCF7L2 gene were also examined for further mapping of the genetic architecture of TCF7L2 variation and kidney phenotypes (Supplemental Table 4).
Table 5 shows that, consistent with the findings for CKD progression in ARIC, the rs7901695 C allele was significantly associated with lower eGFR (P = 0.049; P = 0.013 after further adjustment for fasting glucose available in only 465 participants). The variant rs11196205 was not significantly associated with eGFR in HAPI (P = 0.39; P = 0.12 after adjustment for fasting glucose). The strongest associations were observed for rs7918405, rs7096151, rs7070850, and rs7068695 (smallest P = 0.0003) all within one LD block approximately 200 kb upstream, and rs1904444 and rs10885436 (smallest P = 0.002) in a different LD block approximately 150 kb downstream of the transcribed gene (Table 5, Supplemental Figure 3). The associations remained significant after adjustment for the 14 LD blocks in the region studied for several markers both upstream and downstream of TCF7L2 (P < 0.0036). Neither of these associated regions was in strong LD with rs7901695.
|
| DISCUSSION |
|---|
|
|
|---|
Only one previous study investigated, as a secondary outcome, the association of rs7903146 and renal function among individuals with diabetes. Individuals homozygous for the risk genotype (TT, n = 12) had significantly worse renal function measured by 24-h creatinine clearance.25 This observation is consistent with our findings in ARIC, where white and black individuals with diabetes had 20 to 30% higher risk for CKD progression with each copy of the risk alleles at rs7901346 and rs7901695. Our study extends current knowledge by describing a significant association between TCF7L2 variants and renal phenotypes not only among individuals with diabetes but also in three population-based samples overall and specifically among those without known diabetes.
Several phenomena could explain the associations we observed: (1) The TCF7L2 variants lead to increased diabetes risk, and subsequent hyperglycemia causes the renal phenotypes. The observed association among individuals without diabetes is the result of residual confounding by undetected diabetes and impaired glucose metabolism. (2) The diabetes-susceptibility variants increase the risk for CKD not only through their effect on diabetes but also through other renal-specific mechanisms (pleiotropic effect). (3) Other variants in the TCF7L2 region in LD with the diabetes-susceptibility variants affect susceptibility to CKD (allelic heterogeneity).
Sensitivity analyses in both ARIC and FHS suggest that the impact of frank diabetes does not fully account for the observed association between TCF7L2 SNPs and renal phenotypes; however, the role of IFG as a potential mediator is less clear. In white ARIC participants without IFG (fasting glucose <100 mg/dl), free of baseline or incident diabetes, and accounting for time-dependent fasting glucose, we observed an approximately 40% attenuation of the multivariable adjusted HR for both rs7903146 and rs7901695. Similarly, sensitivity analyses in FHS excluding individuals with IFG showed an approximately 30% attenuation of the effect size. Although the sample sizes for these sensitivity analyses in both ARIC and FHS were reduced by almost 50%, this suggests that IFG contributes partly to the significant association between TCF7L2 variants and kidney function, even among those without known diabetes or those who will not develop diabetes for years.
Conversely, the effect size for rs7901695 among black ARIC participants remained unchanged in these sensitivity analyses. Also, a significant association of rs7901695 and eGFR was observed in HAPI, where >95% of participants with available measurements had fasting glucose values <100 mg/dl, possibly suggesting an association between rs7901695 and renal phenotypes independent of IFG.
A possible mechanism underlying an association of TCF7L2 variants and kidney function independent of diabetes could be pleiotropy. Pleiotropy has been reported from previous genetic studies of people with diabetes as one of the phenotypes studied.26–28 The TCF7L2 gene product TCF4 is a downstream effector in the canonical Wnt signaling pathway, which is essential in the developing kidney.29 Mutations in components of this pathway have been identified as the cause of inherited kidney disease.19,20 Although our study may suggest a mechanistic link between Wnt signaling and kidney function, the underlying molecular mechanisms clearly need to be further studied. Future genome-wide association studies may help to identify additional important components of this pathway.
Regardless of the causal mechanism, however, the consistently observed significant association between common TCF7L2 variants and kidney disease and markers of kidney function in three population-based studies among individuals without diagnosed diabetes is of importance. If independently associated, carriers of the risk alleles at rs7901695 and rs7903146, accounting for 30 to 45% of the population-based samples examined here, might be more susceptible to CKD progression or reduced kidney function independent of hyperglycemia. If mediated by IFG, then our study draws attention to the measurable, significant effect that glucose levels between 100 and 125 mg/dl might have on renal phenotypes among individuals without diabetes in the general population. This is of particular importance because an association of TCF7L2 variants and albuminuria, which usually brings diabetic kidney damage to clinical attention later during the course of the disease, was not observed in our study. This lack of association between TCF7L2 variants and albuminuria in ARIC and FHS might be explained by different underlying biologic processes of increased albuminuria and reduced eGFR.
The results presented here must be viewed with several limitations in mind. First, in contrast to both FHS and HAPI, we did not observe a significant cross-sectional association between TCF7L2 variants and continuous eGFR in ARIC, which might be explained by the apparent U-shape association between allele frequency and eGFR (Supplemental Figure 2). Further research in additional studies therefore seems warranted. The stronger association of TCF7L2 variants with cystatin C compared with eGFR as markers of kidney function in FHS could be due to better sensitivity of cystatin C to changes in eGFR at higher levels.30
Second, our study could not examine identical variants or identical phenotypes across all three of the studies. Instead, two identical variants were available for both ARIC and FHS with similar, significant associations observed in white ARIC and FHS participants, as well as another two variants in both ARIC and HAPI, one of which was significantly associated with renal phenotypes in both studies. It should be noted, however, that rs7903146 and rs7901695 are in strong LD (pair-wise r2 = 0.87 in white ARIC participants), therefore providing indirect evidence for the expected observation with the variant not genotyped. Despite the differences across the studies, the consistency of the direction and magnitude of effects between the studies lends confidence to the results we observed. Although in FHS and HAPI no information was available for a prospective analysis comparable to ARIC, multivariable adjusted analyses accounted for the same covariates, and the rigorous conduct and covariate collection of the FHS and HAPI studies lends additional confidence to the results. The 95% CI of the odds ratios for eGFR <60 ml/min per 1.73 m2 in FHS included the estimates of the HR in ARIC. In fact, the point estimates were almost identical for the two markers typed in both studies, rs12255372 and rs7903146, although the results were not statistically significant in FHS as a result of the smaller sample size. In HAPI, rs7901695 was significantly associated with continuous eGFR, whereas rs11196205 was not. A potential explanation is the limited power in the HAPI sample, allowing only replication of the most significantly associated SNP in ARIC.
Third, our study is further limited because we cannot rule out the possibility that one or more variants in LD with rs7901695 and/or rs7903146 are the causal variant affecting kidney function. Sequencing of all TCF7L2 exons in a previous study did not reveal any nonsynonymous substitutions.8 One previous publication described a SNP upstream of TCF7L2, rs10885390, as more strongly associated with diabetes than rs7903146.15 Because only three HAPI participants had diabetes, we could not assess whether the same variants outside the coding region we observed associated with lower eGFR are also associated with diabetes, which could be addressed in future studies.
Another potential limitation to inferences from association studies is confounding by undetected population stratification. This potential threat is reduced in our study, because we observed the association between TCF7L2 variants and renal traits in several independent study samples of European ancestry (ARIC white participants, FHS, and HAPI) as well as among black ARIC participants. Moreover, sensitivity analyses among black ARIC participants adjusting for percentage of European ancestry did not change the results we observed. Finally, in previous diabetes studies across multiple populations of European ancestry, associations with variants at the TCF7L2 locus have not been found to be due to underlying population stratification.12
Last, our study is limited by eGFR's being an imperfect measure of kidney function. Serum creatinine levels are strongly influenced by muscle mass, and the estimation of GFR using serum creatinine is most accurate in the low ranges of eGFR.30 In FHS, however, we observed significant associations between TCF7L2 variants and serum cystatin C, which is not as influenced by muscle mass.
In summary, variants in TCF7L2 were consistently associated with renal traits in three population-based studies. Risk allele carriers overall as well as among individuals without known diabetes are at increased risk for CKD progression and reduced kidney function. Further studies are needed to replicate our findings and to investigate the underlying mechanisms of this association in more detail.
| CONCISE METHODS |
|---|
|
|
|---|
1.8 mg/dl (
159 µmol/L, women) and
2.0 mg/dl (
177 µmol/L, men; n = 39) were excluded from the CKD progression analyses to ensure that the results were not driven by a few individuals with advanced kidney disease.31 Individuals with missing covariates were excluded from all regression analyses for a final study sample of 10,786 white and 3599 black participants in these analyses.
Framingham Heart Offspring Study.
The FHS is an observational community-based cohort of cardiovascular disease and its risk factors that was initiated in 1948 with the recruitment of 5209 women and men into the original cohort. In 1971, 5124 children or spouses of the original cohort were enrolled into the offspring cohort; participants were examined approximately every 4 to 8 yr.32 This investigation is composed of offspring participants who attended the seventh examination cycle (1998 through 2001). Of a total of 3537 participants who attended, 2468 had genotype data available for TCF7L2, including 2359 for rs12255372 and 2368 for rs7903146. Of these, 2376 had serum creatinine measured, 2330 had cystatin C measured, and 2306 (2323) had a complete covariate profile and made up the final sample size for rs12255372 (rs7901346). Of these participants, 1982 also attended the sixth examination cycle (1995 through 1998) and had a measurement of ACR available. The Framingham offspring study protocol is approved by Boston University Medical Center institutional review board, and all participants signed the informed consent.
HAPI Heart Study.
The HAPI Heart Study consists of 868 related and unrelated Old Order Amish people aged 20 to 80 yr from Lancaster, PA, and is described in detail elsewhere.23 The Old Order Amish are a genetically homogeneous closed founder population of European descent. This generally healthy population lives an active lifestyle with substantial environmental homogeneity. Most of the HAPI participants were recruited on the basis of their participation in two previous studies.33,34 Between 2003 and 2006, HAPI participants underwent two clinic examinations 2 wk apart. For this analysis, data were available for 861 HAPI participants. The protocol was approved by the institutional review board of the University of Maryland, and informed content was obtained.
Assessment of Baseline Characteristics
ARIC participants provided detailed information on demographic, socioeconomic, health behavior, risk factor control, and medical history variables as described previously.31 Racial affiliation was self-reported using the terms "black" or "white" used throughout this report. Details regarding the methods of risk factor measurement have been described previously.31,35 Participants were considered as having diabetes at baseline and each follow-up visit when they reported a physician diagnosis of diabetes or current intake of diabetes medication or when they had a fasting glucose of
126 mg/dl (
7 mmol/L) or nonfasting glucose of
200 mg/dl (
11.1 mmol/L).
In FHS, covariate data were obtained at the seventh examination cycle, and standard definitions were applied. Details regarding the methods of risk factor measurement have been previously described.36 Participants were considered as having diabetes when they were treated with diabetes medication or when they had a fasting glucose of
126 mg/dl (
7 mmol/L).
In HAPI, diabetes was defined as a physician-documented diagnosis of diabetes or the current use of antidiabetic medication. Additional clinical measurements were obtained including fasting glucose (available for 466 [54%] of HAPI participants) and BP.
Genotyping
Genotyping of TCF7L2 polymorphisms rs12255372, rs7903146, rs7901695, rs11196205, and rs7895340 in ARIC was performed separately using the TaqMan assay (Applied Biosystems, Foster City, CA) by the ARIC Central DNA Laboratory. The average call rate of the five SNPs was 93.3%. Hardy-Weinberg expectation was met for all five SNPs among white and black individuals overall (Supplemental Table 1). The
statistic for at least 705 replicates per SNP was
0.92 for each SNP. Missing genotype data were nondifferential with respect to the development of CKD progression. Results for SNP rs7895340 were not further investigated in this study because of the strong LD between rs11196205 and rs7895340 (r2 = 0.97 and 0.98 in white and black ARIC participants, respectively).
In FHS, rs12255372 and rs7903146 were genotyped by allele-specific multiplex primer extension of PCR-amplified products with detection by matrix-assisted laser desorption ionization-time of flight mass spectroscopy using the Sequenom iPLEX platform as described previously.37 The call rates for these two SNPs were 98.6 and 98.0%, respectively, and concordance was 100% among 254 individuals genotyped twice. Hardy-Weinberg expectation was evaluated among unrelated individuals (n = 1689).
HAPI participants were genotyped using the Affymetrix GeneChip Human Mapping 500K Array. The GeneChip Genotyping Analysis Software (GTYPE 4.0) was used for automated genotype calling, and the GTYPE-generated chip files were re-analyzed using the BRLMM genotype calling algorithm (http://www.affymetrix.com). A confidence threshold for call quality of 0.5 was used for this analysis. The resulting mean call rate across of the 861 samples was 97.5%. SNPs that were monomorphic, had a minor allele frequency <2%, or deviated from Hardy-Weinberg equilibrium (P < 0.001) were excluded. The average call rate of the 66 SNPs in the TCF7L2 region was 98.7%.
Definition of the Study Outcome
Study outcomes were categorized as primary and secondary to address the testing of multiple, albeit correlated, outcomes. For the prospective analysis of the primary outcome in ARIC, CKD progression over an average of 13.2 yr of follow-up was defined as (1) an increase in serum creatinine levels
0.4 mg/dl above baseline or (2) a hospitalization discharge or death coded for chronic renal disease (International Classification of Disease, Ninth Revision [ICD-9] codes 581 to 583 or 585 to 588), hypertensive renal disease (ICD-9 code 403), hypertensive heart and renal disease (ICD-9 code 404), unspecified disorder of kidney and ureter (ICD-9 code 593.9), diabetes with renal manifestations (ICD-9 code 250.4), kidney transplant, renal dialysis or adjustment/fitting of catheter (ICD-9 codes V42.0, V45.1, or V56), or either hemodialysis (ICD-9-Clinical Modification 39.95) or peritoneal dialysis procedure codes (ICD-9-Clinical Modification 54.98) without acute renal failure (ICD-9 codes 584, 586, 788.9, and 958.5).31 Assessment of short-term variability within ARIC participants revealed that 0.18 mg/dl was the minimal change in creatinine at which 95% confidence existed that a true change had occurred.38 A rise in serum creatinine was defined as a change of at least twice this amount (0.4 mg/dl)
Secondary kidney traits included eGFR at visits 1 and 4 and the urinary ACR at ARIC Study visit 4. Kidney function was assessed by calculation of eGFR using the abbreviated Modification of Diet in Renal Disease (MDRD) Study equation39: eGFR (ml/min per 1.73 m2) = 186.3 * serum creatinine (mg/dl)–1.154 * age–0.203 * 0.742 (if female) * 1.21 (if black). Serum creatinine was measured using a modified kinetic Jaffe reaction, and creatinine values were standardized and calibrated using regression to the Cleveland Clinic laboratory as described previously.40,41 Albuminuria was measured as described previously.42
In FHS, serum creatinine was measured by the modified Jaffe method and statistically calibrated as described previously40,43 before transformation to eGFR via the MDRD equation.24 Stage 3 CKD of higher was defined according to Kidney Disease Outcomes Quality Initiative (K/DOQI) guidelines as an eGFR <60 ml/min per 1.73 m2.24 Cystatin C concentrations were measured on previously frozen samples (Dade Behring Diagnostic, Marburg Germany); the intra- and interassay coefficients of variation were 2.4 and 3.3%, respectively. Urinary albumin concentrations were measured by immunoturbimetry (Tina-quant Albumin assay; Roche Diagnostics, Indianapolis, IN) and indexed to urinary creatinine.
In HAPI, a blood sample was drawn, centrifuged, and the serum was frozen and shipped to Quest Diagnostics (Horsham, PA), where creatinine was measured using the modified kinetic Jaffe method. eGFR was calculated using the abbreviated MDRD Study equation.39 The urinary ACR was not measured in this study.
Statistical Analysis
In ARIC, all analyses were conducted separately by self-reported race. Differences in baseline characteristics by genotype were assessed using
2 tests and ANOVA, as applicable. The primary analysis evaluated associations of TCF7L2 variants with time to CKD progression. Follow-up time was defined as time from age 45 when participants became eligible for ARIC enrollment with staggered entries at the baseline ARIC visit date, to the ARIC visit date at which the increase in serum creatinine occurred, date of CKD hospitalization discharge or death date, or the earlier of the date of last contact or January 1, 2003. Variants in TCF7L2 as independent predictors of CKD incidence were evaluated in Cox proportional hazards models. An additive genetic model specified a priori was used to evaluate the genetic effect, and previously described risk variants for diabetes were coded as the risk alleles.8,14 The proportionality assumption for Cox models of each variant was assessed by inspection of the complementary log(–log[survival function]) curves and was found not to be violated in the analysis of Schoenfeld residuals. The fully multivariable-adjusted Cox model used in analyses included adjustment for baseline CKD risk factors determined a priori (gender, diabetes, fasting glucose, systolic BP, antihypertensive medication intake, body mass index, smoking) and study center in white participants as covariates. Sensitivity analyses were conducted among black ARIC participants: A set of 1536 ancestry informative markers genotyped among 3965 black ARIC participants was used to derive mean percentage of European ancestry using the software Ancestrymap,44 and this variable was included as a covariate in the regression analyses. The role of diabetes as an effect modifier was evaluated by stratifying by baseline diabetes status and evaluating the statistical interaction between SNP and baseline diabetes status. Sensitivity analyses among participants without baseline diabetes consisted of (1) incorporating diabetes and fasting glucose as time-dependent variables into the multivariable adjusted Cox model; (2) treating CKD progression and incident diabetes as competing events, removing individuals from the analysis at the date of incident diabetes when this date was before or the same as the date of CKD progression and accounting for time-dependent fasting glucose; and (3) additionally excluding individuals with baseline fasting glucose
100 mg/dl from the previous model. Time to incident diabetes was defined as described previously.45 A similar analytical approach was used for the secondary kidney phenotypes of eGFR and ACR (general linear regression) and eGFR <60 ml/min per 1.73 m2 (logistic regression). For normalization of the distribution and emphasize differences at lower ranges of eGFR, where the estimation of GFR is more accurate,30 the transformation ln(eGFR) was analyzed. Similarly, ln(ACR) was analyzed to normalize the distribution. Analyses of eGFR as the outcome included adjustment for age and gender as covariates for two reasons: (1) To account for the relationship of these covariates to the exposure as well as to the outcome, and (2) to account for the relation of age and gender to glomerular filtration beyond their relation to serum creatinine concentrations. All analyses were conducted using Stata 9.2 statistical software (College Station, TX). The LD plots in Figure 1 were generated using the program Haploview 4.0.46 Given the a priori study hypothesis, the investigation of candidate variants in TCF7L2,8 and the moderate to high correlation of the SNP (r2 range 0.44 to 0.87), P < 0.05 for the primary outcomes and P < 0.01 for secondary outcomes was considered statistically significant in all studies.
In FHS, analyses were performed using the generalized estimating equation procedure to account for the relatedness of observations among family members. Models were run with the following cross-sectional outcomes: (2) Stage 3 CKD or higher as a dichotomous trait, (2) ln(eGFR), (3) ln(cystatin C), (4) ln(ACR). Primary outcomes were eGFR and cystatin C levels; albuminuria was evaluated as a secondary outcome. Analogous to ARIC, data were analyzed assuming an additive model and adjusting for age, gender, body mass index, systolic BP, hypertension treatment, diabetes, fasting glucose, and smoking. Secondary analyses were performed after exclusion of individuals with diabetes at examination cycle 7. All analyses were conducted using SAS 8.1 (SAS Institute, Cary, NC).
In HAPI, association analyses were performed using variance component analyses to account for relatedness in individuals from large Amish pedigrees using the SOLAR software.47 Continuous eGFR was analyzed as a quantitative trait using a measured genotype approach, and parameters estimates were obtained by maximum likelihood methods. Likelihood ratio tests were used to determine the significance of associations under an additive model. Analyses were performed using the residuals of eGFR adjusted for age, gender, systolic BP, body mass index, and waist circumference, as well as for fasting glucose in participants with available fasting glucose measurements. Although a similar analytical approach was taken, continuous eGFR without logarithmic transformation was the primary outcome because of the healthy nature of this group of individuals, and the multivariate analysis did not account for diabetes because only three individuals had diabetes. Because analyses in HAPI were conducted to confirm an a priori hypothesis, statistical significance for candidate gene variants was defined by P < 0.05 for rs7901695 and rs11196205 in HAPI participants. Because the remaining 64 variants are not independent of each other, a Bonferroni correction was applied for the 14 LD blocks present in the HAPI data (statistical significance at P < 0.0036), even though several of these blocks contain a marker or a proxy of a marker previously described as associated with diabetes.
| DISCLOSURES |
|---|
|
|
|---|
| Acknowledgments |
|---|
We are indebted to the staff and participants in the ARIC Study, the FHS, and the HAPI Heart Study for important contributions.
| Footnotes |
|---|
S.-J.H. and E.R. contributed equally to this work
C.S.F. and W.H.L.K. share senior authorship.
See related editorial, "Appraising the Value of Genomic Association Studies," on pages 1840–1842.
Supplemental information for this article is available online at http://www.jasn.org/.
| REFERENCES |
|---|
|
|
|---|
Related Articles
This article has been cited by other articles:
![]() |
F. Pesce and F. P. Schena Genome-wide association studies in kidney diseases: Quo Vadis? Nephrol. Dial. Transplant., December 1, 2009; 24(12): 3589 - 3592. [Full Text] [PDF] |
||||
![]() |
A. L. George Jr. Appraising the Value of Genomic Association Studies J. Am. Soc. Nephrol., October 1, 2008; 19(10): 1840 - 1842. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
HOME
CURRENT ISSUE
ARCHIVES
JASN Express
ONLINE SUBMISSION
AUTHOR INFO
EDITORIAL BOARD SUBSCRIBE FEEDBACK ALERTS HELP |