Epigenome-wide association study of DNA methylation in maternal blood leukocytes with BMI in pregnancy and gestational weight gain

OBJECTIVES

We aimed to discover CpG sites with differential DNA methylation in peripheral blood leukocytes associated with body mass index (BMI) in pregnancy and gestational weight gain (GWG) in women of European and South Asian ancestry. Furthermore, we aimed to investigate how the identified sites were associated with methylation quantitative trait loci, gene ontology, and cardiometabolic parameters.

METHODS

In the Epigenetics in pregnancy (EPIPREG) sample we quantified maternal DNA methylation in peripheral blood leukocytes in gestational week 28 with Illumina's MethylationEPIC BeadChip. In women with European (n = 303) and South Asian (n = 164) ancestry, we performed an epigenome-wide association study of BMI in gestational week 28 and GWG between gestational weeks 15 and 28 using a meta-analysis approach. Replication was performed in the Norwegian Mother, Father, and Child Cohort Study, the Study of Assisted Reproductive Technologies (MoBa-START) (n = 877, mainly European/Norwegian).

RESULTS

We identified one CpG site significantly associated with GWG (p 5.8 × 10-8) and five CpG sites associated with BMI at gestational week 28 (p from 4.0 × 10-8 to 2.1 × 10-10). Of these, we were able to replicate three in MoBa-START; cg02786370, cg19758958 and cg10472537. Two sites are located in genes previously associated with blood pressure and BMI. DNA methylation at the three replicated CpG sites were associated with levels of blood pressure, lipids and glucose in EPIPREG (p from 1.2 × 10-8 to 0.04).

CONCLUSIONS

We identified five CpG sites associated with BMI at gestational week 28, and one with GWG. Three of the sites were replicated in an independent cohort. Several genetic variants were associated with DNA methylation at cg02786379 and cg16733643 suggesting a genetic component influencing differential methylation. The identified CpG sites were associated with cardiometabolic traits.

CLINICALTRIALS

GOV REGISTRATION NO

Not applicable.

Human Genetic Variation at rs10071329 Correlates With Adiposity-Related Traits, Modulates PPARGC1B Expression, and Alters Brown Adipocyte Function

Human genetic variation in PPARGC1B has been associated with adiposity, but the genetic variants that affect PPARGC1B expression have not been experimentally determined. Here, guided by previous observational data, we used clustered regularly interspaced short palindromic repeats/CRISPR associated protein 9 (CRISPR/Cas9) to scarlessly edit the alleles of the candidate causal genetic variant rs10071329 in a human brown adipocyte cell line. Switching the rs10071329 genotype from A/A to G/G enhanced PPARGC1B expression throughout the adipogenic differentiation, identifying rs10071329 as a cis-expression quantitative trait loci (eQTL). The higher PPARGC1B expression in G/G cells coincided with greater accumulation of triglycerides and higher expression of mitochondria-encoded genes, but without significant effects on adipogenic marker expression. Furthermore, G/G cells had improved basal- and norepinephrine-stimulated mitochondrial respiration, possibly relating to enhanced mitochondrial gene expression. The G/G cells also exhibited increased norepinephrine-stimulated glycerol release, indicating improved lipolysis. Altogether, our results showed that rs10071329 is a cis-eQTL, with the G/G genotype conferring enhanced PPARGC1B expression, with consequent improved mitochondrial function and response to norepinephrine in brown adipocytes. This genetic variant, and as yet undetermined eQTLs, at PPARGC1B could prove useful in genotype-based precision medicine for obesity treatment.

ARTICLE HIGHLIGHTS

DNA methylation risk score for type 2 diabetes is associated with gestational diabetes

BACKGROUND

Gestational diabetes mellitus (GDM) and type 2 diabetes mellitus (T2DM) share many pathophysiological factors including genetics, but whether epigenetic marks are shared is unknown. We aimed to test whether a DNA methylation risk score (MRS) for T2DM was associated with GDM across ancestry and GDM criteria.

METHODS

In two independent pregnancy cohorts, EPIPREG (n = 480) and EPIDG (n = 32), DNA methylation in peripheral blood leukocytes was measured at a gestational age of 28 ± 2. We constructed an MRS in EPIPREG and EPIDG based on CpG hits from a published epigenome-wide association study (EWAS) of T2DM.

RESULTS

With mixed models logistic regression of EPIPREG and EPIDG, MRS for T2DM was associated with GDM: odd ratio (OR)[95% CI]: 1.3 [1.1-1.8], P = 0.002 for the unadjusted model, and 1.4 [1.1-1.7], P = 0.00014 for a model adjusted by age, pre-pregnant BMI, family history of diabetes and smoking status. Also, we found 6 CpGs through a meta-analysis (cg14020176, cg22650271, cg14870271, cg27243685, cg06378491, cg25130381) associated with GDM, and some of their methylation quantitative loci (mQTLs) were related to T2DM and GDM.

CONCLUSION

For the first time, we show that DNA methylation marks for T2DM are also associated with GDM, suggesting shared epigenetic mechanisms between GDM and T2DM.

Digital lifestyle treatment improves long-term metabolic control in type 2 diabetes with different effects in pathophysiological and genetic subgroups

To address the unmet need for scalable solutions for lifestyle treatment, we developed a new digital method to promote behavioral change. Here we report that patients with type-2 diabetes in Sweden (n = 331) exposed to the intervention have significantly improved HbA1c during a median follow-up of 1038 days (4 mmol/mol compared with matched controls; P = 0.009). This is paralleled by reduced body weight, ameliorated insulin secretion, increased physical activity, and cognitive eating restraints. Participants with high BMI and insulin resistance have an even larger response, as have non-risk allele carriers for the FTO gene. The findings open a new avenue for scalable lifestyle management with sustained efficacy and highlight a previously unrecognized opportunity for digital precision treatment based on genetics and individual pathophysiology. ClinicalTrials.gov NCT04624321.

Polygenic scores of diabetes-related traits in subgroups of type 2 diabetes in India: a cohort study

Background

A machine-learning approach identified five subgroups of diabetes in Europeans which included severe autoimmune diabetes (SAID), severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD), mild obesity-related diabetes (MOD) and mild age-related diabetes (MARD) with partially distinct genetic aetiologies. We previously validated four of the non-autoimmune subgroups in people with young-onset type 2 diabetes (T2D) from the Indian WellGen study. Here, we aimed to apply European-derived centroids and genetic risk scores (GRSs) to the unselected (for age) WellGen to test their applicability and investigate the genetic aetiology of the Indian T2D subgroups.

Methods

We applied European derived centroids and GRSs to T2D participants of Indian ancestry (WellGen, n = 2217, 821 genotyped) and compared them with normal glucose tolerant controls (Pune Maternal Nutrition Study, n = 461).

Findings

SIDD was the predominant subgroup followed by MOD, whereas SIRD and MARD were less frequent. Weighted-GRS for T2D, obesity and lipid-related traits associated with T2D. We replicated some of the previous associations of GRS for T2D, insulin secretion, and BMI with SIDD and MOD. Unique to Indian subgroups was the association of GRS for (a) proinsulin with MOD and MARD, (b) liver-lipids with SIDD, SIRD and MOD, and (c) opposite effect of beta-cell GRS with SIDD and MARD, obesity GRS with MARD compared to Europeans. Genetic variants of fucosyltransferases were associated with T2D and MOD in Indians but not Europeans.

Interpretation

The similarities emphasise the applicability of some of the European-derived GRSs to T2D and its subgroups in India while the differences highlight the need for large-scale studies to identify aetiologies in diverse ancestries. The data provide robust evidence for genetically distinct aetiologies for the T2D subgroups and at least partly mirror those seen in Europeans.

Funding

Vetenskapsrådet, Diabetes Wellness, and Hjärt-Lungfonden (Sweden), DST (India), Wellcome Trust, Crafoord Foundation and Albert Påhlsson Foundation.

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.

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.

