Pancreatic islet enhancer clusters enriched in type 2 diabetes risk-associated variants

Tissue and cancer-specific expression of DIEXF is epigenetically mediated by an Alu repeat

Alu repeats constitute a major fraction of human genome and for a small subset of them a role in gene regulation has been described. The number of studies focused on the functional characterization of particular Alu elements is very limited. Most Alu elements are DNA methylated and then assumed to lie in repressed chromatin domains. We hypothesize that Alu elements with low or variable DNA methylation are candidates for a functional role. In a genome-wide study in normal and cancer tissues, we pinpointed an Alu repeat (AluSq2) with differential methylation located upstream of the promoter region of the DIEXF gene. DIEXF encodes a highly conserved factor essential for the development of zebrafish digestive tract. To characterize the contribution of the Alu element to the regulation of DIEXF we analysed the epigenetic landscapes of the gene promoter and flanking regions in different cell types and cancers. Alternate epigenetic profiles (DNA methylation and histone modifications) of the AluSq2 element were associated with DIEXF transcript diversity as well as protein levels, while the epigenetic profile of the CpG island associated with the DIEXF promoter remained unchanged. These results suggest that AluSq2 might directly contribute to the regulation of DIEXF transcription and protein expression. Moreover, AluSq2 was DNA hypomethylated in different cancer types, pointing out its putative contribution to DIEXF alteration in cancer and its potential as tumoural biomarker.

Liver gene regulatory networks: Contributing factors to nonalcoholic fatty liver disease

Metabolic diseases such as nonalcoholic fatty liver disease (NAFLD) result from complex interactions between intrinsic and extrinsic factors, including genetics and exposure to obesogenic environments. These risk factors converge in aberrant gene expression patterns in the liver, which are underlined by altered cis-regulatory networks. In homeostasis and in disease states, liver cis-regulatory networks are established by coordinated action of liver-enriched transcription factors (TFs), which define enhancer landscapes, activating broad gene programs with spatiotemporal resolution. Recent advances in DNA sequencing have dramatically expanded our ability to map active transcripts, enhancers and TF cistromes, and to define the 3D chromatin topology that contains these elements. Deployment of these technologies has allowed investigation of the molecular processes that regulate liver development and metabolic homeostasis. Moreover, genomic studies of NAFLD patients and NAFLD models have demonstrated that the liver undergoes pervasive regulatory rewiring in NAFLD, which is reflected by aberrant gene expression profiles. We have therefore achieved an unprecedented level of detail in the understanding of liver cis-regulatory networks, particularly in physiological conditions. Future studies should aim to map active regulatory elements with added levels of resolution, addressing how the chromatin landscapes of different cell lineages contribute to and are altered in NAFLD and NAFLD-associated metabolic states. Such efforts would provide additional clues into the molecular factors that trigger this disease. This article is categorized under: Biological Mechanisms > Metabolism Biological Mechanisms > Regulatory Biology Laboratory Methods and Technologies > Genetic/Genomic Methods.

Single-cell ATAC-Seq in human pancreatic islets and deep learning upscaling of rare cells reveals cell-specific type 2 diabetes regulatory signatures

OBJECTIVE

Type 2 diabetes (T2D) is a complex disease characterized by pancreatic islet dysfunction, insulin resistance, and disruption of blood glucose levels. Genome-wide association studies (GWAS) have identified > 400 independent signals that encode genetic predisposition. More than 90% of associated single-nucleotide polymorphisms (SNPs) localize to non-coding regions and are enriched in chromatin-defined islet enhancer elements, indicating a strong transcriptional regulatory component to disease susceptibility. Pancreatic islets are a mixture of cell types that express distinct hormonal programs, so each cell type may contribute differentially to the underlying regulatory processes that modulate T2D-associated transcriptional circuits. Existing chromatin profiling methods such as ATAC-seq and DNase-seq, applied to islets in bulk, produce aggregate profiles that mask important cellular and regulatory heterogeneity.

METHODS

We present genome-wide single-cell chromatin accessibility profiles in >1,600 cells derived from a human pancreatic islet sample using single-cell combinatorial indexing ATAC-seq (sci-ATAC-seq). We also developed a deep learning model based on U-Net architecture to accurately predict open chromatin peak calls in rare cell populations.

RESULTS

We show that sci-ATAC-seq profiles allow us to deconvolve alpha, beta, and delta cell populations and identify cell-type-specific regulatory signatures underlying T2D. Particularly, T2D GWAS SNPs are significantly enriched in beta cell-specific and across cell-type shared islet open chromatin, but not in alpha or delta cell-specific open chromatin. We also demonstrate, using less abundant delta cells, that deep learning models can improve signal recovery and feature reconstruction of rarer cell populations. Finally, we use co-accessibility measures to nominate the cell-specific target genes at 104 non-coding T2D GWAS signals.

CONCLUSIONS

Collectively, we identify the islet cell type of action across genetic signals of T2D predisposition and provide higher-resolution mechanistic insights into genetically encoded risk pathways.

Deep learning models predict regulatory variants in pancreatic islets and refine type 2 diabetes association signals

Genome-wide association analyses have uncovered multiple genomic regions associated with T2D, but identification of the causal variants at these remains a challenge. There is growing interest in the potential of deep learning models - which predict epigenome features from DNA sequence - to support inference concerning the regulatory effects of disease-associated variants. Here, we evaluate the advantages of training convolutional neural network (CNN) models on a broad set of epigenomic features collected in a single disease-relevant tissue - pancreatic islets in the case of type 2 diabetes (T2D) - as opposed to models trained on multiple human tissues. We report convergence of CNN-based metrics of regulatory function with conventional approaches to variant prioritization - genetic fine-mapping and regulatory annotation enrichment. We demonstrate that CNN-based analyses can refine association signals at T2D-associated loci and provide experimental validation for one such signal. We anticipate that these approaches will become routine in downstream analyses of GWAS.

Computational and functional analyses of T2D GWAS SNPs for transcription factor binding

Genome-wide association studies (GWASs) have successfully identified numerous non-coding genetic variants for type 2 diabetes (T2D), but the functional roles underlying these non-coding variants remain largely unknown. The effects of T2D GWAS lead SNPs on transcriptional factors binding motifs were firstly analyzed via JASPAR, followed by functional validations including dual-luciferase reporter assays, biotin-based DNA pull-down assays, real-time quantitative PCR, and western blotting. The results showed that GWAS SNP rs4430796 conferred T allele specific transcriptional enhancer activity via a PAX6 binding element, and upregulated the expression of HNF1B. GWAS SNP rs4607103 showed a bidirectional modulation of ADAMTS9-AS2 and ADAMTS9 by TCF7L2 in a T allele-specific manner. GWAS SNP rs849135 conferred C allele-specific bidirectional transcriptional enhancer activity via a CREB1 binding element. Our findings have uncovered the functional mechanisms of three T2D GWAS SNPs via affecting the binding of transcription factors, providing new insights into the genetics and molecular pathogenesis of T2D.

Functionally Annotating Regulatory Elements in the Equine Genome Using Histone Mark ChIP-Seq

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein–protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse.

Circadian Entrainment Triggers Maturation of Human In Vitro Islets

Stem-cell-derived tissues could transform disease research and therapy, yet most methods generate functionally immature products. We investigate how human pluripotent stem cells (hPSCs) differentiate into pancreatic islets in vitro by profiling DNA methylation, chromatin accessibility, and histone modification changes. We find that enhancer potential is reset upon lineage commitment and show how pervasive epigenetic priming steers endocrine cell fates. Modeling islet differentiation and maturation regulatory circuits reveals genes critical for generating endocrine cells and identifies circadian control as limiting for in vitro islet function. Entrainment to circadian feeding/fasting cycles triggers islet metabolic maturation by inducing cyclic synthesis of energy metabolism and insulin secretion effectors, including antiphasic insulin and glucagon pulses. Following entrainment, hPSC-derived islets gain persistent chromatin changes and rhythmic insulin responses with a raised glucose threshold, a hallmark of functional maturity, and function within days of transplantation. Thus, hPSC-derived tissues are amenable to functional improvement by circadian modulation.

