Functional Cis-regulatory Elements Research Articles

Quantitative analysis of cis-regulatory elements in transcription with KAS-ATAC-seq

Cis-regulatory elements (CREs) are pivotal in orchestrating gene expression throughout diverse biological systems. Accurate identification and in-depth characterization of functional CREs are crucial for decoding gene regulation networks during cellular processes. In this study, we develop Kethoxal-Assisted Single-stranded DNA Assay for Transposase-Accessible Chromatin with Sequencing (KAS-ATAC-seq) to quantitatively analyze the transcriptional activity of CREs. A main advantage of KAS-ATAC-seq lies in its precise measurement of ssDNA levels within both proximal and distal ATAC-seq peaks, enabling the identification of transcriptional regulatory sequences. This feature is particularly adept at defining Single-Stranded Transcribing Enhancers (SSTEs). SSTEs are highly enriched with nascent RNAs and specific transcription factors (TFs) binding sites that define cellular identity. Moreover, KAS-ATAC-seq provides a detailed characterization and functional implications of various SSTE subtypes. Our analysis of CREs during mouse neural differentiation demonstrates that KAS-ATAC-seq can effectively identify immediate-early activated CREs in response to retinoic acid (RA) treatment. Our findings indicate that KAS-ATAC-seq provides more precise annotation of functional CREs in transcription. Future applications of KAS-ATAC-seq would help elucidate the intricate dynamics of gene regulation in diverse biological processes.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq–based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

Genome-Wide Identification of Open Chromatin in Plants Using MH-Seq.

Functional cis-regulatory elements (CREs) act as precise transcriptional switches for fine-tuning gene transcription. Identification of CREs is critical for understanding regulatory mechanisms of gene expression associated with various biological processes in eukaryotes. It is well known that CREs reside in open chromatin that exhibits hypersensitivity to enzyme cleavage and physical shearing. Currently, high-throughput methodologies, such as DNase-seq, ATAC-seq, and FAIRE-seq, have been widely applied in mapping open chromatin in various eukaryotic genomes. More recently, differential MNase (micrococcal nuclease) treatment has been successfully employed to map open chromatin in addition to profiling nucleosome landscape in both mammalian and plant species. We have developed a MNase hypersensitivity sequencing (MH-seq) technique in plants. The MH-seq procedure includes plant nuclei fixation and purification, differential treatments of purified nuclei with MNase, specific recovery of MNase-trimmed small DNA fragments within 20~100bp in length, and MH-seq library construction followed by Illumina sequencing and data analysis. MH-seq has been successfully applied for global identification of open chromatin in both Arabidopsis thaliana and maize. It has been proven to be an attractive alternative for profiling open chromatin. Thus, MH-seq is expected to be valuable in probing chromatin accessibility on a genome-wide scale for other plants with sequenced genomes. Moreover, MHS data allow to implement footprinting assays to unveil binding sites of transcription factors.

Genomic features underlie the co-option of SVA transposons as cis-regulatory elements in human pluripotent stem cells.

Domestication of transposable elements (TEs) into functional cis-regulatory elements is a widespread phenomenon. However, the mechanisms behind why some TEs are co-opted as functional enhancers while others are not are underappreciated. SINE-VNTR-Alus (SVAs) are the youngest group of transposons in the human genome, where ~3,700 copies are annotated, nearly half of which are human-specific. Many studies indicate that SVAs are among the most frequently co-opted TEs in human gene regulation, but the mechanisms underlying such processes have not yet been thoroughly investigated. Here, we leveraged CRISPR-interference (CRISPRi), computational and functional genomics to elucidate the genomic features that underlie SVA domestication into human stem-cell gene regulation. We found that ~750 SVAs are co-opted as functional cis-regulatory elements in human induced pluripotent stem cells. These SVAs are significantly closer to genes and harbor more transcription factor binding sites than non-co-opted SVAs. We show that a long DNA motif composed of flanking YY1/2 and OCT4 binding sites is enriched in the co-opted SVAs and that these two transcription factors bind consecutively on the TE sequence. We used CRISPRi to epigenetically repress active SVAs in stem cell-like NCCIT cells. Epigenetic perturbation of active SVAs strongly attenuated YY1/OCT4 binding and influenced neighboring gene expression. Ultimately, SVA repression resulted in ~3,000 differentially expressed genes, 131 of which were the nearest gene to an annotated SVA. In summary, we demonstrated that SVAs modulate human gene expression, and uncovered that location and sequence composition contribute to SVA domestication into gene regulatory networks.

