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.

The Genetic Landscape of Renal Complications in Type 1 Diabetes.

Diabetes is the leading cause of ESRD. Despite evidence for a substantial heritability of diabetic kidney disease, efforts to identify genetic susceptibility variants have had limited success. We extended previous efforts in three dimensions, examining a more comprehensive set of genetic variants in larger numbers of subjects with type 1 diabetes characterized for a wider range of cross-sectional diabetic kidney disease phenotypes. In 2843 subjects, we estimated that the heritability of diabetic kidney disease was 35% (P=6.4×10-3). Genome-wide association analysis and replication in 12,540 individuals identified no single variants reaching stringent levels of significance and, despite excellent power, provided little independent confirmation of previously published associated variants. Whole-exome sequencing in 997 subjects failed to identify any large-effect coding alleles of lower frequency influencing the risk of diabetic kidney disease. However, sets of alleles increasing body mass index (P=2.2×10-5) and the risk of type 2 diabetes (P=6.1×10-4) associated with the risk of diabetic kidney disease. We also found genome-wide genetic correlation between diabetic kidney disease and failure at smoking cessation (P=1.1×10-4). Pathway analysis implicated ascorbate and aldarate metabolism (P=9.0×10-6), and pentose and glucuronate interconversions (P=3.0×10-6) in pathogenesis of diabetic kidney disease. These data provide further evidence for the role of genetic factors influencing diabetic kidney disease in those with type 1 diabetes and highlight some key pathways that may be responsible. Altogether these results reveal important biology behind the major cause of kidney disease.

New susceptibility loci associated with kidney disease in type 1 diabetes.

Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ~2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P = 1.2 × 10(-8)) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P = 2.0 × 10(-9)). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-ß1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P = 2.1 × 10(-7)), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.

Genome-wide association study of urinary albumin excretion rate in patients with type 1 diabetes.

AIMS/HYPOTHESIS: An abnormal urinary albumin excretion rate (AER) is often the first clinically detectable manifestation of diabetic nephropathy. Our aim was to estimate the heritability and to detect genetic variation associated with elevated AER in patients with type 1 diabetes. METHODS: The discovery phase genome-wide association study (GWAS) included 1,925 patients with type 1 diabetes and with data on 24 h AER. AER was analysed as a continuous trait and the analysis was stratified by the use of antihypertensive medication. Signals with a p value <10(-4) were followed up in 3,750 additional patients with type 1 diabetes from seven studies. RESULTS: The narrow-sense heritability, captured with our genotyping platform, was estimated to explain 27.3% of the total AER variability, and 37.6% after adjustment for covariates. In the discovery stage, five single nucleotide polymorphisms in the GLRA3 gene were strongly associated with albuminuria (p < 5 × 10(-8)). In the replication group, a nominally significant association (p = 0.035) was observed between albuminuria and rs1564939 in GLRA3, but this was in the opposite direction. Sequencing of the surrounding genetic region in 48 Finnish and 48 UK individuals supported the possibility that population-specific rare variants contribute to the synthetic association observed at the common variants in GLRA3. The strongest replication (p = 0.026) was obtained for rs2410601 between the PSD3 and SH2D4A genes. Pathway analysis highlighted natural killer cell mediated immunity processes. CONCLUSIONS/INTERPRETATION: This study suggests novel pathways and molecular mechanisms for the pathogenesis of albuminuria in type 1 diabetes.

Novel genetic susceptibility loci for diabetic end-stage renal disease identified through robust naive Bayes classification.

AIMS/HYPOTHESIS: Diabetic nephropathy is a major diabetic complication, and diabetes is the leading cause of end-stage renal disease (ESRD). Family studies suggest a hereditary component for diabetic nephropathy. However, only a few genes have been associated with diabetic nephropathy or ESRD in diabetic patients. Our aim was to detect novel genetic variants associated with diabetic nephropathy and ESRD. METHODS: We exploited a novel algorithm, 'Bag of Naive Bayes', whose marker selection strategy is complementary to that of conventional genome-wide association models based on univariate association tests. The analysis was performed on a genome-wide association study of 3,464 patients with type 1 diabetes from the Finnish Diabetic Nephropathy (FinnDiane) Study and subsequently replicated with 4,263 type 1 diabetes patients from the Steno Diabetes Centre, the All Ireland-Warren 3-Genetics of Kidneys in Diabetes UK collection (UK-Republic of Ireland) and the Genetics of Kidneys in Diabetes US Study (GoKinD US). RESULTS: Five genetic loci (WNT4/ZBTB40-rs12137135, RGMA/MCTP2-rs17709344, MAPRE1P2-rs1670754, SEMA6D/SLC24A5-rs12917114 and SIK1-rs2838302) were associated with ESRD in the FinnDiane study. An association between ESRD and rs17709344, tagging the previously identified rs12437854 and located between the RGMA and MCTP2 genes, was replicated in independent case-control cohorts. rs12917114 near SEMA6D was associated with ESRD in the replication cohorts under the genotypic model (p < 0.05), and rs12137135 upstream of WNT4 was associated with ESRD in Steno. CONCLUSIONS/INTERPRETATION: This study supports the previously identified findings on the RGMA/MCTP2 region and suggests novel susceptibility loci for ESRD. This highlights the importance of applying complementary statistical methods to detect novel genetic variants in diabetic nephropathy and, in general, in complex diseases.