MafB-dependent neurotransmitter signaling promotes β cell migration in the developing pancreas

Hormone secretion from pancreatic islets is essential for glucose homeostasis and loss or dysfunction of islet cells is a hallmark of type 2 diabetes. Maf transcription factors are critical for establishing and maintaining adult endocrine cell function. However, during pancreas development, MafB is not only expressed in insulin- and glucagon-producing cells, but also Neurog3+ endocrine progenitor cells suggesting additional functions in cell differentiation and islet formation. Here we report that MafB deficiency impairs β cell clustering and islet formation, but also coincides with loss of neurotransmitter and axon guidance receptor gene expression. Moreover, the observed loss of nicotinic receptor gene expression in human and mouse β cells implied that signaling through these receptors contributes to islet cell migration/formation. Inhibition of nicotinic receptor activity resulted in reduced β cell migration towards autonomic nerves and impaired β cell clustering. These findings highlight a novel function of MafB in controlling neuronal-directed signaling events required for islet formation.

Cross-Ancestry DNA Methylation Marks of Insulin Resistance in Pregnancy: An Integrative Epigenome-Wide Association Study

Although there are some epigenome-wide association studies (EWAS) of insulin resistance, for most of them authors did not replicate their findings, and most are focused on populations of European ancestry, limiting the generalizability. In the Epigenetics in Pregnancy (EPIPREG; n = 294 Europeans and 162 South Asians) study, we conducted an EWAS of insulin resistance in maternal peripheral blood leukocytes, with replication in the Born in Bradford (n = 879; n = 430 Europeans and 449 South Asians), Methyl Epigenome Network Association (MENA) (n = 320), and Botnia (n = 56) cohorts. In EPIPREG, we identified six CpG sites inversely associated with insulin resistance across ancestry, of which five were replicated in independent cohorts (cg02988288, cg19693031, and cg26974062 in TXNIP; cg06690548 in SLC7A11; and cg04861640 in ZSCAN26). From methylation quantitative trait loci analysis in EPIPREG, we identified gene variants related to all five replicated cross-ancestry CpG sites, which were associated with several cardiometabolic phenotypes. Mediation analyses suggested that the gene variants regulate insulin resistance through DNA methylation. To conclude, our cross-ancestry EWAS identified five CpG sites related to lower insulin resistance.

Type 2 diabetes candidate genes, including PAX5, cause impaired insulin secretion in human pancreatic islets

Type 2 diabetes (T2D) is caused by insufficient insulin secretion from pancreatic β cells. To identify candidate genes contributing to T2D pathophysiology, we studied human pancreatic islets from approximately 300 individuals. We found 395 differentially expressed genes (DEGs) in islets from individuals with T2D, including, to our knowledge, novel (OPRD1, PAX5, TET1) and previously identified (CHL1, GLRA1, IAPP) candidates. A third of the identified expression changes in islets may predispose to diabetes, as expression of these genes associated with HbA1c in individuals not previously diagnosed with T2D. Most DEGs were expressed in human β cells, based on single-cell RNA-Seq data. Additionally, DEGs displayed alterations in open chromatin and associated with T2D SNPs. Mouse KO strains demonstrated that the identified T2D-associated candidate genes regulate glucose homeostasis and body composition in vivo. Functional validation showed that mimicking T2D-associated changes for OPRD1, PAX5, and SLC2A2 impaired insulin secretion. Impairments in Pax5-overexpressing β cells were due to severe mitochondrial dysfunction. Finally, we discovered PAX5 as a potential transcriptional regulator of many T2D-associated DEGs in human islets. Overall, we have identified molecular alterations in human pancreatic islets that contribute to β cell dysfunction in T2D pathophysiology.

Epigenome-wide association study of serum folate in maternal peripheral blood leukocytes

Aim: To perform an epigenome-wide association study (EWAS) of serum folate in maternal blood. Methods: Cross-ancestry (Europeans = 302, South Asians = 161) and ancestry-specific EWAS in the EPIPREG cohort were performed, followed by methyl quantitative trait loci analysis and association with cardiometabolic phenotypes. Replication was attempted using maternal folate intake and blood methylation data from the MoBa study and verified if the findings were significant in a previous EWAS of maternal serum folate in cord blood. Results & conclusion: cg19888088 (cross-ancestry) in EBF3, cg01952260 (Europeans) and cg07077240 (South Asians) in HERC3 were associated with serum folate. cg19888088 and cg01952260 were associated with diastolic blood pressure. cg07077240 was associated with variants in CASC15. The findings were not replicated and were not significant in cord blood.

Multi-ancestry genome-wide association study of gestational diabetes mellitus highlights genetic links with type 2 diabetes

Gestational diabetes mellitus (GDM) is associated with increased risk of pregnancy complications and adverse perinatal outcomes. GDM often reoccurs and is associated with increased risk of subsequent diagnosis of type 2 diabetes (T2D). To improve our understanding of the aetiological factors and molecular processes driving the occurrence of GDM, including the extent to which these overlap with T2D pathophysiology, the GENetics of Diabetes In Pregnancy Consortium assembled genome-wide association studies of diverse ancestry in a total of 5485 women with GDM and 347 856 without GDM. Through multi-ancestry meta-analysis, we identified five loci with genome-wide significant association (P < 5 × 10-8) with GDM, mapping to/near MTNR1B (P = 4.3 × 10-54), TCF7L2 (P = 4.0 × 10-16), CDKAL1 (P = 1.6 × 10-14), CDKN2A-CDKN2B (P = 4.1 × 10-9) and HKDC1 (P = 2.9 × 10-8). Multiple lines of evidence pointed to the shared pathophysiology of GDM and T2D: (i) four of the five GDM loci (not HKDC1) have been previously reported at genome-wide significance for T2D; (ii) significant enrichment for associations with GDM at previously reported T2D loci; (iii) strong genetic correlation between GDM and T2D and (iv) enrichment of GDM associations mapping to genomic annotations in diabetes-relevant tissues and transcription factor binding sites. Mendelian randomization analyses demonstrated significant causal association (5% false discovery rate) of higher body mass index on increased GDM risk. Our results provide support for the hypothesis that GDM and T2D are part of the same underlying pathology but that, as exemplified by the HKDC1 locus, there are genetic determinants of GDM that are specific to glucose regulation in pregnancy.

Genetic regulation of RNA splicing in human pancreatic islets

BACKGROUND

Non-coding genetic variants that influence gene transcription in pancreatic islets play a major role in the susceptibility to type 2 diabetes (T2D), and likely also contribute to type 1 diabetes (T1D) risk. For many loci, however, the mechanisms through which non-coding variants influence diabetes susceptibility are unknown.

RESULTS

We examine splicing QTLs (sQTLs) in pancreatic islets from 399 human donors and observe that common genetic variation has a widespread influence on the splicing of genes with established roles in islet biology and diabetes. In parallel, we profile expression QTLs (eQTLs) and use transcriptome-wide association as well as genetic co-localization studies to assign islet sQTLs or eQTLs to T2D and T1D susceptibility signals, many of which lack candidate effector genes. This analysis reveals biologically plausible mechanisms, including the association of T2D with an sQTL that creates a nonsense isoform in ERO1B, a regulator of ER-stress and proinsulin biosynthesis. The expanded list of T2D risk effector genes reveals overrepresented pathways, including regulators of G-protein-mediated cAMP production. The analysis of sQTLs also reveals candidate effector genes for T1D susceptibility such as DCLRE1B, a senescence regulator, and lncRNA MEG3.

CONCLUSIONS

These data expose widespread effects of common genetic variants on RNA splicing in pancreatic islets. The results support a role for splicing variation in diabetes susceptibility, and offer a new set of genetic targets with potential therapeutic benefit.

Islet Gene View-a tool to facilitate islet research