Prediction and Experimental Validation of a New Salinity-Responsive Cis-Regulatory Element (CRE) in a Tilapia Cell Line.

Transcriptional regulation is a major mechanism by which organisms integrate gene x environment interactions. It can be achieved by coordinated interplay between cis-regulatory elements (CREs) and transcription factors (TFs). Euryhaline tilapia (Oreochromis mossambicus) tolerate a wide range of salinity and thus are an appropriate model to examine transcriptional regulatory mechanisms during salinity stress in fish. Quantitative proteomics in combination with the transcription inhibitor actinomycin D revealed 19 proteins that are transcriptionally upregulated by hyperosmolality in tilapia brain (OmB) cells. We searched the extended proximal promoter up to intron1 of each corresponding gene for common motifs using motif discovery tools. The top-ranked motif identified (STREME1) represents a binding site for the Forkhead box TF L1 (FoxL1). STREME1 function during hyperosmolality was experimentally validated by choosing two of the 19 genes, chloride intracellular channel 2 (clic2) and uridine phosphorylase 1 (upp1), that are enriched in STREME1 in their extended promoters. Transcriptional induction of these genes during hyperosmolality requires STREME1, as evidenced by motif mutagenesis. We conclude that STREME1 represents a new functional CRE that contributes to gene x environment interactions during salinity stress in tilapia. Moreover, our results indicate that FoxL1 family TFs are contribute to hyperosmotic induction of genes in euryhaline fish.

Investigating chromatin accessibility during development and differentiation by ATAC-sequencing to guide the identification of cis-regulatory elements.

Mapping accessible chromatin across time scales can give insights into its dynamic nature, for example during cellular differentiation and tissue or organism development. Analysis of such data can be utilised to identify functional cis-regulatory elements (CRE) and transcription factor binding sites and, when combined with transcriptomics, can reveal gene regulatory networks (GRNs) of expressed genes. Chromatin accessibility mapping is a powerful approach and can be performed using ATAC-sequencing (ATAC-seq), whereby Tn5 transposase inserts sequencing adaptors into genomic DNA to identify differentially accessible regions of chromatin in different cell populations. It requires low sample input and can be performed and analysed relatively quickly compared with other methods. The data generated from ATAC-seq, along with other genomic approaches, can help uncover chromatin packaging and potential cis-regulatory elements that may be responsible for gene expression. Here, we describe the ATAC-seq approach and give examples from mainly vertebrate embryonic development, where such datasets have identified the highly dynamic nature of chromatin, with differing landscapes between cellular precursors for different lineages.

Structural and functional analysis of the human cone-rod homeobox transcription factor.

The cone‐rod homeobox (CRX) protein is a critical K50 homeodomain transcription factor responsible for the differentiation and maintenance of photoreceptor neurons in the vertebrate retina. Mutant alleles in the human gene encoding CRX result in a variety of distinct blinding retinopathies, including retinitis pigmentosa, cone‐rod dystrophy, and Leber congenital amaurosis. Despite the success of using in vitro biochemistry, animal models, and genomics approaches to study this clinically relevant transcription factor over the past 25 years since its initial characterization, there are no high‐resolution structures in the published literature for the CRX protein. In this study, we use bioinformatic approaches and small‐angle X‐ray scattering (SAXS) structural analysis to further understand the biochemical complexity of the human CRX homeodomain (CRX‐HD). We find that the CRX‐HD is a compact, globular monomer in solution that can specifically bind functional cis‐regulatory elements encoded upstream of retina‐specific genes. This study presents the first structural analysis of CRX, paving the way for a new approach to studying the biochemistry of this protein and its disease‐causing mutant protein variants.

Coupling high-throughput mapping with proteomics analysis delineates cis-regulatory elements at high resolution.

Growing evidence suggests that functional cis-regulatory elements (cis-REs) not only exist in epigenetically marked but also in unmarked sites of the human genome. While it is already difficult to identify cis-REs in the epigenetically marked sites, interrogating cis-REs residing within the unmarked sites is even more challenging. Here, we report adapting Reel-seq, an in vitro high-throughput (HTP) technique, to fine-map cis-REs at high resolution over a large region of the human genome in a systematic and continuous manner. Using Reel-seq, as a proof-of-principle, we identified 408 candidate cis-REs by mapping a 58 kb core region on the aging-related CDKN2A/B locus that harbors p16INK4a. By coupling Reel-seq with FREP-MS, a proteomics analysis technique, we characterized two cis-REs, one in an epigenetically marked site and the other in an epigenetically unmarked site. These elements are shown to regulate the p16INK4a expression over an ∼100 kb distance by recruiting the poly(A) binding protein PABPC1 and the transcription factor FOXC2. Downregulation of either PABPC1 or FOXC2 in human endothelial cells (ECs) can induce the p16INK4a-dependent cellular senescence. Thus, we confirmed the utility of Reel-seq and FREP-MS analyses for the systematic identification of cis-REs at high resolution over a large region of the human genome.

