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

Single-housing-induced islet epigenomic changes are related to polymorphisms in diabetic KK mice

A lack of social relationships is increasingly recognized as a type 2 diabetes (T2D) risk. To investigate the underlying mechanism, we used male KK mice, an inbred strain with spontaneous diabetes. Given the association between living alone and T2D risk in humans, we divided the non-diabetic mice into singly housed (KK-SH) and group-housed control mice. Around the onset of diabetes in KK-SH mice, we compared H3K27ac ChIP-Seq with RNA-Seq using pancreatic islets derived from each experimental group, revealing a positive correlation between single-housing-induced changes in H3K27ac and gene expression levels. In particular, single-housing-induced H3K27ac decreases revealed a significant association with islet cell functions and GWAS loci for T2D and related diseases, with significant enrichment of binding motifs for transcription factors representative of human diabetes. Although these H3K27ac regions were preferentially localized to a polymorphic genomic background, SNVs and indels did not cause sequence disruption of enriched transcription factor motifs in most of these elements. These results suggest alternative roles of genetic variants in environment-dependent epigenomic changes and provide insights into the complex mode of disease inheritance.

Untangling the genetics of beta cell dysfunction and death in type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is a complex multi-system disease which arises from both environmental and genetic factors, resulting in the destruction of insulin-producing pancreatic beta cells. Over the past two decades, human genetic studies have provided new insight into the etiology of T1D, including an appreciation for the role of beta cells in their own demise.

SCOPE OF REVIEW

Here, we outline models supported by human genetic data for the role of beta cell dysfunction and death in T1D. We highlight the importance of strong evidence linking T1D genetic associations to bona fide candidate genes for mechanistic and therapeutic consideration. To guide rigorous interpretation of genetic associations, we describe molecular profiling approaches, genomic resources, and disease models that may be used to construct variant-to-gene links and to investigate candidate genes and their role in T1D.

MAJOR CONCLUSIONS

We profile advances in understanding the genetic causes of beta cell dysfunction and death at individual T1D risk loci. We discuss how genetic risk prediction models can be used to address disease heterogeneity. Further, we present areas where investment will be critical for the future use of genetics to address open questions in the development of new treatment and prevention strategies for T1D.

Super-enhancer omics in stem cell

The hallmarks of stem cells, such as proliferation, self-renewal, development, differentiation, and regeneration, are critical to maintain stem cell identity which is sustained by genetic and epigenetic factors. Super-enhancers (SEs), which consist of clusters of active enhancers, play a central role in maintaining stemness hallmarks by specifically transcriptional model. The SE-navigated transcriptional complex, including SEs, non-coding RNAs, master transcriptional factors, Mediators and other co-activators, forms phase-separated condensates, which offers a toggle for directing diverse stem cell fate. With the burgeoning technologies of multiple-omics applied to examine different aspects of SE, we firstly raise the concept of "super-enhancer omics", inextricably linking to Pan-omics. In the review, we discuss the spatiotemporal organization and concepts of SEs, and describe links between SE-navigated transcriptional complex and stem cell features, such as stem cell identity, self-renewal, pluripotency, differentiation and development. We also elucidate the mechanism of stemness and oncogenic SEs modulating cancer stem cells via genomic and epigenetic alterations hijack in cancer stem cell. Additionally, we discuss the potential of targeting components of the SE complex using small molecule compounds, genome editing, and antisense oligonucleotides to treat SE-associated organ dysfunction and diseases, including cancer. This review also provides insights into the future of stem cell research through the paradigm of SEs.

Atypical KCNQ1/Kv7 channel function in a neonatal diabetes patient: Hypersecretion preceded the failure of pancreatic β-cells

KCNQ1/Kv7, a low-voltage-gated K+ channel, regulates cardiac rhythm and glucose homeostasis. While KCNQ1 mutations are associated with long-QT syndrome and type2 diabetes, its function in human pancreatic cells remains controversial. We identified a homozygous KCNQ1 mutation (R397W) in an individual with permanent neonatal diabetes melitus (PNDM) without cardiovascular symptoms. To decipher the potential mechanism(s), we introduced the mutation into human embryonic stem cells and generated islet-like organoids (SC-islets) using CRISPR-mediated homology-repair. The mutation did not affect pancreatic differentiation, but affected channel function by increasing spike frequency and Ca2+ flux, leading to insulin hypersecretion. With prolonged culturing, the mutant islets decreased their secretion and gradually deteriorated, modeling a diabetic state, which accelerated by high glucose levels. The molecular basis was the downregulated expression of voltage-activated Ca2+ channels and oxidative phosphorylation. Our study provides a better understanding of the role of KCNQ1 in regulating insulin secretion and β-cell survival in hereditary diabetes pathology.

MethNet: a robust approach to identify regulatory hubs and their distal targets from cancer data

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover cis regulatory elements (CREs) in a 1 Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPR interference based single-cell RNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

Ca2+ signaling and metabolic stress-induced pancreatic β-cell failure

Early in the development of Type 2 diabetes (T2D), metabolic stress brought on by insulin resistance and nutrient overload causes β-cell hyperstimulation. Herein we summarize recent studies that have explored the premise that an increase in the intracellular Ca2+ concentration ([Ca2+]i), brought on by persistent metabolic stimulation of β-cells, causes β-cell dysfunction and failure by adversely affecting β-cell function, structure, and identity. This mini-review builds on several recent reviews that also describe how excess [Ca2+]i impairs β-cell function.

Enhancer-promoter hubs organize transcriptional networks promoting oncogenesis and drug resistance

Recent advances in high-resolution mapping of spatial interactions among regulatory elements support the existence of complex topological assemblies of enhancers and promoters known as enhancer-promoter hubs or cliques. Yet, organization principles of these multi-interacting enhancer-promoter hubs and their potential role in regulating gene expression in cancer remains unclear. Here, we systematically identified enhancer-promoter hubs in breast cancer, lymphoma, and leukemia. We found that highly interacting enhancer-promoter hubs form at key oncogenes and lineage-associated transcription factors potentially promoting oncogenesis of these diverse cancer types. Genomic and optical mapping of interactions among enhancer and promoter elements further showed that topological alterations in hubs coincide with transcriptional changes underlying acquired resistance to targeted therapy in T cell leukemia and B cell lymphoma. Together, our findings suggest that enhancer-promoter hubs are dynamic and heterogeneous topological assemblies with the potential to control gene expression circuits promoting oncogenesis and drug resistance.

