Effects on the transcriptome upon deletion of a distal element cannot be predicted by the size of the H3K27Ac peak in human cells

Application of CRISPR-Cas9 gene editing technology in basic research, diagnosis and treatment of colon cancer

Colon cancer is the fourth leading cause of cancer death worldwide, and its progression is accompanied by a complex array of genetic variations. CRISPR/Cas9 can identify new drug-resistant or sensitive mutations in colon cancer, and can use gene editing technology to develop new therapeutic targets and provide personalized treatments, thereby significantly improving the treatment of colon cancer patients. CRISPR/Cas9 systems are driving advances in biotechnology. RNA-directed Cas enzymes have accelerated the pace of basic research and led to clinical breakthroughs. This article reviews the rapid development of CRISPR/Cas in colon cancer, from gene editing to transcription regulation, gene knockout, genome-wide CRISPR tools, therapeutic targets, stem cell genomics, immunotherapy, metabolism-related genes and inflammatory bowel disease. In addition, the limitations and future development of CRISPR/Cas9 in colon cancer studies are reviewed. In conclusion, this article reviews the application of CRISPR-Cas9 gene editing technology in basic research, diagnosis and treatment of colon cancer.

Characterizing chromatin interactions of regulatory elements and nucleosome positions, using Hi-C, Micro-C, and promoter capture Micro-C

BACKGROUND

Regulatory elements such as promoters, enhancers, and insulators interact each other to mediate molecular processes. To capture chromatin interactions of regulatory elements, 3C-derived methods such as Hi-C and Micro-C are developed. Here, we generated and analyzed Hi-C, Micro-C, and promoter capture Micro-C datasets with different sequencing depths to study chromatin interactions of regulatory elements and nucleosome positions in human prostate cancer cells.

RESULTS

Compared to Hi-C, Micro-C identifies more high-resolution loops, including ones around structural variants. By evaluating the effect of sequencing depth, we revealed that more than 2 billion reads of Micro-C are needed to detect chromatin interactions at 1 kb resolution. Moreover, we found that deep-sequencing identifies additional long-range loops that are longer than 1 Mb in distance. Furthermore, we found that more than 50% of the loops are involved in insulators while less than 10% of the loops are promoter-enhancer loops. To comprehensively capture chromatin interactions that promoters are involved in, we performed promoter capture Micro-C. Promoter capture Micro-C identifies loops near promoters with a lower amount of sequencing reads. Sequencing of 160 million reads of promoter capture Micro-C resulted in reaching a plateau of identifying loops. However, there was still a subset of promoters that are not involved in loops even after deep-sequencing. By integrating Micro-C with NOMe-seq and ChIP-seq, we found that active promoters involved in loops have a more accessible region with lower levels of DNA methylation and more highly phased nucleosomes, compared to active promoters that are not involved in loops.

CONCLUSION

We determined the required sequencing depth for Micro-C and promoter capture Micro-C to generate high-resolution chromatin interaction maps and loops. We also investigated the effect of sequencing coverage of Hi-C, Micro-C, and promoter capture Micro-C on detecting chromatin loops. Our analyses suggest the presence of distinct regulatory element groups, which are differently involved in nucleosome positions and chromatin interactions. This study does not only provide valuable insights on understanding chromatin interactions of regulatory elements, but also present guidelines for designing research projects on chromatin interactions among regulatory elements.

Targeting of noncoding RNAs encoded by a novel MYC enhancers inhibits the proliferation of human hepatic carcinoma cells in vitro

The proto-oncogene MYC is important for development and cell growth, however, its abnormal regulation causes cancer. Recent studies identified distinct enhancers of MYC in various cancers, but any MYC enhancer(s) in hepatocellular carcinoma (HCC) remain(s) elusive. By analyzing H3K27ac enrichment and enhancer RNA (eRNA) expression in cultured HCC cells, we identified six putative MYC enhancer regions. Amongst these, two highly active enhancers, located ~ 800 kb downstream of the MYC gene, were identified by qRT-PCR and reporter assays. We functionally confirmed these enhancers by demonstrating a significantly reduced MYC expression and cell proliferation upon CRISPR/Cas9-based deletion and/or antisense oligonucleotide (ASO)-mediated inhibition. In conclusion, we identified potential MYC enhancers of HCC and propose that the associated eRNAs may be suitable targets for HCC treatment.

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

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.