Recovery of viable endocrine-specific cells and transcriptomes from human pancreatic islet-engrafted mice

Human pancreatic islets engrafted into immunodeficient mice serve as an important model for in vivo human diabetes studies. Following engraftment, islet function can be monitored in vivo by measuring circulating glucose and human insulin; however, it will be important to recover viable cells for more complex graft analyses. Moreover, RNA analyses of dissected grafts have not distinguished which hormone-specific cell types contribute to gene expression. We developed a method for recovering live cells suitable for fluorescence-activated cell sorting from human islets engrafted in mice. Although yields of recovered islet cells were relatively low, the ratios of bulk-sorted β, α, and δ cells and their respective hormone-specific RNA-Seq transcriptomes are comparable pretransplant and posttransplant, suggesting that the cellular characteristics of islet grafts posttransplant closely mirror the original donor islets. Single-cell RNA-Seq transcriptome analysis confirms the presence of appropriate β, α, and δ cell subsets. In addition, ex vivo perifusion of recovered human islet grafts demonstrated glucose-stimulated insulin secretion. Viable cells suitable for patch-clamp analysis were recovered from transplanted human embryonic stem cell-derived β cells. Together, our functional and hormone-specific transcriptome analyses document the broad applicability of this system for longitudinal examination of human islet cells undergoing developmental/metabolic/pharmacogenetic manipulation in vivo and may facilitate the discovery of treatments for diabetes.

UHPLC/QTOF-MS-based metabolomics reveal the effect of Melastoma dodecandrum extract in type 2 diabetic rats

Context: Melastoma dodecandrum Lour. (Melastomataceae) is a traditional Chinese medicine. This is the first study to report a protective effect of the ethanol extract from M. dodecandrum (MDE) in type 2 diabetic (T2DM) rats.Objective: To investigate the therapeutic mechanism of MDE in T2DM rats.Materials and methods: Sprague-Dawley rats were fed a high-fat diet for 6 consecutive weeks, followed by intraperitoneal injection of streptozotocin (STZ) (30 mg/kg) to induce diabetes. T2DM rats were divided into untreated diabetic, metformin-treated and MDE-treated groups. Additionally, normal rats without treatment served as the control group (n = 6). Metformin (250 mg/kg) and MDE (600 mg/kg) were intragastrically administered to T2DM rats for 5 consecutive weeks. Serum samples were evaluated via ultra-performance liquid chromatography quadrupole time-of-flight mass spectrometry (UHPLC/QTOF-MS), followed by principal components analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA).Results: The 17 identified potential biomarkers were mainly involved in lipid, amino acid, arachidonic acid, taurine and nicotinic acid metabolism. MDE also significantly reduced the level of fasting blood glucose (FBG), oral glucose tolerance, insulin, total cholesterol (TC), triglyceride (TG), low-density lipoprotein (LDL), malondialdehyde (MDA), aspartate aminotransferase (AST), alanine aminotransferase (ALT) and urea nitrogen (BUN) in T2DM rats. The high-density lipoprotein (HDL), serum creatinine (Scr), superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GSH-Px) levels were elevated in MDE-treated group.Discussion and conclusion: MDE possesses substantial antidiabetic activity, especially in lipid disorder regulation. This suggests that the use of MDE can be generalized to broader pharmacological studies, such as obesity and hyperlipidaemia.

Claudin 4 in pancreatic β cells is involved in regulating the functional state of adult islets

The functional state (FS) of adult pancreatic islets is regulated by a large array of regulatory molecules including numerous transcription factors. Whether any islet structural molecules play such a role has not been well understood. Here, multiple technologies including bioinformatics analyses were used to explore such molecules. The tight junction family molecule claudin 4 (Cldn4) was the highest enriched amongst over 140 structural genes analysed. Cldn4 expression was ~75-fold higher in adult islets than in exocrine tissues and was mostly up-regulated during functional maturation of developing islet cells. Cldn4 was progressively down-regulated in functionally compromised, dedifferentiating insulin-secreting β cells and in db/db type 2 diabetic islets. Furthermore, the genetic deletion of Cldn4 impaired significantly the FS without apparently affecting pancreas morphology, islet architectural structure and cellular distribution, and secretion of enteroendocrine hormones. Thus, we suggest a previously unidentified role for Cldn4 in regulating the FS of islets, with implications in translational research for better diabetes therapies.

Pancreatic Islet Transcriptional Enhancers and Diabetes

PURPOSE OF REVIEW

Common genetic variants that associate with type 2 diabetes risk are markedly enriched in pancreatic islet transcriptional enhancers. This review discusses current advances in the annotation of islet enhancer variants and their target genes.

RECENT FINDINGS

Recent methodological advances now allow genetic and functional mapping of diabetes causal variants at unprecedented resolution. Mapping of enhancer-promoter interactions in human islets has provided a unique appreciation of the complexity of islet gene regulatory processes and enabled direct association of noncoding diabetes risk variants to their target genes. The recently improved human islet enhancer annotations constitute a framework for the interpretation of diabetes genetic signals in the context of pancreatic islet gene regulation. In the future, integration of existing and yet to come regulatory maps with genetic fine-mapping efforts and in-depth functional characterization will foster the discovery of novel diabetes molecular risk mechanisms.

The impact of proinflammatory cytokines on the β-cell regulatory landscape provides insights into the genetics of type 1 diabetes

Early stages of type 1 diabetes (T1D) are characterized by local autoimmune inflammation and progressive loss of insulin-producing pancreatic β cells. We show here that exposure to pro-inflammatory cytokines unmasks a marked plasticity of the β-cell regulatory landscape. We expand the repertoire of human islet regulatory elements by mapping stimulus-responsive enhancers linked to changes in the β-cell transcriptome, proteome and 3D chromatin structure. Our data indicate that the β cell response to cytokines is mediated by the induction of new regulatory regions as well as the activation of primed regulatory elements prebound by islet-specific transcription factors. We find that T1D-associated loci are enriched of the newly mapped cis-regulatory regions and identify T1D-associated variants disrupting cytokine-responsive enhancer activity in human β cells. Our study illustrates how β cells respond to a pro-inflammatory environment and implicate a role for stimulus-response islet enhancers in T1D.

Human stem cell models: lessons for pancreatic development and disease

A comprehensive understanding of mechanisms that underlie the development and function of human cells requires human cell models. For the pancreatic lineage, protocols have been developed to differentiate human pluripotent stem cells (hPSCs) into pancreatic endocrine and exocrine cells through intermediates resembling in vivo development. In recent years, this differentiation system has been employed to decipher mechanisms of pancreatic development, congenital defects of the pancreas, as well as genetic forms of diabetes and exocrine diseases. In this review, we summarize recent insights gained from studies of pancreatic hPSC models. We discuss how genome-scale analyses of the differentiation system have helped elucidate roles of chromatin state, transcription factors, and noncoding RNAs in pancreatic development and how the analysis of cells with disease-relevant mutations has provided insight into the molecular underpinnings of genetically determined diseases of the pancreas.

Identification of atrial fibrillation associated genes and functional non-coding variants

Disease-associated genetic variants that lie in non-coding regions found by genome-wide association studies are thought to alter the functionality of transcription regulatory elements and target gene expression. To uncover causal genetic variants, variant regulatory elements and their target genes, here we cross-reference human transcriptomic, epigenomic and chromatin conformation datasets. Of 104 genetic variant regions associated with atrial fibrillation candidate target genes are prioritized. We optimize EMERGE enhancer prediction and use accessible chromatin profiles of human atrial cardiomyocytes to more accurately predict cardiac regulatory elements and identify hundreds of sub-threshold variants that co-localize with regulatory elements. Removal of mouse homologues of atrial fibrillation-associated regions in vivo uncovers a distal regulatory region involved in Gja1 (Cx43) expression. Our analyses provide a shortlist of genes likely affected by atrial fibrillation-associated variants and provide variant regulatory elements in each region that link genetic variation and target gene regulation, helping to focus future investigations.

