Whole-Genome Sequencing of Finnish Type 1 Diabetic Siblings Discordant for Kidney Disease Reveals DNA Variants associated with Diabetic Nephropathy

BACKGROUND

Several genetic susceptibility loci associated with diabetic nephropathy have been documented, but no causative variants implying novel pathogenetic mechanisms have been elucidated.

METHODS

We carried out whole-genome sequencing of a discovery cohort of Finnish siblings with type 1 diabetes who were discordant for the presence (case) or absence (control) of diabetic nephropathy. Controls had diabetes without complications for 15-37 years. We analyzed and annotated variants at genome, gene, and single-nucleotide variant levels. We then replicated the associated variants, genes, and regions in a replication cohort from the Finnish Diabetic Nephropathy study that included 3531 unrelated Finns with type 1 diabetes.

RESULTS

We observed protein-altering variants and an enrichment of variants in regions associated with the presence or absence of diabetic nephropathy. The replication cohort confirmed variants in both regulatory and protein-coding regions. We also observed that diabetic nephropathy-associated variants, when clustered at the gene level, are enriched in a core protein-interaction network representing proteins essential for podocyte function. These genes include protein kinases (protein kinase C isoforms ε and ι) and protein tyrosine kinase 2.

CONCLUSIONS

Our comprehensive analysis of a diabetic nephropathy cohort of siblings with type 1 diabetes who were discordant for kidney disease points to variants and genes that are potentially causative or protective for diabetic nephropathy. This includes variants in two isoforms of the protein kinase C family not previously linked to diabetic nephropathy, adding support to previous hypotheses that the protein kinase C family members play a role in diabetic nephropathy and might be attractive therapeutic targets.

A remote cis-acting variant at 3q links glomerular NCK1 to diabetic nephropathy

We have previously reported genetic association of a single nucleotide polymorphism (SNP), rs1866813, at 3q locus with increased risk of diabetic nephropathy (DN). The SNP is located approximately 70 kb downstream of a cluster of four genes. This raises a question how the remote noncoding polymorphism affects the risk of DN. In this study, we tested a long-range regulatory potential of this variant by a series of experiments. In a luciferase assay, two alleles of the SNP showed differential effects on the luciferase activity in transfected cells in vitro. Using transgenic zebrafish, we further demonstrated in vivo that two alleles of the SNP differentially regulated GFP expression in zebrafish podocytes. Immunofluorescence staining and Western blotting verified that only Nck1 of the four nearby genes was predominantly expressed in mouse glomeruli as well as in podocytes. Furthermore, genotypes of the SNP rs1866813 were correlated with NCK1 expression in immortalized lymphocytes from diabetic patients. The risk allele was associated with increased NCK1 expression compared to the non-risk allele, consistent with the results of the reporter-based studies. Interestingly, differential expression of glomerular Nck1 between mouse strains carrying the nephropathy-prone 129/Sv allele and nephropathy-resistant C57BL/6 allele was also observed. Our results show that the DN-associated SNP rs1866813 is a remote cis-acting variant differentially regulating glomerular NCK1 expression. This finding implicates an important role for glomerular NCK1 in DN pathogenesis under hyperglycemia.

Association testing of previously reported variants in a large case-control meta-analysis of diabetic nephropathy

We formed the GEnetics of Nephropathy-an International Effort (GENIE) consortium to examine previously reported genetic associations with diabetic nephropathy (DN) in type 1 diabetes. GENIE consists of 6,366 similarly ascertained participants of European ancestry with type 1 diabetes, with and without DN, from the All Ireland-Warren 3-Genetics of Kidneys in Diabetes U.K. and Republic of Ireland (U.K.-R.O.I.) collection and the Finnish Diabetic Nephropathy Study (FinnDiane), combined with reanalyzed data from the Genetics of Kidneys in Diabetes U.S. Study (U.S. GoKinD). We found little evidence for the association of the EPO promoter polymorphism, rs161740, with the combined phenotype of proliferative retinopathy and end-stage renal disease in U.K.-R.O.I. (odds ratio [OR] 1.14, P = 0.19) or FinnDiane (OR 1.06, P = 0.60). However, a fixed-effects meta-analysis that included the previously reported cohorts retained a genome-wide significant association with that phenotype (OR 1.31, P = 2 × 10(-9)). An expanded investigation of the ELMO1 locus and genetic regions reported to be associated with DN in the U.S. GoKinD yielded only nominal statistical significance for these loci. Finally, top candidates identified in a recent meta-analysis failed to reach genome-wide significance. In conclusion, we were unable to replicate most of the previously reported genetic associations for DN, and significance for the EPO promoter association was attenuated.

Association of genetic variants at 3q22 with nephropathy in patients with type 1 diabetes mellitus

Diabetic nephropathy (DN) is the primary cause of morbidity and mortality in patients with type 1 diabetes mellitus (T1DM) and affects about 30% of these patients. We have previously localized a DN locus on chromosome 3q with suggestive linkage in Finnish individuals. Linkage to this region has also been reported earlier by several other groups. To fine map this locus, we conducted a multistage case-control association study in T1DM patients, comprising 1822 cases with nephropathy and 1874 T1DM patients free of nephropathy, from Finland, Iceland, and the British Isles. At the screening stage, we genotyped 3072 tag SNPs, spanning a 28 Mb region, in 234 patients and 215 controls from Finland. SNPs that met the significance threshold of p < 0.01 at this stage were followed up by a series of sample sets. A genetic variant, rs1866813, in the noncoding region at 3q22 was associated with increased risk of DN (overall p = 7.07 x 10(-6), combined odds ratio [OR] of the allele = 1.33). The estimated genotypic ORs of this variant in all Finnish samples suggested a codominant effect, resulting in significant association, with a p value of 4.7 x 10(-5) (OR = 1.38; 95% confidence interval = 1.18-1.62). Additionally, an 11 kb segment flanked by rs62408925 and rs1866813, two strongly correlated variants (r(2) = 0.95), contains three elements highly conserved across multiple species. Independent replication will clarify the role of the associated variants at 3q22 in influencing the risk of DN.