Characterization of gene expression in pancreatic islets and its alteration in type 2 diabetes (T2D) are vital in understanding islet function and T2D pathogenesis. We leveraged RNA sequencing and genome-wide genotyping in islets from 188 donors to create the Islet Gene View (IGW) platform to make this information easily accessible to the scientific community. Expression data were related to islet phenotypes, diabetes status, other islet-expressed genes, islet hormone-encoding genes and for expression in insulin target tissues. The IGW web application produces output graphs for a particular gene of interest. In IGW, 284 differentially expressed genes (DEGs) were identified in T2D donor islets compared with controls. Forty percent of DEGs showed cell-type enrichment and a large proportion significantly co-expressed with islet hormone-encoding genes; glucagon (<i>GCG</i>, 56%), amylin (<i>IAPP</i>, 52%), insulin (<i>INS</i>, 44%), and somatostatin (<i>SST</i>, 24%). Inhibition of two DEGs, <i>UNC5D</i> and <i>SERPINE2</i>, impaired glucose-stimulated insulin secretion and impacted cell survival in a human β-cell model. The exploratory use of IGW could help designing more comprehensive functional follow-up studies and serve to identify therapeutic targets in T2D.

Neuronal Dysfunction Is Linked to the Famine-Associated Risk of Proliferative Retinopathy in Patients With Type 2 Diabetes

Persons with type 2 diabetes born in the regions of famine exposures have disproportionally elevated risk of vision-threatening proliferative diabetic retinopathy (PDR) in adulthood. However, the underlying mechanisms are not known. In the present study, we aimed to investigate the plausible molecular factors underlying progression to PDR. To study the association of genetic variants with PDR under the intrauterine famine exposure, we analyzed single nucleotide polymorphisms (SNPs) that were previously reported to be associated with type 2 diabetes, glucose, and pharmacogenetics. Analyses were performed in the population from northern Ukraine with a history of exposure to the Great Ukrainian Holodomor famine [the Diagnostic Optimization and Treatment of Diabetes and its Complications in the Chernihiv Region (DOLCE study), n = 3,583]. A validation of the top genetic findings was performed in the Hong Kong diabetes registry (HKDR, n = 730) with a history of famine as a consequence of the Japanese invasion during WWII. In DOLCE, the genetic risk for PDR was elevated for the variants in ADRA2A, PCSK9, and CYP2C19*2 loci, but reduced at PROX1 locus. The association of ADRA2A loci with the risk of advanced diabetic retinopathy in famine-exposed group was further replicated in HKDR. The exposure of embryonic retinal cells to starvation for glucose, mimicking the perinatal exposure to famine, resulted in sustained increased expression of Adra2a and Pcsk9, but decreased Prox1. The exposure to starvation exhibited a lasting inhibitory effects on neurite outgrowth, as determined by neurite length. In conclusion, a consistent genetic findings on the famine-linked risk of ADRA2A with PDR indicate that the nerves may likely to be responsible for communicating the effects of perinatal exposure to famine on the elevated risk of advanced stages of diabetic retinopathy in adults. These results suggest the possibility of utilizing neuroprotective drugs for the prevention and treatment of PDR.

Mapping the Cord Blood Transcriptome of Pregnancies Affected by Early Maternal Anemia to Identify Signatures of Fetal Programming

CONTEXT

Anemia during early pregnancy (EP) is common in developing countries and is associated with adverse health consequences for both mothers and children. Offspring of women with EP anemia often have low birth weight, which increases risk for cardiometabolic diseases, including type 2 diabetes (T2D), later in life.

OBJECTIVE

We aimed to elucidate mechanisms underlying developmental programming of adult cardiometabolic disease, including epigenetic and transcriptional alterations potentially detectable in umbilical cord blood (UCB) at time of birth.

METHODS

We leveraged global transcriptome- and accompanying epigenome-wide changes in 48 UCB from newborns of EP anemic Tanzanian mothers and 50 controls to identify differentially expressed genes (DEGs) in UCB exposed to maternal EP anemia. DEGs were assessed for association with neonatal anthropometry and cord insulin levels. These genes were further studied in expression data from human fetal pancreas and adult islets to understand their role in beta-cell development and/or function.

RESULTS

The expression of 137 genes was altered in UCB of newborns exposed to maternal EP anemia. These putative signatures of fetal programming, which included the birth weight locus LCORL, were potentially mediated by epigenetic changes in 27 genes and associated with neonatal anthropometry. Among the DEGs were P2RX7, PIK3C2B, and NUMBL, which potentially influence beta-cell development. Insulin levels were lower in EP anemia-exposed UCB, supporting the notion of developmental programming of pancreatic beta-cell dysfunction and subsequently increased risk of T2D in offspring of mothers with EP anemia.

CONCLUSIONS

Our data provide proof-of-concept on distinct transcriptional and epigenetic changes detectable in UCB from newborns exposed to maternal EP anemia.

Human pancreatic islet miRNA-mRNA networks of altered miRNAs due to glycemic status

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene expression via mRNA targeting, playing important roles in the pancreatic islets. We aimed to identify molecular pathways and genomic regulatory regions associated with altered miRNA expression due to glycemic status, which could contribute to the development of type 2 diabetes (T2D). To this end, miRNAs were identified by a combination of differential miRNA expression and correlation analysis in human islet samples from donors with normal and elevated blood glucose levels. Analysis and clustering of highly correlated, experimentally validated gene targets of these miRNAs revealed two islet-specific clusters, which were associated with key aspects of islet functions and included a high number of T2D-related genes. Finally, cis-eQTLs and public GWAS data integration uncovered suggestive genomic signals of association with insulin secretion and T2D. The miRNA-driven network-based approach presented in this study contributes to a better understanding of impaired insulin secretion in T2D pathogenesis.

•

Network-based islet miRNA approaches provide a better insight into T2D pathogenesis

•

Islets from healthy and hyperglycemic donors have distinct miRNA profiles

•

Differentially expressed miRNAs regulate gene networks that affect islet function

•

Specific miRNA eQTLs overlap with insulin secretion and T2D genetic markers

Perinatal famine is associated with excess risk of proliferative retinopathy in patients with type 2 diabetes

PURPOSE

Intrauterine undernutrition is associated with increased risk of type 2 diabetes. Children born premature or small for gestational age were reported to have abnormal retinal vascularization. However, whether intrauterine famine act as a trigger for diabetes complications, including retinopathy, is unknown. The aim of the current study was to evaluate long-term effects of perinatal famine on the risk of proliferative diabetic retinopathy (PDR).

METHODS

We studied the risk for PDR among type 2 diabetes patients exposed to perinatal famine in two independent cohorts: the Ukrainian National Diabetes Registry (UNDR) and the Hong Kong Diabetes Registry (HKDR). We analysed individuals born during the Great Famine (the Holodomor, 1932-1933) and the WWII (1941-1945) famine in 101 095 (3601 had PDR) UNDR participants. Among 3021 (251 had PDR) HKDR participants, we studied type 2 diabetes patients exposed to perinatal famine during the WWII Japanese invasion in 1942-1945.

RESULTS

During the Holodomor and WWII, perinatal famine was associated with a 1.76-fold (p = 0.019) and 3.02-fold (p = 0.001) increased risk of severe PDR in the UNDR. The risk for PDR was 1.66-fold elevated among individuals born in 1942 in the HKDR (p < 0.05). The associations between perinatal famine and PDR remained statistically significant after corrections for HbA1c in available 18 507 UNDR (p_{additive interaction}  < 0.001) and in 3021 HKDR type 2 diabetes patients (p < 0.05).

CONCLUSION

In conclusion, type 2 diabetes patients, exposed to perinatal famine, have increased risk of PDR compared to those without perinatal famine exposure. Further studies are needed to understand the underlying mechanisms and to extend this finding to other diabetes complications.

MAFA and MAFB regulate exocytosis-related genes in human β-cells

AIMS

Reduced expression of exocytotic genes is associated with functional defects in insulin exocytosis contributing to impaired insulin secretion and type 2 diabetes (T2D) development. MAFA and MAFB transcription factors regulate β-cell physiology, and their gene expression is reduced in T2D β cells. We investigate if loss of MAFA and MAFB in human β cells contributes to T2D progression by regulating genes required for insulin exocytosis.

