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

Evidence from oyster suggests an ancient role for Pdx in regulating insulin gene expression in animals

Hox and ParaHox genes encode transcription factors with similar expression patterns in divergent animals. The Pdx (Xlox) homeobox gene, for example, is expressed in a sharp spatial domain in the endodermal cell layer of the gut in chordates, echinoderms, annelids and molluscs. The significance of comparable gene expression patterns is unclear because it is not known if downstream transcriptional targets are also conserved. Here, we report evidence indicating that a classic transcriptional target of Pdx1 in vertebrates, the insulin gene, is a likely direct target of Pdx in Pacific oyster adults. We show that one insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue. Transcriptomic comparison suggests that this tissue plays a similar role to the vertebrate pancreas. Using ATAC-seq and ChIP, we identify an upstream regulatory element of the cgILP gene which shows binding interaction with cgPdx protein in oyster hepatopancreas and demonstrate, using a cell culture assay, that the oyster Pdx can act as a transcriptional activator through this site, possibly in synergy with NeuroD. These data argue that a classic homeodomain-target gene interaction dates back to the origin of Bilateria.

In vertebrates insulin is a direct transcriptional target of Pdx: the same is true in Pacific oysters and the authors show insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue, showing this gene interaction dates back to the origin of Bilateria.

Glucocorticoid signaling in pancreatic islets modulates gene regulatory programs and genetic risk of type 2 diabetes

Glucocorticoids are key regulators of glucose homeostasis and pancreatic islet function, but the gene regulatory programs driving responses to glucocorticoid signaling in islets and the contribution of these programs to diabetes risk are unknown. In this study we used ATAC-seq and RNA-seq to map chromatin accessibility and gene expression from eleven primary human islet samples cultured in vitro with the glucocorticoid dexamethasone at multiple doses and durations. We identified thousands of accessible chromatin sites and genes with significant changes in activity in response to glucocorticoids. Chromatin sites up-regulated in glucocorticoid signaling were prominently enriched for glucocorticoid receptor binding sites and up-regulated genes were enriched for ion transport and lipid metabolism, whereas down-regulated chromatin sites and genes were enriched for inflammatory, stress response and proliferative processes. Genetic variants associated with glucose levels and T2D risk were enriched in glucocorticoid-responsive chromatin sites, including fine-mapped variants at 51 known signals. Among fine-mapped variants in glucocorticoid-responsive chromatin, a likely casual variant at the 2p21 locus had glucocorticoid-dependent allelic effects on beta cell enhancer activity and affected SIX2 and SIX3 expression. Our results provide a comprehensive map of islet regulatory programs in response to glucocorticoids through which we uncover a role for islet glucocorticoid signaling in mediating genetic risk of T2D.

A Conserved Notochord Enhancer Controls Pancreas Development in Vertebrates

The notochord is an evolutionary novelty in vertebrates that functions as an important signaling center during development. Notochord ablation in chicken has demonstrated that it is crucial for pancreas development; however, the molecular mechanism has not been fully described. Here, we show that in zebrafish, the loss of function of nog2, a Bmp antagonist expressed in the notochord, impairs β cell differentiation, compatible with the antagonistic role of Bmp in β cell differentiation. In addition, we show that nog2 expression in the notochord is induced by at least one notochord enhancer and its loss of function reduces the number of pancreatic progenitors and impairs β cell differentiation. Tracing Nog2 diffusion, we show that Nog2 emanates from the notochord to the pancreas progenitor domain. Finally, we find a notochord enhancer in human and mice Nog genomic landscapes, suggesting that the acquisition of a Nog notochord enhancer occurred early in the vertebrate phylogeny and contributes to the development of complex organs like the pancreas.

CRISPR-based genome editing in primary human pancreatic islet cells

Gene targeting studies in primary human islets could advance our understanding of mechanisms driving diabetes pathogenesis. Here, we demonstrate successful genome editing in primary human islets using clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated protein 9 (Cas9). CRISPR-based targeting efficiently mutated protein-coding exons, resulting in acute loss of islet β-cell regulators, like the transcription factor PDX1 and the KATP channel subunit KIR6.2, accompanied by impaired β-cell regulation and function. CRISPR targeting of non-coding DNA harboring type 2 diabetes (T2D) risk variants revealed changes in ABCC8, SIX2 and SIX3 expression, and impaired β-cell function, thereby linking regulatory elements in these target genes to T2D genetic susceptibility. Advances here establish a paradigm for genetic studies in human islet cells, and reveal regulatory and genetic mechanisms linking non-coding variants to human diabetes risk.

The impact of non-additive genetic associations on age-related complex diseases

Genome-wide association studies (GWAS) are not fully comprehensive, as current strategies typically test only the additive model, exclude the X chromosome, and use only one reference panel for genotype imputation. We implement an extensive GWAS strategy, GUIDANCE, which improves genotype imputation by using multiple reference panels and includes the analysis of the X chromosome and non-additive models to test for association. We apply this methodology to 62,281 subjects across 22 age-related diseases and identify 94 genome-wide associated loci, including 26 previously unreported. Moreover, we observe that 27.7% of the 94 loci are missed if we use standard imputation strategies with a single reference panel, such as HRC, and only test the additive model. Among the new findings, we identify three novel low-frequency recessive variants with odds ratios larger than 4, which need at least a three-fold larger sample size to be detected under the additive model. This study highlights the benefits of applying innovative strategies to better uncover the genetic architecture of complex diseases.

Scrt1, a transcriptional regulator of β-cell proliferation identified by differential chromatin accessibility during islet maturation

Glucose-induced insulin secretion, a hallmark of mature β-cells, is achieved after birth and is preceded by a phase of intense proliferation. These events occurring in the neonatal period are decisive for establishing an appropriate functional β-cell mass that provides the required insulin throughout life. However, key regulators of gene expression involved in functional maturation of β-cells remain to be elucidated. Here, we addressed this issue by mapping open chromatin regions in newborn versus adult rat islets using the ATAC-seq assay. We obtained a genome-wide picture of chromatin accessible sites (~ 100,000) among which 20% were differentially accessible during maturation. An enrichment analysis of transcription factor binding sites identified a group of transcription factors that could explain these changes. Among them, Scrt1 was found to act as a transcriptional repressor and to control β-cell proliferation. Interestingly, Scrt1 expression was controlled by the transcriptional repressor RE-1 silencing transcription factor (REST) and was increased in an in vitro reprogramming system of pancreatic exocrine cells to β-like cells. Overall, this study led to the identification of several known and unforeseen key transcriptional events occurring during β-cell maturation. These findings will help defining new strategies to induce the functional maturation of surrogate insulin-producing cells.

Defining super-enhancer landscape in triple-negative breast cancer by multiomic profiling

