Enhancers: the abundance and function of regulatory sequences beyond promoters

STAT3 drives the expression of HIF1alpha in cancer cells through a novel super-enhancer

Aerobic glycolysis is one of the major hallmarks of malignant tumors. This metabolic reprogramming benefits the rapid proliferation of cancer cells, facilitates the formation of tumor microenvironment to support their growth and survival, and impairs the efficacy of various tumor therapies. Therefore, the elucidation of the mechanisms driving aerobic glycolysis in tumors represents a pivotal breakthrough in developing therapeutic strategies for solid tumors. HIF1α serves as a central regulator of aerobic glycolysis with elevated mRNA and protein expression across multiple tumor types. However, the mechanisms contributing to this upregulation remain elusive. This study reports the identification of a novel HIF1α super enhancer (HSE) in multiple cancer cells using bioinformatics analysis, chromosome conformation capture (3C), chromatin immunoprecipitation (ChIP), and CRISPR/Cas9 genome editing techniques. Deletion of HSE in cancer cells significantly reduces the expression of HIF1α, glycolysis, cell proliferation, colony and tumor formation ability, confirming the role of HSE as the enhancer of HIF1α in cancer cells. Particularly, we demonstrated that STAT3 promotes the expression of HIF1α by binding to HSE. The discovery of HSE will help elucidate the pathways driving tumor aerobic glycolysis, offering new therapeutic targets and potentially resolving the bottleneck in solid tumor treatment.

HPV16-miRNAs exert oncogenic effects through enhancers in human cervical cancer

BACKGROUND

Cervical cancer is a human papillomavirus (HPV)-related disease. HPV type 16 (HPV16), which is the predominant cause of cervical cancer, can encode miRNAs (HPV16-miRNAs). However, the role of HPV16-miRNAs in the pathogenesis of cervical cancer remains unclear.

METHODS

Human cervical cancer cell lines SiHa (HPV16-positive) and C33A (HPV-negative), and cervical cancer tissues were collected to investigate the expression levels of two HPV16-miRNAs (HPV16-miR-H1 and HPV16-miR-H6). The overexpression and knockdown of HPV16-miR-H1 and HPV16-miR-H6 were performed using the lentiviral vector system and miRNA inhibitors, respectively. RNA-sequencing (RNA-seq) analysis and H3K27ac chromatin immunoprecipitation and sequencing (CHIP-seq) experiments were utilized to explore the roles of HPV16-miR-H1 and HPV16-miR-H6 facilitated by enhancers. CCK8, EdU, transwell, and wound healing assays were performed to verify the effects of HPV16-miR-H1 and HPV16-miR-H6 on cell proliferation and migration.

RESULTS

HPV16-miR-H1 and HPV16-miR-H6 were highly expressed in both SiHa cells and tissue samples from HPV16-positive cervical cancer patients. RNA-seq analysis showed that HPV16-miR-H1 and HPV16-miR-H6 induced the upregulation of numerous tumor progression-associated genes. H3K27ac CHIP-seq experiments further revealed that HPV16-miR-H1 and HPV16-miR-H6 modulated the expression of critical genes by regulating their enhancer activity. The functional study demonstrated that HPV16-miR-H1 and HPV16-miR-H6 increased the migratory capacity of SiHa cells.

CONCLUSIONS

Our data shed light on the role of HPV16-encoded miRNAs in cervical cancer, particularly emphasizing their involvement in the miRNA-enhancer-target gene system. This novel regulatory mechanism of HPV16-miRNAs provides new insights and approaches for the development of therapeutic strategies by targeting HPV16-positive cervical cancer.

Mechanosensitive super-enhancers regulate genes linked to atherosclerosis in endothelial cells

Vascular homeostasis and pathophysiology are tightly regulated by mechanical forces generated by hemodynamics. Vascular disorders such as atherosclerotic diseases largely occur at curvatures and bifurcations where disturbed blood flow activates endothelial cells while unidirectional flow at the straight part of vessels promotes endothelial health. Integrated analysis of the endothelial transcriptome, the 3D epigenome, and human genetics systematically identified the SNP-enriched cistrome in vascular endothelium subjected to well-defined atherosclerosis-prone disturbed flow or atherosclerosis-protective unidirectional flow. Our results characterized the endothelial typical- and super-enhancers and underscored the critical regulatory role of flow-sensitive endothelial super-enhancers. CRISPR interference and activation validated the function of a previously unrecognized unidirectional flow-induced super-enhancer that upregulates antioxidant genes NQO1, CYB5B, and WWP2, and a disturbed flow-induced super-enhancer in endothelium which drives prothrombotic genes EDN1 and HIVEP in vascular endothelium. Our results employing multiomics identify the cis-regulatory architecture of the flow-sensitive endothelial epigenome related to atherosclerosis and highlight the regulatory role of super-enhancers in mechanotransduction mechanisms.

PorcineAI-Enhancer: Prediction of Pig Enhancer Sequences Using Convolutional Neural Networks

Understanding the mechanisms of gene expression regulation is crucial in animal breeding. Cis-regulatory DNA sequences, such as enhancers, play a key role in regulating gene expression. Identifying enhancers is challenging, despite the use of experimental techniques and computational methods. Enhancer prediction in the pig genome is particularly significant due to the costliness of high-throughput experimental techniques. The study constructed a high-quality database of pig enhancers by integrating information from multiple sources. A deep learning prediction framework called PorcineAI-enhancer was developed for the prediction of pig enhancers. This framework employs convolutional neural networks for feature extraction and classification. PorcineAI-enhancer showed excellent performance in predicting pig enhancers, validated on an independent test dataset. The model demonstrated reliable prediction capability for unknown enhancer sequences and performed remarkably well on tissue-specific enhancer sequences.The study developed a deep learning prediction framework, PorcineAI-enhancer, for predicting pig enhancers. The model demonstrated significant predictive performance and potential for tissue-specific enhancers. This research provides valuable resources for future studies on gene expression regulation in pigs.

Identification of the enhancer RNAs related to tumorgenesis of pituitary neuroendocrine tumors

Background

Pituitary neuroendocrine tumors (PitNETs), which originate from the pituitary gland, account for 10%-15% of all intracranial neoplasms. Recent studies have indicated that enhancer RNAs (eRNAs) exert regulatory effects on tumor growth. However, the mechanisms underlying the eRNA-mediated tumorigenesis of PitNETs have not been elucidated.

Methods

Normal pituitary and PitNETs tissues were used to identify the differentially expressed eRNAs (DEEs). Immune gene sets and hallmarks of cancer gene sets were quantified based on single sample gene set enrichment analysis (ssGSEA) algorithm using GSVA. The perspective of immune cells among all samples was calculated by the CIBERSORT algorithm. Moreover, the regulatory network composed of key DEEs, target genes of eRNAs, hallmarks of cancer gene sets, differentially expressed TF, immune cells and immune gene sets were constructed by Pearson correlation analysis. Small molecular anti-PitNETs drugs were explored by CMap analysis and the accuracy of the study was verified by in vitro and in vivo experiments, ATAC-seq and ChIP-seq.

Results

In this study, data of 134 PitNETs and 107 non-tumorous pituitary samples were retrieved from a public database to identify differentially expressed genes. In total, 1128 differentially expressed eRNAs (DEEs) (494 upregulated eRNAs and 634 downregulated eRNAs) were identified. Next, the correlation of DEEs with cancer-related and immune-related gene signatures was examined to establish a co-expression regulatory network comprising 18 DEEs, 50 potential target genes of DEEs, 5 cancer hallmark gene sets, 2 differentially expressed transcription factors, 4 immune cell types, and 4 immune gene sets. Based on this network, the following four therapeutics for PitNETs were identified using Connectivity Map analysis: ciclopirox, bepridil, clomipramine, and alexidine. The growth-inhibitory effects of these therapeutics were validated using in vitro experiments. Ciclopirox exerted potential growth-inhibitory effects on PitNETs. Among the DEEs, GNLY, HOXB7, MRPL33, PRDM16, TCF7, and ZNF26 were determined to be potential diagnostic and therapeutic biomarkers for PitNETs.

