E2F transcription factors promote tumorigenicity in pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with limited treatment options, illustrating an urgent need to identify new drugable targets in PDACs. OBJECTIVE: Using the similarities between tumor development and normal embryonic development, which is accompanied by rapid cell expansion, we aimed to identify and characterize embryonic signaling pathways that were reinitiated during tumor formation and expansion. METHODS AND RESULTS: Here, we report that the transcription factors E2F1 and E2F8 are potential key regulators in PDAC. E2F1 and E2F8 RNA expression is mainly localized in proliferating cells in the developing pancreas and in malignant ductal cells in PDAC. Silencing of E2F1 and E2F8 in PANC-1 pancreatic tumor cells inhibited cell proliferation and impaired cell spreading and migration. Moreover, loss of E2F1 also affected cell viability and apoptosis with E2F expression in PDAC tissues correlating with expression of apoptosis and mitosis pathway genes, suggesting that E2F factors promote cell cycle regulation and tumorigenesis in PDAC cells. CONCLUSION: Our findings illustrate that E2F1 and E2F8 transcription factors are expressed in pancreatic progenitor and PDAC cells, where they contribute to tumor cell expansion by regulation of cell proliferation, viability, and cell migration making these genes attractive therapeutic targets and potential prognostic markers for pancreatic cancer.

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.

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.

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 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.

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.

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 crucial 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 in 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.

Cryo-EM structure supports a role of AQP7 as a junction protein.

Aquaglyceroporin 7 (AQP7) facilitates glycerol flux across the plasma membrane with a critical physiological role linked to metabolism, obesity, and associated diseases. Here, we present the single-particle cryo-EM structure of AQP7 determined at 2.55 Å resolution adopting two adhering tetramers, stabilized by extracellularly exposed loops, in a configuration like that of the well-characterized interaction of AQP0 tetramers. The central pore, in-between the four monomers, displays well-defined densities restricted by two leucine filters. Gas chromatography mass spectrometry (GC/MS) results show that the AQP7 sample contains glycerol 3-phosphate (Gro3P), which is compatible with the identified features in the central pore. AQP7 is shown to be highly expressed in human pancreatic α- and ß- cells suggesting that the identified AQP7 octamer assembly, in addition to its function as glycerol channel, may serve as junction proteins within the endocrine pancreas.

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.

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.

The human batokine EPDR1 regulates ß-cell metabolism and function.

OBJECTIVE: Ependymin-Related Protein 1 (EPDR1) was recently identified as a secreted human batokine regulating mitochondrial respiration linked to thermogenesis in brown fat. Despite that EPDR1 is expressed in human pancreatic ß-cells and that glucose-stimulated mitochondrial metabolism is critical for stimulus-secretion coupling in ß-cells, the role of EPDR1 in ß-cell metabolism and function has not been investigated. METHODS: EPDR1 mRNA levels in human pancreatic islets from non-diabetic (ND) and type 2 diabetes (T2D) subjects were assessed. Human islets, EndoC-ßH1 and INS1 832/13 cells were transfected with scramble (control) and EPDR1 siRNAs (EPDR1-KD) or treated with human EPDR1 protein, and glucose-stimulated insulin secretion (GSIS) assessed by ELISA. Mitochondrial metabolism was investigated by extracellular flux analyzer, confocal microscopy and mass spectrometry-based metabolomics analysis. RESULTS: EPDR1 mRNA expression was upregulated in human islets from T2D and obese donors and positively correlated to BMI of donors. In T2D donors, EPDR1 mRNA levels negatively correlated with HbA1c and positively correlated with GSIS. EPDR1 silencing in human islets and ß-cell lines reduced GSIS whereas treatment with human EPDR1 protein increased GSIS. Epdr1 silencing in INS1 832/13 cells reduced glucose- and pyruvate- but not K+-stimulated insulin secretion. Metabolomics analysis in Epdr1-KD INS1 832/13 cells suggests diversion of glucose-derived pyruvate to lactate production and decreased malate-aspartate shuttle and the tricarboxylic acid (TCA) cycle activity. The glucose-stimulated rise in mitochondrial respiration and ATP/ADP-ratio was impaired in Epdr1-deficient cells. CONCLUSION: These results suggests that to maintain glucose homeostasis in obese people, upregulation of EPDR1 may improve ß-cell function via channelling glycolysis-derived pyruvate to the mitochondrial TCA cycle.

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.

