Family-Based Association Analysis Confirms the Role of the Chromosome 9q21.32 Locus in the Susceptibility of Diabetic Nephropathy

Abstract

A genome-wide association scan of type 1 diabetic patients from the GoKinD collections previously identified four novel diabetic nephropathy susceptibility loci that have subsequently been shown to be associated with diabetic nephropathy in unrelated patients with type 2 diabetes. To expand these findings, we examined whether single nucleotide polymorphisms (SNPs) at these susceptibility loci were associated with diabetic nephropathy in patients from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection. Six SNPs across the four loci identified in the GoKinD collections and 7 haplotype tagging SNPs, were genotyped in 66 extended families of European ancestry. Pedigrees from this collection contained an average of 18.5 members, including 2 to 14 members with type 2 diabetes. Among diabetic family members, the 9q21.32 locus approached statistical significance with advanced diabetic nephropathy (P = 0.037 [adjusted P = 0.222]). When we expanded our definition of diabetic nephropathy to include individuals with high microalbuminuria, the strength of this association improved significantly (P = 1.42×10⁻³ [adjusted P = 0.009]). This same locus also trended toward statistical significance with variation in urinary albumin excretion in family members with type 2 diabetes (P = 0.032 [adjusted P = 0.192]) and in analyses expanded to include all relatives (P = 0.019 [adjusted P = 0.114]). These data increase support that SNPs identified in the GoKinD collections on chromosome 9q21.32 are true diabetic nephropathy susceptibility loci.

About the Authors:

Marcus G. Pezzolesi

* E-mail: [email protected]

Affiliations Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America, Department of Medicine, Harvard Medical School, Boston, Massachusetts, United States of America

Jackson Jeong

Affiliation: Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America

Adam M. Smiles

Affiliation: Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America

Jan Skupien

Affiliations Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America, Department of Medicine, Harvard Medical School, Boston, Massachusetts, United States of America

Josyf C. Mychaleckyj

Affiliation: Center for Public Health Genomics, University of Virginia School of Medicine, Charlottesville, Virginia, United States of America

Stephen S. Rich

Affiliation: Center for Public Health Genomics, University of Virginia School of Medicine, Charlottesville, Virginia, United States of America

James H. Warram

Affiliation: Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America

Andrzej S. Krolewski

Affiliations Research Division, Joslin Diabetes Center, Boston, Massachusetts, United States of America, Department of Medicine, Harvard Medical School, Boston, Massachusetts, United States of America

Introduction

Increased urinary albumin excretion, in both the microalbuminuric and proteinuric ranges, is a hallmark of diabetic nephropathy (DN) [1], [2]. Clinically, DN is a progressive disease that advances through characteristic stages. For many diabetic patients, elevated urinary albumin excretion is coupled with declining renal function. In a large proportion of these individuals, renal function continues to deteriorate until end-stage renal disease (ESRD) is reached.

Despite a large body of evidence that favors a genetic basis for susceptibility to DN, identification of the genetic factors that contribute to its risk has proven challenging [3]–[12]. While no single, major DN susceptibility gene has yet been identified, growing support for several loci identified though genome-wide surveys of common genetic variants has recently begun to emerge [13]–[22]. Included among the studies contributing to this success is our recent genome-wide association (GWA) scan of unrelated type 1 diabetic (T1D) subjects from the Genetics of Kidneys in Diabetes (GoKinD) study collections [18]. In this report, we identified strong associations at several common single nucleotide polymorphisms (SNPs, minor allele frequencies >5%) located across four distinct chromosomal regions. Three of these loci, located on chromosome 9q21.32 near the FRMD3 gene, chromosome 11p15.4 at the CARS gene, and chromosome 13q33.3 at the MYO16/IRS2 locus, have since been confirmed in multiple diverse collections of unrelated T1D or type 2 diabetic (T2D) patients [18], [20], [22]. A more recent meta-analysis of T1D nephropathy, defined as end-stage renal disease (ESRD), in European-derived populations, however, failed to confirm these, as well as several other, previously reported genetic associations; reinforcing the need for further investigation of these and other loci to truly understand their role in the genetic basis of DN [23].

