Comparative analyses of super-enhancers reveal conserved elements in vertebrate genomes

Evolutionary conservation of VSX2 super-enhancer modules in retinal development

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. We recently identified developmental stage- and cell type-specific modules within the murine Vsx2 SE. Here, we show that the human VSX2 SE modules have similar developmental stage- and cell type-specific activity in reporter gene assays. By inserting the human sequence of one VSX2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the function of these SE modules during human retinal development, we deleted individual modules in human embryonic stem cells and generated retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module led to disruptions in bipolar neuron development. This prototypical SE serves as a model for understanding developmental stage- and cell type-specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

Computational identification of key transcription factors for embryonic and postnatal Sox2+ dental epithelial stem cell

While many reptiles can replace their tooth throughout life, human loss the tooth replacement capability after formation of the permanent teeth. It was thought that the difference in tooth regeneration capability depends on the persistence of a specialized dental epithelial structure, the dental lamina that contains dental epithelial stem cells (DESC). Currently, we know very little about DESC such as what genes are expressed or its chromatin accessibility profile. Multiple markers of DESC have been proposed such as Sox2 and Lgr5 . Few single cell RNA-seq experiments have been performed previously, but no obvious DESC cluster was identified in these scRNA-seq datasets, possible due to that the expression level of DESC markers such as Sox2 and Lgr5 is too low or the percentage of DESC is too low in whole tooth. We utilize a mouse line Sox2-GFP to enrich Sox2+ DESC and use Smart-Seq2 protocol and ATAC-seq protocol to generate transcriptome profile and chromatin accessibility profile of P2 Sox2+ DESC. Additionally, we generate transcriptome profile and chromatin accessibility profile of E11.5 Sox2+ dental lamina cells. With transcriptome profile and chromatin accessibility profile, we systematically identify potential key transcription factors for E11.5 Sox2+ cells and P2 Sox2+ cells. We identified transcription factors including Pitx2, Id3, Pitx1, Tbx1, Trp63, Nkx2-3, Grhl3, Dlx2, Runx1, Nfix, Zfp536 , etc potentially formed the core transcriptional regulatory networks of Sox2+ DESC in both embryonic and postnatal stages.

Identification of Evolutionarily Conserved VSX2 Enhancers in Retinal Development

Super-enhancers (SEs) are expansive regions of genomic DNA that regulate the expression of genes involved in cell identity and cell fate. Recently, we found that distinct modules within a murine SE regulate gene expression of master regulatory transcription factor Vsx2 in a developmental stage- and cell-type specific manner. Vsx2 is expressed in retinal progenitor cells as well as differentiated bipolar neurons and Müller glia. Mutations in VSX2 in humans and mice lead to microphthalmia due to a defect in retinal progenitor cell proliferation. Deletion of a single module within the Vsx2 SE leads to microphthalmia. Deletion of a separate module within the SE leads to a complete loss of bipolar neurons, yet the remainder of the retina develops normally. Furthermore, the Vsx2 SE is evolutionarily conserved in vertebrates, suggesting that these modules are important for retinal development across species. In the present study, we examine the ability of these modules to drive retinal development between species. By inserting the human build of one Vsx2 SE module into a mouse with microphthalmia, eye size was rescued. To understand the implications of these SE modules in a model of human development, we generated human retinal organoids. Deleting one module results in small organoids, recapitulating the small-eyed phenotype of mice with microphthalmia, while deletion of the other module leads to a complete loss of ON cone bipolar neurons. This prototypical SE serves as a model for uncoupling developmental stage- and cell-type specific effects of neurogenic transcription factors with complex expression patterns. Moreover, by elucidating the gene regulatory mechanisms, we can begin to examine how dysregulation of these mechanisms contributes to phenotypic diversity and disease.

Zebrafish regulatory genomic resources for disease modelling and regeneration

