Reduction of ZFX levels decreases histone H4 acetylation and increases Pol2 pausing at target promoters.

The ZFX transcriptional activator binds to CpG island promoters, with a major peak at ∼200-250 bp downstream from transcription start sites. Because ZFX binds within the transcribed region, we investigated whether it regulates transcriptional elongation. We used GRO-seq to show that loss or reduction of ZFX increased Pol2 pausing at ZFX-regulated promoters. To further investigate the mechanisms by which ZFX regulates transcription, we determined regions of the protein needed for transactivation and for recruitment to the chromatin. Interestingly, although ZFX has 13 grouped zinc fingers, deletion of the first 11 fingers produces a protein that can still bind to chromatin and activate transcription. We next used TurboID-MS to detect ZFX-interacting proteins, identifying ZNF593, as well as proteins that interact with the N-terminal transactivation domain (which included histone modifying proteins), and proteins that interact with ZFX when it is bound to the chromatin (which included TAFs and other histone modifying proteins). Our studies support a model in which ZFX enhances elongation at target promoters by recruiting H4 acetylation complexes and reducing pausing.

m6A epitranscriptome analysis reveals differentially methylated transcripts that drive early chemoresistance in bladder cancer.

N 6-Methyladenosine (m6A) RNA modifications dynamically regulate messenger RNA processing, differentiation and cell fate. Given these functions, we hypothesized that m6A modifications play a role in the transition to chemoresistance. To test this, we took an agnostic discovery approach anchored directly to chemoresistance rather than to any particular m6A effector protein. Specifically, we used methyl-RNA immunoprecipitation followed by sequencing (MeRIP-seq) in parallel with RNA sequencing to identify gene transcripts that were both differentially methylated and differentially expressed between cisplatin-sensitive and cisplatin-resistant bladder cancer (BC) cells. We filtered and prioritized these genes using clinical and functional database tools, and then validated several of the top candidates via targeted quantitative polymerase chain reaction (qPCR) and MeRIP-PCR. In cisplatin-resistant cells, SLC7A11 transcripts had decreased methylation associated with decreased m6A reader YTHDF3 binding, prolonged RNA stability, and increased RNA and protein levels, leading to reduced ferroptosis and increased survival. Consistent with this, cisplatin-sensitive BC cell lines and patient-derived organoids exposed to cisplatin for as little as 48 h exhibited similar mechanisms of SLC7A11 upregulation and chemoresistance, trends that were also reflected in public cancer survival databases. Collectively, these findings highlight epitranscriptomic plasticity as a mechanism of rapid chemoresistance and a potential therapeutic target.

Large-scale manipulation of promoter DNA methylation reveals context-specific transcriptional responses and stability.

BACKGROUND: Cytosine DNA methylation is widely described as a transcriptional repressive mark with the capacity to silence promoters. Epigenome engineering techniques enable direct testing of the effect of induced DNA methylation on endogenous promoters; however, the downstream effects have not yet been comprehensively assessed. RESULTS: Here, we simultaneously induce methylation at thousands of promoters in human cells using an engineered zinc finger-DNMT3A fusion protein, enabling us to test the effect of forced DNA methylation upon transcription, chromatin accessibility, histone modifications, and DNA methylation persistence after the removal of the fusion protein. We find that transcriptional responses to DNA methylation are highly context-specific, including lack of repression, as well as cases of increased gene expression, which appears to be driven by the eviction of methyl-sensitive transcriptional repressors. Furthermore, we find that some regulatory networks can override DNA methylation and that promoter methylation can cause alternative promoter usage. DNA methylation deposited at promoter and distal regulatory regions is rapidly erased after removal of the zinc finger-DNMT3A fusion protein, in a process combining passive and TET-mediated demethylation. Finally, we demonstrate that induced DNA methylation can exist simultaneously on promoter nucleosomes that possess the active histone modification H3K4me3, or DNA bound by the initiated form of RNA polymerase II. CONCLUSIONS: These findings have important implications for epigenome engineering and demonstrate that the response of promoters to DNA methylation is more complex than previously appreciated.

FOXC1 Binds Enhancers and Promotes Cisplatin Resistance in Bladder Cancer.

