A systems genetics approach identifies genes and pathways for type 2 diabetes in human islets

Characterization of eQTLs associated with androstenone by RNA sequencing in porcine testis

Characterization of genetic variants affecting genome-wide gene expression levels (expression quantitative trait loci or eQTLs) in pig testes may improve our understanding of genetic architecture of boar taint (an animal welfare trait) and helps in genome-assisted or genomic selection programs. The aims of this study were to identify eQTLs associated with androstenone, to find candidate eQTLs for low androstenone, and to validate the top eQTL by reverse transcriptase quantitative PCR (RT-qPCR). Gene expression profiles were obtained by RNA sequencing in testis from Danish cross-bred pigs and genotype data by 80K single nucleotide polymorphism panel. A total of 262 eQTLs [false discovery rate (FDR) < 0.05] were identified by using two software packages: Matrix eQTL and Krux eQTL. Of these, 149 cis-acting eQTLs were significantly associated with androstenone concentrations and gene expression (FDR < 0.05). The eQTLs were associated with several genes of boar taint relevance including CYP1A2, CYB5D1, and SPHK2. One eQTL gene, AMPH, was differentially expressed (FDR < 0.05) and affected by chicory. Five candidate eQTLs associated with low androstenone concentrations were discovered, including the top eQTL associated with CYP1A2. RT-qPCR confirmed target gene expression to be significantly (P < 0.05) different based on eQTL genotypes. Furthermore, eQTLs were enriched as QTLs for 15 boar taint related traits from the PigQTLdb. This is the first study to report eQTLs in testes of commercial crossbred pigs used in pork production and to reveal genetic architecture of boar taint. Potential applications include development of a DNA test and in advanced genomic selection models for boar taint.

Contribution of Oxidative Stress and Impaired Biogenesis of Pancreatic β-Cells to Type 2 Diabetes

Significance: Type 2 diabetes development involves multiple changes in β-cells, related to the oxidative stress and impaired redox signaling, beginning frequently by sustained overfeeding due to the resulting lipotoxicity and glucotoxicity. Uncovering relationships among the dysregulated metabolism, impaired β-cell "well-being," biogenesis, or cross talk with peripheral insulin resistance is required for elucidation of type 2 diabetes etiology. Recent Advances: It has been recognized that the oxidative stress, lipotoxicity, and glucotoxicity cannot be separated from numerous other cell pathology events, such as the attempted compensation of β-cell for the increased insulin demand and dynamics of β-cell biogenesis and its "reversal" at dedifferentiation, that is, from the concomitantly decreasing islet β-cell mass (also due to transdifferentiation) and low-grade islet or systemic inflammation. Critical Issues: At prediabetes, the compensation responses of β-cells, attempting to delay the pathology progression-when exaggerated-set a new state, in which a self-checking redox signaling related to the expression of Ins gene expression is impaired. The resulting altered redox signaling, diminished insulin secretion responses to various secretagogues including glucose, may lead to excretion of cytokines or chemokines by β-cells or excretion of endosomes. They could substantiate putative stress signals to the periphery. Subsequent changes and lasting glucolipotoxicity promote islet inflammatory responses and further pathology spiral. Future Directions: Should bring an understanding of the β-cell self-checking and related redox signaling, including the putative stress signal to periphery. Strategies to cure or prevent type 2 diabetes could be based on the substitution of the "wrong" signal by the "correct" self-checking signal.

Free Fatty Acid Receptors in Health and Disease

Fatty acids are metabolized and synthesized as energy substrates during biological responses. Long- and medium-chain fatty acids derived mainly from dietary triglycerides, and short-chain fatty acids (SCFAs) produced by gut microbial fermentation of the otherwise indigestible dietary fiber, constitute the major sources of free fatty acids (FFAs) in the metabolic network. Recently, increasing evidence indicates that FFAs serve not only as energy sources but also as natural ligands for a group of orphan G protein-coupled receptors (GPCRs) termed free fatty acid receptors (FFARs), essentially intertwining metabolism and immunity in multiple ways, such as via inflammation regulation and secretion of peptide hormones. To date, several FFARs that are activated by the FFAs of various chain lengths have been identified and characterized. In particular, FFAR1 (GPR40) and FFAR4 (GPR120) are activated by long-chain saturated and unsaturated fatty acids, while FFAR3 (GPR41) and FFAR2 (GPR43) are activated by SCFAs, mainly acetate, butyrate, and propionate. In this review, we discuss the recent reports on the key physiological functions of the FFAR-mediated signaling transduction pathways in the regulation of metabolism and immune responses. We also attempt to reveal future research opportunities for developing therapeutics for metabolic and immune disorders.

Targeting BCAA Catabolism to Treat Obesity-Associated Insulin Resistance

Recent studies implicate a strong association between elevated plasma branched-chain amino acids (BCAAs) and insulin resistance (IR). However, a causal relationship and whether interrupted BCAA homeostasis can serve as a therapeutic target for diabetes remain to be established experimentally. In this study, unbiased integrative pathway analyses identified a unique genetic link between obesity-associated IR and BCAA catabolic gene expression at the pathway level in human and mouse populations. In genetically obese (ob/ob) mice, rate-limiting branched-chain α-keto acid (BCKA) dehydrogenase deficiency (i.e., BCAA and BCKA accumulation), a metabolic feature, accompanied the systemic suppression of BCAA catabolic genes. Restoring BCAA catabolic flux with a pharmacological inhibitor of BCKA dehydrogenase kinase (BCKDK) ( a suppressor of BCKA dehydrogenase) reduced the abundance of BCAA and BCKA and markedly attenuated IR in ob/ob mice. Similar outcomes were achieved by reducing protein (and thus BCAA) intake, whereas increasing BCAA intake did the opposite; this corroborates the pathogenic roles of BCAAs and BCKAs in IR in ob/ob mice. Like BCAAs, BCKAs also suppressed insulin signaling via activation of mammalian target of rapamycin complex 1. Finally, the small-molecule BCKDK inhibitor significantly attenuated IR in high-fat diet-induced obese mice. Collectively, these data demonstrate a pivotal causal role of a BCAA catabolic defect and elevated abundance of BCAAs and BCKAs in obesity-associated IR and provide proof-of-concept evidence for the therapeutic validity of manipulating BCAA metabolism for treating diabetes.

Gluconeogenic Enzymes in β-Cells: Pharmacological Targets for Improving Insulin Secretion

Pancreatic β-cells express the gluconeogenic enzymes glucose 6-phosphatase (G6Pase), fructose 1,6-bisphosphatase (FBP), and phosphoenolpyruvate (PEP) carboxykinase (PCK), which modulate glucose-stimulated insulin secretion (GSIS) through their ability to reverse otherwise irreversible glycolytic steps. Here, we review current knowledge about the expression and regulation of these enzymes in the context of manipulating them to improve insulin secretion in diabetics. Because the regulation of gluconeogenic enzymes in β-cells is so poorly understood, we propose novel research avenues to study these enzymes as modulators of insulin secretion and β-cell dysfunction, with especial attention to FBP, which constitutes an attractive target with an inhibitor under clinical evaluation at present.

GNAS gene is an important regulator of insulin secretory capacity in pancreatic β-cells

BACKGROUND

Type 2 diabetes (T2D) is a complex polygenic disease with unclear mechanism. In an attempt to identify novel genes involved in β-cell function, we harness a bioinformatics method called Loss-of-function tool (LoFtool) gene score.

METHODS

RNA-sequencing data from human islets were used to cross-reference genes within the 1st quartile of most intolerant LoFtool score with the 100th most expressed genes in human islets. Out of these genes, GNAS and EEF1A1 genes were selected for further investigation in diabetic islets, metabolic tissues along with their correlation with diabetic phenotypes. The influence of GNAS and EEF1A1 on insulin secretion and β-cell function were validated in INS-1 cells.

RESULTS

A comparatively higher expression level of GNAS and EEF1A1 was observed in human islets than fat, liver and muscle tissues. Furthermore, diabetic islets displayed a reduced expression of GNAS, but not of EEF1A, compared to non-diabetic islets. The expression of GNAS was positively correlated with insulin secretory index, GLP1R, GIPR and inversely correlated with HbA_1c. Diabetic human islets displayed a reduced cAMP generation and insulin secretory capacity in response to glucose. Moreover, siRNA silencing of GNAS in INS-1 cells reduced insulin secretion, insulin content, and cAMP production. In addition, the expression of Insulin, PDX1, and MAFA was significantly down-regulated in GNAS-silenced cells. However, cell viability and apoptosis rate were unaffected.