In the past decades, the zebrafish has become a disease model with increasing popularity owing to its advantages that include fast development, easy genetic manipulation, simplicity for imaging, and sharing conserved disease-associated genes and pathways with those of human. In parallel, studies of disease mechanisms are increasingly focusing on non-coding mutations, which require genome annotation maps of regulatory elements, such as enhancers and promoters. In line with this, genomic resources for zebrafish research are expanding, producing a variety of genomic data that help in defining regulatory elements and their conservation between zebrafish and humans. Here, we discuss recent developments in generating functional annotation maps for regulatory elements of the zebrafish genome and how this can be applied to human diseases. We highlight community-driven developments, such as DANIO-CODE, in generating a centralised and standardised catalogue of zebrafish genomics data and functional annotations; consider the advantages and limitations of current annotation maps; and offer considerations for interpreting and integrating existing maps with comparative genomics tools. We also discuss the need for developing standardised genomics protocols and bioinformatic pipelines and provide suggestions for the development of analysis and visualisation tools that will integrate various multiomic bulk sequencing data together with fast-expanding data on single-cell methods, such as single-cell assay for transposase-accessible chromatin with sequencing. Such integration tools are essential to exploit the multiomic chromatin characterisation offered by bulk genomics together with the cell-type resolution offered by emerging single-cell methods. Together, these advances will build an expansive toolkit for interrogating the mechanisms of human disease in zebrafish.

Experimental Validation and Prediction of Super-Enhancers: Advances and Challenges

Super-enhancers (SEs) are cis-regulatory elements of the human genome that have been widely discussed since the discovery and origin of the term. Super-enhancers have been shown to be strongly associated with the expression of genes crucial for cell differentiation, cell stability maintenance, and tumorigenesis. Our goal was to systematize research studies dedicated to the investigation of structure and functions of super-enhancers as well as to define further perspectives of the field in various applications, such as drug development and clinical use. We overviewed the fundamental studies which provided experimental data on various pathologies and their associations with particular super-enhancers. The analysis of mainstream approaches for SE search and prediction allowed us to accumulate existing data and propose directions for further algorithmic improvements of SEs' reliability levels and efficiency. Thus, here we provide the description of the most robust algorithms such as ROSE, imPROSE, and DEEPSEN and suggest their further use for various research and development tasks. The most promising research direction, which is based on topic and number of published studies, are cancer-associated super-enhancers and prospective SE-targeted therapy strategies, most of which are discussed in this review.

Polyglutamine-expanded ATXN7 alters a specific epigenetic signature underlying photoreceptor identity gene expression in SCA7 mouse retinopathy

BACKGROUND

Spinocerebellar ataxia type 7 (SCA7) is a neurodegenerative disorder that primarily affects the cerebellum and retina. SCA7 is caused by a polyglutamine expansion in the ATXN7 protein, a subunit of the transcriptional coactivator SAGA that acetylates histone H3 to deposit narrow H3K9ac mark at DNA regulatory elements of active genes. Defective histone acetylation has been presented as a possible cause for gene deregulation in SCA7 mouse models. However, the topography of acetylation defects at the whole genome level and its relationship to changes in gene expression remain to be determined.

METHODS

We performed deep RNA-sequencing and chromatin immunoprecipitation coupled to high-throughput sequencing to examine the genome-wide correlation between gene deregulation and alteration of the active transcription marks, e.g. SAGA-related H3K9ac, CBP-related H3K27ac and RNA polymerase II (RNAPII), in a SCA7 mouse retinopathy model.

RESULTS

Our analyses revealed that active transcription marks are reduced at most gene promoters in SCA7 retina, while a limited number of genes show changes in expression. We found that SCA7 retinopathy is caused by preferential downregulation of hundreds of highly expressed genes that define morphological and physiological identities of mature photoreceptors. We further uncovered that these photoreceptor genes harbor unusually broad H3K9ac profiles spanning the entire gene bodies and have a low RNAPII pausing. This broad H3K9ac signature co-occurs with other features that delineate superenhancers, including broad H3K27ac, binding sites for photoreceptor specific transcription factors and expression of enhancer-related non-coding RNAs (eRNAs). In SCA7 retina, downregulated photoreceptor genes show decreased H3K9 and H3K27 acetylation and eRNA expression as well as increased RNAPII pausing, suggesting that superenhancer-related features are altered.