Conclusion

This study illustrated the significant influence of eRNAs on the occurrence and development of PitNETs. By constructing the co-expression regulation network, GNLY, HOXB6, MRPL33, PRDM16, TCF7, and ZNF26 were identified as relatively significant DEEs which were considered as the novel biomarkers of diagnosis and treatment of PitNETs. This study demonstrated the roles of eRNAs in the occurrence and development of PitNETs and revealed that ciclopirox was a potential therapeutic for pituitary adenomas.

"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes

Regulatory elements activate promoters by recruiting transcription factors (TFs) to specific motifs. Notably, TF-DNA interactions often depend on cooperativity with colocalized partners, suggesting an underlying cis-regulatory syntax. To explore TF cooperativity in mammals, we analyze ∼500 mouse and human primary cells by combining an atlas of TF motifs, footprints, ChIP-seq, transcriptomes, and accessibility. We uncover two TF groups that colocalize with most expressed factors, forming stripes in hierarchical clustering maps. The first group includes lineage-determining factors that occupy DNA elements broadly, consistent with their key role in tissue-specific transcription. The second one, dubbed universal stripe factors (USFs), comprises ∼30 SP, KLF, EGR, and ZBTB family members that recognize overlapping GC-rich sequences in all tissues analyzed. Knockouts and single-molecule tracking reveal that USFs impart accessibility to colocalized partners and increase their residence time. Mammalian cells have thus evolved a TF superfamily with overlapping DNA binding that facilitate chromatin accessibility.

Prognostic Ability of Enhancer RNAs in Metastasis of Non-Small Cell Lung Cancer

(1) Background: Non-small cell lung cancer (NSCLC) is the most common lung cancer. Enhancer RNA (eRNA) has potential utility in the diagnosis, prognosis and treatment of cancer, but the role of eRNAs in NSCLC metastasis is not clear; (2) Methods: Differentially expressed transcription factors (DETFs), enhancer RNAs (DEEs), and target genes (DETGs) between primary NSCLC and metastatic NSCLC were identified. Prognostic DEEs (PDEEs) were screened by Cox regression analyses and a predicting model for metastatic NSCLC was constructed. We identified DEE interactions with DETFs, DETGs, reverse phase protein arrays (RPPA) protein chips, immunocytes, and pathways to construct a regulation network using Pearson correlation. Finally, the mechanisms and clinical significance were explained using multi-dimensional validation unambiguously; (3) Results: A total of 255 DEEs were identified, and 24 PDEEs were selected into the multivariate Cox regression model (AUC = 0.699). Additionally, the NSCLC metastasis-specific regulation network was constructed, and six key PDEEs were defined (ANXA8L1, CASTOR2, CYP4B1, GTF2H2C, PSMF1 and TNS4); (4) Conclusions: This study focused on the exploration of the prognostic value of eRNAs in the metastasis of NSCLC. Finally, six eRNAs were identified as potential markers for the prediction of metastasis of NSCLC.

Enhancer RNA: What we know and what we can achieve

Enhancers are important cis-acting elements that can regulate gene transcription and cell fate alongside promoters. In fact, many human cancers and diseases are associated with the malfunction of enhancers. Recent studies have shown that enhancers can produce enhancer RNAs (eRNAs) by RNA polymerase II. In this review, we discuss eRNA production, characteristics, functions and mechanics. eRNAs can determine chromatin accessibility, histone modification and gene expression by constructing a 'chromatin loop', thereby bringing enhancers to their target gene. eRNA can also be involved in the phase separation with enhancers and other proteins. eRNAs are abundant, and importantly, tissue-specific in tumours, various diseases and stem cells; thus, eRNAs can be a potential target for disease diagnosis and treatment. As eRNA is produced from the active transcription of enhancers and is involved in the regulation of cell fate, its manipulation will influence cell function, and therefore, it can be a new target for biological therapy.

CT-FOCS: a novel method for inferring cell type-specific enhancer–promoter maps

Spatiotemporal gene expression patterns are governed to a large extent by the activity of enhancer elements, which engage in physical contacts with their target genes. Identification of enhancer–promoter (EP) links that are functional only in a specific subset of cell types is a key challenge in understanding gene regulation. We introduce CT-FOCS (cell type FOCS), a statistical inference method that uses linear mixed effect models to infer EP links that show marked activity only in a single or a small subset of cell types out of a large panel of probed cell types. Analyzing 808 samples from FANTOM5, covering 472 cell lines, primary cells and tissues, CT-FOCS inferred such EP links more accurately than recent state-of-the-art methods. Furthermore, we show that strictly cell type-specific EP links are very uncommon in the human genome.

A deep generative model for multi-view profiling of single-cell RNA-seq and ATAC-seq data

Here, we present a multi-modal deep generative model, the single-cell Multi-View Profiler (scMVP), which is designed for handling sequencing data that simultaneously measure gene expression and chromatin accessibility in the same cell, including SNARE-seq, sci-CAR, Paired-seq, SHARE-seq, and Multiome from 10X Genomics. scMVP generates common latent representations for dimensionality reduction, cell clustering, and developmental trajectory inference and generates separate imputations for differential analysis and cis-regulatory element identification. scMVP can help mitigate data sparsity issues with imputation and accurately identify cell groups for different joint profiling techniques with common latent embedding, and we demonstrate its advantages on several realistic datasets.

The Stochastic Genome and Its Role in Gene Expression

Mammalian genomes have distinct levels of spatial organization and structure that have been hypothesized to play important roles in transcription regulation. Although much has been learned about these architectural features with ensemble techniques, single-cell studies are showing a new universal trend: Genomes are stochastic and dynamic at every level of organization. Stochastic gene expression, on the other hand, has been studied for years. In this review, we probe whether there is a causative link between the two phenomena. We specifically discuss the functionality of chromatin state, topologically associating domains (TADs), and enhancer biology in light of their stochastic nature and their specific roles in stochastic gene expression. We highlight persistent fundamental questions in this area of research.

Identification and Functional Characterization of Two Noncoding RNAs Transcribed from Putative Active Enhancers in Hepatocellular Carcinoma

Enhancers have been conventionally perceived as cis-acting elements that provide binding sites for trans-acting factors. However, recent studies have shown that enhancers are transcribed and that these transcripts, called enhancer RNAs (eRNAs), have a regulatory function. Here, we identified putative eRNAs by profiling and determining the overlap between noncoding RNA expression loci and eRNA-associated histone marks such as H3K27ac and H3K4me1 in hepatocellular carcinoma (HCC) cell lines. Of the 132 HCC-derived noncoding RNAs, 74 overlapped with the eRNA loci defined by the FANTOM consortium, and 65 were located in the proximal regions of genes differentially expressed between normal and tumor tissues in TCGA dataset. Interestingly, knockdown of two selected putative eRNAs, THUMPD3-AS1 and LINC01572, led to downregulation of their target mRNAs and to a reduction in the proliferation and migration of HCC cells. Additionally, the expression of these two noncoding RNAs and target mRNAs was elevated in tumor samples in the TCGA dataset, and high expression was associated with poor survival of patients. Collectively, our study suggests that noncoding RNAs such as THUMPD3-AS1 and LINC01572 (i.e., putative eRNAs) can promote the transcription of genes involved in cell proliferation and differentiation and that the dysregulation of these noncoding RNAs can cause cancers such as HCC.

New insights into the proteins interacting with the promoters of silkworm fibroin genes