To address this need, we chose to extend our focused evaluation of the loci identified in GoKinD to a family-based association study of patients from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection. In addition to dichotomized comparisons of DN status, we investigated whether any of these loci were associated with quantitative variation in urinary albumin in this collection.

Materials and Methods

Study Patients and Ethics Statement

The present study investigated 1,221 individuals (798 with direct genotype and phenotype information) from 66 extended families of European ancestry from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection. The protocols and informed consent procedures used in this study were approved by the Committee on Human Subjects of the Joslin Diabetes Center.

The recruitment of the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection has previously been described [9], [11], [12], [24]. Briefly, between 1993 and 2003, families with an apparent autosomal dominant mode of inheritance of T2D, irrespective of their nephropathy status, were recruited to the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection through T2D probands receiving medical care at the Joslin Clinic. After obtaining informed written consent, trained recruiters administered previously described study protocols that included a structured interview, seated blood pressure measurements, and the collections of blood and urine samples. ESRD status for members of this collection was updated as of August 2008 through the United States Renal Data System.

Classification of Nephropathy Status

Methods for measuring albumin and creatinine in a random urine sample for determination of the albumin-to-creatinine ratio (ACR) and defining normoalbuminuria, microalbuminuria, or proteinuria were described previously [25]. ACR values were used to assign albuminuria status to all individuals included in our analysis; individuals with ACR values less than 30 µg/mg, between 30 µg/mg and 300 µg/mg, between 100 µg/mg and 300 µg/mg, and above 300 µg/mg were considered normoalbuminuric, microalbuminuric, high microalbuminuric, and proteinuric, respectively. Individuals with ESRD were assigned ACR values of 3500 µg/mg. For quantitative trait analyses, a log transformation was applied to the measured/assigned ACR values.

Genotyping

Six SNPs across the four loci identified in the GoKinD collections were selected for inclusion in the present study; including rs39075 on chromosome 7p14.3, rs1888747 and rs10868025 on chromosome 9q21.32, rs451041 on chromosome 11p15.4, and rs1411766 and rs9521445 on chromosome 13q33.3. Seven additional haplotype tagging SNPs (two on chromosomes 7p14.3, 11p15.4, and 13q33.3 and one on chromosome 9q21.32) were selected using Haploview [26] to capture the major haplotypes (haplotype frequencies ≥0.05) for the linkage disequilibrium (LD) blocks containing the SNPs identified in GoKinD. All thirteen SNPs were genotyped using Taqman (Applied Biosystems, Foster City, CA) technology by the Genetics Core of the Diabetes and Endocrinology Research Center at the Joslin Diabetes Center in accordance with the manufacturer's protocols.

Statistical Analysis

Each SNP was tested for deviation from Hardy-Weinberg equilibrium using a chi-square goodness-of-fit test. Family-based single-marker association tests were performed using the FBAT software under an additive model using a conservative empirical variance estimator to test the null hypothesis of no linkage and no association [27]. For all dichotomous trait analyses, allele transmissions from parent to affected and unaffected family members were contrasted by weighting their contribution to the FBAT test statistic using the estimated population prevalence of DN among diabetic individuals (i.e., 30%; ‘Affected and Unaffected’ analyses). Family-based association testing of allele transmission from parents to only affected offsprings was also performed (i.e., ‘Affecteds Only’ analyses). The HBAT procedure in FBAT was used to estimate haplotype frequencies and perform haplotype-specific and global tests of association. P-values <8.33×10⁻³ (0.05/6) were considered statistically significant.

Results

A total of 1,221 individuals from 66 extended families of European ancestry from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection, including 798 (382 non-diabetic and 416 diabetic) members with direct genotype and phenotype information, were included in the present study. Pedigrees from these families included an average of 18.5 members, ranging in size from 6 to 39 members, and formed a total of 318 nuclear families. Four to 26 individuals from each family, including 2 to 14 members with diabetes, had DNA available for genotyping. The mean age of diabetes diagnosis within these families, was 43.4±16.8 years.