Chemotherapy resistance is traditionally attributed to DNA mutations that confer a survival advantage under drug selection pressure. However, in bladder cancer and other malignancies, we and others have previously reported that cancer cells can convert spontaneously to an aggressive drug-resistant phenotype without prior drug selection or mutational events. In the current work, we explored possible epigenetic mechanisms behind this phenotypic plasticity. Using Hoechst dye exclusion and flow cytometry, we isolated the aggressive drug-resistant cells and analyzed their chromatin accessibility at regulatory elements. Compared to the rest of the cancer cell population, the aggressive drug-resistant cells exhibited enhancer accessibility changes. In particular, we found that differentially accessible enhancers were enriched for the FOXC1 transcription factor motif, and that FOXC1 was the most significantly overexpressed gene in aggressive drug-resistant cells. ChIP-seq analysis revealed that differentially accessible enhancers in aggressive drug-resistant cells had a higher FOXC1 binding, which regulated the expression of adjacent cancer-relevant genes like ABCB1 and ID3. Accordingly, cisplatin treatment of bladder cancer cells led to an increased FOXC1 expression, which mediated cell survival and conversion to a drug-resistant phenotype. Collectively, these findings suggest that FOXC1 contributes to phenotypic plasticity by binding enhancers and promoting a mutation-independent shift towards cisplatin resistance in bladder cancer.

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.

Expanded encyclopaedias of DNA elements in the human and mouse genomes.

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

Characterization of the ZFX family of transcription factors that bind downstream of the start site of CpG island promoters.

Our study focuses on a family of ubiquitously expressed human C2H2 zinc finger proteins comprised of ZFX, ZFY and ZNF711. Although their protein structure suggests that ZFX, ZFY and ZNF711 are transcriptional regulators, the mechanisms by which they influence transcription have not yet been elucidated. We used CRISPR-mediated deletion to create bi-allelic knockouts of ZFX and/or ZNF711 in female HEK293T cells (which naturally lack ZFY). We found that loss of either ZFX or ZNF711 reduced cell growth and that the double knockout cells have major defects in proliferation. RNA-seq analysis revealed that thousands of genes showed altered expression in the double knockout clones, suggesting that these TFs are critical regulators of the transcriptome. To gain insight into how these TFs regulate transcription, we created mutant ZFX proteins and analyzed them for DNA binding and transactivation capability. We found that zinc fingers 11-13 are necessary and sufficient for DNA binding and, in combination with the N terminal region, constitute a functional transactivator. Our functional analyses of the ZFX family provides important new insights into transcriptional regulation in human cells by members of the large, but under-studied family of C2H2 zinc finger proteins.

Genome-wide analysis of HOXC4 and HOXC6 regulated genes and binding sites in prostate cancer cells.

Aberrant expression of HOXC6 and HOXC4 is commonly detected in prostate cancer. The high expression of these transcription factors is associated with aggressive prostate cancer and can predict cancer recurrence after treatment. Thus, HOXC4 and HOXC6 are clinically relevant biomarkers of aggressive prostate cancer. However, the molecular mechanisms by which these HOXC genes contribute to prostate cancer is not yet understood. To begin to address the role of HOXC4 and HOXC6 in prostate cancer, we performed RNA-seq analyses before and after siRNA-mediated knockdown of HOXC4 and/or HOXC6 and also performed ChIP-seq to identify genomic binding sites for both of these transcription factors. Our studies demonstrate that HOXC4 and HOXC6 co-localize with HOXB13, FOXA1 and AR, three transcription factors previously shown to contribute to the development of prostate cancer. We suggest that the aberrantly upregulated HOXC4 and HOXC6 proteins may compete with HOXB13 for binding sites, thus altering the prostate transcriptome. This competition model may be applicable to many different human cancers that display increased expression of a HOX transcription factor.

A high-resolution 3D epigenomic map reveals insights into the creation of the prostate cancer transcriptome.

To better understand the impact of chromatin structure on regulation of the prostate cancer transcriptome, we develop high-resolution chromatin interaction maps in normal and prostate cancer cells using in situ Hi-C. By combining the in situ Hi-C data with active and repressive histone marks, CTCF binding sites, nucleosome-depleted regions, and transcriptome profiling, we identify topologically associating domains (TADs) that change in size and epigenetic states between normal and prostate cancer cells. Moreover, we identify normal and prostate cancer-specific enhancer-promoter loops and involved transcription factors. For example, we show that FOXA1 is enriched in prostate cancer-specific enhancer-promoter loop anchors. We also find that the chromatin structure surrounding the androgen receptor (AR) locus is altered in the prostate cancer cells with many cancer-specific enhancer-promoter loops. This creation of 3D epigenomic maps enables a better understanding of prostate cancer biology and mechanisms of gene regulation.

