Residual insulin secretion in individuals with type 1 diabetes in Finland: longitudinal and cross-sectional analyses

BACKGROUND

Contrary to the presumption that type 1 diabetes leads to an absolute insulin deficiency, many individuals with type 1 diabetes have circulating C-peptide years after the diagnosis. We studied factors affecting random serum C-peptide concentration in individuals with type 1 diabetes and the association with diabetic complications.

METHODS

Our longitudinal analysis included individuals newly diagnosed with type 1 diabetes from Helsinki University Hospital (Helsinki, Finland) with repeated random serum C-peptide and concomitant glucose measurements from within 3 months of diagnosis and at least once later. The long-term cross-sectional analysis included data from participants from 57 centres in Finland who had type 1 diabetes diagnosed after 5 years of age, initiation of insulin treatment within 1 year from diagnosis, and a C-peptide concentration of less than 1·0 nmol/L (FinnDiane study) and patients with type 1 diabetes from the DIREVA study. We tested the association of random serum C-peptide concentrations and polygenic risk scores with one-way ANOVA, and association of random serum C-peptide concentrations, polygenic risk scores, and clinical factors with logistic regression.

FINDINGS

The longitudinal analysis included 847 participants younger than 16 years and 110 aged 16 years or older. In the longitudinal analysis, age at diagnosis strongly correlated with the decline in C-peptide secretion. The cross-sectional analysis included 3984 participants from FinnDiane and 645 from DIREVA. In the cross-sectional analysis, at a median duration of 21·6 years (IQR 12·5-31·2), 776 (19·4%) of 3984 FinnDiane participants had residual random serum C-peptide secretion (>0·02 nmol/L), which was associated with lower type 1 diabetes polygenic risk compared with participants without random serum C-peptide (p<0·0001). Random serum C-peptide was inversely associated with hypertension, HbA1c, and cholesterol, but also independently with microvascular complications (adjusted OR 0·61 [95% CI 0·38-0·96], p=0·033, for nephropathy; 0·55 [0·34-0·89], p=0·014, for retinopathy).

INTERPRETATION

Although children with multiple autoantibodies and HLA risk genotypes progressed to absolute insulin deficiency rapidly, many adolescents and adults had residual random serum C-peptide decades after the diagnosis. Polygenic risk of type 1 and type 2 diabetes affected residual random serum C-peptide. Even low residual random serum C-peptide concentrations seemed to be associated with a beneficial complications profile.

FUNDING

Folkhälsan Research Foundation; Academy of Finland; University of Helsinki and Helsinki University Hospital; Medical Society of Finland; the Sigrid Juselius Foundation; the "Liv and Hälsa" Society; Novo Nordisk Foundation; and State Research Funding via the Helsinki University Hospital, the Vasa Hospital District, Turku University Hospital, Vasa Central Hospital, Jakobstadsnejdens Heart Foundation, and the Medical Foundation of Vaasa.

Genome-wide mRNA profiling in urinary extracellular vesicles reveals stress gene signature for diabetic kidney disease

Urinary extracellular vesicles (uEV) are a largely unexplored source of kidney-derived mRNAs with potential to serve as a liquid kidney biopsy. We assessed ∼200 uEV mRNA samples from clinical studies by genome-wide sequencing to discover mechanisms and candidate biomarkers of diabetic kidney disease (DKD) in Type 1 diabetes (T1D) with replication in Type 1 and 2 diabetes. Sequencing reproducibly showed >10,000 mRNAs with similarity to kidney transcriptome. T1D DKD groups showed 13 upregulated genes prevalently expressed in proximal tubules, correlated with hyperglycemia and involved in cellular/oxidative stress homeostasis. We used six of them (GPX3, NOX4, MSRB, MSRA, HRSP12, and CRYAB) to construct a transcriptional "stress score" that reflected long-term decline of kidney function and could even identify normoalbuminuric individuals showing early decline. We thus provide workflow and web resource for studying uEV transcriptomes in clinical urine samples and stress-linked DKD markers as potential early non-invasive biomarkers or drug targets.

Distribution and characteristics of newly-defined subgroups of type 2 diabetes in randomised clinical trials: Post hoc cluster assignment analysis of over 12,000 study participants