Breast cancer is a heterogeneous disease, affecting over 3.5 million women worldwide, yet the functional role of cis-regulatory elements including super-enhancers in different breast cancer subtypes remains poorly characterized. Triple-negative breast cancer (TNBC) is an aggressive subtype of breast cancer with a poor prognosis. Here we apply integrated epigenomic and transcriptomic profiling to uncover super-enhancer heterogeneity between breast cancer subtypes, and provide clinically relevant biological insights towards TNBC. Using CRISPR/Cas9-mediated gene editing, we identify genes that are specifically regulated by TNBC-specific super-enhancers, including FOXC1 and MET, thereby unveiling a mechanism for specific overexpression of the key oncogenes in TNBC. We also identify ANLN as a TNBC-specific gene regulated by super-enhancer. Our studies reveal a TNBC-specific epigenomic landscape, contributing to the dysregulated oncogene expression in breast tumorigenesis.

Neonatal diabetes mutations disrupt a chromatin pioneering function that activates the human insulin gene

Despite the central role of chromosomal context in gene transcription, human noncoding DNA variants are generally studied outside of their genomic location. This limits our understanding of disease-causing regulatory variants. INS promoter mutations cause recessive neonatal diabetes. We show that all INS promoter point mutations in 60 patients disrupt a CC dinucleotide, whereas none affect other elements important for episomal promoter function. To model CC mutations, we humanized an ∼3.1-kb region of the mouse Ins2 gene. This recapitulated developmental chromatin states and cell-specific transcription. A CC mutant allele, however, abrogated active chromatin formation during pancreas development. A search for transcription factors acting through this element revealed that another neonatal diabetes gene product, GLIS3, has a pioneer-like ability to derepress INS chromatin, which is hampered by the CC mutation. Our in vivo analysis, therefore, connects two human genetic defects in an essential mechanism for developmental activation of the INS gene.

Single-cell chromatin accessibility identifies pancreatic islet cell type- and state-specific regulatory programs of diabetes risk

Single-nucleus assay for transposase-accessible chromatin using sequencing (snATAC-seq) creates new opportunities to dissect cell type-specific mechanisms of complex diseases. Since pancreatic islets are central to type 2 diabetes (T2D), we profiled 15,298 islet cells by using combinatorial barcoding snATAC-seq and identified 12 clusters, including multiple alpha, beta and delta cell states. We cataloged 228,873 accessible chromatin sites and identified transcription factors underlying lineage- and state-specific regulation. We observed state-specific enrichment of fasting glucose and T2D genome-wide association studies for beta cells and enrichment for other endocrine cell types. At T2D signals localized to islet-accessible chromatin, we prioritized variants with predicted regulatory function and co-accessibility with target genes. A causal T2D variant rs231361 at the KCNQ1 locus had predicted effects on a beta cell enhancer co-accessible with INS and genome editing in embryonic stem cell-derived beta cells affected INS levels. Together our findings demonstrate the power of single-cell epigenomics for interpreting complex disease genetics.

ZFAND3 Overexpression in the Mouse Liver Improves Glucose Tolerance and Hepatic Insulin Resistance

Genome-wide association studies have identified more than 300 loci associated with type 2 diabetes mellitus; however, the mechanisms underlying their role in type 2 diabetes mellitus susceptibility remain largely unknown. Zinc finger AN1-type domain 3 (ZFAND3), known as testis-expressed sequence 27, is a type 2 diabetes mellitus-susceptibility gene. Limited information is available regarding the physiological role of ZFAND3 in vivo. This study aimed to investigate the association between ZFAND3 and type 2 diabetes mellitus. ZFAND3 was significantly upregulated in the liver of diabetic mice compared to wild-type mice. To overexpress ZFAND3, we generated a ZFAND3-expressing adenovirus (Ad) vector using an improved Ad vector exhibiting significantly lower hepatotoxicity (Ad-ZFAND3). Glucose tolerance was significantly improved in Ad-ZFAND3-treated mice compared to the control Ad-treated mice. ZFAND3 overexpression in the mouse liver also improved insulin resistance. Furthermore, gluconeogenic gene expression was significantly lower in primary mouse hepatocytes transduced with Ad-ZFAND3 than those transduced with the control Ad vector. The present results suggest that ZFAND3 improves glucose tolerance by improving insulin resistance and suppressing gluconeogenesis, serving as a potential novel therapeutic target for type 2 diabetes mellitus.

SeqEnhDL: sequence-based classification of cell type-specific enhancers using deep learning models

Objective

To address the challenge of computational identification of cell type-specific regulatory elements on a genome-wide scale.

Results

We propose SeqEnhDL, a deep learning framework for classifying cell type-specific enhancers based on sequence features. DNA sequences of “strong enhancer” chromatin states in nine cell types from the ENCODE project were retrieved to build and test enhancer classifiers. For any DNA sequence, positional k-mer (k = 5, 7, 9 and 11) fold changes relative to randomly selected non-coding sequences across each nucleotide position were used as features for deep learning models. Three deep learning models were implemented, including multi-layer perceptron (MLP), Convolutional Neural Network (CNN) and Recurrent Neural Network (RNN). All models in SeqEnhDL outperform state-of-the-art enhancer classifiers (including gkm-SVM and DanQ) in distinguishing cell type-specific enhancers from randomly selected non-coding sequences. Moreover, SeqEnhDL can directly discriminate enhancers from different cell types, which has not been achieved by other enhancer classifiers. Our analysis suggests that both enhancers and their tissue-specificity can be accurately identified based on their sequence features. SeqEnhDL is publicly available at https://github.com/wyp1125/SeqEnhDL.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05518-7.

Long non-coding RNA-regulated pathways in pancreatic β cells: Their role in diabetes

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides that have not coding potential, but act as gene expression regulators through several molecular mechanisms. Several studies have identified tons of lncRNAs that are expressed in pancreatic β cells and many of them have been shown to have β cell-specific expression, suggesting a potential role in the regulation of basal β cell functions. Indeed, accumulating evidence based on numerous studies, has highlighted the implication of lncRNAs in the regulation of pancreatic β cell differentiation and proliferation, insulin synthesis and secretion, and apoptosis. In addition, several lncRNAs have shown to be implicated in pancreatic β cell dysfunction linked to different types of diabetes, including type 1 and type 2 diabetes, and monogenic forms of the disease. Pathogenic conditions linked to diabetes (inflammation or lipoglucotoxicity, for example) dysregulate the expression of several lncRNAs, suggesting that changes in lncRNA may alter potentially important pathways for β cell function, and eventually leading to β cell dysfunction and diabetes development. In this sense, functional characterization of some lncRNAs has demonstrated that these non-coding molecules participate in the regulation of several crucial pathways at the pancreatic β cell level, and dysregulation of these pathways leads to pathogenic phenotypes. In this review, we provide an overview of the action mechanisms of functionally characterized lncRNAs in healthy β cells and describe the contribution of some diabetes-associated lncRNAs to pancreatic β cell failure.

Elucidation of disease etiology by trans-layer omics analysis

