Mice lacking a Myc enhancer that includes human SNP rs6983267 are resistant to intestinal tumors

CCAT1 and CCAT2 long noncoding RNAs, located within the 8q.24.21 'gene desert', serve as important prognostic biomarkers in colorectal cancer

Background

8q24.21 is a frequently amplified genomic region in colorectal cancer (CRC). This region is often referred to as a 'gene desert' due to lack of any important protein-coding genes, highlighting the potential role of noncoding RNAs, including long noncoding RNAs (lncRNAs) located around the proto-oncogene MYC. In this study, we have firstly evaluated the clinical significance of altered expression of lncRNAs mapped to this genomic locus in CRC.

Patients and methods

A total of 300 tissues, including 280 CRC and 20 adjacent normal mucosa specimens were evaluated for the expression of 12 lncRNAs using qRT-PCR assays. We analyzed the associations between lncRNA expression and various clinicopathological features, as well as with recurrence free survival (RFS) and overall survival (OS) in two independent cohorts.

Results

The expression of CCAT1, CCAT1-L, CCAT2, PVT1, and CASC19 were elevated in cancer tissues (P = 0.039, <0.001, 0.018, <0.001, 0.002, respectively). Among these, high expression of CCAT1 and CCAT2 was significantly associated with poor RFS (P = 0.049 and 0.022, respectively) and OS (P = 0.028 and 0.015, respectively). These results were validated in an independent patient cohort, in which combined expression of CCAT1 and CCAT2 expression was significantly associated with a poor RFS (HR:2.60, 95% confidence interval [CI]: 1.04-6.06, P = 0.042) and a poor OS (HR:8.38, 95%CI: 2.68-37.0, P < 0.001). We established a RFS prediction model which revealed that combined expression of CCAT1, CCAT2, and carcinoembryonic antigen was a significant determinant for efficiently predicting RFS in stage II (P = 0.034) and stage III (P = 0.001) CRC patients.

Conclusions

Several lncRNAs located in 8q24.21 locus are highly over-expressed in CRC. High expression of CCAT1 and CCAT2 significantly associates with poor RFS and OS. The expression of these two lncRNAs independently, or in combination, serves as important prognostic biomarkers in CRC.

Emerging applications of genome-editing technology to examine functionality of GWAS-associated variants for complex traits

Over the last decade, genome-wide association studies (GWAS) have propelled the discovery of thousands of loci associated with complex diseases. The focus is now turning toward the function of these association signals, determining the causal variant(s) among those in strong linkage disequilibrium, and identifying their underlying mechanisms, such as long-range gene regulation. Genome-editing techniques utilizing zinc-finger nucleases (ZFN), transcription activator-like effector nucleases (TALENs), and clustered regularly-interspaced short palindromic repeats with Cas9 nuclease (CRISPR-Cas9) are becoming the tools of choice to establish functionality for these variants, due to the ability to assess effects of single variants in vivo. This review will discuss examples of how these technologies have begun to aid functional analysis of GWAS loci for complex traits such as cardiovascular disease, Type 2 diabetes, cancer, obesity, and autoimmune disease. We focus on analysis of variants occurring within noncoding genomic regions, as these comprise the majority of GWAS variants, providing the greatest challenges to determining functionality, and compare editing strategies that provide different levels of evidence for variant functionality. The review describes molecular insights into some of these potentially causal variants and how these may relate to the pathology of the trait and look toward future directions for these technologies in post-GWAS analysis, such as base-editing.

Common Genetic Variation and Susceptibility to Ovarian Cancer: Current Insights and Future Directions

In this review, we summarize current progress in the genetic epidemiology of epithelial ovarian cancer (EOC), focusing exclusively on elucidating the role of common germline genetic variation in conferring susceptibility to EOC. We provide an overview of the more than 30 EOC risk loci identified to date by genome-wide association studies (GWAS) and describe the contribution of large-scale, cross-cancer type, custom genotyping projects, such as the OncoArray and the Collaborative Oncological Gene-Environment Study, to locus discovery and replication. We discuss the histotype-specific nature of these EOC risk loci, pleiotropy, or overlapping genetic effects between EOC and other hormone-related cancer types, and the application of findings to polygenic risk prediction for EOC. The second part of the article offers a concise review of primarily laboratory-based studies that have led to the identification of several putative EOC susceptibility genes using common variants at the known EOC risk loci as starting points. More global biological insights emerging from network- and pathway-based analyses of GWAS for EOC susceptibility are also highlighted. Finally, we delve into potential future directions, including the need to identify EOC risk loci in non-European populations and the next generation of GWAS functional studies that are likely to involve genome editing to establish the cell type-specific carcinogenic effects of EOC risk variants Cancer Epidemiol Biomarkers Prev; 27(4); 395-404. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Pan-cancer Alterations of the MYC Oncogene and Its Proximal Network across the Cancer Genome Atlas

Although the MYC oncogene has been implicated in cancer, a systematic assessment of alterations of MYC, related transcription factors, and co-regulatory proteins, forming the proximal MYC network (PMN), across human cancers is lacking. Using computational approaches, we define genomic and proteomic features associated with MYC and the PMN across the 33 cancers of The Cancer Genome Atlas. Pan-cancer, 28% of all samples had at least one of the MYC paralogs amplified. In contrast, the MYC antagonists MGA and MNT were the most frequently mutated or deleted members, proposing a role as tumor suppressors. MYC alterations were mutually exclusive with PIK3CA, PTEN, APC, or BRAF alterations, suggesting that MYC is a distinct oncogenic driver. Expression analysis revealed MYC-associated pathways in tumor subtypes, such as immune response and growth factor signaling; chromatin, translation, and DNA replication/repair were conserved pan-cancer. This analysis reveals insights into MYC biology and is a reference for biomarkers and therapeutics for cancers with alterations of MYC or the PMN.

Deregulation of polycomb repressor complex 1 modifier AUTS2 in T-cell leukemia