A polymorphism in the angiotensin II type 1 receptor gene has different effects on the risk of diabetic nephropathy in men and women.

BACKGROUND: The etiology of diabetic nephropathy depends partly on genetic factors. Elevated systemic and intraglomerular blood pressure and glomerular filtration rate, partly regulated by the renin-angiotensin system, increase the risk of diabetic nephropathy. METHODS: The present case-control study investigated the association of the rs5186 polymorphism, in the angiotensin II type 1 receptor gene (AGTR1), with diabetic nephropathy. The study included 3561 patients with type 1 diabetes from Denmark, Finland, France and Sweden. Microalbuminuria was defined as albumin excretion rate (AER) ≥20 to <200 µg/min or albumin concentration ≥30 to <300 mg/l (n=707), macroalbuminuria was defined as AER≥200 µg/min or ≥300 mg/l (n=1546), and patients with renal replacement therapy were also included in this group. The controls had >15 years diabetes duration, AER <20 µg/min or <30 mg/l, and no antihypertensive treatment (n=1308). RESULTS: AA genotype of the rs5186 polymorphism significantly increased the risk of diabetic nephropathy in male patients, OR=1.27 (95% CI=1.02-1.58), P=0.03, adjusted for age at diabetes onset, HbA1c, diabetes duration, smoking and country of origin. Among the women, there were no significant associations between rs5186 and diabetic nephropathy, OR=0.89 (0.71-1.11), P=0.30. CONCLUSION: We conclude that the AGTR1 gene may be associated with increased risk of diabetic nephropathy in men with type 1 diabetes.

1H NMR metabonomics approach to the disease continuum of diabetic complications and premature death.

Subtle metabolic changes precede and accompany chronic vascular complications, which are the primary causes of premature death in diabetes. To obtain a multimetabolite characterization of these high-risk individuals, we measured proton nuclear magnetic resonance (1H NMR) data from the serum of 613 patients with type I diabetes and a diverse spread of complications. We developed a new metabonomics framework to visualize and interpret the data and to link the metabolic profiles to the underlying diagnostic and biochemical variables. Our results indicate complex interactions between diabetic kidney disease, insulin resistance and the metabolic syndrome. We illustrate how a single 1H NMR protocol is able to identify the polydiagnostic metabolite manifold of type I diabetes and how its alterations translate to clinical phenotypes, clustering of micro- and macrovascular complications, and mortality during several years of follow-up. This work demonstrates the diffuse nature of complex vascular diseases and the limitations of single diagnostic biomarkers. However, it also promises cost-effective solutions through high-throughput analytics and advanced computational methods, as applied here in a case that is representative of the real clinical situation.

Genetic Determinants of Glycated Hemoglobin in Type 1 Diabetes.

Glycated hemoglobin (HbA1c) is an important measure of glycemia in diabetes. HbA1c is influenced by environmental and genetic factors both in people with and in people without diabetes. We performed a genome-wide association study (GWAS) for HbA1c in a Finnish type 1 diabetes (T1D) cohort, FinnDiane. Top results were examined for replication in T1D cohorts DCCT/EDIC, WESDR, CACTI, EDC, and RASS, and a meta-analysis was performed. Three SNPs in high linkage disequilibrium on chromosome 13 near relaxin family peptide receptor 2 (RXFP2) were associated with HbA1c in FinnDiane at genome-wide significance (P < 5 × 10-8). The minor alleles of rs2085277 and rs1360072 were associated with higher HbA1c also in the meta-analysis with RASS (P < 5 × 10-8), where these variants had minor allele frequencies ≥1%. Furthermore, these SNPs were associated with HbA1c in an East Asian population without diabetes (P ≤ 0.013). A weighted genetic risk score created from 55 HbA1c-associated variants from the literature was associated with HbA1c in FinnDiane but explained only a small amount of variation. Understanding the genetic basis of glycemic control and HbA1c may lead to better prevention of diabetes complications.

A quality assessment survey of SNP genotyping laboratories.