Chromosome 20p11.2 deletions cause congenital hyperinsulinism via the loss of FOXA2 or its regulatory elements

Persistent congenital hyperinsulinism (HI) is a rare genetically heterogeneous condition characterised by dysregulated insulin secretion leading to life-threatening hypoglycaemia. For up to 50% of affected individuals screening of the known HI genes does not identify a disease-causing variant. Large deletions have previously been used to identify novel regulatory regions causing HI. Here, we used genome sequencing to search for novel large (>1 Mb) deletions in 180 probands with HI of unknown cause and replicated our findings in a large cohort of 883 genetically unsolved individuals with HI using off-target copy number variant calling from targeted gene panels. We identified overlapping heterozygous deletions in five individuals (range 3-8 Mb) spanning chromosome 20p11.2. The pancreatic beta-cell transcription factor gene, FOXA2, a known cause of HI was deleted in two of the five individuals. In the remaining three, we found a minimal deleted region of 2.4 Mb adjacent to FOXA2 that encompasses multiple non-coding regulatory elements that are in conformational contact with FOXA2. Our data suggests that the deletions in these three children may cause disease through the dysregulation of FOXA2 expression. These findings provide new insights into the regulation of FOXA2 in the beta-cell and confirm an aetiological role for chromosome 20p11.2 deletions in syndromic HI.

Senescence of human pancreatic beta cells enhances functional maturation through chromatin reorganization and promotes interferon responsiveness

Senescent cells can influence the function of tissues in which they reside, and their propensity for disease. A portion of adult human pancreatic beta cells express the senescence marker p16, yet it is unclear whether they are in a senescent state, and how this affects insulin secretion. We analyzed single-cell transcriptome datasets of adult human beta cells, and found that p16-positive cells express senescence gene signatures, as well as elevated levels of beta-cell maturation genes, consistent with enhanced functionality. Senescent human beta-like cells in culture undergo chromatin reorganization that leads to activation of enhancers regulating functional maturation genes and acquisition of glucose-stimulated insulin secretion capacity. Strikingly, Interferon-stimulated genes are elevated in senescent human beta cells, but genes encoding senescence-associated secretory phenotype (SASP) cytokines are not. Senescent beta cells in culture and in human tissue show elevated levels of cytoplasmic DNA, contributing to their increased interferon responsiveness. Human beta-cell senescence thus involves chromatin-driven upregulation of a functional-maturation program, and increased responsiveness of interferon-stimulated genes, changes that could increase both insulin secretion and immune reactivity.

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic β-cell identity and function. Elimination of Lkb1 from the β-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of β cell- selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from β-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors important for β-cell identity, such as FOXA, MAFA or RFX6 and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates β-cell identity and function.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic β-cell survival and function

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, we examined multiple variants that influence SLC30A8 allele-specific expression. Epigenomic mapping has previously identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighboring genes. Here, we show that deletion of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-βH3 cells lowers the expression of SLC30A8 and several neighboring genes and improves glucose-stimulated insulin secretion. While downregulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21, or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Constructing eRNA-mediated gene regulatory networks to explore the genetic basis of muscle and fat-relevant traits in pigs

BACKGROUND

Enhancer RNAs (eRNAs) play a crucial role in transcriptional regulation. While significant progress has been made in understanding epigenetic regulation mediated by eRNAs, research on the construction of eRNA-mediated gene regulatory networks (eGRN) and the identification of critical network components that influence complex traits is lacking.

RESULTS

Here, employing the pig as a model, we conducted a comprehensive study using H3K27ac histone ChIP-seq and RNA-seq data to construct eRNA expression profiles from multiple tissues of two distinct pig breeds, namely Enshi Black (ES) and Duroc. In addition to revealing the regulatory landscape of eRNAs at the tissue level, we developed an innovative network construction and refinement method by integrating RNA-seq, ChIP-seq, genome-wide association study (GWAS) signals and enhancer-modulating effects of single nucleotide polymorphisms (SNPs) measured by self-transcribing active regulatory region sequencing (STARR-seq) experiments. Using this approach, we unraveled eGRN that significantly influence the growth and development of muscle and fat tissues, and identified several novel genes that affect adipocyte differentiation in a cell line model.

CONCLUSIONS

Our work not only provides novel insights into the genetic basis of economic pig traits, but also offers a generalizable approach to elucidate the eRNA-mediated transcriptional regulation underlying a wide spectrum of complex traits for diverse organisms.

RFX6 Maintains Gene Expression and Function of Adult Human Islet α-Cells

Mutations in the gene encoding the transcription factor regulatory factor X-box binding 6 (RFX6) are associated with human diabetes. Within pancreatic islets, RFX6 expression is most abundant in islet α-cells, and α-cell RFX6 expression is altered in diabetes. However, the roles of RFX6 in regulating gene expression, glucagon output, and other crucial human adult α-cell functions are not yet understood. We developed a method for selective genetic targeting of human α-cells and assessed RFX6-dependent α-cell function. RFX6 suppression with RNA interference led to impaired α-cell exocytosis and dysregulated glucagon secretion in vitro and in vivo. By contrast, these phenotypes were not observed with RFX6 suppression across all islet cells. Transcriptomics in α-cells revealed RFX6-dependent expression of genes governing nutrient sensing, hormone processing, and secretion, with some of these exclusively expressed in human α-cells. Mapping of RFX6 DNA-binding sites in primary human islet cells identified a subset of direct RFX6 target genes. Together, these data unveil RFX6-dependent genetic targets and mechanisms crucial for regulating adult human α-cell function.

ARTICLE HIGHLIGHTS

Lessons and Applications of Omics Research in Diabetes Epidemiology

PURPOSE OF REVIEW