The silkworm, Bombyx mori, is a silk-producing insect that has contributed greatly to human society. The silk gland of B. mori is a specialized organ responsible for synthesizing silk fibroin and sericin proteins under control of numerous factors. However, which factors are involved in direct silk protein synthesis regulation remains largely unknown. We report the identification of promoter-interacting proteins (PIPs) necessary for the regulation of genes encoding fibroin proteins, including the fibroin heavy chain (fibH), fibroin light chain (fibL), and a 25-kD polypeptide protein (P25). In the fourth larval molting stage (M4) or day 5 fifth-instar larvae (L5D5), a total of 198, 292, and 247 or 330, 305, and 460 proteins interacting with the promoter region of fibH, fibL and P25, respectively, were identified from the posterior silk gland by DNA pull-down combined with mass spectrometry. Many PIPs were particularly involved in ribosome- and metabolism-related pathways. Additionally, 135 and 212 proteins were identified as common PIPs of fibH, fibL and P25 in M4 and L5D5, respectively. Among all PIPs, we identified 31 potential transcription factors, such as Y-box and poly A-binding proteins, which play roles in nucleotide binding, ATP binding, or protein folding. This study provides the first in-depth profile of proteins interacting with fibroin gene promoters and contributes to a better understanding of silk protein synthesis regulation.

Reactivation of tumour suppressor in breast cancer by enhancer switching through NamiRNA network

Dysfunction of Tumour Suppressor Genes (TSGs) is a common feature in carcinogenesis. Epigenetic abnormalities including DNA hypermethylation or aberrant histone modifications in promoter regions have been described for interpreting TSG inactivation. However, in many instances, how TSGs are silenced in tumours are largely unknown. Given that miRNA with low expression in tumours is another recognized signature, we hypothesize that low expression of miRNA may reduce the activity of TSG related enhancers and further lead to inactivation of TSG during cancer development. Here, we reported that low expression of miRNA in cancer as a recognized signature leads to loss of function of TSGs in breast cancer. In 157 paired breast cancer and adjacent normal samples, tumour suppressor gene GPER1 and miR-339 are both downregulated in Luminal A/B and Triple Negative Breast Cancer subtypes. Mechanistic investigations revealed that miR-339 upregulates GPER1 expression in breast cancer cells by switching on the GPER1 enhancer, which can be blocked by enhancer deletion through the CRISPR/Cas9 system. Collectively, our findings reveal novel mechanistic insights into TSG dysfunction in cancer development, and provide evidence that reactivation of TSG by enhancer switching may be a promising alternative strategy for clinical breast cancer treatment.

Epigenetic landscapes of intracranial aneurysm risk haplotypes implicate enhancer function of endothelial cells and fibroblasts in dysregulated gene expression

BACKGROUND

Genome-wide association studies have identified many single nucleotide polymorphisms (SNPs) associated with increased risk for intracranial aneurysm (IA). However, how such variants affect gene expression within IA is poorly understood. We used publicly-available ChIP-Seq data to study chromatin landscapes surrounding risk loci to determine whether IA-associated SNPs affect functional elements that regulate gene expression in cell types comprising IA tissue.

METHODS

We mapped 16 significant IA-associated SNPs to linkage disequilibrium (LD) blocks within human genome. Using ChIP-Seq data, we examined these regions for presence of H3K4me1, H3K27ac, and H3K9ac histone marks (typically associated with latent/active enhancers). This analysis was conducted in several cell types that are present in IA tissue (endothelial cells, smooth muscle cells, fibroblasts, macrophages, monocytes, neutrophils, T cells, B cells, NK cells). In cell types with significant histone enrichment, we used HiC data to investigate topologically associated domains (TADs) encompassing the LD blocks to identify genes that may be affected by IA-associated variants. Bioinformatics were performed to determine the biological significance of these genes. Genes within HiC-defined TADs were also compared to differentially expressed genes from RNA-seq/microarray studies of IA tissues.

RESULTS

We found that endothelial cells and fibroblasts, rather than smooth muscle or immune cells, have significant enrichment for enhancer marks on IA risk haplotypes (p < 0.05). Bioinformatics demonstrated that genes within TADs subsuming these regions are associated with structural extracellular matrix components and enzymatic activity. The majority of histone marked TADs (83% fibroblasts [IMR90], 77% HUVEC) encompassed at least one differentially expressed gene from IA tissue studies.

CONCLUSIONS

These findings provide evidence that genetic variants associated with IA risk act on endothelial cells and fibroblasts. There is strong circumstantial evidence that this may be mediated through altered enhancer function, as genes in TADs encompassing enhancer marks have also been shown to be differentially expressed in IA tissue. These genes are largely related to organization and regulation of the extracellular matrix. This study builds upon our previous (Poppenberg et al., BMC Med Genomics, 2019) by including a more diverse set of data from additional cell types and by identifying potential affected genes (i.e. those in TADs).

CellWalker integrates single-cell and bulk data to resolve regulatory elements across cell types in complex tissues

Single-cell and bulk genomics assays have complementary strengths and weaknesses, and alone neither strategy can fully capture regulatory elements across the diversity of cells in complex tissues. We present CellWalker, a method that integrates single-cell open chromatin (scATAC-seq) data with gene expression (RNA-seq) and other data types using a network model that simultaneously improves cell labeling in noisy scATAC-seq and annotates cell type-specific regulatory elements in bulk data. We demonstrate CellWalker’s robustness to sparse annotations and noise using simulations and combined RNA-seq and ATAC-seq in individual cells. We then apply CellWalker to the developing brain. We identify cells transitioning between transcriptional states, resolve regulatory elements to cell types, and observe that autism and other neurological traits can be mapped to specific cell types through their regulatory elements.

A comprehensive enhancer screen identifies TRAM2 as a key and novel mediator of YAP oncogenesis

BACKGROUND

Frequent activation of the co-transcriptional factor YAP is observed in a large number of solid tumors. Activated YAP associates with enhancer loci via TEAD4-DNA-binding protein and stimulates cancer aggressiveness. Although thousands of YAP/TEAD4 binding-sites are annotated, their functional importance is unknown. Here, we aim at further identification of enhancer elements that are required for YAP functions.

RESULTS

We first apply genome-wide ChIP profiling of YAP to systematically identify enhancers that are bound by YAP/TEAD4. Next, we implement a genetic approach to uncover functions of YAP/TEAD4-associated enhancers, demonstrate its robustness, and use it to reveal a network of enhancers required for YAP-mediated proliferation. We focus on EnhancerTRAM2, as its target gene TRAM2 shows the strongest expression-correlation with YAP activity in nearly all tumor types. Interestingly, TRAM2 phenocopies the YAP-induced cell proliferation, migration, and invasion phenotypes and correlates with poor patient survival. Mechanistically, we identify FSTL-1 as a major direct client of TRAM2 that is involved in these phenotypes. Thus, TRAM2 is a key novel mediator of YAP-induced oncogenic proliferation and cellular invasiveness.

CONCLUSIONS

YAP is a transcription co-factor that binds to thousands of enhancer loci and stimulates tumor aggressiveness. Using unbiased functional approaches, we dissect YAP enhancer network and characterize TRAM2 as a novel mediator of cellular proliferation, migration, and invasion. Our findings elucidate how YAP induces cancer aggressiveness and may assist diagnosis of cancer metastasis.

The regulatory genome of the malaria vector Anopheles gambiae: integrating chromatin accessibility and gene expression

Anopheles gambiae mosquitoes are primary human malaria vectors, but we know very little about their mechanisms of transcriptional regulation. We profiled chromatin accessibility by the assay for transposase-accessible chromatin by sequencing (ATAC-seq) in laboratory-reared A. gambiae mosquitoes experimentally infected with the human malaria parasite Plasmodium falciparum. By integrating ATAC-seq, RNA-seq and ChIP-seq data, we showed a positive correlation between accessibility at promoters and introns, gene expression and active histone marks. By comparing expression and chromatin structure patterns in different tissues, we were able to infer cis-regulatory elements controlling tissue-specific gene expression and to predict the in vivo binding sites of relevant transcription factors. The ATAC-seq assay also allowed the precise mapping of active regulatory regions, including novel transcription start sites and enhancers that were annotated to mosquito immune-related genes. Not only is this study important for advancing our understanding of mechanisms of transcriptional regulation in the mosquito vector of human malaria, but the information we produced also has great potential for developing new mosquito-control and anti-malaria strategies.