CONCLUSIONS

Our study thus provides evidence that distinctive epigenetic configurations underlying high expression of cell-type specific genes are preferentially impaired in SCA7, resulting in a defect in the maintenance of identity features of mature photoreceptors. Our results also suggest that continuous SAGA-driven acetylation plays a role in preserving post-mitotic neuronal identity.

Regeneration and developmental enhancers are differentially compatible with minimal promoters

Enhancers and promoters are cis-regulatory elements that control gene expression. Enhancers are activated in a cell type-, tissue-, and condition-specific manner to stimulate promoter function and transcription. Zebrafish have emerged as a powerful animal model for examining the activities of enhancers derived from various species through transgenic enhancer assays, in which an enhancer is coupled with a minimal promoter. However, the efficiency of minimal promoters and their compatibility with multiple developmental and regeneration enhancers have not been systematically tested in zebrafish. Thus, we assessed the efficiency of six minimal promoters and comprehensively interrogated the compatibility of the promoters with developmental and regeneration enhancers. We found that the fos minimal promoter and Drosophila synthetic core promoter (DSCP) yielded high rates of leaky expression that may complicate the interpretation of enhancer assays. Notably, the adenovirus E1b promoter, the zebrafish lepb 0.8-kb (P0.8) and lepb 2-kb (P2) promoters, and a new zebrafish synthetic promoter (ZSP) that combines elements of the E1b and P0.8 promoters drove little or no ectopic expression, making them suitable for transgenic assays. We also found significant differences in compatibility among specific combinations of promoters and enhancers, indicating the importance of promoters as key regulatory elements determining the specificity of gene expression. Our study provides guidelines for transgenic enhancer assays in zebrafish to aid in the discovery of functional enhancers regulating development and regeneration.

Identification and functional validation of super-enhancers in Arabidopsis thaliana

Super-enhancers (SEs) are exceptionally large enhancers and are recognized to play prominent roles in cell identity in mammalian species. We surveyed the genomic regions containing large clusters of accessible chromatin regions (ACRs) marked by deoxyribonuclease (DNase) I hypersensitivity in Arabidopsis thaliana. We identified a set of 749 putative SEs, which have a minimum length of 1.5 kilobases and represent the top 2.5% of the largest ACR clusters. We demonstrate that the genomic regions associating with these SEs were more sensitive to DNase I than other nonpromoter ACRs. The SEs were preferentially associated with topologically associating domains. Furthermore, the SEs and their predicted cognate genes were frequently associated with organ development and tissue identity in A. thaliana. Therefore, the A. thaliana SEs and their cognate genes mirror the functional characteristics of those reported in mammalian species. We developed CRISPR/Cas-mediated deletion lines of a 3,578-bp SE associated with the thalianol biosynthetic gene cluster (BGC). Small deletions (131-157 bp) within the SE resulted in distinct phenotypic changes and transcriptional repression of all five thalianol genes. In addition, T-DNA insertions in the SE region resulted in transcriptional alteration of all five thalianol genes. Thus, this SE appears to play a central role in coordinating the operon-like expression pattern of the thalianol BGC.

Super-enhancers conserved within placental mammals maintain stem cell pluripotency