To survey the quality of SNP genotyping, a joint Nordic quality assessment (QA) round was organized between 11 laboratories in the Nordic and Baltic countries. The QA round involved blinded genotyping of 47 DNA samples for 18 or six randomly selected SNPs. The methods used by the participating laboratories included all major platforms for small- to medium-size SNP genotyping. The laboratories used their standard procedures for SNP assay design, genotyping, and quality control. Based on the joint results from all laboratories, a consensus genotype for each DNA sample and SNP was determined by the coordinator of the survey, and the results from each laboratory were compared to this genotype. The overall genotyping accuracy achieved in the survey was excellent. Six laboratories delivered genotype data that were in full agreement with the consensus genotype. The average accuracy per SNP varied from 99.1 to 100% between the laboratories, and it was frequently 100% for the majority of the assays for which SNP genotypes were reported. Lessons from the survey are that special attention should be given to the quality of the DNA samples prior to genotyping, and that a conservative approach for calling the genotypes should be used to achieve a high accuracy.

High-throughput pedigree drawing.

Family trees have long been a valuable visual tool for geneticists in identifying clusters of inherited traits and genotypes. As more data are collected, drawing the graphs by hand becomes impractical and, for this reason, we have developed the pedigree software CraneFoot. It can process any family graph with minimal computational cost by making a pedigree transformation that enables the use of a linear node positioning algorithm. The program is designed for automated drawing to printed media and efficient visual classification of genetically interesting families from large data sets. It also incorporates a robust pedigree topology check with detailed error messages.

Fine mapping of the multiple sclerosis susceptibility locus on 5p14-p12.

Linkage analyses have identified four major MS susceptibility loci in Finns. Here we have fine mapped the region on chromosome 5p in 28 Finnish MS families. Marker D5S416 provided the highest pairwise LOD score, and multipoint and haplotype analyses restrict the critical region to about 5.3 Mb on 5p15 between markers D5S1987 and D5S416. Ascertaining for HLA type and geographical origin indicated that families with and without the HLA DR15 risk haplotype, as well as families within and outside an internal high-risk region, contributed to the linkage to 5p, implying the general significance for this locus in Finnish MS families.

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.

Novel susceptibility locus at 22q11 for diabetic nephropathy in type 1 diabetes.

BACKGROUND: Diabetic nephropathy (DN) affects about 30% of patients with type 1 diabetes (T1D) and contributes to serious morbidity and mortality. So far only the 3q21-q25 region has repeatedly been indicated as a susceptibility region for DN. The aim of this study was to search for new DN susceptibility loci in Finnish, Danish and French T1D families. METHODS AND RESULTS: We performed a genome-wide linkage study using 384 microsatellite markers. A total of 175 T1D families were studied, of which 94 originated from Finland, 46 from Denmark and 35 from France. The whole sample set consisted of 556 individuals including 42 sib-pairs concordant and 84 sib-pairs discordant for DN. Two-point and multi-point non-parametric linkage analyses were performed using the Analyze package and the MERLIN software. A novel DN locus on 22q11 was identified in the joint analysis of the Finnish, Danish and French families by genome-wide multipoint non-parametric linkage analysis using the Kong and Cox linear model (NPL(pairs) LOD score 3.58). Nominal or suggestive evidence of linkage to this locus was also detected when the three populations were analyzed separately. Suggestive evidence of linkage was found to six additional loci in the Finnish and French sample sets. CONCLUSIONS: This study identified a novel DN locus at chromosome 22q11 with significant evidence of linkage to DN. Our results suggest that this locus may be of importance in European populations. In addition, this study supports previously indicated DN loci on 3q21-q25 and 19q13.

Acute hyperglycaemia induces an inflammatory response in young patients with type 1 diabetes.

BACKGROUND: Patients with type 1 diabetes (T1D) are at a substantially increased risk of cardiovascular disease. Stress-induced hyperglycaemia in turn is shown to worsen the prognosis of patients suffering from an acute myocardial infarction. However, the mechanisms behind these findings are incompletely known. AIM: To investigate whether markers of chronic inflammation, and oxidative stress respond to acute hyperglycaemia in patients with T1D. METHODS: The plasma glucose concentration was rapidly raised from 5 to 15 mmol/L in 35 males (22 men with T1D and 13 age-matched non-diabetic volunteers) and maintained for 2 h. All participants were young non-smokers without any signs of diabetic or other complications. Markers of chronic inflammation, and oxidative stress were analysed in serum/plasma samples drawn at base-line and after 120 min of hyperglycaemia. RESULTS: Compared to normoglycaemia, acute hyperglycaemia increased the interleukin (IL)-6 concentrations by 39% in patients with T1D (P<0.01) and 26% in healthy volunteers (P<0.05). During hyperglycaemia the superoxide dismutase concentration was increased by 17% in the healthy volunteers (P<0.01) and 5% in the patients with type 1 diabetes (P=NS). The increase in tumour necrosis factor (TNF)-alpha was larger in patients with type 1 diabetes than in non-diabetic volunteers (35% versus -10%, P<0.05). CONCLUSIONS: This study shows that acute hyperglycaemia induces an inflammatory response in patients with type 1 diabetes.