Characterization of a functional endothelial super-enhancer that regulates ADAMTS18 and angiogenesis

Super-enhancers are clusters of enhancers associated with cell lineage. They can be powerful gene-regulators and may be useful in cell-type specific viral-vector development. Here, we have screened for endothelial super-enhancers and identified an enhancer from within a cluster that conferred 5–70-fold increase in transgene expression. Importantly, CRISPR/Cas9 deletion of enhancers demonstrated regulation of ADAMTS18, corresponding to evidence of chromatin contacts between these genomic regions. Cell division-related pathways were primarily affected by the enhancer deletions, which correlated with significant reduction in cell proliferation. Furthermore, we observed changes in angiogenesis-related genes consistent with the endothelial specificity of this SE. Indeed, deletion of the enhancers affected tube formation, resulting in reduced or shortened sprouts. The super-enhancer angiogenic role is at least partly due to its regulation of ADAMTS18, as siRNA knockdown of ADAMTS18 resulted in significantly shortened endothelial sprouts. Hence, functional characterization of a novel endothelial super-enhancer has revealed substantial downstream effects from single enhancer deletions and led to the discovery of the cis-target gene ADAMTS18 and its role in endothelial function.

Genomic editing of intronic enhancers unveils their role in fine-tuning tissue-specific gene expression in Arabidopsis thaliana

Enhancers located in introns are abundant and play a major role in the regulation of gene expression in mammalian species. By contrast, the functions of intronic enhancers in plants have largely been unexplored and only a handful of plant intronic enhancers have been reported. We performed a genome-wide prediction of intronic enhancers in Arabidopsis thaliana using open chromatin signatures based on DNase I sequencing. We identified 941 candidate intronic enhancers associated with 806 genes in seedling tissue and 1,271 intronic enhancers associated with 1,069 genes in floral tissue. We validated the function of 15 of 21 (71%) of the predicted intronic enhancers in transgenic assays using a reporter gene. We also created deletion lines of three intronic enhancers associated with two different genes using CRISPR/Cas. Deletion of these enhancers, which span key transcription factor binding sites, did not abolish gene expression but caused varying levels of transcriptional repression of their cognate genes. Remarkably, the transcriptional repression of the deletion lines occurred at specific developmental stages and resulted in distinct phenotypic effects on plant morphology and development. Clearly, these three intronic enhancers are important in fine-tuning tissue- and development-specific expression of their cognate genes.

CoBRA: Containerized Bioinformatics Workflow for Reproducible ChIP/ATAC-seq Analysis

Chromatin immunoprecipitation sequencing (ChIP-seq) and the Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) have become essential technologies to effectively measure protein–DNA interactions and chromatin accessibility. However, there is a need for a scalable and reproducible pipeline that incorporates proper normalization between samples, correction of copy number variations, and integration of new downstream analysis tools. Here we present Containerized Bioinformatics workflow for Reproducible ChIP/ATAC-seq Analysis (CoBRA), a modularized computational workflow which quantifies ChIP-seq and ATAC-seq peak regions and performs unsupervised and supervised analyses. CoBRA provides a comprehensive state-of-the-art ChIP-seq and ATAC-seq analysis pipeline that can be used by scientists with limited computational experience. This enables researchers to gain rapid insight into protein–DNA interactions and chromatin accessibility through sample clustering, differential peak calling, motif enrichment, comparison of sites to a reference database, and pathway analysis. CoBRA is publicly available online at https://bitbucket.org/cfce/cobra

Validation of prostate cancer risk variants rs10993994 and rs7098889 by CRISPR/Cas9 mediated genome editing