METHODS

Three approaches were performed: (1) RNAseq analysis with the focus on exocytosis-related genes in MafA^-/- mouse islets, (2) correlational analysis between MAFA, MAFB and exocytosis-related genes in human islets and (3) MAFA and MAFB silencing in human islets and EndoC-βH1 cells followed by functional in vitro studies.

RESULTS

The expression of 30 exocytosis-related genes was significantly downregulated in MafA^-/- mouse islets. In human islets, the expression of 29 exocytosis-related genes correlated positively with MAFA and MAFB. Eight exocytosis-related genes were downregulated in MafA^-/- mouse islets and positively correlated with MAFA and MAFB in human islets. From this analysis, the expression of RAB3A, STXBP1, UNC13A, VAMP2, NAPA, NSF, STX1A and SYT7 was quantified after acute MAFA or MAFB silencing in EndoC-βH1 cells and human islets. MAFA and MAFB silencing resulted in impaired insulin secretion and reduced STX1A, SYT7 and STXBP1 (EndoC-βH1) and STX1A (human islets) mRNA expression. STX1A and STXBP1 protein expression was also impaired in islets from T2D donors which lack MAFA expression.

CONCLUSION

Our data indicate that STXBP1 and STX1A are important MAFA/B-regulated exocytosis genes which may contribute to insulin exocytosis defects observed in MAFA-deficient human T2D β cells.

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.

100 YEARS OF INSULIN: Towards improved precision and a new classification of diabetes mellitus

Type 2 diabetes (T2D) is one of the fastest increasing diseases worldwide. Although it is defined by a single metabolite, glucose, it is increasingly recognized as a highly heterogeneous disease with varying clinical manifestations. Identification of different subtypes at an early stage of disease when complications might still be prevented could hopefully allow for more personalized medicine. An important step toward precision medicine would be to target the right resources to the right patients, thereby improving patient health and reducing health costs for the society. More well-defined disease populations also offer increased power in experimental, genetic and clinical studies. In a recent study, we used six clinical variables (glutamate decarboxylase autoantibodies, age at onset of diabetes, glycated hemoglobin, BMI and simple measures of insulin resistance and insulin secretion (so called HOMA estimates) to cluster adult-onset diabetes patients into five subgroups. These subgroups have been robustly reproduced in several populations worldwide and are associated with different risks of diabetic complications and responses to treatment. Importantly, the group with severe insulin-deficient diabetes had increased risk of retinopathy and neuropathy, whereas the severe insulin-resistant diabetes group has the highest risk for diabetic kidney disease (DKD) and fatty liver. This emphasizes the key role of insulin resistance in the pathogenesis of DKD and fatty liver in T2D. In conclusion, this novel subclassification, breaking down T2D in clinically meaningful subgroups, provides the prerequisite framework for expanded personalized medicine in diabetes beyond what is already available for monogenic and to some extent type 1 diabetes.

Association of single nucleotide polymorphisms with insulin secretion, insulin sensitivity, and diabetes in women with a history of gestational diabetes mellitus

Background

This study investigated whether single nucleotide polymorphisms (SNPs) reported by previous genome-wide association studies (GWAS) to be associated with impaired insulin secretion, insulin resistance, and/or type 2 diabetes are associated with disposition index, the homeostasis model assessment of insulin resistance (HOMA-IR), and/or development of diabetes following a pregnancy complicated by gestational diabetes mellitus (GDM).

Methods

Seventy-two SNPs were genotyped in 374 women with previous GDM from Southern Sweden. An oral glucose tolerance test was performed 1–2 years postpartum, although data on the diagnosis of diabetes were accessible up to 5 years postpartum. HOMA-IR and disposition index were used to measure insulin resistance and secretion, respectively.

Results

The risk A-allele in the rs11708067 polymorphism of the adenylate cyclase 5 gene (ADCY5) was associated with decreased disposition index (beta = − 0.90, SE 0.38, p = 0.019). This polymorphism was an expression quantitative trait loci (eQTL) in islets for both ADCY5 and its antisense transcript. The risk C-allele in the rs2943641 polymorphism, near the insulin receptor substrate 1 gene (IRS1), showed a trend towards association with increased HOMA-IR (beta = 0.36, SE 0.18, p = 0.050), and the T-allele of the rs4607103 polymorphism, near the ADAM metallopeptidase with thrombospondin type 1 motif 9 gene (ADAMTS9), was associated with postpartum diabetes (OR = 2.12, SE 0.22, p = 0.00055). The genetic risk score (GRS) of the top four SNPs tested for association with the disposition index using equal weights was associated with the disposition index (beta = − 0.31, SE = 0.29, p = 0.00096). In addition, the GRS of the four SNPs studied for association with HOMA-IR using equal weights showed an association with HOMA-IR (beta = 1.13, SE = 0.48, p = 9.72874e−11). All analyses were adjusted for age, body mass index, and ethnicity.

Conclusions

This study demonstrated the genetic susceptibility of women with a history of GDM to impaired insulin secretion and sensitivity and, ultimately, to diabetes development.

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.

Cohort profile: Epigenetics in Pregnancy (EPIPREG) - population-based sample of European and South Asian pregnant women with epigenome-wide DNA methylation (850k) in peripheral blood leukocytes

Pregnancy is a valuable model to study the association between DNA methylation and several cardiometabolic traits, due to its direct potential to influence mother's and child's health. Epigenetics in Pregnancy (EPIPREG) is a population-based sample with the aim to study associations between DNA-methylation in pregnancy and cardiometabolic traits in South Asian and European pregnant women and their offspring. This cohort profile paper aims to present our sample with genetic and epigenetic data and invite researchers with similar cohorts to collaborative projects, such as replication of ours or their results and meta-analysis. In EPIPREG we have quantified epigenome-wide DNA methylation in maternal peripheral blood leukocytes in gestational week 28±1 in Europeans (n = 312) and South Asians (n = 168) that participated in the population-based cohort STORK Groruddalen, in Norway. DNA methylation was measured with Infinium MethylationEPIC BeadChip (850k sites), with technical validation of four CpG sites using bisulphite pyrosequencing in a subset (n = 30). The sample is well characterized with few missing data on e.g. genotype, universal screening for gestational diabetes, objectively measured physical activity, bioelectrical impedance, anthropometrics, biochemical measurements, and a biobank with maternal serum and plasma, urine, placenta tissue. In the offspring, we have repeated ultrasounds during pregnancy, cord blood, and anthropometrics up to 4 years of age. We have quantified DNA methylation in peripheral blood leukocytes in nearly all eligible women from the STORK Groruddalen study, to minimize the risk of selection bias. Genetic principal components distinctly separated Europeans and South Asian women, which fully corresponded with the self-reported ethnicity. Technical validation of 4 CpG sites from the methylation bead chip showed good agreement with bisulfite pyrosequencing. We plan to study associations between DNA methylation and cardiometabolic traits and outcomes.

The MafA-target gene PPP1R1A regulates GLP1R-mediated amplification of glucose-stimulated insulin secretion in β-cells

The amplification of glucose-stimulated insulin secretion (GSIS) through incretin signaling is critical for maintaining physiological glucose levels. Incretins, like glucagon-like peptide 1 (GLP1), are a target of type 2 diabetes drugs aiming to enhance insulin secretion. Here we show that the protein phosphatase 1 inhibitor protein 1A (PPP1R1A), is expressed in β-cells and that its expression is reduced in dysfunctional β-cells lacking MafA and upon acute MafA knock down. MafA is a central regulator of GSIS and β-cell function. We observed a strong correlation of MAFA and PPP1R1A mRNA levels in human islets, moreover, PPP1R1A mRNA levels were reduced in type 2 diabetic islets and positively correlated with GLP1-mediated GSIS amplification. PPP1R1A silencing in INS1 (832/13) β-cells impaired GSIS amplification, PKA-target protein phosphorylation, mitochondrial coupling efficiency and also the expression of critical β-cell marker genes like MafA, Pdx1, NeuroD1 and Pax6. Our results demonstrate that the β-cell transcription factor MafA is required for PPP1R1A expression and that reduced β-cell PPP1R1A levels impaired β-cell function and contributed to β-cell dedifferentiation during type 2 diabetes. Loss of PPP1R1A in type 2 diabetic β-cells may explains the unresponsiveness of type 2 diabetic patients to GLP1R-based treatments.

SCRT1 is a novel beta cell transcription factor with insulin regulatory properties

Here we show that scratch family transcriptional repressor 1 (SCRT1), a zinc finger transcriptional regulator, is a novel regulator of beta cell function. SCRT1 was found to be expressed in beta cells in rodent and human islets. In human islets, expression of SCRT1 correlated with insulin secretion capacity and the expression of the insulin (INS) gene. Furthermore, SCRT1 mRNA expression was lower in beta cells from T2D patients. siRNA-mediated Scrt1 silencing in INS-1832/13 cells, mouse- and human islets resulted in impaired glucose-stimulated insulin secretion and decreased expression of the insulin gene. This is most likely due to binding of SCRT1 to E-boxes of the Ins1 gene as shown with ChIP. Scrt1 silencing also reduced the expression of several key beta cell transcription factors. Moreover, Scrt1 mRNA expression was reduced by glucose and SCRT1 protein was found to translocate between the nucleus and the cytosol in a glucose-dependent fashion in INS-1832/13 cells as well as in a rodent model of T2D. SCRT1 was also regulated by a GSK3β-dependent SCRT1-serine phosphorylation. Taken together, SCRT1 is a novel beta cell transcription factor that regulates insulin secretion and is affected in T2D.

Robustness and lethality in multilayer biological molecular networks

Robustness is a prominent feature of most biological systems. Most previous related studies have been focused on homogeneous molecular networks. Here we propose a comprehensive framework for understanding how the interactions between genes, proteins and metabolites contribute to the determinants of robustness in a heterogeneous biological network. We integrate heterogeneous sources of data to construct a multilayer interaction network composed of a gene regulatory layer, a protein-protein interaction layer, and a metabolic layer. We design a simulated perturbation process to characterize the contribution of each gene to the overall system's robustness, and find that influential genes are enriched in essential and cancer genes. We show that the proposed mechanism predicts a higher vulnerability of the metabolic layer to perturbations applied to genes associated with metabolic diseases. Furthermore, we find that the real network is comparably or more robust than expected in multiple random realizations. Finally, we analytically derive the expected robustness of multilayer biological networks starting from the degree distributions within and between layers. These results provide insights into the non-trivial dynamics occurring in the cell after a genetic perturbation is applied, confirming the importance of including the coupling between different layers of interaction in models of complex biological systems.

Blood-based epigenetic estimators of chronological age in human adults using DNA methylation data from the Illumina MethylationEPIC array

BACKGROUND

Epigenetic clocks have been recognized for their precise prediction of chronological age, age-related diseases, and all-cause mortality. Existing epigenetic clocks are based on CpGs from the Illumina HumanMethylation450 BeadChip (450 K) which has now been replaced by the latest platform, Illumina MethylationEPIC BeadChip (EPIC). Thus, it remains unclear to what extent EPIC contributes to increased precision and accuracy in the prediction of chronological age.

RESULTS

We developed three blood-based epigenetic clocks for human adults using EPIC-based DNA methylation (DNAm) data from the Norwegian Mother, Father and Child Cohort Study (MoBa) and the Gene Expression Omnibus (GEO) public repository: 1) an Adult Blood-based EPIC Clock (ABEC) trained on DNAm data from MoBa (n = 1592, age-span: 19 to 59 years), 2) an extended ABEC (eABEC) trained on DNAm data from MoBa and GEO (n = 2227, age-span: 18 to 88 years), and 3) a common ABEC (cABEC) trained on the same training set as eABEC but restricted to CpGs common to 450 K and EPIC. Our clocks showed high precision (Pearson correlation between chronological and epigenetic age (r) > 0.94) in independent cohorts, including GSE111165 (n = 15), GSE115278 (n = 108), GSE132203 (n = 795), and the Epigenetics in Pregnancy (EPIPREG) study of the STORK Groruddalen Cohort (n = 470). This high precision is unlikely due to the use of EPIC, but rather due to the large sample size of the training set.