CONCLUSION

LoFtool is a powerful tool to identify genes associated with pancreatic islets dysfunction. GNAS is a crucial gene for the β-cell insulin secretory capacity.

Bioinformatic prediction of critical genes and pathways involved in longevity in Drosophila melanogaster

The pursuit of longevity has been the goal of humanity since ancient times. Genetic alterations have been demonstrated to affect lifespan. As increasing numbers of pro-longevity genes and anti-longevity genes have been discovered in Drosophila, screening for functionally important genes among the large number of genes has become difficult. The aim of the present study was to explore critical genes and pathways affecting longevity in Drosophila melanogaster. In this study, 168 genes associated with longevity in D. melanogaster were collected from the Human Ageing Genomic Resources (HAGR) database. Network clustering analysis, network topological analysis, and pathway analysis were integrated to identify key genes and pathways. Quantitative real-time PCR (qRT-PCR) was applied to verify the expression of genes in representative pathways and of predicted genes derived from the gene-gene sub-network. Our results revealed that six key pathways might be associated with longevity, including the longevity-regulating pathway, the peroxisome pathway, the mTOR-signalling pathway, the FOXO-signalling pathway, the AGE-RAGE-signalling pathway in diabetic complications, and the TGF-beta-signalling pathway. Moreover, the results revealed that six key genes in representative pathways, including Cat, Ry, S6k, Sod, Tor, and Tsc1, and the predicted genes Jra, Kay, and Rheb exhibited significant expression changes in ageing D. melanogaster strain w¹¹¹⁸ compared to young ones. Overall, our results revealed that six pathways and six key genes might play pivotal roles in regulating longevity, and three interacting genes might be implicated in longevity. The results will not only provide new insight into the mechanisms of longevity, but also provide novel ideas for network-based approaches for longevity-related research.

JAZF1, a relevant metabolic regulator in type 2 diabetes

Excessive adiposity and metabolic inflammation are the key risk factors of type 2 diabetes mellitus (T2DM). Juxtaposed with another zinc finger gene 1 (JAZF1) has been identified as a novel transcriptional cofactor, with function of regulating glucose and lipid homeostasis and inflammation. JAZF1 is involved in metabolic process of T2DM via interaction with several nuclear receptors and protein kinases. Additionally, increasing evidence from genome-wide association studies (GWAS) has shown that JAZF1 polymorphisms are closely associated with T2DM. In this review, we have updated the latest research advances on JAZF1 and discussed its regulatory network in T2DM. The association between JAZF1 polymorphisms and T2DM is discussed as well. The information provided is of importance for guiding future studies as well as for the design of JAZF1-based T2DM therapy.

Fine-scale haplotype mapping of MUT, AACS, SLC6A15 and PRKCA genes indicates association with insulin resistance of metabolic syndrome and relationship with branched chain amino acid metabolism or regulation

Branched chain amino acids (BCAA) are essential elements of the human diet, which display increased plasma levels in obesity and regained particular interest as potential biomarkers for development of diabetes. To define determinants of insulin resistance (IR) we investigated 73 genes involved in BCAA metabolism or regulation by fine-scale haplotype mapping in two European populations with metabolic syndrome. French and Romanians (n = 465) were genotyped for SNPs (Affymetrix) and enriched by imputation (BEAGLE 4.1) at 1000 genome project density. Initial association hits detected by sliding window were refined (HAPLOVIEW 3.1 and PHASE 2.1) and correlated to homeostasis model assessment (HOMAIR) index, in vivo insulin sensitivity (SI) and BCAA plasma levels (ANOVA). Four genomic regions were associated with IR located downstream of MUT, AACS, SLC6A15 and PRKCA genes (P between 9.3 and 3.7 x 10-5). Inferred haplotypes up to 13 SNPs length were associated with IR (e.g. MUT gene with P < 4.9 x 10-5; Bonferroni 1.3 x 10-3) and synergistic to HOMAIR. SNPs in the same regions were also associated with one order of magnitude lower P values (e.g. rs20167284 in the MUT gene with P < 1.27 x 10-4) and replicated in Mediterranean samples (n = 832). In French, influential haplotypes (OR > 2.0) were correlated with in vivo insulin sensitivity (1/SI) except for SLC6A15 gene. Association of these genes with BCAA levels was variable, but influential haplotypes confirmed implication of MUT from BCAA metabolism as well as a role of regulatory genes (AACS and PRKCA) and suggested potential changes in transcriptional activity. These data drive attention towards new regulatory regions involved in IR in relation with BCAA and show the ability of haplotypes in phased DNA to detect signals complimentary to SNPs, which may be useful in designing genetic markers for clinical applications in ethnic populations.

Aberrant fatty acid profile and FFAR4 signaling confer endocrine resistance in breast cancer

BACKGROUND

Evidence suggests that fatty acid receptor FFAR4 plays a tumor-promoting role in adipose tissue-adjacent malignancies, but its clinical relevance remains unexplored. Here, we investigated the clinical significance and underlying mechanisms of FFAR4 in hormone receptor-positive breast cancer (HRPBC).

METHODS

FFAR4 expression was assessed by immunohistochemistry in an exploration cohort of 307 breast cancer cases collected from two independent institutes. Two public breast cancer microarray datasets served as validation cohorts. Gas chromatography-mass spectrometry was employed to identify FFAR4 ligands in normal and cancerous breast tissues. Survival analyses were performed in all cohorts and designated molecular subgroups. Mechanistic studies were performed in vitro in hormone receptor-positive breast cancer cell lines MCF-7 and T-47D.

RESULTS

Aberrant FFAR4 expression and endogenous FFAR4 ligands were identified in breast cancer tissues, five FFAR4 ligands showed significantly elevated proportions in cancerous versus normal tissues. In the exploration cohort, FFAR4 was demonstrated as an independent prognostic factor for recurrences (HR: 2.183, 95% CI: 1.360-3.504, P = 0.001) and breast cancer-specific deaths (HR: 2.102, 95% CI: 1.173-3.766, P = 0.013) in HRPBC cases. In contrast, FFAR4 expression was not associated with prognosis in hormone receptor-negative cases. In the validation cohorts, FFAR4 mRNA levels were also observed to be associated with disease recurrence in estrogen receptor-positive cases, but not so in estrogen receptor-negative cases. FFAR4 activation by endogenous ligands and a synthetic ligand TUG891 significantly dampened tamoxifen's efficacy on HRPBC cells, whereas FFAR4 knockdown or antagonist AH7614 abrogated this effect. Furthermore, FFAR4-induced tamoxifen resistance was dependent on ERK and AKT pathways in HRPBC.

CONCLUSIONS

Our results establish a novel role of FFAR4 and its ligands in the complicated interactions between tissue lipid profile and cancer biology. FFAR4 signaling confers tamoxifen resistance in HRPBC cell line and FFAR4 expression can serve as a prognostic biomarker for tamoxifen-treated HRPBC patients. FFAR4 may serve as a potential target for anti-breast cancer therapies, especially in endocrine resistant cases.

GPR120 protects lipotoxicity-induced pancreatic β-cell dysfunction through regulation of PDX1 expression and inhibition of islet inflammation

G-protein coupled receptor 120 (GPR120) has been shown to act as an omega-3 unsaturated fatty acid sensor and is involved in insulin secretion. However, the underlying mechanism in pancreatic β cells remains unclear. To explore the potential link between GPR120 and β-cell function, its agonists docosahexaenoic acid (DHA) and GSK137647A were used in palmitic acid (PA)-induced pancreatic β-cell dysfunction, coupled with GPR120 knockdown (KD) in MIN6 cells and GPR120 knockout (KO) mice to identify the underlying signaling pathways. In vitro and ex vivo treatments of MIN6 cells and islets isolated from wild-type (WT) mice with DHA and GSK137647A restored pancreatic duodenal homeobox-1 (PDX1) expression levels and β-cell function via inhibiting PA-induced elevation of proinflammatory chemokines and activation of nuclear factor κB, c-Jun amino (N)-terminal kinases1/2 and p38MAPK signaling pathways. On the contrary, these GPR120 agonism-mediated protective effects were abolished in GPR120 KD cells and islets isolated from GPR120 KO mice. Furthermore, GPR120 KO mice displayed glucose intolerance and insulin resistance relative to WT littermates, and β-cell functional related genes were decreased while inflammation was exacerbated in islets with increased macrophages in pancreas from GPR120 KO mice. DHA and GSK137647A supplementation ameliorated glucose tolerance and insulin sensitivity, as well as improved Pdx1 expression and islet inflammation in diet-induced obese WT mice, but not in GPR120 KO mice. These findings indicate that GPR120 activation is protective against lipotoxicity-induced pancreatic β-cell dysfunction, via the mediation of PDX1 expression and inhibition of islet inflammation, and that GPR120 activation may serve as a preventative and therapeutic target for obesity and diabetes.

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.