Recent advances in genomic technology and molecular techniques have greatly facilitated the identification of disease biomarkers, advanced understanding of pathogenesis of different common diseases, and heralded the dawn of precision medicine. Much of these advances in the area of diabetes have been made possible through deep phenotyping of epidemiological cohorts, and analysis of the different omics data in relation to detailed clinical information. In this review, we aim to provide an overview on how omics research could be incorporated into the design of current and future epidemiological studies.

RECENT FINDINGS

We provide an up-to-date review of the current understanding in the area of genetic, epigenetic, proteomic and metabolomic markers for diabetes and related outcomes, including polygenic risk scores. We have drawn on key examples from the literature, as well as our own experience of conducting omics research using the Hong Kong Diabetes Register and Hong Kong Diabetes Biobank, as well as other cohorts, to illustrate the potential of omics research in diabetes. Recent studies highlight the opportunity, as well as potential benefit, to incorporate molecular profiling in the design and set-up of diabetes epidemiology studies, which can also advance understanding on the heterogeneity of diabetes. Learnings from these examples should facilitate other researchers to consider incorporating research on omics technologies into their work to advance the field and our understanding of diabetes and its related co-morbidities. Insights from these studies would be important for future development of precision medicine in diabetes.

Priority index for critical Covid-19 identifies clinically actionable targets and drugs

While genome-wide studies have identified genomic loci in hosts associated with life-threatening Covid-19 (critical Covid-19), the challenge of resolving these loci hinders further identification of clinically actionable targets and drugs. Building upon our previous success, we here present a priority index solution designed to address this challenge, generating the target and drug resource that consists of two indexes: the target index and the drug index. The primary purpose of the target index is to identify clinically actionable targets by prioritising genes associated with Covid-19. We illustrate the validity of the target index by demonstrating its ability to identify pre-existing Covid-19 phase-III drug targets, with the majority of these targets being found at the leading prioritisation (leading targets). These leading targets have their evolutionary origins in Amniota ('four-leg vertebrates') and are predominantly involved in cytokine-cytokine receptor interactions and JAK-STAT signaling. The drug index highlights opportunities for repurposing clinically approved JAK-STAT inhibitors, either individually or in combination. This proposed strategic focus on the JAK-STAT pathway is supported by the active pursuit of therapeutic agents targeting this pathway in ongoing phase-II/III clinical trials for Covid-19.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying β cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR-Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR-Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for Insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired β cell PCSK1 regulation and Insulin secretion. Multiplex CRISPR-Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

Repression of latent NF-κB enhancers by PDX1 regulates β cell functional heterogeneity

Interactions between lineage-determining and activity-dependent transcription factors determine single-cell identity and function within multicellular tissues through incompletely known mechanisms. By assembling a single-cell atlas of chromatin state within human islets, we identified β cell subtypes governed by either high or low activity of the lineage-determining factor pancreatic duodenal homeobox-1 (PDX1). β cells with reduced PDX1 activity displayed increased chromatin accessibility at latent nuclear factor κB (NF-κB) enhancers. Pdx1 hypomorphic mice exhibited de-repression of NF-κB and impaired glucose tolerance at night. Three-dimensional analyses in tandem with chromatin immunoprecipitation (ChIP) sequencing revealed that PDX1 silences NF-κB at circadian and inflammatory enhancers through long-range chromatin contacts involving SIN3A. Conversely, Bmal1 ablation in β cells disrupted genome-wide PDX1 and NF-κB DNA binding. Finally, antagonizing the interleukin (IL)-1β receptor, an NF-κB target, improved insulin secretion in Pdx1 hypomorphic islets. Our studies reveal functional subtypes of single β cells defined by a gradient in PDX1 activity and identify NF-κB as a target for insulinotropic therapy.

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

Genes with epigenetic alterations in human pancreatic islets impact mitochondrial function, insulin secretion, and type 2 diabetes

Epigenetic dysregulation may influence disease progression. Here we explore whether epigenetic alterations in human pancreatic islets impact insulin secretion and type 2 diabetes (T2D). In islets, 5,584 DNA methylation sites exhibit alterations in T2D cases versus controls and are associated with HbA1c in individuals not diagnosed with T2D. T2D-associated methylation changes are found in enhancers and regions bound by β-cell-specific transcription factors and associated with reduced expression of e.g. CABLES1, FOXP1, GABRA2, GLR1A, RHOT1, and TBC1D4. We find RHOT1 (MIRO1) to be a key regulator of insulin secretion in human islets. Rhot1-deficiency in β-cells leads to reduced insulin secretion, ATP/ADP ratio, mitochondrial mass, Ca2+, and respiration. Regulators of mitochondrial dynamics and metabolites, including L-proline, glycine, GABA, and carnitines, are altered in Rhot1-deficient β-cells. Islets from diabetic GK rats present Rhot1-deficiency. Finally, RHOT1methylation in blood is associated with future T2D. Together, individuals with T2D exhibit epigenetic alterations linked to mitochondrial dysfunction in pancreatic islets.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

TFEB and TFE3 control glucose homeostasis by regulating insulin gene expression

To fulfill their function, pancreatic beta cells require precise nutrient-sensing mechanisms that control insulin production. Transcription factor EB (TFEB) and its homolog TFE3 have emerged as crucial regulators of the adaptive response of cell metabolism to environmental cues. Here, we show that TFEB and TFE3 regulate beta-cell function and insulin gene expression in response to variations in nutrient availability. We found that nutrient deprivation in beta cells promoted TFEB/TFE3 activation, which resulted in suppression of insulin gene expression. TFEB overexpression was sufficient to inhibit insulin transcription, whereas beta cells depleted of both TFEB and TFE3 failed to suppress insulin gene expression in response to amino acid deprivation. Interestingly, ChIP-seq analysis showed binding of TFEB to super-enhancer regions that regulate insulin transcription. Conditional, beta-cell-specific, Tfeb-overexpressing, and Tfeb/Tfe3 double-KO mice showed severe alteration of insulin transcription, secretion, and glucose tolerance, indicating that TFEB and TFE3 are important physiological mediators of pancreatic function. Our findings reveal a nutrient-controlled transcriptional mechanism that regulates insulin production, thus playing a key role in glucose homeostasis at both cellular and organismal levels.

An Enhancer Within Abcb11 Regulates G6pc2 in C57BL/6 Mouse Pancreatic Islets