Ninety-seven (23.3%) diabetic individuals were considered to have advanced DN (proteinuria, n = 40, or ESRD, n = 57) while 312 (75.0%) were classified as non-DN controls. To improve power to detect significant associations between DN and SNPs identified in the GoKinD collections, we also expanded our nephropathy phenotype to include individuals with less severe nephropathy. For these comparisons, dichotomized cases included 28 additional individuals with high microalbuminuria. Seven diabetic individuals did not have ACR data available for classification of their nephropathy status or for the quantitative analysis of this trait; these individuals were excluded from all analyses. Proteinuria developed in 3 (0.8%) non-diabetic individuals. Clinical characteristics for examined members of the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection included in this study are provided in Table 1.

[Figure omitted. See PDF.]

Table 1. Clinical characteristics of 798 examined members from 66 families from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection.

https://doi.org/10.1371/journal.pone.0060301.t001

The distribution of relative pairs in the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection based on their relationship and DN status are provided in Table 2. These 66 extended families generated a total of 6,421 relative pairs; 1,026 of whom are concordant for diabetes. The diabetic relative pairs include 53 relative pairs that are concordant for advanced DN, 239 relative pairs that are disconcordant for advanced DN, 88 relative pairs concordant for DN, and 329 relative pairs disconcordant for DN. Additionally, there are 15 relative pairs concordant for ESRD and 134 relative pairs disconcordant for this phenotype.

[Figure omitted. See PDF.]

Table 2. Summary of the relative pairs in the 66 families from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection according to diabetes and nephropathy status.

https://doi.org/10.1371/journal.pone.0060301.t002

No deviation from Hardy-Weinberg equilibrium was observed among the 13 SNPs included in this study in the entire collection or in analyses performed separately in affected and unaffected individuals (P>0.05).

Family-based association analyses of 6 SNPs across the four loci identified in the GoKinD GWA scan were performed in diabetic relatives as well as in all relatives combined. Among diabetic family members, rs1888747 on chromosome 9q21.32 showed evidence of association with advanced nephropathy among diabetic family members (affecteds only: P = 0.029 [adjusted P = 0.174], Z = 2.18; affecteds and unaffecteds: P = 0.037 [adjusted P = 0.222], Z = 2.08, Table 3). When we expanded our definition of DN to include individuals with high microalbuminuria, the strength of this association improved significantly (affecteds only: P = 1.74×10⁻³ [adjusted P = 0.010], Z = 3.13; affecteds and unaffecteds: P = 1.42×10⁻³ [adjusted P = 0.009], Z = 3.19, Table 4). In both comparisons, the associations on 9q21.32 were in the same direction as initially reported in the GoKinD collections. Although no other SNPs achieved statistical significance, variants on 7p14.3 and 13q33.3 shared directionality with associations identified in GoKinD (Table 3 and Table 4). Among all family members, only rs1888747 was nominally associated with nephropathy in analyses that included high microalbuminurics (affecteds only: P = 0.026 [adjusted P = 0.156], Z = 2.23; affecteds and unaffecteds: P = 0.017 [adjusted P = 0.102], Z = 2.38, Table S1 and Table S2). Lastly, dichotomous analyses restricted to diabetic family members with ESRD showed modest evidence of a statistical association with rs1888747 (affecteds only: P = 0.036, Z = 2.09; affecteds and unaffecteds: P = 0.046, Z = 2.00, Table 5); this SNP, however, did not achieve statistical significance when a conservative Bonferroni correction was applied (adjusted P>0.05).

[Figure omitted. See PDF.]

Table 3. Family-based association results between DN-associated SNPs and advanced nephropathy (normoalbuminuria vs. proteinuria/ESRD) among diabetic family members.

https://doi.org/10.1371/journal.pone.0060301.t003

[Figure omitted. See PDF.]

Table 4. Family-based association results between DN-associated SNPs and nephropathy (normoalbuminuria vs. high microalbuminuria/proteinuria/ESRD) among diabetic family members.