Despite pluripotent stem cells sharing key transcription factors, their maintenance involves distinct genetic inputs. Emerging evidence suggests that super-enhancers (SEs) can function as master regulatory hubs to control cell identity and pluripotency in humans and mice. However, whether pluripotency-associated SEs share an evolutionary origin in mammals remains elusive. Here, we performed comprehensive comparative epigenomic and transcription factor binding analyses among pigs, humans, and mice to identify pluripotency-associated SEs. Like typical enhancers, SEs displayed rapid evolution in mammals. We showed that BRD4 is an essential and conserved activator for mammalian pluripotency-associated SEs. Comparative motif enrichment analysis revealed 30 shared transcription factor binding motifs among the three species. The majority of transcriptional factors that bind to identified motifs are known regulators associated with pluripotency. Further, we discovered three pluripotency-associated SEs (SE-SOX2, SE-PIM1, and SE-FGFR1) that displayed remarkable conservation in placental mammals and were sufficient to drive reporter gene expression in a pluripotency-dependent manner. Disruption of these conserved SEs through the CRISPR-Cas9 approach severely impaired stem cell pluripotency. Our study provides insights into the understanding of conserved regulatory mechanisms underlying the maintenance of pluripotency as well as species-specific modulation of the pluripotency-associated regulatory networks in mammals.

Senescent cells limit p53 activity via multiple mechanisms to remain viable

Super-enhancers regulate genes with important functions in processes that are cell type-specific or define cell identity. Mouse embryonic fibroblasts establish 40 senescence-associated super-enhancers regardless of how they become senescent, with 50 activated genes located in the vicinity of these enhancers. Here we show, through gene knockdown and analysis of three core biological properties of senescent cells that a relatively large number of senescence-associated super-enhancer-regulated genes promote survival of senescent mouse embryonic fibroblasts. Of these, Mdm2, Rnase4, and Ang act by suppressing p53-mediated apoptosis through various mechanisms that are also engaged in response to DNA damage. MDM2 and RNASE4 transcription is also elevated in human senescent fibroblasts to restrain p53 and promote survival. These insights identify key survival mechanisms of senescent cells and provide molecular entry points for the development of targeted therapeutics that eliminate senescent cells at sites of pathology.

Role of CNSs Conserved Distal Cis-Regulatory Elements in CD4 + T Cell Development and Differentiation

Naïve CD4⁺ T cells differentiate into diverse subsets of effector cells and perform various homeostatic and immune functions. The differentiation and maintenance of these different subsets are controlled through the upregulation and silencing of master genes. Mechanistic studies of the regulation of these master genes identified conserved and distal intronic regulatory elements, which are accessible subsets of conserved non-coding sequences (CNSs), acting as cis-regulatory elements in a lineage-specific manner that controls the function of CD4⁺ T cells. Abnormal CNS activity is associated with incorrect expression of master genes and development of autoimmune diseases or immune suppression. Here, we describe the function of several conserved, distal cis-regulatory elements at the Foxp3, Rorc, Il-4, Il-10 and Il-17 gene locus were shown to play important roles in CD4⁺ T cells differentiation. Together, this review briefly outlines currently known CNSs, with a focus on their regulations and functions in complexes modulating the differentiation and maintenance of various CD4⁺ T cells subsets, in health and disease contexts, as well as during the conversion of T regulatory cells to T helper 17 cells. This article will provide a comprehensive view of CNSs conserved distal cis-regulatory elements at a few loci that control aspects of CD4⁺ T cells function.

Dynamic 3D genome reorganization during development and metabolic stress of the porcine liver

Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

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.

Dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis

Folliculogenesis is a complex biological process involving a central oocyte and its surrounding somatic cells. Three-dimensional chromatin architecture is an important transcription regulator; however, little is known about its dynamics and role in transcriptional regulation of granulosa cells during chicken folliculogenesis. We investigate the transcriptomic dynamics of chicken granulosa cells over ten follicular stages and assess the chromatin architecture dynamics and how it influences gene expression in granulosa cells at three key stages: the prehierarchical small white follicles, the first largest preovulatory follicles, and the postovulatory follicles. Our results demonstrate the consistency between the global reprogramming of chromatin architecture and the transcriptomic divergence during folliculogenesis, providing ample evidence for compartmentalization rearrangement, variable organization of topologically associating domains, and rewiring of the long-range interaction between promoter and enhancers. These results provide key insights into avian reproductive biology and provide a foundational dataset for the future in-depth functional characterization of granulosa cells.

