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.

Systematic analysis and functional annotation of variations in the genome of an Indian individual.

Whole genome sequencing of personal genomes has revealed a large repertoire of genomic variations and has provided a rich template for identification of common and rare variants in genomes in addition to understanding the genetic basis of diseases. The widespread application of personal genome sequencing in clinical settings for predictive and preventive medicine has been limited due to the lack of comprehensive computational analysis pipelines. We have used next-generation sequencing technology to sequence the whole genome of a self-declared healthy male of Indian origin. We have generated around 28X of the reference human genome with over 99% coverage. Analysis revealed over 3 million single nucleotide variations and about 490,000 small insertion-deletion events including several novel variants. Using this dataset as a template, we designed a comprehensive computational analysis pipeline for the systematic analysis and annotation of functionally relevant variants in the genome. This study follows a systematic and intuitive data analysis workflow to annotate genome variations and its potential functional effects. Moreover, we integrate predictive analysis of pharmacogenomic traits with emphasis on drugs for which pharmacogenomic testing has been recommended. This study thus provides the template for genome-scale analysis of personal genomes for personalized medicine.

Common variants of SLAMF1 and ITLN1 on 1q21 are associated with type 2 diabetes in Indian population.

Though multiple studies link chromosomal regions 1q21-q23 and 20q13 with type 2 diabetes, fine mapping of these regions is yet to confirm gene(s) explaining the linkages. These candidate regions remain unexplored in Indians, which is a high-risk population for type 2 diabetes. Hypothesizing regulatory regions to have a more important role in complex disorders, we examined association of 207 common variants in proximal promoter and untranslated regions of genes on 1q21-23 and 20q13 with type 2 diabetes in 2115 North Indians. Further, top signals were replicated in an independent group of 2085 North Indians. Variants-rs11265455-SLAMF1 (odds ratios (OR)=1.32, P=1.1 × 10(-3)), rs1062827-F11R (OR=1.36, P=1.7 × 10(-3)) and rs12565932-F11R (OR=1.35, P=1.8 × 10(-3)) were top signals for association with type 2 diabetes whereas rs1333062-ITLN1 (OR=1.28, P=3.4 × 10(-3)) showed strongest association in body mass index-stratified analysis. Replication of these four variants confirmed associations of rs11265455-SLAMF1 (OR=1.27, P=9.1 × 10(-3)) and rs1333062-ITLN1 (OR=1.25, P=1.1 × 10(-3)) with type 2 diabetes. Meta-analysis further corroborated the association of rs11265455-SLAMF1 (OR random effect=1.29, P random effect=3.9 × 10(-5)) and rs1333062-ITLN1 (OR random effect=1.19, P random effect=1.8 × 10(-4)). In conclusion, the study demonstrates that variants of SLAMF1 and ITLN1, both implicated in inflammation, are associated with type 2 diabetes in Indians.

The type 1 diabetes gene TYK2 regulates ß-cell development and its responses to interferon-α.

Type 1 diabetes (T1D) is an autoimmune disease that results in the destruction of insulin producing pancreatic ß-cells. One of the genes associated with T1D is TYK2, which encodes a Janus kinase with critical roles in type-Ι interferon (IFN-Ι) mediated intracellular signalling. To study the role of TYK2 in ß-cell development and response to IFNα, we generated TYK2 knockout human iPSCs and directed them into the pancreatic endocrine lineage. Here we show that loss of TYK2 compromises the emergence of endocrine precursors by regulating KRAS expression, while mature stem cell-islets (SC-islets) function is not affected. In the SC-islets, the loss or inhibition of TYK2 prevents IFNα-induced antigen processing and presentation, including MHC Class Ι and Class ΙΙ expression, enhancing their survival against CD8+ T-cell cytotoxicity. These results identify an unsuspected role for TYK2 in ß-cell development and support TYK2 inhibition in adult ß-cells as a potent therapeutic target to halt T1D progression.

No association of TNFRSF1B variants with type 2 diabetes in Indians of Indo-European origin.