To date, genome-wide association studies (GWASs) have successfully identified thousands of associations between genetic polymorphisms and human traits. However, the pathways between the associated genotype and phenotype are often poorly understood. The transcriptome, proteome, and metabolome, the omics, are positioned along the pathway and can provide useful information to translate from genotype to phenotype. This review shows useful data resources for connecting each omics and describes how they are combined into a cohesive analysis. Quantitative trait loci (QTL) are useful information for connecting the genome and other omics. QTL represent how much genetic variants have effects on other omics and give us clues to how GWAS risk SNPs affect biological mechanisms. Integration of each omics provides a robust analytical framework for estimating disease causality, discovering drug targets, and identifying disease-associated tissues. Technological advances and the rise of consortia and biobanks have facilitated the analyses of unprecedented data, improving both the quality and quantity of research. Proficient management of these valuable datasets allows discovering novel insights into the genetic background and etiology of complex human diseases and contributing to personalized medicine.

Pancreatic progenitor epigenome maps prioritize type 2 diabetes risk genes with roles in development

Genetic variants associated with type 2 diabetes (T2D) risk affect gene regulation in metabolically relevant tissues, such as pancreatic islets. Here, we investigated contributions of regulatory programs active during pancreatic development to T2D risk. Generation of chromatin maps from developmental precursors throughout pancreatic differentiation of human embryonic stem cells (hESCs) identifies enrichment of T2D variants in pancreatic progenitor-specific stretch enhancers that are not active in islets. Genes associated with progenitor-specific stretch enhancers are predicted to regulate developmental processes, most notably tissue morphogenesis. Through gene editing in hESCs, we demonstrate that progenitor-specific enhancers harboring T2D-associated variants regulate cell polarity genes LAMA1 and CRB2. Knockdown of lama1 or crb2 in zebrafish embryos causes a defect in pancreas morphogenesis and impairs islet cell development. Together, our findings reveal that a subset of T2D risk variants specifically affects pancreatic developmental programs, suggesting that dysregulation of developmental processes can predispose to T2D.

Chromatin 3D interaction analysis of the STARD10 locus unveils FCHSD2 as a regulator of insulin secretion

Using chromatin conformation capture, we show that an enhancer cluster in the STARD10 type 2 diabetes (T2D) locus forms a defined 3-dimensional (3D) chromatin domain. A 4.1-kb region within this locus, carrying 5 T2D-associated variants, physically interacts with CTCF-binding regions and with an enhancer possessing strong transcriptional activity. Analysis of human islet 3D chromatin interaction maps identifies the FCHSD2 gene as an additional target of the enhancer cluster. CRISPR-Cas9-mediated deletion of the variant region, or of the associated enhancer, from human pancreas-derived EndoC-βH1 cells impairs glucose-stimulated insulin secretion. Expression of both STARD10 and FCHSD2 is reduced in cells harboring CRISPR deletions, and lower expression of STARD10 and FCHSD2 is associated, the latter nominally, with the possession of risk variant alleles in human islets. Finally, CRISPR-Cas9-mediated loss of STARD10 or FCHSD2, but not ARAP1, impairs regulated insulin secretion. Thus, multiple genes at the STARD10 locus influence β cell function.

Isoform-Specific Roles of Mutant p63 in Human Diseases

Simple Summary

The protein p63 belongs to the family of the p53 tumor suppressor. Mouse models have, however, shown that it is not a classical tumor suppressor but instead involved in developmental processes. Mutations in the p63 gene cause several developmental defects in human patients characterized by limb deformation, cleft lip/palate, and ectodermal dysplasia due to p63’s role as a master regulator of epidermal development. In addition, p63 plays a key role as a quality control factor in oocytes and p63 mutations can result either in compromised genetic quality control or premature cell death of all oocytes.

Abstract

The p63 gene encodes a master regulator of epidermal commitment, development, and differentiation. Heterozygous mutations in the DNA binding domain cause Ectrodactyly, Ectodermal Dysplasia, characterized by limb deformation, cleft lip/palate, and ectodermal dysplasia while mutations in in the C-terminal domain of the α-isoform cause Ankyloblepharon-Ectodermal defects-Cleft lip/palate (AEC) syndrome, a life-threatening disorder characterized by skin fragility, severe, long-lasting skin erosions, and cleft lip/palate. The molecular disease mechanisms of these syndromes have recently become elucidated and have enhanced our understanding of the role of p63 in epidermal development. Here we review the molecular cause and functional consequences of these p63-mutations for skin development and discuss the consequences of p63 mutations for female fertility.

PDX1LOW MAFALOW β-cells contribute to islet function and insulin release

Transcriptionally mature and immature β-cells co-exist within the adult islet. How such diversity contributes to insulin release remains poorly understood. Here we show that subtle differences in β-cell maturity, defined using PDX1 and MAFA expression, contribute to islet operation. Functional mapping of rodent and human islets containing proportionally more PDX1^HIGH and MAFA^HIGH β-cells reveals defects in metabolism, ionic fluxes and insulin secretion. At the transcriptomic level, the presence of increased numbers of PDX1^HIGH and MAFA^HIGH β-cells leads to dysregulation of gene pathways involved in metabolic processes. Using a chemogenetic disruption strategy, differences in PDX1 and MAFA expression are shown to depend on islet Ca²⁺ signaling patterns. During metabolic stress, islet function can be restored by redressing the balance between PDX1 and MAFA levels across the β-cell population. Thus, preserving heterogeneity in PDX1 and MAFA expression, and more widely in β-cell maturity, might be important for the maintenance of islet function.

Mechanisms controlling pancreatic islet cell function in insulin secretion

Metabolic homeostasis in mammals is tightly regulated by the complementary actions of insulin and glucagon. The secretion of these hormones from pancreatic β-cells and α-cells, respectively, is controlled by metabolic, endocrine, and paracrine regulatory mechanisms and is essential for the control of blood levels of glucose. The deregulation of these mechanisms leads to various pathologies, most notably type 2 diabetes, which is driven by the combined lesions of impaired insulin action and a loss of the normal insulin secretion response to glucose. Glucose stimulates insulin secretion from β-cells in a bi-modal fashion, and new insights about the underlying mechanisms, particularly relating to the second or amplifying phase of this secretory response, have been recently gained. Other recent work highlights the importance of α-cell-produced proglucagon-derived peptides, incretin hormones from the gastrointestinal tract and other dietary components, including certain amino acids and fatty acids, in priming and potentiation of the β-cell glucose response. These advances provide a new perspective for the understanding of the β-cell failure that triggers type 2 diabetes.

The β-Cell Genomic Landscape in T1D: Implications for Disease Pathogenesis

PURPOSE OF REVIEW

Type 1 diabetes (T1D) develops as a consequence of a combination of genetic predisposition and environmental factors. Combined, these events trigger an autoimmune disease that results in progressive loss of pancreatic β cells, leading to insulin deficiency. This article reviews the current knowledge on the genetics of T1D with a specific focus on genetic variation in pancreatic islet regulatory networks and its implication to T1D risk and disease development.

RECENT FINDINGS

Accumulating evidence suggest an active role of β cells in T1D pathogenesis. Based on such observation several studies aimed in mapping T1D risk variants acting at the β cell level. Such studies unravel T1D risk loci shared with type 2 diabetes (T2D) and T1D risk variants potentially interfering with β-cell responses to external stimuli. The characterization of regulatory genomics maps of disease-relevant states and cell types can be used to elucidate the mechanistic role of β cells in the pathogenesis of T1D.