Sex-Specific Alterations in Cardiac DNA Methylation in Adult Mice by Perinatal Lead Exposure

Environmental factors play an important role in the etiology of cardiovascular diseases. Cardiovascular diseases exhibit marked sexual dimorphism; however, the sex-specific effects of environmental exposures on cardiac health are incompletely understood. Perinatal and adult exposures to the metal lead (Pb) are linked to several adverse cardiovascular outcomes, but the sex-specific effects of this toxicant on the heart have received little attention. Perinatal environmental exposures can lead to disease through disruption of the normal epigenetic programming that occurs during early development. Using a mouse model of human-relevant perinatal environmental exposure, we investigated the effects of exposure to Pb during gestation and lactation on DNA methylation in the hearts of adult offspring mice (n = 6 per sex). Two weeks prior to mating, dams were assigned to control or Pb acetate (32 ppm) water, and exposure continued until offspring were weaned at three weeks of age. Enhanced reduced-representation bisulfite sequencing was used to measure DNA methylation in the hearts of offspring at five months of age. Although Pb exposure stopped at three weeks of age, we discovered hundreds of differentially methylated cytosines (DMCs) and regions (DMRs) in males and females at five months of age. DMCs/DMRs and their associated genes were sex-specific, with a small, but statistically significant subset overlapping between sexes. Pathway analysis revealed altered methylation of genes important for cardiac and other tissue development in males, and histone demethylation in females. Together, these data demonstrate that perinatal exposure to Pb induces sex-specific changes in cardiac DNA methylation that are present long after cessation of exposure, and highlight the importance of considering sex in environmental epigenetics and mechanistic toxicology studies.

Generation of Transgenic Fluorescent Reporter Lines for Studying Hematopoietic Development in the Mouse

Hematopoiesis in the mouse and other mammals occurs in several waves and arises from distinct anatomic sites. Transgenic mice expressing fluorescent reporter proteins at various points in the hematopoietic hierarchy, from hematopoietic stem cell to more restricted progenitors to each of the final differentiated cell types, have provided valuable tools for tagging, tracking, and isolating these cells. In this chapter, we discuss general considerations in designing a transgene, survey available fluorescent probes, and describe methods for confirming and analyzing transgene expression in the hematopoietic tissues of the embryo, fetus, and postnatal/adult animal.

Transcriptionally active enhancers in human cancer cells

The growth of human cancer cells is driven by aberrant enhancer and gene transcription activity. Here, we use transient transcriptome sequencing (TT-seq) to map thousands of transcriptionally active putative enhancers in fourteen human cancer cell lines covering seven types of cancer. These enhancers were associated with cell type-specific gene expression, enriched for genetic variants that predispose to cancer, and included functionally verified enhancers. Enhancer-promoter (E-P) pairing by correlation of transcription activity revealed ~ 40,000 putative E-P pairs, which were depleted for housekeeping genes and enriched for transcription factors, cancer-associated genes, and 3D conformational proximity. The cell type specificity and transcription activity of target genes increased with the number of paired putative enhancers. Our results represent a rich resource for future studies of gene regulation by enhancers and their role in driving cancerous cell growth.

Mediator's Kinase Module: A Modular Regulator of Cell Fate

Selective gene expression is crucial in maintaining the self-renewing and multipotent properties of stem cells. Mediator is a large, evolutionarily conserved, multi-subunit protein complex that modulates gene expression by relaying signals from cell type-specific transcription factors to RNA polymerase II. In humans, this complex consists of 30 subunits arranged in four modules. One critical module of the Mediator complex is the kinase module consisting of four subunits: MED12, MED13, CDK8, and CCNC. The kinase module exists in variable association with the 26-subunit Mediator core and affects transcription through phosphorylation of transcription factors and by controlling Mediator structure and function. Many studies have shown the kinase module to be a key player in the maintenance of stem cells that is distinct from a general role in transcription. Genetic studies have revealed that dysregulation of this kinase subunit contributes to the development of many human diseases. In this review, we discuss the importance of the Mediator kinase module by examining how this module functions with the more recently identified transcriptional super-enhancers, how changes in the kinase module and its activity can lead to the development of human disease, and the role of this unique module in directing and maintaining cell state. As we look to use stem cells to understand human development and treat human disease through both cell-based therapies and tissue engineering, we need to remain aware of the on-going research and address critical gaps in knowledge related to the molecular mechanisms that control cell fate.

Detection of CRISPR-mediated genome modifications through altered methylation patterns of CpG islands

BACKGROUND

The development and application of CRISPR technologies for the modification of the genome are rapidly expanding. Advances in the field describe new CRISPR components that are strategically engineered to improve the precision and reliability of CRISPR editing within the genome sequence. Genome modification using induced genome breaks that are targeted and mediated by CRISPR components leverage cellular mechanisms for repair like homology directed repair (HDR) to incorporate genomic edits with increased precision.

RESULTS

In this report, we describe the gain of methylation at typically hypomethylated CpG island (CGI) locations affected by the CRISPR-mediated incorporation of donor DNA using HDR mechanisms. With characterization of CpG methylation patterns using whole genome bisulfite sequencing, these CGI methylation disruptions trace the insertion of the donor DNA during the genomic edit. These insertions mediated by homology-directed recombination disrupt the generational methylation pattern stability of the edited CGI within the cells and their cellular lineage within the animal strain, persisting across generations. Our approach describes a statistically based workflow for indicating locations of modified CGIs and provides a mechanism for evaluating the directed modification of the methylome of the affected CGI at the CpG-level.

CONCLUSIONS

With advances in genome modification technology comes the need to detect the level and persistence of methylation change that modifications to the genomic sequence impose upon the collaterally edited methylome. Any modification of the methylome of somatic or germline cells could have implications for gene regulation mechanisms governed by the methylation patterns of CGI regions in the application of therapeutic edits of more sensitively regulated genomic regions. The method described here locates the directed modification of the mouse epigenome that persists over generations. While this observance would require supporting molecular observations such as direct sequence changes or gene expression changes, the observation of epigenetic modification provides an indicator that intentionally directed genomic edits can lead to collateral, unintentional epigenomic changes post modification with generational persistence.

Spirits in the Material World: Enhancer RNAs in Transcriptional Regulation

Responses to developmental and environmental cues depend on precise spatiotemporal control of gene transcription. Enhancers, which comprise DNA elements bound by regulatory proteins, can activate target genes in response to these external signals. Recent studies have shown that enhancers are transcribed to produce enhancer RNAs (eRNAs). Do eRNAs play a functional role in activating gene expression or are they non-functional byproducts of nearby transcription machinery? The unstable nature of eRNAs and over-reliance on knockdown approaches have made elucidating the possible functions of eRNAs challenging. We focus here on studies using cloned eRNAs to study their function as transcripts, revealing roles for eRNAs in enhancer-promoter looping, recruiting transcriptional machinery, and facilitating RNA polymerase pause-release to regulate gene expression.

TENET 2.0: Identification of key transcriptional regulators and enhancers in lung adenocarcinoma

