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

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.

Loci for insulin processing and secretion provide insight into type 2 diabetes risk

Insulin secretion is critical for glucose homeostasis, and increased levels of the precursor proinsulin relative to insulin indicate pancreatic islet beta-cell stress and insufficient insulin secretory capacity in the setting of insulin resistance. We conducted meta-analyses of genome-wide association results for fasting proinsulin from 16 European-ancestry studies in 45,861 individuals. We found 36 independent signals at 30 loci (p value < 5 × 10-8), which validated 12 previously reported loci for proinsulin and ten additional loci previously identified for another glycemic trait. Half of the alleles associated with higher proinsulin showed higher rather than lower effects on glucose levels, corresponding to different mechanisms. Proinsulin loci included genes that affect prohormone convertases, beta-cell dysfunction, vesicle trafficking, beta-cell transcriptional regulation, and lysosomes/autophagy processes. We colocalized 11 proinsulin signals with islet expression quantitative trait locus (eQTL) data, suggesting candidate genes, including ARSG, WIPI1, SLC7A14, and SIX3. The NKX6-3/ANK1 proinsulin signal colocalized with a T2D signal and an adipose ANK1 eQTL signal but not the islet NKX6-3 eQTL. Signals were enriched for islet enhancers, and we showed a plausible islet regulatory mechanism for the lead signal in the MADD locus. These results show how detailed genetic studies of an intermediate phenotype can elucidate mechanisms that may predispose one to disease.

Pancreatic microexons regulate islet function and glucose homeostasis

Pancreatic islets control glucose homeostasis by the balanced secretion of insulin and other hormones, and their abnormal function causes diabetes or hypoglycaemia. Here we uncover a conserved programme of alternative microexons included in mRNAs of islet cells, particularly in genes involved in vesicle transport and exocytosis. Islet microexons (IsletMICs) are regulated by the RNA binding protein SRRM3 and represent a subset of the larger neural programme that are particularly sensitive to SRRM3 levels. Both SRRM3 and IsletMICs are induced by elevated glucose levels, and depletion of SRRM3 in human and rat beta cell lines and mouse islets, or repression of particular IsletMICs using antisense oligonucleotides, leads to inappropriate insulin secretion. Consistently, mice harbouring mutations in Srrm3 display defects in islet cell identity and function, leading to hyperinsulinaemic hypoglycaemia. Importantly, human genetic variants that influence SRRM3 expression and IsletMIC inclusion in islets are associated with fasting glucose variation and type 2 diabetes risk. Taken together, our data identify a conserved microexon programme that regulates glucose homeostasis.

3D genome alterations in T cells associated with disease activity of systemic lupus erythematosus

OBJECTIVES

Three-dimensional (3D) genome alterations can dysregulate gene expression by rewiring physical interactions within chromosomes in a tissue-specific or cell-specific manner and lead to diseases. We aimed to elucidate the 3D genome structure and its role in gene expression networks dysregulated in systemic lupus erythematosus (SLE).

METHODS

We performed Hi-C experiments using CD4+ T cells from 7 patients with SLE and 5 age-matched and sex-matched healthy controls (HCs) combined with RNA sequencing analysis. Further integrative analyses, including transcription factor motif enrichment, SPI1 knockdown and histone modifications (H3K27ac, H3K4me1, H3K4me3), were performed for altered loop-associated gene loci in SLE.

RESULTS

We deciphered the 3D chromosome organisation in T cells of patients with SLE and found it was clearly distinct from that of HCs and closely associated with the disease activity of SLE. Importantly, we identified loops within chromosomes associated with the disease activity of SLE and differentially expressed genes and found some key histone modifications close to these loops. Moreover, we demonstrated the contribution of the transcription factor SPI1, whose motif is located in the altered loop in SLE, to the overexpression of interferon pathway gene. In addition, we identified the potential influences of genetic variations in 3D genome alterations in SLE.

CONCLUSIONS

Our results highlight the 3D genome structure alterations associated with SLE development and provide a foundation for future interrogation of the relationships between chromosome structure and gene expression control in SLE.

Low input capture Hi-C (liCHi-C) identifies promoter-enhancer interactions at high-resolution

Long-range interactions between regulatory elements and promoters are key in gene transcriptional control; however, their study requires large amounts of starting material, which is not compatible with clinical scenarios nor the study of rare cell populations. Here we introduce low input capture Hi-C (liCHi-C) as a cost-effective, flexible method to map and robustly compare promoter interactomes at high resolution. As proof of its broad applicability, we implement liCHi-C to study normal and malignant human hematopoietic hierarchy in clinical samples. We demonstrate that the dynamic promoter architecture identifies developmental trajectories and orchestrates transcriptional transitions during cell-state commitment. Moreover, liCHi-C enables the identification of disease-relevant cell types, genes and pathways potentially deregulated by non-coding alterations at distal regulatory elements. Finally, we show that liCHi-C can be harnessed to uncover genome-wide structural variants, resolve their breakpoints and infer their pathogenic effects. Collectively, our optimized liCHi-C method expands the study of 3D chromatin organization to unique, low-abundance cell populations, and offers an opportunity to uncover factors and regulatory networks involved in disease pathogenesis.

Insights from multi-omics integration in complex disease primary tissues

Genome-wide association studies (GWAS) have provided insights into the genetic basis of complex diseases. In the next step, integrative multi-omics approaches can characterize molecular profiles in relevant primary tissues to reveal the mechanisms that underlie disease development. Here, we highlight recent progress in four examples of complex diseases generated by integrative studies: type 2 diabetes (T2D), osteoarthritis, Alzheimer's disease (AD), and systemic lupus erythematosus (SLE). High-resolution methodologies such as single-cell and spatial omics techniques will become even more important in the future. Furthermore, we emphasize the urgent need to include as yet understudied cell types and increase the diversity of studied populations.

A genome-wide CRISPR screen identifies CALCOCO2 as a regulator of beta cell function influencing type 2 diabetes risk

Identification of the genes and processes mediating genetic association signals for complex diseases represents a major challenge. As many of the genetic signals for type 2 diabetes (T2D) exert their effects through pancreatic islet-cell dysfunction, we performed a genome-wide pooled CRISPR loss-of-function screen in a human pancreatic beta cell line. We assessed the regulation of insulin content as a disease-relevant readout of beta cell function and identified 580 genes influencing this phenotype. Integration with genetic and genomic data provided experimental support for 20 candidate T2D effector transcripts including the autophagy receptor CALCOCO2. Loss of CALCOCO2 was associated with distorted mitochondria, less proinsulin-containing immature granules and accumulation of autophagosomes upon inhibition of late-stage autophagy. Carriers of T2D-associated variants at the CALCOCO2 locus further displayed altered insulin secretion. Our study highlights how cellular screens can augment existing multi-omic efforts to support mechanistic understanding and provide evidence for causal effects at genome-wide association studies loci.