Recently, we identified deregulated expression of the B-cell specific transcription factor MEF2C in T-cell acute lymphoid leukemia (T-ALL). Here, we performed sequence analysis of a regulatory upstream section of MEF2C in T-ALL cell lines which, however, proved devoid of mutations. Unexpectedly, we found strong conservation between the regulatory upstream region of MEF2C (located at chromosomal band 5q14) and an intergenic stretch at 7q11 located between STAG3L4 and AUTS2, covering nearly 20 kb. While the non-coding gene STAG3L4 was inconspicuously expressed, AUTS2 was aberrantly upregulated in 6% of T-ALL patients (public dataset GSE42038) and in 3/24 T-ALL cell lines, two of which represented very immature differentiation stages. AUTS2 expression was higher in normal B-cells than in T-cells, indicating lineage-specific activity in lymphopoiesis. While excluding chromosomal aberrations, examinations of AUTS2 transcriptional regulation in T-ALL cells revealed activation by IL7-IL7R-STAT5-signalling and MEF2C. AUTS2 protein has been shown to interact with polycomb repressor complex 1 subtype 5 (PRC1.5), transforming this particular complex into an activator. Accordingly, expression profiling and functional analyses demonstrated that AUTS2 activated while PCGF5 repressed transcription of NKL homeobox gene MSX1 in T-ALL cells. Forced expression and pharmacological inhibition of EZH2 in addition to H3K27me3 analysis indicated that PRC2 repressed MSX1 as well. Taken together, we found that AUTS2 and MEF2C, despite lying on different chromosomes, share strikingly similar regulatory upstream regions and aberrant expression in T-ALL subsets. Our data implicate chromatin complexes PRC1/AUTS2 and PRC2 in a gene network in T-ALL regulating early lymphoid differentiation.

Wnt signalling is a bi-directional vulnerability of cancer cells

Wnt signalling is involved in the formation, metastasis and relapse of a wide array of cancers. However, there is ongoing debate as to whether activation or inhibition of the pathway holds the most promise as a therapeutic treatment for cancer, with conflicting evidence from a variety of tumour types. We show that Wnt/β-catenin signalling is a bi-directional vulnerability of neuroblastoma, malignant melanoma and colorectal cancer, with hyper-activation or repression of the pathway both representing a promising therapeutic strategy, even within the same cancer type. Hyper-activation directs cancer cells to undergo apoptosis, even in cells oncogenically driven by β-catenin. Wnt inhibition blocks proliferation of cancer cells and promotes neuroblastoma differentiation. Wnt and retinoic acid co-treatments synergise, representing a promising combination treatment for MYCN-amplified neuroblastoma. Additionally, we report novel cross-talks between MYCN and β-catenin signalling, which repress normal β-catenin mediated transcriptional regulation. A β-catenin target gene signature could predict patient outcome, as could the expression level of its DNA binding partners, the TCF/LEFs. This β-catenin signature provides a tool to identify neuroblastoma patients likely to benefit from Wnt-directed therapy. Taken together, we show that Wnt/β-catenin signalling is a bi-directional vulnerability of a number of cancer entities, and potentially a more broadly conserved feature of malignant cells.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Molecular and physiological roles of the adaptor protein CARD9 in immunity

CARD9 is a caspase recruitment domain-containing signaling protein that plays a critical role in innate and adaptive immunity. It has been widely demonstrated that CARD9 adaptor allows pattern recognition receptors to induce NF-κB and MAPK activation, which initiates a “downstream” inflammation cytokine cascade and provides effective protection against microbial invasion, especially fungal infection. Here our aim is to update existing paradigms and summarize the most recent findings on the CARD9 signaling pathway, revealing significant mechanistic insights into the pathogenesis of CARD9 deficiency. We also discuss the effect of CARD9 genetic mutations on the in vivo immune response, and highlight clinical advances in non-infection inflammation.

A Myc enhancer cluster regulates normal and leukaemic haematopoietic stem cell hierarchies

The transcription factor Myc is essential for the regulation of haematopoietic stem cells and progenitors and has a critical function in haematopoietic malignancies. Here we show that an evolutionarily conserved region located 1.7 megabases downstream of the Myc gene that has previously been labelled as a 'super-enhancer' is essential for the regulation of Myc expression levels in both normal haematopoietic and leukaemic stem cell hierarchies in mice and humans. Deletion of this region in mice leads to a complete loss of Myc expression in haematopoietic stem cells and progenitors. This caused an accumulation of differentiation-arrested multipotent progenitors and loss of myeloid and B cells, mimicking the phenotype caused by Mx1-Cre-mediated conditional deletion of the Myc gene in haematopoietic stem cells. This super-enhancer comprises multiple enhancer modules with selective activity that recruits a compendium of transcription factors, including GFI1b, RUNX1 and MYB. Analysis of mice carrying deletions of individual enhancer modules suggests that specific Myc expression levels throughout most of the haematopoietic hierarchy are controlled by the combinatorial and additive activity of individual enhancer modules, which collectively function as a 'blood enhancer cluster' (BENC). We show that BENC is also essential for the maintenance of MLL-AF9-driven leukaemia in mice. Furthermore, a BENC module, which controls Myc expression in mouse haematopoietic stem cells and progenitors, shows increased chromatin accessibility in human acute myeloid leukaemia stem cells compared to blasts. This difference correlates with MYC expression and patient outcome. We propose that clusters of enhancers, such as BENC, form highly combinatorial systems that allow precise control of gene expression across normal cellular hierarchies and which also can be hijacked in malignancies.

Association of rs6983267 at 8q24, HULC rs7763881 polymorphisms and serum lncRNAs CCAT2 and HULC with colorectal cancer in Egyptian patients

The impact of HULC rs7763881 on colorectal cancer (CRC) susceptibility is not yet known. Also, the biological function of the cancer-related rs6983267 remains unclear. We investigated the association of these SNPs with the risk of CRC and adenomatous polyps (AP), their correlation with CCAT2 and HULC expression, and the potential of serum CCAT2 and HULC as biomarkers for CRC. 120 CRC patients, 30 AP patients, and 96 healthy controls were included. Genotyping and serum lncRNAs were assayed by qPCR. Studied SNPs were not associated with AP susceptibility. rs6983267 GG was associated with increased CRC risk, whereas rs7763881 AC was protective. rs7763881 and rs6983267 CT haplotype was protective. Serum CCAT2 and HULC were upregulated in CRC and AP patients versus controls and discriminated these groups by ROC analysis. rs6983267 GG and rs7763881 AA patients demonstrated higher serum CCAT2 and HULC compared with GT/TT and AC, respectively. rs6983267 and serum HULC predicted CRC diagnosis among non-CRC groups (AP + controls) by multivariate analysis. Studied SNPs or serum long noncoding RNAs weren't correlated with nodal or distant metastasis. In conclusion, rs6983267 and rs7763881 are potential genetic markers of CRC predisposition and correlate with serum CCAT2 and HULC, two novel potential non-invasive diagnostic biomarkers for CRC.

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.