The prostate cancer risk variant rs55958994 regulates multiple gene expression through extreme long-range chromatin interaction to control tumor progression.

Genome-wide association studies identified single-nucleotide polymorphism (SNP) rs55958994 as a significant variant associated with increased susceptibility to prostate cancer. However, the mechanisms by which this SNP mediates increased risk to cancer are still unknown. In this study, we show that this variant is located in an enhancer active in prostate cancer cells. Deletion of this enhancer from prostate tumor cells resulted in decreased tumor initiation, tumor growth, and invasive migration, as well as a loss of stem-like cells. Using a combination of capture chromosome conformation capture (Capture-C) and RNA sequencing, we identified genes on the same and different chromosomes as targets regulated by the enhancer. Furthermore, we show that expression of individual candidate target genes in an enhancer-deleted cell line rescued different aspects of tumorigenesis. Our data suggest that the rs55958994-associated enhancer affects prostate cancer progression by influencing expression of multiple genes via long-range chromatin interactions.

Ezh2-dCas9 and KRAB-dCas9 enable engineering of epigenetic memory in a context-dependent manner.

BACKGROUND: Rewriting of the epigenome has risen as a promising alternative to gene editing for precision medicine. In nature, epigenetic silencing can result in complete attenuation of target gene expression over multiple mitotic divisions. However, persistent repression has been difficult to achieve in a predictable manner using targeted systems. RESULTS: Here, we report that persistent epigenetic memory required both a DNA methyltransferase (DNMT3A-dCas9) and a histone methyltransferase (Ezh2-dCas9 or KRAB-dCas9). We demonstrate that the histone methyltransferase requirement can be locus specific. Co-targeting Ezh2-dCas9, but not KRAB-dCas9, with DNMT3A-dCas9 and DNMT3L induced long-term HER2 repression over at least 50 days (approximately 57 cell divisions) and triggered an epigenetic switch to a heterochromatic environment. An increase in H3K27 trimethylation and DNA methylation was stably maintained and accompanied by a sustained loss of H3K27 acetylation. Interestingly, substitution of Ezh2-dCas9 with KRAB-dCas9 enabled long-term repression at some target genes (e.g., SNURF) but not at HER2, at which H3K9me3 and DNA methylation were transiently acquired and subsequently lost. Off-target DNA hypermethylation occurred at many individual CpG sites but rarely at multiple CpGs in a single promoter, consistent with no detectable effect on transcription at the off-target loci tested. Conversely, robust hypermethylation was observed at HER2. We further demonstrated that Ezh2-dCas9 required full-length DNMT3L for maximal activity and that co-targeting DNMT3L was sufficient for persistent repression by Ezh2-dCas9 or KRAB-dCas9. CONCLUSIONS: These data demonstrate that targeting different combinations of histone and DNA methyltransferases is required to achieve maximal repression at different loci. Fine-tuning of targeting tools is a necessity to engineer epigenetic memory at any given locus in any given cell type.

The Enigmatic HOX Genes: Can We Crack Their Code?

Homeobox genes (HOX) are a large family of transcription factors that direct the formation of many body structures during early embryonic development. There are 39 genes in the subgroup of homeobox genes that constitute the human HOX gene family. Correct embryonic development of flies and vertebrates is, in part, mediated by the unique and highly regulated expression pattern of the HOX genes. Disruptions in these fine-tuned regulatory mechanisms can lead to developmental problems and to human diseases such as cancer. Unfortunately, the molecular mechanisms of action of the HOX family of transcription factors are severely under-studied, likely due to idiosyncratic details of their structure, expression, and function. We suggest that a concerted and collaborative effort to identify interacting protein partners, produce genome-wide binding profiles, and develop HOX network inhibitors in a variety of human cell types will lead to a deeper understanding of human development and disease. Within, we review the technological challenges and possible approaches needed to achieve this goal.