https://doi.org/10.1371/journal.pone.0060301.t004

[Figure omitted. See PDF.]

Table 5. Family-based association results between DN-associated SNPs and ESRD among diabetic family members.

https://doi.org/10.1371/journal.pone.0060301.t005

In quantitative trait analyses, rs1888747 on chromosome 9q21.32 was similarly shown to be associated with logACR among diabetic family members (P = 0.030 [adjusted P = 0.180], Z = 2.17, Table 6) and in analyses extended to all family members in the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection (P = 0.017 [adjusted P = 0.102], Z = 2.39, Table S3). Controlling for reported ACE inhibitor treatment weakened the association between rs1888747 and albuminuria among diabetic family members and all family members (P = 0.217 [adjusted P = 1.00], Z = 1.24 and P = 0.121 [adjusted P = 0.726], Z = 1.55, respectively) while associations at rs10868025 (also located on chromosome 9q21.32) improved (diabetic family members: P = 0.012 [adjusted P = 0.072], Z = 2.50; all family members: P = 0.033 [adjusted P = 0.198], Z = 2.13).

[Figure omitted. See PDF.]

Table 6. Family-based association analysis between DN-associated SNPs and logACR among diabetic family members.

https://doi.org/10.1371/journal.pone.0060301.t006

No additional associations were observed between the remaining GoKinD SNPs and logACR in analyses either restricted to diabetic individuals or in those extended to include all available family members.

To further examine the susceptibility loci identified in GoKinD, we genotyped haplotype tagging SNPs across each locus and performed family-based multi-marker analyses for each nephropathy phenotype using the HBAT procedure in FBAT. Haplotypes formed by genotyped SNPs on chromosomes 7p14.3, 11p15.4, and 13q33.3 were not associated with any of the examined nephropathy phenotypes in analyses of diabetic individuals or in those performed in all family members (global P>0.05, data not shown). In contrast, haplotypes on chromosome 9q21.32 showed evidence of association with advanced nephropathy, nephropathy and logACR among diabetic family members (Table 7 and Table 8). While no haplotype on 9q21.32 was more strongly associated with nephropathy or logACR than the individual SNPs at this locus, the GTA haplotype, which is comprised of the risk alleles for both rs1888747 and rs10868025, is more strongly associated with an increased risk of advanced nephropathy than any individual SNP at this locus (affecteds only: P = 0.012, Z = 2.50; affecteds and unaffecteds: P = 0.018, Z = 2.37, Table 7). None of the genotyped haplotype tagging SNPs were associated with nephropathy in single marker analyses (Tables S4 through S9).

[Figure omitted. See PDF.]

Table 7. Family-based haplotype analysis between chromosome 9q21.32^* haplotypes and advanced nephropathy (normoalbuminuria vs. proteinuria/ESRD) and nephropathy (normoalbuminuria vs. high microalbuminuria/proteinuria/ESRD) among diabetic family members.

https://doi.org/10.1371/journal.pone.0060301.t007

[Figure omitted. See PDF.]

Table 8. Family-based haplotype analysis between chromosome 9q21.32^* haplotypes and logACR among diabetic family members.

https://doi.org/10.1371/journal.pone.0060301.t008

Discussion

Our GWA scan of the GoKinD collections identified strong association at four distinct chromosomal regions, including loci on chromosomes 9q21.32, 11p15.4, and 13q33.3 that have since been confirmed in multiple collections comprised of unrelated T1D or T2D subjects [18], [20]–[22]. In the present report, we extend these findings further by providing additional support for associations at chromosome 9q21.32 in a large collection of related T2D patients from the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection. In this study, a statistically significant association was observed between this locus and the risk of high microalbuminuria, proteinuria, and ESRD among diabetic individuals in these families.