GWAS have identified numerous SNPs associated with prostate cancer risk. One such SNP is rs10993994. It is located in the β-microseminoprotein (MSMB) promoter region, mediates MSMB prostate secretion levels, and is linked to mRNA expression changes in both MSMB and the adjacent gene NCOA4. In addition, our previous work showed a second SNP, rs7098889, is in positive linkage disequilibrium with rs10993994 and associated with MSMB expression independent of rs10993994. Here, we generate a series of clones with single alleles removed by double guide RNA (gRNA) mediated CRISPR/Cas9 deletions, through which we demonstrate that each of these SNPs independently and greatly alters MSMB expression in an allele-specific manner. We further show that these SNPs have no substantial effect on the expression of NCOA4. These data demonstrate that a single SNP can have a large effect on gene expression and illustrate the importance of functional validation studies to deconvolute observed correlations. The method we have developed is generally applicable to test any SNP for which a relevant heterozygous cell line is available. AUTHOR SUMMARY: In pursuing the underlying biological mechanism of prostate cancer pathogenesis, scientists utilized the existence of common single nucleotide polymorphisms (SNPs) in the human genome as genetic markers to perform large scale genome wide association studies (GWAS) and have so far identified more than a hundred prostate cancer risk variants. Such variants provide an unbiased and systematic new venue to study the disease mechanism, and the next big challenge is to translate these genetic associations to the causal role of altered gene function in oncogenesis. The majority of these variants are waiting to be studied and lots of them may act in oncogenesis through gene expression regulation. To prove the concept, we took rs10993994 and its linked rs7098889 as an example and engineered single cell clones by allelic-specific CRISPR/Cas9 deletion to separate the effect of each allele. We observed that a single nucleotide difference would lead to surprisingly high level of MSMB gene expression change in a gene specific and cell-type specific manner. Our study strongly supports the notion that differential level of gene expression caused by risk variants and their associated genetic locus play a major role in oncogenesis and also highlights the importance of studying the function of MSMB encoded β-MSP in prostate cancer pathogenesis.

Signatures of co-deregulated genes and their transcriptional regulators in colorectal cancer

The deregulated genes in colorectal cancer (CRC) vary significantly across different studies. Thus, a systems biology approach is needed to identify the co-deregulated genes (co-DEGs), explore their molecular networks, and spot the major hub proteins within these networks. We reanalyzed 19 GEO gene expression profiles to identify and annotate CRC versus normal signatures, single-gene perturbation, and single-drug perturbation signatures. We identified the co-DEGs across different studies, their upstream regulating kinases and transcription factors (TFs). Connectivity Map was used to identify likely repurposing drugs against CRC within each group. The functional changes of the co-upregulated genes in the first category were mainly associated with negative regulation of transforming growth factor β production and glomerular epithelial cell differentiation; whereas the co-downregulated genes were enriched in cotranslational protein targeting to the membrane. We identified 17 hub proteins across the co-upregulated genes and 18 hub proteins across the co-downregulated genes, composed of well-known TFs (MYC, TCF3, PML) and kinases (CSNK2A1, CDK1/4, MAPK14), and validated most of them using GEPIA2 and HPA, but also through two signature gene lists composed of the co-up and co-downregulated genes. We further identified a list of repurposing drugs that can potentially target the co-DEGs in CRC, including camptothecin, neostigmine bromide, emetine, remoxipride, cephaeline, thioridazine, and omeprazole. Similar analyses were performed in the co-DEG signatures in single-gene or drug perturbation experiments in CRC. MYC, PML, CDKs, CSNK2A1, and MAPKs were common hub proteins among all studies. Overall, we identified the critical genes in CRC and we propose repurposing drugs that could be used against them.

Total Functional Score of Enhancer Elements Identifies Lineage-Specific Enhancers That Drive Differentiation of Pancreatic Cells