CONCLUSIONS

Our ABECs predicted adults' chronological age precisely in independent cohorts. As EPIC is now the dominant platform for measuring DNAm, these clocks will be useful in further predictions of chronological age, age-related diseases, and mortality.

Subtypes of Type 2 Diabetes Determined From Clinical Parameters

Type 2 diabetes (T2D) is defined by a single metabolite, glucose, but is increasingly recognized as a highly heterogeneous disease, including individuals with varying clinical characteristics, disease progression, drug response, and risk of complications. Identification of subtypes with differing risk profiles and disease etiologies at diagnosis could open up avenues for personalized medicine and allow clinical resources to be focused to the patients who would be most likely to develop diabetic complications, thereby both improving patient health and reducing costs for the health sector. More homogeneous populations also offer increased power in experimental, genetic, and clinical studies. Clinical parameters are easily available and reflect relevant disease pathways, including the effects of both genetic and environmental exposures. We used six clinical parameters (GAD autoantibodies, age at diabetes onset, HbA_1c, BMI, and measures of insulin resistance and insulin secretion) to cluster adult-onset diabetes patients into five subtypes. These subtypes have been robustly reproduced in several populations and associated with different risks of complications, comorbidities, genetics, and response to treatment. Importantly, the group with severe insulin-deficient diabetes (SIDD) had increased risk of retinopathy and neuropathy, whereas the severe insulin-resistant diabetes (SIRD) group had the highest risk for diabetic kidney disease (DKD) and fatty liver, emphasizing the importance of insulin resistance for DKD and hepatosteatosis in T2D. In conclusion, we believe that subclassification using these highly relevant parameters could provide a framework for personalized medicine in diabetes.

Genetic variant effects on gene expression in human pancreatic islets and their implications for T2D

Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, key tissues and cell-types required for functional inference are absent from large-scale resources. Here we explore the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using data from 420 donors. We find: (a) 7741 cis-eQTLs in islets with a replication rate across 44 GTEx tissues between 40% and 73%; (b) marked overlap between islet cis-eQTL signals and active regulatory sequences in islets, with reduced eQTL effect size observed in the stretch enhancers most strongly implicated in GWAS signal location; (c) enrichment of islet cis-eQTL signals with T2D risk variants identified in genome-wide association studies; and (d) colocalization between 47 islet cis-eQTLs and variants influencing T2D or glycemic traits, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in disease relevant tissues.

Ghrelin suppresses insulin secretion in human islets and type 2 diabetes patients have diminished islet ghrelin cell number and lower plasma ghrelin levels

It is not known how ghrelin affects insulin secretion in human islets from patients with type 2 diabetes (T2D) or whether islet ghrelin expression or circulating ghrelin levels are altered in T2D. Here we sought out to identify the effect of ghrelin on insulin secretion in human islets and the impact of T2D on circulating ghrelin levels and on islet ghrelin cells. The effect of ghrelin on insulin secretion was assessed in human T2D and non-T2D islets. Ghrelin expression was assessed with RNA-sequencing (n = 191) and immunohistochemistry (n = 21). Plasma ghrelin was measured with ELISA in 40 T2D and 40 non-T2D subjects. Ghrelin exerted a glucose-dependent insulin-suppressing effect in islets from both T2D and non-T2D donors. Compared with non-T2D donors, T2D donors had reduced ghrelin mRNA expression and 75% less islet ghrelin cells, and ghrelin mRNA expression correlated negatively with HbA1c. T2D subjects had 25% lower fasting plasma ghrelin levels than matched controls. Thus, ghrelin has direct insulin-suppressing effects in human islets and T2D patients have lower fasting ghrelin levels, likely as a result of reduced number of islet ghrelin cells. These findings support inhibition of ghrelin signaling as a potential therapeutic avenue for stimulation of insulin secretion in T2D patients.