Reduced Expression of PLCXD3 Associates With Disruption of Glucose Sensing and Insulin Signaling in Pancreatic β-Cells

Previous work has shown that reduced expression of PLCXD3, a member of the phosphoinositide-specific phospholipases (PI-PLC) family, impaired insulin secretion with an unclear mechanism. In the current study, we aim to investigate the mechanism underlying this effect using human islets and rat INS-1 (832/13) cells. Microarray and RNA sequencing data showed that PLCXD3 is among the highly expressed PI-PLCs in human islets and INS-1 (832/13) cells. Expression of PLCXD3 was reduced in human diabetic islets, correlated positively with Insulin and GLP1R expression and inversely with the donor's body mass index (BMI) and glycated hemoglobin (HbA_1c). Expression silencing of PLCXD3 in INS-1 (832/13) cells was found to reduce glucose-stimulated insulin secretion (GSIS) and insulin content. In addition, the expression of Insulin, NEUROD1, GLUT2, GCK, INSR, IRS2, and AKT was downregulated. Cell viability and apoptosis rate were unaffected. In conclusion, our data suggest that low expression of PLCXD3 in pancreatic β-cells associates with downregulation of the key insulin signaling and insulin biosynthesis genes as well as reduction in glucose sensing.

Chromatin interactions and expression quantitative trait loci reveal genetic drivers of multimorbidities

Clinical studies of non-communicable diseases identify multimorbidities that suggest a common set of predisposing factors. Despite the fact that humans have ~24,000 genes, we do not understand the genetic pathways that contribute to the development of multimorbid non-communicable disease. Here we create a multimorbidity atlas of traits based on pleiotropy of spatially regulated genes. Using chromatin interaction and expression Quantitative Trait Loci (eQTL) data, we analyse 20,782 variants (p < 5 × 10-6) associated with 1351 phenotypes to identify 16,248 putative spatial eQTL-eGene pairs that are involved in 76,013 short- and long-range regulatory interactions (FDR < 0.05) in different human tissues. Convex biclustering of spatial eGenes that are shared among phenotypes identifies complex interrelationships between nominally different phenotype-associated SNPs. Our approach enables the simultaneous elucidation of variant interactions with target genes that are drivers of multimorbidity, and those that contribute to unique phenotype associated characteristics.

Multi-omic analysis of signalling factors in inflammatory comorbidities

BACKGROUND

Inflammation is a core element of many different, systemic and chronic diseases that usually involve an important autoimmune component. The clinical phase of inflammatory diseases is often the culmination of a long series of pathologic events that started years before. The systemic characteristics and related mechanisms could be investigated through the multi-omic comparative analysis of many inflammatory diseases. Therefore, it is important to use molecular data to study the genesis of the diseases. Here we propose a new methodology to study the relationships between inflammatory diseases and signalling molecules whose dysregulation at molecular levels could lead to systemic pathological events observed in inflammatory diseases.

RESULTS

We first perform an exploratory analysis of gene expression data of a number of diseases that involve a strong inflammatory component. The comparison of gene expression between disease and healthy samples reveals the importance of members of gene families coding for signalling factors. Next, we focus on interested signalling gene families and a subset of inflammation related diseases with multi-omic features including both gene expression and DNA methylation. We introduce a phylogenetic-based multi-omic method to study the relationships between multi-omic features of inflammation related diseases by integrating gene expression, DNA methylation through sequence based phylogeny of the signalling gene families. The models of adaptations between gene expression and DNA methylation can be inferred from pre-estimated evolutionary relationship of a gene family. Members of the gene family whose expression or methylation levels significantly deviate from the model are considered as the potential disease associated genes.

CONCLUSIONS

Applying the methodology to four gene families (the chemokine receptor family, the TNF receptor family, the TGF- β gene family, the IL-17 gene family) in nine inflammation related diseases, we identify disease associated genes which exhibit significant dysregulation in gene expression or DNA methylation in the inflammation related diseases, which provides clues for functional associations between the diseases.

The adhesion receptor GPR56 is activated by extracellular matrix collagen III to improve β-cell function

AIMS

G-protein coupled receptor 56 (GPR56) is the most abundant islet-expressed G-protein coupled receptor, suggesting a potential role in islet function. This study evaluated islet expression of GPR56 and its endogenous ligand collagen III, and their effects on β-cell function.

METHODS

GPR56 and collagen III expression in mouse and human pancreas sections was determined by fluorescence immunohistochemistry. Effects of collagen III on β-cell proliferation, apoptosis, intracellular calcium ([Ca2+]i) and insulin secretion were determined by cellular BrdU incorporation, caspase 3/7 activities, microfluorimetry and radioimmunoassay, respectively. The role of GPR56 in islet vascularisation and innervation was evaluated by immunohistochemical staining for CD31 and TUJ1, respectively, in pancreases from wildtype (WT) and Gpr56-/- mice, and the requirement of GPR56 for normal glucose homeostasis was determined by glucose tolerance tests in WT and Gpr56-/- mice.

RESULTS

Immunostaining of mouse and human pancreases revealed that GPR56 was expressed by islet β-cells while collagen III was confined to the peri-islet basement membrane and islet capillaries. Collagen III protected β-cells from cytokine-induced apoptosis, triggered increases in [Ca2+]i and potentiated glucose-induced insulin secretion from WT islets but not from Gpr56-/- islets. Deletion of GPR56 did not affect glucose-induced insulin secretion in vitro and it did not impair glucose tolerance in adult mice. GPR56 was not required for normal islet vascularisation or innervation.

CONCLUSION

We have demonstrated that collagen III improves islet function by increasing insulin secretion and protecting against apoptosis. Our data suggest that collagen III may be effective in optimising islet function to improve islet transplantation outcomes, and GPR56 may be a target for the treatment of type 2 diabetes.

A type 2 diabetes disease module with a high collective influence for Cdk2 and PTPLAD1 is localized in endosomes

Despite the identification of many susceptibility genes our knowledge of the underlying mechanisms responsible for complex disease remains limited. Here, we identified a type 2 diabetes disease module in endosomes, and validate it for functional relevance on selected nodes. Using hepatic Golgi/endosomes fractions, we established a proteome of insulin receptor-containing endosomes that allowed the study of physical protein interaction networks on a type 2 diabetes background. The resulting collated network is formed by 313 nodes and 1147 edges with a topology organized around a few major hubs with Cdk2 displaying the highest collective influence. Overall, 88% of the nodes are associated with the type 2 diabetes genetic risk, including 101 new candidates. The Type 2 diabetes module is enriched with cytoskeleton and luminal acidification–dependent processes that are shared with secretion-related mechanisms. We identified new signaling pathways driven by Cdk2 and PTPLAD1 whose expression affects the association of the insulin receptor with TUBA, TUBB, the actin component ACTB and the endosomal sorting markers Rab5c and Rab11a. Therefore, the interactome of internalized insulin receptors reveals the presence of a type 2 diabetes disease module enriched in new layers of feedback loops required for insulin signaling, clearance and islet biology.

Reduced expression of stearoyl-CoA desaturase-1, but not free fatty acid receptor 2 or 4 in subcutaneous adipose tissue of patients with newly diagnosed type 2 diabetes mellitus

Background

In subcutaneous adipose tissue (SAT), higher stearoyl-CoA desaturase-1 (SCD1) expression has been related to improved insulin sensitivity in thiazolidinedione-treated type 2 diabetes mellitus patients. In animal models, deficiency of the free fatty acid receptor (FFAR) 2 associated with higher and FFAR4-deficiency with lower insulin sensitivity. We hypothesized that increased FFAR2 expression and reductions in FFAR4 and SCD1 expression in SAT of type 2 diabetes mellitus patients associate positively with insulin resistance and impaired beta cell function.

