Genomewide Linkage Scan for Diabetic Renal Failure and Albuminuria: The FIND StudyIgo, Jr. R.P.a · Iyengar S.K.a · Nicholas S.B.c · Goddard K.A.B.d · Langefeld C.D.e · Hanson R.L.f · Duggirala R.g · Divers J.e · Abboud H.h · Adler S.G.i · Arar N.H.h · Horvath A.a · Elston R.C.a · Bowden D.W.e · Guo X.j · Ipp E.i · Kao W.H.L.k · Kimmel P.L.l · Knowler W.C.f · Meoni L.A.k · Molineros J.m · Nelson R.G.f · Pahl M.V.n · Parekh R.S.k · Rasooly R.S.l · Schelling J.R.b · Shah V.O.o · Smith M.W.p · Winkler C.A.q · Zager P.G.o · Sedor J.R.b · Freedman B.I.e · on behalf of the Family Investigation of Nephropathy and Diabetes Research Group
aGenetic Analysis and Data Coordinating Center at bCase Western Reserve University, Cleveland, Ohio, cUniversity of California, Los Angeles, Calif., dCenter for Health Research, Kaiser Permanente Northwest, Portland, Ore., eWake Forest School of Medicine, Winston-Salem, N.C., fNational Institute of Diabetes and Digestive and Kidney Diseases (NIDDK), Phoenix, Ariz., gTexas Biomedical Research Institute, and hUniversity of Texas Health Science Center at San Antonio, San Antonio, Tex., iHarbor-University of California Los Angeles Medical Center, Torrance, Calif., jCedars-Sinai Medical Center, Los Angeles, Calif., kJohns Hopkins University, Baltimore, Md., lNIDDK, NIH, Bethesda, Md., mOklahoma Medical Research Foundation, Oklahoma City, Okla., nUniversity of California, Irvine, Calif., oUniversity of New Mexico, Albuquerque,N. Mex., and pGenetics and Genomics Group, Advanced Technology Program, andqBasic Research Laboratory, Center for Cancer Research, National Cancer Institute-Frederick, SAIC-Frederick, Frederick, Md., USA
Background: Diabetic nephropathy (DN) is a leading cause of mortality and morbidity in patients with type 1 and type 2 diabetes. The multicenter FIND consortium aims to identify genes for DN and its associated quantitative traits, e.g. the urine albumin:creatinine ratio (ACR). Herein, the results of whole-genome linkage analysis and a sparse association scan for ACR and a dichotomous DN phenotype are reported in diabetic individuals. Methods: A genomewide scan comprising more than 5,500 autosomal single nucleotide polymorphism markers (average spacing of 0.6 cM) was performed on 1,235 nuclear and extended pedigrees (3,972 diabetic participants) ascertained for DN from African-American (AA), American-Indian (AI), European-American (EA) and Mexican-American (MA) populations. Results: Strong evidence for linkage to DN was detected on chromosome 6p (p = 8.0 × 10–5, LOD = 3.09) in EA families as well as suggestive evidence for linkage to chromosome 7p in AI families. Regions on chromosomes 3p in AA, 7q in EA, 16q in AA and 22q in MA displayed suggestive evidence of linkage for urine ACR. The linkage peak on chromosome 22q overlaps the MYH9/APOL1 gene region, previously implicated in AA diabetic and nondiabetic nephropathies. Conclusion: These results strengthen the evidence for previously identified genomic regions and implicate several novel loci potentially involved in the pathogenesis of DN.
Copyright © 2011 S. Karger AG, Basel
Diabetic nephropathy (DN) associated with type 1 and type 2 diabetes mellitus (DM) remains the leading cause of chronic kidney disease in the US, contributing to approximately 50% of incident cases of end-stage renal disease (ESRD) . DN imposes a significant personal and socioeconomic burden on patients and their families, society, and healthcare systems due, in part, to its contribution to cardiovascular disease [2,3].
Genetic factors contribute to risk of DN in all ethnic groups [4,5]. Genomewide linkage scans of kidney function in the presence of DM implicated regions on chromosomes 3q, 7p, 7q, 9, 10q and 18q (online suppl. table S1; see www.karger.com/doi/10.1159/000326763 for all online suppl. material), supporting the view that genetic factors contribute to DN and albuminuria. To date, four genomewide association analyses have been reported for DN [6,7,8,9]. These and other studies have repeatedly implicated a number of genes in DN, including ELMO1, PVT1, ACACB and SLC12A3, as well as specific candidate genes ACE, CNDP1, NOS3, SOD2, APOE and PRKCB [4,10]. The estimated sibling risk ratio for DN is about 2.3 , but only a modest proportion of that risk is explained by the above loci. Herein, the Family Investigation of Nephropathy and Diabetes (FIND) reports the largest genomewide linkage study for DN and urine albumin:creatinine ratio (ACR) in African Americans (AA), Southwest American Indians (AI), European Americans (EA) and Mexican Americans (MA) .
Materials and Methods
The FIND study design has been reported . Families of probands with DN with a diabetic sibling with or without nephropathy were recruited from eleven participating investigative centers. Living parents and other relatives (i.e. avuncular, cousin, half-sibling and grandparental affected pairs) were recruited when available. Recruitment was performed according to the principles of the Declaration of Helsinki. Written informed consent was obtained from all subjects with approval of the institutional review board at each center and the FIND Genetic Analysis and Data Coordinating Center.
Measurement of DN-related phenotypes in FIND has been described [12,13] (online suppl. material: Methods). The total available sample, with DN scored as affected or unaffected, included 2,616 individuals in 1,235 pedigrees across all ethnic groups (table 1), comprising 1,435 full-sib pairs and 188 half-sib pairs. Whether an individual was considered affected or unaffected with DN in linkage analysis depended on the participant category . Probands met the FIND criteria for DM and had either biopsy-proven DN, ESRD attributable to DN or chronic kidney disease attributable to DN. To be considered ‘unaffected’, individuals were required to have DM for at least 10 years without evidence of kidney disease, as ascertained through history, estimated glomerular filtration rate and urine ACR.
|Table 1. Description of the FIND study sample|
Urine was not collected from ESRD patients and a urine ACR of 3 g/g was assigned to these subjects. ACR values greater than 3 g/g in non-ESRD cases with DN were Winsorized to 3 g/g.
The binary DN trait was adjusted for sex, and the quantitative ACR trait for sex and age at DM diagnosis in all analyses. We did not adjust for use of angiotensin-converting enzyme inhibitors or angiotensin II receptor blockers, as information on these medications was not available for all participants.
We conducted the linkage and association analyses using the Illumina IV panel of approximately 5,500 autosomal single nucleotide polymorphisms (SNPs) that passed initial quality control criteria (online suppl. table S2). The Illumina IV markers were selected for uniform spacing across the genome  and high information content across multiple ethnic groups . The Haseman-Elston regression approach , implemented in SIBPAL (part of the S.A.G.E. software package), was the primary analysis method. The Haseman-Elston linkage test was performed using multipoint IBD sharing estimates, which combine genetic information across multiple linked SNPs, obtained from reduced sets of SNPs designed to maximize information in the absence of linkage disequilibrium. The W4 weighting scheme was used to adjust for dependence of sib-pairs within a sibship and of squared trait sums and differences as previously described , except that half-sib pairs were included in the analysis where possible. Linkage analysis was first attempted for each study sample including half-sib pairs. The analysis was repeated without half-sibs if model instability was encountered.
Linkage analyses were computed separately in each ethnic group. In addition to asymptotic p values, statistical significance was estimated via permutation testing, constructing an approximate empirical distribution of the test statistic. Because SIBPAL estimates empirical p values on full-sib pairs only, we did not estimate empirical p values for analyses including half-sib pairs. Fisher’s method was used to combine p values across ethnic groups. Following the linkage scans, we performed association testing using all SNPs that passed quality control criteria and the approach implemented in ASSOC, adjusting for familial relationships within a linear mixed model framework (online suppl. material: Methods). Because of differences in allele frequencies and linkage disequilibrium structure, the sets of SNPs contributing to linkage and association analyses differed among ethnic groups (online suppl. table S2).
The FIND sample comprised 2,616 participants in 1,235 pedigrees qualifying as affected or unaffected for DN, and 3,089 subjects with measures of urine ACR. Although the entire FIND sample did not measure autoantibodies diagnostic for type 1 diabetes, roughly 90–95% of DN-affected individuals were believed to have type 2 diabetes based on clinical criteria and a type 1 diabetes prevalence of 4% observed in a subsample of 857 individuals . Significant differences in BMI, age at DM diagnosis and DM duration prior to recruitment were present between probands and affected family members (p ≤ 2 × 10–4; table 1). Relatives of probands affected and unaffected for DN had significant differences in sex ratio, BMI and age at DM diagnosis. Significant differences in age were not detected between probands and DN-affected relatives, nor between DN-affected and -unaffected relatives (p > 0.05).
The FIND sample included 1,435 full-sib pairs and 113 half-sib pairs contributing to the genomewide linkage analysis for DN (online suppl. table S3), and 2,201 full-sib pairs and 268 half-sib pairs in the linkage analysis for urine ACR (online suppl. table S4). Subsets of these samples were included in preliminary scans for these traits .
The major ethnic-group-specific linkage signals for the DN binary trait, adjusted for sex, occurred on chromosomes 6p, 7p, 7q and 11p (fig. 1; table 2). The EA sample provided the strongest overall evidence for linkage on chromosome 6p (asymptotic p = 0.00015 at 24.9 cM, equivalent to a LOD score of 2.84; fig. 2a). Permutation analysis yielded a slightly more significant empirical p value of 8.0 × 10–5, LOD = 3.09. The most significant peak on 7p in AI was nearly as strong (p = 0.00016 at 15.5 cM, LOD = 2.81; fig. 2b). Several other regions showed suggestive evidence for linkage in AI, including 1q, 15q and 18p (fig. 2c). Top results in AA and MA were less robust. No peak in AA had an asymptotic p value ≤0.001; whereas chromosome 11p showed suggestive evidence for linkage in MA (fig. 1). The most significant evidence for linkage across all four populations, on chromosome 7p (p = 0.0019, LOD = 1.81), was driven primarily by the AI group, with a smaller contribution from AA (fig. 2c).
|Table 2. Major linkage peaks for DN and urine ACR|
|Fig. 1. Autosomal genomewide model-free linkage scan for DN. From top to bottom, panels display results for AA, AI, EA and MA samples, as well as pooled Fisher p values from all study samples (All).|
|Fig. 2. Detailed linkage and association results for DN for segments of chromosomes 6 (a), 7 (b), 18 (c, e) and 22 (d, f, g). Heavy lines trace linkage results for AA, AI, EA and MA samples, as well as combined Fisher p values (alternating dots and dashes). Tick marks at the top show the location of Illumina IV SNPs. In panel e, results from association analysis for AI (·) are shown along with the linkage profile from panel c. In panels f and g, association results for AI and EA are shown (·), respectively, with the linkage profiles from panel d.|
The linkage analysis was followed by a sparse genomewide association scan. Association testing for DN detected SNPs with p < 0.0001 for at least one ethnic group in two locations on chromosome 18, and SNPs with p < 0.001 on several chromosomes (fig. 2; online suppl. fig. 3; online suppl. table S5). No SNP with p < 0.001 in the association analysis coincided with a region showing evidence for linkage in the same ethnic group. In the ethnicity-combined association analysis, Fisher p values below 0.001 were observed for SNPs on chromosomes 5, 17 and 18 (online suppl. fig. S3; online suppl. table S5). All of these results were driven by a single population.
The most significant association result for AI (rs1241893 on chromosome 18; 22.56 cM, 6.87 Mb, p = 3.0 × 10–5) lies between two regions with modest evidence for linkage in AI (fig. 2e). Two SNPs spaced 2.7 kb apart on chromosome 18 and exhibiting low p values for AA (rs1662910 and rs948438; 57.98 cM and 33.10 Mb) did not correspond to linkage peaks in any ethnic group. In contrast, chromosome 22 marker rs5769116 (62.92 cM, 45.44 Mb) revealed moderate evidence for both linkage (p = 0.0018, LOD = 1.83) and association (p = 0.0019) in EA (fig. 2g) some 10 Mb distal to MYH9/APOL1 – genes associated with several nondiabetic kidney diseases [9,18,19,20].
In the genomewide linkage analysis for the quantitative ACR trait, suggestive evidence for linkage was obtained on chromosomes 3, 7, 16 and 22 (fig. 3; table 2). As for DN, the strongest population-specific linkage peak for urine ACR was seen in EA – in this case on chromosome 7q (p = 0.00011 at 99.3 cM, LOD = 2.96; fig. 4b). This signal overlaps the broad peak of linkage for DN in AA, but evidence for linkage to DN in EA was minimal (fig. 2b). The AA sample yielded strong evidence for linkage on chromosome 3p (p = 0.00018 at 95.5 cM, LOD = 2.76; fig. 4a), and weaker evidence on chromosomes 16q and 21q (fig. 4c). The most significant evidence for linkage in the MA sample was found on chromosome 22 (40.00 cM; empirical p = 0.00058, LOD = 2.29). The peak occurred approximately 1 Mb proximal to MYH9 (fig. 4d); the asymptotic and empirical p values for linkage at rs735853, within the MYH9 gene, were 0.0031 and 0.0014 (LOD = 1.92), respectively. Nonetheless, no test of association for chromosome 22 SNPs for the MA sample yielded a p value less than 0.01 (online suppl. fig. 4; online suppl. table S6).
|Fig. 3. Genomewide linkage scan for ACR, adjusted for sex and age at DM diagnosis. See legend of fig. 1 for explanation of plots and symbols.|
|Fig. 4. Detailed linkage and association results for ACR, on segments of chromosomes 3 (a), 7 (b), 21 (c, e) and 22 (d). Axes and lines are as in figure 2. In panel e, results from association analysis for the AA sample are shown (·) with the AA linkage profile from panel c.|
Weak evidence of linkage for urine ACR in AI (p = 0.02) reinforced the AA chromosome 21q linkage peak (rs220271; 42.36 Mb), providing the strongest evidence for linkage in the ethnicity-combined sample (p = 0.0010 at 59.28 cM, LOD = 2.06; fig. 4c). This peak coincided with the second most significant association result in AA, but there were no other strong association results in the region (fig. 4e). The most significant population-specific p value for association with urine ACR, found on chromosome 4p in AI (p = 1.3 × 10–5 at 38.65 Mb, 58.04 cM), was of borderline genomewide significance (p = 0.06 after Bonferroni correction for 4,925 tests; online suppl. fig. 4; online suppl. table S6). Combined association results across all four populations revealed p values ≤0.001 on chromosomes 6, 8, 9, 11 and 21. All signals were primarily driven by a single ethnic group, except for chromosomes 11 and 21 (online suppl. table S6). On chromosome 11, both the EA and MA samples gave significant results at rs722317 (15.88 Mb, 24.27 cM; p = 0.00046 and 0.0026, respectively) to yield a Fisher p value of 7.3 × 10–5, the most significant combined association result for ACR. rs2250226 on chromosome 21 yielded p values of 0.0025 and 0.0013, respectively, in AA and MA (Fisher p = 0.00028).
The FIND results from genomewide linkage and sparse association scans in four American ethnic groups with DN confirm previous linkage findings and suggest new genetic loci contributing to DN and related phenotypes. Suggestive evidence for linkage was observed on chromosomes 6p and 7p for DN, and on 3p, 7q, 16q and 22q for urine ACR in one or more ethnic groups. No region displayed consistent evidence for linkage across all 4 ethnic groups, although the modest urine ACR linkage and association peak on chromosome 21 in AA received weak support from AI. Significant evidence for linkage was detected in overlapping regions of chromosome 7q for both DN and urine ACR in AA and EA families, respectively.
Expanding the FIND family study samples from the partial cohort previously examined , altered the pattern of linkage across the autosomal genome. Newly detected linkage peaks with equivalent LOD scores ≧2.0 were observed on chromosomes 1q, 7p, 10p, 11p and 15q for DN, and on 2q, 3p, 16q and 22q for urine ACR. The full FIND sample provided greater power to detect genetic influences, as did modification of SIBPAL to include half-sib pairs in the Haseman-Elston regression. Moreover, use of a relatively dense set of 6,000 SNPs as a linkage panel, as opposed to a set of 400 microsatellite markers, ensures that a high level of information about linkage is available uniformly across each chromosome. The linkage peaks previously detected in a subset of these cases on chromosome 14q23 for DN and on 15q26 for urine ACR, now less prominent, were likely false positives, although genetic heterogeneity could also influence results. Similarly, there was no substantial overlap between the linkage findings reported here and in a linkage analysis of estimated glomerular filtration rate on the partial cohort .
In general, these results did not closely overlap genomic regions implicated in previous non-FIND linkage studies for DN  (online suppl. table S1). Our linkage signals on chromosome 7p for DN in AI and 22q12.3 for urine ACR in MA, coincided with those for glomerular filtration rate and ACR, respectively, in type 2 diabetes from 63 mostly Caucasian extended pedigrees [22,23]. The strongest evidence for linkage in this study occurred on chromosome 6p, with an empirical p value for linkage of 8.0 × 10–5 (LOD = 3.09) in EA for the binary DN trait.
The current report is the first of a linkage signal of such magnitude near the p terminus of chromosome 6 in a genomewide linkage study of DN or nondiabetic kidney disease. However, a genomewide scan for insulin resistance reported modest evidence for linkage in this location for acute insulin response to glucose in Hispanic families (LOD = 1.72 at 15 cM) . Suggestive evidence for linkage to urine ACR was found on chromosome 22 in MA (p = 5.9 × 10–4, LOD = 2.29), with peak linkage near the nonmuscle myosin heavy chain 9/apolipoprotein L1 (MYH9/APOL1) gene region implicated in diabetic and nondiabetic ESRD and focal segmental glomerulosclerosis in AA [9,18,19,25,26]. Additional supporting evidence was observed in EA using the binary DN trait. However, no association signal was detected, which is reasonable since none of the reported MYH9/APOL1 kidney disease risk variants were included in our genotyping panels.
Several results from the FIND analyses occur in regions previously associated with DN phenotypes. Two SNPs near the candidate gene carnosinase 1 (CNDP1) showed evidence for association in at least one FIND cohort. The first SNP, rs999647, is 49 kb distal to CNDP1, with association specific to DN in AI families (p = 0.0057) and weak association with urine ACR in EA (p = 0.031; Fisher combined p = 0.028 over all cohorts). A neighboring SNP, rs872994, 919 kb distal to CNDP1 (at position 73,171,838, human genome build 18), displayed weak evidence for association in three ethnic groups, AI (p = 0.053), EA (p = 0.017) and MA (p = 0.033), with a combined Fisher p of 0.0043.
The chromosome 3p DN association in the AI population occurred at rs892605, approximately 400 kb proximal to the ghrelin (GHRL) gene, for which association with kidney function in DM, but not DM itself, has been reported [27,28,29].
Genomewide linkage scans for DN and related traits differ in many ways, including pedigree structures recruited, ascertainment criteria, definitions of DN and measurement of kidney function, diabetes type and duration. In contrast to other publications, FIND is a multiethnic sample of DN-affected and -unaffected subjects with extremes of phenotypic severity . This genomewide linkage and sparse association scan in the full FIND study sample provides evidence for both novel (e.g. on chromosomes 3p and 6p) and previously detected (chromosomes 7p and 22) genetic loci predisposing to renal impairment and DN. Genomewide association studies have the most power to detect genetic variants with relatively modest effect, but require that causal alleles be well captured (tagged) by the marker panel, which is true in practice mostly for common variants. The linkage study presented here only has power to detect variants with relatively large effect, but can do so regardless of the population allele frequency and number of tightly linked variants involved. These results should provide guidance for further exploration of the genetic basis for diabetic kidney disease using novel technologies.
We thank all FIND participants. This study was supported by grants U01DK57292, U01DK57329, U01DK057300, U01DK057298, U01DK057249, U01DK57295, U01DK070657, U01DK057303, U01DK070657, U01DK57304 and DK57292-05 from the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) and, in part, by the Intramural Research Program of the NIDDK. This project has been funded in whole or in part with federal funds from the National Cancer Institute, National Institutes of Health (NIH), under contract N01-CO-12400 and the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research. This work was also supported by the National Center for Research Resources for the General Clinical Research Center grants: Case Western Reserve University, M01-RR-000080; Wake Forest University, M01-RR-07122; Harbor-University of California, Los Angeles Medical Center, M01-RR-00425; College of Medicine, University of California, Irvine, M01-RR-00827–29; University of New Mexico, HSC M01-RR-00997; and Frederic C. Bartter, M01-RR-01346. Genotyping was performed by the Center for Inherited Disease Research, which is fully funded through a federal contract from the NIH to Johns Hopkins University (N01-HG-65403). The results of this analysis were obtained using the S.A.G.E. package of genetic epidemiology software, which is supported by a U.S. Public Health Service Resource Grant (RR03655) from the National Center for Research Resources.
Supplementary information is available at the American Journal of Nephrology web site.
Dr. Sudha Iyengar
Department of Epidemiology and Biostatistics, Case Western Reserve University
Wolstein Research Building, Rm # 1300, 10900 Euclid Avenue
Cleveland, OH 44106-7281 (USA)
Tel. +1 216 368 5636, E-Mail firstname.lastname@example.org
Received: January 4, 2011
Accepted: February 23, 2011
Published online: March 31, 2011
Number of Print Pages : 9
Number of Figures : 4, Number of Tables : 2, Number of References : 29
Additional supplementary material is available online - Number of Parts : 1
American Journal of Nephrology
Vol. 33, No. 5, Year 2011 (Cover Date: May 2011)
Journal Editor: Bakris G. (Chicago, Ill.)
ISSN: 0250-8095 (Print), eISSN: 1421-9670 (Online)
For additional information: http://www.karger.com/AJN