Three-dimensional analysis reveals altered chromatin interaction by enhancer inhibitors harbors TCF7L2-regulated cancer gene signature.

Distal regulatory elements influence the activity of gene promoters through chromatin looping. Chromosome conformation capture (3C) methods permit identification of chromatin contacts across different regions of the genome. However, due to limitations in the resolution of these methods, the detection of functional chromatin interactions remains a challenge. In the current study, we employ an integrated approach to define and characterize the functional chromatin contacts of human pancreatic cancer cells. We applied tethered chromatin capture to define classes of chromatin domains on a genome-wide scale. We identified three types of structural domains (topologically associated, boundary, and gap) and investigated the functional relationships of these domains with respect to chromatin state and gene expression. We uncovered six distinct sub-domains associated with epigenetic states. Interestingly, specific epigenetically active domains are sensitive to treatment with histone acetyltransferase (HAT) inhibitors and decrease in H3K27 acetylation levels. To examine whether the subdomains that change upon drug treatment are functionally linked to transcription factor regulation, we compared TCF7L2 chromatin binding and gene regulation to HAT inhibition. We identified a subset of coding RNA genes that together can stratify pancreatic cancer patients into distinct survival groups. Overall, this study describes a process to evaluate the functional features of chromosome architecture and reveals the impact of epigenetic inhibitors on chromosome architecture and identifies genes that may provide insight into disease outcome.

Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation.

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium.

CRISPR-mediated deletion of prostate cancer risk-associated CTCF loop anchors identifies repressive chromatin loops.

BACKGROUND: Recent genome-wide association studies (GWAS) have identified more than 100 loci associated with increased risk of prostate cancer, most of which are in non-coding regions of the genome. Understanding the function of these non-coding risk loci is critical to elucidate the genetic susceptibility to prostate cancer. RESULTS: We generate genome-wide regulatory element maps and performed genome-wide chromosome confirmation capture assays (in situ Hi-C) in normal and tumorigenic prostate cells. Using this information, we annotate the regulatory potential of 2,181 fine-mapped prostate cancer risk-associated SNPs and predict a set of target genes that are regulated by prostate cancer risk-related H3K27Ac-mediated loops. We next identify prostate cancer risk-associated CTCF sites involved in long-range chromatin loops. We use CRISPR-mediated deletion to remove prostate cancer risk-associated CTCF anchor regions and the CTCF anchor regions looped to the prostate cancer risk-associated CTCF sites, and we observe up to 100-fold increases in expression of genes within the loops when the prostate cancer risk-associated CTCF anchor regions are deleted. CONCLUSIONS: We identify GWAS risk loci involved in long-range loops that function to repress gene expression within chromatin loops. Our studies provide new insights into the genetic susceptibility to prostate cancer.

Defining Regulatory Elements in the Human Genome Using Nucleosome Occupancy and Methylome Sequencing (NOMe-Seq).

NOMe-seq (nucleosome occupancy and methylome sequencing) identifies nucleosome-depleted regions that correspond to promoters, enhancers, and insulators. The NOMe-seq method is based on the treatment of chromatin with the M.CviPI methyltransferase, which methylates GpC dinucleotides that are not protected by nucleosomes or other proteins that are tightly bound to the chromatin (GpCm does not occur in the human genome and therefore there is no endogenous background of GpCm). Following bisulfite treatment of the M.CviPI-methylated chromatin (which converts unmethylated Cs to Ts and thus allows the distinction of GpC from GpCm) and subsequent genomic sequencing, nucleosome-depleted regions can be ascertained on a genome-wide scale. The bisulfite treatment also allows the distinction of CpG from CmpG (most endogenous methylation occurs at CpG dinucleotides) and thus the endogenous methylation status of the genome can also be obtained in the same sequencing reaction. Importantly, open chromatin is expected to have high levels of GpCm but low levels of CmpG; thus, each of the two separate methylation analyses serve as independent (but opposite) measures which provide matching chromatin designations for each regulatory element.NOMe-seq has advantages over ChIP-seq for identification of regulatory elements because it is not reliant upon knowing the exact modifications on the surrounding nucleosomes. Also, NOMe-seq has advantages over DHS (DNase hypersensitive site)-seq, FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements)-seq, and ATAC (Assay for Transposase-Accessible Chromatin)-seq because it also gives positioning information for several nucleosomes on either side of each open regulatory element. Here, we provide a detailed protocol for NOMe-seq that begins with the isolation of chromatin, followed by methylation of GpCs with M.CviPI and treatment with bisulfite, and ending with the creation of next generation sequencing libraries. We also include sequencing QC analysis metrics and bioinformatics steps that can be used to identify nucleosome-depleted regions throughout the genome.