The differentiation of embryonic stem cells into various lineages is highly dependent on the chromatin state of the genome and patterns of gene expression. To identify lineage-specific enhancers driving the differentiation of progenitors into pancreatic cells, we used a previously described computational framework called Total Functional Score of Enhancer Elements (TFSEE), which integrates multiple genomic assays that probe both transcriptional and epigenomic states. First, we evaluated and compared TFSEE as an enhancer-calling algorithm with enhancers called using GRO-seq-defined enhancer transcripts (method 1) versus enhancers called using histone modification ChIP-seq data (method 2). Second, we used TFSEE to define the enhancer landscape and identify transcription factors (TFs) that maintain the multipotency of a subpopulation of endodermal stem cells during differentiation into pancreatic lineages. Collectively, our results demonstrate that TFSEE is a robust enhancer-calling algorithm that can be used to perform multilayer genomic data integration to uncover cell type-specific TFs that control lineage-specific enhancers.

H3K27ac-activated LINC00519 promotes lung squamous cell carcinoma progression by targeting miR-450b-5p/miR-515-5p/YAP1 axis

Objectives

Long non‐coding RNAs (lncRNAs) are extensively reported as participants in the biological process of diverse malignancies, including lung squamous cell carcinoma (LUSC). Long intergenic non‐protein coding RNA 519 (LINC00519) is identified as a novel lncRNA which has not yet been studied in cancers.

Materials and Methods

LINC00519 expression was detected by qRT‐PCR. The effect of LINC00519 on LUSC cellular activities was determined by in vitro and in vivo assays. Subcellular fractionation and FISH assays were conducted to identify the localization of LINC00519. The interaction between miR‐450b‐5p/miR‐515‐5p and LINC00519/YAP1 was verified by RIP, RNA pull‐down and luciferase reporter assays.

Results

Elevated level of LINC00519 was identified in LUSC tissues and cell lines. High LINC00519 level predicted unsatisfactory prognosis. Then, loss‐of‐function assays suggested the inhibitive role of silenced LINC00519 in cell proliferation, migration, invasion and tumour growth and promoting effect on cell apoptosis in LUSC. Mechanically, LINC00519 was activated by H3K27 acetylation (H3K27ac). Moreover, LINC00519 sponged miR‐450b‐5p and miR‐515‐5p to up‐regulate Yes1 associated transcriptional regulator (YAP1). Additionally, miR‐450b‐5p and miR‐515‐5p elicited anti‐carcinogenic effects in LUSC. Finally, rescue assays validated the effect of LINC00519‐miR‐450b‐5p‐miR‐515‐5p‐YAP1 axis in LUSC.

Conclusions

H3K27ac‐activated LINC00519 acts as a competing endogenous RNA (ceRNA) to promote LUSC progression by targeting miR‐450b‐5p/miR‐515‐5p/YAP1 axis.

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.

Functional characterization of two enhancers located downstream FOXP2

BACKGROUND

Mutations in the coding region of FOXP2 are known to cause speech and language impairment. However, it is not clear how dysregulation of the gene contributes to language deficit. Interestingly, microdeletions of the region downstream the gene have been associated with cognitive deficits.

METHODS

Here, we investigate changes in FOXP2 expression in the SK-N-MC neuroblastoma human cell line after deletion by CRISPR-Cas9 of two enhancers located downstream of the gene.

RESULTS

Deletion of any of these two functional enhancers downregulates FOXP2, but also upregulates the closest 3' gene MDFIC. Because this effect is not statistically significant in a HEK 293 cell line, derived from the human kidney, both enhancers might confer a tissue specific regulation to both genes. We have also found that the deletion of any of these enhancers downregulates six well-known FOXP2 target genes in the SK-N-MC cell line.

CONCLUSIONS

We expect these findings contribute to a deeper understanding of how FOXP2 and MDFIC are regulated to pace neuronal development supporting cognition, speech and language.

lncRNA PLAC2 activated by H3K27 acetylation promotes cell proliferation and invasion via the activation of Wnt/β‑catenin pathway in oral squamous cell carcinoma