KDM4B promotes acute myeloid leukemia associated with AML1-ETO by regulating chromatin accessibility.

Epigenetic alterations of chromatin structure affect chromatin accessibility and collaborate with genetic alterations in the development of cancer. Lysine demethylase 4B (KDM4B) has been identified as a JmjC domain‐containing epigenetic modifier that possesses histone demethylase activity. Although recent studies have demonstrated that KDM4B positively regulates the pathogenesis of multiple types of solid tumors, the tissue specificity and context dependency have not been fully elucidated. In this study, we investigated gene expression profiles established from clinical samples and found that KDM4B is elevated specifically in acute myeloid leukemia (AML) associated with chromosomal translocation 8;21 [t(8;21)], which results in a fusion of the AML1 and the eight‐twenty‐one (ETO) genes to generate a leukemia oncogene, AML1‐ETO fusion transcription factor. Short hairpin RNA‐mediated KDM4B silencing significantly reduced cell proliferation in t(8;21)‐positive AML cell lines. Meanwhile, KDM4B silencing suppressed the expression of AML1‐ETO‐inducible genes, and consistently perturbed chromatin accessibility of AML1‐ETO‐binding sites involving altered active enhancer marks and functional cis‐regulatory elements. Notably, transduction of murine KDM4B orthologue mutants followed by KDM4B silencing demonstrated a requirement of methylated‐histone binding modules for a proliferative surge. To address the role of KDM4B in leukemia development, we further generated and analyzed Kdm4b conditional knockout mice. As a result, Kdm4b deficiency attenuated clonogenic potential mediated by AML1‐ETO and delayed leukemia progression in vivo. Thus, our results highlight a tumor‐promoting role of KDM4B in AML associated with t(8;21).

The cis-Regulatory Element of SNCA Intron 4 Modulates Susceptibility to Parkinson's Disease in Han Chinese.

Objective: A novel functional cis-regulatory element (CRE) located at SNCA intron 4 has recently been identified in association with Parkinson’s disease (PD) risk in European descendants. We aimed to investigate whether this CRE is associated with PD in Han Chinese ethnicity.Methods: A Chinese cohort comprising 513 sporadic PD patients and 517 controls was recruited. CRE variants were identified by sequencing and then analyzed.Results: A total of nine variants were detected, namely eight single nucleotide variants and one new insertion variant. Two variants, rs17016188 and rs7684892, had minor allele frequency greater than 5%. A difference of rs17016188 was observed in males with the C allele serving as a recessive risk factor (p = 0.001, OR = 2.349, 95% CI = 1.414–3.901) following Bonferroni correction. Haplotypes of rs17016188 and rs7684892 showed distribution differences in the total and the male populations (p = 0.002 and 4.08 × 10−5, respectively). Among the haplotypes, rs17016188/T-rs7684892/G was associated with a reduced risk for PD (p = 4.8 × 10−4, OR = 0.731, 95% CI = 0.614–0.872).Conclusions: Our results provide insight into how the SNCA intron 4 CRE harbors variants and its contribution to PD risk in Chinese ethnicity.

Co-optation of Tandem DNA Repeats for the Maintenance of Mesenchymal Identity

Tandem repeats (TRs) are generated by DNA replication errors and retain a high level of instability, which in principle would make them unsuitable for integration into gene regulatory networks. However, the appearance of DNA sequence motifs recognized by transcription factors may turn TRs into functional cis-regulatory elements, thus favoring their stabilization in genomes. Here, we show that, in human cells, the transcriptional repressor ZEB1, which promotes the maintenance of mesenchymal features largely by suppressing epithelial genes and microRNAs, occupies TRs harboring dozens of copies of its DNA-binding motif within genomic loci relevant for maintenance of epithelial identity. The deletion of one such TR caused quasi-mesenchymal cancer cells to reacquire epithelial features, partially recapitulating the effects of ZEB1 gene deletion. These data demonstrate that the high density of identical motifs in TRs can make them suitable platforms for recruitment of transcriptional repressors, thus promoting their exaptation into pre-existing cis-regulatory networks.