DN is well-recognized to be a complex disease, characterized by both abnormalities in urinary albumin excretion and declining renal function. While most patients with DN exhibit some degree of elevated urinary albumin excretion, for some, this abnormality is concurrent with impaired renal function; a subset of these individuals eventually require renal replacement therapy. This phenotypic heterogeneity suggests that multiple factors, some genetic and some non-genetic, contribute to the distinct stages of DN. In the present study, we report evidence of association at the 9q21.32 locus with advanced DN (i.e., proteinuria and ESRD). The strength of these associations improved significantly when we expanded our nephropathy phenotype to include individuals with less severe DN (i.e., patients with high microalbuminuria, ACR values >100 µg/mg). Based on these findings, we hypothesize that variation at this locus contributes to the early stages of nephropathy in diabetes. Furthermore, we hypothesize that other genetic factors are likely involved in the progression of DN and the decline in renal function that accompanies the latter stages of this disease process.

Associations at common variants on the 9q21.32 locus have now been confirmed in four distinct collections; T1D patients from the GoKinD collections, the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection, Caucasian participants from the Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) study [18], and African Americans with ESRD due to T2D [21]. In addition to these studies, a strong protective affect of the rs1888747 DN-risk allele observed in populations of European ancestry has recently been reported in unrelated Japanese T2D patients with microalbuminuria [20]; a notable finding that may be attributed to allelic heterogeneity resulting from the ancestral backgrounds of these two ethnic groups. Despite being quite underpowered, in this study Maeda et al. identified a strong association at rs1888747 in a comparison of only 32 microalbuminuric patients who progressed to overt proteinuria and 168 individuals who remained microalbuminuric (i.e., non-progressors).

Additional support for a role of this locus in nephropathy comes from a genetic study of albuminuria quantitative trait loci (QTL) performed in an intercross of albuminuria resistant and susceptible mouse strains [28]. Using this approach, Sheehan et al. were able to localize a QTL associated with increased urinary albumin on mouse chromosome 4, a region homologous to the DN-associated locus on chromosome 9q21.32 seen in human populations. This striking concordance between human and mouse suggests that a common disease mechanism may link these renal damage phenotypes.

A major challenge in dissecting the genetic basis of complex traits, including DN, is that many of the common variants that have been reproducibly shown to be associated with disease explain only a modest proportion of the overall risk of disease. Although less powerful than population-based designs, family-based approaches such as the one employed in our study will prove crucial in uncovering variants that have much larger contributions to the genetic basis of disease as rare variants that are expected to explain a larger proportion of the heritability of a given phenotype are enrich in related individuals.

The present study has modest power (<80%) to detect similar effects as those observed by Pezzolesi et al. [18]. Despite this limitation, our analysis increases support that previously reported associations at the 9q21.32 locus are genuine diabetic nephropathy susceptibility loci and, given the limited power of our study, suggest that the true effect sizes attributed to variants at this locus may in fact be larger than previously estimated. More specifically, the associations identified in the GoKinD collections at the 9q21.32 locus and confirmed in the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection occur at common SNPs (risk allele frequency ≈ 60–70%) of modest effect (odds ratio = 1.40–1.50). It remains possible that this association is due to the presence of a ‘synthetic association’ caused by one or more rare variants located some distance from the observed associations [29]. Carriers of the rs1888744 risk allele in our families are potential future candidates for targeted next-generation sequencing of this locus to identify rare variants that may explain a substantial proportion of the heritability of DN risk observed in these families.