Novel biological insights revealed from the investigation of multiscale genome architecture

Gene expression and cell fate determination require precise and coordinated epigenetic regulation. The complex three-dimensional (3D) genome organization plays a critical role in transcription in myriad biological processes. A wide range of architectural features of the 3D genome, including chromatin loops, topologically associated domains (TADs), chromatin compartments, and phase separation, together regulate the chromatin state and transcriptional activity at multiple levels. With the help of 3D genome informatics, recent biochemistry and imaging approaches based on different strategies have revealed functional interactions among biomacromolecules, even at the single-cell level. Here, we review the occurrence, mechanistic basis, and functional implications of dynamic genome organization, and outline recent experimental and computational approaches for profiling multiscale genome architecture to provide robust tools for studying the 3D genome.

Using single-cell multi-omics screening of human fetal pancreas to identify novel players in human beta cell development

Islet transplantation from organ donors can considerably improve glucose homeostasis and well-being in individuals with type 1 diabetes, where the beta cells are destroyed by the autoimmune attack, but there are insufficient donor islets to make this a widespread therapy. Strategies are therefore being developed to generate unlimited amounts of insulin-producing beta cells from pluripotent stem cells, with the aim that they will be transplanted to treat diabetes. Whilst much progress has been made in recent years in the directed differentiation of pluripotent stem cells to beta-like cells, essential gaps still exist in generating stem cell-derived beta cells that are fully functional in vitro. This short review provides details of recent multi-'omics' studies of the human fetal pancreas, which are revealing granular information on the various cell types in the developing pancreas. It is anticipated that this fine mapping of the pancreatic cells at single-cell resolution will provide additional insights that can be utilised to reproducibly produce human beta cells in vitro that have the functional characteristics of beta cells within native human islets.

The HASTER lncRNA promoter is a cis-acting transcriptional stabilizer of HNF1A

The biological purpose of long non-coding RNAs (lncRNAs) is poorly understood. Haploinsufficient mutations in HNF1A homeobox A (HNF1A), encoding a homeodomain transcription factor, cause diabetes mellitus. Here, we examine HASTER, the promoter of an lncRNA antisense to HNF1A. Using mouse and human models, we show that HASTER maintains cell-specific physiological HNF1A concentrations through positive and negative feedback loops. Pancreatic β cells from Haster mutant mice consequently showed variegated HNF1A silencing or overexpression, resulting in hyperglycaemia. HASTER-dependent negative feedback was essential to prevent HNF1A binding to inappropriate genomic regions. We demonstrate that the HASTER promoter DNA, rather than the lncRNA, modulates HNF1A promoter-enhancer interactions in cis and thereby regulates HNF1A transcription. Our studies expose a cis-regulatory element that is unlike classic enhancers or silencers, it stabilizes the transcription of its target gene and ensures the fidelity of a cell-specific transcription factor program. They also show that disruption of a mammalian lncRNA promoter can cause diabetes mellitus.

How dysregulation of the immune system promotes diabetes mellitus and cardiovascular risk complications

Diabetes mellitus (DM) is a chronic metabolic disorder characterized by persistent hyperglycemia due to insulin resistance or failure to produce insulin. Patients with DM develop microvascular complications that include chronic kidney disease and retinopathy, and macrovascular complications that mainly consist in an accelerated and more severe atherosclerosis compared to the general population, increasing the risk of cardiovascular (CV) events, such as stroke or myocardial infarction by 2- to 4-fold. DM is commonly associated with a low-grade chronic inflammation that is a known causal factor in its development and its complications. Moreover, it is now well-established that inflammation and immune cells play a major role in both atherosclerosis genesis and progression, as well as in CV event occurrence. In this review, after a brief presentation of DM physiopathology and its macrovascular complications, we will describe the immune system dysregulation present in patients with type 1 or type 2 diabetes and discuss its role in DM cardiovascular complications development. More specifically, we will review the metabolic changes and aberrant activation that occur in the immune cells driving the chronic inflammation through cytokine and chemokine secretion, thus promoting atherosclerosis onset and progression in a DM context. Finally, we will discuss how genetics and recent systemic approaches bring new insights into the mechanisms behind these inflammatory dysregulations and pave the way toward precision medicine.

Genetic regulation of RNA splicing in human pancreatic islets

BACKGROUND

Non-coding genetic variants that influence gene transcription in pancreatic islets play a major role in the susceptibility to type 2 diabetes (T2D), and likely also contribute to type 1 diabetes (T1D) risk. For many loci, however, the mechanisms through which non-coding variants influence diabetes susceptibility are unknown.

RESULTS

We examine splicing QTLs (sQTLs) in pancreatic islets from 399 human donors and observe that common genetic variation has a widespread influence on the splicing of genes with established roles in islet biology and diabetes. In parallel, we profile expression QTLs (eQTLs) and use transcriptome-wide association as well as genetic co-localization studies to assign islet sQTLs or eQTLs to T2D and T1D susceptibility signals, many of which lack candidate effector genes. This analysis reveals biologically plausible mechanisms, including the association of T2D with an sQTL that creates a nonsense isoform in ERO1B, a regulator of ER-stress and proinsulin biosynthesis. The expanded list of T2D risk effector genes reveals overrepresented pathways, including regulators of G-protein-mediated cAMP production. The analysis of sQTLs also reveals candidate effector genes for T1D susceptibility such as DCLRE1B, a senescence regulator, and lncRNA MEG3.

CONCLUSIONS

These data expose widespread effects of common genetic variants on RNA splicing in pancreatic islets. The results support a role for splicing variation in diabetes susceptibility, and offer a new set of genetic targets with potential therapeutic benefit.

Multiple parameters shape the 3D chromatin structure of single nuclei at the doc locus in Drosophila

The spatial organization of chromatin at the scale of topologically associating domains (TADs) and below displays large cell-to-cell variations. Up until now, how this heterogeneity in chromatin conformation is shaped by chromatin condensation, TAD insulation, and transcription has remained mostly elusive. Here, we used Hi-M, a multiplexed DNA-FISH imaging technique providing developmental timing and transcriptional status, to show that the emergence of TADs at the ensemble level partially segregates the conformational space explored by single nuclei during the early development of Drosophila embryos. Surprisingly, a substantial fraction of nuclei display strong insulation even before TADs emerge. Moreover, active transcription within a TAD leads to minor changes to the local inter- and intra-TAD chromatin conformation in single nuclei and only weakly affects insulation to the neighboring TAD. Overall, our results indicate that multiple parameters contribute to shaping the chromatin architecture of single nuclei at the TAD scale.