Lung cancer is the leading cause of cancer-related death and lung adenocarcinoma is its most common subtype. Although genetic alterations have been identified as drivers in subsets of lung adenocarcinoma, they do not fully explain tumor development. Epigenetic alterations have been implicated in the pathogenesis of tumors. To identify epigenetic alterations driving lung adenocarcinoma, we used an improved version of the Tracing Enhancer Networks using Epigenetic Traits method (TENET 2.0) in primary normal lung and lung adenocarcinoma cells. We found over 32,000 enhancers that appear differentially activated between normal lung and lung adenocarcinoma. Among the identified transcriptional regulators inactivated in lung adenocarcinoma vs. normal lung, NKX2-1 was linked to a large number of silenced enhancers. Among the activated transcriptional regulators identified, CENPA, FOXM1, and MYBL2 were linked to numerous cancer-specific enhancers. High expression of CENPA, FOXM1, and MYBL2 is particularly observed in a subgroup of lung adenocarcinomas and is associated with poor patient survival. Notably, CENPA, FOXM1, and MYBL2 are also key regulators of cancer-specific enhancers in breast adenocarcinoma of the basal subtype, but they are associated with distinct sets of activated enhancers. We identified individual lung adenocarcinoma enhancers linked to CENPA, FOXM1, or MYBL2 that were associated with poor patient survival. Knockdown experiments of FOXM1 and MYBL2 suggest that these factors regulate genes involved in controlling cell cycle progression and cell division. For example, we found that expression of TK1, a potential target gene of a MYBL2-linked enhancer, is associated with poor patient survival. Identification and characterization of key transcriptional regulators and associated enhancers in lung adenocarcinoma provides important insights into the deregulation of lung adenocarcinoma epigenomes, highlighting novel potential targets for clinical intervention.

Multiomics data integration unveils core transcriptional regulatory networks governing cell-type identity

A plethora of computational approaches have been proposed for reconstructing gene regulatory networks (GRNs) from gene expression data. However, gene regulatory processes are often too complex to predict from the transcriptome alone. Here, we present a computational method, Moni, that systematically integrates epigenetics, transcriptomics, and protein–protein interactions to reconstruct GRNs among core transcription factors and their co-factors governing cell identity. We applied Moni to 57 datasets of human cell types and lines and demonstrate that it can accurately infer GRNs, thereby outperforming state-of-the-art methods.

Pervasive and CpG-dependent promoter-like characteristics of transcribed enhancers

The temporal and spatial expression of genes is controlled by promoters and enhancers. Findings obtained over the last decade that not only promoters but also enhancers are characterized by bidirectional, divergent transcription have challenged the traditional notion that promoters and enhancers represent distinct classes of regulatory elements. Over half of human promoters are associated with CpG islands (CGIs), relatively CpG-rich stretches of generally several hundred nucleotides that are often associated with housekeeping genes. Only about 6% of transcribed enhancers defined by CAGE-tag analysis are associated with CGIs. Here, we present an analysis of enhancer and promoter characteristics and relate them to the presence or absence of CGIs. We show that transcribed enhancers share a number of CGI-dependent characteristics with promoters, including statistically significant local overrepresentation of core promoter elements. CGI-associated enhancers are longer, display higher directionality of transcription, greater expression, a lesser degree of tissue specificity, and a higher frequency of transcription-factor binding events than non-CGI-associated enhancers. Genes putatively regulated by CGI-associated enhancers are enriched for transcription regulator activity. Our findings show that CGI-associated transcribed enhancers display a series of characteristics related to sequence, expression and function that distinguish them from enhancers not associated with CGIs.

Rare or Overlooked? Structural Disruption of Regulatory Domains in Human Neurocristopathies

In the last few years, the role of non-coding regulatory elements and their involvement in human disease have received great attention. Among the non-coding regulatory sequences, enhancers are particularly important for the proper establishment of cell type-specific gene-expression programs. Furthermore, the disruption of enhancers can lead to human disease through two main mechanisms: (i) Mutations or copy number variants can directly alter the enhancer sequences and thereby affect expression of their target genes; (ii) structural variants can provoke changes in 3-D chromatin organization that alter neither the enhancers nor their target genes, but rather the physical communication between them. In this review, these pathomechanisms are mostly discussed in the context of neurocristopathies, congenital disorders caused by defects that occur during neural crest development. We highlight why, due to its contribution to multiple tissues and organs, the neural crest represents an important, yet understudied, cell type involved in multiple congenital disorders. Moreover, we discuss currently available resources and experimental models for the study of human neurocristopathies. Last, we provide some practical guidelines that can be followed when investigating human neurocristopathies caused by structural variants. Importantly, these guidelines can be useful not only to uncover the etiology of human neurocristopathies, but also of other human congenital disorders in which enhancer disruption is involved.

Regulation of inflammation in diabetes: From genetics to epigenomics evidence

Background

Diabetes is one of the greatest public health challenges worldwide, and we still lack complementary approaches to significantly enhance the efficacy of preventive and therapeutic approaches. Genetic and environmental factors are the culprits involved in diabetes risk. Evidence from the last decade has highlighted that deregulation in the immune and inflammatory responses increase susceptibility to type 1 and type 2 diabetes. Spatiotemporal patterns of gene expression involved in immune cell polarisation depend on genomic enhancer elements in response to inflammatory and metabolic cues. Several studies have reported that most regulatory genetic variants are located in the non-protein coding regions of the genome and particularly in enhancer regions. The progress of high-throughput technologies has permitted the characterisation of enhancer chromatin properties. These advances support the concept that genetic alteration of enhancers may influence the immune and inflammatory responses in relation to diabetes.

Scope of review

Results from genome-wide association studies (GWAS) combined with functional and integrative analyses have elucidated the impacts of some diabetes risk-associated variants that are involved in the regulation of the immune system. Additionally, genetic variant mapping to enhancer regions may alter enhancer status, which in turn leads to aberrant expression of inflammatory genes associated with diabetes susceptibility. The focus of this review was to provide an overview of the current indications that inflammatory processes are regulated at the genetic and epigenomic levels in diabetes, along with perspectives on future research avenues that may improve understanding of the disease.

Major conclusions

In this review, we provide genetic evidence in support of a deregulated immune response as a risk factor in diabetes. We also argue about the importance of enhancer regions in the regulation of immune cell polarisation and how the recent advances using genome-wide methods for enhancer identification have enabled the determination of the impact of enhancer genetic variation on diabetes onset and phenotype. This could eventually lead to better management plans and improved treatment responses in human diabetes.

RegulationSpotter: annotation and interpretation of extratranscriptic DNA variants

RegulationSpotter is a web-based tool for the user-friendly annotation and interpretation of DNA variants located outside of protein-coding transcripts (extratranscriptic variants). It is designed for clinicians and researchers who wish to assess the potential impact of the considerable number of non-coding variants found in Whole Genome Sequencing runs. It annotates individual variants with underlying regulatory features in an intuitive way by assessing over 100 genome-wide annotations. Additionally, it calculates a score, which reflects the regulatory potential of the variant region. Its dichotomous classifications, ‘functional’ or ‘non-functional’, and a human-readable presentation of the underlying evidence allow a biologically meaningful interpretation of the score. The output shows key aspects of every variant and allows rapid access to more detailed information about its possible role in gene regulation. RegulationSpotter can either analyse single variants or complete VCF files. Variants located within protein-coding transcripts are automatically assessed by MutationTaster as well as by RegulationSpotter to account for possible intragenic regulatory effects. RegulationSpotter offers the possibility of using phenotypic data to focus on known disease genes or genomic elements interacting with them. RegulationSpotter is freely available at https://www.regulationspotter.org.

Expression vector cassette engineering for recombinant therapeutic production in mammalian cell systems

Human tissue plasminogen activator was the first recombinant therapy protein that successfully produced in Chinese hamster ovary cells in 1986 and approved for clinical use. Since then, more and more therapeutic proteins are being manufactured in mammalian cells, and the technologies for recombinant protein production in this expression system have developed rapidly, with the optimization of both upstream and downstream processes. One of the most promising strategies is expression vector cassette optimization based on the expression vector cassette. In this review paper, these approaches and developments are summarized, and the future strategy on the utilizing of expression cassettes for the production of recombinant therapeutic proteins in mammalian cells is discussed.

A CRISPR-Cas9 screen identifies essential CTCF anchor sites for estrogen receptor-driven breast cancer cell proliferation