In contrast to studies of unrelated cases and controls, family-based designs, such as the one used in the present study, are robust to population admixture and stratification. Additionally, family members also tend to have more homogeneity of environmental factors that could confound genetic associations with the phenotype of interest. Despite these advantages, we acknowledge that the present study is not without its limitations. Patients in Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection were recruited irrespective of their nephropathy status and their phenotypic characteristics were primarily derived from baseline data taken at the time of enrollment. Understanding the natural history of diabetic nephropathy and recognizing the limitations of this cross-sectional assessment of kidney phenotypes, we chose to use the USRDS database to track individuals that progressed to ESRD. Individuals who might have progressed to ESRD but refused renal replacement therapy are not represented in USRDS and, depending on their renal status at baseline, may have misclassification of their renal phenotype. In lieu of longitudinal follow-up of all members of the Joslin Study of Genetics of Nephropathy in Type 2 Diabetes Family Collection, we recognize that the potential misclassification of these individuals is a limitation of our study. Importantly, however, we note that any phenotypic misclassification due to our inability to track participants who refused renal replacement therapy is independent of their carrier status of DN risk and/or non-risk alleles; the resulting non-differential misclassification could only have biased our findings toward the null. A second limitation worthy of discussion is the potential competing risk of death prior to the development of ESRD. Because of the high rate of mortality experienced by diabetic patients with ESRD, variants associated with this phenotype may alternatively be associated with survival in the presence of severe kidney disease. Although the present study lacks the power to formally assess whether the effects of variants at the 9q21.32 locus differ according to duration of ESRD, we have previously shown that the odds ratios for these variants are consistent across tertiles of ESRD duration in patients with T1D and ESRD, suggesting that these associations are not due to survival bias [18].

rs1888747 lies approximately 2 kilo-basepair upstream of FRMD3, a gene that we have previously shown to be expressed in both human kidney mesangial cells and proximal tubular cells [18]. Additional work has further demonstrated that FRMD3 is also expressed in human podocytes and that its protein product, the 4.1O protein, interacts with nephrin, podocin, and actin, suggesting this protein is involved in maintaining the function and integrity of the slit diaphragm (unpublished data). Most recently, we have shown that rs1888747's risk allele generates a transcription factor binding site in a repressive promote module that is shared by multiple members of the bone morphogenetic protein (BMP) signaling pathway; a pathway that has previously been implicated in the development of DN [30], [31]. Hierarchical clustering of expression data for FRMD3 and its co-expressed transcripts suggests that these genes are linked to early progression in DN [30]. Coupled with the strong association we observed at this locus in diabetic family members with less severe DN in the present study, we hypothesize that the 9q21.32 locus contributes to glomerular injury early in DN's pathogenesis.

In summary, our study provides further evidence that the 9q21.32 region is a susceptibility locus for DN. Coupled with its proximity to the association at this locus, FRMD3 appears to be both a strong positional and biological candidate gene for DN.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Family-based association results between DN-associated SNPs and advanced nephropathy (normoalbuminuria vs. proteinuria/ESRD) among all family members.

https://doi.org/10.1371/journal.pone.0060301.s001

(DOC)

Table S2.

Family-based association results between DN-associated SNPs and nephropathy (normoalbuminuria vs. high microalbuminuria/proteinuria/ESRD) among all family members.

https://doi.org/10.1371/journal.pone.0060301.s002

(DOC)

Table S3.

Family-based association analysis between DN-associated SNPs and logACR among all family members.

https://doi.org/10.1371/journal.pone.0060301.s003

(DOC)

Table S4.

Single marker family-based association analyses between haplotype tagging SNPs across the four GoKinD loci and advanced nephropathy among diabetic family members. Affecteds and unaffecteds analyses are presented.

https://doi.org/10.1371/journal.pone.0060301.s004

(DOC)

Table S5.

Single marker family-based association analyses between haplotype tagging SNPs across the four GoKinD loci and nephropathy among diabetic family members. Affecteds and unaffecteds analyses are presented.

https://doi.org/10.1371/journal.pone.0060301.s005

(DOC)

Table S6.

Single marker family-based association analyses between haplotype tagging SNPs across the four GoKinD loci and logACR among diabetic family members. Affecteds and unaffecteds analyses are presented.

https://doi.org/10.1371/journal.pone.0060301.s006

(DOC)

Table S7.

Single marker family-based association analyses between haplotype tagging SNPs across the four GoKinD loci and advanced nephropathy among all family members. Affecteds and unaffecteds analyses are presented.

https://doi.org/10.1371/journal.pone.0060301.s007

(DOC)

Table S8.

Single marker family-based association analyses between haplotype tagging SNPs across the four GoKinD loci and nephropathy among all family members. Affecteds and unaffecteds analyses are presented.

https://doi.org/10.1371/journal.pone.0060301.s008

(DOC)

Table S9.