Genome-wide association studies of cancer: current insights and future perspectives

Genome-wide association studies (GWAS) provide an agnostic approach for investigating the genetic basis of complex diseases. In oncology, GWAS of nearly all common malignancies have been performed, and over 450 genetic variants associated with increased risks have been identified. As well as revealing novel pathways important in carcinogenesis, these studies have shown that common genetic variation contributes substantially to the heritable risk of many common cancers. The clinical application of GWAS is starting to provide opportunities for drug discovery and repositioning as well as for cancer prevention. However, deciphering the functional and biological basis of associations is challenging and is in part a barrier to fully unlocking the potential of GWAS.

To Wnt or Lose: The Missing Non-Coding Linc in Colorectal Cancer

Colorectal cancer (CRC) is the third most frequent cancer and one of the leading causes for cancer-related mortality. Aberrant activation of the Wnt signaling is an essential initiating factor in colon carcinogenesis, and a driving force of CRC progression. Recently, long non-coding RNAs (lncRNAs) have emerged as significant players in CRC pathogenesis through diversified mechanisms. Although both Wnt signaling and lncRNAs represent interesting research areas for CRC, an effort of directly connecting these two areas is lacking. To fill in the knowledge gap, we focus on the reported findings of lncRNAs that regulate Wnt signaling or essential Wnt signaling targets. These include several newly discovered lncRNAs originated from the amplified cancer-associated chromosome 8q24 region that surrounds the essential Wnt target MYC gene, lncRNAs reported to be involved in CRC stem cells, and several individual lncRNAs connected to Wnt signaling through other mechanisms. This review will provide essential information that assists in understanding the missing link of lncRNAs to the classical Wnt signaling in CRC.

Contrasting effects of an Mdm2 functional polymorphism on tumor phenotypes

MDM2, an E3 ubiquitin ligase, is a potent inhibitor of the p53 tumor suppressor and is elevated in many human cancers that retain wild-type p53. MDM2 SNP309G is a functional polymorphism that results in elevated levels of MDM2 (due to enhanced SP1 binding to the MDM2 promoter) thus decreasing p53 activity. Mdm2^SNP309G/G mice are more prone to spontaneous tumor formation than Mdm2^SNP309T/T mice, providing direct evidence for the impact of this SNP in tumor development. We asked whether environmental factors impact SNP309G function and show that SNP309G cooperates with ionizing radiation to exacerbate tumor development. Surprisingly, ultraviolet B light or Benzo(a)pyrene exposure of skin shows that SNP309G allele actually protects against squamous cell carcinoma susceptibility. These contrasting differences led us to interrogate the mechanism by which Mdm2 SNP309 regulates tumor susceptibility in a tissue-specific manner. Although basal Mdm2 levels were significantly higher in most tissues in Mdm2^SNP309G/G mice compared with Mdm2^SNP309T/T mice, they were significantly lower in Mdm2^SNP309G/G keratinocytes, the cell-type susceptible to squamous cell carcinoma. The assessment of potential transcriptional regulators in ENCODE ChIP-seq database identified transcriptional repressor E2F6 as a possible negative regulator of MDM2 expression. Our data show that E2F6 suppresses Mdm2 expression in cells harboring the SNP309G allele but not the SNP309T allele. Thus, Mdm2 SNP309G exhibits tissue-specific regulation and differentially impacts cancer risk.

A Genetic Variant Associated with Five Vascular Diseases Is a Distal Regulator of Endothelin-1 Gene Expression

Genome-wide association studies (GWASs) implicate the PHACTR1 locus (6p24) in risk for five vascular diseases, including coronary artery disease, migraine headache, cervical artery dissection, fibromuscular dysplasia, and hypertension. Through genetic fine mapping, we prioritized rs9349379, a common SNP in the third intron of the PHACTR1 gene, as the putative causal variant. Epigenomic data from human tissue revealed an enhancer signature at rs9349379 exclusively in aorta, suggesting a regulatory function for this SNP in the vasculature. CRISPR-edited stem cell-derived endothelial cells demonstrate rs9349379 regulates expression of endothelin 1 (EDN1), a gene located 600 kb upstream of PHACTR1. The known physiologic effects of EDN1 on the vasculature may explain the pattern of risk for the five associated diseases. Overall, these data illustrate the integration of genetic, phenotypic, and epigenetic analysis to identify the biologic mechanism by which a common, non-coding variant can distally regulate a gene and contribute to the pathogenesis of multiple vascular diseases.

Association between 8q24 rs6983267 polymorphism and cancer susceptibility: a meta-analysis involving 170,737 subjects

Published data on the association between 8q24 rs6983267 polymorphism and cancer risk are inconsistent. Thus, we conducted a meta-analysis to evaluate the relationship between rs6983267 polymorphism and cancer risk. We searched on PubMed, EMBASE, Web of Science and China National Knowledge Infrastructure (CNKI) up to November 1, 2016 for relevant studies. Odds ratios (ORs) and 95% confidence intervals (CIs) were used to estimate the strength of this association. We included 78 case-control studies with a total of 73,996 cases and 96,741 controls in this meta-analysis. The pooled results showed that rs6983267 polymorphism was significantly associated with increased risk of overall cancer in all genetic models (dominant model: OR = 1.19, 95% CI = 1.13–1.26; recessive model: OR = 1.19, 95% CI = 1.14–1.25; homozygous model: OR= 1.31, 95% CI = 1.23–1.40; heterozygous model: OR = 1.14, 95% CI = 1.10–1.19; allelic model: OR = 1.14, 95% CI = 1.11–1.18). Stratified analyses indicated that rs6983267 significantly increased the risk of colorectal cancer in Caucasians, prostate cancer in Caucasians and Asians, thyroid cancer in Caucasians and lung cancer in Asians. When studies were stratified by study quality, source of controls and genotyping method, significant associations were especially found in the high quality studies, the publication-based studies, the hospital-based studies, and the PCR-RFLP studies. Additional well-designed studies with large samples should be performed to validate our results.

MYC Modulation around the CDK2/p27/SKP2 Axis

