scholarly journals Integrative epigenomic and high-throughput functional enhancer profiling reveals determinants of enhancer heterogeneity in gastric cancer

2021 ◽  
Vol 13 (1) ◽  
Author(s):  
Taotao Sheng ◽  
Shamaine Wei Ting Ho ◽  
Wen Fong Ooi ◽  
Chang Xu ◽  
Manjie Xing ◽  
...  

Abstract Background Enhancers are distal cis-regulatory elements required for cell-specific gene expression and cell fate determination. In cancer, enhancer variation has been proposed as a major cause of inter-patient heterogeneity—however, most predicted enhancer regions remain to be functionally tested. Methods We analyzed 132 epigenomic histone modification profiles of 18 primary gastric cancer (GC) samples, 18 normal gastric tissues, and 28 GC cell lines using Nano-ChIP-seq technology. We applied Capture-based Self-Transcribing Active Regulatory Region sequencing (CapSTARR-seq) to assess functional enhancer activity. An Activity-by-contact (ABC) model was employed to explore the effects of histone acetylation and CapSTARR-seq levels on enhancer-promoter interactions. Results We report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are GC-associated in vivo (> 50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying CapSTARR-seq to the enhancer catalog, we observed significant correlations between CapSTARR-seq functional activity and H3K27ac/H3K4me1 levels. Super-enhancer regions exhibited increased CapSTARR-seq signals compared to regular enhancers, even when decoupled from native chromatin contexture. We show that combining histone modification and CapSTARR-seq functional enhancer data improves the prediction of enhancer-promoter interactions and pinpointing of germline single nucleotide polymorphisms (SNPs), somatic copy number alterations (SCNAs), and trans-acting TFs involved in GC expression. We identified cancer-relevant genes (ING1, ARL4C) whose expression between patients is influenced by enhancer differences in genomic copy number and germline SNPs, and HNF4α as a master trans-acting factor associated with GC enhancer heterogeneity. Conclusions Our results indicate that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, providing insights into the relative contribution of genetic (cis) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.


2021 ◽  
Author(s):  
Taotao Sheng ◽  
Shamaine Wei Ting Ho ◽  
Wen Fong Ooi ◽  
Chang Xu ◽  
Manjie Xing ◽  
...  

Background: Enhancers are distal cis-regulatory elements required for cell-specific gene expression and cell fate determination. In cancer, enhancer variation has been proposed as a major cause of inter-patient heterogeneity – however, most predicted enhancer regions remain to be functionally tested. Results: Analyzing 128 epigenomic histone modification profiles of primary GC samples, normal gastric tissues, and GC cell lines, we report a comprehensive catalog of 75,730 recurrent predicted enhancers, the majority of which are tumor-associated in vivo (>50,000) and associated with lower somatic mutation rates inferred by whole-genome sequencing. Applying Capture-based Self-Transcribing Active Regulatory Region sequencing (CapSTARR-seq) to the enhancer catalog, we observed significant correlations between CapSTARR-seq functional activity and H3K27ac/H3K4me1 levels. Super-enhancer regions exhibited increased CapSTARR-seq signals compared to regular enhancers even when decoupled from native chromatin contexture. We show that combining histone modification and CapSTARR-seq functional enhancer data improves the prediction of enhancer-promoter interactions and pinpointing of germline single nucleotide polymorphisms (SNPs), somatic copy number alterations (SCNAs), and trans-acting TFs involved in GC expression. Specifically, we identified cancer-relevant genes (e.g. ING1, ARL4C) whose expression between patients is influenced by enhancer differences in genomic copy number and germline SNPs, and HNF4α as a master trans-acting factor associated with GC enhancer heterogeneity. Conclusions: Our study indicates that combining histone modification and functional assay data may provide a more accurate metric to assess enhancer activity than either platform individually, and provides insights into the relative contribution of genetic (cis ) and regulatory (trans) mechanisms to GC enhancer functional heterogeneity.



2019 ◽  
Author(s):  
Kailong Li ◽  
Yuxuan Liu ◽  
Hui Cao ◽  
Yuannyu Zhang ◽  
Zhimin Gu ◽  
...  

ABSTRACTTissue-specific gene expression requires coordinated control of gene-proximal and -distal cis-regulatory elements (CREs), yet functional analysis of gene-distal CREs such as enhancers remains challenging. Here we describe enhanced CRISPR/dCas9-based epigenetic editing systems, enCRISPRa and enCRISPRi, for multiplexed analysis of enhancer function in situ and in vivo. Using dual effectors capable of re-writing enhancer-associated chromatin modifications, we show that enCRISPRa and enCRISPRi modulate gene transcription by remodeling local epigenetic landscapes at sgRNA-targeted enhancers and associated genes. Comparing with existing methods, the new systems display more robust perturbation of enhancer activity and gene transcription with minimal off-targets. Allele-specific targeting of enCRISPRa to oncogenic TAL1 super-enhancer modulates TAL1 expression and cancer progression in xenotransplants. Multiplexed perturbations of lineage-specific enhancers using an enCRISPRi knock-in mouse establish in vivo evidence for lineage-restricted essentiality of developmental enhancers during hematopoietic lineage specification. Hence, enhanced CRSIPR epigenetic editing provides opportunities for interrogating enhancer function in native biological contexts.



2002 ◽  
Vol 22 (12) ◽  
pp. 4293-4308 ◽  
Author(s):  
Charles K. Kaufman ◽  
Satrajit Sinha ◽  
Diana Bolotin ◽  
Jie Fan ◽  
Elaine Fuchs

ABSTRACT In this report, we explored the mechanisms underlying keratinocyte-specific and differentiation-specific gene expression in the skin. We have identified five keratinocyte-specific, open chromatin regions that exist within the 6 kb of 5′ upstream regulatory sequence known to faithfully recapitulate the strong endogenous keratin 5 (K5) promoter and/or enhancer activity. One of these, DNase I-hypersensitive site (HSs) 4, was unique in that it acted independently to drive abundant and keratinocyte-specific reporter gene activity in culture and in transgenic mice, despite the fact that it was not essential for K5 enhancer activity. We have identified evolutionarily conserved regulatory elements and a number of their associated proteins that bind to this compact and complex enhancer element. The 125-bp 3′ half of this element (referred to as 4.2) is by far the smallest known strong enhancer element possessing keratinocyte-specific activity in vivo. Interestingly, its activity is restricted to a subset of progeny of K5-expressing cells located within the sebaceous gland. The other half of HSs 4 (termed 4.1) possesses activity to suppress sebocyte-specific expression and induce expression in the channel (inner root sheath) cells surrounding the hair shaft. Our findings lead us to a view of keratinocyte gene expression which is determined by multiple regulatory modules, many of which contain AP-2 and/or Sp1/Sp3 binding sites for enhancing expression in skin epithelium, but which also harbor one or more unique sites for the binding of factors which determine specificity. Through mixing and matching of these modules, additional levels of specificity are obtained, indicating that both transcriptional repressors and activators govern the specificity.



2021 ◽  
Vol 49 (7) ◽  
pp. 3856-3875
Author(s):  
Marina Kulik ◽  
Melissa Bothe ◽  
Gözde Kibar ◽  
Alisa Fuchs ◽  
Stefanie Schöne ◽  
...  

Abstract The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied them in an equivalent cellular context. Analysis of chromatin and sequence suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the result of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.



Author(s):  
Fernanda M Bosada ◽  
Mathilde R Rivaud ◽  
Jae-Sun Uhm ◽  
Sander Verheule ◽  
Karel van Duijvenboden ◽  
...  

Rationale: Atrial Fibrillation (AF) is the most common cardiac arrhythmia diagnosed in clinical practice. Genome-wide association studies have identified AF-associated common variants across 100+ genomic loci, but the mechanism underlying the impact of these variant loci on AF susceptibility in vivo has remained largely undefined. One such variant region, highly associated with AF, is found at 1q24, close to PRRX1, encoding the Paired Related Homeobox 1 transcription factor. Objective: To identify the mechanistic link between the variant region at 1q24 and AF predisposition. Methods and Results: The mouse orthologue of the noncoding variant genomic region (R1A) at 1q24 was deleted using CRISPR genome editing. Among the genes sharing the topologically associated domain with the deleted R1A region (Kifap3, Prrx1, Fmo2, Prrc2c), only the broadly expressed gene Prrx1 was downregulated in mutants, and only in cardiomyocytes. Expression and epigenetic profiling revealed that a cardiomyocyte lineage-specific gene program (Mhrt, Myh6, Rbm20, Tnnt2, Ttn, Ckm) was upregulated in R1A-/- atrial cardiomyocytes, and that Mef2 binding motifs were significantly enriched at differentially accessible chromatin sites. Consistently, Prrx1 suppressed Mef2-activated enhancer activity in HL-1 cells. Mice heterozygous or homozygous for the R1A deletion were susceptible to atrial arrhythmia induction, had atrial conduction slowing and more irregular RR intervals. Isolated R1A-/- mouse left atrial cardiomyocytes showed lower action potential upstroke velocities and sodium current, as well as increased systolic and diastolic calcium concentrations compared to controls. Conclusions: The noncoding AF variant region at 1q24 modulates Prrx1 expression in cardiomyocytes. Cardiomyocyte-specific reduction of Prrx1 expression upon deletion of the noncoding region leads to a profound induction of a cardiac lineage-specific gene program and to propensity for AF. These data indicate that AF-associated variants in humans may exert AF predisposition through reduced PRRX1 expression in cardiomyocytes.



Development ◽  
1997 ◽  
Vol 124 (24) ◽  
pp. 4971-4982 ◽  
Author(s):  
Z. Yin ◽  
X.L. Xu ◽  
M. Frasch

The Drosophila tinman homeobox gene has a major role in early mesoderm patterning and determines the formation of visceral mesoderm, heart progenitors, specific somatic muscle precursors and glia-like mesodermal cells. These functions of tinman are reflected in its dynamic pattern of expression, which is characterized by initial widespread expression in the trunk mesoderm, then refinement to a broad dorsal mesodermal domain, and finally restricted expression in heart progenitors. Here we show that each of these phases of expression is driven by a discrete enhancer element, the first being active in the early mesoderm, the second in the dorsal mesoderm and the third in cardioblasts. We provide evidence that the early-active enhancer element is a direct target of twist, a gene encoding a basic helix-loop-helix (bHLH) protein, which is necessary for tinman activation. This 180 bp enhancer includes three E-box sequences which bind Twist protein in vitro and are essential for enhancer activity in vivo. Ectodermal misexpression of twist causes ectopic activation of this enhancer in ectodermal cells, indicating that twist is the only mesoderm-specific activator of early tinman expression. We further show that the 180 bp enhancer also includes negatively acting sequences. Binding of Even-skipped to these sequences appears to reduce twist-dependent activation in a periodic fashion, thus producing a striped tinman pattern in the early mesoderm. In addition, these sequences prevent activation of tinman by twist in a defined portion of the head mesoderm that gives rise to hemocytes. We find that this repression requires the function of buttonhead, a head-patterning gene, and that buttonhead is necessary for normal activation of the hematopoietic differentiation gene serpent in the same area. Together, our results show that tinman is controlled by an array of discrete enhancer elements that are activated successively by differential genetic inputs, as well as by closely linked activator and repressor binding sites within an early-acting enhancer, which restrict twist activity to specific areas within the twist expression domain.



1996 ◽  
Vol 16 (6) ◽  
pp. 2777-2786 ◽  
Author(s):  
V Gailus-Durner ◽  
J Xie ◽  
C Chintamaneni ◽  
A K Vershon

The meiosis-specific gene HOP1, which encodes a component of the synaptonemal complex, is controlled through two regulatory elements, UASH and URS1H. Sites similar to URS1H have been identified in the promoter region of virtually every early meiosis-specific gene, as well as in many promoters of nonmeiotic genes, and it has been shown that the proteins that bind to this site function to regulate meiotic and nonmeiotic transcription. Sites similar to the UASH site have been found in a number of meiotic and nonmeiotic genes as well. Since it has been shown that UASH functions as an activator site in vegetative haploid cells, it seemed likely that the factors binding to this site regulate both meiotic and nonmeiotic transcription. We purified the factor binding to the UASH element of the HOP1 promoter. Sequence analysis identified the protein as Abf1 (autonomously replicating sequence-binding factor 1), a multifunctional protein involved in DNA replication, silencing, and transcriptional regulation. We show by mutational analysis of the UASH site, that positions outside of the proposed UASH consensus sequence (TNTGN[A/T]GT) are required for DNA binding in vitro and transcriptional activation in vivo. A new UASH consensus sequence derived from this mutational analysis closely matches a consensus Abf1 binding site. We also show that an Abf1 site from a nonmeiotic gene can replace the function of the UASH site in the HOP1 promoter. Taken together, these results show that Abf1 functions to regulate meiotic gene expression.



Blood ◽  
2001 ◽  
Vol 98 (3) ◽  
pp. 610-617 ◽  
Author(s):  
Chunyan Zhang ◽  
Michael A. Thornton ◽  
M. Anna Kowalska ◽  
Bruce S. Sachis ◽  
Michael Feldman ◽  
...  

Abstract The genes for the related human (h) chemokines, PBP (platelet basic protein) and PF4 (platelet factor 4), are within 5.3 kilobases (kb) of each other and form a megakaryocyte-specific gene locus. The hypothesis was considered that the PBP and PF4 genes share a common distal regulatory region(s) that leads to their high-level megakaryocyte-specific expression in vivo. This study examined PBP and PF4 expression in transgenic mice using 4 distinct humanPBP/PF4 gene locus constructs. These studies showed that within the region studied there was sufficient information to regulate tissue-specific expression of both hPBP and hPF4. Indeed this region contained sufficient DNA information to lead to expression levels of PBP and PF4 comparable to the homologous mouse genes in a position-independent, copy number–dependent fashion. These studies also indicated that the DNA domains that led to this expression were distinct for the 2 genes; hPBP expression is regulated by a region that is 1.5 to 4.4 kb upstream of that gene. Expression of hPF4 is regulated by a region that is either intergenic between the 2 genes or immediately downstream of the hPF4 gene. Comparison of the available human and mouse sequences shows conserved flanking region domains containing potential megakaryocyte-related transcriptional factor DNA-binding sites. Further analysis of these regulatory regions may identify enhancer domains involved in megakaryopoiesis that may be useful in the selective expression of other genes in megakaryocytes and platelets as a strategy for regulating hemostasis, thrombosis, and inflammation.



1995 ◽  
Vol 15 (2) ◽  
pp. 642-652 ◽  
Author(s):  
A Pierani ◽  
C Pouponnot ◽  
G Calothy

Developmental control of gene expression often results from the coupling of growth arrest with the establishment of differentiation programs. QR1 is a gene specifically expressed in retinas during the late phase of embryogenesis. At this stage neuroectodermal precursors have reached terminal mitosis and are undergoing differentiation into distinct cell types. Transcription of the QR1 gene is tightly regulated during retinal development: this gene is expressed between embryonic day 9 (ED9) and ED17 and is completely repressed at hatching in quail. Moreover, QR1 transcription is downregulated when postmitotic neural retina cells are induced to proliferate by pp60v-src. We studied the stage-dependent transcriptional control of this gene during quail neural retina (QNR) cell development. Transient transfection experiments with QR1/CAT constructs at various stages of development showed that a region located between -935 and -1265 bp upstream of the transcription start site is necessary to promote transcription in retina cells during the late phase of embryonal development (QNR9, corresponding to ED9). By in vivo footprinting assays we identified at least two elements that are occupied by DNA-protein complexes in QNR cells: the A and B boxes. The A box allows formation of several biochemically distinct complexes: C1, C2, C3, and C4. Formation of the C2 complex mainly during early stages (ED7) and of C2, C3, and C4 complexes during postnatal life correlates with repression of QR1 transcription, whereas the C1 complex is strongly induced at ED11 when the QR1 gene is expressed. We previously showed that C1 was involved in downregulation of QR1 transcription by pp60v-src. Several complexes are also formed on the B box. We show that these complexes are exclusively present in neural tissues and that they involve members of the POU family of transcription factors. Mutations of each one of the two regions which abolish the binding of the C1 factor(s) on the A box and of the POU factor(s) on the B box also prevent stimulation of QR1 transcription in QNR9. Therefore, both elements appear to be required for the stage-specific transcription of the QR1 gene. We also show that the regulatory region from position -1265 to position -935 is able to confer stage-specific transcription upon a heterologous promoter (thymidine kinase). Indeed, this region stimulates transcription in differentiating retinas (QNR9) and represses transcription in terminally differentiated retinas (QNR17, corresponding to postnatal life). Our results suggest that cell growth regulation and developmental control are coordinated through the A and B boxes in regulating QR1 transcription during retinal differentiation.



2020 ◽  
Vol 29 (10) ◽  
pp. 1689-1699
Author(s):  
Harrison Pantera ◽  
Bo Hu ◽  
Daniel Moiseev ◽  
Chris Dunham ◽  
Jibraan Rashid ◽  
...  

Abstract Copy number variation of the peripheral nerve myelin gene Peripheral Myelin Protein 22 (PMP22) causes multiple forms of inherited peripheral neuropathy. The duplication of a 1.4 Mb segment surrounding this gene in chromosome 17p12 (c17p12) causes the most common form of Charcot-Marie-Tooth disease type 1A, whereas the reciprocal deletion of this gene causes a separate neuropathy termed hereditary neuropathy with liability to pressure palsies (HNPP). PMP22 is robustly induced in Schwann cells in early postnatal development, and several transcription factors and their cognate regulatory elements have been implicated in coordinating the gene’s proper expression. We previously found that a distal super-enhancer domain was important for Pmp22 expression in vitro, with particular impact on a Schwann cell-specific alternative promoter. Here, we investigate the consequences of deleting this super-enhancer in vivo. We find that loss of the super-enhancer in mice reduces Pmp22 expression throughout development and into adulthood, with greater impact on the Schwann cell-specific promoter. Additionally, these mice display tomacula formed by excessive myelin folding, a pathological hallmark of HNPP, as have been previously observed in heterozygous Pmp22 mice as well as sural biopsies from patients with HNPP. Our findings demonstrate a mechanism by which smaller copy number variations, not including the Pmp22 gene, are sufficient to reduce gene expression and phenocopy a peripheral neuropathy caused by the HNPP-associated deletion encompassing PMP22.



Sign in / Sign up

Export Citation Format

Share Document