BACKGROUND: There has been no systematic evaluation of the association between genetic variants of type 2 receptor for TNFα (TNFR2) and type 2 diabetes, despite strong biological evidence for the role of this receptor in the pathogenesis of this complex disorder. In view of this, we performed a comprehensive association analysis of TNFRSF1B variants with type 2 diabetes in 4,200 Indo-European subjects from North India. METHODS: The initial phase evaluated association of seven SNPs viz. rs652625, rs496888, rs6697733, rs945439, rs235249, rs17883432 and rs17884213 with type 2 diabetes in 2,115 participants (1,073 type 2 diabetes patients and 1,042 control subjects). Further, we conducted replication analysis of three associated SNPs in 2,085 subjects (1,047 type 2 diabetes patients and 1,038 control subjects). RESULTS: We observed nominal association of rs945439, rs235249 and rs17884213 with type 2 diabetes (P < 0.05) in the initial phase. Haplotype CC of rs945439 and rs235249 conferred increased susceptibility for type 2 diabetes [OR = 1.19 (95%CI 1.03-1.37), P = 0.019/Pperm = 0.076] whereas, TG haplotype of rs235249 and rs17884213 provided protection against type 2 diabetes [OR = 0.83 (95%CI 0.72-0.95, P = 7.2 × 10-3/Pperm = 0.019]. We also observed suggestive association of rs496888 with plasma hsCRP levels [P = 0.042]. However, the association of rs945439, rs235249 and rs17884213 with type 2 diabetes was not replicated in the second study population. Meta-analysis of the two studies also failed to detect any association with type 2 diabetes. CONCLUSIONS: Our two-stage association analysis suggests that TNFRSF1B variants are not the determinants of genetic risk of type 2 diabetes in North Indians.

Association of variants in genes involved in pancreatic ß-cell development and function with type 2 diabetes in North Indians.

Variants in genes involved in pancreatic ß-cell development and function are known to cause monogenic forms of type 2 diabetes and are also associated with complex form. In this study, we studied the genetic association of polymorphisms in such important genes with type 2 diabetes in the high-risk Indians. We genotyped 91 polymorphisms in 19 genes (ABCC8, HNF1A, HNF1B, HNF4A, INS, INSM1, ISL1, KCNJ11, MAFA, MNX1, NEUROD1, NEUROG3, NKX2.2, NKX6.1, PAX4, PAX6, PDX1, USF1 and WFS1) in 2025 unrelated North Indians of Indo-European ethnicity comprising of 1019 diabetic and 1006 non-diabetic subjects. HNF4A promoter P2 polymorphisms rs1884613 and rs2144908, which are in high linkage disequilibrium, showed significant association with type 2 diabetes (odds ratio (OR)=1.37 (95% confidence interval (CI) 1.19-1.57), P=9.4 × 10(-6) for rs1884613 and OR=1.37 (95%CI 1.20-1.57), P=6.0 × 10(-6) for rs2144908), as previously shown in other populations. We observed body mass index-dependent association of these variants with type 2 diabetes in normal-weight/lean subjects. Variants in USF1, ABCC8, ISL1 and KCNJ11 showed nominal association, while haplotypes in these genes were significantly associated. rs3812704 upstream of NEUROG3 significantly increased risk for type 2 diabetes in normal-weight/lean subjects (OR=1.68 (95%CI 1.25-2.24), P=4.9 × 10(-4)). Thus, pancreatic ß-cell development and function genes contribute to susceptibility to type 2 diabetes in North Indians.

Common variants of FTO and the risk of obesity and type 2 diabetes in Indians.

Common variants of fat mass and obesity-associated gene (FTO, fat mass- and obesity-associated gene) have been shown to be associated with obesity and type 2 diabetes in population of European and non-European ethnicity. However, studies in Indian population have provided inconsistent results. Here, we examined association of eight FTO variants (rs1421085, rs8050136, rs9939609, rs9930506, rs1861867, rs9926180, rs2540769 and rs708277) with obesity and type 2 diabetes in 5364 North Indians (2474 type 2 diabetes patients and 2890 non-diabetic controls) in two stages. None of the variants including previously reported intron 1 variants (rs1421085, rs8050136, rs9939609 and rs9930506) showed body mass index (BMI)-dependent/independent association with type 2 diabetes. However, rs1421085, rs8050136 and rs9939609 were associated with obesity status and measures of obesity (BMI, waist circumference and waist-to-hip ratio) in stage 2 and combined study population. Meta-analysis of the two study population results also revealed that rs1421085, rs8050136 and rs9939609 were significantly associated with BMI both under the random- and fixed-effect models (P (random/fixed)=0.02/0.0001, 0.004/0.0006 and 0.01/0.01, respectively). In conclusion, common variants of FTO were associated with obesity, but not with type 2 diabetes in North Indian population.

Urinary extracellular vesicles: Assessment of pre-analytical variables and development of a quality control with focus on transcriptomic biomarker research.

Urinary extracellular vesicles (uEV) are a topical source of non-invasive biomarkers for health and diseases of the urogenital system. However, several challenges have become evident in the standardization of uEV pipelines from collection of urine to biomarker analysis. Here, we studied the effect of pre-analytical variables and developed means of quality control for uEV isolates to be used in transcriptomic biomarker research. We included urine samples from healthy controls and individuals with type 1 or type 2 diabetes and normo-, micro- or macroalbuminuria and isolated uEV by ultracentrifugation. We studied the effect of storage temperature (-20°C vs. -80°C), time (up to 4 years) and storage format (urine or isolated uEV) on quality of uEV by nanoparticle tracking analysis, electron microscopy, Western blotting and qPCR. Urinary EV RNA was compared in terms of quantity, quality, and by mRNA or miRNA sequencing. To study the stability of miRNA levels in samples isolated by different methods, we created and tested a list of miRNAs commonly enriched in uEV isolates. uEV and their transcriptome were preserved in urine or as isolated uEV even after long-term storage at -80°C. However, storage at -20°C degraded particularly the GC-rich part of the transcriptome and EV protein markers. Transcriptome was preserved in RNA samples extracted with and without DNAse, but read distributions still showed some differences in e.g. intergenic and intronic reads. MiRNAs commonly enriched in uEV isolates were stable and concordant between different EV isolation methods. Analysis of never frozen uEV helped to identify surface characteristics of particles by EM. In addition to uEV, qPCR assays demonstrated that uEV isolates commonly contained polyoma viruses. Based on our results, we present recommendations how to store and handle uEV isolates for transcriptomics studies that may help to expedite standardization of the EV biomarker field.

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.

Evaluation of DOK5 as a susceptibility gene for type 2 diabetes and obesity in North Indian population.

BACKGROUND: Type 2 diabetes is a complex metabolic disorder with obesity being a major contributing factor in its development. Susceptibility loci for type 2 diabetes and obesity have been localized on different chromosomal regions by various genome-wide linkage scans. Of these chromosomal regions, 20q13 is one of the strongest linked regions for type 2 diabetes as well as obesity. On 20q13 lies DOK5 that seems to be a strong functional and positional candidate for type 2 diabetes and obesity because of its involvement in insulin signaling and immune responses. Hence, for the first time, we explored DOK5 as a potential type 2 diabetes and obesity susceptibility gene. METHODS: We sequenced 43 subjects for polymorphisms in functionally relevant regions of DOK5. A total of 10 SNPs that included 5 that were identified by sequencing and 5 additional SNPs from NCBI Variation Database were genotyped in 2,115 participants comprising of 1,073 patients with type 2 diabetes and 1,042 controls of Indo-European ethnicity from North India. RESULTS: We identified a novel variant in intron 7 referred to as DK176673. We found nominal association of three SNPs-rs6064099 (OR = 0.75, P = 0.019), rs873079 (OR = 0.76, P = 0.036) and DK176673 (OR = 1.55, P = 0.037) with type 2 diabetes among normal-weight subjects [BMI < 23 kg/m2]. The haplotype GGC harboring rs6068916, rs6064099 and rs873079 showed strong association with type 2 diabetes among normal-weight subjects (OR = 1.37, P/Pperm = 5.8 x 10-3/0.037). Association analysis with obesity revealed that rs6064099 is associated with reduced susceptibility for obesity (OR = 0.48, P = 6.8 x 10-3). Also, haplotype GGC conferred increased susceptibility for obesity (OR = 1.27, P/Pperm = 9.0 x 10-3/0.039). Also, rs6064099 was significantly associated with reduced BMI [median(IQR) = 24.0(20.7-27.1) vs 23.9(20.2-26.8) vs 21.8(19.2-24.7) for GG vs GC vs CC, P = 7.0 x 10-3]. CONCLUSIONS: We identified DOK5 as a novel susceptibility gene for obesity and type 2 diabetes in North Indian subjects. Association of DOK5 variants both with obesity and type 2 diabetes suggests that these variants might modulate type 2 diabetes susceptibility through obesity.

Comparison of urinary extracellular vesicle isolation methods for transcriptomic biomarker research in diabetic kidney disease.

Urinary Extracellular Vesicles (uEV) have emerged as a source for biomarkers of kidney damage, holding potential to replace the conventional invasive techniques including kidney biopsy. However, comprehensive studies characterizing uEV isolation methods with patient samples are rare. Here we compared performance of three established uEV isolation workflows for their subsequent use in transcriptomics analysis for biomarker discovery in diabetic kidney disease. We collected urine samples from individuals with type 1 diabetes with macroalbuminuria and healthy controls. We isolated uEV by Hydrostatic Filtration Dialysis (HFD), ultracentrifugation (UC), and a commercial kit- based isolation method (NG), each with different established urine clearing steps. Purified EVs were analysed by electron microscopy, nanoparticle tracking analysis, and Western blotting. Isolated RNAs were subjected to miRNA and RNA sequencing. HFD and UC samples showed close similarities based on mRNA sequencing data. NG samples had a lower number of reads and different mRNA content compared to HFD or UC. For miRNA sequencing data, satisfactory miRNA counts were obtained by all methods, but miRNA contents differed slightly. This suggests that the isolation workflows enrich specific subpopulations of miRNA-rich uEV preparation components. Our data shows that HFD,UC and the kit-based method are suitable methods to isolate uEV for miRNA-seq. However, only HFD and UC were suitable for mRNA-seq in our settings.