The domestic chicken Gallus gallus domesticus is a classic model for the study of folliculogenesis. Here the authors integrate multi-omics analyses characterizing the dynamic transcriptome and chromatin architecture in granulosa cells during chicken folliculogenesis.

Genome-wide chromatin contacts of super-enhancer-associated lncRNA identify LINC01013 as a regulator of fibrosis in the aortic valve

Calcific aortic valve disease (CAVD) is characterized by a fibrocalcific process. The regulatory mechanisms that drive the fibrotic response in the aortic valve (AV) are poorly understood. Long noncoding RNAs derived from super-enhancers (lncRNA-SE) control gene expression and cell fate. Herein, multidimensional profiling including chromatin immunoprecipitation and sequencing, transposase-accessible chromatin sequencing, genome-wide 3D chromatin contacts of enhancer-promoter identified LINC01013 as an overexpressed lncRNA-SE during CAVD. LINC01013 is within a loop anchor, which has contact with the promoter of CCN2 (CTGF) located at ~180 kb upstream. Investigation showed that LINC01013 acts as a decoy factor for the negative transcription elongation factor E (NELF-E), whereby it controls the expression of CCN2. LINC01013-CCN2 is part of a transforming growth factor beta 1 (TGFB1) network and exerts a control over fibrogenesis. These findings illustrate a novel mechanism whereby a dysregulated lncRNA-SE controls, through a looping process, the expression of CCN2 and fibrogenesis of the AV.

Calcific aortic valve disease is the most common heart valve disorder characterized by a thickening of the aortic valve resulting from fibrotic and calcific processes. Because the aortic valve replacement is currently the only therapeutic option, the identification of key molecular processes that control the progression of the disease could lead to the development of novel noninvasive therapies. Growing evidence suggests that long noncoding RNAs (lncRNAs) fine tune gene expression in health and disease states. By using a multidimensional profiling including genome-wide 3D enhancer-promoter looping data, we identified LINC01013, a lncRNA, as a regulator of fibrogenesis. Specifically, we found that LINC01013 is located in a cluster of distant enhancers (super-enhancer) in aortic valve interstitial cells and has significant long-range looping with the promoter of CCN2, a gene that orchestrates fibrogenesis. We discovered that LINC01013 is acting as a decoy factor for a negative transcription elongation factor, whereby it controls the transcription of CCN2. In turn, higher expression of LINC01013 during calcific aortic valve disease promoted the expression of CCN2 and a fibrogenic program. These findings provide evidence that LINC01013 is a key regulator of fibrogenesis in CAVD.

A Comparative Analysis of Super-Enhancers and Broad H3K4me3 Domains in Pig, Human, and Mouse Tissues

Super-enhancers (SEs) and broad H3K4me3 domains (BDs) are crucial regulators in the control of tissue identity in human and mouse. However, their features in pig remain largely unknown. In this study, by integrative computational analyses of epigenomic and transcriptomic data, we have characterized SEs and BDs in six pig tissues and analyzed their conservation in comparison with human and mouse tissues. Similar to human and mouse, pig SEs and BDs display higher tissue specificity than their typical counterparts. Genes proximal to SEs and BDs are associated with tissue identity in most tissues. About 55-182 SEs (5-17% in total) and 99-309 BDs (8-16% in total) across pig tissues are considered as functionally conserved elements because they have orthologous SEs and BDs in human and mouse. However, these elements do not necessarily exhibit sequence conservation. The functionally conserved SEs are correlated to tissue identity in majority of pig tissues, while those conserved BDs are linked to tissue identity in a few tissues. Our study provides resources for future gene regulatory studies in pig. It highlights that SEs are more effective in defining tissue identity than BDs, which is contrasting to a previous study. It also provides novel insights on understanding the sequence features of functionally conserved elements.

Computational biology approaches for mapping transcriptional regulatory networks