As a new group of important effector molecules involved in multiple cancer types, including breast cancer, lung cancer and oral squamous cell carcinoma, long noncoding RNAs (lncRNAs) have attracted considerable attention recently. However, the underlying cause that induces the dysregulated lncRNAs in cancer remains poorly understood. In the present study, the regulatory model of the lncRNA placenta-specific protein 2 (PLAC2) upregulation in oral squamous cell carcinoma (OSCC) was investigated and its biological functions in OSCC malignant progression was identified. A reverse transcription-quantitative polymerase chain reaction assay identified that PLAC2 is upregulated in OSCC cell lines and primary tissue samples. Furthermore, bioinformatic analysis followed by chromatin immunoprecipitation verified an enriched histone H3 on lysine 27 (H3K27) acetylation (H3K27ac) at the promoter region of the PLAC2 gene. Knockdown of cAMP-response element binding protein-binding protein (CBP) significantly reduced the enrichment level of H3K27ac, and thereby induced a decreased expression of PLAC2. Functionally, overexpression of PLAC2 promotes OSCC cell proliferation, migration and invasion, whereas knockdown of PLAC2 exerted an opposite effect. Furthermore, the Wnt/β-catenin signaling pathway was activated by PLAC2 and mediated the PLAC2-induced malignant progress of OSCC. In conclusion, the present results indicated that lncRNA PLAC2 is transcriptionally activated by H3K27ac modification at the promoter region in OSCC, and promotes cell growth and metastasis via activating Wnt/β-catenin signaling pathway. Therefore, PLAC2 may serve as a promising biomarker for OSCC prognosis and therapy.

Regulatory SNPs and their widespread effects on the transcriptome

Currently, it is generally accepted that the cis-acting effects of noncoding variants on gene expression are a major factor for phenotypic variation in complex traits and disease susceptibility. Meanwhile, the protein products of many target genes for the identified cis-regulatory variants (rSNPs) are regulatory molecules themselves (transcription factors, effectors, components of signal transduction pathways, etc.), which implies dramatic downstream effects of these variations on complex gene networks. Here, we brief the results of recent most comprehensive studies on the role of rSNPs in transcriptional regulation across the genome.

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.

Development of Toolboxes for Precision Genome/Epigenome Editing and Imaging of Epigenetics

Zinc finger (ZF) proteins are composed of repeated ββα modules and coordinate a zinc ion. ZF domains recognizing specific DNA target sequences can be substituted for the binding domains of various DNA-modifying enzymes to create designer nucleases, recombinases, and methyltransferases with programmable sequence specificity. Enzymatic genome editing and modification can be applied to many fields of basic research and medicine. The recent development of new platforms using transcription activator-like effector (TALE) proteins or the CRISPR-Cas9 system has expanded the range of possibilities for genome-editing technologies. In addition, these DNA binding domains can also be utilized to build a toolbox for epigenetic controls by fusing them with protein- or DNA-modifying enzymes. Here, our research on epigenome editing including the development of artificial zinc finger recombinase (ZFR), split DNA methyltransferase, and fluorescence imaging of histone proteins by ZIP tag-probe system is introduced. Advances in the ZF, TALE, and CRISPR-Cas9 platforms have paved the way for the next generation of genome/epigenome engineering approaches.

NF-κB-Dependent Lymphoid Enhancer Co-option Promotes Renal Carcinoma Metastasis

Metastases, the spread of cancer cells to distant organs, cause the majority of cancer-related deaths. Few metastasis-specific driver mutations have been identified, suggesting aberrant gene regulation as a source of metastatic traits. However, how metastatic gene expression programs arise is poorly understood. Here, using human-derived metastasis models of renal cancer, we identify transcriptional enhancers that promote metastatic carcinoma progression. Specific enhancers and enhancer clusters are activated in metastatic cancer cell populations, and the associated gene expression patterns are predictive of poor patient outcome in clinical samples. We find that the renal cancer metastasis-associated enhancer complement consists of multiple coactivated tissue-specific enhancer modules. Specifically, we identify and functionally characterize a coregulatory enhancer cluster, activated by the renal cancer driver HIF2A and an NF-κB-driven lymphoid element, as a mediator of metastasis in vivo We conclude that oncogenic pathways can acquire metastatic phenotypes through cross-lineage co-option of physiologic epigenetic enhancer states.Significance: Renal cancer is associated with significant mortality due to metastasis. We show that in metastatic renal cancer, functionally important metastasis genes are activated via co-option of gene regulatory enhancer modules from distant developmental lineages, thus providing clues to the origins of metastatic cancer. Cancer Discov; 8(7); 850-65. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 781.