Beta-cell dysfunction induced by non-cytotoxic concentrations of Interleukin-1β is associated with changes in expression of beta-cell maturity genes and associated histone modifications

Decreased insulin secretory capacity in Type 2 diabetes mellitus is associated with beta-cell dedifferentiation and inflammation. We hypothesize that prolonged exposure of beta-cells to low concentrations of IL-1β induce beta-cell dedifferentiation characterized by impaired glucose-stimulated insulin secretion, reduced expression of key beta-cell genes and changes in histone modifications at gene loci known to affect beta-cell function. Ten days exposure to IL-1β at non-cytotoxic concentrations reduced insulin secretion and beta-cell proliferation and decreased expression of key beta-cell identity genes, including MafA and Ucn3 and decreased H3K27ac at the gene loci, suggesting that inflammatory cytokines directly affects the epigenome. Following removal of IL-1β, beta-cell function was normalized and mRNA expression of beta-cell identity genes, such as insulin and Ucn3 returned to pre-stimulation levels. Our findings indicate that prolonged exposure to low concentrations of IL-1β induces epigenetic changes associated with loss of beta-cell identity as observed in Type 2 diabetes.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

Identification of C2CD4A as a human diabetes susceptibility gene with a role in β cell insulin secretion

Fine mapping and validation of genes causing β cell failure from susceptibility loci identified in type 2 diabetes genome-wide association studies (GWAS) poses a significant challenge. The VPS13C-C2CD4A-C2CD4B locus on chromosome 15 confers diabetes susceptibility in every ethnic group studied to date. However, the causative gene is unknown. FoxO1 is involved in the pathogenesis of β cell dysfunction, but its link to human diabetes GWAS has not been explored. Here we generated a genome-wide map of FoxO1 superenhancers in chemically identified β cells using 2-photon live-cell imaging to monitor FoxO1 localization. When parsed against human superenhancers and GWAS-derived diabetes susceptibility alleles, this map revealed a conserved superenhancer in C2CD4A, a gene encoding a β cell/stomach-enriched nuclear protein of unknown function. Genetic ablation of C2cd4a in β cells of mice phenocopied the metabolic abnormalities of human carriers of C2CD4A-linked polymorphisms, resulting in impaired insulin secretion during glucose tolerance tests as well as hyperglycemic clamps. C2CD4A regulates glycolytic genes, and notably represses key β cell "disallowed" genes, such as lactate dehydrogenase A We propose that C2CD4A is a transcriptional coregulator of the glycolytic pathway whose dysfunction accounts for the diabetes susceptibility associated with the chromosome 15 GWAS locus.

DiseaseEnhancer: a resource of human disease-associated enhancer catalog