Genome-wide identification of novel intergenic enhancer-like elements: implications in the regulation of transcription in Plasmodium falciparum

BackgroundThe molecular mechanisms of transcriptional regulation are poorly understood in Plasmodium falciparum. In addition, most of the genes in Plasmodium falciparum are transcriptionally poised and only a handful of cis-regulatory elements are known to operate in transcriptional regulation. Here, we employed an epigenetic signature based approach to identify significance of previously uncharacterised intergenic regions enriched with histone modification marks leading to discovery of enhancer-like elements.ResultsWe found that enhancer-like elements are significantly enriched with H3K4me1, generate unique non-coding bi-directional RNAs and majority of them can function as cis-regulators. Furthermore, functional enhancer reporter assay demonstrates that the enhancer-like elements regulate transcription of target genes in Plasmodium falciparum. Our study also suggests that the Plasmodium genome segregates functionally related genes into discrete housekeeping and pathogenicity/virulence clusters, presumably for robust transcriptional control of virulence/pathogenicity genes.ConclusionsThis report contributes to the understanding of parasite regulatory genomics by identification of enhancer-like elements, defining their epigenetic and transcriptional features and provides a resource of functional cis-regulatory elements that may give insights into the virulence/pathogenicity of Plasmodium falciparum.

Emergence of the Noncoding Cancer Genome: A Target of Genetic and Epigenetic Alterations.

The majority of genetic and epigenetic alterations accumulate in the noncoding genome throughout oncogenesis. Discriminating driver from passenger events is a challenge that holds great promise to improve our understanding of the etiology of different cancer types. Advancing our understanding of the noncoding cancer genome may thus identify new therapeutic opportunities and accelerate our capacity to find improved biomarkers to monitor various stages of cancer development. Cancer Discov; 6(11); 1215-29. ©2016 AACR.

In Silico Expression Analysis.

Information on the specificity of cis-sequences enables the design of functional synthetic plant promoters that are responsive to specific stresses. Potential cis-sequences may be experimentally tested, however, correlation of genomic sequence with gene expression data enables an in silico expression analysis approach to bioinformatically assess the stress specificity of candidate cis-sequences prior to experimental verification. The present chapter demonstrates an example for the in silico validation of a potential cis-regulatory sequence responsive to cold stress. The described online tool can be applied for the bioinformatic assessment of cis-sequences responsive to most abiotic and biotic stresses of plants. Furthermore, a method is presented based on a reverted in silico expression analysis approach that predicts highly specific potentially functional cis-regulatory elements for a given stress.

Muscle fiber type specific activation of the slow myosin heavy chain 2 promoter by a non-canonical E-box

Different mechanisms control skeletal muscle fiber type gene expression at specific times in vertebrate development. Embryonic myogenesis leading to formation of primary muscle fibers in avian species is largely directed by myoblast cell commitment to the formation of diverse fiber types. In contrast, development of different secondary fiber types during fetal myogenesis is partly determined by neural influences. In both primary and secondary chicken muscle fibers, differential expression of the slow myosin heavy chain 2 (MyHC2) gene distinguishes fast from fast/slow muscle fibers. This study focused on the transcriptional regulation of the slow MyHC2 gene in primary myotubes formed from distinct fast/slow and fast myogenic cell lineages. Promoter deletion analyses identified a discrete 86 bp promoter segment that conferred fiber type, lineage-specific gene expression in fast/slow versus fast myoblast derived primary myotubes. Sequence analysis and promoter activity assays determined that this segment contains two functional cis-regulatory elements. One element is a non-canonical E-box, and electromobility shift assays demonstrated that both cis-elements interacted with the E-protein, E47. The results indicate that primary muscle fiber type specific expression of the slow MyHC2 gene is controlled by a novel mechanism involving a transcriptional complex that includes E47 at a non-canonical E-box.

Functional Assessment of Disease-Associated Regulatory Variants In Vivo Using a Versatile Dual Colour Transgenesis Strategy in Zebrafish