Methods

Twenty-five type 2 diabetes mellitus patients and 25 glucose-tolerant humans (CON) matched for sex, age, and BMI underwent mixed-meal tests to assess insulin sensitivity (OGIS) and beta cell function (ΔAUC(C-peptide)_{0–180 min}/ΔAUC(glucose)_{0–180 min}) in a cross-sectional study. Gene and protein expression of SCD1 and FFAR2/4 were quantified in SAT biopsies.

Results

Insulin sensitivity was 14% and beta cell function 71% (both p < 0.001) lower in type 2 diabetes mellitus patients. In type 2 diabetes mellitus, SCD1 mRNA was fivefold (p < 0.001) and protein expression twofold (p < 0.01) lower. While FFAR2/4 mRNA and protein expression did not differ between groups, FFAR2 protein levels correlated negatively with beta cell function only in CON (r = −0.74, p < 0.01). However, neither SCD1 nor FFAR2/4 protein expression correlated with insulin sensitivity in both groups.

Conclusions

Type 2 diabetes patients have lower SCD1, which does not associate with insulin resistance. Only in non-diabetic humans, FFAR2 associated with impaired beta cell function.

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 of branched chain amino acids related variant rs1440581 with risk of incident diabetes and longitudinal changes in insulin resistance in Chinese

AIMS

Previous genome-wide association studies reported rs1440581 was significantly associated with circulating branched chain amino acids (BCAAs) levels in Europeans. We aimed to investigate association of BCAAs related variant rs1440581 with incident T2D risk and longitudinal changes in glucose-related metabolic traits in a community-based prospective cohort of Chinese.

METHODS

6043 non-diabetic participants aged ≥ 40 years from a community-based population at baseline were included and followed-up for 5 years. The BCAAs related variant rs1440581 was genotyped. Incident T2D was defined as fasting plasma glucose (FPG) ≥ 7.0 mmol/L or taking anti-diabetic therapy. Anthropometry and biochemical measurements were evaluated at both baseline and follow-up.

RESULTS

576 (9.5%) participants developed T2D during the 5-year follow-up. Each C-allele was associated with a 20% higher risk of incident T2D (odds ratio = 1.20, 95% confidence interval [1.05, 1.36]) after adjustments for the confounders. We did not find a main effect of the variant on increase in fasting serum insulin (FSI) level or insulin resistance (IR). However, we found rs1440581 significantly modified effect of weight gain on increase in FSI and HOMA-IR. In the C-allele carriers, body mass index increase was associated with greater increase in Log₁₀_FSI (β ± SE 0.027 ± 0.002) and Log₁₀_HOMA-IR (0.030 ± 0.003), as compared to T-allele (both P for interaction = 0.003).

CONCLUSIONS

BCAAs related genetic variant rs1440581 was associated with an increased risk of incident T2D in a Chinese population. This variant might modify effect of weight gain on development in IR.

Molecular Biology and Physiology of Volume-Regulated Anion Channel (VRAC)

The Volume-Regulated Anion Channel (VRAC) is activated by cell swelling and plays a key role in cell volume regulation. VRAC is ubiquitously expressed in vertebrate cells and also implicated in many other physiological and cellular processes including fluid secretion, glutamate release, membrane potential regulation, cell proliferation, migration, and apoptosis. Although its biophysical properties have been well characterized, the molecular identity of VRAC remained a mystery for almost three decades. The field was transformed by recent discoveries showing that the leucine-rich repeat-containing protein 8A (LRRC8A, also named SWELL1) and its four other homologs form heteromeric VRAC channels. The composition of LRRC8 subunits determines channel properties and substrate selectivity of a large variety of different VRACs. Incorporating purified SWELL1-containing protein complexes into lipid bilayers is sufficient to reconstitute channel activities, a finding that supports the decrease in intracellular ionic strength as the mechanism of VRAC activation during cell swelling. Characterization of Swell1 knockout mice uncovers the important role of VRAC in T cell development, pancreatic β-cell glucose-stimulated insulin secretion, and adipocyte metabolic function. The ability to permeate organic osmolytes and metabolites is a major feature of VRAC. The list of VRAC substrates is expected to grow, now also including some cancer drugs and antibiotics even under non-cell swelling conditions. Therefore, a critical role of VRAC in drug resistance and cell-cell communication is emerging. This review summarizes the exciting recent progress on the structure-function relationship and physiology of VRAC and discusses key future questions to be solved.

Protective role of the ELOVL2/docosahexaenoic acid axis in glucolipotoxicity-induced apoptosis in rodent beta cells and human islets

AIMS/HYPOTHESIS

Dietary n-3 polyunsaturated fatty acids, especially docosahexaenoic acid (DHA), are known to influence glucose homeostasis. We recently showed that Elovl2 expression in beta cells, which regulates synthesis of endogenous DHA, was associated with glucose tolerance and played a key role in insulin secretion. The present study aimed to examine the role of the very long chain fatty acid elongase 2 (ELOVL2)/DHA axis on the adverse effects of palmitate with high glucose, a condition defined as glucolipotoxicity, on beta cells.

METHODS

We detected ELOVL2 in INS-1 beta cells and mouse and human islets using quantitative PCR and western blotting. Downregulation and adenoviral overexpression of Elovl2 was carried out in beta cells. Ceramide and diacylglycerol levels were determined by radio-enzymatic assay and lipidomics. Apoptosis was quantified using caspase-3 assays and poly (ADP-ribose) polymerase cleavage. Palmitate oxidation and esterification were determined by [U-¹⁴C]palmitate labelling.

RESULTS

We found that glucolipotoxicity decreased ELOVL2 content in rodent and human beta cells. Downregulation of ELOVL2 drastically potentiated beta cell apoptosis induced by glucolipotoxicity, whereas adenoviral Elovl2 overexpression and supplementation with DHA partially inhibited glucolipotoxicity-induced cell death in rodent and human beta cells. Inhibition of beta cell apoptosis by the ELOVL2/DHA axis was associated with a decrease in ceramide accumulation. However, the ELOVL2/DHA axis was unable to directly alter ceramide synthesis or metabolism. By contrast, DHA increased palmitate oxidation but did not affect its esterification. Pharmacological inhibition of AMP-activated protein kinase and etomoxir, an inhibitor of carnitine palmitoyltransferase 1 (CPT1), the rate-limiting enzyme in fatty acid β-oxidation, attenuated the protective effect of the ELOVL2/DHA axis during glucolipotoxicity. Downregulation of CPT1 also counteracted the anti-apoptotic action of the ELOVL2/DHA axis. By contrast, a mutated active form of Cpt1 inhibited glucolipotoxicity-induced beta cell apoptosis when ELOVL2 was downregulated.

CONCLUSIONS/INTERPRETATION

Our results identify ELOVL2 as a critical pro-survival enzyme for preventing beta cell death and dysfunction induced by glucolipotoxicity, notably by favouring palmitate oxidation in mitochondria through a CPT1-dependent mechanism.

ANRIL: A lncRNA at the CDKN2A/B Locus With Roles in Cancer and Metabolic Disease

The CDKN2A/B genomic locus is associated with risk of human cancers and metabolic disease. Although the locus contains several important protein-coding genes, studies suggest disease roles for a lesser-known antisense lncRNA encoded at this locus, called ANRIL. ANRIL is a complex gene containing at least 21 exons in simians, with many reported linear and circular isoforms. Like other genes, abundance of ANRIL is regulated by epigenetics, classic transcription regulation, splicing, and post-transcriptional influences such as RNA stability and microRNAs. Known molecular functions of ANRIL include in cis and in trans gene regulation through chromatin modification complexes, and influence over microRNA signaling networks. Polymorphisms at the ANRIL gene are linked to risk for many different cancers, as well as risk of atherosclerotic cardiovascular disease, bone mass, obesity and type 2 diabetes. A broad array of variable reported impacts of polymorphisms on ANRIL abundance, splicing and function suggests that ANRIL has cell-type and context-dependent regulation and actions. In cancer cells, ANRIL gain of function increases proliferation, metastasis, cell survival and epithelial-mesenchymal transformation, whereas ANRIL loss of function decreases tumor size and growth, invasion and metastasis, and increases apoptosis and senescence. In metabolic disease, polymorphisms at the ANRIL gene are linked to risk of type 2 diabetes, coronary artery disease, coronary artery calcium score, myocardial infarction, and stroke. Intriguingly, with the exception of one polymorphism in exon 2 of ANRIL, the single nucleotide polymorphisms (SNPs) associated with atherosclerosis and diabetes are non-overlapping. Evidence suggests that ANRIL gain of function increases atherosclerosis; in diabetes, a risk-SNP reduced the pancreatic beta cell proliferation index. Studies are limited by the uncertain relevance of rodent models to ANRIL studies, since most ANRIL exons do not exist in mouse. Diverse cell-type-dependent results suggest it is necessary to perform studies in the relevant primary human tissue for each disease. Much remains to be learned about the biology of ANRIL in human health and disease; this research area may lead to insight into disease mechanisms and therapeutic approaches.

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.