G6PC2 is predominantly expressed in pancreatic islet β-cells where it encodes a glucose-6-phosphatase catalytic subunit that modulates the sensitivity of insulin secretion to glucose by opposing the action of glucokinase, thereby regulating fasting blood glucose (FBG). Prior studies have shown that the G6pc2 promoter alone is unable to confer sustained islet-specific gene expression in mice, suggesting the existence of distal enhancers that regulate G6pc2 expression. Using information from both mice and humans and knowledge that single nucleotide polymorphisms (SNPs) both within and near G6PC2 are associated with variations in FBG in humans, we identified several putative enhancers 3' of G6pc2. One region, herein referred to as enhancer I, resides in the 25th intron of Abcb11 and binds multiple islet-enriched transcription factors. CRISPR-mediated deletion of enhancer I in C57BL/6 mice had selective effects on the expression of genes near the G6pc2 locus. In isolated islets, G6pc2 and Spc25 expression were reduced ∼50%, and Gm13613 expression was abolished, whereas Cers6 and nostrin expression were unaffected. This partial reduction in G6pc2 expression enhanced islet insulin secretion at basal glucose concentrations but did not affect FBG or glucose tolerance in vivo, consistent with the absence of a phenotype in G6pc2 heterozygous C57BL/6 mice.

ARTICLE HIGHLIGHTS

Dietary flavonoids prevent diabetes through epigenetic regulation: advance and challenge

The pathophysiology of diabetes has been studied extensively in various countries, but effective prevention and treatment methods are still insufficient. In recent years, epigenetics has received increasing attention from researchers in exploring the etiology and treatment of diabetes. DNA methylation, histone modifications, and non-coding RNAs play critical roles in the occurrence, maintenance, and progression of diabetes and its complications. Therefore, preventing or reversing the epigenetic alterations that occur during the development of diabetes may reduce the individual and societal burden of the disease. Dietary flavonoids serve as natural epigenetic modulators for the discovery of biomarkers for diabetes prevention and the development of alternative therapies. However, there is limited knowledge about the potential beneficial effects of flavonoids on the epigenetics of diabetes. In this review, the multidimensional epigenetic effects of different flavonoid subtypes in diabetes were summarized. Furthermore, it was discussed that parental flavonoid diets might reduce diabetes incidence in offspring, which represent a promising opportunity to prevent diabetes in the future. Future work will depend on exploring anti-diabetic effects of different flavonoids with different epigenetic regulation mechanisms and clinical trials.Highlights• "Epigenetic therapy" could reduce the burden of diabetic patients• "Epigenetic diet" ameliorates diabetes• Targeting epigenetic regulations by dietary flavonoids in the diabetes prevention• Dietary flavonoids prevent diabetes via transgenerational epigenetic inheritance.

Multiple genetic variants at the SLC30A8 locus affect local super-enhancer activity and influence pancreatic β-cell survival and function

Variants at the SLC30A8 locus are associated with type 2 diabetes (T2D) risk. The lead variant, rs13266634, encodes an amino acid change, Arg325Trp (R325W), at the C-terminus of the secretory granule-enriched zinc transporter, ZnT8. Although this protein-coding variant was previously thought to be the sole driver of T2D risk at this locus, recent studies have provided evidence for lowered expression of SLC30A8 mRNA in protective allele carriers. In the present study, combined allele-specific expression (cASE) analysis in human islets revealed multiple variants that influence SLC30A8 expression. Epigenomic mapping identified an islet-selective enhancer cluster at the SLC30A8 locus, hosting multiple T2D risk and cASE associations, which is spatially associated with the SLC30A8 promoter and additional neighbouring genes. Deletions of variant-bearing enhancer regions using CRISPR-Cas9 in human-derived EndoC-βH3 cells lowered the expression of SLC30A8 and several neighbouring genes, and improved insulin secretion. Whilst down-regulation of SLC30A8 had no effect on beta cell survival, loss of UTP23, RAD21 or MED30 markedly reduced cell viability. Although eQTL or cASE analyses in human islets did not support the association between these additional genes and diabetes risk, the transcriptional regulator JQ1 lowered the expression of multiple genes at the SLC30A8 locus and enhanced stimulated insulin secretion.

Regulation of CTCF loop formation during pancreatic cell differentiation

Transcription reprogramming during cell differentiation involves targeting enhancers to genes responsible for establishment of cell fates. To understand the contribution of CTCF-mediated chromatin organization to cell lineage commitment, we analyzed 3D chromatin architecture during the differentiation of human embryonic stem cells into pancreatic islet organoids. We find that CTCF loops are formed and disassembled at different stages of the differentiation process by either recruitment of CTCF to new anchor sites or use of pre-existing sites not previously involved in loop formation. Recruitment of CTCF to new sites in the genome involves demethylation of H3K9me3 to H3K9me2, demethylation of DNA, recruitment of pioneer factors, and positioning of nucleosomes flanking the new CTCF sites. Existing CTCF sites not involved in loop formation become functional loop anchors via the establishment of new cohesin loading sites containing NIPBL and YY1 at sites between the new anchors. In both cases, formation of new CTCF loops leads to strengthening of enhancer promoter interactions and increased transcription of genes adjacent to loop anchors. These results suggest an important role for CTCF and cohesin in controlling gene expression during cell differentiation.

Three-Dimensional Gene Regulation Network in Glioblastoma Ferroptosis

Ferroptosis is an iron-dependent form of cell death, which is reported to be associated with glioma progression and drug sensitivity. Targeting ferroptosis is a potential therapeutic approach for glioma. However, the molecular mechanism of glioma cell ferroptosis is not clear. In this study, we profile the change of 3D chromatin structure in glioblastoma ferroptosis by using HiChIP and study the 3D gene regulation network in glioblastoma ferroptosis. A combination of an analysis of HiChIP and RNA-seq data suggests that change of chromatin loops mediated by 3D chromatin structure regulates gene expressions in glioblastoma ferroptosis. Genes that are regulated by 3D chromatin structures include genes that were reported to function in ferroptosis, like HDM2 and TXNRD1. We propose a new regulatory mechanism governing glioblastoma cell ferroptosis by 3D chromatin structure.