Single marker family-based association analyses between haplotype tagging SNPs across the four GoKinD loci and logACR among all family members. Affecteds and unaffecteds analyses are presented.

https://doi.org/10.1371/journal.pone.0060301.s009

(DOC)

Author Contributions

Conceived and designed the experiments: MGP ASK. Performed the experiments: MGP. Analyzed the data: MGP ASK. Contributed reagents/materials/analysis tools: MGP ASK JJ AMS JS JCM SSR JHW. Wrote the paper: MGP.

Citation: Pezzolesi MG, Jeong J, Smiles AM, Skupien J, Mychaleckyj JC, Rich SS, et al. (2013) Family-Based Association Analysis Confirms the Role of the Chromosome 9q21.32 Locus in the Susceptibility of Diabetic Nephropathy. PLoS ONE 8(3): e60301. https://doi.org/10.1371/journal.pone.0060301

References

1. Krolewski AS, Warram JH (1997) Clinical Features and Epidemiology of Diabetic Nephropathy. In: Pickup JC WG, editor. Textbook of Diabetes.2nd ed.Oxford: Blackwell Scientific Publications. pp. 53.51–53.13.

2. Parving HH, Mauer M, Ritz E (2004) Diabetic Nephropathy. In: BM B, editor. Brenner and Rector's The Kidney.7th ed. Philadelphia: Elsevier.pp. 1777–1818.

3. Seaquist ER, Goetz FC, Rich S, Barbosa J (1989) Familial clustering of diabetic kidney disease. Evidence for genetic susceptibility to diabetic nephropathy. N Engl J Med 320: 1161–1165.

4. Quinn M, Angelico MC, Warram JH, Krolewski AS (1996) Familial factors determine the development of diabetic nephropathy in patients with IDDM. Diabetologia 39: 940–945.

5. Faronato PP, Maioli M, Tonolo G, Brocco E, Noventa F, et al. (1997) Clustering of albumin excretion rate abnormalities in Caucasian patients with NIDDM. The Italian NIDDM Nephropathy Study Group. Diabetologia 40: 816–823.

6. Canani LH, Gerchman F, Gross JL (1999) Familial clustering of diabetic nephropathy in Brazilian type 2 diabetic patients. Diabetes 48: 909–913.

7. Forsblom CM, Kanninen T, Lehtovirta M, Saloranta C, Groop LC (1999) Heritability of albumin excretion rate in families of patients with Type II diabetes. Diabetologia 42: 1359–1366.

8. Fogarty DG, Hanna LS, Wantman M, Warram JH, Krolewski AS, et al. (2000) Segregation analysis of urinary albumin excretion in families with type 2 diabetes. Diabetes 49: 1057–1063.

9. Fogarty DG, Rich SS, Hanna L, Warram JH, Krolewski AS (2000) Urinary albumin excretion in families with type 2 diabetes is heritable and genetically correlated to blood pressure. Kidney Int 57: 250–257.

10. Langefeld CD, Beck SR, Bowden DW, Rich SS, Wagenknecht LE, et al. (2004) Heritability of GFR and albuminuria in Caucasians with type 2 diabetes mellitus. Am J Kidney Dis 43: 796–800.

11. Placha G, Canani LH, Warram JH, Krolewski AS (2005) Evidence for different susceptibility genes for proteinuria and ESRD in type 2 diabetes. Adv Chronic Kidney Dis 12: 155–169.

12. Krolewski AS, Poznik GD, Placha G, Canani L, Dunn J, et al. (2006) A genome-wide linkage scan for genes controlling variation in urinary albumin excretion in type II diabetes. Kidney Int 69: 129–136.

13. Shimazaki A, Kawamura Y, Kanazawa A, Sekine A, Saito S, et al. (2005) Genetic variations in the gene encoding ELMO1 are associated with susceptibility to diabetic nephropathy. Diabetes 54: 1171–1178.

14. Hanson RL, Craig DW, Millis MP, Yeatts KA, Kobes S, et al. (2007) Identification of PVT1 as a candidate gene for end-stage renal disease in type 2 diabetes using a pooling-based genome-wide single nucleotide polymorphism association study. Diabetes 56: 975–983.