Estrogen receptor α (ERα) is an enhancer activating transcription factor, a key driver of breast cancer and a main target for cancer therapy. ERα-mediated gene regulation requires proper chromatin-conformation to facilitate interactions between ERα-bound enhancers and their target promoters. A major determinant of chromatin structure is the CCCTC-binding factor (CTCF), that dimerizes and together with cohesin stabilizes chromatin loops and forms the boundaries of topologically associated domains. However, whether CTCF-binding elements (CBEs) are essential for ERα-driven cell proliferation is unknown. To address this question in a global manner, we implemented a CRISPR-based functional genetic screen targeting CBEs located in the vicinity of ERα-bound enhancers. We identified four functional CBEs and demonstrated the role of one of them in inducing chromatin conformation changes in favor of activation of PREX1, a key ERα target gene in breast cancer. Indeed, high PREX1 expression is a bona-fide marker of ERα-dependency in cell lines, and is associated with good outcome after anti-hormonal treatment. Altogether, our data show that distinct CTCF-mediated chromatin structures are required for ERα- driven breast cancer cell proliferation.

Identification of gene specific cis-regulatory elements during differentiation of mouse embryonic stem cells: An integrative approach using high-throughput datasets

Gene expression governs cell fate, and is regulated via a complex interplay of transcription factors and molecules that change chromatin structure. Advances in sequencing-based assays have enabled investigation of these processes genome-wide, leading to large datasets that combine information on the dynamics of gene expression, transcription factor binding and chromatin structure as cells differentiate. While numerous studies focus on the effects of these features on broader gene regulation, less work has been done on the mechanisms of gene-specific transcriptional control. In this study, we have focussed on the latter by integrating gene expression data for the in vitro differentiation of murine ES cells to macrophages and cardiomyocytes, with dynamic data on chromatin structure, epigenetics and transcription factor binding. Combining a novel strategy to identify communities of related control elements with a penalized regression approach, we developed individual models to identify the potential control elements predictive of the expression of each gene. Our models were compared to an existing method and evaluated using the existing literature and new experimental data from embryonic stem cell differentiation reporter assays. Our method is able to identify transcriptional control elements in a gene specific manner that reflect known regulatory relationships and to generate useful hypotheses for further testing.

Genetic Variation in Long-Range Enhancers

Cis-regulatory elements (CREs), including insulators, promoters, and enhancers, play critical roles in the establishment and maintenance of normal cellular function. Within each cell, the 3D structure of chromatin is arranged in specific patterns to expose the CREs required for optimal spatiotemporal regulation of gene expression. CREs can act over large distances along the linear genome, facilitated by looping of the intervening chromatin to allow direct interaction between distal regulatory elements and their target genes. A number of pathologies are associated with dysregulation of CRE function, including developmental disorders, cancers, and neuropsychiatric disease. A majority of known neuropsychiatric disease risk loci are noncoding, and increasing evidence suggests that they contribute to disease through disruption of CREs. As such, rather than directly altering the amino acid content of proteins, these variants are instead thought to affect where, when, and to what extent a given gene is expressed. The distances over which CREs can operate often render their target genes difficult to identify. Furthermore, as many risk loci contain multiple variants in high linkage disequilibrium, identification of the causative single nucleotide polymorphism(s) therein is not straightforward. Thus, deciphering the genetic etiology of complex neuropsychiatric disorders presents a significant challenge.

Insights into active intragenic enhancers

Enhancers can regulate gene transcription from afar. Many enhancers are located in genes. Although the regulatory roles of several individual intragenic enhancers have been elaborated, a genome-wide insight into intragenic enhancers remains to be elucidated. We found that active intragenic enhancers have a preference for being located in expressed genes. Unlike intergenic enhancers, active intragenic enhancers are enriched of H3K79me2 epigenetic signal, and depleted of variant histone H2A.Z. Moreover, eRNAs of active intragenic enhancers show lower degradation rates than those of the other enhancers. Our findings will have implications in understanding functions of intragenic enhancers.

Identifying lncRNA-mediated regulatory modules via ChIA-PET network analysis

Background

Although several studies have provided insights into the role of long non-coding RNAs (lncRNAs), the majority of them have unknown function. Recent evidence has shown the importance of both lncRNAs and chromatin interactions in transcriptional regulation. Although network-based methods, mainly exploiting gene-lncRNA co-expression, have been applied to characterize lncRNA of unknown function by means of ’guilt-by-association’, no strategy exists so far which identifies mRNA-lncRNA functional modules based on the 3D chromatin interaction graph.

Results

To better understand the function of chromatin interactions in the context of lncRNA-mediated gene regulation, we have developed a multi-step graph analysis approach to examine the RNA polymerase II ChIA-PET chromatin interaction network in the K562 human cell line. We have annotated the network with gene and lncRNA coordinates, and chromatin states from the ENCODE project. We used centrality measures, as well as an adaptation of our previously developed Markov State Models (MSM) clustering method, to gain a better understanding of lncRNAs in transcriptional regulation. The novelty of our approach resides in the detection of fuzzy regulatory modules based on network properties and their optimization based on co-expression analysis between genes and gene-lncRNA pairs. This results in our method returning more bona fide regulatory modules than other state-of-the art approaches for clustering on graphs.

Conclusions

Interestingly, we find that lncRNA network hubs tend to be significantly enriched in evolutionary conserved lncRNAs and enhancer-like functions. We validated regulatory functions for well known lncRNAs, such as MALAT1 and the enhancer-like lncRNA FALEC. In addition, by investigating the modular structure of bigger components we mine putative regulatory functions for uncharacterized lncRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2900-8) contains supplementary material, which is available to authorized users.

Live-cell imaging reveals enhancer-dependent Sox2 transcription in the absence of enhancer proximity

Enhancers are important regulatory elements that can control gene activity across vast genetic distances. However, the underlying nature of this regulation remains obscured because it has been difficult to observe in living cells. Here, we visualize the spatial organization and transcriptional output of the key pluripotency regulator Sox2 and its essential enhancer Sox2 Control Region (SCR) in living embryonic stem cells (ESCs). We find that Sox2 and SCR show no evidence of enhanced spatial proximity and that spatial dynamics of this pair is limited over tens of minutes. Sox2 transcription occurs in short, intermittent bursts in ESCs and, intriguingly, we find this activity demonstrates no association with enhancer proximity, suggesting that direct enhancer-promoter contacts do not drive contemporaneous Sox2 transcription. Our study establishes a framework for interrogation of enhancer function in living cells and supports an unexpected mechanism for enhancer control of Sox2 expression that uncouples transcription from enhancer proximity.

Arid3a regulates nephric tubule regeneration via evolutionarily conserved regeneration signal-response enhancers

Amphibians and fish have the ability to regenerate numerous tissues, whereas mammals have a limited regenerative capacity. Despite numerous developmental genes becoming reactivated during regeneration, an extensive analysis is yet to be performed on whether highly regenerative animals utilize unique cis-regulatory elements for the reactivation of genes during regeneration and how such cis-regulatory elements become activated. Here, we screened regeneration signal-response enhancers at the lhx1 locus using Xenopus and found that the noncoding elements conserved from fish to human function as enhancers in the regenerating nephric tubules. A DNA-binding motif of Arid3a, a component of H3K9me3 demethylases, was commonly found in RSREs. Arid3a binds to RSREs and reduces the H3K9me3 levels. It promotes cell cycle progression and causes the outgrowth of nephric tubules, whereas the conditional knockdown of arid3a using photo-morpholino inhibits regeneration. These results suggest that Arid3a contributes to the regeneration of nephric tubules by decreasing H3K9me3 on RSREs.

Inferring the Molecular Mechanisms of Noncoding Alzheimer's Disease-Associated Genetic Variants