Methods for Scarless, Selection-Free Generation of Human Cells and Allele-Specific Functional Analysis of Disease-Associated SNPs and Variants of Uncertain Significance

With the continued emergence of risk loci from Genome-Wide Association studies and variants of uncertain significance identified from patient sequencing, better methods are required to translate these human genetic findings into improvements in public health. Here we combine CRISPR/Cas9 gene editing with an innovative high-throughput genotyping pipeline utilizing KASP (Kompetitive Allele-Specific PCR) genotyping technology to create scarless isogenic cell models of cancer variants in ~1 month. We successfully modeled two novel variants previously identified by our lab in the PALB2 gene in HEK239 cells, resulting in isogenic cells representing all three genotypes for both variants. We also modeled a known functional risk SNP of colorectal cancer, rs6983267, in HCT-116 cells. Cells with extremely low levels of gene editing could still be identified and isolated using this approach. We also introduce a novel molecular assay, ChIPnQASO (Chromatin Immunoprecipitation and Quantitative Allele-Specific Occupation), which uses the same technology to reveal allele-specific function of these variants at the DNA-protein interaction level. We demonstrated preferential binding of the transcription factor TCF7L2 to the rs6983267 risk allele over the non-risk. Our pipeline provides a platform for functional variant discovery and validation that is accessible and broadly applicable for the progression of efforts towards precision medicine.

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.

Identification of activated enhancers and linked transcription factors in breast, prostate, and kidney tumors by tracing enhancer networks using epigenetic traits

Background

Although technological advances now allow increased tumor profiling, a detailed understanding of the mechanisms leading to the development of different cancers remains elusive. Our approach toward understanding the molecular events that lead to cancer is to characterize changes in transcriptional regulatory networks between normal and tumor tissue. Because enhancer activity is thought to be critical in regulating cell fate decisions, we have focused our studies on distal regulatory elements and transcription factors that bind to these elements.

Results

Using DNA methylation data, we identified more than 25,000 enhancers that are differentially activated in breast, prostate, and kidney tumor tissues, as compared to normal tissues. We then developed an analytical approach called Tracing Enhancer Networks using Epigenetic Traits that correlates DNA methylation levels at enhancers with gene expression to identify more than 800,000 genome-wide links from enhancers to genes and from genes to enhancers. We found more than 1200 transcription factors to be involved in these tumor-specific enhancer networks. We further characterized several transcription factors linked to a large number of enhancers in each tumor type, including GATA3 in non-basal breast tumors, HOXC6 and DLX1 in prostate tumors, and ZNF395 in kidney tumors. We showed that HOXC6 and DLX1 are associated with different clusters of prostate tumor-specific enhancers and confer distinct transcriptomic changes upon knockdown in C42B prostate cancer cells. We also discovered de novo motifs enriched in enhancers linked to ZNF395 in kidney tumors.

Conclusions

Our studies characterized tumor-specific enhancers and revealed key transcription factors involved in enhancer networks for specific tumor types and subgroups. Our findings, which include a large set of identified enhancers and transcription factors linked to those enhancers in breast, prostate, and kidney cancers, will facilitate understanding of enhancer networks and mechanisms leading to the development of these cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0102-4) contains supplementary material, which is available to authorized users.

The role of enhancers in cancer

Enhancer elements function as the logic gates of the genetic regulatory circuitry. One of their most important functions is the integration of extracellular signals with intracellular cell fate information to generate cell type-specific transcriptional responses. Mutations occurring in cancer often misregulate enhancers that normally control the signal-dependent expression of growth-related genes. This misregulation can result from trans-acting mechanisms, such as activation of the transcription factors or epigenetic regulators that control enhancer activity, or can be caused in cis by direct mutations that alter the activity of the enhancer or its target gene specificity. These processes can generate tumour type-specific super-enhancers and establish a 'locked' gene regulatory state that drives the uncontrolled proliferation of cancer cells. Here, we review the role of enhancers in cancer, and their potential as therapeutic targets.