Fatty Acid-Stimulated Insulin Secretion vs. Lipotoxicity

Fatty acid (FA)-stimulated insulin secretion (FASIS) is reviewed here in contrast to type 2 diabetes etiology, resulting from FA overload, oxidative stress, intermediate hyperinsulinemia, and inflammation, all converging into insulin resistance. Focusing on pancreatic islet β-cells, we compare the physiological FA roles with the pathological ones. Considering FAs not as mere amplifiers of glucose-stimulated insulin secretion (GSIS), but as parallel insulin granule exocytosis inductors, partly independent of the K_ATP channel closure, we describe the FA initiating roles in the prediabetic state that is induced by retardations in the glycerol-3-phosphate (glucose)-promoted glycerol/FA cycle and by the impaired GPR40/FFA1 (free FA1) receptor pathway, specifically in its amplification by the redox-activated mitochondrial phospholipase, iPLA2γ. Also, excessive dietary FAs stimulate intestine enterocyte incretin secretion, further elevating GSIS, even at low glucose levels, thus contributing to diabetic hyperinsulinemia. With overnutrition and obesity, the FA overload causes impaired GSIS by metabolic dysbalance, paralleled by oxidative and metabolic stress, endoplasmic reticulum stress and numerous pro-apoptotic signaling, all leading to decreased β-cell survival. Lipotoxicity is exerted by saturated FAs, whereas ω-3 polyunsaturated FAs frequently exert antilipotoxic effects. FA-facilitated inflammation upon the recruitment of excess M1 macrophages into islets (over resolving M2 type), amplified by cytokine and chemokine secretion by β-cells, leads to an inevitable failure of pancreatic β-cells.

Impaired β-cell glucokinase as an underlying mechanism in diet-induced diabetes

High-fat diet (HFD)-fed mouse models have been widely used to study early type 2 diabetes. Decreased β-cell glucokinase (GCK) expression has been observed in HFD-induced diabetes. However, owing to its crucial roles in glucose metabolism in the liver and in islet β-cells, the contribution of decreased GCK expression to the development of HFD-induced diabetes is unclear. Here, we employed a β-cell-targeted gene transfer vector and determined the impact of β-cell-specific increase in GCK expression on β-cell function and glucose handling in vitro and in vivo Overexpression of GCK enhanced glycolytic flux, ATP-sensitive potassium channel activation and membrane depolarization, and increased proliferation in Min6 cells. β-cell-targeted GCK transduction did not change glucose handling in chow-fed C57BL/6 mice. Although adult mice fed a HFD showed reduced islet GCK expression, impaired glucose tolerance and decreased glucose-stimulated insulin secretion (GSIS), β-cell-targeted GCK transduction improved glucose tolerance and restored GSIS. Islet perifusion experiments verified restored GSIS in isolated HFD islets by GCK transduction. Thus, our data identify impaired β-cell GCK expression as an underlying mechanism for dysregulated β-cell function and glycemic control in HFD-induced diabetes. Our data also imply an etiological role of GCK in diet-induced diabetes.This article has an associated First Person interview with the first author of the paper.

Altered expression of WFS1 and NOTCH2 genes associated with diabetic nephropathy in T2DM patients

AIM

The increased incidence of type 2 diabetes mellitus (T2DM) and the importance of early identification and management of its complications, especially diabetic nephropathy (DN), have spotted the light on genetic factors that increase risk of T2DM and its related nephropathy. The present study aimed at investigating expression of (KCNJ11, ABCC8, JAZF1, WFS1, PPARG, NOTCH2 and EXOSC4) genes in peripheral blood of T2DM patients.

METHOD

The study included 30 non-complicated T2DM patients, 30 patients with DN and 40 healthy controls. Quantitative Real Time PCR Array was used to study gene expression.

RESULTS

NOTCH2 showed higher expression while KCNJ11, JAZF1, WFS1 and PPARG genes showed lower expression in DN patients compared to non-complicated patients. KCNJ11, JAZF1, WFS1, PPARG, and EXOSC4 expression showed significant negative correlation with microalbumin, while NOTCH2 expression was significantly positively correlated with microalbumin. AS regard HbA1c and studied genes expression, there was significant negative correlation between WFS1 expression and HbA1c, while NOTCH2, KCNJ11, JAZF1, PPARG, EXOSC4 expression didn't show significant correlation with HbA1c. Risk ratio of studied genes expression showed that WFS1 and NOTCH2 had highest risk ratio (30) and highest sensitivity and specificity, in relation to DN and they were the best predictors in the group of studied genes at cut off value of ≤0.861 for WFS1 and ≥0.678 for NOTCH2.

CONCLUSION

Altered expression of WFS1 and NOTCH2 genes may play a role in pathogenesis and development of DN in patients with T2DM. These results may contribute in early identification and management of DN.

Proteasomal degradation of the histone acetyl transferase p300 contributes to beta-cell injury in a diabetes environment

In type 2 diabetes, amyloid oligomers, chronic hyperglycemia, lipotoxicity, and pro-inflammatory cytokines are detrimental to beta-cells, causing apoptosis and impaired insulin secretion. The histone acetyl transferase p300, involved in remodeling of chromatin structure by epigenetic mechanisms, is a key ubiquitous activator of the transcriptional machinery. In this study, we report that loss of p300 acetyl transferase activity and expression leads to beta-cell apoptosis, and most importantly, that stress situations known to be associated with diabetes alter p300 levels and functional integrity. We found that proteasomal degradation is the mechanism subserving p300 loss in beta-cells exposed to hyperglycemia or pro-inflammatory cytokines. We also report that melatonin, a hormone produced in the pineal gland and known to play key roles in beta-cell health, preserves p300 levels altered by these toxic conditions. Collectively, these data imply an important role for p300 in the pathophysiology of diabetes.

Increased secreted frizzled-related protein 4 and ficolin-3 levels in gestational diabetes mellitus women

By biochemical and epidemiological similarity with type 2 diabetes mellitus (T2DM), gestational diabetes mellitus (GDM) has some overlap between prediction markers and risk factors of T2DM. The present study aimed to establish that secreted frizzled-related protein 4 (SFRP4) and ficolin-3 levels, which have been linked to insulin resistance and the development of T2DM, are elevated in GDM women. A longitudinal prospective cohort study of 86 GDM and 273 normal glucose tolerant (NGT) pregnant women was performed. The clinical parameters, lipid profiles, and serum SFRP4 and ficolin-3 levels were tested during the early and late second-trimester and third-trimester of pregnancy. Both SFRP4 and ficolin-3 levels were significantly higher in GDM women as compared to the NGT participants at three test points (p < 0.01). Spearman's correlation analysis showed that serum SFRP4 levels were significantly positively correlated with ficolin-3 during the early and late second-trimester and third-trimester of pregnancy. The elevated SFRP4 and ficolin-3 concentrations at 16-18 weeks gestation significantly associated with GDM were conformed using binary logistic regression analysis after controlling for other variables [odds ratios (OR) with 95% confidence intervals (CI) for SFRP4: 2.84 (1.78-4.53), p < 0.01; for ficolin-3: 2.45 (1.55-3.88), p < 0.01]. In Conclusions, increased SFRP4 and ficolin-3 levels are significantly associated with GDM development and might be important risk factors for this pregnancy complication.

Polymorphism of the GLIS3 gene in a Caucasian population and among individuals with carbohydrate metabolism disorders in Russia

Objective

Earlier, GLIS3 gene polymorphisms have been shown to be associated with the development of maturity onset diabetes of the young (MODY). We screened GLIS3 gene sequences among patients with MODY to identify probably pathogenic variants by whole-exome sequencing. We estimated frequency of rare single-nucleotide variants in the coding region of GLIS3 in a Caucasian population and among individuals with carbohydrate metabolism disorders in Russia.