Human Pluripotent Stem Cells to Model Islet Defects in Diabetes

Diabetes mellitus is characterized by elevated levels of blood glucose and is ultimately caused by insufficient insulin production from pancreatic beta cells. Different research models have been utilized to unravel the molecular mechanisms leading to the onset of diabetes. The generation of pancreatic endocrine cells from human pluripotent stem cells constitutes an approach to study genetic defects leading to impaired beta cell development and function. Here, we review the recent progress in generating and characterizing functional stem cell-derived beta cells. We summarize the diabetes disease modeling possibilities that stem cells offer and the challenges that lie ahead to further improve these models.

A Multi-omic Integrative Scheme Characterizes Tissues of Action at Loci Associated with Type 2 Diabetes

Resolving the molecular processes that mediate genetic risk remains a challenge because most disease-associated variants are non-coding and functional characterization of these signals requires knowledge of the specific tissues and cell-types in which they operate. To address this challenge, we developed a framework for integrating tissue-specific gene expression and epigenomic maps to obtain "tissue-of-action" (TOA) scores for each association signal by systematically partitioning posterior probabilities from Bayesian fine-mapping. We applied this scheme to credible set variants for 380 association signals from a recent GWAS meta-analysis of type 2 diabetes (T2D) in Europeans. The resulting tissue profiles underscored a predominant role for pancreatic islets and, to a lesser extent, adipose and liver, particularly among signals with greater fine-mapping resolution. We incorporated resulting TOA scores into a rule-based classifier and validated the tissue assignments through comparison with data from cis-eQTL enrichment, functional fine-mapping, RNA co-expression, and patterns of physiological association. In addition to implicating signals with a single TOA, we found evidence for signals with shared effects in multiple tissues as well as distinct tissue profiles between independent signals within heterogeneous loci. Lastly, we demonstrated that TOA scores can be directly coupled with eQTL colocalization to further resolve effector transcripts at T2D signals. This framework guides mechanistic inference by directing functional validation studies to the most relevant tissues and can gain power as fine-mapping resolution and cell-specific annotations become richer. This method is generalizable to all complex traits with relevant annotation data and is made available as an R package.

In Vivo Reporter Assays Uncover Changes in Enhancer Activity Caused by Type 2 Diabetes-Associated Single Nucleotide Polymorphisms

Many single nucleotide polymorphisms (SNPs) associated with type 2 diabetes overlap with putative endocrine pancreatic enhancers, suggesting that these SNPs modulate enhancer activity and, consequently, gene expression. We performed in vivo mosaic transgenesis assays in zebrafish to quantitatively test the enhancer activity of type 2 diabetes-associated loci. Six out of 10 tested sequences are endocrine pancreatic enhancers. The risk variant of two sequences decreased enhancer activity, while in another two incremented it. One of the latter (rs13266634) locates in an SLC30A8 exon, encoding a tryptophan-to-arginine substitution that decreases SLC30A8 function, which is the canonical explanation for type 2 diabetes risk association. However, other type 2 diabetes-associated SNPs that truncate SLC30A8 confer protection from this disease, contradicting this explanation. Here, we clarify this incongruence, showing that rs13266634 boosts the activity of an overlapping enhancer and suggesting an SLC30A8 gain of function as the cause for the increased risk for the disease. We further dissected the functionality of this enhancer, finding a single nucleotide mutation sufficient to impair its activity. Overall, this work assesses in vivo the importance of disease-associated SNPs in the activity of endocrine pancreatic enhancers, including a poorly explored case where a coding SNP modulates the activity of an enhancer.

Epigenetics in β-cell adaptation and type 2 diabetes

Healthy pancreatic β-cells adapt to systemic insulin resistance to maintain normal blood glucose levels, and a failure of this adaptation leads to type 2 diabetes in humans. While genome-wide association studies have uncovered genetic variants that are associated with type 2 diabetes, it is still insufficient to explain the high prevalence of this disease. Epigenetics is the study of gene expression changes that do not involve DNA sequence alterations such as DNA methylation, histone modification, and non-coding RNAs. Over the last decade, a large number of studies have reported on the role of epigenetics in β-cell biology. In this review, we summarize the epigenetic mechanisms in β-cell adaptation and type 2 diabetes, including alterations in three-dimensional chromatin structure and RNA modifications.

Long noncoding RNAs as novel biomarkers for Type 2 diabetes

Type 2 diabetes (T2D) is a metabolic disease characterized by disordered glucagon secretion, insulin resistance in target tissues, and decreased islet β-cell mass and function. The routine diagnosis was based on measurements of metabolic markers, while genetic risk factors have been considered to increase the probability of predicting the development of the disease. Recent evidence suggests that long noncoding RNAs (lncRNAs) regulate gene expression in various physiological and pathological processes. As increasing lncRNAs are identified in β cells, understanding the regulatory roles of lncRNAs in T2D becomes indispensable. In this review, we discuss the potential role of lncRNAs contributing to β-cell identity and T2D susceptibility, which provide a perspective insight into the development of novel diagnosis biomarkers for T2D.

Transcription factor GLIS3: Critical roles in thyroid hormone biosynthesis, hypothyroidism, pancreatic beta cells and diabetes

GLI-Similar 3 (GLIS3) is a member of the GLIS subfamily of Krüppel-like zinc finger transcription factors that functions as an activator or repressor of gene expression. Study of GLIS3-deficiency in mice and humans revealed that GLIS3 plays a critical role in the regulation of several biological processes and is implicated in the development of various diseases, including hypothyroidism and diabetes. This was supported by genome-wide association studies that identified significant associations of common variants in GLIS3 with increased risk of these pathologies. To obtain insights into the causal mechanisms underlying these diseases, it is imperative to understand the mechanisms by which this protein regulates the development of these pathologies. Recent studies of genes regulated by GLIS3 led to the identification of a number of target genes and have provided important molecular insights by which GLIS3 controls cellular processes. These studies revealed that GLIS3 is essential for thyroid hormone biosynthesis and identified a critical function for GLIS3 in the generation of pancreatic β cells and insulin gene transcription. These observations raised the possibility that the GLIS3 signaling pathway might provide a potential therapeutic target in the management of diabetes, hypothyroidism, and other diseases. To develop such strategies, it will be critical to understand the upstream signaling pathways that regulate the activity, expression and function of GLIS3. Here, we review the recent progress on the molecular mechanisms by which GLIS3 controls key functions in thyroid follicular and pancreatic β cells and how this causally relates to the development of hypothyroidism and diabetes.

Aberrant methylation underlies insulin gene expression in human insulinoma