Large-scale sequencing studies discovered substantial genetic variants occurring in enhancers which regulate genes via long range chromatin interactions. Importantly, such variants could affect enhancer regulation by changing transcription factor bindings or enhancer hijacking, and in turn, make an essential contribution to disease progression. To facilitate better usage of published data and exploring enhancer deregulation in various human diseases, we created DiseaseEnhancer (http://biocc.hrbmu.edu.cn/DiseaseEnhancer/), a manually curated database for disease-associated enhancers. As of July 2017, DiseaseEnhancer includes 847 disease-associated enhancers in 143 human diseases. Database features include basic enhancer information (i.e. genomic location and target genes); disease types; associated variants on the enhancer and their mediated phenotypes (i.e. gain/loss of enhancer and the alterations of transcription factor bindings). We also include a feature on our website to export any query results into a file and download the full database. DiseaseEnhancer provides a promising avenue for researchers to facilitate the understanding of enhancer deregulation in disease pathogenesis, and identify new biomarkers for disease diagnosis and therapy.

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

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

The long noncoding RNA MALAT1 predicts human pancreatic islet isolation quality

Human islet isolation is a cost-/resource-intensive program generating islets for cell therapy in Type 1 diabetes. However, only a third of cadaveric pancreas get to clinical transplantation due to low quality/number of islets. There is a need to identify biomarker(s) that predict the quality of islets, prior to initiating their isolation. Here, we sequenced transcriptome from 18 human islet preparations stratified into three groups (Gr.1: Best quality/transplantable islets, Gr.2: Intermediary quality, Gr.3: Inferior quality/non-transplantable islets) based on routine measurements including islet purity/viability. Machine-learning algorithms involving penalized regression analyses identified 10 long-non-coding(lnc)RNAs significantly different across all group-wise comparisons (Gr1VsGr2, Gr2vsGr3, Gr1vsGr3). Two variants of Metastasis-Associated Lung Adenocarcinoma Transcript-1(MALAT1) lncRNA were common across all comparisons. We confirmed RNA-seq findings in a "validation set" of 75 human islet preparations. Finally, in 19 pancreas samples, we demonstrate that assessing the levels of MALAT1 variants alone (ROC curve AUC: 0.83) offers highest specificity in predicting post-isolation islet quality and improves the predictive potential for clinical islet transplantation when combined with Edmonton Donor Points/Body Mass Index(BMI)/North American Islet Donor Score(NAIDS). We present this resource of islet-quality-stratified lncRNA transcriptome data and identify MALAT1 as a biomarker that significantly enhances current selection methods for clinical (GMP)-grade islet isolation.

Gene loci associated with insulin secretion in islets from non-diabetic mice

Genetic susceptibility to type 2 diabetes is primarily due to β-cell dysfunction. However, a genetic study to directly interrogate β-cell function ex vivo has never been previously performed. We isolated 233,447 islets from 483 Diversity Outbred (DO) mice maintained on a Western-style diet, and measured insulin secretion in response to a variety of secretagogues. Insulin secretion from DO islets ranged >1,000-fold even though none of the mice were diabetic. The insulin secretory response to each secretagogue had a unique genetic architecture; some of the loci were specific for one condition, whereas others overlapped. Human loci that are syntenic to many of the insulin secretion QTL from mouse are associated with diabetes-related SNPs in human genome-wide association studies. We report on three genes, Ptpn18, Hunk and Zfp148, where the phenotype predictions from the genetic screen were fulfilled in our studies of transgenic mouse models. These three genes encode a non-receptor type protein tyrosine phosphatase, a serine/threonine protein kinase, and a Krϋppel-type zinc-finger transcription factor, respectively. Our results demonstrate that genetic variation in insulin secretion that can lead to type 2 diabetes is discoverable in non-diabetic individuals.

The Genetic Epidemiology of Type 2 Diabetes: Opportunities for Health Translation

PURPOSE OF REVIEW

Genome-wide association studies have delineated the genetic architecture of type 2 diabetes. While functional studies to identify target transcripts are ongoing, new genetic knowledge can be translated directly to health applications. The review covers several translation directions but focuses on genomic polygenic scores for screening and prevention.

RECENT FINDINGS

Over 400 genomic variants associated with T2D and its related quantitative traits are now known. Genetic scores comprising dozens to millions of associated variants can predict incident T2D. However, measurement of body mass index is more efficient than genetic scores to detect T2D risk groups, and knowledge of T2D genetic risk alone seems insufficient to improve health. Genetically determined metabolic sub-phenotypes can be identified by clustering variants associated with physiological axes like insulin resistance. Genetic sub-phenotyping may be a way forward to identify specific individual phenotypes for prevention and treatment. Genomic polygenic scores for T2D can predict incident diabetes but may not be useful to improve health overall. Genetic detection of T2D sub-phenotypes could be useful to personalize screening and care.

GLIS3 binds pancreatic beta cell regulatory regions alongside other islet transcription factors

The Krüppel-like zinc finger transcription factor Gli-similar 3 (GLIS3) plays a critical role in the regulation of pancreatic beta cells, with global Glis3 knockout mice suffering from severe hyperglycemia and dying by post-natal day 11. In addition, GLIS3 has been shown to directly regulate the early endocrine marker Ngn3, as well as Ins2 gene expression in mature beta cells. We hypothesize that GLIS3 regulates several other genes critical to beta cell function, in addition to Ins2, by directly binding to regulatory regions. We therefore generated a pancreas-specific Glis3 deletion mouse model (Glis3Δpanc) using a Pdx1-driven Cre mouse line. Roughly 20% of these mice develop hyperglycemia by 8-weeks and lose most of their insulin expression. However, this did not appear to be due to loss of the beta cells themselves, as no change in cell death was observed. Indeed, presumptive beta cells appeared to persist as PDX1+/INS-/MAFA-/GLUT2- cells. Islet RNA-seq analysis combined with GLIS3 ChIP-seq analysis revealed apparent direct regulation of a variety of diabetes related genes, including Slc2a2 and Mafa. GLIS3 binding near these genes coincided with binding for other islet-enriched transcription factors, indicating these are distinct regulatory hubs. Our data indicates that GLIS3 not only regulates insulin expression, but several additional genes critical for beta cell function.

Enhancer signatures stratify and predict outcomes of non-functional pancreatic neuroendocrine tumors

Most pancreatic neuroendocrine tumors (PNETs) do not produce excess hormones and are therefore considered 'non-functional'^1-3. As clinical behaviors vary widely and distant metastases are eventually lethal^2,4, biological classifications might guide treatment. Using enhancer maps to infer gene regulatory programs, we find that non-functional PNETs fall into two major subtypes, with epigenomes and transcriptomes that partially resemble islet α- and β-cells. Transcription factors ARX and PDX1 specify these normal cells, respectively^5,6, and 84% of 142 non-functional PNETs expressed one or the other factor, occasionally both. Among 103 cases, distant relapses occurred almost exclusively in patients with ARX⁺PDX1^- tumors and, within this subtype, in cases with alternative lengthening of telomeres. These markedly different outcomes belied similar clinical presentations and histology and, in one cohort, occurred irrespective of MEN1 mutation. This robust molecular stratification provides insight into cell lineage correlates of non-functional PNETs, accurately predicts disease course and can inform postoperative clinical decisions.

Human pancreatic islet three-dimensional chromatin architecture provides insights into the genetics of type 2 diabetes

Genetic studies promise to provide insight into the molecular mechanisms underlying type 2 diabetes (T2D). Variants associated with T2D are often located in tissue-specific enhancer clusters or super-enhancers. So far, such domains have been defined through clustering of enhancers in linear genome maps rather than in three-dimensional (3D) space. Furthermore, their target genes are often unknown. We have created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated enhancers to their target genes, often located hundreds of kilobases away. It also revealed >1,300 groups of islet enhancers, super-enhancers and active promoters that form 3D hubs, some of which show coordinated glucose-dependent activity. We demonstrate that genetic variation in hubs impacts insulin secretion heritability, and show that hub annotations can be used for polygenic scores that predict T2D risk driven by islet regulatory variants. Human islet 3D chromatin architecture, therefore, provides a framework for interpretation of T2D genome-wide association study (GWAS) signals.

Genetic associations of breast and prostate cancer are enriched for regulatory elements identified in disease-related tissues

Although genome-wide association studies (GWAS) have identified hundreds of risk loci for breast and prostate cancer, only a few studies have characterized the GWAS association signals across functional genomic annotations with a particular focus on single nucleotide polymorphisms (SNPs) located in DNA regulatory elements. In this study, we investigated the enrichment pattern of GWAS signals for breast and prostate cancer in genomic functional regions located in normal tissue and cancer cell lines. We quantified the overall enrichment of SNPs with breast and prostate cancer association p values < 1 × 10^-8 across regulatory categories. We then obtained annotations for DNaseI hypersensitive sites (DHS), typical enhancers, and super enhancers across multiple tissue types, to assess if significant GWAS signals were selectively enriched in annotations found in disease-related tissue. Finally, we quantified the enrichment of breast and prostate cancer SNP heritability in regulatory regions, and compared the enrichment pattern of SNP heritability with GWAS signals. DHS, typical enhancers, and super enhancers identified in the breast cancer cell line MCF-7 were observed with the highest enrichment of genome-wide significant variants for breast cancer. For prostate cancer, GWAS signals were mostly enriched in DHS and typical enhancers identified in the prostate cancer cell line LNCaP. With progressively stringent GWAS p value thresholds, an increasing trend of enrichment was observed for both diseases in DHS, typical enhancers, and super enhancers located in disease-related tissue. Results from heritability enrichment analysis supported the selective enrichment pattern of functional genomic regions in disease-related cell lines for both breast and prostate cancer. Our results suggest the importance of studying functional annotations identified in disease-related tissues when characterizing GWAS results, and further demonstrate the role of germline DNA regulatory elements from disease-related tissue in breast and prostate carcinogenesis.

Genome-wide enhancer annotations differ significantly in genomic distribution, evolution, and function

Background

Non-coding gene regulatory enhancers are essential to transcription in mammalian cells. As a result, a large variety of experimental and computational strategies have been developed to identify cis-regulatory enhancer sequences. Given the differences in the biological signals assayed, some variation in the enhancers identified by different methods is expected; however, the concordance of enhancers identified by different methods has not been comprehensively evaluated. This is critically needed, since in practice, most studies consider enhancers identified by only a single method. Here, we compare enhancer sets from eleven representative strategies in four biological contexts.

Results

All sets we evaluated overlap significantly more than expected by chance; however, there is significant dissimilarity in their genomic, evolutionary, and functional characteristics, both at the element and base-pair level, within each context. The disagreement is sufficient to influence interpretation of candidate SNPs from GWAS studies, and to lead to disparate conclusions about enhancer and disease mechanisms. Most regions identified as enhancers are supported by only one method, and we find limited evidence that regions identified by multiple methods are better candidates than those identified by a single method. As a result, we cannot recommend the use of any single enhancer identification strategy in all settings.

Conclusions

Our results highlight the inherent complexity of enhancer biology and identify an important challenge to mapping the genetic architecture of complex disease. Greater appreciation of how the diverse enhancer identification strategies in use today relate to the dynamic activity of gene regulatory regions is needed to enable robust and reproducible results.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5779-x) contains supplementary material, which is available to authorized users.

ATAC-seq reveals alterations in open chromatin in pancreatic islets from subjects with type 2 diabetes

Impaired insulin secretion from pancreatic islets is a hallmark of type 2 diabetes (T2D). Altered chromatin structure may contribute to the disease. We therefore studied the impact of T2D on open chromatin in human pancreatic islets. We used assay for transposase-accessible chromatin using sequencing (ATAC-seq) to profile open chromatin in islets from T2D and non-diabetic donors. We identified 57,105 and 53,284 ATAC-seq peaks representing open chromatin regions in islets of non-diabetic and diabetic donors, respectively. The majority of ATAC-seq peaks mapped near transcription start sites. Additionally, peaks were enriched in enhancer regions and in regions where islet-specific transcription factors (TFs), e.g. FOXA2, MAFB, NKX2.2, NKX6.1 and PDX1, bind. Islet ATAC-seq peaks overlap with 13 SNPs associated with T2D (e.g. rs7903146, rs2237897, rs757209, rs11708067 and rs878521 near TCF7L2, KCNQ1, HNF1B, ADCY5 and GCK, respectively) and with additional 67 SNPs in LD with known T2D SNPs (e.g. SNPs annotated to GIPR, KCNJ11, GLIS3, IGF2BP2, FTO and PPARG). There was enrichment of open chromatin regions near highly expressed genes in human islets. Moreover, 1,078 open chromatin peaks, annotated to 898 genes, differed in prevalence between diabetic and non-diabetic islet donors. Some of these peaks are annotated to candidate genes for T2D and islet dysfunction (e.g. HHEX, HMGA2, GLIS3, MTNR1B and PARK2) and some overlap with SNPs associated with T2D (e.g. rs3821943 near WFS1 and rs508419 near ANK1). Enhancer regions and motifs specific to key TFs including BACH2, FOXO1, FOXA2, NEUROD1, MAFA and PDX1 were enriched in differential islet ATAC-seq peaks of T2D versus non-diabetic donors. Our study provides new understanding into how T2D alters the chromatin landscape, and thereby accessibility for TFs and gene expression, in human pancreatic islets.