AIMS

Newly-defined subgroups of type 2 diabetes mellitus (T2DM) have been reported from real-world cohorts but not in detail from randomised clinical trials (RCTs).

METHODS

T2DM participants, uncontrolled on different pre-study therapies (n = 12.738; 82 % Caucasian; 44 % with diabetes duration > 10 years) from 14 RCTs, were assigned to new subgroups according to age at onset of diabetes, HbA1c, BMI, and fasting C-peptide using the nearest centroid approach. Subgroup distribution, characteristics and influencing factors were analysed.

RESULTS

In both, pooled and single RCTs, "mild-obesity related diabetes" predominated (45 %) with mean BMI of 35 kg/m². "Severe insulin-resistant diabetes" was found least often (4.6 %) and prevalence of "mild age-related diabetes" (23.9 %) was mainly influenced by age at onset of diabetes and age cut-offs. Subgroup characteristics were widely comparable to those from real-world cohorts, but all subgroups showed higher frequencies of diabetes-related complications which were associated with longer diabetes duration. A high proportion of "severe insulin-deficient diabetes" (25.4 %) was identified with poor pre-study glycaemic control.

CONCLUSIONS

Classification of RCT participants into newly-defined diabetes subgroups revealed the existence of a heterogeneous population of T2DM. For future RCTs, subgroup-based randomisation of T2DM will better define the target population and relevance of the outcomes by avoiding clinical heterogeneity.

Subgroups of patients with young-onset type 2 diabetes in India reveal insulin deficiency as a major driver

AIM/HYPOTHESIS

Five subgroups were described in European diabetes patients using a data driven machine learning approach on commonly measured variables. We aimed to test the applicability of this phenotyping in Indian individuals with young-onset type 2 diabetes.

METHODS

We applied the European-derived centroids to Indian individuals with type 2 diabetes diagnosed before 45 years of age from the WellGen cohort (n = 1612). We also applied de novo k-means clustering to the WellGen cohort to validate the subgroups. We then compared clinical and metabolic-endocrine characteristics and the complication rates between the subgroups. We also compared characteristics of the WellGen subgroups with those of two young European cohorts, ANDIS (n = 962) and DIREVA (n = 420). Subgroups were also assessed in two other Indian cohorts, Ahmedabad (n = 187) and PHENOEINDY-2 (n = 205).

RESULTS

Both Indian and European young-onset type 2 diabetes patients were predominantly classified into severe insulin-deficient (SIDD) and mild obesity-related (MOD) subgroups, while the severe insulin-resistant (SIRD) and mild age-related (MARD) subgroups were rare. In WellGen, SIDD (53%) was more common than MOD (38%), contrary to findings in Europeans (Swedish 26% vs 68%, Finnish 24% vs 71%, respectively). A higher proportion of SIDD compared with MOD was also seen in Ahmedabad (57% vs 33%) and in PHENOEINDY-2 (67% vs 23%). Both in Indians and Europeans, the SIDD subgroup was characterised by insulin deficiency and hyperglycaemia, MOD by obesity, SIRD by severe insulin resistance and MARD by mild metabolic-endocrine disturbances. In WellGen, nephropathy and retinopathy were more prevalent in SIDD compared with MOD while the latter had higher prevalence of neuropathy.

CONCLUSIONS /INTERPRETATION

Our data identified insulin deficiency as the major driver of type 2 diabetes in young Indians, unlike in young European individuals in whom obesity and insulin resistance predominate. Our results provide useful clues to pathophysiological mechanisms and susceptibility to complications in type 2 diabetes in the young Indian population and suggest a need to review management strategies.

Genome-wide association analyses highlight etiological differences underlying newly defined subtypes of diabetes

Type 2 diabetes has been reproducibly clustered into five subtypes with different disease progression and risk of complications; however, etiological differences are unknown. We used genome-wide association and genetic risk score (GRS) analysis to compare the underlying genetic drivers. Individuals from the Swedish ANDIS (All New Diabetics In Scania) study were compared to individuals without diabetes; the Finnish DIREVA (Diabetes register in Vasa) and Botnia studies were used for replication. We show that subtypes differ with regard to family history of diabetes and association with GRS for diabetes-related traits. The severe insulin-resistant subtype was uniquely associated with GRS for fasting insulin but not with variants in the TCF7L2 locus or GRS reflecting insulin secretion. Further, an SNP (rs10824307) near LRMDA was uniquely associated with mild obesity-related diabetes. Therefore, we conclude that the subtypes have partially distinct genetic backgrounds indicating etiological differences.