Novel regulators of islet function identified from genetic variation in mouse islet Ca2+ oscillations

Insufficient insulin secretion to meet metabolic demand results in diabetes. The intracellular flux of Ca2+ into β-cells triggers insulin release. Since genetics strongly influences variation in islet secretory responses, we surveyed islet Ca2+ dynamics in eight genetically diverse mouse strains. We found high strain variation in response to four conditions: (1) 8 mM glucose; (2) 8 mM glucose plus amino acids; (3) 8 mM glucose, amino acids, plus 10 nM glucose-dependent insulinotropic polypeptide (GIP); and (4) 2 mM glucose. These stimuli interrogate β-cell function, α- to β-cell signaling, and incretin responses. We then correlated components of the Ca2+ waveforms to islet protein abundances in the same strains used for the Ca2+ measurements. To focus on proteins relevant to human islet function, we identified human orthologues of correlated mouse proteins that are proximal to glycemic-associated single-nucleotide polymorphisms in human genome-wide association studies. Several orthologues have previously been shown to regulate insulin secretion (e.g. ABCC8, PCSK1, and GCK), supporting our mouse-to-human integration as a discovery platform. By integrating these data, we nominate novel regulators of islet Ca2+ oscillations and insulin secretion with potential relevance for human islet function. We also provide a resource for identifying appropriate mouse strains in which to study these regulators.

Rare variant association analysis in 51,256 type 2 diabetes cases and 370,487 controls informs the spectrum of pathogenicity of monogenic diabetes genes

We meta-analyzed array data imputed with the TOPMed reference panel and whole-genome sequence (WGS) datasets and performed the largest, rare variant (minor allele frequency as low as 5×10-5) GWAS meta-analysis of type 2 diabetes (T2D) comprising 51,256 cases and 370,487 controls. We identified 52 novel variants at genome-wide significance (p<5 × 10-8), including 8 novel variants that were either rare or ancestry-specific. Among them, we identified a rare missense variant in HNF4A p.Arg114Trp (OR=8.2, 95% confidence interval [CI]=4.6-14.0, p = 1.08×10-13), previously reported as a variant implicated in Maturity Onset Diabetes of the Young (MODY) with incomplete penetrance. We demonstrated that the diabetes risk in carriers of this variant was modulated by a T2D common variant polygenic risk score (cvPRS) (carriers in the top PRS tertile [OR=18.3, 95%CI=7.2-46.9, p=1.2×10-9] vs carriers in the bottom PRS tertile [OR=2.6, 95% CI=0.97-7.09, p = 0.06]. Association results identified eight variants of intermediate penetrance (OR>5) in monogenic diabetes (MD), which in aggregate as a rare variant PRS were associated with T2D in an independent WGS dataset (OR=4.7, 95% CI=1.86-11.77], p = 0.001). Our data also provided support evidence for 21% of the variants reported in ClinVar in these MD genes as benign based on lack of association with T2D. Our work provides a framework for using rare variant imputation and WGS analyses in large-scale population-based association studies to identify large-effect rare variants and provide evidence for informing variant pathogenicity.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying β cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR/Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR/Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired β cell PCSK1 regulation and insulin secretion. Multiplex CRISPR/Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

Modeling islet enhancers using deep learning identifies candidate causal variants at loci associated with T2D and glycemic traits

Genetic association studies have identified hundreds of independent signals associated with type 2 diabetes (T2D) and related traits. Despite these successes, the identification of specific causal variants underlying a genetic association signal remains challenging. In this study, we describe a deep learning (DL) method to analyze the impact of sequence variants on enhancers. Focusing on pancreatic islets, a T2D relevant tissue, we show that our model learns islet-specific transcription factor (TF) regulatory patterns and can be used to prioritize candidate causal variants. At 101 genetic signals associated with T2D and related glycemic traits where multiple variants occur in linkage disequilibrium, our method nominates a single causal variant for each association signal, including three variants previously shown to alter reporter activity in islet-relevant cell types. For another signal associated with blood glucose levels, we biochemically test all candidate causal variants from statistical fine-mapping using a pancreatic islet beta cell line and show biochemical evidence of allelic effects on TF binding for the model-prioritized variant. To aid in future research, we publicly distribute our model and islet enhancer perturbation scores across ~67 million genetic variants. We anticipate that DL methods like the one presented in this study will enhance the prioritization of candidate causal variants for functional studies.

Pancreatic islet α cell function and proliferation requires the arginine transporter SLC7A2

Interrupting glucagon signaling decreases gluconeogenesis and the fractional extraction of amino acids by liver from blood resulting in lower glycemia. The resulting hyperaminoacidemia stimulates α cell proliferation and glucagon secretion via a liver-α cell axis. We hypothesized that α cells detect and respond to circulating amino acids levels via a unique amino acid transporter repertoire. We found that Slc7a2ISLC7A2 is the most highly expressed cationic amino acid transporter in α cells with its expression being three-fold greater in α than β cells in both mouse and human. Employing cell culture, zebrafish, and knockout mouse models, we found that the cationic amino acid arginine and SLC7A2 are required for α cell proliferation in response to interrupted glucagon signaling. Ex vivo and in vivo assessment of islet function in Slc7a2-/- mice showed decreased arginine-stimulated glucagon and insulin secretion. We found that arginine activation of mTOR signaling and induction of the glutamine transporter SLC38A5 was dependent on SLC7A2, showing that both's role in α cell proliferation is dependent on arginine transport and SLC7A2. Finally, we identified single nucleotide polymorphisms in SLC7A2 associated with HbA1c. Together, these data indicate a central role for SLC7A2 in amino acid-stimulated α cell proliferation and islet hormone secretion.

Inferring causal genes at type 2 diabetes GWAS loci through chromosome interactions in islet cells