Induction of α cell-restricted Gc in dedifferentiating β cells contributes to stress-induced β-cell dysfunction

Diabetic β cell failure is associated with β cell dedifferentiation. To identify effector genes of dedifferentiation, we integrated analyses of histone methylation as a surrogate of gene activation status and RNA expression in β cells sorted from mice with multiparity-induced diabetes. Interestingly, only a narrow subset of genes demonstrated concordant changes to histone methylation and RNA levels in dedifferentiating β cells. Notable among them was the α cell signature gene Gc, encoding a vitamin D-binding protein. While diabetes was associated with Gc induction, Gc-deficient islets did not induce β cell dedifferentiation markers and maintained normal ex vivo insulin secretion in the face of metabolic challenge. Moreover, Gc-deficient mice exhibited a more robust insulin secretory response than normal controls during hyperglycemic clamps. The data are consistent with a functional role of Gc activation in β cell dysfunction, and indicate that multiparity-induced diabetes is associated with altered β cell fate.

Mechanistic Characterization of RASGRP1 Variants Identifies an hnRNP-K-Regulated Transcriptional Enhancer Contributing to SLE Susceptibility

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component. We recently identified a novel SLE susceptibility locus near RASGRP1, which governs the ERK/MAPK kinase cascade and B-/T-cell differentiation and development. However, precise causal RASGRP1 functional variant(s) and their mechanisms of action in SLE pathogenesis remain undefined. Our goal was to fine-map this locus, prioritize genetic variants likely to be functional, experimentally validate their biochemical mechanisms, and determine the contribution of these SNPs to SLE risk. We performed a meta-analysis across six Asian and European cohorts (9,529 cases; 22,462 controls), followed by in silico bioinformatic and epigenetic analyses to prioritize potentially functional SNPs. We experimentally validated the functional significance and mechanism of action of three SNPs in cultured T-cells. Meta-analysis identified 18 genome-wide significant (p < 5 × 10⁻⁸) SNPs, mostly concentrated in two haplotype blocks, one intronic and the other intergenic. Epigenetic fine-mapping, allelic, eQTL, and imbalance analyses predicted three transcriptional regulatory regions with four SNPs (rs7170151, rs11631591-rs7173565, and rs9920715) prioritized for functional validation. Luciferase reporter assays indicated significant allele-specific enhancer activity for intronic rs7170151 and rs11631591-rs7173565 in T-lymphoid (Jurkat) cells, but not in HEK293 cells. Following up with EMSA, mass spectrometry, and ChIP-qPCR, we detected allele-dependent interactions between heterogeneous nuclear ribonucleoprotein K (hnRNP-K) and rs11631591. Furthermore, inhibition of hnRNP-K in Jurkat and primary T-cells downregulated RASGRP1 and ERK/MAPK signaling. Comprehensive association, bioinformatics, and epigenetic analyses yielded putative functional variants of RASGRP1, which were experimentally validated. Notably, intronic variant (rs11631591) is located in a cell type-specific enhancer sequence, where its risk allele binds to the hnRNP-K protein and modulates RASGRP1 expression in Jurkat and primary T-cells. As risk allele dosage of rs11631591 correlates with increased RASGRP1 expression and ERK activity, we suggest that this SNP may underlie SLE risk at this locus.

Pancreatic islet chromatin accessibility and conformation reveals distal enhancer networks of type 2 diabetes risk

Genetic variants affecting pancreatic islet enhancers are central to T2D risk, but the gene targets of islet enhancer activity are largely unknown. We generate a high-resolution map of islet chromatin loops using Hi-C assays in three islet samples and use loops to annotate target genes of islet enhancers defined using ATAC-seq and published ChIP-seq data. We identify candidate target genes for thousands of islet enhancers, and find that enhancer looping is correlated with islet-specific gene expression. We fine-map T2D risk variants affecting islet enhancers, and find that candidate target genes of these variants defined using chromatin looping and eQTL mapping are enriched in protein transport and secretion pathways. At IGF2BP2, a fine-mapped T2D variant reduces islet enhancer activity and IGF2BP2 expression, and conditional inactivation of IGF2BP2 in mouse islets impairs glucose-stimulated insulin secretion. Our findings provide a resource for studying islet enhancer function and identifying genes involved in T2D risk.

Risk loci for type 2 diabetes (T2D) reside in pancreatic islet enhancers. Here, the authors generate high-resolution maps of islet chromatin conformation using Hi-C which they combine with ATAC-seq and ChIP-seq data to annotate candidate target genes of enhancers and validate IGF2BP2 activity in mouse islets.

Epigenetics in Human Obesity and Type 2 Diabetes

Epigenetic mechanisms control gene activity and the development of an organism. The epigenome includes DNA methylation, histone modifications, and RNA-mediated processes, and disruption of this balance may cause several pathologies and contribute to obesity and type 2 diabetes (T2D). This Review summarizes epigenetic signatures obtained from human tissues of relevance for metabolism-i.e., adipose tissue, skeletal muscle, pancreatic islets, liver, and blood-in relation to obesity and T2D. Although this research field is still young, these comprehensive data support not only a role for epigenetics in disease development, but also epigenetic alterations as a response to disease. Genetic predisposition, as well as aging, contribute to epigenetic variability, and several environmental factors, including exercise and diet, further interact with the human epigenome. The reversible nature of epigenetic modifications holds promise for future therapeutic strategies in obesity and T2D.

A transcriptome-wide association study of high-grade serous epithelial ovarian cancer identifies new susceptibility genes and splice variants

We sought to identify susceptibility genes for high-grade serous ovarian cancer (HGSOC) by performing a transcriptome-wide association study of gene expression and splice junction usage in HGSOC-relevant tissue types (N = 2,169) and the largest genome-wide association study available for HGSOC (N = 13,037 cases and 40,941 controls). We identified 25 transcriptome-wide association study significant genes, 7 at the junction level only, including LRRC46 at 19q21.32, (P = 1 × 10-9), CHMP4C at 8q21 (P = 2 × 10-11) and a PRC1 junction at 15q26 (P = 7 × 10-9). In vitro assays for CHMP4C showed that the associated variant induces allele-specific exon inclusion (P = 0.0024). Functional screens in HGSOC cell lines found evidence of essentiality for three of the new genes we identified: HAUS6, KANSL1 and PRC1, with the latter comparable to MYC. Our study implicates at least one target gene for 6 out of 13 distinct genome-wide association study regions, identifying 23 new candidate susceptibility genes for HGSOC.

Exploring the underlying biology of intrinsic cardiorespiratory fitness through integrative analysis of genomic variants and muscle gene expression profiling