First Genome-Wide Association Study of Latent Autoimmune Diabetes in Adults Reveals Novel Insights Linking Immune and Metabolic Diabetes

OBJECTIVE

Latent autoimmune diabetes in adults (LADA) shares clinical features with both type 1 and type 2 diabetes; however, there is ongoing debate regarding the precise definition of LADA. Understanding its genetic basis is one potential strategy to gain insight into appropriate classification of this diabetes subtype.

RESEARCH DESIGN AND METHODS

We performed the first genome-wide association study of LADA in case subjects of European ancestry versus population control subjects (n = 2,634 vs. 5,947) and compared against both case subjects with type 1 diabetes (n = 2,454 vs. 968) and type 2 diabetes (n = 2,779 vs. 10,396).

RESULTS

The leading genetic signals were principally shared with type 1 diabetes, although we observed positive genetic correlations genome-wide with both type 1 and type 2 diabetes. Additionally, we observed a novel independent signal at the known type 1 diabetes locus harboring PFKFB3, encoding a regulator of glycolysis and insulin signaling in type 2 diabetes and inflammation and autophagy in autoimmune disease, as well as an attenuation of key type 1-associated HLA haplotype frequencies in LADA, suggesting that these are factors that distinguish childhood-onset type 1 diabetes from adult autoimmune diabetes.

CONCLUSIONS

Our results support the need for further investigations of the genetic factors that distinguish forms of autoimmune diabetes as well as more precise classification strategies.

A Genome-Wide Association Study of Diabetic Kidney Disease in Subjects With Type 2 Diabetes

Identification of sequence variants robustly associated with predisposition to diabetic kidney disease (DKD) has the potential to provide insights into the pathophysiological mechanisms responsible. We conducted a genome-wide association study (GWAS) of DKD in type 2 diabetes (T2D) using eight complementary dichotomous and quantitative DKD phenotypes: the principal dichotomous analysis involved 5,717 T2D subjects, 3,345 with DKD. Promising association signals were evaluated in up to 26,827 subjects with T2D (12,710 with DKD). A combined T1D+T2D GWAS was performed using complementary data available for subjects with T1D, which, with replication samples, involved up to 40,340 subjects with diabetes (18,582 with DKD). Analysis of specific DKD phenotypes identified a novel signal near GABRR1 (rs9942471, P = 4.5 × 10^-8) associated with microalbuminuria in European T2D case subjects. However, no replication of this signal was observed in Asian subjects with T2D or in the equivalent T1D analysis. There was only limited support, in this substantially enlarged analysis, for association at previously reported DKD signals, except for those at UMOD and PRKAG2, both associated with estimated glomerular filtration rate. We conclude that, despite challenges in addressing phenotypic heterogeneity, access to increased sample sizes will continue to provide more robust inference regarding risk variant discovery for DKD.

Novel subgroups of adult-onset diabetes and their association with outcomes: a data-driven cluster analysis of six variables

BACKGROUND

Diabetes is presently classified into two main forms, type 1 and type 2 diabetes, but type 2 diabetes in particular is highly heterogeneous. A refined classification could provide a powerful tool to individualise treatment regimens and identify individuals with increased risk of complications at diagnosis.

METHODS

We did data-driven cluster analysis (k-means and hierarchical clustering) in patients with newly diagnosed diabetes (n=8980) from the Swedish All New Diabetics in Scania cohort. Clusters were based on six variables (glutamate decarboxylase antibodies, age at diagnosis, BMI, HbA_1c, and homoeostatic model assessment 2 estimates of β-cell function and insulin resistance), and were related to prospective data from patient records on development of complications and prescription of medication. Replication was done in three independent cohorts: the Scania Diabetes Registry (n=1466), All New Diabetics in Uppsala (n=844), and Diabetes Registry Vaasa (n=3485). Cox regression and logistic regression were used to compare time to medication, time to reaching the treatment goal, and risk of diabetic complications and genetic associations.