3D chromatin maps of the human pancreas reveal lineage-specific regulatory architecture of T2D risk

Three-dimensional (3D) chromatin organization maps help dissect cell-type-specific gene regulatory programs. Furthermore, 3D chromatin maps contribute to elucidating the pathogenesis of complex genetic diseases by connecting distal regulatory regions and genetic risk variants to their respective target genes. To understand the cell-type-specific regulatory architecture of diabetes risk, we generated transcriptomic and 3D epigenomic profiles of human pancreatic acinar, alpha, and beta cells using single-cell RNA-seq, single-cell ATAC-seq, and high-resolution Hi-C of sorted cells. Comparisons of these profiles revealed differential A/B (open/closed) chromatin compartmentalization, chromatin looping, and transcriptional factor-mediated control of cell-type-specific gene regulatory programs. We identified a total of 4,750 putative causal-variant-to-target-gene pairs at 194 type 2 diabetes GWAS signals using pancreatic 3D chromatin maps. We found that the connections between candidate causal variants and their putative target effector genes are cell-type stratified and emphasize previously underappreciated roles for alpha and acinar cells in diabetes pathogenesis.

Pancreas agenesis mutations disrupt a lead enhancer controlling a developmental enhancer cluster

Sequence variants in cis-acting enhancers are important for polygenic disease, but their role in Mendelian disease is poorly understood. Redundancy between enhancers that regulate the same gene is thought to mitigate the pathogenic impact of enhancer mutations. Recent findings, however, have shown that loss-of-function mutations in a single enhancer near PTF1A cause pancreas agenesis and neonatal diabetes. Using mouse and human genetic models, we show that this enhancer activates an entire PTF1A enhancer cluster in early pancreatic multipotent progenitors. This leading role, therefore, precludes functional redundancy. We further demonstrate that transient expression of PTF1A in multipotent progenitors sets in motion an epigenetic cascade that is required for duct and endocrine differentiation. These findings shed insights into the genome regulatory mechanisms that drive pancreas differentiation. Furthermore, they reveal an enhancer that acts as a regulatory master key and is thus vulnerable to pathogenic loss-of-function mutations.

•

The pancreas agenesis enhancer (EnhP) activates PTF1A in early pancreatic progenitors

•

EnhP also activates other progenitor PTF1A enhancers

•

This master key function explains why EnhP is vulnerable to loss-of-function mutations

•

Transient PTF1A expression in progenitors controls pancreas growth and endocrinogenesis

Thymocyte regulatory variant alters transcription factor binding and protects from type 1 diabetes in infants

We recently mapped a genetic susceptibility locus on chromosome 6q22.33 for type 1 diabetes (T1D) diagnosed below the age of 7 years between the PTPRK and thymocyte-selection-associated (THEMIS) genes. As the thymus plays a central role in shaping the T cell repertoire, we aimed to identify the most likely causal genetic factors behind this association using thymocyte genomic data. In four thymocyte populations, we identified 253 DNA sequence motifs underlying histone modifications. The G insertion allele of rs138300818, associated with protection from diabetes, created thymocyte motifs for multiple histone modifications and thymocyte types. In a parallel approach to identifying variants that alter transcription factor binding motifs, the same variant disrupted a predicted motif for Rfx7, which is abundantly expressed in the thymus. Chromatin state and RNA sequencing data suggested strong transcription overlapping rs138300818 in fetal thymus, while expression quantitative trait locus and chromatin conformation data associate the insertion with lower THEMIS expression. Extending the analysis to other T1D loci further highlighted rs66733041 affecting the GATA3 transcription factor binding in the AFF3 locus. Taken together, our results support a role for thymic THEMIS gene expression and the rs138300818 variant in promoting the development of early-onset T1D.

TCF7L2 plays a complex role in human adipose progenitor biology, which might contribute to genetic susceptibility to type 2 diabetes

INTRODUCTION

Non-coding genetic variation at TCF7L2 is the strongest genetic determinant of type 2 diabetes (T2D) risk in humans. TCF7L2 encodes a transcription factor mediating the nuclear effects of WNT signaling in adipose tissue (AT). In vivo studies in transgenic mice have highlighted important roles for TCF7L2 in adipose tissue biology and systemic metabolism.

OBJECTIVE

To map the expression of TCF7L2 in human AT, examine its role in human adipose cell biology in vitro, and investigate the effects of the fine-mapped T2D-risk allele at rs7903146 on AT morphology and TCF7L2 expression.

METHODS

Ex vivo gene expression studies of TCF7L2 in whole and fractionated human AT. In vitro TCF7L2 gain- and/or loss-of-function studies in primary and immortalized human adipose progenitor cells (APCs) and mature adipocytes (mADs). AT phenotyping of rs7903146 T2D-risk variant carriers and matched controls.

RESULTS

Adipose progenitors (APs) exhibited the highest TCF7L2 mRNA abundance compared to mature adipocytes and adipose-derived endothelial cells. Obesity was associated with reduced TCF7L2 transcript levels in whole subcutaneous abdominal AT but paradoxically increased expression in APs. In functional studies, TCF7L2 knockdown (KD) in abdominal APs led to dose-dependent activation of WNT/β-catenin signaling, impaired proliferation and dose-dependent effects on adipogenesis. Whilst partial KD enhanced adipocyte differentiation, near-total KD impaired lipid accumulation and adipogenic gene expression. Over-expression of TCF7L2 accelerated adipogenesis. In contrast, TCF7L2-KD in gluteal APs dose-dependently enhanced lipid accumulation. Transcriptome-wide profiling revealed that TCF7L2 might modulate multiple aspects of AP biology including extracellular matrix secretion, immune signaling and apoptosis. The T2D-risk allele at rs7903146 was associated with reduced AP TCF7L2 expression and enhanced AT insulin sensitivity.

CONCLUSIONS

TCF7L2 plays a complex role in AP biology and has both dose- and depot-dependent effects on adipogenesis. In addition to regulating pancreatic insulin secretion, genetic variation at TCF7L2 might also influence T2D risk by modulating AP function.

Recruitment of the CoREST transcription repressor complexes by Nerve Growth factor IB-like receptor (Nurr1/NR4A2) mediates silencing of HIV in microglial cells