MYC is a pleiotropic transcription factor that controls a number of fundamental cellular processes required for the proliferation and survival of normal and malignant cells, including the cell cycle. MYC interacts with several central cell cycle regulators that control the balance between cell cycle progression and temporary or permanent cell cycle arrest (cellular senescence). Among these are the cyclin E/A/cyclin-dependent kinase 2 (CDK2) complexes, the CDK inhibitor p27^KIP1 (p27) and the E3 ubiquitin ligase component S-phase kinase-associated protein 2 (SKP2), which control each other by forming a triangular network. MYC is engaged in bidirectional crosstalk with each of these players; while MYC regulates their expression and/or activity, these factors in turn modulate MYC through protein interactions and post-translational modifications including phosphorylation and ubiquitylation, impacting on MYC's transcriptional output on genes involved in cell cycle progression and senescence. Here we elaborate on these network interactions with MYC and their impact on transcription, cell cycle, replication and stress signaling, and on the role of other players interconnected to this network, such as CDK1, the retinoblastoma protein (pRB), protein phosphatase 2A (PP2A), the F-box proteins FBXW7 and FBXO28, the RAS oncoprotein and the ubiquitin/proteasome system. Finally, we describe how the MYC/CDK2/p27/SKP2 axis impacts on tumor development and discuss possible ways to interfere therapeutically with this system to improve cancer treatment.

Genomic Insight into the Role of lncRNA in Cancer Susceptibility

With the development of advanced genomic methods, a large amount of long non-coding RNAs (lncRNAs) has been found to be important for cancer initiation and progression. Given that most of the genome-wide association study (GWAS)-identified cancer risk SNPs are located in the noncoding region, the expression and function of lncRNAs are more likely to be affected by the SNPs. The SNPs may affect the expression of lncRNAs directly through disrupting the binding of transcription factors or indirectly by affecting the expression of regulatory factors. Moreover, SNPs may disrupt the interaction between lncRNAs and other RNAs or proteins. Unveiling the relationship of lncRNA, protein-coding genes, transcription factors and miRNAs from the angle of genomics will improve the accuracy of disease prediction and help find new therapeutic targets.

Enhancers not required

Laboratory mice with over half a megabase of DNA upstream of their Myc gene removed still thrive in the absence of stress.

Mice deficient of Myc super-enhancer region reveal differential control mechanism between normal and pathological growth

The gene desert upstream of the MYC oncogene on chromosome 8q24 contains susceptibility loci for several major forms of human cancer. The region shows high conservation between human and mouse and contains multiple MYC enhancers that are activated in tumor cells. However, the role of this region in normal development has not been addressed. Here we show that a 538 kb deletion of the entire MYC upstream super-enhancer region in mice results in 50% to 80% decrease in Myc expression in multiple tissues. The mice are viable and show no overt phenotype. However, they are resistant to tumorigenesis, and most normal cells isolated from them grow slowly in culture. These results reveal that only cells whose MYC activity is increased by serum or oncogenic driver mutations depend on the 8q24 super-enhancer region, and indicate that targeting the activity of this element is a promising strategy of cancer chemoprevention and therapy.

Wnt target genes and where to find them

Wnt/β-catenin signaling is highly conserved throughout metazoans, is required for numerous essential events in development, and serves as a stem cell niche signal in many contexts. Misregulation of the pathway is linked to several human pathologies, most notably cancer. Wnt stimulation results in stabilization and nuclear import of β-catenin, which then acts as a transcriptional co-activator. Transcription factors of the T-cell family (TCF) are the best-characterized nuclear binding partners of β-catenin and mediators of Wnt gene regulation. This review provides an update on what is known about the transcriptional activation of Wnt target genes, highlighting recent work that modifies the conventional model. Wnt/β-catenin signaling regulates genes in a highly context-dependent manner, and the role of other signaling pathways and TCF co-factors in this process will be discussed. Understanding Wnt gene regulation has served to elucidate many biological roles of the pathway, and we will use examples from stem cell biology, metabolism, and evolution to illustrate some of the rich Wnt biology that has been uncovered.

Altered TERT promoter and other genomic regulatory elements: occurrence and impact

Study of genetic alterations, inherited or acquired, that increase the risk or drive cancers and many other diseases had remained mostly confined to coding sequences of the human genome. Data from genome wide associations studies, development of the Encyclopedia of DNA Elements (ENCODE), and a spurt in detection of driver somatic mutations have shifted focus towards noncoding regions of the human genome. The majority of genetic variants robustly associated with cancers and other syndromes identified through genome wide studies are located within noncoding regulatory regions of the genome. Genome wide techniques have put an emphasis on the role of three-dimensional chromosomal structures and cis-acting elements in regulations of different genes. The variants within noncoding genomic regions can potentially alter a number of regulatory elements including promoters, enhancers, insulators, noncoding long RNAs and others that affect cancers and various diseases through altered expression of critical genes. With effect of genetic alterations within regulatory elements dependent on other partner molecules like transcription factors and histone marks, an understanding of such modifications can potentially identify extended therapeutic targets. That concept has been augmented by the detection of driver somatic noncoding mutations within the promoter region of the telomerase reverse transcriptase (TERT) gene in different cancers. The acquired somatic noncoding mutations within different regulatory elements are now being reported in different cancers with an increased regularity. In this review we discuss the occurrence and impact of germline and somatic alterations within the TERT promoter and other genomic regulatory elements.

Small genomic insertions form enhancers that misregulate oncogenes

The non-coding regions of tumour cell genomes harbour a considerable fraction of total DNA sequence variation, but the functional contribution of these variants to tumorigenesis is ill-defined. Among these non-coding variants, somatic insertions are among the least well characterized due to challenges with interpreting short-read DNA sequences. Here, using a combination of Chip-seq to enrich enhancer DNA and a computational approach with multiple DNA alignment procedures, we identify enhancer-associated small insertion variants. Among the 102 tumour cell genomes we analyse, small insertions are frequently observed in enhancer DNA sequences near known oncogenes. Further study of one insertion, somatically acquired in primary leukaemia tumour genomes, reveals that it nucleates formation of an active enhancer that drives expression of the LMO2 oncogene. The approach described here to identify enhancer-associated small insertion variants provides a foundation for further study of these abnormalities across human cancers.