Association between the rs1544410 polymorphism in the vitamin D receptor (VDR) gene and insulin secretion after gestational diabetes mellitus

Background and aims

Genetic variants involved in vitamin D metabolism have been associated with diabetes and related syndromes/diseases. We wanted to investigate possible associations of polymorphisms in genes involved in vitamin D metabolism with indices of insulin resistance and insulin secretion, and also with development of diabetes after gestational diabetes mellitus (GDM).

Materials and methods

We have studied 376 women with previous GDM. Eight single nucleotide polymorphisms (SNPs) in the genes for vitamin D receptor (VDR) [rs731236, rs7975232, rs10735810, and rs1544410], vitamin D binding protein (DBP) [rs7041 and rs4588], and cytochrome P450 family 27 subfamily B member 1 (CYP27B1) [rs10877012 and rs4646536] were genotyped by TaqMan Allelic Discrimination Assay using the Quantstudio 7 Flex system. A 75-g oral glucose tolerance test (OGTT) was performed 1–2 years postpartum. The homeostasis model assessment of insulin resistance (HOMA-IR) and the disposition index [(insulinogenic index: I30/G30)/HOMA-IR] were used to calculate insulin resistance and insulin secretion, respectively. Serum samples for determination of 25(OH)D3 were collected at the time of the OGTT. Manifestation of diabetes was followed up to five years postpartum.

Results

After adjustment for BMI, age, and ethnicity, the A-allele of the VDR rs1544410 polymorphism was found to be associated with increased disposition index (difference per allele = 3.56, 95% CI: 0.4567–6.674; p = 0.03). The A-allele of the DBP rs7041 polymorphism was found to be associated with 25(OH)D3 levels (difference [in nmol/L] per allele = −5.478, 95% CI: -8.315 to −2.641; p = 0.0002), as was the T-allele of the DBP rs4588 polymorphism (OR = −6.319, 95% CI: −9.466 to −3.171; p = 0.0001). None of the SNPs were significantly associated with HOMA-IR or postpartum diabetes.

Conclusions

This study provides evidence that the rs1544410 polymorphism of the VDR gene may be associated with increased insulin secretion in women after pregnancy complicated by GDM. Further studies in other populations are needed to confirm the results.

Glucose-dependent insulinotropic peptide and risk of cardiovascular events and mortality: a prospective study

AIMS/HYPOTHESIS

Evidence that glucose-dependent insulinotropic peptide (GIP) and/or the GIP receptor (GIPR) are involved in cardiovascular biology is emerging. We hypothesised that GIP has untoward effects on cardiovascular biology, in contrast to glucagon-like peptide 1 (GLP-1), and therefore investigated the effects of GIP and GLP-1 concentrations on cardiovascular disease (CVD) and mortality risk.

METHODS

GIP concentrations were successfully measured during OGTTs in two independent populations (Malmö Diet Cancer-Cardiovascular Cohort [MDC-CC] and Prevalence, Prediction and Prevention of Diabetes in Botnia [PPP-Botnia]) in a total of 8044 subjects. GLP-1 (n = 3625) was measured in MDC-CC. The incidence of CVD and mortality was assessed via national/regional registers or questionnaires. Further, a two-sample Mendelian randomisation (2SMR) analysis between the GIP pathway and outcomes (coronary artery disease [CAD] and myocardial infarction) was carried out using a GIP-associated genetic variant, rs1800437, as instrumental variable. An additional reverse 2SMR was performed with CAD as exposure variable and GIP as outcome variable, with the instrumental variables constructed from 114 known genetic risk variants for CAD.

RESULTS

In meta-analyses, higher fasting levels of GIP were associated with risk of higher total mortality (HR[95% CI] = 1.22 [1.11, 1.35]; p = 4.5 × 10^-5) and death from CVD (HR[95% CI] 1.30 [1.11, 1.52]; p = 0.001). In accordance, 2SMR analysis revealed that increasing GIP concentrations were associated with CAD and myocardial infarction, and an additional reverse 2SMR revealed no significant effect of CAD on GIP levels, thus confirming a possible effect solely of GIP on CAD.

CONCLUSIONS/INTERPRETATION

In two prospective, community-based studies, elevated levels of GIP were associated with greater risk of all-cause and cardiovascular mortality within 5-9 years of follow-up, whereas GLP-1 levels were not associated with excess risk. Further studies are warranted to determine the cardiovascular effects of GIP per se.

Life-long impairment of glucose homeostasis upon prenatal exposure to psychostimulants

Maternal drug abuse during pregnancy is a rapidly escalating societal problem. Psychostimulants, including amphetamine, cocaine, and methamphetamine, are amongst the illicit drugs most commonly consumed by pregnant women. Neuropharmacology concepts posit that psychostimulants affect monoamine signaling in the nervous system by their affinities to neurotransmitter reuptake and vesicular transporters to heighten neurotransmitter availability extracellularly. Exacerbated dopamine signaling is particularly considered as a key determinant of psychostimulant action. Much less is known about possible adverse effects of these drugs on peripheral organs, and if in utero exposure induces lifelong pathologies. Here, we addressed this question by combining human RNA-seq data with cellular and mouse models of neuroendocrine development. We show that episodic maternal exposure to psychostimulants during pregnancy coincident with the intrauterine specification of pancreatic β cells permanently impairs their ability of insulin production, leading to glucose intolerance in adult female but not male offspring. We link psychostimulant action specifically to serotonin signaling and implicate the sex-specific epigenetic reprogramming of serotonin-related gene regulatory networks upstream from the transcription factor Pet1/Fev as determinants of reduced insulin production.

Loss of ZnT8 function protects against diabetes by enhanced insulin secretion

A rare loss-of-function allele p.Arg138* in SLC30A8 encoding the zinc transporter 8 (ZnT8), which is enriched in Western Finland, protects against type 2 diabetes (T2D). We recruited relatives of the identified carriers and showed that protection was associated with better insulin secretion due to enhanced glucose responsiveness and proinsulin conversion, particularly when compared with individuals matched for the genotype of a common T2D-risk allele in SLC30A8, p.Arg325. In genome-edited human induced pluripotent stem cell (iPSC)-derived β-like cells, we establish that the p.Arg138* allele results in reduced SLC30A8 expression due to haploinsufficiency. In human β cells, loss of SLC30A8 leads to increased glucose responsiveness and reduced KATP channel function similar to isolated islets from carriers of the T2D-protective allele p.Trp325. These data position ZnT8 as an appealing target for treatment aimed at maintaining insulin secretion capacity in T2D.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

Genome editing of human pancreatic beta cell models: problems, possibilities and outlook

Understanding the molecular mechanisms behind beta cell dysfunction is essential for the development of effective and specific approaches for diabetes care and prevention. Physiological human beta cell models are needed for this work. We review the possibilities and limitations of currently available human beta cell models and how they can be dramatically enhanced using genome-editing technologies. In addition to the gold standard, primary isolated islets, other models now include immortalised human beta cell lines and pluripotent stem cell-derived islet-like cells. The scarcity of human primary islet samples limits their use, but valuable gene expression and functional data from large collections of human islets have been made available to the scientific community. The possibilities for studying beta cell physiology using immortalised human beta cell lines and stem cell-derived islets are rapidly evolving. However, the functional immaturity of these cells is still a significant limitation. CRISPR-Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR-associated protein 9) has enabled precise engineering of specific genetic variants, targeted transcriptional modulation and genome-wide genetic screening. These approaches can now be exploited to gain understanding of the mechanisms behind coding and non-coding diabetes-associated genetic variants, allowing more precise evaluation of their contribution to diabetes pathogenesis. Despite all the progress, genome editing in primary pancreatic islets remains difficult to achieve, an important limitation requiring further technological development.