Results

We identified 15 single-nucleotide variants in GLIS3. Three rare variants (minor allele frequency < 1%) rs806052, rs143051164, and rs149840771 were genotyped in 126 cases of MODY, in 188 patients with type 2 diabetes mellitus (DM2), and 564 randomly selected Caucasian individuals in Russia. A heterozygous rs806052 variant was identified in one patient with DM2; c.1270T frequency was 0.003. Prevalence of rs143051164 c.844G was 0.003 in the control population and 0.004 and 0.003 in MODY and DM2 samples, respectively. Prevalence of rs149840771 c.2096A was 0.003 and 0.004 in the control population and among MODY patients, respectively. In DM2 patients, rs149840771 c.2096A was not identified. We did not detect any associations of rs806052, rs143051164, and rs149840771 with carbohydrate metabolism disorders among patients with MODY and DM2 in Russia.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-3338-1) contains supplementary material, which is available to authorized users.

BRANCHED CHAIN AMINO ACIDS AT THE EDGE BETWEEN MENDELIAN AND COMPLEX DISORDERS

Branched chained amino acids (BCAA) are essential components of the human diet and important nutrient signals, which regain particular interest in recent years with the avenue of metabolomics studies suggesting their potential role as biomarkers. There is now compelling evidence for predictive role of BCAA in progression of diabetes, but causality relationship is still debated concerning insulin resistance and genetic versus non-genetic pathogenesis. Mendelian randomization studies in large cohorts of diabetes indicated pathogenic role of PPM1K (protein phosphatase Mg2+/Mn2+ dependent 1K) on Chr 4q22.1 gene, encoding for a phosphatase that activates BCKDH (branched chain keto acid dehydrogenase) complex. Recent studies indicated that insulin rapidly and dose-dependently regulates gene expression of the same complex, but the relationship with systemic insulin resistance and glucose levels is complex. Rare genetic syndromes due to Mendelian mutations in key genes in BCAA catabolism may be good models to understand potential role of gene of BCAA catabolism. However, in studying complex disorders geneticists are faced to complete new aspects of metabolic regulation complicating understanding genetics of obesity, diabetes or metabolic syndrome. A review of genetic syndromes of BCAA metabolism suggests that insulin resistance is not present, except rare cases of methylmalonic aciduria due to MUT (methylmalonyl-coA mutase) gene on Chr 6p12.3. Another aspect that complicates understanding is the new role of central nervous system (CNS) in insulin resistance. For a long time the hypothalamic hunger/satiety neuronal system was considered a key site of nutrient regulation. Genes may also affect the brain rewarding system (BRS) that would regulate food intake by modulating the motivation to obtain food and considering hedonic properties. Nutrigenomic and nutrigenetic investigations taking into account concurrently BCAA intake, metabolic regulation and gene variation have large perspectives to merge genetic and nutritional understanding in complex disorders.

Genetic Drivers of Pancreatic Islet Function

The majority of gene loci that have been associated with type 2 diabetes play a role in pancreatic islet function. To evaluate the role of islet gene expression in the etiology of diabetes, we sensitized a genetically diverse mouse population with a Western diet high in fat (45% kcal) and sucrose (34%) and carried out genome-wide association mapping of diabetes-related phenotypes. We quantified mRNA abundance in the islets and identified 18,820 expression QTL. We applied mediation analysis to identify candidate causal driver genes at loci that affect the abundance of numerous transcripts. These include two genes previously associated with monogenic diabetes (PDX1 and HNF4A), as well as three genes with nominal association with diabetes-related traits in humans (FAM83E, IL6ST, and SAT2). We grouped transcripts into gene modules and mapped regulatory loci for modules enriched with transcripts specific for α-cells, and another specific for δ-cells. However, no single module enriched for β-cell-specific transcripts, suggesting heterogeneity of gene expression patterns within the β-cell population. A module enriched in transcripts associated with branched-chain amino acid metabolism was the most strongly correlated with physiological traits that reflect insulin resistance. Although the mice in this study were not overtly diabetic, the analysis of pancreatic islet gene expression under dietary-induced stress enabled us to identify correlated variation in groups of genes that are functionally linked to diabetes-associated physiological traits. Our analysis suggests an expected degree of concordance between diabetes-associated loci in the mouse and those found in human populations, and demonstrates how the mouse can provide evidence to support nominal associations found in human genome-wide association mapping.

Systems biology of the IMIDIA biobank from organ donors and pancreatectomised patients defines a novel transcriptomic signature of islets from individuals with type 2 diabetes

AIMS/HYPOTHESIS

Pancreatic islet beta cell failure causes type 2 diabetes in humans. To identify transcriptomic changes in type 2 diabetic islets, the Innovative Medicines Initiative for Diabetes: Improving beta-cell function and identification of diagnostic biomarkers for treatment monitoring in Diabetes (IMIDIA) consortium ( www.imidia.org ) established a comprehensive, unique multicentre biobank of human islets and pancreas tissues from organ donors and metabolically phenotyped pancreatectomised patients (PPP).

METHODS

Affymetrix microarrays were used to assess the islet transcriptome of islets isolated either by enzymatic digestion from 103 organ donors (OD), including 84 non-diabetic and 19 type 2 diabetic individuals, or by laser capture microdissection (LCM) from surgical specimens of 103 PPP, including 32 non-diabetic, 36 with type 2 diabetes, 15 with impaired glucose tolerance (IGT) and 20 with recent-onset diabetes (<1 year), conceivably secondary to the pancreatic disorder leading to surgery (type 3c diabetes). Bioinformatics tools were used to (1) compare the islet transcriptome of type 2 diabetic vs non-diabetic OD and PPP as well as vs IGT and type 3c diabetes within the PPP group; and (2) identify transcription factors driving gene co-expression modules correlated with insulin secretion ex vivo and glucose tolerance in vivo. Selected genes of interest were validated for their expression and function in beta cells.

RESULTS

Comparative transcriptomic analysis identified 19 genes differentially expressed (false discovery rate ≤0.05, fold change ≥1.5) in type 2 diabetic vs non-diabetic islets from OD and PPP. Nine out of these 19 dysregulated genes were not previously reported to be dysregulated in type 2 diabetic islets. Signature genes included TMEM37, which inhibited Ca²⁺-influx and insulin secretion in beta cells, and ARG2 and PPP1R1A, which promoted insulin secretion. Systems biology approaches identified HNF1A, PDX1 and REST as drivers of gene co-expression modules correlated with impaired insulin secretion or glucose tolerance, and 14 out of 19 differentially expressed type 2 diabetic islet signature genes were enriched in these modules. None of these signature genes was significantly dysregulated in islets of PPP with impaired glucose tolerance or type 3c diabetes.

CONCLUSIONS/INTERPRETATION

These studies enabled the stringent definition of a novel transcriptomic signature of type 2 diabetic islets, regardless of islet source and isolation procedure. Lack of this signature in islets from PPP with IGT or type 3c diabetes indicates differences possibly due to peculiarities of these hyperglycaemic conditions and/or a role for duration and severity of hyperglycaemia. Alternatively, these transcriptomic changes capture, but may not precede, beta cell failure.

Transcriptional Regulatory Mechanisms in Adipose and Muscle Tissue Associated with Composite Glucometabolic Phenotypes

OBJECTIVE

Tissue-specific gene expression is associated with individual metabolic measures. However, these measures may not reflect the true but latent underlying biological phenotype. This study reports gene expression associations with multidimensional glucometabolic characterizations of obesity, glucose homeostasis, and lipid traits.

METHODS

Factor analysis was computed by using orthogonal rotation to construct composite phenotypes (CPs) from 23 traits in 256 African Americans without diabetes. Genome-wide transcript expression data from adipose and muscle were tested for association with CPs, and expression quantitative trait loci (eQTLs) were identified by associations between cis-acting single-nucleotide polymorphisms (SNPs) and gene expression.

RESULTS

The factor analysis identified six CPs. CPs 1 through 6 individually explained 34%, 12%, 9%, 8%, 6%, and 5% of the variation in 23 glucometabolic traits studied. There were 3,994 and 929 CP-associated transcripts identified in adipose and muscle tissue, respectively; CP2 had the largest number of associated transcripts. Pathway analysis identified multiple canonical pathways from the CP-associated transcripts. In adipose and muscle, significant cis-eQTLs were identified for 558 and 164 CP-associated transcripts (q-value < 0.01), respectively.