Human insulinomas are rare, benign, slowly proliferating, insulin-producing beta cell tumors that provide a molecular "recipe" or "roadmap" for pathways that control human beta cell regeneration. An earlier study revealed abnormal methylation in the imprinted p15.5-p15.4 region of chromosome 11, known to be abnormally methylated in another disorder of expanded beta cell mass and function: the focal variant of congenital hyperinsulinism. Here, we compare deep DNA methylome sequencing on 19 human insulinomas, and five sets of normal beta cells. We find a remarkably consistent, abnormal methylation pattern in insulinomas. The findings suggest that abnormal insulin (INS) promoter methylation and altered transcription factor expression create alternative drivers of INS expression, replacing canonical PDX1-driven beta cell specification with a pathological, looping, distal enhancer-based form of transcriptional regulation. Finally, NFaT transcription factors, rather than the canonical PDX1 enhancer complex, are predicted to drive INS transactivation.

Functional Genomics in Pancreatic β Cells: Recent Advances in Gene Deletion and Genome Editing Technologies for Diabetes Research

The inheritance of variants that lead to coding changes in, or the mis-expression of, genes critical to pancreatic beta cell function can lead to alterations in insulin secretion and increase the risk of both type 1 and type 2 diabetes. Recently developed clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) gene editing tools provide a powerful means of understanding the impact of identified variants on cell function, growth, and survival and might ultimately provide a means, most likely after the transplantation of genetically "corrected" cells, of treating the disease. Here, we review some of the disease-associated genes and variants whose roles have been probed up to now. Next, we survey recent exciting developments in CRISPR/Cas9 technology and their possible exploitation for β cell functional genomics. Finally, we will provide a perspective as to how CRISPR/Cas9 technology may find clinical application in patients with diabetes.

Identification of a novel lncRNA (G3R1) regulated by GLIS3 in pancreatic β-cells

Recent advances in high throughput RNA sequencing have revealed that, in addition to messenger RNAs (mRNAs), long non-coding RNAs (lncRNAs) play an important role in the regulation of many cell functions and of organ development. While a number of lncRNAs have been identified in pancreatic islets, their function remains largely undetermined. Here, we identify a novel long ncRNA regulated by the transcription factor GLIS3, which we refer to as GLIS3 regulated 1 (G3R1). This lncRNA was identified for its significant loss of expression in GLIS3 knockout mouse pancreatic islets. G3R1 appears to be specifically expressed in mouse pancreatic β-cells and in a β-cell line (βTC-6). ChIP-seq analysis indicated that GLIS3 and other islet-enriched transcription factors bind near the G3R1 gene, suggesting they directly regulate G3R1 transcription. Similarly, an apparent human homolog of G3R1 displays a similar expression pattern, with additional expression seen in human brain. In order to determine the function of G3R1 in mouse pancreatic β-cells, we utilized CRISPR to develop a knockout mouse where ~80% of G3R1 sequence is deleted. Phenotypic analysis of these mice did not reveal any impairment in β-cell function or glucose regulation, indicating the complexity underlying the study of lncRNA function.

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

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

Straightforward protocol for allele-specific chromatin conformation capture

A key advance in our understanding of gene regulation came with the finding that the genome undergoes three-dimensional nuclear folding in a genetically determined process. This 3D conformation directly influences the association between enhancers and their target promoters. This complex interplay has been proven to be essential for gene regulation, and genetic variants affecting this process have been associated to human diseases. The development of new technologies that quantify these DNA interactions represented a revolution in the field. High throughput techniques like HiC provide a general picture of chromatin topology. However, they often lack resolution to evidence subtle effects that single nucleotide polymorphisms exert over the contacts between cis-regulatory regions and target promoters. Here we propose a cost-efficient approach to perform allele-specific chromatin conformation analysis. As a proof of concept, we analyzed the impact of a common deletion mapping between SIRPB1 promoter and one of its downstream enhancers.

TrxG Complex Catalytic and Non-catalytic Activity Play Distinct Roles in Pancreas Progenitor Specification and Differentiation

Appropriate regulation of genes that coordinate pancreas progenitor proliferation and differentiation is required for pancreas development. Here, we explore the role of H3K4 methylation and the Trithorax group (TrxG) complexes in mediating gene expression during pancreas development. Disruption of TrxG complex assembly, but not catalytic activity, prevented endocrine cell differentiation in pancreas progenitor spheroids. In vivo loss of TrxG catalytic activity in PDX1⁺ cells increased apoptosis and the fraction of progenitors in the G1 phase of the cell cycle. Pancreas progenitors were reallocated to the acinar lineage, primarily at the expense of NEUROG3⁺ endocrine progenitors. Later in development, acinar and endocrine cell numbers were decreased, and increased gene expression variance and reduced terminal marker activation in acinar cells led to their incomplete differentiation. These findings demonstrate that TrxG co-activator activity is required for gene induction, whereas TrxG catalytic activity and H3K4 methylation help maintain transcriptional stability.

Capturing functional epigenomes for insight into metabolic diseases

BACKGROUND

Metabolic diseases such as obesity are known to be driven by both environmental and genetic factors. Although genome-wide association studies of common variants and their impact on complex traits have provided some biological insight into disease etiology, identified genetic variants have been found to contribute only a small proportion to disease heritability, and to map mainly to non-coding regions of the genome. To link variants to function, association studies of cellular traits, such as epigenetic marks, in disease-relevant tissues are commonly applied.

SCOPE OF THE REVIEW

We review large-scale efforts to generate genome-wide maps of coordinated epigenetic marks and their utility in complex disease dissection with a focus on DNA methylation. We contrast DNA methylation profiling methods and discuss the advantages of using targeted methods for single-base resolution assessments of methylation levels across tissue-specific regulatory regions to deepen our understanding of contributing factors leading to complex diseases.

MAJOR CONCLUSIONS

Large-scale assessments of DNA methylation patterns in metabolic disease-linked study cohorts have provided insight into the impact of variable epigenetic variants in disease etiology. In-depth profiling of epigenetic marks at regulatory regions, particularly at tissue-specific elements, will be key to dissect the genetic and environmental components contributing to metabolic disease onset and progression.

A Long Non-coding RNA, LOC157273, Is an Effector Transcript at the Chromosome 8p23.1-PPP1R3B Metabolic Traits and Type 2 Diabetes Risk Locus

AIMS

Causal transcripts at genomic loci associated with type 2 diabetes (T2D) are mostly unknown. The chr8p23.1 variant rs4841132, associated with an insulin-resistant diabetes risk phenotype, lies in the second exon of a long non-coding RNA (lncRNA) gene, LOC157273, located 175 kilobases from PPP1R3B, which encodes a key protein regulating insulin-mediated hepatic glycogen storage in humans. We hypothesized that LOC157273 regulates expression of PPP1R3B in human hepatocytes.

METHODS

We tested our hypothesis using Stellaris fluorescent in situ hybridization to assess subcellular localization of LOC157273; small interfering RNA (siRNA) knockdown of LOC157273, followed by RT-PCR to quantify LOC157273 and PPP1R3B expression; RNA-seq to quantify the whole-transcriptome gene expression response to LOC157273 knockdown; and an insulin-stimulated assay to measure hepatocyte glycogen deposition before and after knockdown.

RESULTS