Metabolic phenotypes, vascular complications, and premature deaths in a population of 4,197 patients with type 1 diabetes.

OBJECTIVE: Poor glycemic control, elevated triglycerides, and albuminuria are associated with vascular complications in diabetes. However, few studies have investigated combined associations between metabolic markers, diabetic kidney disease, retinopathy, hypertension, obesity, and mortality. Here, the goal was to reveal previously undetected association patterns between clinical diagnoses and biochemistry in the FinnDiane dataset. RESEARCH DESIGN AND METHODS: At baseline, clinical records, serum, and 24-h urine samples of 2,173 men and 2,024 women with type 1 diabetes were collected. The data were analyzed by the self-organizing map, which is an unsupervised pattern recognition algorithm that produces a two-dimensional layout of the patients based on their multivariate biochemical profiles. At follow-up, the results were compared against all-cause mortality during 6.5 years (295 deaths). RESULTS: The highest mortality was associated with advanced kidney disease. Other risk factors included 1) a profile of insulin resistance, abdominal obesity, high cholesterol, triglycerides, and low HDL(2) cholesterol, and 2) high adiponectin and high LDL cholesterol for older patients. The highest population-adjusted risk of death was 10.1-fold (95% CI 7.3-13.1) for men and 10.7-fold (7.9-13.7) for women. Nonsignificant risk was observed for a profile with good glycemic control and high HDL(2) cholesterol and for a low cholesterol profile with a short diabetes duration. CONCLUSIONS: The self-organizing map analysis enabled detailed risk estimates, described the associations between known risk factors and complications, and uncovered statistical patterns difficult to detect by classical methods. The results also suggest that diabetes per se, without an adverse metabolic phenotype, does not contribute to increased mortality.

Diagnosing diabetic nephropathy by 1H NMR metabonomics of serum.

OBJECT: The most severe complication of type 1 diabetes (T1DM) is diabetic nephropathy. It is associated with a high risk of cardiovascular complications and premature death and requires early detection to be efficiently treated. The clinical practice to diagnose diabetic nephropathy is also a non-optimal and tedious set up based on albumin excretion rate in multiple overnight or 24h urine samples. Conversely, in this study, these independent diagnostic data are used to provide a realistic testing case for applying (1)H NMR metabonomics of serum in a diagnostic fashion. MATERIALS AND METHODS: 182 T1DM and 21 non-diabetic (non-T1DM) individuals were studied. The (1)H NMR of serum at 500 MHz was targeted at two molecular windows: lipoprotein lipids and low-molecular-weight metabolites. RESULTS: T1DM and non-T1DM individuals were exclusively separated by (1)H NMR. For diabetic nephropathy diagnosis in the T1DM patients, (1)H NMR data (and clinical biochemistry data) gave a sensitivity of 87.1% (83.9%) and a specificity of 87.7% (95.9%). The predictive values of positive and negative tests were 89.0% (95.5%) and 83.6% (79.2%), respectively. CONCLUSIONS: (1)H NMR metabonomics clearly distinguishes metabolic characteristics of T1DM and appears approximately as good a means to diagnose diabetic nephropathy from serum as an advanced set of biochemical variables.

Association analysis of the AIRE and insulin genes in Finnish type 1 diabetic patients.

Mutations in the autoimmune regulator (AIRE) gene cause a recessive Mendelian disorder autoimmune polyendocrinopathy syndrome type 1 (APS-1 or autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy). APS-1 patients develop multiorgan autoimmune diseases including type 1 diabetes (prevalence 12%). The AIRE protein controls the central tolerance induction in the thymus by regulating the expression levels of tissue-specific peripheral antigens, such as insulin. We hypothesized that the insulin gene (INS) polymorphisms together with the AIRE variations may predispose individuals to diabetes. The role of the AIRE gene was tested both independently and on the condition of the INS risk genotype in the Finnish type 1 diabetes sample. A total of 733 type 1 diabetic cases and 735 age- and sex-matched healthy controls were used in the analysis. Five common single nucleotide polymorphisms (SNPs) in the AIRE gene were selected from the public database (dbSNP). The -23HphI polymorphism was used as a surrogate marker for the INS gene promoter repeat. The five genotyped SNPs in the AIRE gene showed no evidence of association with type 1 diabetes. As expected, the INS gene polymorphism -23HphI was significantly associated with susceptibility to type 1 diabetes (P=6.8 x 10(-12), chi(2) test). When the subclass of patients carrying the homozygote genotype of the INS gene was used in the analysis, the AIRE polymorphisms showed no association with the disease. In conclusion, the AIRE gene does not seem to contribute to disease susceptibility in Finnish type 1 diabetic patients, whereas the insulin gene represents a notable risk factor for disease in this population.