Transcriptional Regulatory Networks (TRNs) are mainly responsible for the cell-type- or cell-state-specific expression of gene sets from the same DNA sequence. However, so far there are no precise maps of TRNs available for each cell-type or cell-state, and no ideal tool to map those networks clearly and in full from biological samples. In this review, major approaches and tools to map TRNs from high-throughput data are presented, depending on the type of methods or data used to infer them, and their advantages and limitations are discussed. After summarizing the main principles defining the topology and structure–function relationships in TRNs, an overview of the extensive work done to map TRNs from bulk transcriptomic data will be presented by type of methodological approach. Most recent modellings of TRNs using other types of molecular data or integrating different data types, including single-cell RNA-sequencing and chromatin information, will then be discussed, before briefly concluding with improvements expected to come in the field.

RNA polymerase II clusters form in line with surface condensation on regulatory chromatin

It is essential for cells to control which genes are transcribed into RNA. In eukaryotes, two major control points are recruitment of RNA polymerase II (Pol II) into a paused state, and subsequent pause release toward transcription. Pol II recruitment and pause release occur in association with macromolecular clusters, which were proposed to be formed by a liquid-liquid phase separation mechanism. How such a phase separation mechanism relates to the interaction of Pol II with DNA during recruitment and transcription, however, remains poorly understood. Here, we use live and super-resolution microscopy in zebrafish embryos to reveal Pol II clusters with a large variety of shapes, which can be explained by a theoretical model in which regulatory chromatin regions provide surfaces for liquid-phase condensation at concentrations that are too low for canonical liquid-liquid phase separation. Model simulations and chemical perturbation experiments indicate that recruited Pol II contributes to the formation of these surface-associated condensates, whereas elongating Pol II is excluded from these condensates and thereby drives their unfolding.

Reorganization of chromatin architecture during prenatal development of porcine skeletal muscle

Myofibres (primary and secondary myofibre) are the basic structure of muscle and the determinant of muscle mass. To explore the skeletal muscle developmental processes from primary myofibres to secondary myofibres in pigs, we conducted an integrative three-dimensional structure of genome and transcriptomic characterization of longissimus dorsi muscle of pig from primary myofibre formation stage [embryonic Day 35 (E35)] to secondary myofibre formation stage (E80). In the hierarchical genomic structure, we found that 11.43% of genome switched compartment A/B status, 14.53% of topologically associating domains are changed intradomain interactions (D-scores) and 2,730 genes with differential promoter–enhancer interactions and (or) enhancer activity from E35 to E80. The alterations of genome architecture were found to correlate with expression of genes that play significant roles in neuromuscular junction, embryonic morphogenesis, skeletal muscle development or metabolism, typically, NEFL, MuSK, SLN, Mef2D and GCK. Significantly, Sox6 and MATN2 play important roles in the process of primary to secondary myofibres formation and increase the regulatory potential score and genes expression in it. In brief, we reveal the genomic reorganization from E35 to E80 and construct genome-wide high-resolution interaction maps that provide a resource for studying long-range control of gene expression from E35 to E80.

Chromatin architecture transitions from zebrafish sperm through early embryogenesis

Chromatin architecture mapping in 3D formats has increased our understanding of how regulatory sequences and gene expression are connected and regulated in a genome. The 3D chromatin genome shows extensive remodeling during embryonic development, and although the cleavage-stage embryos of most species lack structure before zygotic genome activation (pre-ZGA), zebrafish has been reported to have structure. Here, we aimed to determine the chromosomal architecture in paternal/sperm zebrafish gamete cells to discern whether it either resembles or informs early pre-ZGA zebrafish embryo chromatin architecture. First, we assessed the higher-order architecture through advanced low-cell in situ Hi-C. The structure of zebrafish sperm, packaged by histones, lacks topological associated domains and instead displays “hinge-like” domains of ∼150 kb that repeat every 1–2 Mbs, suggesting a condensed repeating structure resembling mitotic chromosomes. The pre-ZGA embryos lacked chromosomal structure, in contrast to prior work, and only developed structure post-ZGA. During post-ZGA, we find chromatin architecture beginning to form at small contact domains of a median length of ∼90 kb. These small contact domains are established at enhancers, including super-enhancers, and chemical inhibition of Ep300a (p300) and Crebbpa (CBP) activity, lowering histone H3K27ac, but not transcription inhibition, diminishes these contacts. Together, this study reveals hinge-like domains in histone-packaged zebrafish sperm chromatin and determines that the initial formation of high-order chromatin architecture in zebrafish embryos occurs after ZGA primarily at enhancers bearing high H3K27ac.

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

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