Most of the loci identified by genome-wide association studies (GWAS) for late-onset Alzheimer's disease (LOAD) are in strong linkage disequilibrium (LD) with nearby variants all of which could be the actual functional variants, often in non-protein-coding regions and implicating underlying gene regulatory mechanisms. We set out to characterize the causal variants, regulatory mechanisms, tissue contexts, and target genes underlying these associations. We applied our INFERNO algorithm to the top 19 non-APOE loci from the IGAP GWAS study. INFERNO annotated all LD-expanded variants at each locus with tissue-specific regulatory activity. Bayesian co-localization analysis of summary statistics and eQTL data was performed to identify tissue-specific target genes. INFERNO identified enhancer dysregulation in all 19 tag regions analyzed, significant enrichments of enhancer overlaps in the immune-related blood category, and co-localized eQTL signals overlapping enhancers from the matching tissue class in ten regions (ABCA7, BIN1, CASS4, CD2AP, CD33, CELF1, CLU, EPHA1, FERMT2, ZCWPW1). In several cases, we identified dysregulation of long noncoding RNA (lncRNA) transcripts and applied the lncRNA target identification algorithm from INFERNO to characterize their downstream biological effects. We also validated the allele-specific effects of several variants on enhancer function using luciferase expression assays. By integrating functional genomics with GWAS signals, our analysis yielded insights into the regulatory mechanisms, tissue contexts, genes, and biological processes affected by noncoding genetic variation associated with LOAD risk.

Differential analysis of chromatin accessibility and histone modifications for predicting mouse developmental enhancers

Enhancers are distal cis-regulatory elements that modulate gene expression. They are depleted of nucleosomes and enriched in specific histone modifications; thus, calling DNase-seq and histone mark ChIP-seq peaks can predict enhancers. We evaluated nine peak-calling algorithms for predicting enhancers validated by transgenic mouse assays. DNase and H3K27ac peaks were consistently more predictive than H3K4me1/2/3 and H3K9ac peaks. DFilter and Hotspot2 were the best DNase peak callers, while HOMER, MUSIC, MACS2, DFilter and F-seq were the best H3K27ac peak callers. We observed that the differential DNase or H3K27ac signals between two distant tissues increased the area under the precision-recall curve (PR-AUC) of DNase peaks by 17.5-166.7% and that of H3K27ac peaks by 7.1-22.2%. We further improved this differential signal method using multiple contrast tissues. Evaluated using a blind test, the differential H3K27ac signal method substantially improved PR-AUC from 0.48 to 0.75 for predicting heart enhancers. We further validated our approach using postnatal retina and cerebral cortex enhancers identified by massively parallel reporter assays, and observed improvements for both tissues. In summary, we compared nine peak callers and devised a superior method for predicting tissue-specific mouse developmental enhancers by reranking the called peaks.

coTRaCTE predicts co-occurring transcription factors within cell-type specific enhancers

Cell-type specific gene expression is regulated by the combinatorial action of transcription factors (TFs). In this study, we predict transcription factor (TF) combinations that cooperatively bind in a cell-type specific manner. We first divide DNase hypersensitive sites into cell-type specifically open vs. ubiquitously open sites in 64 cell types to describe possible cell-type specific enhancers. Based on the pattern contrast between these two groups of sequences we develop "co-occurring TF predictor on Cell-Type specific Enhancers" (coTRaCTE) - a novel statistical method to determine regulatory TF co-occurrences. Contrasting the co-binding of TF pairs between cell-type specific and ubiquitously open chromatin guarantees the high cell-type specificity of the predictions. coTRaCTE predicts more than 2000 co-occurring TF pairs in 64 cell types. The large majority (70%) of these TF pairs is highly cell-type specific and overlaps in TF pair co-occurrence are highly consistent among related cell types. Furthermore, independently validated co-occurring and directly interacting TFs are significantly enriched in our predictions. Focusing on the regulatory network derived from the predicted co-occurring TF pairs in embryonic stem cells (ESCs) we find that it consists of three subnetworks with distinct functions: maintenance of pluripotency governed by OCT4, SOX2 and NANOG, regulation of early development governed by KLF4, STAT3, ZIC3 and ZNF148 and general functions governed by MYC, TCF3 and YY1. In summary, coTRaCTE predicts highly cell-type specific co-occurring TFs which reveal new insights into transcriptional regulatory mechanisms.

CRISPR/Cas9-Based Engineering of the Epigenome

Determining causal relationships between distinct chromatin features and gene expression, and ultimately cell behavior, remains a major challenge. Recent developments in targetable epigenome-editing tools enable us to assign direct transcriptional and functional consequences to locus-specific chromatin modifications. This Protocol Review discusses the unprecedented opportunity that CRISPR/Cas9 technology offers for investigating and manipulating the epigenome to facilitate further understanding of stem cell biology and engineering of stem cells for therapeutic applications. We also provide technical considerations for standardization and further improvement of the CRISPR/Cas9-based tools to engineer the epigenome.

Characterization of a morphogenetic furrow specific Gal4 driver in the developing Drosophila eye

The ability to express a gene of interest in a spatio-temporal manner using Gal4-UAS system has allowed the use of Drosophila model to study various biological phenomenon. During Drosophila eye development, a synchronous wave of differentiation called Morphogenetic furrow (MF) initiates at the posterior margin resulting in differentiation of retinal neurons. This synchronous differentiation is also observed in the differentiating retina of vertebrates. Since MF is highly dynamic, it can serve as an excellent model to study patterning and differentiation. However, there are not any Gal4 drivers available to observe the gain- of- function or loss- of- function of a gene specifically along the dynamic MF. The decapentaplegic (dpp) gene encodes a secreted protein of the transforming growth factor-beta (TGF-beta) superfamily that expresses at the posterior margin and then moves with the MF. However, unlike the MF associated pattern of dpp gene expression, the targeted dpp-Gal4 driver expression is restricted to the posterior margin of the developing eye disc. We screened GMR lines harboring regulatory regions of dpp fused with Gal4 coding region to identify MF specific enhancer of dpp using a GFP reporter gene. We employed immuno-histochemical approaches to detect gene expression. The rationale was that GFP reporter expression will correspond to the dpp expression domain in the developing eye. We identified two new dpp-Gal4 lines, viz., GMR17E04-Gal4 and GMR18D08-Gal4 that carry sequences from first intron region of dpp gene. GMR17E04-Gal4 drives expression along the MF during development and later in the entire pupal retina whereas GMR18D08-Gal4 drives expression of GFP transgene in the entire developing eye disc, which later drives expression only in the ventral half of the pupal retina. Thus, GMR18D08-Gal4 will serve as a new reagent for targeting gene expression in the ventral half of the pupal retina. We compared misexpression phenotypes of Wg, a negative regulator of eye development, using GMR17E04-Gal4, GMR18D08-Gal4 with existing dpp-Gal4 driver. The eye phenotypes generated by using our newly identified MF specific driver are not similar to the ones generated by existing dpp-Gal4 driver. It suggests that misexpression studies along MF needs revisiting using the new Gal4 drivers generated in our studies.

Differential tissue specific, temporal and spatial expression patterns of the Aggrecan gene is modulated by independent enhancer elements

The transcriptional mechanism through which chondrocytes control the spatial and temporal composition of the cartilage tissue has remained largely elusive. The central aim of this study was to identify whether transcriptional enhancers played a role in the organisation of the chondrocytes in cartilaginous tissue. We focused on the Aggrecan gene (Acan) as it is essential for the normal structure and function of cartilage and it is expressed developmentally in different stages of chondrocyte maturation. Using transgenic reporter studies in mice we identified four elements, two of which showed individual chondrocyte developmental stage specificity. In particular, one enhancer (-80) distinguishes itself from the others by being predominantly active in adult cartilage. Furthermore, the -62 element uniquely drove reporter activity in early chondrocytes. The remaining chondrocyte specific enhancers, +28 and -30, showed no preference to chondrocyte type. The transcription factor SOX9 interacted with all the enhancers in vitro and mutation of SOX9 binding sites in one of the enhancers (-30) resulted in a loss of its chondrocyte specificity and ectopic enhancer reporter activity. Thus, the Acan enhancers orchestrate the precise spatiotemporal expression of this gene in cartilage types at different stages of development and adulthood.

MuSERA: Multiple Sample Enriched Region Assessment

