Loci for insulin processing and secretion provide insight into type 2 diabetes risk.

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

RFX6 haploinsufficiency predisposes to diabetes through impaired beta cell function.

AIMS/HYPOTHESIS: Regulatory factor X 6 (RFX6) is crucial for pancreatic endocrine development and differentiation. The RFX6 variant p.His293LeufsTer7 is significantly enriched in the Finnish population, with almost 1:250 individuals as a carrier. Importantly, the FinnGen study indicates a high predisposition for heterozygous carriers to develop type 2 and gestational diabetes. However, the precise mechanism of this predisposition remains unknown. METHODS: To understand the role of this variant in beta cell development and function, we used CRISPR technology to generate allelic series of pluripotent stem cells. We created two isogenic stem cell models: a human embryonic stem cell model; and a patient-derived stem cell model. Both were differentiated into pancreatic islet lineages (stem-cell-derived islets, SC-islets), followed by implantation in immunocompromised NOD-SCID-Gamma mice. RESULTS: Stem cell models of the homozygous variant RFX6-/- predictably failed to generate insulin-secreting pancreatic beta cells, mirroring the phenotype observed in Mitchell-Riley syndrome. Notably, at the pancreatic endocrine stage, there was an upregulation of precursor markers NEUROG3 and SOX9, accompanied by increased apoptosis. Intriguingly, heterozygous RFX6+/- SC-islets exhibited RFX6 haploinsufficiency (54.2% reduction in protein expression), associated with reduced beta cell maturation markers, altered calcium signalling and impaired insulin secretion (62% and 54% reduction in basal and high glucose conditions, respectively). However, RFX6 haploinsufficiency did not have an impact on beta cell number or insulin content. The reduced insulin secretion persisted after in vivo implantation in mice, aligning with the increased risk of variant carriers to develop diabetes. CONCLUSIONS/INTERPRETATION: Our allelic series isogenic SC-islet models represent a powerful tool to elucidate specific aetiologies of diabetes in humans, enabling the sensitive detection of aberrations in both beta cell development and function. We highlight the critical role of RFX6 in augmenting and maintaining the pancreatic progenitor pool, with an endocrine roadblock and increased cell death upon its loss. We demonstrate that RFX6 haploinsufficiency does not affect beta cell number or insulin content but does impair function, predisposing heterozygous carriers of loss-of-function variants to diabetes. DATA AVAILABILITY: Ultra-deep bulk RNA-seq data for pancreatic differentiation stages 3, 5 and 7 of H1 RFX6 genotypes are deposited in the Gene Expression Omnibus database with accession code GSE234289. Original western blot images are deposited at Mendeley ( https://data.mendeley.com/datasets/g75drr3mgw/2 ).

Identification of monogenic variants in more than ten per cent of children without type 1 diabetes-related autoantibodies at diagnosis in the Finnish Pediatric Diabetes Register.

AIMS/HYPOTHESIS: Monogenic forms of diabetes (MODY, neonatal diabetes mellitus and syndromic forms) are rare, and affected individuals may be misclassified and treated suboptimally. The prevalence of type 1 diabetes is high in Finnish children but systematic screening for monogenic diabetes has not been conducted. We assessed the prevalence and clinical manifestations of monogenic diabetes in children initially registered with type 1 diabetes in the Finnish Pediatric Diabetes Register (FPDR) but who had no type 1 diabetes-related autoantibodies (AABs) or had only low-titre islet cell autoantibodies (ICAs) at diagnosis. METHODS: The FPDR, covering approximately 90% of newly diagnosed diabetic individuals aged ≤15 years in Finland starting from 2002, includes data on diabetes-associated HLA genotypes and AAB data (ICA, and autoantibodies against insulin, GAD, islet antigen 2 and zinc transporter 8) at diagnosis. A next generation sequencing gene panel including 42 genes was used to identify monogenic diabetes. We interpreted the variants in HNF1A by using the gene-specific standardised criteria and reported pathogenic and likely pathogenic findings only. For other genes, we also reported variants of unknown significance if an individual's phenotype suggested monogenic diabetes. RESULTS: Out of 6482 participants, we sequenced DNA for 152 (2.3%) testing negative for all AABs and 49 (0.8%) positive only for low-titre ICAs (ICAlow). A monogenic form of diabetes was revealed in 19 (12.5%) of the AAB-negative patients (14 [9.2%] had pathogenic or likely pathogenic variants) and two (4.1%) of the ICAlow group. None had ketoacidosis at diagnosis or carried HLA genotypes conferring high risk for type 1 diabetes. The affected genes were GCK, HNF1A, HNF4A, HNF1B, INS, KCNJ11, RFX6, LMNA and WFS1. A switch from insulin to oral medication was successful in four of five patients with variants in HNF1A, HNF4A or KCNJ11. CONCLUSIONS/INTERPRETATION: More than 10% of AAB-negative children with newly diagnosed diabetes had a genetic finding associated with monogenic diabetes. Because the genetic diagnosis can lead to major changes in treatment, we recommend referring all AAB-negative paediatric patients with diabetes for genetic testing. Low-titre ICAs in the absence of other AABs does not always indicate a diagnosis of type 1 diabetes.

A multigenerational study on phenotypic consequences of the most common causal variant of HNF1A-MODY.

AIMS/HYPOTHESIS: Systematic studies on the phenotypic consequences of variants causal of HNF1A-MODY are rare. Our aim was to assess the phenotype of carriers of a single HNF1A variant and genetic and clinical factors affecting the clinical spectrum. METHODS: We conducted a family-based multigenerational study by comparing heterozygous carriers of the HNF1A p.(Gly292fs) variant with the non-carrier relatives irrespective of diabetes status. During more than two decades, 145 carriers and 131 non-carriers from 12 families participated in the study, and 208 underwent an OGTT at least once. We assessed the polygenic risk score for type 2 diabetes, age at onset of diabetes and measures of body composition, as well as plasma glucose, serum insulin, proinsulin, C-peptide, glucagon and NEFA response during the OGTT. RESULTS: Half of the carriers remained free of diabetes at 23 years, one-third at 33 years and 13% even at 50 years. The median age at diagnosis was 21 years (IQR 17-35). We could not identify clinical factors affecting the age at conversion; sex, BMI, insulin sensitivity or parental carrier status had no significant effect. However, for 1 SD unit increase of a polygenic risk score for type 2 diabetes, the predicted age at diagnosis decreased by 3.2 years. During the OGTT, the carriers had higher levels of plasma glucose and lower levels of serum insulin and C-peptide than the non-carriers. The carriers were also leaner than the non-carriers (by 5.0 kg, p=0.012, and by 2.1 kg/m2 units of BMI, p=2.2 × 10-4, using the first adult measurements) and, possibly as a result of insulin deficiency, demonstrated higher lipolytic activity (with medians of NEFA at fasting 621 vs 441 µmol/l, p=0.0039; at 120 min during an OGTT 117 vs 64 µmol/l, p=3.1 × 10-5). CONCLUSIONS/INTERPRETATION: The most common causal variant of HNF1A-MODY, p.(Gly292fs), presents not only with hyperglycaemia and insulin deficiency, but also with increased lipolysis and markedly lower adult BMI. Serum insulin was more discriminative than C-peptide between carriers and non-carriers. A considerable proportion of carriers develop diabetes after young adulthood. Even among individuals with a monogenic form of diabetes, polygenic risk of diabetes modifies the age at onset of diabetes.

Impaired hepatic lipid synthesis from polyunsaturated fatty acids in TM6SF2 E167K variant carriers with NAFLD.

BACKGROUND: Carriers of the transmembrane 6 superfamily member 2 E167K gene variant (TM6SF2EK/KK) have decreased expression of the TM6SF2 gene and increased risk of NAFLD and NASH. Unlike common 'obese/metabolic' NAFLD, these subjects lack hypertriglyceridemia and have lower risk of cardiovascular disease. In animals, phosphatidylcholine (PC) deficiency results in a similar phenotype. PCs surround the core of VLDL consisting of triglycerides (TGs) and cholesteryl-esters (CEs). We determined the effect of the TM6SF2 E167K on these lipids in the human liver and serum and on hepatic gene expression and studied the effect of TM6SF2 knockdown on hepatocyte handling of these lipids. METHODS: Liver biopsies were taken from subjects characterized with respect to the TM6SF2 genotype, serum and liver lipidome, gene expression and histology. In vitro, after TM6SF2 knockdown in HuH-7 cells, we compared incorporation of different fatty acids into TGs, CEs, and PCs. RESULTS: The TM6SF2EK/KK and TM6SF2EE groups had similar age, gender, BMI and HOMA-IR. Liver TGs and CEs were higher and liver PCs lower in the TM6SF2EK/KK than the TM6SF2EE group (p<0.05). Polyunsaturated fatty acids (PUFA) were deficient in liver and serum TGs and liver PCs but hepatic free fatty acids were relatively enriched in PUFA (p<0.05). Incorporation of PUFA into TGs and PCs in TM6SF2 knockdown hepatocytes was decreased (p<0.05). Hepatic expression of TM6SF2 was decreased in variant carriers, and was co-expressed with genes regulated by PUFAs. CONCLUSIONS: Hepatic lipid synthesis from PUFAs is impaired and could contribute to deficiency in PCs and increased intrahepatic TG in TM6SF2 E167K variant carriers.

Strong influence of variants near MC4R on adiposity in children and adults: a cross-sectional study in Indian population.

Common variants near melanocortin 4 receptor (MC4R) gene are shown to be associated with adiposity but have varied effects in different age groups. Among Indians, studies have shown association of these variants with obesity in adults, but their association in children is yet to be confirmed. We evaluated association of rs17782313 and rs12970134 near MC4R with adiposity and related traits in Indians including 1362 children and 4077 adults (consisting of 2049 diabetic and 2028 nondiabetic adult subjects). Both variants rs17782313 and rs12970134 showed strong association with adiposity measures (weight, body mass index and waist circumference) in children (P-range 7.6 × 10(-5)-3.8 × 10(-12)) and nominal association in nondiabetic adults (P-range 0.05-0.003). Effect sizes on adiposity measures in children (ß range 0.22-0.26 Z-score) were ~3-fold higher compared with adults (ß range 0.06-0.08). The minor alleles of both variants showed borderline association (P-range 0.08-0.04) with risk of type 2 diabetes in adults. Meta-analysis of rs12970134 in >12 000 Indian adults corroborated its association with adiposity (P≤2.2 × 10(-9)), homeostasis model assessment-estimated insulin resistance (P=4.0 × 10(-5)) and type 2 diabetes (P=0.003) with only moderate heterogeneity, suggesting similar effect on adult Indians residing in different geographical regions. In conclusion, the study demonstrates association of variants near MC4R with obesity and related traits in Indian children and adults, with higher impact during childhood.

Capturing the Kidney Transcriptome by Urinary Extracellular Vesicles-From Pre-Analytical Obstacles to Biomarker Research.

Urinary extracellular vesicles (uEV) hold non-invasive RNA biomarkers for genitourinary tract diseases. However, missing knowledge about reference genes and effects of preanalytical choices hinder biomarker studies. We aimed to assess how preanalytical variables (urine storage temperature, isolation workflow) affect diabetic kidney disease (DKD)-linked miRNAs or kidney-linked miRNAs and mRNAs (kidney-RNAs) in uEV isolates and to discover stable reference mRNAs across diverse uEV datasets. We studied nine raw and normalized sequencing datasets including healthy controls and individuals with prostate cancer or type 1 diabetes with or without albuminuria. We focused on kidney-RNAs reviewing literature for DKD-linked miRNAs from kidney tissue, cell culture and uEV/urine experiments. RNAs were analyzed by expression heatmaps, hierarchical clustering and selecting stable mRNAs with normalized counts (>200) and minimal coefficient of variation. Kidney-RNAs were decreased after urine storage at -20 °C vs. -80 °C. Isolation workflows captured kidney-RNAs with different efficiencies. Ultracentrifugation captured DKD -linked miRNAs that separated healthy and diabetic macroalbuminuria groups. Eleven mRNAs were stably expressed across the datasets. Hence, pre-analytical choices had variable effects on kidney-RNAs-analyzing kidney-RNAs complemented global correlation, which could fade differences in some relevant RNAs. Replicating prior DKD-marker results and discovery of candidate reference mRNAs encourages further uEV biomarker studies.

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.

Genome-wide association study and functional characterization identifies candidate genes for insulin-stimulated glucose uptake.

Distinct tissue-specific mechanisms mediate insulin action in fasting and postprandial states. Previous genetic studies have largely focused on insulin resistance in the fasting state, where hepatic insulin action dominates. Here we studied genetic variants influencing insulin levels measured 2 h after a glucose challenge in >55,000 participants from three ancestry groups. We identified ten new loci (P < 5 × 10-8) not previously associated with postchallenge insulin resistance, eight of which were shown to share their genetic architecture with type 2 diabetes in colocalization analyses. We investigated candidate genes at a subset of associated loci in cultured cells and identified nine candidate genes newly implicated in the expression or trafficking of GLUT4, the key glucose transporter in postprandial glucose uptake in muscle and fat. By focusing on postprandial insulin resistance, we highlighted the mechanisms of action at type 2 diabetes loci that are not adequately captured by studies of fasting glycemic traits.