CONCLUSIONS

Adipose and muscle transcripts comprehensively define pathways involved in regulating glucometabolic disorders. Cis-eQTLs for CP-associated genes may act as primary causal determinants of glucometabolic phenotypes by regulating transcription of key genes.

Transcriptomic Analysis Reveals Novel Mechanisms Mediating Islet Dysfunction in the Intrauterine Growth-Restricted Rat

Intrauterine growth restriction (IUGR) increases the risk of type 2 diabetes developing in adulthood. In previous studies that used bilateral uterine artery ligation in a rat model of IUGR, age-associated decline in glucose homeostasis and islet function was revealed. To elucidate mechanisms contributing to IUGR pathogenesis, the islet transcriptome was sequenced from 2-week-old rats, when in vivo glucose tolerance is mildly impaired, and at 10 weeks of age, when rats are hyperglycemic and have reduced β-cell mass. RNA sequencing and functional annotation with Ingenuity Pathway Analysis revealed temporal changes in IUGR islets. For instance, gene expression involving amino acid metabolism was significantly reduced primarily at 2 weeks of age, but ion channel expression, specifically that involved in cell-volume regulation, was more disrupted in adult IUGR islets. Additionally, we observed alterations in the microenvironment of IUGR islets with extracellular matrix genes being significantly increased at 2 weeks of age and significantly decreased at 10 weeks. Specifically, hyaluronan synthase 2 expression and hyaluronan staining were increased in IUGR islets at 2 weeks of age (P < 0.05). Mesenchymal stromal cell-derived factors that have been shown to preserve islet allograft function, such as Anxa1, Cxcl12, and others, also were increased at 2 weeks and decreased in adult islets. Finally, comparisons of differentially expressed genes with those of type 2 diabetic human islets support a role for these pathways in human patients with diabetes. Together, these data point to new mechanisms in the pathogenesis of IUGR-mediated islet dysfunction in type 2 diabetes.

Mechanisms of β-cell dedifferentiation in diabetes: recent findings and future research directions

Like all the cells of an organism, pancreatic β-cells originate from embryonic stem cells through a complex cellular process termed differentiation. Differentiation involves the coordinated and tightly controlled activation/repression of specific effectors and gene clusters in a time-dependent fashion thereby giving rise to particular morphological and functional cellular features. Interestingly, cellular differentiation is not a unidirectional process. Indeed, growing evidence suggests that under certain conditions, mature β-cells can lose, to various degrees, their differentiated phenotype and cellular identity and regress to a less differentiated or a precursor-like state. This concept is termed dedifferentiation and has been proposed, besides cell death, as a contributing factor to the loss of functional β-cell mass in diabetes. β-cell dedifferentiation involves: (1) the downregulation of β-cell-enriched genes, including key transcription factors, insulin, glucose metabolism genes, protein processing and secretory pathway genes; (2) the concomitant upregulation of genes suppressed or expressed at very low levels in normal β-cells, the β-cell forbidden genes; and (3) the likely upregulation of progenitor cell genes. These alterations lead to phenotypic reconfiguration of β-cells and ultimately defective insulin secretion. While the major role of glucotoxicity in β-cell dedifferentiation is well established, the precise mechanisms involved are still under investigation. This review highlights the identified molecular mechanisms implicated in β-cell dedifferentiation including oxidative stress, endoplasmic reticulum (ER) stress, inflammation and hypoxia. It discusses the role of Foxo1, Myc and inhibitor of differentiation proteins and underscores the emerging role of non-coding RNAs. Finally, it proposes a novel hypothesis of β-cell dedifferentiation as a potential adaptive mechanism to escape cell death under stress conditions.

Integrative network analysis highlights biological processes underlying GLP-1 stimulated insulin secretion: A DIRECT study

Glucagon-like peptide 1 (GLP-1) stimulated insulin secretion has a considerable heritable component as estimated from twin studies, yet few genetic variants influencing this phenotype have been identified. We performed the first genome-wide association study (GWAS) of GLP-1 stimulated insulin secretion in non-diabetic individuals from the Netherlands Twin register (n = 126). This GWAS was enhanced using a tissue-specific protein-protein interaction network approach. We identified a beta-cell protein-protein interaction module that was significantly enriched for low gene scores based on the GWAS P-values and found support at the network level in an independent cohort from Tübingen, Germany (n = 100). Additionally, a polygenic risk score based on SNPs prioritized from the network was associated (P < 0.05) with glucose-stimulated insulin secretion phenotypes in up to 5,318 individuals in MAGIC cohorts. The network contains both known and novel genes in the context of insulin secretion and is enriched for members of the focal adhesion, extracellular-matrix receptor interaction, actin cytoskeleton regulation, Rap1 and PI3K-Akt signaling pathways. Adipose tissue is, like the beta-cell, one of the target tissues of GLP-1 and we thus hypothesized that similar networks might be functional in both tissues. In order to verify peripheral effects of GLP-1 stimulation, we compared the transcriptome profiling of ob/ob mice treated with liraglutide, a clinically used GLP-1 receptor agonist, versus baseline controls. Some of the upstream regulators of differentially expressed genes in the white adipose tissue of ob/ob mice were also detected in the human beta-cell network of genes associated with GLP-1 stimulated insulin secretion. The findings provide biological insight into the mechanisms through which the effects of GLP-1 may be modulated and highlight a potential role of the beta-cell expressed genes RYR2, GDI2, KIAA0232, COL4A1 and COL4A2 in GLP-1 stimulated insulin secretion.

Future Roadmaps for Precision Medicine Applied to Diabetes: Rising to the Challenge of Heterogeneity

Precision medicine, the concept that specific treatments can be targeted to groups of individuals with specific genetic, cellular, or molecular features, is a key aspect of modern healthcare, and its use is rapidly expanding. In diabetes, the application of precision medicine has been demonstrated in monogenic disease, where sulphonylureas are used to treat patients with neonatal diabetes due to mutations in ATP-dependent potassium (K_ATP) channel genes. However, diabetes is highly heterogeneous, both between and within polygenic and monogenic subtypes. Making the correct diagnosis and using the correct treatment from diagnosis can be challenging for clinicians, but it is crucial to prevent long-term morbidity and mortality. To facilitate precision medicine in diabetes, research is needed to develop a better understanding of disease heterogeneity and its impact on potential treatments for specific subtypes. Animal models have been used in diabetes research, but they are not translatable to humans in the majority of cases. Advances in molecular genetics and functional laboratory techniques and availability and sharing of large population data provide exciting opportunities for human studies. This review will map the key elements of future diabetes research in humans and its potential for clinical translation to promote precision medicine in all diabetes subtypes.

HbA1c is associated with altered expression in blood of cell cycle- and immune response-related genes

AIMS/HYPOTHESIS

Individuals with type 2 diabetes are heterogeneous in their glycaemic control as tracked by blood HbA_1c levels. Here, we investigated the extent to which gene expression levels in blood reflect current and future HbA_1c levels.

METHODS

HbA_1c levels at baseline and 1 and 2 year follow-up were compared with gene expression levels in 391 individuals with type 2 diabetes from the Hoorn Diabetes Care System Cohort (15,564 genes, RNA sequencing). The functions of associated baseline genes were investigated further using pathway enrichment analysis. Using publicly available data, we investigated whether the genes identified are also associated with HbA_1c in the target tissues, muscle and pancreas.

RESULTS

At baseline, 220 genes (1.4%) were associated with baseline HbA_1c. Identified genes were enriched for cell cycle and complement system activation pathways. The association of 15 genes extended to the target tissues, muscle (n = 113) and pancreatic islets (n = 115). At follow-up, expression of 25 genes (0.16%) associated with 1 year HbA_1c and nine genes (0.06%) with 2 year HbA_1c. Five genes overlapped across all time points, and 18 additional genes between baseline and 1 year follow-up. After adjustment for baseline HbA_1c, the number of significant genes at 1 and 2 years markedly decreased, suggesting that gene expression levels in whole blood reflect the current glycaemic state and but not necessarily the future glycaemic state.

CONCLUSIONS/INTERPRETATION

HbA_1c levels in individuals with type 2 diabetes are associated with expression levels of genes that link to the cell cycle and complement system activation.

Pancreatic β-Cell Electrical Activity and Insulin Secretion: Of Mice and Men