A map of cis-regulatory elements and 3D genome structures in zebrafish

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species¹. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Chromatin accessibility established by Pou5f3, Sox19b and Nanog primes genes for activity during zebrafish genome activation

In many organisms, early embryonic development is driven by maternally provided factors until the controlled onset of transcription during zygotic genome activation. The regulation of chromatin accessibility and its relationship to gene activity during this transition remain poorly understood. Here, we generated chromatin accessibility maps with ATAC-seq from genome activation until the onset of lineage specification. During this period, chromatin accessibility increases at regulatory elements. This increase is independent of RNA polymerase II-mediated transcription, with the exception of the hypertranscribed miR-430 locus. Instead, accessibility often precedes the transcription of associated genes. Loss of the maternal transcription factors Pou5f3, Sox19b, and Nanog, which are known to be required for zebrafish genome activation, results in decreased accessibility at regulatory elements. Importantly, the accessibility of regulatory regions, especially when established by Pou5f3, Sox19b and Nanog, is predictive for future transcription. Our results show that the maternally provided transcription factors Pou5f3, Sox19b, and Nanog open up chromatin and prime genes for activity during zygotic genome activation in zebrafish.

Decoding regulatory structures and features from epigenomics profiles: A Roadmap-ENCODE Variational Auto-Encoder (RE-VAE) model

The development of chromatin immunoprecipitation (ChIP) with massively parallel DNA sequencing (ChIP-seq) technologies has promoted generation of large-scale epigenomics data, providing us unprecedented opportunities to explore the landscape of epigenomic profiles at scales across both histone marks and tissue types. In addition to many tools directly for data analysis, advanced computational approaches, such as deep learning, have recently become promising to deeply mine the data structures and identify important regulators from complex functional genomics data. We implemented a neural network framework, a Variational Auto-Encoder (VAE) model, to explore the epigenomic data from the Roadmap Epigenomics Project and the Encyclopedia of DNA Elements (ENCODE) project. Our model is applied to 935 reference samples, covering 28 tissues and 12 histone marks. We used the enhancer and promoter regions as the annotation features and ChIP-seq signal values in these regions as the feature values. Through a parameter sweep process, we identified the suitable hyperparameter values and built a VAE model to represent the epigenomics data and to further explore the biological regulation. The resultant Roadmap-ENCODE VAE (RE-VAE) model contained data compression and feature representation. Using the compressed data in the latent space, we found that the majority of histone marks were well clustered but not for tissues or cell types. Tissue or cell specificity was observed only in some histone marks (e.g., H3K4me3 and H3K27ac) and could be characterized when the number of tissue samples is large (e.g., blood and brain). In blood, the contributive regions and genes identified by RE-VAE model were confirmed by tissue-specificity enrichment analysis with an independent tissue expression panel. Finally, we demonstrated that RE-VAE model could detect cancer cell lines with similar epigenomics profiles. In conclusion, we introduced and implemented a VAE model to represent large-scale epigenomics data. The model could be used to explore classifications of histone modifications and tissue/cell specificity and to classify new data with unknown sources.

Strategies for genetic inactivation of long noncoding RNAs in zebrafish