Background: Resolving causal genes for type 2 diabetes at loci implicated by genome-wide association studies (GWAS) requires integrating functional genomic data from relevant cell types. Chromatin features in endocrine cells of the pancreatic islet are particularly informative and recent studies leveraging chromosome conformation capture (3C) with Hi-C based methods have elucidated regulatory mechanisms in human islets. However, these genome-wide approaches are less sensitive and afford lower resolution than methods that target specific loci. Methods: To gauge the extent to which targeted 3C further resolves chromatin-mediated regulatory mechanisms at GWAS loci, we generated interaction profiles at 23 loci using next-generation (NG) capture-C in a human beta cell model (EndoC-βH1) and contrasted these maps with Hi-C maps in EndoC-βH1 cells and human islets and a promoter capture Hi-C map in human islets. Results: We found improvements in assay sensitivity of up to 33-fold and resolved ~3.6X more chromatin interactions. At a subset of 18 loci with 25 co-localised GWAS and eQTL signals, NG Capture-C interactions implicated effector transcripts at five additional genetic signals relative to promoter capture Hi-C through physical contact with gene promoters. Conclusions: High resolution chromatin interaction profiles at selectively targeted loci can complement genome- and promoter-wide maps.

Priority index for asthma (PIA): In silico discovery of shared and distinct drug targets for adult- and childhood-onset disease

Asthma is a chronic disease that is caused by a combination of genetic risks and environmental triggers and can affect both adults and children. Genome-wide association studies have revealed partly distinct genetic architectures for its two age-of-onset subtypes (namely, adult-onset and childhood-onset). We reason that identifying shared and distinct drug targets between these subtypes may inform the development of subtype-specific therapeutic strategies. In attempting this, we here introduce Priority Index for Asthma or PIA, a genetics-led and network-driven drug target prioritisation tool for asthma. We demonstrate the validity of the tool in improving drug target prioritisation for asthma compared to the status quo methods, as well as in capturing the underlying etiology and existing therapeutics for the disease. We also illustrate how PIA can be used to prioritise drug targets for adult- and childhood-onset asthma, as well as to identify shared and distinct pathway crosstalk genes. Shared crosstalk genes are mostly involved in JAK-STAT signaling, with clinical evidence supporting that targeting this pathway may be a promising drug repurposing opportunity for both subtypes. Crosstalk genes specific to childhood-onset asthma are enriched for PI3K-AKT-mTOR signaling, and we identify genes that are already targeted by licensed medications as repurposed drug candidates for this subtype. We make all our results accessible and reproducible at http://www.genetictargets.com/PIA. Collectively, our study has significant implications for asthma computational medicine research and can guide the future development of subtype-specific therapeutic strategies for the disease.

RELA governs a network of islet-specific metabolic genes necessary for beta cell function

AIMS/HYPOTHESIS

NF-κB activation unites metabolic and inflammatory responses in many diseases yet less is known about the role that NF-κB plays in normal metabolism. In this study we investigated how RELA impacts the beta cell transcriptional landscape and provides network control over glucoregulation.

METHODS

We generated novel mouse lines harbouring beta cell-specific deletion of either the Rela gene, encoding the canonical NF-κB transcription factor p65 (βp65KO mice), or the Ikbkg gene, encoding the NF-κB essential modulator NEMO (βNEMOKO mice), as well as βA20Tg mice that carry beta cell-specific and forced transgenic expression of the NF-κB-negative regulator gene Tnfaip3, which encodes the A20 protein. Mouse studies were complemented by bioinformatics analysis of human islet chromatin accessibility (assay for transposase-accessible chromatin with sequencing [ATAC-seq]), promoter capture Hi-C (pcHi-C) and p65 binding (chromatin immunoprecipitation-sequencing [ChIP-seq]) data to investigate genome-wide control of the human beta cell metabolic programme.

RESULTS

Rela deficiency resulted in complete loss of stimulus-dependent inflammatory gene upregulation, consistent with its known role in governing inflammation. However, Rela deletion also rendered mice glucose intolerant because of functional loss of insulin secretion. Glucose intolerance was intrinsic to beta cells as βp65KO islets failed to secrete insulin ex vivo in response to a glucose challenge and were unable to restore metabolic control when transplanted into secondary chemical-induced hyperglycaemic recipients. Maintenance of glucose tolerance required Rela but was independent of classical NF-κB inflammatory cascades, as blocking NF-κB signalling in vivo by beta cell knockout of Ikbkg (NEMO), or beta cell overexpression of Tnfaip3 (A20), did not cause severe glucose intolerance. Thus, basal p65 activity has an essential and islet-intrinsic role in maintaining normal glucose homeostasis. Genome-wide bioinformatic mapping revealed the presence of p65 binding sites in the promoter regions of specific metabolic genes and in the majority of islet enhancer hubs (~70% of ~1300 hubs), which are responsible for shaping beta cell type-specific gene expression programmes. Indeed, the islet-specific metabolic genes Slc2a2, Capn9 and Pfkm identified within the large network of islet enhancer hub genes showed dysregulated expression in βp65KO islets.

CONCLUSIONS/INTERPRETATION

These data demonstrate an unappreciated role for RELA as a regulator of islet-specific transcriptional programmes necessary for the maintenance of healthy glucose metabolism. These findings have clinical implications for the use of anti-inflammatories, which influence NF-κB activation and are associated with diabetes.

MethNet: a robust approach to identify regulatory hubs and their distal targets in cancer

Aberrations in the capacity of DNA/chromatin modifiers and transcription factors to bind non-coding regions can lead to changes in gene regulation and impact disease phenotypes. However, identifying distal regulatory elements and connecting them with their target genes remains challenging. Here, we present MethNet, a pipeline that integrates large-scale DNA methylation and gene expression data across multiple cancers, to uncover novel cis regulatory elements (CREs) in a 1Mb region around every promoter in the genome. MethNet identifies clusters of highly ranked CREs, referred to as 'hubs', which contribute to the regulation of multiple genes and significantly affect patient survival. Promoter-capture Hi-C confirmed that highly ranked associations involve physical interactions between CREs and their gene targets, and CRISPRi based scRNA Perturb-seq validated the functional impact of CREs. Thus, MethNet-identified CREs represent a valuable resource for unraveling complex mechanisms underlying gene expression, and for prioritizing the verification of predicted non-coding disease hotspots.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that are irresponsive to topological changes. In each lineage, a high degree of connectivity or "hubness" positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages, compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

Interactive enhancer hubs (iHUBs) mediate transcriptional reprogramming and adaptive resistance in pancreatic cancer

OBJECTIVE