Sequencing initiatives have detected multiple types of mutations in cancer. Here the authors, analysing enhancer-targeting sequence data, show that small insertions in transcriptional enhancers are frequently found near oncogenes, and demonstrate how one mutation deregulates expression of LMO2 in leukemia cells.

ASPirin Intervention for the REDuction of colorectal cancer risk (ASPIRED): a study protocol for a randomized controlled trial

BACKGROUND

Although aspirin is recommended for the prevention of colorectal cancer, the specific individuals for whom the benefits outweigh the risks are not clearly defined. Moreover, the precise mechanisms by which aspirin reduces the risk of cancer are unclear. We recently launched the ASPirin Intervention for the REDuction of colorectal cancer risk (ASPIRED) trial to address these uncertainties.

METHODS/DESIGN

ASPIRED is a prospective, double-blind, multidose, placebo-controlled, biomarker clinical trial of aspirin use in individuals previously diagnosed with colorectal adenoma. Individuals (n = 180) will be randomized in a 1:1:1 ratio to low-dose (81 mg/day) or standard-dose (325 mg/day) aspirin or placebo. At two study visits, participants will provide lifestyle, dietary and biometric data in addition to urine, saliva and blood specimens. Stool, grossly normal colorectal mucosal biopsies and cytology brushings will be collected during a flexible sigmoidoscopy without bowel preparation. The study will examine the effect of aspirin on urinary prostaglandin metabolites (PGE-M; primary endpoint), plasma inflammatory markers (macrophage inhibitory cytokine-1 (MIC-1)), colonic expression of transcription factor binding (transcription factor 7-like 2 (TCF7L2)), colonocyte gene expression, including hydroxyprostaglandin dehydrogenase 15-(NAD) (HPGD) and those that encode Wnt signaling proteins, colonic cellular nanocytology and oral and gut microbial composition and function.

DISCUSSION

Aspirin may prevent colorectal cancer through multiple, interrelated mechanisms. The ASPIRED trial will scrutinize these pathways and investigate putative mechanistically based risk-stratification biomarkers.

TRIAL REGISTRATION

This protocol is registered with the U.S. National Institutes of Health trial registry, ClinicalTrials.gov, under the identifier NCT02394769 . Registered on 16 March 2015.

The NOTCH1-MYC highway toward T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is a highly proliferative hematologic malignancy that results from the transformation of immature T-cell progenitors. Aberrant cell growth and proliferation in T-ALL lymphoblasts are sustained by activation of strong oncogenic drivers promoting cell anabolism and cell cycle progression. Oncogenic NOTCH signaling, which is activated in more than 65% of T-ALL patients by activating mutations in the NOTCH1 gene, has emerged as a major regulator of leukemia cell growth and metabolism. T-ALL NOTCH1 mutations result in ligand-independent and sustained NOTCH1-receptor signaling, which translates into activation of a broad transcriptional program dominated by upregulation of genes involved in anabolic pathways. Among these, the MYC oncogene plays a major role in NOTCH1-induced transformation. As result, the oncogenic activity of NOTCH1 in T-ALL is strictly dependent on MYC upregulation, which makes the NOTCH1-MYC regulatory circuit an attractive therapeutic target for the treatment of T-ALL.

The Shh Topological Domain Facilitates the Action of Remote Enhancers by Reducing the Effects of Genomic Distances

Gene expression often requires interaction between promoters and distant enhancers, which occur within the context of highly organized topologically associating domains (TADs). Using a series of engineered chromosomal rearrangements at the Shh locus, we carried out an extensive fine-scale characterization of the factors that govern the long-range regulatory interactions controlling Shh expression. We show that Shh enhancers act pervasively, yet not uniformly, throughout the TAD. Importantly, changing intra-TAD distances had no impact on Shh expression. In contrast, inversions disrupting the TAD altered global folding of the region and prevented regulatory contacts in a distance-dependent manner. Our data indicate that the Shh TAD promotes distance-independent contacts between distant regions that would otherwise interact only sporadically, enabling functional communication between them. In large genomes where genomic distances per se can limit regulatory interactions, this function of TADs could be as essential for gene expression as the formation of insulated neighborhoods.

MLL3/MLL4/COMPASS Family on Epigenetic Regulation of Enhancer Function and Cancer

During development, precise spatiotemporal patterns of gene expression are coordinately controlled by cis-regulatory modules known as enhancers. Their crucial role in development helped spur numerous studies aiming to elucidate the functional properties of enhancers within their physiological and disease contexts. In recent years, the role of enhancer malfunction in tissue-specific tumorigenesis is increasingly investigated. Here, we direct our focus to two primary players in enhancer regulation and their role in cancer pathogenesis: MLL3 and MLL4, members of the COMPASS family of histone H3 lysine 4 (H3K4) methyltransferases, and their complex-specific subunit UTX, a histone H3 lysine 27 (H3K27) demethylase. We review the most recent evidence on the underlying roles of MLL3/MLL4 and UTX in cancer and highlight key outstanding questions to help drive future research and contribute to our fundamental understanding of cancer and facilitate identification of therapeutic opportunities.

Genetic variation at the 8q24.21 renal cancer susceptibility locus affects HIF binding to a MYC enhancer

Clear cell renal cell carcinoma (ccRCC) is characterized by loss of function of the von Hippel-Lindau tumour suppressor (VHL) and unrestrained activation of hypoxia-inducible transcription factors (HIFs). Genetic and epigenetic determinants have an impact on HIF pathways. A recent genome-wide association study on renal cancer susceptibility identified single-nucleotide polymorphisms (SNPs) in an intergenic region located between the oncogenes MYC and PVT1. Here using assays of chromatin conformation, allele-specific chromatin immunoprecipitation and genome editing, we show that HIF binding to this regulatory element is necessary to trans-activate MYC and PVT1 expression specifically in cells of renal tubular origins. Moreover, we demonstrate that the risk-associated polymorphisms increase chromatin accessibility and activity as well as HIF binding to the enhancer. These findings provide further evidence that genetic variation at HIF-binding sites modulates the oncogenic transcriptional output of the VHL-HIF axis and provide a functional explanation for the disease-associated effects of SNPs in ccRCC.

3D genome organization in health and disease: emerging opportunities in cancer translational medicine

