Swedish Research Council, European Research Council, Vinnova, Academy of Finland, Novo Nordisk Foundation, Scania University Hospital, Sigrid Juselius Foundation, Innovative Medicines Initiative 2 Joint Undertaking, Vasa Hospital district, Jakobstadsnejden Heart Foundation, Folkhälsan Research Foundation, Ollqvist Foundation, and Swedish Foundation for Strategic Research.

Genetic variation can modulate gene expression, and thereby phenotypic variation and susceptibility to complex diseases such as type 2 diabetes (T2D). Here we harnessed the potential of DNA and RNA sequencing in human pancreatic islets from 89 deceased donors to identify genes of potential importance in the pathogenesis of T2D. We present a catalog of genetic variants regulating gene expression (eQTL) and exon use (sQTL), including many long noncoding RNAs, which are enriched in known T2D-associated loci. Of 35 eQTL genes, whose expression differed between normoglycemic and hyperglycemic individuals, siRNA of tetraspanin 33 (TSPAN33), 5′-nucleotidase, ecto (NT5E), transmembrane emp24 protein transport domain containing 6 (TMED6), and p21 protein activated kinase 7 (PAK7) in INS1 cells resulted in reduced glucose-stimulated insulin secretion. In addition, we provide a genome-wide catalog of allelic expression imbalance, which is also enriched in known T2D-associated loci. Notably, allelic imbalance in paternally expressed gene 3 (PEG3) was associated with its promoter methylation and T2D status. Finally, RNA editing events were less common in islets than previously suggested in other tissues. Taken together, this study provides new insights into the complexity of gene regulation in human pancreatic islets and better understanding of how genetic variation can influence glucose metabolism.T ype 2 diabetes (T2D) is an increasing global health problem (1). Although genome-wide association studies (GWAS) have yielded more than 70 loci associated with T2D or related traits (2, 3), they have not provided the expected breakthrough in our understanding of the pathogenesis of the disease. They have nonetheless pointed at a central role of the pancreatic islets and β-cell dysfunction in the development of the disease (4, 5). It therefore seems pertinent to focus on human pancreatic islets to obtain insights into the molecular mechanisms causing the disease (6, 7). Given that most SNPs associated with T2D lie in noncoding regions, the majority of causal variants are likely to regulate gene expression rather than protein function per se. Therefore, combination of DNA and RNA sequencing in the same individuals may help to disentangle the role these SNPs play in the pathogenesis of the disease (8). Although the human pancreatic islet transcriptome has been previously described (6, 9-18), using microarrays or RNA sequencing of a limited number of nondiabetic individuals, this has not allowed a more global analysis of the complexity of the islet transcriptome in T2D. Here we combined genotypic imputation, expression microarrays, and exome and RNA sequencing (ExomeSeq and RNA-Seq) in a large number of human pancreatic islets from deceased donors with and without T2D. This study identified a number of novel genes, including long intergenic noncoding RNAs (lincRNAs), whose expression and/or splicing influences insulin secretion and is associated with glycemia. In addition, we provide a catalog of RNA editing and allele-specific expr...

We performed fine-mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in/near KCNQ1. “Credible sets” of variants most likely to drive each distinct signal mapped predominantly to non-coding sequence, implying that T2D association is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine-mapping implicated rs10830963 as driving T2D association. We confirmed that this T2D-risk allele increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D-risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease.

Type 2 diabetes (T2D) is a global pandemic. Genome-wide association studies (GWASs) have identified >100 genetic variants associated with the disease, including a common variant in the melatonin receptor 1 b gene (MTNR1B). Here, we demonstrate increased MTNR1B expression in human islets from risk G-allele carriers, which likely leads to a reduction in insulin release, increasing T2D risk. Accordingly, in insulin-secreting cells, melatonin reduced cAMP levels, and MTNR1B overexpression exaggerated the inhibition of insulin release exerted by melatonin. Conversely, mice with a disruption of the receptor secreted more insulin. Melatonin treatment in a human recall-by-genotype study reduced insulin secretion and raised glucose levels more extensively in risk G-allele carriers. Thus, our data support a model where enhanced melatonin signaling in islets reduces insulin secretion, leading to hyperglycemia and greater future risk of T2D. The findings also imply that melatonin physiologically serves to inhibit nocturnal insulin release.