Disruption of gene regulation by sequence variation in non-coding regions of the genome is now recognised as a significant cause of human disease and disease susceptibility. Sequence variants in cis-regulatory elements (CREs), the primary determinants of spatio-temporal gene regulation, can alter transcription factor binding sites. While technological advances have led to easy identification of disease-associated CRE variants, robust methods for discerning functional CRE variants from background variation are lacking. Here we describe an efficient dual-colour reporter transgenesis approach in zebrafish, simultaneously allowing detailed in vivo comparison of spatio-temporal differences in regulatory activity between putative CRE variants and assessment of altered transcription factor binding potential of the variant. We validate the method on known disease-associated elements regulating SHH, PAX6 and IRF6 and subsequently characterise novel, ultra-long-range SOX9 enhancers implicated in the craniofacial abnormality Pierre Robin Sequence. The method provides a highly cost-effective, fast and robust approach for simultaneously unravelling in a single assay whether, where and when in embryonic development a disease-associated CRE-variant is affecting its regulatory function.

A novel pairwise comparison method for in silico discovery of statistically significant cis-regulatory elements in eukaryotic promoter regions: application to Arabidopsis.

Cis regulatory elements (CREs), located within promoter regions, play a significant role in the blueprint for transcriptional regulation of genes. There is a growing interest to study the combinatorial nature of CREs including presence or absence of CREs, the number of occurrences of each CRE, as well as of their order and location relative to their target genes. Comparative promoter analysis has been shown to be a reliable strategy to test the significance of each component of promoter architecture. However, it remains unclear what level of difference in the number of occurrences of each CRE is of statistical significance in order to explain different expression patterns of two genes. In this study, we present a novel statistical approach for pairwise comparison of promoters of Arabidopsis genes in the context of number of occurrences of each CRE within the promoters. First, using the sample of 1000 Arabidopsis promoters, the results of the goodness of fit test and non-parametric analysis revealed that the number of occurrences of CREs in a promoter sequence is Poisson distributed. As a promoter sequence contained functional and non-functional CREs, we addressed the issue of the statistical distribution of functional CREs by analyzing the ChIP-seq datasets. The results showed that the number of occurrences of functional CREs over the genomic regions was determined as being Poisson distributed. In accordance with the obtained distribution of CREs occurrences, we suggested the Audic and Claverie (AC) test to compare two promoters based on the number of occurrences for the CREs. Superiority of the AC test over Chi-square (2×2) and Fisher's exact tests was also shown, as the AC test was able to detect a higher number of significant CREs. The two case studies on the Arabidopsis genes were performed in order to biologically verify the pairwise test for promoter comparison. Consequently, a number of CREs with significantly different occurrences was identified between the promoters. The results of the pairwise comparative analysis together with the expression data for the studied genes revealed the biological significance of the identified CREs.

TAF4, a subunit of transcription factor II D, directs promoter occupancy of nuclear receptor HNF4A during post-natal hepatocyte differentiation.

The functions of the TAF subunits of mammalian TFIID in physiological processes remain poorly characterised. In this study, we describe a novel function of TAFs in directing genomic occupancy of a transcriptional activator. Using liver-specific inactivation in mice, we show that the TAF4 subunit of TFIID is required for post-natal hepatocyte maturation. TAF4 promotes pre-initiation complex (PIC) formation at post-natal expressed liver function genes and down-regulates a subset of embryonic expressed genes by increased RNA polymerase II pausing. The TAF4-TAF12 heterodimer interacts directly with HNF4A and in vivo TAF4 is necessary to maintain HNF4A-directed embryonic gene expression at post-natal stages and promotes HNF4A occupancy of functional cis-regulatory elements adjacent to the transcription start sites of post-natal expressed genes. Stable HNF4A occupancy of these regulatory elements requires TAF4-dependent PIC formation highlighting that these are mutually dependent events. Local promoter-proximal HNF4A-TFIID interactions therefore act as instructive signals for post-natal hepatocyte differentiation.

Genome-wide identification and characterization of functional neuronal activity-dependent enhancers.

Experience-dependent gene transcription is required for nervous system development and function. However, the DNA regulatory elements that control this program of gene expression are not well defined. Here we characterize the enhancers that function across the genome to mediate activity-dependent transcription in mouse cortical neurons. We find that the subset of enhancers enriched for monomethylation of histone H3 Lys4 (H3K4me1) and binding of the transcriptional coactivator CREBBP (also called CBP) that shows increased acetylation of histone H3 Lys27 (H3K27ac) after membrane depolarization of cortical neurons functions to regulate activity-dependent transcription. A subset of these enhancers appears to require binding of FOS, which was previously thought to bind primarily to promoters. These findings suggest that FOS functions at enhancers to control activity-dependent gene programs that are critical for nervous system function and provide a resource of functional cis-regulatory elements that may give insight into the genetic variants that contribute to brain development and disease.