Intrinsic cardiorespiratory fitness (CRF) is defined as the level of CRF in the sedentary state. There are large individual differences in intrinsic CRF among sedentary adults. The physiology of variability in CRF has received much attention, but little is known about the genetic and molecular mechanisms that impact intrinsic CRF. These issues were explored in the present study by interrogating intrinsic CRF-associated DNA sequence variation and skeletal muscle gene expression data from the HERITAGE Family Study through an integrative bioinformatics guided approach. A combined analytic strategy involving genetic association, pathway enrichment, tissue-specific network structure, cis-regulatory genome effects, and expression quantitative trait loci was used to select and rank genes through a variation-adjusted weighted ranking scheme. Prioritized genes were further interrogated for corroborative evidence from knockout mouse phenotypes and relevant physiological traits from the HERITAGE cohort. The mean intrinsic V̇o2max was 33.1 ml O2·kg-1·min-1 (SD = 8.8) for the sample of 493 sedentary adults. Suggestive evidence was found for gene loci related to cardiovascular physiology (ATE1, CASQ2, NOTO, and SGCG), hematopoiesis (PICALM, SSB, CA9, and CASQ2), skeletal muscle phenotypes (SGCG, DMRT2, ADARB1, and CASQ2), and metabolism (ATE1, PICALM, RAB11FIP5, GBA2, SGCG, PRADC1, ARL6IP5, and CASQ2). Supportive evidence for a role of several of these loci was uncovered via association between DNA variants and muscle gene expression levels with exercise cardiovascular and muscle physiological traits. This initial effort to define the underlying molecular substrates of intrinsic CRF warrants further studies based on appropriate cohorts and study designs, complemented by functional investigations. NEW & NOTEWORTHY Intrinsic cardiorespiratory fitness (CRF) is measured in the sedentary state and is highly variable among sedentary adults. The physiology of variability in intrinsic cardiorespiratory fitness has received much attention, but little is known about the genetic and molecular mechanisms that impact intrinsic CRF. These issues were explored computationally in the present study, with further corroborative evidence obtained from analysis of phenotype data from knockout mouse models and human cardiovascular and skeletal muscle measurements.

Genomic annotation of disease-associated variants reveals shared functional contexts

Variation in non-coding DNA, encompassing gene regulatory regions such as enhancers and promoters, contributes to risk for complex disorders, including type 2 diabetes. While genome-wide association studies have successfully identified hundreds of type 2 diabetes loci throughout the genome, the vast majority of these reside in non-coding DNA, which complicates the process of determining their functional significance and level of priority for further study. Here we review the methods used to experimentally annotate these non-coding variants, to nominate causal variants and to link them to diabetes pathophysiology. In recent years, chromatin profiling, massively parallel sequencing, high-throughput reporter assays and CRISPR gene editing technologies have rapidly become indispensable tools. Rather than treating individual variants in isolation, we discuss the importance of accounting for context, both genetic (such as flanking DNA sequence) and environmental (such as cellular state or environmental exposure). Incorporating these features shows promise in terms of revealing biologically convergent molecular signatures across distant and seemingly unrelated loci. Studying regulatory elements in the proper context will be crucial for interpreting the functional significance of disease-associated variants and applying the resulting knowledge to improve patient care.

The Transcription Factor ERG Regulates Super-Enhancers Associated With an Endothelial-Specific Gene Expression Program

RATIONALE

The ETS (E-26 transformation-specific) transcription factor ERG (ETS-related gene) is essential for endothelial homeostasis, driving expression of lineage genes and repressing proinflammatory genes. Loss of ERG expression is associated with diseases including atherosclerosis. ERG's homeostatic function is lineage-specific, because aberrant ERG expression in cancer is oncogenic. The molecular basis for ERG lineage-specific activity is unknown. Transcriptional regulation of lineage specificity is linked to enhancer clusters (super-enhancers).

OBJECTIVE

To investigate whether ERG regulates endothelial-specific gene expression via super-enhancers.

METHODS AND RESULTS

Chromatin immunoprecipitation with high-throughput sequencing in human umbilical vein endothelial cells showed that ERG binds 93% of super-enhancers ranked according to H3K27ac, a mark of active chromatin. These were associated with endothelial genes such as DLL4 (Delta-like protein 4), CLDN5 (claudin-5), VWF (von Willebrand factor), and CDH5 (VE-cadherin). Comparison between human umbilical vein endothelial cell and prostate cancer TMPRSS2 (transmembrane protease, serine-2):ERG fusion-positive human prostate epithelial cancer cell line (VCaP) cells revealed distinctive lineage-specific transcriptome and super-enhancer profiles. At a subset of endothelial super-enhancers (including DLL4 and CLDN5), loss of ERG results in significant reduction in gene expression which correlates with decreased enrichment of H3K27ac and MED (Mediator complex subunit)-1, and reduced recruitment of acetyltransferase p300. At these super-enhancers, co-occupancy of GATA2 (GATA-binding protein 2) and AP-1 (activator protein 1) is significantly lower compared with super-enhancers that remained constant following ERG inhibition. These data suggest distinct mechanisms of super-enhancer regulation in endothelial cells and highlight the unique role of ERG in controlling a core subset of super-enhancers. Most disease-associated single nucleotide polymorphisms from genome-wide association studies lie within noncoding regions and perturb transcription factor recognition sequences in relevant cell types. Analysis of genome-wide association studies data shows significant enrichment of risk variants for cardiovascular disease and other diseases, at ERG endothelial enhancers and super-enhancers.

CONCLUSIONS

The transcription factor ERG promotes endothelial homeostasis via regulation of lineage-specific enhancers and super-enhancers. Enrichment of cardiovascular disease-associated single nucleotide polymorphisms at ERG super-enhancers suggests that ERG-dependent transcription modulates disease risk.

Functional characteristics of novel pancreatic Pax6 regulatory elements

Complex diseases, such as diabetes, are influenced by comprehensive transcriptional networks. Genome-wide association studies have revealed that variants located in regulatory elements for pancreatic transcription factors are linked to diabetes, including those functionally linked to the paired box transcription factor Pax6. Pax6 deletions in adult mice cause rapid onset of classic diabetes, but the full spectrum of pancreatic Pax6 regulators is unknown. Using a regulatory element discovery approach, we identified two novel Pax6 pancreatic cis-regulatory elements in a poorly characterized regulatory desert. Both new elements, Pax6 pancreas cis-regulatory element 3 (PE3) and PE4, are located 50 and 100 kb upstream and interact with different parts of the Pax6 promoter and nearby non-coding RNAs. They drive expression in the developing pancreas and brain and code for multiple pancreas-related transcription factor-binding sites. PE3 binds CCCTC-binding factor (CTCF) and is marked by stem cell identity markers in embryonic stem cells, whilst a common variant located in the PE4 element affects binding of Pax4, a known pancreatic regulator, altering Pax6 gene expression. To determine the ability of these elements to regulate gene expression, synthetic transcriptional activators and repressors were targeted to PE3 and PE4, modulating Pax6 gene expression, as well as influencing neighbouring genes and long non-coding RNAs, implicating the Pax6 locus in pancreas function and diabetes.

Laser capture microdissection of human pancreatic islets reveals novel eQTLs associated with type 2 diabetes

OBJECTIVE

Genome wide association studies (GWAS) for type 2 diabetes (T2D) have identified genetic loci that often localise in non-coding regions of the genome, suggesting gene regulation effects. We combined genetic and transcriptomic analysis from human islets obtained from brain-dead organ donors or surgical patients to detect expression quantitative trait loci (eQTLs) and shed light into the regulatory mechanisms of these genes.

METHODS

Pancreatic islets were isolated either by laser capture microdissection (LCM) from surgical specimens of 103 metabolically phenotyped pancreatectomized patients (PPP) or by collagenase digestion of pancreas from 100 brain-dead organ donors (OD). Genotyping (> 8.7 million single nucleotide polymorphisms) and expression (> 47,000 transcripts and splice variants) analyses were combined to generate cis-eQTLs.

RESULTS

After applying genome-wide false discovery rate significance thresholds, we identified 1,173 and 1,021 eQTLs in samples of OD and PPP, respectively. Among the strongest eQTLs shared between OD and PPP were CHURC1 (OD p-value=1.71 × 10^-24; PPP p-value = 3.64 × 10^-24) and PSPH (OD p-value = 3.92 × 10^-26; PPP p-value = 3.64 × 10^-24). We identified eQTLs in linkage-disequilibrium with GWAS loci T2D and associated traits, including TTLL6, MLX and KIF9 loci, which do not implicate the nearest gene. We found in the PPP datasets 11 eQTL genes, which were differentially expressed in T2D and two genes (CYP4V2 and TSEN2) associated with HbA1c but none in the OD samples.