15. Millis MP, Bowen D, Kingsley C, Watanabe RM, Wolford JK (2007) Variants in the plasmacytoma variant translocation gene (PVT1) are associated with end-stage renal disease attributed to type 1 diabetes. Diabetes 56: 3027–3032.

16. Leak TS, Perlegas PS, Smith SG, Keene KL, Hicks PJ, et al. (2009) Variants in intron 13 of the ELMO1 gene are associated with diabetic nephropathy in African Americans. Ann Hum Genet 73: 152–159.

17. Pezzolesi MG, Katavetin P, Kure M, Poznik GD, Skupien J, et al. (2009) Confirmation of genetic associations at ELMO1 in the GoKinD collection supports its role as a susceptibility gene in diabetic nephropathy. Diabetes 58: 2698–2702.

18. Pezzolesi MG, Poznik GD, Mychaleckyj JC, Paterson AD, Barati MT, et al. (2009) Genome-wide association scan for diabetic nephropathy susceptibility genes in type 1 diabetes. Diabetes 58: 1403–1410.

19. Hanson RL, Millis MP, Young NJ, Kobes S, Nelson RG, et al. (2010) ELMO1 variants and susceptibility to diabetic nephropathy in American Indians. Mol Genet Metab 101: 383–390.

20. Maeda S, Araki S, Babazono T, Toyoda M, Umezono T, et al. (2010) Replication study for the association between 4 loci identified by a genomewide association study on European American subjects with type 1 diabetes and susceptibility to diabetic nephropathy in Japanese subjects with type 2 diabetes. Diabetes In Press.

21. Freedman BI, Langefeld CD, Lu L, Divers J, Comeau ME, et al. (2011) Differential effects of MYH9 and APOL1 risk variants on FRMD3 Association with Diabetic ESRD in African Americans. PLoS Genet 7: e1002150.

22. Pezzolesi MG, Poznik GD, Skupien J, Smiles AM, Mychaleckyj JC, et al. (2011) An intergenic region on chromosome 13q33.3 is associated with the susceptibility to kidney disease in type 1 and 2 diabetes. Kidney Int 80: 105–111.

23. Williams WW, Salem RM, McKnight AJ, Sandholm N, Forsblom C, et al. (2012) Association testing of previously reported variants in a large case-control meta-analysis of diabetic nephropathy. Diabetes 61: 2187–2194.

24. Klupa T, Warram JH, Antonellis A, Pezzolesi M, Nam M, et al. (2002) Determinants of the development of diabetes (maturity-onset diabetes of the young-3) in carriers of HNF-1alpha mutations: evidence for parent-of-origin effect. Diabetes Care 25: 2292–2301.

25. Warram JH, Gearin G, Laffel L, Krolewski AS (1996) Effect of duration of type I diabetes on the prevalence of stages of diabetic nephropathy defined by urinary albumin/creatinine ratio. J Am Soc Nephrol 7: 930–937.

26. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.

27. Laird NM, Horvath S, Xu X (2000) Implementing a unified approach to family-based tests of association. Genet Epidemiol 19 Suppl 1S36–42.

28. Sheehan S, Tsaih SW, King BL, Stanton C, Churchill GA, et al. (2007) Genetic analysis of albuminuria in a cross between C57BL/6J and DBA/2J mice. Am J Physiol Renal Physiol 293: F1649–1656.

29. Dickson SP, Wang K, Krantz I, Hakonarson H, Goldstein DB (2010) Rare variants create synthetic genome-wide associations. PLoS Biol 8: e1000294.

30. Martini S, Nair V, Patel SR, Eichinger F, Nelson RG, et al. (2013) From SNP to transcriptional mechanism: A model for FRMD3 in diabetic nephropathy. Diabetes In press.

31. Dolan V, Hensey C, Brady HR (2003) Diabetic nephropathy: renal development gone awry? Pediatr Nephrol 18: 75–84.