Organizing the DNA to fit inside a spatially constrained nucleus is a challenging problem that has attracted the attention of scientists across all disciplines of science. Increasing evidence has demonstrated the importance of genome geometry in several cellular contexts that affect human health. Among several approaches, the application of sequencing technologies has substantially increased our understanding of this intricate organization, also known as chromatin interactions. These structures are involved in transcriptional control of gene expression by connecting distal regulatory elements with their target genes and regulating co-transcriptional splicing. In addition, chromatin interactions play pivotal roles in the organization of the genome, the formation of structural variants, recombination, DNA replication and cell division. Mutations in factors that regulate chromatin interactions lead to the development of pathological conditions, for example, cancer. In this review, we discuss key findings that have shed light on the importance of these structures in the context of cancers, and highlight the applicability of chromatin interactions as potential biomarkers in molecular medicine as well as therapeutic implications of chromatin interactions.

Systematic mapping of functional enhancer-promoter connections with CRISPR interference

Gene expression in mammals is regulated by noncoding elements that can affect physiology and disease, yet the functions and target genes of most noncoding elements remain unknown. We present a high-throughput approach that uses clustered regularly interspaced short palindromic repeats (CRISPR) interference (CRISPRi) to discover regulatory elements and identify their target genes. We assess >1 megabase of sequence in the vicinity of two essential transcription factors, MYC and GATA1, and identify nine distal enhancers that control gene expression and cellular proliferation. Quantitative features of chromatin state and chromosome conformation distinguish the seven enhancers that regulate MYC from other elements that do not, suggesting a strategy for predicting enhancer-promoter connectivity. This CRISPRi-based approach can be applied to dissect transcriptional networks and interpret the contributions of noncoding genetic variation to human disease.

BNC2 is a putative tumor suppressor gene in high-grade serous ovarian carcinoma and impacts cell survival after oxidative stress

Rs3814113 is the single-nucleotide polymorphism (SNP) showing the strongest association with high-grade serous ovarian carcinoma (HGSOC) incidence and is located in an intergenic region about 44 kb downstream of basonuclin 2 (BNC2) gene. Lifetime number of ovulations is associated with increased risk to develop HGSOC, probably because of cell damage of extrauterine Müllerian epithelium by ovulation-induced oxidative stress. However, the impact of low-penetrance HGSOC risk alleles (e.g. rs3814113) on the damage induced by oxidative stress remains unclear. Therefore, the purpose of this study was to investigate whether rs3814113 genetic interval regulates BNC2 expression and whether BNC2 expression levels impact on cell survival after oxidative stress. To do this, we analyzed gene expression levels of BNC2 first in HGSOC data sets and then in an isogenic cell line that we engineered to carry a 5 kb deletion around rs3814113. Finally, we silenced BNC2 and measured surviving cells after hydrogen peroxide (H2O2) treatment to simulate oxidative stress after ovulation. In this paper, we describe that BNC2 expression levels are reduced in HGSOC samples compared with control samples, and that BNC2 expression levels decrease following oxidative stress and ovulation in vitro and in vivo, respectively. Moreover, deletion of 5 kb surrounding rs3814113 decreases BNC2 expression levels in an isogenic cell line, and silencing of BNC2 expression levels increases cell survival after H2O2 treatment. Altogether, our findings suggest that the intergenic region located around rs3814113 regulates BNC2 expression, which in turn affects cell survival after oxidative stress response. Indeed, HGSOC samples present lower BNC2 expression levels that probably, in the initial phases of oncogenic transformation, conferred resistance to oxidative stress and ultimately reduced the clearance of cells with oxidative-induced damages.

Enhancer deregulation in cancer and other diseases

Mutations in enhancer-associated chromatin-modifying components and genomic alterations in non-coding regions of the genome occur frequently in cancer, and other diseases pointing to the importance of enhancer fidelity to ensure proper tissue homeostasis. In this review, I will use specific examples to discuss how mutations in chromatin-modifying factors might affect enhancer activity of disease-relevant genes. I will then consider direct evidence from single nucleotide polymorphisms, small insertions, or deletions but also larger genomic rearrangements such as duplications, deletions, translocations, and inversions of specific enhancers to demonstrate how they have the ability to impact enhancer activity of disease genes including oncogenes and tumor suppressor genes. Considering that the scientific community only fairly recently has begun to focus its attention on "enhancer malfunction" in disease, I propose that multiple new enhancer-regulated and disease-relevant processes will be uncovered in the near future that will constitute the mechanistic basis for novel therapeutic avenues.

The rs6983267 SNP and long non-coding RNA CARLo-5 are associated with endometrial carcinoma

The single nucleotide polymorphism (SNP) rs6983267 and cancer-associated region long non-coding RNA (CARLo-5) are associated with various human cancers. This study aimed to investigate the expression of CARLo-5 in endometrial carcinoma (EC) and its relationship with clinicopathological features and patient survival. The association of the rs6983267 SNP with EC risk and its involvement in the regulation of CARLo-5 expression in EC were investigated. The rs6983267 SNP was genotyped by polymerase chain reaction (PCR) and ligase detection reaction in 543 EC patients and 584 controls. The expression of CARLo-5 in 108 EC tissues and 66 normal endometrial tissues (NETs) was determined using quantitative real-time PCR. The genotype and allele distributions of the rs6983267 SNP differed significantly between patients and controls. There was a significant correlation between the rs6983267 genotypes and lymph node metastasis of EC patients (P = 0.026). CARLo-5 expression was significantly higher in EC tissues than in NETs (P < 0.001) and significantly associated with FIGO stage (P = 0.029) and lymph node metastasis (P = 0.030). Patients with high CARLo-5 expression had significantly shorter overall survival than those with low CARLo-5 expression (P = 0.003). The rs6983267 genotype was significantly correlated with CARLo-5 expression (P < 0.05). In conclusion, CARLo-5 was identified as a pro-oncogenic lncRNA that may play an important role in EC progression and represent a prognostic marker for EC. The expression of CARLo-5 was significantly correlated with the rs6983267 genotype associated with increased susceptibility to EC. Environ. Mol. Mutagen. 57:508-515, 2016. © 2016 Wiley Periodicals, Inc.

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.