CONCLUSIONS

eQTL analysis of LCM islets from PPP led us to identify novel genes which had not been previously linked to islet biology and T2D. The understanding gained from eQTL approaches, especially using surgical samples of living patients, provides a more accurate 3-dimensional representation than those from genetic studies alone.

Combined Inhibition of DYRK1A, SMAD, and Trithorax Pathways Synergizes to Induce Robust Replication in Adult Human Beta Cells

Small-molecule inhibitors of dual-specificity tyrosine-regulated kinase 1A (DYRK1A) induce human beta cells to proliferate, generating a labeling index of 1.5%-3%. Here, we demonstrate that combined pharmacologic inhibition of DYRK1A and transforming growth factor beta superfamily (TGFβSF)/SMAD signaling generates remarkable further synergistic increases in human beta cell proliferation (average labeling index, 5%-8%, and as high as 15%-18%), and increases in both mouse and human beta cell numbers. This synergy reflects activation of cyclins and cdks by DYRK1A inhibition, accompanied by simultaneous reductions in key cell-cycle inhibitors (CDKN1C and CDKN1A). The latter results from interference with the basal Trithorax- and SMAD-mediated transactivation of CDKN1C and CDKN1A. Notably, combined DYRK1A and TGFβ inhibition allows preservation of beta cell differentiated function. These beneficial effects extend from normal human beta cells and stem cell-derived human beta cells to those from people with type 2 diabetes, and occur both in vitro and in vivo.

Integrative modeling reveals key chromatin and sequence signatures predicting super-enhancers

Super-enhancers (SEs) are clusters of transcriptional enhancers which control the expression of cell identity and disease-associated genes. Current studies demonstrated the role of multiple factors in SE formation; however, a systematic analysis to assess the relative predictive importance of chromatin and sequence features of SEs and their constituents is lacking. In addition, a predictive model that integrates various types of data to predict SEs has not been established. Here, we integrated diverse types of genomic and epigenomic datasets to identify key signatures of SEs and investigated their predictive importance. Through integrative modeling, we found Cdk8, Cdk9, and Smad3 as new features of SEs, which can define known and new SEs in mouse embryonic stem cells and pro-B cells. We compared six state-of-the-art machine learning models to predict SEs and showed that non-parametric ensemble models performed better as compared to parametric. We validated these models using cross-validation and also independent datasets in four human cell-types. Taken together, our systematic analysis and ranking of features can be used as a platform to define and understand the biology of SEs in other cell-types.

The Continuing Evolution of Precision Health in Type 2 Diabetes: Achievements and Challenges

PURPOSE OF REVIEW

The purpose of this review was to summarize recent advances in the genomics of type 2 diabetes (T2D) and to highlight current initiatives to advance precision health.

RECENT FINDINGS

Generation of multi-omic data to measure each of the "biologic layers," developments in describing genomic function and annotation in T2D relevant tissue, along with the increasing recognition that T2D is a heterogeneous disease, and large-scale collaborations have all contributed to advancing our understanding of the molecular basis of T2D. Substantial advances have been made in understanding the molecular basis of T2D pathogenesis, such that precision health diabetes is increasingly becoming a reality. For precision diabetes to become a routine in clinical and public health, additional large-scale multi-omic initiatives are needed along with better assessment of our environment to delineate an individual's diabetes subtype for improved detection and management.

Long non-coding RNA repertoire and open chromatin regions constitute midbrain dopaminergic neuron - specific molecular signatures

Midbrain dopaminergic (DA) neurons are involved in diverse neurological functions, including control of movements, emotions or reward. In turn, their dysfunctions cause severe clinical manifestations in humans, such as the appearance of motor and cognitive symptoms in Parkinson’s Disease. The physiology and pathophysiology of these neurons are widely studied, mostly with respect to molecular mechanisms implicating protein-coding genes. In contrast, the contribution of non-coding elements of the genome to DA neuron function is poorly investigated. In this study, we isolated DA neurons from E14.5 ventral mesencephalons in mice, and used RNA-seq and ATAC-seq to establish and describe repertoires of long non-coding RNAs (lncRNAs) and putative DNA regulatory regions specific to this neuronal population. We identified 1,294 lncRNAs constituting the repertoire of DA neurons, among which 939 were novel. Most of them were not found in hindbrain serotonergic (5-HT) neurons, indicating a high degree of cell-specificity. This feature was also observed regarding open chromatin regions, as 39% of the ATAC-seq peaks from the DA repertoire were not detected in the 5-HT neurons. Our work provides for the first time DA-specific catalogues of non-coding elements of the genome that will undoubtedly participate in deepening our knowledge regarding DA neuronal development and dysfunctions.

Cell Specificity of Human Regulatory Annotations and Their Genetic Effects on Gene Expression

Epigenomic signatures from histone marks and transcription factor (TF)-binding sites have been used to annotate putative gene regulatory regions. However, a direct comparison of these diverse annotations is missing, and it is unclear how genetic variation within these annotations affects gene expression. Here, we compare five widely used annotations of active regulatory elements that represent high densities of one or more relevant epigenomic marks-"super" and "typical" (nonsuper) enhancers, stretch enhancers, high-occupancy target (HOT) regions, and broad domains-across the four matched human cell types for which they are available. We observe that stretch and super enhancers cover cell type-specific enhancer "chromatin states," whereas HOT regions and broad domains comprise more ubiquitous promoter states. Expression quantitative trait loci (eQTL) in stretch enhancers have significantly smaller effect sizes compared to those in HOT regions. Strikingly, chromatin accessibility QTL in stretch enhancers have significantly larger effect sizes compared to those in HOT regions. These observations suggest that stretch enhancers could harbor genetically primed chromatin to enable changes in TF binding, possibly to drive cell type-specific responses to environmental stimuli. Our results suggest that current eQTL studies are relatively underpowered or could lack the appropriate environmental context to detect genetic effects in the most cell type-specific "regulatory annotations," which likely contributes to infrequent colocalization of eQTL with genome-wide association study signals.

Prediction of drug response and adverse drug reactions: From twin studies to Next Generation Sequencing

Understanding and predicting inter-individual differences related to the success of drug therapy is of tremendous importance, both during drug development and for clinical applications. Importantly, while seminal twin studies indicate that the majority of inter-individual differences in drug disposition are driven by hereditary factors, common genetic polymorphisms explain only less than half of this genetically encoded variability. Recent progress in Next Generation Sequencing (NGS) technologies has for the first time allowed to comprehensively map the genetic landscape of human pharmacogenes. Importantly, these projects have unveiled vast numbers of rare genetic variants, which are estimated to contribute substantially to the missing heritability of drug metabolism phenotypes. However, functional interpretation of these rare variants remains challenging and constitutes one of the important frontiers of contemporary pharmacogenomics. Furthermore, NGS technologies face challenges in the interrogation of genes residing in complex genomic regions, such as CYP2D6 and HLA genes. We here provide an update of the implementation of pharmacogenomic variations in the clinical setting and present emerging strategies that facilitate the translation of NGS data into clinically useful information. Importantly, we anticipate that these developments will soon result in a paradigm shift of pre-emptive genotyping away from the interrogation to candidate variants and towards the comprehensive profiling of an individuals genotype, thus allowing for a true individualization of patient drug treatment regimens.

Inhibition of Notch signaling attenuates pituitary adenoma growth in Nude mice