We found that siRNA knockdown decreased LOC157273 transcript levels by approximately 80%, increased PPP1R3B mRNA levels by 1.7-fold, and increased glycogen deposition by >50% in primary human hepatocytes. An A/G heterozygous carrier (vs. three G/G carriers) had reduced LOC157273 abundance due to reduced transcription of the A allele and increased PPP1R3B expression and glycogen deposition.

CONCLUSION

We show that the lncRNA LOC157273 is a negative regulator of PPP1R3B expression and glycogen deposition in human hepatocytes and a causal transcript at an insulin-resistant T2D risk locus.

Roles of Noncoding RNAs in Islet Biology

The discovery that most mammalian genome sequences are transcribed to ribonucleic acids (RNA) has revolutionized our understanding of the mechanisms governing key cellular processes and of the causes of human diseases, including diabetes mellitus. Pancreatic islet cells were found to contain thousands of noncoding RNAs (ncRNAs), including micro-RNAs (miRNAs), PIWI-associated RNAs, small nucleolar RNAs, tRNA-derived fragments, long non-coding RNAs, and circular RNAs. While the involvement of miRNAs in islet function and in the etiology of diabetes is now well documented, there is emerging evidence indicating that other classes of ncRNAs are also participating in different aspects of islet physiology. The aim of this article will be to provide a comprehensive and updated view of the studies carried out in human samples and rodent models over the past 15 years on the role of ncRNAs in the control of α- and β-cell development and function and to highlight the recent discoveries in the field. We not only describe the role of ncRNAs in the control of insulin and glucagon secretion but also address the contribution of these regulatory molecules in the proliferation and survival of islet cells under physiological and pathological conditions. It is now well established that most cells release part of their ncRNAs inside small extracellular vesicles, allowing the delivery of genetic material to neighboring or distantly located target cells. The role of these secreted RNAs in cell-to-cell communication between β-cells and other metabolic tissues as well as their potential use as diabetes biomarkers will be discussed. © 2020 American Physiological Society. Compr Physiol 10:893-932, 2020.

The New Frontier of Functional Genomics: From Chromatin Architecture and Noncoding RNAs to Therapeutic Targets

Common diseases are complex, multifactorial disorders whose pathogenesis is influenced by the interplay of genetic predisposition and environmental factors. Genome-wide association studies have interrogated genetic polymorphisms across genomes of individuals to test associations between genotype and susceptibility to specific disorders, providing insights into the genetic architecture of several complex disorders. However, genetic variants associated with the susceptibility to common diseases are often located in noncoding regions of the genome, such as tissue-specific enhancers or long noncoding RNAs, suggesting that regulatory elements might play a relevant role in human diseases. Enhancers are cis-regulatory genomic sequences that act in concert with promoters to regulate gene expression in a precise spatiotemporal manner. They can be located at a considerable distance from their cognate target promoters, increasing the difficulty of their identification. Genomes are organized in domains of chromatin folding, namely topologically associating domains (TADs). Identification of enhancer-promoter interactions within TADs has revealed principles of cell-type specificity across several organisms and tissues. The vast majority of mammalian genomes are pervasively transcribed, accounting for a previously unappreciated complexity of the noncoding RNA fraction. Particularly, long noncoding RNAs have emerged as key players for the establishment of chromatin architecture and regulation of gene expression. In this perspective, we describe the new advances in the fields of transcriptomics and genome organization, focusing on the role of noncoding genomic variants in the predisposition of common diseases. Finally, we propose a new framework for the identification of the next generation of pharmacological targets for common human diseases.

Pancreatic β-cells in type 1 and type 2 diabetes mellitus: different pathways to failure

Loss of functional β-cell mass is the key mechanism leading to the two main forms of diabetes mellitus - type 1 diabetes mellitus (T1DM) and type 2 diabetes mellitus (T2DM). Understanding the mechanisms behind β-cell failure is critical to prevent or revert disease. Basic pathogenic differences exist in the two forms of diabetes mellitus; T1DM is immune mediated and T2DM is mediated by metabolic mechanisms. These mechanisms differentially affect early β-cell dysfunction and eventual fate. Over the past decade, major advances have been made in the field, mostly delivered by studies on β-cells in human disease. These advances include studies of islet morphology and human β-cell gene expression in T1DM and T2DM, the identification and characterization of the role of T1DM and T2DM candidate genes at the β-cell level and the endoplasmic reticulum stress signalling that contributes to β-cell failure in T1DM (mostly IRE1 driven) and T2DM (mostly PERK-eIF2α dependent). Here, we review these new findings, focusing on studies performed on human β-cells or on samples obtained from patients with diabetes mellitus.

(Epi)genomic heterogeneity of pancreatic islet function and failure in type 2 diabetes

Background

Pancreatic Islets of Langerhans are heterogeneous tissues consisting of multiple endocrine cell types that carry out distinct yet coordinated roles to regulate blood glucose homeostasis. Islet dysfunction and specifically failure of the beta cells to secrete adequate insulin are known precursors to type 2 diabetes (T2D) onset. However, the exact genetic, (epi)genomic, and environmental mechanisms that contribute to islet failure, and ultimately to T2D pathogenesis, require further elucidation.

Scope of review

This review summarizes efforts and advances in dissection of the complex genetic underpinnings of islet function and resilience in T2D pathogenesis. In this review, we will highlight results of the latest T2D genome-wide association study (GWAS) and discuss how these data are being combined with clinical measures in patients to uncover putative T2D subtypes and with functional (epi)genomic studies in islets to understand the genetic programming of islet cell identity, function, and adaptation. Finally, we discuss new and important opportunities to address major knowledge gaps in our understanding of islet (dys)function in T2D risk and progression.

Major conclusions

Genetic variation exerts clear effects on the islet epigenome, regulatory element usage, and gene expression. Future (epi)genomic comparative analyses between T2D and normal islets should incorporate genetics to distinguish patient-specific from disease-specific differences. Incorporating genotype information into future analyses and studies will also enable more precise insights into the molecular genetics of islet deficiency and failure in T2D risk, and should ultimately contribute to a stratified view of T2D and more precise treatment strategies. Islet cellular heterogeneity continues to remain a challenge for understanding the associations between islet failure and T2D development. Further efforts to obtain purified islet cell type populations and determine the specific genetic and environmental effects on each will help address this. Beyond observation of islets at steady state conditions, more research of islet stress and stimulation responses are needed to understand the transition of these tissues from a healthy to diseased state. Together, focusing on these objectives will provide more opportunities to prevent, treat, and manage T2D.

Polymorphic Inversions Underlie the Shared Genetic Susceptibility of Obesity-Related Diseases

The burden of several common diseases including obesity, diabetes, hypertension, asthma, and depression is increasing in most world populations. However, the mechanisms underlying the numerous epidemiological and genetic correlations among these disorders remain largely unknown. We investigated whether common polymorphic inversions underlie the shared genetic influence of these disorders. We performed an inversion association analysis including 21 inversions and 25 obesity-related traits on a total of 408,898 Europeans and validated the results in 67,299 independent individuals. Seven inversions were associated with multiple diseases while inversions at 8p23.1, 16p11.2, and 11q13.2 were strongly associated with the co-occurrence of obesity with other common diseases. Transcriptome analysis across numerous tissues revealed strong candidate genes for obesity-related traits. Analyses in human pancreatic islets indicated the potential mechanism of inversions in the susceptibility of diabetes by disrupting the cis-regulatory effect of SNPs from their target genes. Our data underscore the role of inversions as major genetic contributors to the joint susceptibility to common complex diseases.