The pancreatic β-cell plays a key role in glucose homeostasis by secreting insulin, the only hormone capable of lowering the blood glucose concentration. Impaired insulin secretion results in the chronic hyperglycemia that characterizes type 2 diabetes (T2DM), which currently afflicts >450 million people worldwide. The healthy β-cell acts as a glucose sensor matching its output to the circulating glucose concentration. It does so via metabolically induced changes in electrical activity, which culminate in an increase in the cytoplasmic Ca2+ concentration and initiation of Ca2+-dependent exocytosis of insulin-containing secretory granules. Here, we review recent advances in our understanding of the β-cell transcriptome, electrical activity, and insulin exocytosis. We highlight salient differences between mouse and human β-cells, provide models of how the different ion channels contribute to their electrical activity and insulin secretion, and conclude by discussing how these processes become perturbed in T2DM.

Identification of key gene pathways and coexpression networks of islets in human type 2 diabetes

PURPOSE

The number of people with type 2 diabetes (T2D) is growing rapidly worldwide. Islet β-cell dysfunction and failure are the main causes of T2D pathological processes. The aim of this study was to elucidate the underlying pathways and coexpression networks in T2D islets.

MATERIALS AND METHODS

We analyzed the differentially expressed genes (DEGs) in the data set GSE41762, which contained 57 nondiabetic and 20 diabetic samples, and developed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein-protein interaction (PPI) network, the modules from the PPI network, and the gene annotation enrichment of modules were analyzed as well. Moreover, a weighted correlation network analysis (WGCNA) was applied to screen critical gene modules and coexpression networks and explore the biological significance.

RESULTS

We filtered 957 DEGs in T2D islets. Then GO and KEGG analyses identified that key pathways like inflammatory response, type B pancreatic cell differentiation, and calcium ion-dependent exocytosis were involved in human T2D. Three significant modules were filtered from the PPI network. Ribosome biogenesis, extrinsic apoptotic signaling pathway, and membrane depolarization during action potential were associated with the modules, respectively. Furthermore, coexpression network analysis by WGCNA identified 13 distinct gene modules of T2D islets and revealed four modules, which were strongly correlated with T2D and T2D biomarker hemoglobin A1c (HbA1c). Functional annotation showed that these modules mainly enriched KEGG pathways such as NF-kappa B signaling pathway, tumor necrosis factor signaling pathway, cyclic adenosine monophosphate signaling pathway, and peroxisome proliferators-activated receptor signaling pathway.

CONCLUSION

The results provide potential gene pathways and underlying molecular mechanisms for the prevention, diagnosis, and treatment of T2D.

Polymorphisms in FFAR4 (GPR120) Gene Modulate Insulin Levels and Sensitivity after Fish Oil Supplementation

The objective was to test whether FFAR4 single nucleotide polymorphisms (SNPs) are associated with glycemic control-related traits in humans following fish oil supplementation. A total of 210 participants were given 3 g/day of omega-3 (n-3) fatty acids (FA) (1.9–2.2 g of eicosapentaenoic acid (EPA) and 1.1 g of docosahexaenoic acid (DHA)) during six weeks. Biochemical parameters were taken before and after the supplementation. Using the HapMap database and the tagger procedure in Haploview, 12 tagging SNPs in FFAR4 were selected and then genotyped using TaqMan technology. Transcript expression levels were measured for 30 participants in peripheral mononuclear blood cells. DNA methylation levels were measured for 35 participants in leukocytes. In silico analyses were also performed. Four gene–diet interactions on fasting insulin levels and homeostatic model assessment of insulin resistance (HOMA-IR) index values were found. rs17108973 explained a significant proportion of the variance of insulin levels (3.0%) and HOMA-IR (2.03%) index values. Splice site prediction was different depending on the allele for rs11187527. rs17108973 and rs17484310 had different affinity for transcription factors depending on the allele. n-3 FAs effectively improve insulin-related traits for major allele homozygotes of four FFAR4 SNPs as opposed to carriers of the minor alleles.

FFA4 (GPR120) as a fatty acid sensor involved in appetite control, insulin sensitivity and inflammation regulation

Unsaturated long-chain fatty acids have been suggested to be beneficial in the context of cardiovascular disorders based in epidemiologic studies conducted in Greenland and Mediterranean. DHA and EPA are omega-3 polyunsaturated fatty acids that are plentiful in fish oil, and oleic acid is an omega-9 monounsaturated fatty acid, rich in olive oil. Dietary intake of these unsaturated long-chain fatty acids have been associated with insulin sensitivity and weight loss, which contrasts with the impairment of insulin sensitivity and weight gain associated with high intakes of saturated long-chain fatty acids. The recent discovery that free fatty acid receptor 4 (FFA4, also known as GPR120) acts as a sensor for unsaturated long-chain fatty acids started to unveil the molecular mechanisms underlying the beneficial functions played by these unsaturated long-chain fatty acids in various physiological processes, which include the secretions of gastrointestinal peptide hormones and glucose homeostasis. In this review, the physiological roles and therapeutic significance of FFA4 in appetite control, insulin sensitization, and inflammation reduction are discussed in relation to obesity and type 2 diabetes from pharmacological viewpoints.

Expression of the ZIP/SLC39A transporters in β-cells: a systematic review and integration of multiple datasets

BACKGROUND

Pancreatic β-cells require a constant supply of zinc to maintain normal insulin secretory function. Following co-exocytosis with insulin, zinc is replenished via the Zrt- and Irt-like (ZIP; SLC39A) family of transporters. However the ZIP paralogues of particular importance for zinc uptake, and associations with β-cell function and Type 2 Diabetes remain largely unexplored. We retrieved and statistically analysed publically available microarray and RNA-seq datasets to perform a systematic review on the expression of β-cell SLC39A paralogues. We complemented results with experimental data on expression profiling of human islets and mouse β-cell derived MIN6 cells, and compared transcriptomic and proteomic sequence conservation between human, mouse and rat.

RESULTS

The 14 ZIP paralogues have 73-98% amino sequence conservation between human and rodents. We identified 18 datasets for β-cell SLC39A analysis, which compared relative expression to non-β-cells, and expression in response to PDX-1 activity, cytokines, glucose and type 2 diabetic status. Published expression data demonstrate enrichment of transcripts for ZIP7 and ZIP9 transporters within rodent β-cells and of ZIP6, ZIP7 and ZIP14 within human β-cells, with ZIP1 most differentially expressed in response to cytokines and PDX-1 within rodent, and ZIP6 in response to diabetic status in human and glucose in rat. Our qPCR expression profiling data indicate that SLC39A6, -9, -13, and - 14 are the highest expressed paralogues in human β-cells and Slc39a6 and -7 in MIN6 cells.

CONCLUSIONS

Our systematic review, expression profiling and sequence alignment reveal similarities and potentially important differences in ZIP complements between human and rodent β-cells. We identify ZIP6, ZIP7, ZIP9, ZIP13 and ZIP14 in human and rodent and ZIP1 in rodent as potentially biologically important for β-cell zinc trafficking. We propose ZIP6 and ZIP7 are key functional orthologues in human and rodent β-cells and highlight these zinc importers as important targets for exploring associations between zinc status and normal physiology of β-cells and their decline in Type 2 Diabetes.

Molecular genetics of the transcription factor GLIS3 identifies its dual function in beta cells and neurons

The GLIS family zinc finger 3 isoform (GLIS3) is a risk gene for Type 1 and Type 2 diabetes, glaucoma and Alzheimer's disease endophenotype. We identified GLIS3 binding sites in insulin secreting cells (INS1) (FDR q<0.05; enrichment range 1.40-9.11 fold) sharing the motif wrGTTCCCArTAGs, which were enriched in genes involved in neuronal function and autophagy and in risk genes for metabolic and neuro-behavioural diseases. We confirmed experimentally Glis3-mediated regulation of the expression of genes involved in autophagy and neuron function in INS1 and neuronal PC12 cells. Naturally-occurring coding polymorphisms in Glis3 in the Goto-Kakizaki rat model of type 2 diabetes were associated with increased insulin production in vitro and in vivo, suggestive alteration of autophagy in PC12 and INS1 and abnormal neurogenesis in hippocampus neurons. Our results support biological pleiotropy of GLIS3 in pathologies affecting β-cells and neurons and underline the existence of trans‑nosology pathways in diabetes and its co-morbidities.