Preclinical and clinical studies support that Notch signaling may play an important oncogenic role in cancer, but there is scarce information for pituitary tumors. We therefore undertook a functional study to evaluate Notch participation in pituitary adenoma growth. Tumors generated in Nude mice by subcutaneous GH3 somatolactotrope cell injection were treated in vivo with DAPT, a γ-secretase inhibitor, thus inactivating Notch signaling. This treatment led to pituitary tumor reduction, lower prolactin and GH tumor content and a decrease in angiogenesis. Furthermore, in silico transcriptomic and epigenomic analyses uncovered several tumor suppressor genes related to Notch signaling in pituitary tissue, namely Btg2, Nr4a1, Men1, Zfp36 and Cnot1. Gene evaluation suggested that Btg2, Nr4a1 and Cnot1 may be possible players in GH3 xenograft growth. Btg2 mRNA expression was lower in GH3 tumors compared to the parental line, and DAPT increased its expression levels in the tumor in parallel with the inhibition of its volume. Cnot1 mRNA levels were also increased in the pituitary xenografts by DAPT treatment. And the Nr4a1 gene was lower in tumors compared to the parental line, though not modified by DAPT. Finally, because DAPT in vivo may also be acting on tumor microenvironment, we determined the direct effect of DAPT on GH3 cells in vitro. We found that DAPT decreases the proliferative, secretory and migration potential of GH3 cells. These results position selective interruption of Notch signaling as a potential therapeutic tool in adjuvant treatments for aggressive or resistant pituitary tumors.

1,25-Dihydroxyvitamin D3 enhances glucose-stimulated insulin secretion in mouse and human islets: a role for transcriptional regulation of voltage-gated calcium channels by the vitamin D receptor

AIM

Vitamin D deficiency in rodents negatively affects glucose-stimulated insulin secretion (GSIS) and human epidemiological studies connect poor vitamin D status with type 2 diabetes. Previous studies performed primarily in rat islets have shown that vitamin D can enhance GSIS. However the molecular pathways linking vitamin D and insulin secretion are currently unknown. Therefore, experiments were undertaken to elucidate the transcriptional role(s) of the vitamin D receptor (VDR) in islet function.

METHODS

Human and mouse islets were cultured with vehicle or 1,25-dihydroxyvitamin-D₃ (1,25D₃) and then subjected to GSIS assays. Insulin expression, insulin content, glucose uptake and glucose-stimulated calcium influx were tested. Microarray analysis was performed. In silico analysis was used to identify VDR response elements (VDRE) within target genes and their activity was tested using reporter assays.

RESULTS

Vdr mRNA is abundant in islets and Vdr expression is glucose-responsive. Preincubation of mouse and human islets with 1,25D₃ enhances GSIS and increases glucose-stimulated calcium influx. Microarray analysis identified the R-type voltage-gated calcium channel (VGCC) gene, Cacna1e, which is highly upregulated by 1,25D₃ in human and mouse islets and contains a conserved VDRE in intron 7. Results from GSIS assays suggest that 1,25D₃ might upregulate a variant of R-type VGCC that is resistant to chemical inhibition.

CONCLUSION

These results suggest that the role of 1,25D₃ in regulating calcium influx acts through the R-Type VGCC during GSIS, thereby modulating the capacity of beta cells to secrete insulin.

Epigenetics of chronic inflammatory diseases

Chronic, noncommunicable, and inflammation-associated diseases remain the largest cause of morbidity and mortality globally and within the United States. This is mainly due to our limited understanding of the molecular mechanisms that underlie these complex pathologies. The available evidence indicates that studies of epigenetics (traditionally defined as the heritable changes to gene expression that are independent of changes to DNA) are significantly advancing our knowledge of these inflammatory conditions. This review will focus on epigenetic studies of three diseases, that are among the most burdensome globally: cardiovascular disease, the number one cause of deaths worldwide, type 2 diabetes and, Alzheimer’s disease. The current status of epigenetic research, including the ability to predict disease risk, and key pathophysiological defects are discussed. The significance of defining the contribution of epigenetic defects to nonresolving inflammation and aging, each associated with these diseases, is highlighted, as these are likely to provide new insights into inflammatory disease pathogenesis.

Genome-wide profiling of histone H3K27 acetylation featured fatty acid signalling in pancreatic beta cells in diet-induced obesity in mice

AIMS/HYPOTHESIS

Epigenetic regulation of gene expression has been implicated in the pathogenesis of obesity and type 2 diabetes. However, detailed information, such as key transcription factors in pancreatic beta cells that mediate environmental effects, is not yet available.

METHODS

To analyse genome-wide cis-regulatory profiles and transcriptome of pancreatic islets derived from a diet-induced obesity (DIO) mouse model, we conducted chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) of histone H3 lysine 27 acetylation (histone H3K27ac) and high-throughput RNA sequencing. Transcription factor-binding motifs enriched in differential H3K27ac regions were examined by de novo motif analysis. For the predicted transcription factors, loss of function experiments were performed by transfecting specific siRNA in INS-1, a rat beta cell line, with and without palmitate treatment. Epigenomic and transcriptional changes of possible target genes were evaluated by ChIP and quantitative RT-PCR.

RESULTS

After long-term feeding with a high-fat diet, C57BL/6J mice were obese and mildly glucose intolerant. Among 39,350 islet cis-regulatory regions, 13,369 and 4610 elements showed increase and decrease in ChIP-Seq signals, respectively, significantly associated with global change in gene expression. Remarkably, increased H3K27ac showed a distinctive genomic localisation, mainly in the proximal-promoter regions, revealing enriched elements for nuclear respiratory factor 1 (NRF1), GA repeat binding protein α (GABPA) and myocyte enhancer factor 2A (MEF2A) by de novo motif analysis, whereas decreased H3K27ac was enriched for v-maf musculoaponeurotic fibrosarcoma oncogene family protein K (MAFK), a known negative regulator of beta cells. By siRNA-mediated knockdown of NRF1, GABPA or MEF2A we found that INS-1 cells exhibited downregulation of fatty acid β-oxidation genes in parallel with decrease in the associated H3K27ac. Furthermore, in line with the epigenome in DIO mice, palmitate treatment caused increase in H3K27ac and induction of β-oxidation genes; these responses were blunted when NRF1, GABPA or MEF2A were suppressed.

CONCLUSIONS/INTERPRETATION

These results suggest novel roles for DNA-binding proteins and fatty acid signalling in obesity-induced epigenomic regulation of beta cell function.

DATA AVAILABILITY

The next-generation sequencing data in the present study were deposited at ArrayExpress. RNA-Seq: Dataset name: ERR2538129 (Control), ERR2538130 (Diet-induced obesity) Repository name and number: E-MTAB-6718 - RNA-Seq of pancreatic islets derived from mice fed a long-term high-fat diet against chow-fed controls. ChIP-Seq: Dataset name: ERR2538131 (Control), ERR2538132 (Diet-induced obesity) Repository name and number: E-MTAB-6719 - H3K27ac ChIP-Seq of pancreatic islets derived from mice fed a long-term high-fat diet (HFD) against chow-fed controls.

Haploinsufficiency of Insm1 Impairs Postnatal Baseline β-Cell Mass

Baseline β-cell mass is established during the early postnatal period when β-cells expand. In this study, we show that heterozygous ablation of Insm1 decreases baseline β-cell mass and subsequently impairs glucose tolerance. When exposed to a high-fat diet or on an ob/ob background, glucose intolerance was more severe in Insm1^+/lacZ mice compared with Insm1^+/+ mice, although no further decrease in the β-cell mass was detected. In islets of early postnatal Insm1^+/lacZ mice, the cell cycle was prolonged in β-cells due to downregulation of the cell cycle gene Ccnd1 Although Insm1 had a low affinity for the Ccnd1 promoter compared with other binding sites, binding affinity was strongly dependent on Insm1 levels. We observed dramatically decreased binding of Insm1 to the Ccnd1 promoter after downregulation of Insm1 expression. Furthermore, downregulation of Ccnd1 resulted in a prolonged cell cycle, and overexpression of Ccnd1 rescued cell cycle abnormalities observed in Insm1-deficient β-cells. We conclude that decreases in Insm1 interfere with β-cell specification during the early postnatal period and impair glucose homeostasis during metabolic stress in adults. Insm1 levels are therefore a factor that can influence the development of diabetes.