Loss of the transcription factor MAFB limits β-cell derivation from human PSCs

Next generation sequencing studies have highlighted discrepancies in β-cells which exist between mice and men. Numerous reports have identified MAF BZIP Transcription Factor B (MAFB) to be present in human β-cells postnatally, while its expression is restricted to embryonic and neo-natal β-cells in mice. Using CRISPR/Cas9-mediated gene editing, coupled with endocrine cell differentiation strategies, we dissect the contribution of MAFB to β-cell development and function specifically in humans. Here we report that MAFB knockout hPSCs have normal pancreatic differentiation capacity up to the progenitor stage, but favor somatostatin- and pancreatic polypeptide–positive cells at the expense of insulin- and glucagon-producing cells during endocrine cell development. Our results describe a requirement for MAFB late in the human pancreatic developmental program and identify it as a distinguishing transcription factor within islet cell subtype specification. We propose that hPSCs represent a powerful tool to model human pancreatic endocrine development and associated disease pathophysiology.

The MAF bZIP transcription factor B (MAFB) is present in postnatal human beta cells but its role is unclear. Here, the authors show that MAFB regulates endocrine pancreatic cell fate specification.

Skeletal muscle enhancer interactions identify genes controlling whole-body metabolism

Obesity and type 2 diabetes (T2D) are metabolic disorders influenced by lifestyle and genetic factors that are characterized by insulin resistance in skeletal muscle, a prominent site of glucose disposal. Numerous genetic variants have been associated with obesity and T2D, of which the majority are located in non-coding DNA regions. This suggests that most variants mediate their effect by altering the activity of gene-regulatory elements, including enhancers. Here, we map skeletal muscle genomic enhancer elements that are dynamically regulated after exposure to the free fatty acid palmitate or the inflammatory cytokine TNFα. By overlapping enhancer positions with the location of disease-associated genetic variants, and resolving long-range chromatin interactions between enhancers and gene promoters, we identify target genes involved in metabolic dysfunction in skeletal muscle. The majority of these genes also associate with altered whole-body metabolic phenotypes in the murine BXD genetic reference population. Thus, our combined genomic investigations identified genes that are involved in skeletal muscle metabolism.

Obesity and type 2 diabetes (T2D) are metabolic disorders characterized by insulin resistance in skeletal muscle. Here, the authors map skeletal muscle enhancer elements dynamically regulated after exposure to free fatty acid palmitate or inflammatory cytokine TNFα and identify target genes involved in metabolic dysfunction in skeletal muscle.

Genomic resources for dissecting the role of non-protein coding variation in gene-environment interactions

The majority of single nucleotide variants (SNVs) identified in Genome Wide Association Studies (GWAS) fall within non-protein coding DNA and have the potential to alter gene expression. Non-protein coding DNA can control gene expression by acting as transcription factor (TF) binding sites or by regulating the organization of DNA into chromatin. SNVs in non-coding DNA sequences can disrupt TF binding and chromatin structure and this can result in pathology. Further, environmental health studies have shown that exposure to xenobiotics can disrupt the ability of TFs to regulate entire gene networks and result in pathology. However, there is a large amount of interindividual variability in exposure-linked health outcomes. One explanation for this heterogeneity is that genetic variation and exposure combine to disrupt gene regulation, and this eventually manifests in disease. Many resources exist that annotate common variants from GWAS and combine them with conservation, functional genomics, and TF binding data. These annotation tools provide clues regarding the biological implications of an SNV, as well as lead to the generation of hypotheses regarding potentially disrupted target genes, epigenetic markers, pathways, and cell types. Collectively this information can be used to predict how SNVs can alter an individual's response to exposure and disease risk. A basic understanding of the regulatory information contained within non-protein coding DNA is needed to predict the biological consequences of SNVs, and to determine how these SNVs impact exposure-related disease. We hope that this review will aid in the characterization of disease-associated genetic variation in the non-protein coding genome.

ENdb: a manually curated database of experimentally supported enhancers for human and mouse

Enhancers are a class of cis-regulatory elements that can increase gene transcription by forming loops in intergenic regions, introns and exons. Enhancers, as well as their associated target genes, and transcription factors (TFs) that bind to them, are highly associated with human disease and biological processes. Although some enhancer databases have been published, most only focus on enhancers identified by high-throughput experimental techniques. Therefore, it is highly desirable to construct a comprehensive resource of manually curated enhancers and their related information based on low-throughput experimental evidences. Here, we established a comprehensive manually-curated enhancer database for human and mouse, which provides a resource for experimentally supported enhancers, and to annotate the detailed information of enhancers. The current release of ENdb documents 737 experimentally validated enhancers and their related information, including 384 target genes, 263 TFs, 110 diseases and 153 functions in human and mouse. Moreover, the enhancer-related information was supported by experimental evidences, such as RNAi, in vitro knockdown, western blotting, qRT-PCR, luciferase reporter assay, chromatin conformation capture (3C) and chromosome conformation capture-on-chip (4C) assays. ENdb provides a user-friendly interface to query, browse and visualize the detailed information of enhancers. The database is available at http://www.licpathway.net/ENdb.

From GWAS to Function: Using Functional Genomics to Identify the Mechanisms Underlying Complex Diseases

Genome-wide association studies (GWAS) have successfully mapped thousands of loci associated with complex traits. These associations could reveal the molecular mechanisms altered in common complex diseases and result in the identification of novel drug targets. However, GWAS have also left a number of outstanding questions. In particular, the majority of disease-associated loci lie in non-coding regions of the genome and, even though they are thought to play a role in gene expression regulation, it is unclear which genes they regulate and in which cell types or physiological contexts this regulation occurs. This has hindered the translation of GWAS findings into clinical interventions. In this review we summarize how these challenges have been addressed over the last decade, with a particular focus on the integration of GWAS results with functional genomics datasets. Firstly, we investigate how the tissues and cell types involved in diseases can be identified using methods that test for enrichment of GWAS variants in genomic annotations. Secondly, we explore how to find the genes regulated by GWAS loci using methods that test for colocalization of GWAS signals with molecular phenotypes such as quantitative trait loci (QTLs). Finally, we highlight potential future research avenues such as integrating GWAS results with single-cell sequencing read-outs, designing functionally informed polygenic risk scores (PRS), and validating disease associated genes using genetic engineering. These tools will be crucial to identify new drug targets for common complex diseases.

The glucocorticoid receptor interferes with progesterone receptor-dependent genomic regulation in breast cancer cells