Fine-Mapping of Common Genetic Variants Associated with Colorectal Tumor Risk Identified Potential Functional Variants

Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) associated with colorectal cancer risk. These SNPs may tag correlated variants with biological importance. Fine-mapping around GWAS loci can facilitate detection of functional candidates and additional independent risk variants. We analyzed 11,900 cases and 14,311 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. To fine-map genomic regions containing all known common risk variants, we imputed high-density genetic data from the 1000 Genomes Project. We tested single-variant associations with colorectal tumor risk for all variants spanning genomic regions 250-kb upstream or downstream of 31 GWAS-identified SNPs (index SNPs). We queried the University of California, Santa Cruz Genome Browser to examine evidence for biological function. Index SNPs did not show the strongest association signals with colorectal tumor risk in their respective genomic regions. Bioinformatics analysis of SNPs showing smaller P-values in each region revealed 21 functional candidates in 12 loci (5q31.1, 8q24, 11q13.4, 11q23, 12p13.32, 12q24.21, 14q22.2, 15q13, 18q21, 19q13.1, 20p12.3, and 20q13.33). We did not observe evidence of additional independent association signals in GWAS-identified regions. Our results support the utility of integrating data from comprehensive fine-mapping with expanding publicly available genomic databases to help clarify GWAS associations and identify functional candidates that warrant more onerous laboratory follow-up. Such efforts may aid the eventual discovery of disease-causing variant(s).

TCGA Workflow: Analyze cancer genomics and epigenomics data using Bioconductor packages

Biotechnological advances in sequencing have led to an explosion of publicly available data via large international consortia such as The Cancer Genome Atlas (TCGA), The Encyclopedia of DNA Elements (ENCODE), and The NIH Roadmap Epigenomics Mapping Consortium (Roadmap). These projects have provided unprecedented opportunities to interrogate the epigenome of cultured cancer cell lines as well as normal and tumor tissues with high genomic resolution. The Bioconductor project offers more than 1,000 open-source software and statistical packages to analyze high-throughput genomic data. However, most packages are designed for specific data types (e.g. expression, epigenetics, genomics) and there is no one comprehensive tool that provides a complete integrative analysis of the resources and data provided by all three public projects. A need to create an integration of these different analyses was recently proposed. In this workflow, we provide a series of biologically focused integrative analyses of different molecular data. We describe how to download, process and prepare TCGA data and by harnessing several key Bioconductor packages, we describe how to extract biologically meaningful genomic and epigenomic data. Using Roadmap and ENCODE data, we provide a work plan to identify biologically relevant functional epigenomic elements associated with cancer. To illustrate our workflow, we analyzed two types of brain tumors: low-grade glioma (LGG) versus high-grade glioma (glioblastoma multiform or GBM). This workflow introduces the following Bioconductor packages: AnnotationHub, ChIPSeeker, ComplexHeatmap, pathview, ELMER, GAIA, MINET, RTCGAToolbox,  TCGAbiolinks.

TCGA Workflow: Analyze cancer genomics and epigenomics data using Bioconductor packages

Biotechnological advances in sequencing have led to an explosion of publicly available data via large international consortia such as The Cancer Genome Atlas (TCGA), The Encyclopedia of DNA Elements (ENCODE), and The NIH Roadmap Epigenomics Mapping Consortium (Roadmap). These projects have provided unprecedented opportunities to interrogate the epigenome of cultured cancer cell lines as well as normal and tumor tissues with high genomic resolution. The Bioconductor project offers more than 1,000 open-source software and statistical packages to analyze high-throughput genomic data. However, most packages are designed for specific data types (e.g. expression, epigenetics, genomics) and there is no one comprehensive tool that provides a complete integrative analysis of the resources and data provided by all three public projects. A need to create an integration of these different analyses was recently proposed. In this workflow, we provide a series of biologically focused integrative analyses of different molecular data. We describe how to download, process and prepare TCGA data and by harnessing several key Bioconductor packages, we describe how to extract biologically meaningful genomic and epigenomic data. Using Roadmap and ENCODE data, we provide a work plan to identify biologically relevant functional epigenomic elements associated with cancer. To illustrate our workflow, we analyzed two types of brain tumors: low-grade glioma (LGG) versus high-grade glioma (glioblastoma multiform or GBM). This workflow introduces the following Bioconductor packages: AnnotationHub, ChIPSeeker, ComplexHeatmap, pathview, ELMER, GAIA, MINET, RTCGAToolbox,  TCGAbiolinks.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

IDH-mutant glioma specific association of rs55705857 located at 8q24.21 involves MYC deregulation

The single nucleotide polymorphism rs55705857, located in a non-coding but evolutionarily conserved region at 8q24.21, is strongly associated with IDH-mutant glioma development and was suggested to be a causal variant. However, the molecular mechanism underlying this association has remained unknown. With a case control study in 285 gliomas, 316 healthy controls, 380 systemic cancers, 31 other CNS-tumors, and 120 IDH-mutant cartilaginous tumors, we identified that the association was specific to IDH-mutant gliomas. Odds-ratios were 9.25 (5.17–16.52; 95% CI) for IDH-mutated gliomas and 12.85 (5.94–27.83; 95% CI) for IDH-mutated, 1p/19q co-deleted gliomas. Decreasing strength with increasing anaplasia implied a modulatory effect. No somatic mutations were noted at this locus in 114 blood-tumor pairs, nor was there a copy number difference between risk-allele and only-ancestral allele carriers. CCDC26 RNA-expression was rare and not different between the two groups. There were only minor subtype-specific differences in common glioma driver genes. RNA sequencing and LC-MS/MS comparisons pointed to significantly altered MYC-signaling. Baseline enhancer activity of the conserved region specifically on the MYC promoter and its further positive modulation by the SNP risk-allele was shown in vitro. Our findings implicate MYC deregulation as the underlying cause of the observed association.

A Common Variant at the 14q32 Endometrial Cancer Risk Locus Activates AKT1 through YY1 Binding

A recent meta-analysis of multiple genome-wide association and follow-up endometrial cancer case-control datasets identified a novel genetic risk locus for this disease at chromosome 14q32.33. To prioritize the functional SNP(s) and target gene(s) at this locus, we employed an in silico fine-mapping approach using genotyped and imputed SNP data for 6,608 endometrial cancer cases and 37,925 controls of European ancestry. Association and functional analyses provide evidence that the best candidate causal SNP is rs2494737. Multiple experimental analyses show that SNP rs2494737 maps to a silencer element located within AKT1, a member of the PI3K/AKT/MTOR intracellular signaling pathway activated in endometrial tumors. The rs2494737 risk A allele creates a YY1 transcription factor-binding site and abrogates the silencer activity in luciferase assays, an effect mimicked by transfection of YY1 siRNA. Our findings suggest YY1 is a positive regulator of AKT1, mediating the stimulatory effects of rs2494737 increasing endometrial cancer risk. Identification of an endometrial cancer risk allele within a member of the PI3K/AKT signaling pathway, more commonly activated in tumors by somatic alterations, raises the possibility that well tolerated inhibitors targeting this pathway could be candidates for evaluation as chemopreventive agents in individuals at high risk of developing endometrial cancer.

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.