Desmoglein-2 is important for islet function and ß-cell survival.

Type 1 diabetes is a complex disease characterized by the lack of endogenous insulin secreted from the pancreatic ß-cells. Although ß-cell targeted autoimmune processes and ß-cell dysfunction are known to occur in type 1 diabetes, a complete understanding of the cell-to-cell interactions that support pancreatic function is still lacking. To characterize the pancreatic endocrine compartment, we studied pancreata from healthy adult donors and investigated a single cell surface adhesion molecule, desmoglein-2 (DSG2). Genetically-modified mice lacking Dsg2 were examined for islet cell mass, insulin production, responses to glucose, susceptibility to a streptozotocin-induced mouse model of hyperglycaemia, and ability to cure diabetes in a syngeneic transplantation model. Herein, we have identified DSG2 as a previously unrecognized adhesion molecule that supports ß-cells. Furthermore, we reveal that DSG2 is within the top 10 percent of all genes expressed by human pancreatic islets and is expressed by the insulin-producing ß-cells but not the somatostatin-producing δ-cells. In a Dsg2 loss-of-function mice (Dsg2lo/lo), we observed a significant reduction in the number of pancreatic islets and islet size, and consequently, there was less total insulin content per islet cluster. Dsg2lo/lo mice also exhibited a reduction in blood vessel barrier integrity, an increased incidence of streptozotocin-induced diabetes, and islets isolated from Dsg2lo/lo mice were more susceptible to cytokine-induced ß-cell apoptosis. Following transplantation into diabetic mice, islets isolated from Dsg2lo/lo mice were less effective than their wildtype counterparts at curing diabetes. In vitro assays using the Beta-TC-6 murine ß-cell line suggest that DSG2 supports the actin cytoskeleton as well as the release of cytokines and chemokines. Taken together, our study suggests that DSG2 is an under-appreciated regulator of ß-cell function in pancreatic islets and that a better understanding of this adhesion molecule may provide new opportunities to combat type 1 diabetes.

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.

Twins in Guinea-Bissau have a 'thin-fat' body composition compared to singletons.

The 'thrifty phenotype' hypothesis proposed that fetal undernutrition increases risk of diabetes in later life. Undernourished low birthweight Indian babies are paradoxically more adipose compared to well-nourished European babies, and are at higher risk of diabetes in later life. Twin pregnancies are an example of in utero growth restrictive environment due to shared maternal nutrition. There are few studies of body composition in twins. We performed secondary analysis of anthropometric body composition of twins and singletons in Guinea-Bissau, an economically deprived African country.Anthropometric data were available on 7-34 year-old twins (n = 209, 97 males) and singletons (n = 182, 86 males) in the Guinea-Bissau Twin Registry at the Bandim Health Project. Twins had lower birthweight (2420 vs 3100 g, p < 0.001); and at follow-up, lower height (HAZ mean Z-score difference, -0.21, p = 0.055), weight (WAZ -0.73, p = 0.024) and BMI (BAZ -0.22, p = 0.079) compared to singletons but higher adiposity (skinfolds: +0.33 SD, p = 0.001). Twins also had higher fasting (+0.38 SD, p < 0.001) and 2-hour OGTT glucose concentrations (+0.29 SD, p < 0.05). Linear mixed-effect model accounting for intrapair correlations and interactions confirmed that twins were thinner but fatter across the age range. Data on maternal morbidity and prematurity were not available in this cohort.African populations are known to have a muscular (less adipose) body composition. Demonstration of a thin-fat phenotype in twins in a low socio-economic African country supports the thesis that it could be a manifestation of early life undernutrition and not exclusive to Indians. This phenotype could increase risk of diabetes and related conditions.

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.

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.