Enriched region (ER) identification is a fundamental step in several next-generation sequencing (NGS) experiment types. Yet, although NGS experimental protocols recommend producing replicate samples for each evaluated condition and their consistency is usually assessed, typically pipelines for ER identification do not consider available NGS replicates. This may alter genome-wide descriptions of ERs, hinder significance of subsequent analyses on detected ERs and eventually preclude biological discoveries that evidence in replicate could support. MuSERA is a broadly useful stand-alone tool for both interactive and batch analysis of combined evidence from ERs in multiple ChIP-seq or DNase-seq replicates. Besides rigorously combining sample replicates to increase statistical significance of detected ERs, it also provides quantitative evaluations and graphical features to assess the biological relevance of each determined ER set within its genomic context; they include genomic annotation of determined ERs, nearest ER distance distribution, global correlation assessment of ERs and an integrated genome browser. We review MuSERA rationale and implementation, and illustrate how sets of significant ERs are expanded by applying MuSERA on replicates for several types of NGS data, including ChIP-seq of transcription factors or histone marks and DNase-seq hypersensitive sites. We show that MuSERA can determine a new, enhanced set of ERs for each sample by locally combining evidence on replicates, and prove how the easy-to-use interactive graphical displays and quantitative evaluations that MuSERA provides effectively support thorough inspection of obtained results and evaluation of their biological content, facilitating their understanding and biological interpretations. MuSERA is freely available at http://www.bioinformatics.deib.polimi.it/MuSERA/.

Scc2/Nipbl hops between chromosomal cohesin rings after loading

The cohesin complex mediates DNA-DNA interactions both between (sister chromatid cohesion) and within chromosomes (DNA looping). It has been suggested that intra-chromosome loops are generated by extrusion of DNAs through the lumen of cohesin's ring. Scc2 (Nipbl) stimulates cohesin's ABC-like ATPase and is essential for loading cohesin onto chromosomes. However, it is possible that the stimulation of cohesin's ATPase by Scc2 also has a post-loading function, for example driving loop extrusion. Using fluorescence recovery after photobleaching (FRAP) and single-molecule tracking in human cells, we show that Scc2 binds dynamically to chromatin, principally through an association with cohesin. Scc2's movement within chromatin is consistent with a 'stop-and-go' or 'hopping' motion. We suggest that a low diffusion coefficient, a low stoichiometry relative to cohesin, and a high affinity for chromosomal cohesin enables Scc2 to move rapidly from one chromosomal cohesin complex to another, performing a function distinct from loading.

Single-cell profiling reveals that eRNA accumulation at enhancer-promoter loops is not required to sustain transcription

Enhancers are intergenic DNA elements that regulate the transcription of target genes in response to signaling pathways by interacting with promoters over large genomic distances. Recent studies have revealed that enhancers are bi-directionally transcribed into enhancer RNAs (eRNAs). Using single-molecule fluorescence in situ hybridization (smFISH), we investigated the eRNA-mediated regulation of transcription during estrogen induction in MCF-7 cells. We demonstrate that eRNAs are localized exclusively in the nucleus and are induced with similar kinetics as target mRNAs. However, eRNAs are mostly nascent at enhancers and their steady-state levels remain lower than those of their cognate mRNAs. Surprisingly, at the single-allele level, eRNAs are rarely co-expressed with their target loci, demonstrating that active gene transcription does not require the continuous transcription of eRNAs or their accumulation at enhancers. When co-expressed, sub-diffraction distance measurements between nascent mRNA and eRNA signals reveal that co-transcription of eRNAs and mRNAs rarely occurs within closed enhancer–promoter loops. Lastly, basal eRNA transcription at enhancers, but not E2-induced transcription, is maintained upon depletion of MLL1 and ERα, suggesting some degree of chromatin accessibility prior to signal-dependent activation of transcription. Together, our findings suggest that eRNA accumulation at enhancer–promoter loops is not required to sustain target gene transcription.

Genetic variation in the developmental regulation of cortical avpr1a among prairie voles

Early experiences can have enduring impacts on brain and behavior, but the strength of these effects can be influenced by genetic variation. In principle, polymorphic CpGs (polyCpGs) may contribute to gene-by-environment interactions (G × E) by altering DNA methylation. In this study, we investigate the influence of polyCpGs on the development of vasopressin receptor 1a abundance in the retrosplenial cortex (RSC-V1aR) of prairie voles (Microtus ochrogaster). Two alternative alleles ('HI'/'LO') predict RSC avpr1a expression, V1aR abundance and sexual fidelity in adulthood; these alleles differ in the frequency of CpG sites and in methylation at a putative intron enhancer. We hypothesized that the elevated CpG abundance in the LO allele would make homozygous LO/LO voles more sensitive to developmental perturbations. We found that genotype differences in RSC-V1aR abundance emerged early in ontogeny and were accompanied by differences in methylation of the putative enhancer. As predicted, postnatal treatment with an oxytocin receptor antagonist (OTA) reduced RSC-V1aR abundance in LO/LO adults but not their HI/HI siblings. Similarly, methylation inhibition by zebularine increased RSC-V1aR in LO/LO adults, but not in HI/HI siblings. These data show a gene-by-environment interaction in RSC-V1aR. Surprisingly, however, neither OTA nor zebularine altered adult methylation of the intronic enhancer, suggesting that differences in sensitivity could not be explained by CpG density at the enhancer alone. Methylated DNA immunoprecipiation-sequencing showed additional differentially methylated regions between HI/HI and LO/LO voles. Future research should examine the role of these regions and other regulatory elements in the ontogeny of RSC-V1aR and its developmentally induced changes.

Systematic tissue-specific functional annotation of the human genome highlights immune-related DNA elements for late-onset Alzheimer's disease

Continuing efforts from large international consortia have made genome-wide epigenomic and transcriptomic annotation data publicly available for a variety of cell and tissue types. However, synthesis of these datasets into effective summary metrics to characterize the functional non-coding genome remains a challenge. Here, we present GenoSkyline-Plus, an extension of our previous work through integration of an expanded set of epigenomic and transcriptomic annotations to produce high-resolution, single tissue annotations. After validating our annotations with a catalog of tissue-specific non-coding elements previously identified in the literature, we apply our method using data from 127 different cell and tissue types to present an atlas of heritability enrichment across 45 different GWAS traits. We show that broader organ system categories (e.g. immune system) increase statistical power in identifying biologically relevant tissue types for complex diseases while annotations of individual cell types (e.g. monocytes or B-cells) provide deeper insights into disease etiology. Additionally, we use our GenoSkyline-Plus annotations in an in-depth case study of late-onset Alzheimer's disease (LOAD). Our analyses suggest a strong connection between LOAD heritability and genetic variants contained in regions of the genome functional in monocytes. Furthermore, we show that LOAD shares a similar localization of SNPs to monocyte-functional regions with Parkinson's disease. Overall, we demonstrate that integrated genome annotations at the single tissue level provide a valuable tool for understanding the etiology of complex human diseases. Our GenoSkyline-Plus annotations are freely available at http://genocanyon.med.yale.edu/GenoSkyline.

Global reorganisation of cis-regulatory units upon lineage commitment of human embryonic stem cells

Long-range cis-regulatory elements such as enhancers coordinate cell-specific transcriptional programmes by engaging in DNA looping interactions with target promoters. Deciphering the interplay between the promoter connectivity and activity of cis-regulatory elements during lineage commitment is crucial for understanding developmental transcriptional control. Here, we use Promoter Capture Hi-C to generate a high-resolution atlas of chromosomal interactions involving ~22,000 gene promoters in human pluripotent and lineage-committed cells, identifying putative target genes for known and predicted enhancer elements. We reveal extensive dynamics of cis-regulatory contacts upon lineage commitment, including the acquisition and loss of promoter interactions. This spatial rewiring occurs preferentially with predicted changes in the activity of cis-regulatory elements and is associated with changes in target gene expression. Our results provide a global and integrated view of promoter interactome dynamics during lineage commitment of human pluripotent cells.