The importance of p53 pathway genetics in inherited and somatic cancer genomes

Decades of research have shown that mutations in the p53 stress response pathway affect the incidence of diverse cancers more than mutations in other pathways. However, most evidence is limited to somatic mutations and rare inherited mutations. Using newly abundant genomic data, we demonstrate that commonly inherited genetic variants in the p53 pathway also affect the incidence of a broad range of cancers more than variants in other pathways. The cancer-associated single nucleotide polymorphisms (SNPs) of the p53 pathway have strikingly similar genetic characteristics to well-studied p53 pathway cancer-causing somatic mutations. Our results enable insights into p53-mediated tumour suppression in humans and into p53 pathway-based cancer surveillance and treatment strategies.

Models of human core transcriptional regulatory circuitries

A small set of core transcription factors (TFs) dominates control of the gene expression program in embryonic stem cells and other well-studied cellular models. These core TFs collectively regulate their own gene expression, thus forming an interconnected auto-regulatory loop that can be considered the core transcriptional regulatory circuitry (CRC) for that cell type. There is limited knowledge of core TFs, and thus models of core regulatory circuitry, for most cell types. We recently discovered that genes encoding known core TFs forming CRCs are driven by super-enhancers, which provides an opportunity to systematically predict CRCs in poorly studied cell types through super-enhancer mapping. Here, we use super-enhancer maps to generate CRC models for 75 human cell and tissue types. These core circuitry models should prove valuable for further investigating cell-type–specific transcriptional regulation in healthy and diseased cells.

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

Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with increased risk for colorectal cancer (CRC). A molecular understanding of the functional consequences of this genetic variation is complicated because most GWAS SNPs are located in non-coding regions. We used epigenomic information to identify H3K27Ac peaks in HCT116 colon cancer cells that harbor SNPs associated with an increased risk for CRC. Employing CRISPR/Cas9 nucleases, we deleted a CRC risk-associated H3K27Ac peak from HCT116 cells and observed large-scale changes in gene expression, resulting in decreased expression of many nearby genes. As a comparison, we showed that deletion of a robust H3K27Ac peak not associated with CRC had minimal effects on the transcriptome. Interestingly, although there is no H3K27Ac peak in HEK293 cells in the E7 region, deletion of this region in HEK293 cells decreased expression of several of the same genes that were downregulated in HCT116 cells, including the MYC oncogene. Accordingly, deletion of E7 causes changes in cell culture assays in HCT116 and HEK293 cells. In summary, we show that effects on the transcriptome upon deletion of a distal regulatory element cannot be predicted by the size or presence of an H3K27Ac peak.

Making sense of GWAS: using epigenomics and genome engineering to understand the functional relevance of SNPs in non-coding regions of the human genome

Considerable progress towards an understanding of complex diseases has been made in recent years due to the development of high-throughput genotyping technologies. Using microarrays that contain millions of single-nucleotide polymorphisms (SNPs), Genome Wide Association Studies (GWASs) have identified SNPs that are associated with many complex diseases or traits. For example, as of February 2015, 2111 association studies have identified 15,396 SNPs for various diseases and traits, with the number of identified SNP-disease/trait associations increasing rapidly in recent years. However, it has been difficult for researchers to understand disease risk from GWAS results. This is because most GWAS-identified SNPs are located in non-coding regions of the genome. It is important to consider that the GWAS-identified SNPs serve only as representatives for all SNPs in the same haplotype block, and it is equally likely that other SNPs in high linkage disequilibrium (LD) with the array-identified SNPs are causal for the disease. Because it was hoped that disease-associated coding variants would be identified if the true casual SNPs were known, investigators have expanded their analyses using LD calculation and fine-mapping. However, such analyses also identified risk-associated SNPs located in non-coding regions. Thus, the GWAS field has been left with the conundrum as to how a single-nucleotide change in a non-coding region could confer increased risk for a specific disease. One possible answer to this puzzle is that the variant SNPs cause changes in gene expression levels rather than causing changes in protein function. This review provides a description of (1) advances in genomic and epigenomic approaches that incorporate functional annotation of regulatory elements to prioritize the disease risk-associated SNPs that are located in non-coding regions of the genome for follow-up studies, (2) various computational tools that aid in identifying gene expression changes caused by the non-coding disease-associated SNPs, and (3) experimental approaches to identify target genes of, and study the biological phenotypes conferred by, non-coding disease-associated SNPs.

A dynamic exchange of TCF3 and TCF4 transcription factors controls MYC expression in colorectal cancer cells

Deregulated Wnt/β-catenin signaling promotes colorectal cancer (CRC) by activating expression of the c-MYC proto-oncogene (MYC). In the nucleus, the β-catenin transcriptional co-activator binds T-cell factor (TCF) transcription factors, and together TCF/β-catenin complexes activate MYC expression through Wnt responsive DNA regulatory elements (WREs). The MYC 3' WRE maps 1.4-kb downstream from the MYC transcription stop site and binds TCF4/β-catenin transcription complexes to activate MYC. However, the underlying mechanisms for how this element operates are not fully understood. Here, we report that the TCF family member, TCF3, plays an important role in regulating MYC expression in CRCs. We demonstrate that TCF3 binds the MYC 3' WRE to repress MYC. When TCF3 is depleted using shRNAs, the MYC 3' WRE is more available to bind TCF4/β-catenin complexes. Stimulating downstream Wnt/β-catenin signaling by inhibiting GSK3β causes an exchange of TCF3 with TCF4/β-catenin complexes to activate MYC. Finally, this transcription factor switch at the MYC 3' WRE controls MYC expression as quiescent cells re-enter the cell cycle and progress to S phase. These results indicate that a dynamic interplay of TCF transcription factors governs MYC gene expression in CRCs.