Phenotypic and genotypic differences between Indian and Scandinavian women with gestational diabetes mellitus

OBJECTIVE

Gestational diabetes mellitus (GDM) is a transient form of diabetes characterized by impaired insulin secretion and action during pregnancy. Population-based differences in prevalence exist which could be explained by phenotypic and genetic differences. The aim of this study was to examine these differences in pregnant women from Punjab, India and Scandinavia.

METHODS

Eighty-five GDM/T2D loci in European and/or Indian populations from previous studies were assessed for association with GDM based on Swedish GDM criteria in 4018 Punjabi Indian and 507 Swedish pregnant women. Selected loci were replicated in Scandinavian cohorts, Radiel (N = 398, Finnish) and STORK/STORK-G (N = 780, Norwegian).

RESULTS

Punjabi Indian women had higher GDM prevalence, lower insulin secretion and better insulin sensitivity than Swedish women. There were significant frequency differences of GDM/T2D risk alleles between both populations. rs7178572 at HMG20A, previously associated with GDM in South Indian and European women, was replicated in North Indian women. The T2D risk SNP rs11605924 in the CRY2 gene was associated with increased GDM risk in Scandinavian but decreased GDM risk in Punjabi Indian women. No other overlap was seen between GDM loci in both populations.

CONCLUSIONS

Gestational diabetes mellitus is more common in Indian than Swedish women, which partially can be attributed to differences in insulin secretion and action. There was marked heterogeneity in the GDM phenotypes between the populations which could only partially be explained by genetic differences.

Loss of MafA and MafB expression promotes islet inflammation

Maf transcription factors are critical regulators of beta-cell function. We have previously shown that reduced MafA expression in human and mouse islets is associated with a pro-inflammatory gene signature. Here, we investigate if the loss of Maf transcription factors induced autoimmune processes in the pancreas. Transcriptomics analysis showed expression of pro-inflammatory as well as immune cell marker genes. However, clusters of CD4+ T and B220+ B cells were associated primarily with adult MafA^−/−MafB^+/−, but not MafA^−/− islets. MafA expression was detected in the thymus, lymph nodes and bone marrow suggesting a novel role of MafA in regulating immune-cell function. Analysis of pancreatic lymph node cells showed activation of CD4+ T cells, but lack of CD8+ T cell activation which also coincided with an enrichment of naïve CD8+ T cells. Further analysis of T cell marker genes revealed a reduction of T cell receptor signaling gene expression in CD8, but not in CD4+ T cells, which was accompanied with a defect in early T cell receptor signaling in mutant CD8+ T cells. These results suggest that loss of MafA impairs both beta- and T cell function affecting the balance of peripheral immune responses against islet autoantigens, resulting in local inflammation in pancreatic islets.

RORB and RORC associate with human islet dysfunction and inhibit insulin secretion in INS-1 cells

Little is known about the expression and function of Retinoic acid-related orphan receptors (RORA, B, and C) in pancreatic β cells. Here in, we utilized cDNA microarray and RNA sequencing approaches to investigate the expression pattern of ROR receptors in normal and diabetic human pancreatic islets. Possible correlations between RORs expression and HbA_1c levels as well as insulin secretory capacity in isolated human islets were evaluated. The impact of RORB and RORC expression on insulin secretion in INS-1 (832/13) cells was validated as well. While RORA was the highest expressed gene among the three RORs in human islet cells, RORC was the highest expressed in INS-1 cells (832/13) and while RORB was the lowest expressed gene in human islet cells, RORA was the highest expressed in INS-1 cells (832/13). The expression of RORB and RORC was significantly lower in diabetic/hyperglycemic donors as compared with non-diabetic counterparts. Furthermore, while the expression of RORB correlated positively with insulin secretion and negatively with HbA_1c, that of RORC correlated negatively with HbA_1c. The expression pattern of RORA did not correlate with either of the two parameters. siRNA silencing of RORB or RORC in INS-1 (832/13) cells resulted in a significant downregulation of insulin mRNA expression and insulin secretion. These findings suggest that RORB and RORC are part of the molecular cascade that regulates insulin secretion in pancreatic β cells; and insight that provides for further work on the potential therapeutic utility of RORB and RORC genes in β cell dysfunction in type 2 diabetes.

Precision medicine in type 2 diabetes

The Precision Medicine Initiative defines precision medicine as 'an emerging approach for disease treatment and prevention that takes into account individual variability in genes, environment and lifestyle for each person'. This approach will facilitate more accurate treatment and prevention strategies in contrast to a one-size-fits-all approach, in which disease treatment and prevention strategies are developed for generalized usage. Diabetes is clearly more heterogeneous than the conventional subclassification into type 1 and type 2 diabetes. Monogenic forms of diabetes like MODY and neonatal diabetes have paved the way for precision medicine in diabetes, as carriers of unique mutations require unique treatment. Diagnosis of diabetes in the past has been dependent upon measuring one metabolite, glucose. By instead including six variables in a clustering analysis, we could break down diabetes into five distinct subgroups, with better prediction of disease progression and outcome. The severe insulin-resistant diabetes (SIRD) cluster showed the highest risk of kidney disease and highest prevalence of nonalcoholic fatty liver disease, whereas patients in the insulin-deficient cluster 2 (SIDD) had the highest risk of retinopathy. In the future, this will certainly be improved and expanded by including genetic, epigenetic and other biomarker to allow better prediction of outcome and choice of more precise treatment.

MafA Expression Preserves Immune Homeostasis in Human and Mouse Islets

Type 1 (T1D) and type 2 (T2D) diabetes are triggered by a combination of environmental and/or genetic factors. Maf transcription factors regulate pancreatic beta (β)-cell function, and have also been implicated in the regulation of immunomodulatory cytokines like interferon-β (IFNβ1). In this study, we assessed MAFA and MAFB co-expression with pro-inflammatory cytokine signaling genes in RNA-seq data from human pancreatic islets. Interestingly, MAFA expression was strongly negatively correlated with cytokine-induced signaling (such as IFNAR1, DDX58) and T1D susceptibility genes (IFIH1), whereas correlation of these genes with MAFB was weaker. In order to evaluate if the loss of MafA altered the immune status of islets, MafA deficient mouse islets (MafA^−/−) were assessed for inherent anti-viral response and susceptibility to enterovirus infection. MafA deficient mouse islets had elevated basal levels of Ifnβ1, Rig1 (DDX58 in humans), and Mda5 (IFIH1) which resulted in reduced virus propagation in response to coxsackievirus B3 (CVB3) infection. Moreover, an acute knockdown of MafA in β-cell lines also enhanced Rig1 and Mda5 protein levels. Our results suggest that precise regulation of MAFA levels is critical for islet cell-specific cytokine production, which is a critical parameter for the inflammatory status of pancreatic islets.

Genetic determinants of glucose levels in pregnancy: genetic risk scores analysis and GWAS in the Norwegian STORK cohort

Objective Hyperglycaemia during pregnancy increases the risk of adverse health outcomes in mother and child, but the genetic aetiology is scarcely studied. Our aims were to (1) assess the overlapping genetic aetiology between the pregnant and non-pregnant population and (2) assess the importance of genome-wide polygenic contributions to glucose traits during pregnancy, by exploring whether genetic risk scores (GRSs) for fasting glucose (FG), 2-h glucose (2hG), type 2 diabetes (T2D) and BMI in non-pregnant individuals were associated with glucose measures in pregnant women. Methods We genotyped 529 Norwegian pregnant women and constructed GRS from known genome-wide significant variants and SNPs weakly associated (p > 5 × 10-8) with FG, 2hG, BMI and T2D from external genome-wide association studies (GWAS) and examined the association between these scores and glucose measures at gestational weeks 14-16 and 30-32. We also performed GWAS of FG, 2hG and shape information from the glucose curve during an oral glucose tolerance test (OGTT). Results GRSFG explained similar variance during pregnancy as in the non-pregnant population (~5%). GRSBMI and GRST2D explained up to 1.3% of the variation in the glucose traits in pregnancy. If we included variants more weakly associated with these traits, GRS2hG and GRST2D explained up to 2.4% of the variation in the glucose traits in pregnancy, highlighting the importance of polygenic contributions. Conclusions Our results suggest overlap in the genetic aetiology of FG in pregnant and non-pregnant individuals. This was less apparent with 2hG, suggesting potential differences in postprandial glucose metabolism inside and outside of pregnancy.