Human immune deficiency virus (HIV) infection in the brain leads to chronic neuroinflammation due to the production of pro-inflammatory cytokines, which in turn promotes HIV transcription in infected microglial cells. However, powerful counteracting silencing mechanisms in microglial cells result in the rapid shutdown of HIV expression after viral reactivation to limit neuronal damage. Here we investigated whether the Nerve Growth Factor IB-like nuclear receptor Nurr1 (NR4A2), which is a repressor of inflammation in the brain, acts directly to restrict HIV expression. HIV silencing following activation by TNF-α, or a variety of toll-like receptor (TLR) agonists, in both immortalized human microglial cells (hμglia) and induced pluripotent stem cells (iPSC)-derived human microglial cells (iMG) was enhanced by Nurr1 agonists. Similarly, overexpression of Nurr1 led to viral suppression, while conversely, knock down (KD) of endogenous Nurr1 blocked HIV silencing. The effect of Nurr1 on HIV silencing is direct: Nurr1 binds directly to the specific consensus binding sites in the U3 region of the HIV LTR and mutation of the Nurr1 DNA binding domain blocked its ability to suppress HIV-1 transcription. Chromatin immunoprecipitation (ChIP) assays also showed that after Nurr1 binding to the LTR, the CoREST/HDAC1/G9a/EZH2 transcription repressor complex is recruited to the HIV provirus. Finally, transcriptomic studies demonstrated that in addition to repressing HIV transcription, Nurr1 also downregulated numerous cellular genes involved in inflammation, cell cycle, and metabolism, further promoting HIV latency and microglial homoeostasis. Nurr1 therefore plays a pivotal role in modulating the cycles of proviral reactivation by potentiating the subsequent proviral transcriptional shutdown. These data highlight the therapeutic potential of Nurr1 agonists for inducing HIV silencing and microglial homeostasis and ultimately for the amelioration of the neuroinflammation associated with HIV-associated neurocognitive disorders (HAND).

β Cell-specific deletion of Zfp148 improves nutrient-stimulated β cell Ca2+ responses

Insulin secretion from pancreatic β cells is essential for glucose homeostasis. An insufficient response to the demand for insulin results in diabetes. We previously showed that β cell-specific deletion of Zfp148 (β-Zfp148KO) improves glucose tolerance and insulin secretion in mice. Here, we performed Ca2+ imaging of islets from β‑Zfp148KO and control mice fed both a chow and a Western-style diet. β-Zfp148KO islets demonstrated improved sensitivity and sustained Ca2+ oscillations in response to elevated glucose levels. β-Zfp148KO islets also exhibited elevated sensitivity to amino acid-induced Ca2+ influx under low glucose conditions, suggesting enhanced mitochondrial phosphoenolpyruvate-dependent (PEP-dependent), ATP-sensitive K+ channel closure, independent of glycolysis. RNA-Seq and proteomics of β-Zfp148KO islets revealed altered levels of enzymes involved in amino acid metabolism (specifically, SLC3A2, SLC7A8, GLS, GLS2, PSPH, PHGDH, and PSAT1) and intermediary metabolism (namely, GOT1 and PCK2), consistent with altered PEP cycling. In agreement with this, β-Zfp148KO islets displayed enhanced insulin secretion in response to l-glutamine and activation of glutamate dehydrogenase. Understanding pathways controlled by ZFP148 may provide promising strategies for improving β cell function that are robust to the metabolic challenge imposed by a Western diet.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

Multidimensional chromatin profiling of zebrafish pancreas to uncover and investigate disease-relevant enhancers

The pancreas is a central organ for human diseases. Most alleles uncovered by genome-wide association studies of pancreatic dysfunction traits overlap with non-coding sequences of DNA. Many contain epigenetic marks of cis-regulatory elements active in pancreatic cells, suggesting that alterations in these sequences contribute to pancreatic diseases. Animal models greatly help to understand the role of non-coding alterations in disease. However, interspecies identification of equivalent cis-regulatory elements faces fundamental challenges, including lack of sequence conservation. Here we combine epigenetic assays with reporter assays in zebrafish and human pancreatic cells to identify interspecies functionally equivalent cis-regulatory elements, regardless of sequence conservation. Among other potential disease-relevant enhancers, we identify a zebrafish ptf1a distal-enhancer whose deletion causes pancreatic agenesis, a phenotype previously found to be induced by mutations in a distal-enhancer of PTF1A in humans, further supporting the causality of this condition in vivo. This approach helps to uncover interspecies functionally equivalent cis-regulatory elements and their potential role in human disease.

Alterations in cis-regulatory elements (CREs) can contribute to pancreatic diseases. Here the authors combine chromatin profiling and interaction points with in vivo reporter assays in zebrafish to uncover functionally equivalent human CREs, helping to predict disease-relevant enhancers.

Lupus enhancer risk variant causes dysregulation of IRF8 through cooperative lncRNA and DNA methylation machinery

Despite strong evidence that human genetic variants affect the expression of many key transcription factors involved in autoimmune diseases, establishing biological links between non-coding risk variants and the gene targets they regulate remains a considerable challenge. Here, we combine genetic, epigenomic, and CRISPR activation approaches to screen for functional variants that regulate IRF8 expression. We demonstrate that the locus containing rs2280381 is a cell-type-specific enhancer for IRF8 that spatially interacts with the IRF8 promoter. Further, rs2280381 mediates IRF8 expression through enhancer RNA AC092723.1, which recruits TET1 to the IRF8 promoter regulating IRF8 expression by affecting methylation levels. The alleles of rs2280381 modulate PU.1 binding and chromatin state to regulate AC092723.1 and IRF8 expression differentially. Our work illustrates an integrative strategy to define functional genetic variants that regulate the expression of critical genes in autoimmune diseases and decipher the mechanisms underlying the dysregulation of IRF8 expression mediated by lupus risk variants.

Epigenome rewiring in human pluripotent stem cells

The epigenome plays a crucial role in modulating the activity of regulatory elements, thereby orchestrating diverse transcriptional programs during embryonic development. Human (h)PSC stepwise differentiation provides an excellent platform for capturing dynamic epigenomic events during lineage transition in human development. Here we discuss how recent technological advances, from epigenomic mapping to targeted perturbation, are providing a more comprehensive appreciation of remodeling of the chromatin landscape during human development with implications for aberrant rewiring in disease. We predict that the continuous innovation of hPSC differentiation methods, epigenome mapping, and CRISPR (clustered regularly interspaced short palindromic repeats) perturbation technologies will allow researchers to build toward not only a comprehensive understanding of the epigenomic mechanisms governing development, but also a highly flexible way to model diseases with opportunities for translation.

Incorporating regulatory interactions into gene-set analyses for GWAS data: A controlled analysis with the MAGMA tool