Pancreatic ductal adenocarcinoma (PDAC) displays a remarkable propensity towards therapy resistance. However, molecular epigenetic and transcriptional mechanisms enabling this are poorly understood. In this study, we aimed to identify novel mechanistic approaches to overcome or prevent resistance in PDAC.

DESIGN

We used in vitro and in vivo models of resistant PDAC and integrated epigenomic, transcriptomic, nascent RNA and chromatin topology data. We identified a JunD-driven subgroup of enhancers, called interactive hubs (iHUBs), which mediate transcriptional reprogramming and chemoresistance in PDAC.

RESULTS

iHUBs display characteristics typical for active enhancers (H3K27ac enrichment) in both therapy sensitive and resistant states but exhibit increased interactions and production of enhancer RNA (eRNA) in the resistant state. Notably, deletion of individual iHUBs was sufficient to decrease transcription of target genes and sensitise resistant cells to chemotherapy. Overlapping motif analysis and transcriptional profiling identified the activator protein 1 (AP1) transcription factor JunD as a master transcription factor of these enhancers. JunD depletion decreased iHUB interaction frequency and transcription of target genes. Moreover, targeting either eRNA production or signaling pathways upstream of iHUB activation using clinically tested small molecule inhibitors decreased eRNA production and interaction frequency, and restored chemotherapy responsiveness in vitro and in vivo. Representative iHUB target genes were found to be more expressed in patients with poor response to chemotherapy compared with responsive patients.

CONCLUSION

Our findings identify an important role for a subgroup of highly connected enhancers (iHUBs) in regulating chemotherapy response and demonstrate targetability in sensitisation to chemotherapy.

Integrative genomic analyses in adipocytes implicate DNA methylation in human obesity and diabetes

DNA methylation variations are prevalent in human obesity but evidence of a causative role in disease pathogenesis is limited. Here, we combine epigenome-wide association and integrative genomics to investigate the impact of adipocyte DNA methylation variations in human obesity. We discover extensive DNA methylation changes that are robustly associated with obesity (N = 190 samples, 691 loci in subcutaneous and 173 loci in visceral adipocytes, P < 1 × 10-7). We connect obesity-associated methylation variations to transcriptomic changes at >500 target genes, and identify putative methylation-transcription factor interactions. Through Mendelian Randomisation, we infer causal effects of methylation on obesity and obesity-induced metabolic disturbances at 59 independent loci. Targeted methylation sequencing, CRISPR-activation and gene silencing in adipocytes, further identifies regional methylation variations, underlying regulatory elements and novel cellular metabolic effects. Our results indicate DNA methylation is an important determinant of human obesity and its metabolic complications, and reveal mechanisms through which altered methylation may impact adipocyte functions.

Divergent acute versus prolonged pharmacological GLP-1R responses in adult β cell-specific β-arrestin 2 knockout mice

The glucagon-like peptide-1 receptor (GLP-1R) is a major type 2 diabetes therapeutic target. Stimulated GLP-1Rs are rapidly desensitized by β-arrestins, scaffolding proteins that not only terminate G protein interactions but also act as independent signaling mediators. Here, we have assessed in vivo glycemic responses to the pharmacological GLP-1R agonist exendin-4 in adult β cell-specific β-arrestin 2 knockout (KO) mice. KOs displayed a sex-dimorphic phenotype consisting of weaker acute responses that improved 6 hours after agonist injection. Similar effects were observed for semaglutide and tirzepatide but not with biased agonist exendin-phe1. Acute cyclic adenosine 5'-monophosphate increases were impaired, but desensitization reduced in KO islets. The former defect was attributed to enhanced β-arrestin 1 and phosphodiesterase 4 activities, while reduced desensitization co-occurred with impaired GLP-1R recycling and lysosomal targeting, increased trans-Golgi network signaling, and reduced GLP-1R ubiquitination. This study has unveiled fundamental aspects of GLP-1R response regulation with direct application to the rational design of GLP-1R-targeting therapeutics.

Regulatory de novo mutations underlying intellectual disability

The genetic aetiology of a major fraction of patients with intellectual disability (ID) remains unknown. De novo mutations (DNMs) in protein-coding genes explain up to 40% of cases, but the potential role of regulatory DNMs is still poorly understood. We sequenced 63 whole genomes from 21 ID probands and their unaffected parents. In addition, we analysed 30 previously sequenced genomes from exome-negative ID probands. We found that regulatory DNMs were selectively enriched in fetal brain-specific enhancers as compared with adult brain enhancers. DNM-containing enhancers were associated with genes that show preferential expression in the prefrontal cortex. Furthermore, we identified recurrently mutated enhancer clusters that regulate genes involved in nervous system development (CSMD1, OLFM1, and POU3F3). Most of the DNMs from ID probands showed allele-specific enhancer activity when tested using luciferase assay. Using CRISPR-mediated mutation and editing of epigenomic marks, we show that DNMs at regulatory elements affect the expression of putative target genes. Our results, therefore, provide new evidence to indicate that DNMs in fetal brain-specific enhancers play an essential role in the aetiology of ID.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

Chromatin accessibility differences between alpha, beta, and delta cells identifies common and cell type-specific enhancers

BACKGROUND

High throughput sequencing has enabled the interrogation of the transcriptomic landscape of glucagon-secreting alpha cells, insulin-secreting beta cells, and somatostatin-secreting delta cells. These approaches have furthered our understanding of expression patterns that define healthy or diseased islet cell types and helped explicate some of the intricacies between major islet cell crosstalk and glucose regulation. All three endocrine cell types derive from a common pancreatic progenitor, yet alpha and beta cells have partially opposing functions, and delta cells modulate and control insulin and glucagon release. While gene expression signatures that define and maintain cellular identity have been widely explored, the underlying epigenetic components are incompletely characterized and understood. However, chromatin accessibility and remodeling is a dynamic attribute that plays a critical role to determine and maintain cellular identity.

RESULTS