The number of annotated long noncoding RNAs (lncRNAs) continues to grow; however, their functional characterization in model organisms has been hampered by the lack of reliable genetic inactivation strategies. While partial or full deletions of lncRNA loci disrupt lncRNA expression, they do not permit the formal association of a phenotype with the encoded transcript. Here, we examined several alternative strategies for generating lncRNA null alleles in zebrafish and found that they often resulted in unpredicted changes to lncRNA expression. Removal of the transcription start sites (TSSs) of lncRNA genes resulted in hypomorphic mutants, due to the usage of either constitutive or tissue-specific alternative TSSs. Deletions of short, highly conserved lncRNA regions can also lead to overexpression of truncated transcripts. In contrast, knock-in of a polyadenylation signal enabled complete inactivation of malat1, the most abundant vertebrate lncRNA. In summary, lncRNA null alleles require extensive in vivo validation, and we propose insertion of transcription termination sequences as the most reliable approach to generate lncRNA-deficient zebrafish.

Long-Range Enhancer Interactions Are Prevalent in Mouse Embryonic Stem Cells and Are Reorganized upon Pluripotent State Transition

Transcriptional enhancers, including super-enhancers (SEs), form physical interactions with promoters to regulate cell-type-specific gene expression. SEs are characterized by high transcription factor occupancy and large domains of active chromatin, and they are commonly assigned to target promoters using computational predictions. How promoter-SE interactions change upon cell state transitions, and whether transcription factors maintain SE interactions, have not been reported. Here, we used promoter-capture Hi-C to identify promoters that interact with SEs in mouse embryonic stem cells (ESCs). We found that SEs form complex, spatial networks in which individual SEs contact multiple promoters, and a rewiring of promoter-SE interactions occurs between pluripotent states. We also show that long-range promoter-SE interactions are more prevalent in ESCs than in epiblast stem cells (EpiSCs) or Nanog-deficient ESCs. We conclude that SEs form cell-type-specific interaction networks that are partly dependent on core transcription factors, thereby providing insights into the gene regulatory organization of pluripotent cells.

•

Promoter-capture Hi-C identifies 3D interactions in mouse pluripotent cells

•

Super-enhancers (SEs) form complex spatial networks contacting multiple promoters

•

Rewiring of promoter-SE interactions between ESC and EpiSC pluripotent states

•

Long-range SE interactions are a hallmark of mouse ESCs

Packaging development: how chromatin controls transcription in zebrafish embryogenesis

How developmental gene expression is activated, co-ordinated and maintained is one of the biggest questions in developmental biology. While transcription factors lead the way in directing developmental gene expression, their accessibility to the correct repertoire of genes can depend on other factors such as DNA methylation, the presence of particular histone variants and post-translational modifications of histones. Collectively, factors that modify DNA or affect its packaging and accessibility contribute to a chromatin landscape that helps to control the timely expression of developmental genes. Zebrafish, perhaps better known for their strength as a model of embryology and organogenesis during development, are coming to the fore as a powerful model for interpreting the role played by chromatin in gene expression. Several recent advances have shown that zebrafish exhibit both similarities and differences to other models (and humans) in the way that they employ chromatin mechanisms of gene regulation. Here, I review how chromatin influences developmental transcriptional programmes during early zebrafish development, patterning and organogenesis. Lastly, I briefly highlight the importance of zebrafish chromatin research towards the understanding of human disease and transgenerational inheritance.

High-resolution comparative analysis of great ape genomes

Genetic studies of human evolution require high-quality contiguous ape genome assemblies that are not guided by the human reference. We coupled long-read sequence assembly and full-length complementary DNA sequencing with a multiplatform scaffolding approach to produce ab initio chimpanzee and orangutan genome assemblies. By comparing these with two long-read de novo human genome assemblies and a gorilla genome assembly, we characterized lineage-specific and shared great ape genetic variation ranging from single- to mega-base pair-sized variants. We identified ~17,000 fixed human-specific structural variants identifying genic and putative regulatory changes that have emerged in humans since divergence from nonhuman apes. Interestingly, these variants are enriched near genes that are down-regulated in human compared to chimpanzee cerebral organoids, particularly in cells analogous to radial glial neural progenitors.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.