To date, genome-wide association studies have identified thousands of statistically-significant associations between genetic variants, and phenotypes related to a myriad of traits and diseases. A key goal for human-genetics research is to translate these associations into functional mechanisms. Popular gene-set analysis tools, like MAGMA, map variants to genes they might affect, and then integrate genome-wide association study data (that is, variant-level associations for a phenotype) to score genes for association with a phenotype. Gene scores are subsequently used in competitive gene-set analyses to identify biological processes that are enriched for phenotype association. By default, variants are mapped to genes in their proximity. However, many variants that affect phenotypes are thought to act at regulatory elements, which can be hundreds of kilobases away from their target genes. Thus, we explored the idea of augmenting a proximity-based mapping scheme with publicly-available datasets of regulatory interactions. We used MAGMA to analyze genome-wide association study data for ten different phenotypes, and evaluated the effects of augmentation by comparing numbers, and identities, of genes and gene sets detected as statistically significant between mappings. We detected several pitfalls and confounders of such "augmented analyses", and introduced ways to control for them. Using these controls, we demonstrated that augmentation with datasets of regulatory interactions only occasionally strengthened the enrichment for phenotype association amongst (biologically-relevant) gene sets for different phenotypes. Still, in such cases, genes and regulatory elements responsible for the improvement could be pinpointed. For instance, using brain regulatory-interactions for augmentation, we were able to implicate two acetylcholine receptor subunits involved in post-synaptic chemical transmission, namely CHRNB2 and CHRNE, in schizophrenia. Collectively, our study presents a critical approach for integrating regulatory interactions into gene-set analyses for genome-wide association study data, by introducing various controls to distinguish genuine results from spurious discoveries.

Multi-ethnic GWAS and fine-mapping of glycaemic traits identify novel loci in the PAGE Study

AIMS/HYPOTHESIS

Type 2 diabetes is a growing global public health challenge. Investigating quantitative traits, including fasting glucose, fasting insulin and HbA1c, that serve as early markers of type 2 diabetes progression may lead to a deeper understanding of the genetic aetiology of type 2 diabetes development. Previous genome-wide association studies (GWAS) have identified over 500 loci associated with type 2 diabetes, glycaemic traits and insulin-related traits. However, most of these findings were based only on populations of European ancestry. To address this research gap, we examined the genetic basis of fasting glucose, fasting insulin and HbA1c in participants of the diverse Population Architecture using Genomics and Epidemiology (PAGE) Study.

METHODS

We conducted a GWAS of fasting glucose (n = 52,267), fasting insulin (n = 48,395) and HbA1c (n = 23,357) in participants without diabetes from the diverse PAGE Study (23% self-reported African American, 46% Hispanic/Latino, 40% European, 4% Asian, 3% Native Hawaiian, 0.8% Native American), performing transethnic and population-specific GWAS meta-analyses, followed by fine-mapping to identify and characterise novel loci and independent secondary signals in known loci.

RESULTS

Four novel associations were identified (p < 5 × 10-9), including three loci associated with fasting insulin, and a novel, low-frequency African American-specific locus associated with fasting glucose. Additionally, seven secondary signals were identified, including novel independent secondary signals for fasting glucose at the known GCK locus and for fasting insulin at the known PPP1R3B locus in transethnic meta-analysis.

CONCLUSIONS/INTERPRETATION

Our findings provide new insights into the genetic architecture of glycaemic traits and highlight the continued importance of conducting genetic studies in diverse populations.

DATA AVAILABILITY