ZFX acts as a transcriptional activator in multiple types of human tumors by binding downstream of transcription start sites at the majority of CpG island promoters.

High expression of the transcription factor ZFX is correlated with proliferation, tumorigenesis, and patient survival in multiple types of human cancers. However, the mechanism by which ZFX influences transcriptional regulation has not been determined. We performed ChIP-seq in four cancer cell lines (representing kidney, colon, prostate, and breast cancers) to identify ZFX binding sites throughout the human genome. We identified ~9,000 ZFX binding sites and found that the majority of the sites are in CpG island promoters. Moreover, genes with promoters bound by ZFX are expressed at higher levels than genes with promoters not bound by ZFX. To determine if ZFX contributes to regulation of the promoters to which it is bound, we performed RNA-seq analysis after knockdown of ZFX by siRNA in prostate and breast cancer cells. Many genes with promoters bound by ZFX were downregulated upon ZFX knockdown, supporting the hypothesis that ZFX acts as a transcriptional activator. Surprisingly, ZFX binds at +240 bp downstream of the TSS of the responsive promoters. Using Nucleosome Occupancy and Methylome Sequencing (NOMe-seq), we show that ZFX binds between the open chromatin region at the TSS and the first downstream nucleosome, suggesting that ZFX may play a critical role in promoter architecture. We have also shown that a closely related zinc finger protein ZNF711 has a similar binding pattern at CpG island promoters, but ZNF711 may play a subordinate role to ZFX. This functional characterization of ZFX provides important new insights into transcription, chromatin structure, and the regulation of the cancer transcriptome.

A Prostate Cancer Risk Element Functions as a Repressive Loop that Regulates HOXA13.

Prostate cancer (PCa) is the leading cancer among men in the United States, with genetic factors contributing to ∼42% of the susceptibility to PCa. We analyzed a PCa risk region located at 7p15.2 to gain insight into the mechanisms by which this noncoding region may affect gene regulation and contribute to PCa risk. We performed Hi-C analysis and demonstrated that this region has long-range interactions with the HOXA locus, located ∼873 kb away. Using the CRISPR/Cas9 system, we deleted a 4-kb region encompassing several PCa risk-associated SNPs and performed RNA-seq to investigate transcriptomic changes in prostate cells lacking the regulatory element. Our results suggest that the risk element affects the expression of HOXA13 and HOTTIP, but not other genes in the HOXA locus, via a repressive loop. Forced expression of HOXA13 was performed to gain further insight into the mechanisms by which this risk element affects PCa risk.

dCas9-based epigenome editing suggests acquisition of histone methylation is not sufficient for target gene repression.

Distinct epigenomic profiles of histone marks have been associated with gene expression, but questions regarding the causal relationship remain. Here we investigated the activity of a broad collection of genomically targeted epigenetic regulators that could write epigenetic marks associated with a repressed chromatin state (G9A, SUV39H1, Krüppel-associated box (KRAB), DNMT3A as well as the first targetable versions of Ezh2 and Friend of GATA-1 (FOG1)). dCas9 fusions produced target gene repression over a range of 0- to 10-fold that varied by locus and cell type. dCpf1 fusions were unable to repress gene expression. The most persistent gene repression required the action of several effector domains; however, KRAB-dCas9 did not contribute to persistence in contrast to previous reports. A 'direct tethering' strategy attaching the Ezh2 methyltransferase enzyme to dCas9, as well as a 'recruitment' strategy attaching the N-terminal 45 residues of FOG1 to dCas9 to recruit the endogenous nucleosome remodeling and deacetylase complex, were both successful in targeted deposition of H3K27me3. Surprisingly, however, repression was not correlated with deposition of either H3K9me3 or H3K27me3. Our results suggest that so-called repressive histone modifications are not sufficient for gene repression. The easily programmable dCas9 toolkit allowed precise control of epigenetic information and dissection of the relationship between the epigenome and gene regulation.