The glucocorticoid and progesterone receptors (GR and PR) are closely related members of the steroid receptor family. Despite sharing similar structural and functional characteristics; the cognate hormones display very distinct physiological responses. In mammary epithelial cells, PR activation is associated with the incidence and progression of breast cancer, whereas the GR is related to growth suppression and differentiation. Despite their pharmacological relevance, only a few studies have compared GR and PR activities in the same system. Using a PR⁺/GR⁺ breast cancer cell line, here we report that either glucocorticoid-free or dexamethasone (DEX)-activated GR inhibits progestin-dependent gene expression associated to epithelial-mesenchymal-transition and cell proliferation. When both receptors are activated with their cognate hormones, PR and GR can form part of the same complex according to co-immunoprecipitation, quantitative microscopy and sequential ChIP experiments. Moreover, genome-wide studies in cells treated with either DEX or R5020, revealed the presence of several regions co-bound by both receptors. Surprisingly, GR also binds novel genomic sites in cells treated with R5020 alone. This progestin-induced GR binding was enriched in REL DNA motifs and located close to genes coding for chromatin remodelers. Understanding GR behavior in the context of progestin-dependent breast cancer could provide new targets for tumor therapy.

Missing heritability in Parkinson's disease: the emerging role of non-coding genetic variation

Parkinson's disease (PD) is a neurodegenerative disorder caused by a complex interplay of genetic and environmental factors. For the stratification of PD patients and the development of advanced clinical trials, including causative treatments, a better understanding of the underlying genetic architecture of PD is required. Despite substantial efforts, genome-wide association studies have not been able to explain most of the observed heritability. The majority of PD-associated genetic variants are located in non-coding regions of the genome. A systematic assessment of their functional role is hampered by our incomplete understanding of genotype-phenotype correlations, for example through differential regulation of gene expression. Here, the recent progress and remaining challenges for the elucidation of the role of non-coding genetic variants is reviewed with a focus on PD as a complex disease with multifactorial origins. The function of gene regulatory elements and the impact of non-coding variants on them, and the means to map these elements on a genome-wide level, will be delineated. Moreover, examples of how the integration of functional genomic annotations can serve to identify disease-associated pathways and to prioritize disease- and cell type-specific regulatory variants will be given. Finally, strategies for functional validation and considerations for suitable model systems are outlined. Together this emphasizes the contribution of rare and common genetic variants to the complex pathogenesis of PD and points to remaining challenges for the dissection of genetic complexity that may allow for better stratification, improved diagnostics and more targeted treatments for PD in the future.

Multiomic Profiling Identifies cis-Regulatory Networks Underlying Human Pancreatic β Cell Identity and Function

EndoC-βH1 is emerging as a critical human β cell model to study the genetic and environmental etiologies of β cell (dys)function and diabetes. Comprehensive knowledge of its molecular landscape is lacking, yet required, for effective use of this model. Here, we report chromosomal (spectral karyotyping), genetic (genotyping), epigenomic (ChIP-seq and ATAC-seq), chromatin interaction (Hi-C and Pol2 ChIA-PET), and transcriptomic (RNA-seq and miRNA-seq) maps of EndoC-βH1. Analyses of these maps define known (e.g., PDX1 and ISL1) and putative (e.g., PCSK1 and mir-375) β cell-specific transcriptional cis-regulatory networks and identify allelic effects on cis-regulatory element use. Importantly, comparison with maps generated in primary human islets and/or β cells indicates preservation of chromatin looping but also highlights chromosomal aberrations and fetal genomic signatures in EndoC-βH1. Together, these maps, and a web application we created for their exploration, provide important tools for the design of experiments to probe and manipulate the genetic programs governing β cell identity and (dys)function in diabetes.

EndoC-βH1 is becoming an important cellular model to study genes and pathways governing human β cell identity and function, but its (epi)genomic similarity to primary human islets is unknown. Lawlor et al. complete and compare extensive EndoC and primary human islet multiomic maps to identify shared and distinct genomic circuitry.

Insights into pancreatic islet cell dysfunction from type 2 diabetes mellitus genetics

Type 2 diabetes mellitus (T2DM) is an increasingly prevalent multifactorial disease that has both genetic and environmental risk factors, resulting in impaired glucose homeostasis. Genome-wide association studies (GWAS) have identified over 400 genetic signals that are associated with altered risk of T2DM. Human physiology and epigenomic data support a central role for the pancreatic islet in the pathogenesis of T2DM. This Review focuses on the promises and challenges of moving from genetic associations to molecular mechanisms and highlights efforts to identify the causal variant and effector transcripts at T2DM GWAS susceptibility loci. In addition, we examine current human models that are used to study both β-cell development and function, including EndoC-β cell lines and human induced pluripotent stem cell-derived β-like cells. We use examples of four T2DM susceptibility loci (CDKAL1, MTNR1B, SLC30A8 and PAM) to emphasize how a holistic approach involving genetics, physiology, and cellular and developmental biology can disentangle disease mechanisms at T2DM GWAS signals.

Characterization of functional transposable element enhancers in acute myeloid leukemia

Transposable elements (TEs) have been shown to have important gene regulatory functions and their alteration could lead to disease phenotypes. Acute myeloid leukemia (AML) develops as a consequence of a series of genetic changes in hematopoietic precursor cells, including mutations in epigenetic factors. Here, we set out to study the gene regulatory role of TEs in AML. We first explored the epigenetic landscape of TEs in AML patients using ATAC-seq data. We show that a large number of TEs in general, and more specifically mammalian-wide interspersed repeats (MIRs), are more enriched in AML cells than in normal blood cells. We obtained a similar finding when analyzing histone modification data in AML patients. Gene Ontology enrichment analysis showed that genes near MIRs in open chromatin regions are involved in leukemogenesis. To functionally validate their regulatory role, we selected 19 MIR regions in AML cells, and tested them for enhancer activity in an AML cell line (Kasumi-1) and a chronic myeloid leukemia (CML) cell line (K562); the results revealed several MIRs to be functional enhancers. Taken together, our results suggest that TEs are potentially involved in myeloid leukemogenesis and highlight these sequences as potential candidates harboring AML-associated variation.

From Genetic Association to Molecular Mechanisms for Islet-cell Dysfunction in Type 2 Diabetes

Genome-wide association studies (GWAS) have identified over 400 signals robustly associated with risk for type 2 diabetes (T2D). At the vast majority of these loci, the lead single nucleotide polymorphisms (SNPs) reside in noncoding regions of the genome, which hampers biological inference and translation of genetic discoveries into disease mechanisms. The study of these T2D risk variants in normoglycemic individuals has revealed that a significant proportion are exerting their disease risk through islet-cell dysfunction. The central role of the islet is also demonstrated by numerous studies, which have shown an enrichment of these signals in islet-specific epigenomic annotations. In recent years the emergence of authentic human beta-cell lines, and advances in genome-editing technologies coupled with improved protocols differentiating human pluripotent stem cells into beta-like cells has opened up new opportunities for T2D disease modeling. Here we review the current understanding on the genetic basis of T2D focusing on approaches, which have facilitated the identification of causal variants and their effector transcripts in human islets. We will present examples of functional studies based on animal and conventional cellular systems and highlight the potential of novel stem cell-based T2D disease models.

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.