The mucosal immune system identifies and fights invading pathogens, while allowing non-pathogenic organisms to persist. Mechanisms of pathogen/non-pathogen discrimination are poorly understood, as is the contribution of human genetic variation in disease susceptibility. We describe here a new, IRF3-dependent signaling pathway that is critical for distinguishing pathogens from normal flora at the mucosal barrier. Following uropathogenic E. coli infection, Irf3−/− mice showed a pathogen-specific increase in acute mortality, bacterial burden, abscess formation and renal damage compared to wild type mice. TLR4 signaling was initiated after ceramide release from glycosphingolipid receptors, through TRAM, CREB, Fos and Jun phosphorylation and p38 MAPK-dependent mechanisms, resulting in nuclear translocation of IRF3 and activation of IRF3/IFNβ-dependent antibacterial effector mechanisms. This TLR4/IRF3 pathway of pathogen discrimination was activated by ceramide and by P-fimbriated E. coli, which use ceramide-anchored glycosphingolipid receptors. Relevance of this pathway for human disease was supported by polymorphic IRF3 promoter sequences, differing between children with severe, symptomatic kidney infection and children who were asymptomatic bacterial carriers. IRF3 promoter activity was reduced by the disease-associated genotype, consistent with the pathology in Irf3−/− mice. Host susceptibility to common infections like UTI may thus be strongly influenced by single gene modifications affecting the innate immune response.

Ca-sensor proteins are generally implicated in insulin release through SNARE interactions. Here, secretagogin, whose expression in human pancreatic islets correlates with their insulin content and the incidence of type 2 diabetes, is shown to orchestrate an unexpectedly distinct mechanism. Single-cell RNA-seq reveals retained expression of the TRP family members in β-cells from diabetic donors. Amongst these, pharmacological probing identifies Ca-permeable transient receptor potential vanilloid type 1 channels (TRPV1) as potent inducers of secretagogin expression through recruitment of Sp1 transcription factors. Accordingly, agonist stimulation of TRPV1s fails to rescue insulin release from pancreatic islets of glucose intolerant secretagogin knock-out() mice. However, instead of merely impinging on the SNARE machinery, reduced insulin availability in secretagogin mice is due to β-cell loss, which is underpinned by the collapse of protein folding and deregulation of secretagogin-dependent USP9X deubiquitinase activity. Therefore, and considering the desensitization of TRPV1s in diabetic pancreata, a TRPV1-to-secretagogin regulatory axis seems critical to maintain the structural integrity and signal competence of β-cells.

Cartilage oligomeric matrix protein (COMP) is a soluble pentameric protein expressed in cartilage and involved in collagen organization. Tissue microarrays derived from two cohorts of patients with breast cancer (n=122 and n=498) were immunostained, revealing varying expression of COMP, both in the tumor cells and surrounding stroma. High levels of COMP in tumor cells correlated, independently of other variables, with poor survival and decreased recurrence-free survival. Breast cancer cells, MDA-MB-231, stably expressing COMP were injected into the mammary fat pad of SCID (CB-17/Icr-Prkdc/Rj) mice. Tumors expressing COMP were significantly larger and were more prone to metastasize as compared with control, mock-transfected, tumors. In vitro experiments confirmed that COMP-expressing cells had a more invasive phenotype, which could in part be attributed to an upregulation of matrix metalloprotease-9. Furthermore, microarray analyses of gene expression in tumors formed in vivo showed that COMP expression induced higher expression of genes protecting against endoplasmic reticulum stress. This observation was confirmed in vitro as COMP-expressing cells showed better survival as well as a higher rate of protein synthesis when treated with brefeldin A, compared with control cells. Further, COMP-expressing cells appeared to undergo a metabolic switch, that is, a Warburg effect. Thus, in vitro measurement of cell respiration indicated decreased mitochondrial metabolism. In conclusion, COMP is a novel biomarker in breast cancer, which contributes to the severity of the disease by metabolic switching and increasing invasiveness and tumor cell viability, leading to reduced survival in animal models and human patients.

No abstract