Full summary statistics from each of the population-specific and transethnic results are available at NHGRI-EBI GWAS catalog ( https://www.ebi.ac.uk/gwas/downloads/summary-statistics ).

An effector index to predict target genes at GWAS loci

Drug development and biological discovery require effective strategies to map existing genetic associations to causal genes. To approach this problem, we selected 12 common diseases and quantitative traits for which highly powered genome-wide association studies (GWAS) were available. For each disease or trait, we systematically curated positive control gene sets from Mendelian forms of the disease and from targets of medicines used for disease treatment. We found that these positive control genes were highly enriched in proximity of GWAS-associated single-nucleotide variants (SNVs). We then performed quantitative assessment of the contribution of commonly used genomic features, including open chromatin maps, expression quantitative trait loci (eQTL), and chromatin conformation data. Using these features, we trained and validated an Effector Index (Ei), to map target genes for these 12 common diseases and traits. Ei demonstrated high predictive performance, both with cross-validation on the training set, and an independently derived set for type 2 diabetes. Key predictive features included coding or transcript-altering SNVs, distance to gene, and open chromatin-based metrics. This work outlines a simple, understandable approach to prioritize genes at GWAS loci for functional follow-up and drug development, and provides a systematic strategy for prioritization of GWAS target genes.

Identification of Novel Regulatory Regions Induced by Intrauterine Growth Restriction in Rat Islets

Intrauterine growth restriction (IUGR) leads to the development of type 2 diabetes in adulthood, and the permanent alterations in gene expression implicate an epigenetic mechanism. Using a rat model of IUGR, we performed TrueSeq-HELP Tagging to assess the association of DNA methylation changes and gene dysregulation in islets. We identified 511 differentially methylated regions (DMRs) and 4377 significantly altered single CpG sites. Integrating the methylome and our published transcriptome data sets resulted in the identification of pathways critical for islet function. The identified DMRs were enriched with transcription factor binding motifs, such as Elk1, Etv1, Foxa1, Foxa2, Pax7, Stat3, Hnf1, and AR. In silico analysis of 3-dimensional chromosomal interactions using human pancreas and islet Hi-C data sets identified interactions between 14 highly conserved DMRs and 35 genes with significant expression changes at an early age, many of which persisted in adult islets. In adult islets, there were far more interactions between DMRs and genes with significant expression changes identified with Hi-C, and most of them were critical to islet metabolism and insulin secretion. The methylome was integrated with our published genome-wide histone modification data sets from IUGR islets, resulting in further characterization of important regulatory regions of the genome altered by IUGR containing both significant changes in DNA methylation and specific histone marks. We identified novel regulatory regions in islets after exposure to IUGR, suggesting that epigenetic changes at key transcription factor binding motifs and other gene regulatory regions may contribute to gene dysregulation and an abnormal islet phenotype in IUGR rats.

Epigenomic and Transcriptomic Prioritization of Candidate Obesity-Risk Regulatory GWAS SNPs

Concern about rising rates of obesity has prompted searches for obesity-related single nucleotide polymorphisms (SNPs) in genome-wide association studies (GWAS). Identifying plausible regulatory SNPs is very difficult partially because of linkage disequilibrium. We used an unusual epigenomic and transcriptomic analysis of obesity GWAS-derived SNPs in adipose versus heterologous tissues. From 50 GWAS and 121,064 expanded SNPs, we prioritized 47 potential causal regulatory SNPs (Tier-1 SNPs) for 14 gene loci. A detailed examination of seven loci revealed that four (CABLES1, PC, PEMT, and FAM13A) had Tier-1 SNPs positioned so that they could regulate use of alternative transcription start sites, resulting in different polypeptides being generated or different amounts of an intronic microRNA gene being expressed. HOXA11 and long noncoding RNA gene RP11-392O17.1 had Tier-1 SNPs in their 3' or promoter region, respectively, and strong preferences for expression in subcutaneous versus visceral adipose tissue. ZBED3-AS1 had two intragenic Tier-1 SNPs, each of which could contribute to mediating obesity risk through modulating long-distance chromatin interactions. Our approach not only revealed especially credible novel regulatory SNPs, but also helped evaluate previously highlighted obesity GWAS SNPs that were candidates for transcription regulation.

Stem Cell-Derived β Cells: A Versatile Research Platform to Interrogate the Genetic Basis of β Cell Dysfunction

Pancreatic β cell dysfunction is a central component of diabetes progression. During the last decades, the genetic basis of several monogenic forms of diabetes has been recognized. Genome-wide association studies (GWAS) have also facilitated the identification of common genetic variants associated with an increased risk of diabetes. These studies highlight the importance of impaired β cell function in all forms of diabetes. However, how most of these risk variants confer disease risk, remains unanswered. Understanding the specific contribution of genetic variants and the precise role of their molecular effectors is the next step toward developing treatments that target β cell dysfunction in the era of personalized medicine. Protocols that allow derivation of β cells from pluripotent stem cells, represent a powerful research tool that allows modeling of human development and versatile experimental designs that can be used to shed some light on diabetes pathophysiology. This article reviews different models to study the genetic basis of β cell dysfunction, focusing on the recent advances made possible by stem cell applications in the field of diabetes research.

Oxidative stress-mediated beta cell death and dysfunction as a target for diabetes management

The biggest drawback of a current diabetes therapy is the treatment of the consequences not the cause of the disease. Regardless of the diabetes type, preservation and recovery of functional pancreatic beta cells stands as the biggest challenge in the treatment of diabetes. Free radicals and oxidative stress are among the major mediators of autoimmune destruction of beta cells in type 1 diabetes (T1D) or beta cell malfunction and death provoked by glucotoxicity and insulin resistance in type 2 diabetes (T2D). Additionally, oxidative stress reduces functionality of beta cells in T2D by stimulating their de-/trans-differentiation through the loss of transcription factors critical for beta cell development, maturity and regeneration. This review summarizes up to date clarified redox-related mechanisms involved in regulating beta cell identity and death, underlining similarities and differences between T1D and T2D. The protective effects of natural antioxidants on the oxidative stress-induced beta cell failure were also discussed. Considering that oxidative stress affects epigenetic regulatory mechanisms involved in the regulation of pancreatic beta cell survival and insulin secretion, this review highlighted huge potential of epigenetic therapy. Special attention was paid on application of the state-of-the-art CRISPR/Cas9 technology, based on targeted epigenome editing with the purpose of changing the differentiation state of different cell types, making them insulin-producing with ability to attenuate diabetes. Clarification of the above-mentioned mechanisms could provide better insight into diabetes etiology and pathogenesis, which would allow development of novel, potentially more efficient therapeutic strategies for the prevention or reversion of beta cell loss.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

BACH2 inhibition reverses β cell failure in type 2 diabetes models

Type 2 diabetes (T2D) is associated with defective insulin secretion and reduced β cell mass. Available treatments provide a temporary reprieve, but secondary failure rates are high, making insulin supplementation necessary. Reversibility of β cell failure is a key translational question. Here, we reverse engineered and interrogated pancreatic islet-specific regulatory networks to discover T2D-specific subpopulations characterized by metabolic inflexibility and endocrine progenitor/stem cell features. Single-cell gain- and loss-of-function and glucose-induced Ca2+ flux analyses of top candidate master regulatory (MR) proteins in islet cells validated transcription factor BACH2 and associated epigenetic effectors as key drivers of T2D cell states. BACH2 knockout in T2D islets reversed cellular features of the disease, restoring a nondiabetic phenotype. BACH2-immunoreactive islet cells increased approximately 4-fold in diabetic patients, confirming the algorithmic prediction of clinically relevant subpopulations. Treatment with a BACH inhibitor lowered glycemia and increased plasma insulin levels in diabetic mice, and restored insulin secretion in diabetic mice and human islets. The findings suggest that T2D-specific populations of failing β cells can be reversed and indicate pathways for pharmacological intervention, including via BACH2 inhibition.

Enhancing Acsl4 in absence of mTORC2/Rictor drove β-cell dedifferentiation via inhibiting FoxO1 and promoting ROS production

Rapamycin insensitive companion of mechanistic target of Rapamycin (Rictor), the key component of mTOR complex 2 (mTORC2), controls both β-cell proliferation and function. We sought to study whether long chain acyl-CoA synthetase 4 (Acsl4) worked downstream of Rictor/mTORC2 to maintain β-cell functional mass. We found Acsl4 was positively regulated by Rictor at transcriptional and posttranslational levels in mouse β-cell. Infecting adenovirus expressing Acsl4 in β-cell-specific-Rictor-knockout (βRicKO) islets and Min6 cells knocking down Rictor with lentivirus-expressing siRNA-oligos targeting Rictor(siRic), recovered the β-cell dysplasia but not dysfunction. Cell bioenergetic experiment performed with Seahorse XF showed that Acsl4 could not rescue the dampened glucose oxidation in Rictor-lacking β-cell, but further promoted lipid oxidation. Transposase-Accessible Chromatin (ATAC) and H3K27Ac chromatin immunoprecipitation (ChIP) sequencing studies reflected the epigenetic elevated molecular signature for β-cell dedifferentiation and mitigated oxidative defense/response. These results were confirmed by the observations of elevated acetylation and ubiquitination of FoxO1, increased protein levels of Gpx1 and Hif1an, excessive reactive oxygen species (ROS) production and diminished MafA in Acsl4 overexpressed Rictor-lacking β-cells. In these cells, antioxidant treatment significantly recovered MafA level and insulin content. Inducing lipid oxidation alone could not mimic the effect of Acsl4 in Rictor lacking β-cell. Our study suggested that Acsl4 function in β-cell was context dependent and might facilitate β-cell dedifferentiation with attenuated Rictor/mTORC2 activity or insulin signaling via posttranslational inhibiting FoxO1 and epigenetically enhancing ROS induced MafA degradation.

Gaining insight into metabolic diseases from human genetic discoveries

Human large-scale genetic association studies have identified sequence variations at thousands of genetic risk loci that are more common in patients with diverse metabolic disease compared with healthy controls. While these genetic associations have been replicated in multiple large cohorts and sometimes can explain up to 50% of heritability, the molecular and cellular mechanisms affected by common genetic variation associated with metabolic disease remains mostly unknown. A variety of new genome-wide data types, in conjunction with novel biostatistical and computational analytical methodologies and foundational experimental technologies, are paving the way for a principled approach to systematic variant-to-function (V2F) studies for metabolic diseases, turning associated regions into causal variants, cell types and states of action, effector genes, and cellular and physiological mechanisms. Identification of new target genes and cellular programs for metabolic risk loci will improve mechanistic understanding of disease biology and identification of novel therapeutic strategies.

Transcriptional changes and the role of ONECUT1 in hPSC pancreatic differentiation

Cell type specification during pancreatic development is tightly controlled by a transcriptional and epigenetic network. The precise role of most transcription factors, however, has been only described in mice. To convey such concepts to human pancreatic development, alternative model systems such as pancreatic in vitro differentiation of human pluripotent stem cells can be employed. Here, we analyzed stage-specific RNA-, ChIP-, and ATAC-sequencing data to dissect transcriptional and regulatory mechanisms during pancreatic development. Transcriptome and open chromatin maps of pancreatic differentiation from human pluripotent stem cells provide a stage-specific pattern of known pancreatic transcription factors and indicate ONECUT1 as a crucial fate regulator in pancreas progenitors. Moreover, our data suggest that ONECUT1 is also involved in preparing pancreatic progenitors for later endocrine specification. The dissection of the transcriptional and regulatory circuitry revealed an important role for ONECUT1 within such network and will serve as resource to study human development and disease.

Mutations and variants of ONECUT1 in diabetes

Genes involved in distinct diabetes types suggest shared disease mechanisms. Here we show that One Cut Homeobox 1 (ONECUT1) mutations cause monogenic recessive syndromic diabetes in two unrelated patients, characterized by intrauterine growth retardation, pancreas hypoplasia and gallbladder agenesis/hypoplasia, and early-onset diabetes in heterozygous relatives. Heterozygous carriers of rare coding variants of ONECUT1 define a distinctive subgroup of diabetic patients with early-onset, nonautoimmune diabetes, who respond well to diabetes treatment. In addition, common regulatory ONECUT1 variants are associated with multifactorial type 2 diabetes. Directed differentiation of human pluripotent stem cells revealed that loss of ONECUT1 impairs pancreatic progenitor formation and a subsequent endocrine program. Loss of ONECUT1 altered transcription factor binding and enhancer activity and NKX2.2/NKX6.1 expression in pancreatic progenitor cells. Collectively, we demonstrate that ONECUT1 controls a transcriptional and epigenetic machinery regulating endocrine development, involved in a spectrum of diabetes, encompassing monogenic (recessive and dominant) as well as multifactorial inheritance. Our findings highlight the broad contribution of ONECUT1 in diabetes pathogenesis, marking an important step toward precision diabetes medicine.

TIGER: The gene expression regulatory variation landscape of human pancreatic islets

Genome-wide association studies (GWASs) identified hundreds of signals associated with type 2 diabetes (T2D). To gain insight into their underlying molecular mechanisms, we have created the translational human pancreatic islet genotype tissue-expression resource (TIGER), aggregating >500 human islet genomic datasets from five cohorts in the Horizon 2020 consortium T2DSystems. We impute genotypes using four reference panels and meta-analyze cohorts to improve the coverage of expression quantitative trait loci (eQTL) and develop a method to combine allele-specific expression across samples (cASE). We identify >1 million islet eQTLs, 53 of which colocalize with T2D signals. Among them, a low-frequency allele that reduces T2D risk by half increases CCND2 expression. We identify eight cASE colocalizations, among which we found a T2D-associated SLC30A8 variant. We make all data available through the TIGER portal (http://tiger.bsc.es), which represents a comprehensive human islet genomic data resource to elucidate how genetic variation affects islet function and translates into therapeutic insight and precision medicine for T2D.

3D disorganization and rearrangement of genome provide insights into pathogenesis of NAFLD by integrated Hi-C, Nanopore, and RNA sequencing

The three-dimensional (3D) conformation of chromatin is integral to the precise regulation of gene expression. The 3D genome and genomic variations in non-alcoholic fatty liver disease (NAFLD) are largely unknown, despite their key roles in cellular function and physiological processes. High-throughput chromosome conformation capture (Hi-C), Nanopore sequencing, and RNA-sequencing (RNA-seq) assays were performed on the liver of normal and NAFLD mice. A high-resolution 3D chromatin interaction map was generated to examine different 3D genome hierarchies including A/B compartments, topologically associated domains (TADs), and chromatin loops by Hi-C, and whole genome sequencing identifying structural variations (SVs) and copy number variations (CNVs) by Nanopore sequencing. We identified variations in thousands of regions across the genome with respect to 3D chromatin organization and genomic rearrangements, between normal and NAFLD mice, and revealed gene dysregulation frequently accompanied by these variations. Candidate target genes were identified in NAFLD, impacted by genetic rearrangements and spatial organization disruption. Our data provide a high-resolution 3D genome interaction resource for NAFLD investigations, revealed the relationship among genetic rearrangements, spatial organization disruption, and gene regulation, and identified candidate genes associated with these variations implicated in the pathogenesis of NAFLD. The newly findings offer insights into novel mechanisms of NAFLD pathogenesis and can provide a new conceptual framework for NAFLD therapy.

Leveraging three-dimensional chromatin architecture for effective reconstruction of enhancer-target gene regulatory interactions

A growing amount of evidence in literature suggests that germline sequence variants and somatic mutations in non-coding distal regulatory elements may be crucial for defining disease risk and prognostic stratification of patients, in genetic disorders as well as in cancer. Their functional interpretation is challenging because genome-wide enhancer-target gene (ETG) pairing is an open problem in genomics. The solutions proposed so far do not account for the hierarchy of structural domains which define chromatin three-dimensional (3D) architecture. Here we introduce a change of perspective based on the definition of multi-scale structural chromatin domains, integrated in a statistical framework to define ETG pairs. In this work (i) we develop a computational and statistical framework to reconstruct a comprehensive map of ETG pairs leveraging functional genomics data; (ii) we demonstrate that the incorporation of chromatin 3D architecture information improves ETG pairing accuracy and (iii) we use multiple experimental datasets to extensively benchmark our method against previous solutions for the genome-wide reconstruction of ETG pairs. This solution will facilitate the annotation and interpretation of sequence variants in distal non-coding regulatory elements. We expect this to be especially helpful in clinically oriented applications of whole genome sequencing in cancer and undiagnosed genetic diseases research.

Functional characterization of T2D-associated SNP effects on baseline and ER stress-responsive β cell transcriptional activation

Genome-wide association studies (GWAS) have linked single nucleotide polymorphisms (SNPs) at >250 loci in the human genome to type 2 diabetes (T2D) risk. For each locus, identifying the functional variant(s) among multiple SNPs in high linkage disequilibrium is critical to understand molecular mechanisms underlying T2D genetic risk. Using massively parallel reporter assays (MPRA), we test the cis-regulatory effects of SNPs associated with T2D and altered in vivo islet chromatin accessibility in MIN6 β cells under steady state and pathophysiologic endoplasmic reticulum (ER) stress conditions. We identify 1,982/6,621 (29.9%) SNP-containing elements that activate transcription in MIN6 and 879 SNP alleles that modulate MPRA activity. Multiple T2D-associated SNPs alter the activity of short interspersed nuclear element (SINE)-containing elements that are strongly induced by ER stress. We identify 220 functional variants at 104 T2D association signals, narrowing 54 signals to a single candidate SNP. Together, this study identifies elements driving β cell steady state and ER stress-responsive transcriptional activation, nominates causal T2D SNPs, and uncovers potential roles for repetitive elements in β cell transcriptional stress response and T2D genetics.

A Comprehensive Toolbox to Analyze Enhancer-Promoter Functions

Knowledge in gene transcription and chromatin regulation has been intensely studied for decades, but thanks to next-generation sequencing (NGS) techniques there has been a major leap forward in the last few years. Historically, identification of specific enhancer elements has led to the identification of master transcription factors (TFs) in the 1990s. Genetic and biochemical experiments have identified the key regulators controlling RNA polymerase II (RNAPII) transcription and structurally analyses have elucidated detailed mechanisms. NGS and the development of chromatin immunoprecipitation (ChIP) have accelerated the gain of knowledge in the recent years. By now, we have a dazzling wealth of techniques that are currently used to put gene expression into a genome-wide context. This book is an attempt to assemble useful protocols for many researchers within and nearby research areas. In general, these innovative techniques focus on enhancer and promoter studies. The techniques should also be of interest for related fields such as DNA repair and replication.

Spatial genome architecture and the emergence of malignancy

Human chromosomes are large spatially and hierarchically structured entities, the integrity of which needs to be preserved throughout the lifespan of the cell and in conjunction with cell cycle progression. Preservation of chromosomal structure is important for proper deployment of cell type-specific gene expression programs. Thus, aberrations in the integrity and structure of chromosomes will predictably lead to disease, including cancer. Here, we provide an updated standpoint with respect to chromatin misfolding and the emergence of various cancer types. We discuss recent studies implicating the disruption of topologically associating domains, switching between active and inactive compartments, rewiring of promoter-enhancer interactions in malignancy as well as the effects of single nucleotide polymorphisms in non-coding regions involved in long-range regulatory interactions. In light of these findings, we argue that chromosome conformation studies may now also be useful for patient diagnosis and drug target discovery.

Enhancers in disease: molecular basis and emerging treatment strategies

Enhancers are genomic sequences that play a key role in regulating tissue-specific gene expression levels. An increasing number of diseases are linked to impaired enhancer function through chromosomal rearrangement, genetic variation within enhancers, or epigenetic modulation. Here, we review how these enhancer disruptions have recently been implicated in congenital disorders, cancers, and common complex diseases and address the implications for diagnosis and treatment. Although further fundamental research into enhancer function, target genes, and context is required, enhancer-targeting drugs and gene editing approaches show great therapeutic promise for a range of diseases.

INFIMA leverages multi-omics model organism data to identify effector genes of human GWAS variants

Genome-wide association studies reveal many non-coding variants associated with complex traits. However, model organism studies largely remain as an untapped resource for unveiling the effector genes of non-coding variants. We develop INFIMA, Integrative Fine-Mapping, to pinpoint causal SNPs for diversity outbred (DO) mice eQTL by integrating founder mice multi-omics data including ATAC-seq, RNA-seq, footprinting, and in silico mutation analysis. We demonstrate INFIMA's superior performance compared to alternatives with human and mouse chromatin conformation capture datasets. We apply INFIMA to identify novel effector genes for GWAS variants associated with diabetes. The results of the application are available at http://www.statlab.wisc.edu/shiny/INFIMA/ .

REST is a major negative regulator of endocrine differentiation during pancreas organogenesis

In this study, Rovira et al. report that inactivation of the transcriptional repressor REST causes a drastic increase in pancreatic endocrine progenitors and endocrine cells, and establish that REST is a major negative regulator of embryonic pancreas endocrine differentiation in mice and zebrafish. Their findings show that REST-dependent inhibition ensures a balanced production of endocrine cells from embryonic pancreatic progenitors.

Multiple transcription factors have been shown to promote pancreatic β-cell differentiation, yet much less is known about negative regulators. Earlier epigenomic studies suggested that the transcriptional repressor REST could be a suppressor of endocrinogenesis in the embryonic pancreas. However, pancreatic Rest knockout mice failed to show abnormal numbers of endocrine cells, suggesting that REST is not a major regulator of endocrine differentiation. Using a different conditional allele that enables profound REST inactivation, we observed a marked increase in pancreatic endocrine cell formation. REST inhibition also promoted endocrinogenesis in zebrafish and mouse early postnatal ducts and induced β-cell-specific genes in human adult duct-derived organoids. We also defined genomic sites that are bound and repressed by REST in the embryonic pancreas. Our findings show that REST-dependent inhibition ensures a balanced production of endocrine cells from embryonic pancreatic progenitors.

The molecular basis, genetic control and pleiotropic effects of local gene co-expression

Nearby genes are often expressed as a group. Yet, the prevalence, molecular mechanisms and genetic control of local gene co-expression are far from being understood. Here, by leveraging gene expression measurements across 49 human tissues and hundreds of individuals, we find that local gene co-expression occurs in 13% to 53% of genes per tissue. By integrating various molecular assays (e.g. ChIP-seq and Hi-C), we estimate the ability of several mechanisms, such as enhancer-gene interactions, in distinguishing gene pairs that are co-expressed from those that are not. Notably, we identify 32,636 expression quantitative trait loci (eQTLs) which associate with co-expressed gene pairs and often overlap enhancer regions. Due to affecting several genes, these eQTLs are more often associated with multiple human traits than other eQTLs. Our study paves the way to comprehend trait pleiotropy and functional interpretation of QTL and GWAS findings. All local gene co-expression identified here is available through a public database ( https://glcoex.unil.ch/ ).