Turning Vice into Virtue: Using Batch-Effects to Detect Errors in Large Genomic Data Sets

It is often unavoidable to combine data from different sequencing centers or sequencing platforms when compiling data sets with a large number of individuals. However, the different data are likely to contain specific systematic errors that will appear as SNPs. Here, we devise a method to detect systematic errors in combined data sets. To measure quality differences between individual genomes, we study pairs of variants that reside on different chromosomes and co-occur in individuals. The abundance of these pairs of variants in different genomes is then used to detect systematic errors due to batch effects. Applying our method to the 1000 Genomes data set, we find that coding regions are enriched for errors, where ∼1% of the higher frequency variants are predicted to be erroneous, whereas errors outside of coding regions are much rarer (<0.001%). As expected, predicted errors are found less often than other variants in a data set that was generated with a different sequencing technology, indicating that many of the candidates are indeed errors. However, predicted 1000 Genomes errors are also found in other large data sets; our observation is thus not specific to the 1000 Genomes data set. Our results show that batch effects can be turned into a virtue by using the resulting variation in large scale data sets to detect systematic errors.

Silencing of the FTO gene inhibits insulin secretion: An in vitro study using GRINCH cells

Expression of fat mass and obesity-associated gene (FTO) and ADP-ribosylation factor-like 15 (ARL15) in human islets is inversely correlated with HbA_1c. However, their impact on insulin secretion is still ambiguous. Here in, we investigated the role of FTO and ARL15 using GRINCH (Glucose-Responsive Insulin-secreting C-peptide-modified Human proinsulin) clonal rat β-cells. GRINCH cells have inserted GFP into the human C-peptide insulin gene. Hence, secreted CpepGFP served to monitor insulin secretion. mRNA silencing of FTO in GRINCH cells showed a significant reduction in glucose but not depolarization-stimulated insulin secretion, whereas ARL15 silencing had no effect. A significant down-regulation of insulin mRNA was observed in FTO knockdown cells. Type-2 Diabetic islets revealed a reduced expression of FTO mRNA. In conclusion, our data suggest that fluorescent CpepGFP released from GRINCH cells may serve as a convenient marker for insulin secretion. Silencing of FTO expression, but not ARL15, inhibits insulin secretion by affecting metabolic signaling.

Association between genetic risk variants and glucose intolerance during pregnancy in north Indian women

Background

Gestational diabetes (GDM) is a more common problem in India than in many other parts of the world but it is not known whether this is due to unique environmental factors or a unique genetic background. To address this question we examined whether the same genetic variants associated with GDM and Type 2 Diabetes (T2D) in Caucasians also were associated with GDM in North Indian women.

Methods

Five thousand one hundred pregnant women of gestational age 24–28 weeks from Punjab were studied by a 75 g oral glucose tolerance test (OGTT). GDM was diagnosed by both WHO1999 and 2013 criteria. 79 single nucleotide polymorphisms (SNPs) previously associated with T2D and glycemic traits (12 of them also with GDM) and 6 SNPs from previous T2D associations based on Indian population (some also with European) were genotyped on a Sequenom platform or using Taqman assays in DNA from 4018 women.

Results

In support of previous findings in Caucasian GDM, SNPs at KCJN11 and GRB14 loci were nominally associated with GDM1999 risk in Indian women (both p = 0.02). Notably, T2D risk alleles of the variant rs1552224 near CENTD2, rs11708067 in ADCY5 and rs11605924 in CRY2 genes associated with protection from GDM regardless of criteria applied (p < 0.025). SNPs rs7607980 near COBLL1 (p = 0.0001), rs13389219 near GRB14 (p = 0.026) and rs10423928 in the GIPR gene (p = 0.012) as well as the genetic risk score (GRS) for these previously shown insulin resistance loci here associated with insulin resistance defined by HOMA2-IR and showed a trend towards GDM. GRS comprised of 3 insulin secretion loci here associated with insulin secretion but not GDM.

Conclusions

GDM in women from Punjab in Northern India shows a genetic component, seemingly driven by insulin resistance and secretion and partly shared with GDM in other parts of the world. Most previous T2D loci discovered in European studies did not associate with GDM in North India, indicative of different genetic etiology or alternately, differences in the linkage disequilibrium (LD) structure between populations in which the associated SNPs were identified and Northern Indian women. Interestingly some T2D risk variants were in fact indicative of being protective for GDM in these Indian women.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0380-8) contains supplementary material, which is available to authorized users.

Controllability in an islet specific regulatory network identifies the transcriptional factor NFATC4, which regulates Type 2 Diabetes associated genes

Probing the dynamic control features of biological networks represents a new frontier in capturing the dysregulated pathways in complex diseases. Here, using patient samples obtained from a pancreatic islet transplantation program, we constructed a tissue-specific gene regulatory network and used the control centrality (Cc) concept to identify the high control centrality (HiCc) pathways, which might serve as key pathobiological pathways for Type 2 Diabetes (T2D). We found that HiCc pathway genes were significantly enriched with modest GWAS p-values in the DIAbetes Genetics Replication And Meta-analysis (DIAGRAM) study. We identified variants regulating gene expression (expression quantitative loci, eQTL) of HiCc pathway genes in islet samples. These eQTL genes showed higher levels of differential expression compared to non-eQTL genes in low, medium, and high glucose concentrations in rat islets. Among genes with highly significant eQTL evidence, NFATC4 belonged to four HiCc pathways. We asked if the expressions of T2D-associated candidate genes from GWAS and literature are regulated by Nfatc4 in rat islets. Extensive in vitro silencing of Nfatc4 in rat islet cells displayed reduced expression of 16, and increased expression of four putative downstream T2D genes. Overall, our approach uncovers the mechanistic connection of NFATC4 with downstream targets including a previously unknown one, TCF7L2, and establishes the HiCc pathways' relationship to T2D.

Genetic inactivation of ANGPTL4 improves glucose homeostasis and is associated with reduced risk of diabetes

Angiopoietin-like 4 (ANGPTL4) is an endogenous inhibitor of lipoprotein lipase that modulates lipid levels, coronary atherosclerosis risk, and nutrient partitioning. We hypothesize that loss of ANGPTL4 function might improve glucose homeostasis and decrease risk of type 2 diabetes (T2D). We investigate protein-altering variants in ANGPTL4 among 58,124 participants in the DiscovEHR human genetics study, with follow-up studies in 82,766 T2D cases and 498,761 controls. Carriers of p.E40K, a variant that abolishes ANGPTL4 ability to inhibit lipoprotein lipase, have lower odds of T2D (odds ratio 0.89, 95% confidence interval 0.85–0.92, p = 6.3 × 10⁻¹⁰), lower fasting glucose, and greater insulin sensitivity. Predicted loss-of-function variants are associated with lower odds of T2D among 32,015 cases and 84,006 controls (odds ratio 0.71, 95% confidence interval 0.49–0.99, p = 0.041). Functional studies in Angptl4-deficient mice confirm improved insulin sensitivity and glucose homeostasis. In conclusion, genetic inactivation of ANGPTL4 is associated with improved glucose homeostasis and reduced risk of T2D.

Genetic variation in ANGPTL4 is associated with lipid traits. Here, the authors find that predicted loss-of-function variants in ANGPTL4 are associated with glucose homeostasis and reduced risk of type 2 diabetes and that Angptl4^−/− mice on a high-fat diet show improved insulin sensitivity.