Here, we compare and contrast the chromatin landscape between mouse alpha, beta, and delta cells using ATAC-Seq to evaluate the significant differences in chromatin accessibility. The similarities and differences in chromatin accessibility between these related islet endocrine cells help define their fate in support of their distinct functional roles. We identify patterns that suggest that both alpha and delta cells are poised, but repressed, from becoming beta-like. We also identify patterns in differentially enriched chromatin that have transcription factor motifs preferentially associated with different regions of the genome. Finally, we not only confirm and visualize previously discovered common endocrine- and cell specific- enhancer regions across differentially enriched chromatin, but identify novel regions as well. We compiled our chromatin accessibility data in a freely accessible database of common endocrine- and cell specific-enhancer regions that can be navigated with minimal bioinformatics expertise.

CONCLUSIONS

Both alpha and delta cells appear poised, but repressed, from becoming beta cells in murine pancreatic islets. These data broadly support earlier findings on the plasticity in identity of non-beta cells under certain circumstances. Furthermore, differential chromatin accessibility shows preferentially enriched distal-intergenic regions in beta cells, when compared to either alpha or delta cells.

Nutrient regulation of the islet epigenome controls adaptive insulin secretion

Adaptation of the islet β cell insulin-secretory response to changing insulin demand is critical for blood glucose homeostasis, yet the mechanisms underlying this adaptation are unknown. Here, we have shown that nutrient-stimulated histone acetylation plays a key role in adapting insulin secretion through regulation of genes involved in β cell nutrient sensing and metabolism. Nutrient regulation of the epigenome occurred at sites occupied by the chromatin-modifying enzyme lysine-specific demethylase 1 (Lsd1) in islets. β Cell-specific deletion of Lsd1 led to insulin hypersecretion, aberrant expression of nutrient-response genes, and histone hyperacetylation. Islets from mice adapted to chronically increased insulin demand exhibited shared epigenetic and transcriptional changes. Moreover, we found that genetic variants associated with type 2 diabetes were enriched at LSD1-bound sites in human islets, suggesting that interpretation of nutrient signals is genetically determined and clinically relevant. Overall, these studies revealed that adaptive insulin secretion involves Lsd1-mediated coupling of nutrient state to regulation of the islet epigenome.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.

Contributions of 3D chromatin structure to cell-type-specific gene regulation

Eukaryotic genomes are organized in 3D in a multiscale manner, and different mechanisms acting at each of these scales can contribute to transcriptional regulation. However, the large single-cell variability in 3D chromatin structures represents a challenge to understand how transcription may be differentially regulated between cell types in a robust and efficient manner. Here, we describe the different mechanisms by which 3D chromatin structure was shown to contribute to cell-type-specific transcriptional regulation. Excitingly, several novel methodologies able to measure 3D chromatin conformation and transcription in single cells in their native tissue context, or to detect the dynamics of cis-regulatory interactions, are starting to allow quantitative dissection of chromatin structure noise and relate it to how transcription may be regulated between different cell types and cell states.

Methylation of histone H3 lysine 4 is required for maintenance of beta cell function in adult mice

AIMS/HYPOTHESIS

Beta cells control glucose homeostasis via regulated production and secretion of insulin. This function arises from a highly specialised gene expression programme that is established during development and then sustained, with limited flexibility, in terminally differentiated cells. Dysregulation of this programme is seen in type 2 diabetes but mechanisms that preserve gene expression or underlie its dysregulation in mature cells are not well resolved. This study investigated whether methylation of histone H3 lysine 4 (H3K4), a marker of gene promoters with unresolved functional importance, is necessary for the maintenance of mature beta cell function.

METHODS

Beta cell function, gene expression and chromatin modifications were analysed in conditional Dpy30 knockout mice, in which H3K4 methyltransferase activity is impaired, and in a mouse model of diabetes.

RESULTS

H3K4 methylation maintains expression of genes that are important for insulin biosynthesis and glucose responsiveness. Deficient methylation of H3K4 leads to a less active and more repressed epigenome profile that locally correlates with gene expression deficits but does not globally reduce gene expression. Instead, developmentally regulated genes and genes in weakly active or suppressed states particularly rely on H3K4 methylation. We further show that H3K4 trimethylation (H3K4me3) is reorganised in islets from the Lepr^db/db mouse model of diabetes in favour of weakly active and disallowed genes at the expense of terminal beta cell markers with broad H3K4me3 peaks.

CONCLUSIONS/INTERPRETATION

Sustained methylation of H3K4 is critical for the maintenance of beta cell function. Redistribution of H3K4me3 is linked to gene expression changes that are implicated in diabetes pathology.

Interaction-integrated linear mixed model reveals 3D-genetic basis underlying Autism

Genetic interactions play critical roles in genotype-phenotype associations. We developed a novel interaction-integrated linear mixed model (ILMM) that integrates a priori knowledge into linear mixed models. ILMM enables statistical integration of genetic interactions upfront and overcomes the problems of searching for combinations. To demonstrate its utility, with 3D genomic interactions (assessed by Hi-C experiments) as a priori, we applied ILMM to whole-genome sequencing data for Autism Spectrum Disorders (ASD) and brain transcriptome data, revealing the 3D-genetic basis of ASD and 3D-expression quantitative loci (3D-eQTLs) for brain tissues. Notably, we reported a potential mechanism involving distal regulation between FOXP2 and DNMT3A, conferring the risk of ASD.

Multiscale 3D genome reorganization during skeletal muscle stem cell lineage progression and aging

Little is known about three-dimensional (3D) genome organization in skeletal muscle stem cells [also called satellite cells (SCs)]. Here, we comprehensively map the 3D genome topology reorganization during mouse SC lineage progression. Specifically, rewiring at the compartment level is most pronounced when SCs become activated. Marked loss in topologically associating domain (TAD) border insulation and chromatin looping also occurs during early activation process. Meanwhile, TADs can form TAD clusters and super-enhancer-containing TAD clusters orchestrate stage-specific gene expression. Furthermore, we uncover that transcription factor PAX7 is pivotal in enhancer-promoter (E-P) loop formation. We also identify cis-regulatory elements that are crucial for local chromatin organization at Pax7 locus and Pax7 expression. Lastly, we unveil that geriatric SC displays a prominent gain in long-range contacts and loss of TAD border insulation. Together, our results uncover that 3D chromatin extensively reorganizes at multiple architectural levels and underpins the transcriptome remodeling during SC lineage development and SC aging.