A CRISPR Path to Engineering New Genetic Mouse Models for Cardiovascular Research

Previous efforts to target the mouse genome for the addition, subtraction, or substitution of biologically informative sequences required complex vector design and a series of arduous steps only a handful of laboratories could master. The facile and inexpensive clustered regularly interspaced short palindromic repeats (CRISPR) method has now superseded traditional means of genome modification such that virtually any laboratory can quickly assemble reagents for developing new mouse models for cardiovascular research. Here, we briefly review the history of CRISPR in prokaryotes, highlighting major discoveries leading to its formulation for genome modification in the animal kingdom. Core components of CRISPR technology are reviewed and updated. Practical pointers for 2-component and 3-component CRISPR editing are summarized with many applications in mice including frameshift mutations, deletion of enhancers and noncoding genes, nucleotide substitution of protein-coding and gene regulatory sequences, incorporation of loxP sites for conditional gene inactivation, and epitope tag integration. Genotyping strategies are presented and topics of genetic mosaicism and inadvertent targeting discussed. Finally, clinical applications and ethical considerations are addressed as the biomedical community eagerly embraces this astonishing innovation in genome editing to tackle previously intractable questions.

The MYC 3' Wnt-Responsive Element Drives Oncogenic MYC Expression in Human Colorectal Cancer Cells

Mutations in components of the Wnt/β-catenin signaling pathway drive colorectal cancer (CRC) by deregulating expression of downstream target genes including the c-MYC proto-oncogene (MYC). The critical regulatory DNA enhancer elements that control oncogenic MYC expression in CRC have yet to be fully elucidated. In previous reports, we correlated T-cell factor (TCF) and β-catenin binding to the MYC 3′ Wnt responsive DNA element (MYC 3′ WRE) with MYC expression in HCT116 cells. Here we used CRISPR/Cas9 to determine whether this element is a critical driver of MYC. We isolated a clonal population of cells that contained a deletion of a single TCF binding element (TBE) within the MYC 3′ WRE. This deletion reduced TCF/β-catenin binding to this regulatory element and decreased MYC expression. Using RNA-Seq analysis, we found altered expression of genes that regulate metabolic processes, many of which are known MYC target genes. We found that 3′ WRE-Mut cells displayed a reduced proliferative capacity, diminished clonogenic growth, and a decreased potential to form tumors in vivo. These findings indicate that the MYC 3′ WRE is a critical driver of oncogenic MYC expression and suggest that this element may serve as a therapeutic target for CRC.

4C-seq revealed long-range interactions of a functional enhancer at the 8q24 prostate cancer risk locus

Genome-wide association studies (GWAS) have identified >100 independent susceptibility loci for prostate cancer, including the hot spot at 8q24. However, how genetic variants at this locus confer disease risk hasn’t been fully characterized. Using circularized chromosome conformation capture (4C) coupled with next-generation sequencing and an enhancer at 8q24 as “bait”, we identified genome-wide partners interacting with this enhancer in cell lines LNCaP and C4-2B. These 4C-identified regions are distributed in open nuclear compartments, featuring active histone marks (H3K4me1, H3K4me2 and H3K27Ac). Transcription factors NKX3-1, FOXA1 and AR (androgen receptor) tend to occupy these 4C regions. We identified genes located at the interacting regions, and found them linked to positive regulation of mesenchymal cell proliferation in LNCaP and C4-2B, and several pathways (TGF beta signaling pathway in LNCaP and p53 pathway in C4-2B). Common genes (e.g. MYC and POU5F1B) were identified in both prostate cancer cell lines. However, each cell line also had exclusive genes (e.g. ELAC2 and PTEN in LNCaP and BRCA2 and ZFHX3 in C4-2B). In addition, BCL-2 identified in C4-2B might contribute to the progression of androgen-refractory prostate cancer. Overall, our work reveals key genes and pathways involved in prostate cancer onset and progression.