FINDINGS

We identified five replicable clusters of patients with diabetes, which had significantly different patient characteristics and risk of diabetic complications. In particular, individuals in cluster 3 (most resistant to insulin) had significantly higher risk of diabetic kidney disease than individuals in clusters 4 and 5, but had been prescribed similar diabetes treatment. Cluster 2 (insulin deficient) had the highest risk of retinopathy. In support of the clustering, genetic associations in the clusters differed from those seen in traditional type 2 diabetes.

INTERPRETATION

We stratified patients into five subgroups with differing disease progression and risk of diabetic complications. This new substratification might eventually help to tailor and target early treatment to patients who would benefit most, thereby representing a first step towards precision medicine in diabetes.

FUNDING

Swedish Research Council, European Research Council, Vinnova, Academy of Finland, Novo Nordisk Foundation, Scania University Hospital, Sigrid Juselius Foundation, Innovative Medicines Initiative 2 Joint Undertaking, Vasa Hospital district, Jakobstadsnejden Heart Foundation, Folkhälsan Research Foundation, Ollqvist Foundation, and Swedish Foundation for Strategic Research.

Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes

We aggregated coding variant data for 81,412 type 2 diabetes cases and 370,832 controls of diverse ancestry, identifying 40 coding variant association signals (P < 2.2 × 10-7); of these, 16 map outside known risk-associated loci. We make two important observations. First, only five of these signals are driven by low-frequency variants: even for these, effect sizes are modest (odds ratio ≤1.29). Second, when we used large-scale genome-wide association data to fine-map the associated variants in their regional context, accounting for the global enrichment of complex trait associations in coding sequence, compelling evidence for coding variant causality was obtained for only 16 signals. At 13 others, the associated coding variants clearly represent 'false leads' with potential to generate erroneous mechanistic inference. Coding variant associations offer a direct route to biological insight for complex diseases and identification of validated therapeutic targets; however, appropriate mechanistic inference requires careful specification of their causal contribution to disease predisposition.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Type 2 diabetes susceptibility gene variants predispose to adult-onset autoimmune diabetes

AIMS/HYPOTHESIS

Latent autoimmune diabetes in adults (LADA) is phenotypically a hybrid of type 1 and type 2 diabetes. Genetically LADA is poorly characterised but does share genetic predisposition with type 1 diabetes. We aimed to improve the genetic characterisation of LADA and hypothesised that type 2 diabetes-associated gene variants also predispose to LADA, and that the associations would be strongest in LADA patients with low levels of GAD autoantibodies (GADA).

METHODS

We assessed 41 type 2 diabetes-associated gene variants in Finnish (phase I) and Swedish (phase II) patients with LADA (n = 911) or type 1 diabetes (n = 406), all diagnosed after the age of 35 years, as well as in non-diabetic control individuals 40 years or older (n = 4,002).

RESULTS

Variants in the ZMIZ1 (rs12571751, p = 4.1 × 10(-5)) and TCF7L2 (rs7903146, p = 5.8 × 10(-4)) loci were strongly associated with LADA. Variants in the KCNQ1 (rs2237895, p = 0.0012), HHEX (rs1111875, p = 0.0024 in Finns) and MTNR1B (rs10830963, p = 0.0039) loci showed the strongest association in patients with low GADA, supporting the hypothesis that the disease in these patients is more like type 2 diabetes. In contrast, variants in the KLHDC5 (rs10842994, p = 9.5 × 10(-4) in Finns), TP53INP1 (rs896854, p = 0.005), CDKAL1 (rs7756992, p = 7.0 × 10(-4); rs7754840, p = 8.8 × 10(-4)) and PROX1 (rs340874, p = 0.003) loci showed the strongest association in patients with high GADA. For type 1 diabetes, a strong association was seen for MTNR1B (rs10830963, p = 3.2 × 10(-6)) and HNF1A (rs2650000, p = 0.0012).

CONCLUSIONS/INTERPRETATION

LADA and adult-onset type 1 diabetes share genetic risk variants with type 2 diabetes, supporting the idea of a hybrid form of diabetes and distinguishing them from patients with classical young-onset type 1 diabetes.