Hotspots of aberrant enhancer activity punctuate the colorectal cancer epigenome.

In addition to mutations in genes, aberrant enhancer element activity at non-coding regions of the genome is a key driver of tumorigenesis. Here, we perform epigenomic enhancer profiling of a cohort of more than forty genetically diverse human colorectal cancer (CRC) specimens. Using normal colonic crypt epithelium as a comparator, we identify enhancers with recurrently gained or lost activity across CRC specimens. Of the enhancers highly recurrently activated in CRC, most are constituents of super enhancers, are occupied by AP-1 and cohesin complex members, and originate from primed chromatin. Many activate known oncogenes, and CRC growth can be mitigated through pharmacologic inhibition or genome editing of these loci. Nearly half of all GWAS CRC risk loci co-localize to recurrently activated enhancers. These findings indicate that the CRC epigenome is defined by highly recurrent epigenetic alterations at enhancers which activate a common, aberrant transcriptional programme critical for CRC growth and survival.

Transposon mutagenesis identifies candidate genes that cooperate with loss of transforming growth factor-beta signaling in mouse intestinal neoplasms.

Colorectal cancer (CRC) results from the accumulation of gene mutations and epigenetic alterations in colon epithelial cells, which promotes CRC formation through deregulating signaling pathways. One of the most commonly deregulated signaling pathways in CRC is the transforming growth factor ß (TGF-ß) pathway. Importantly, the effects of TGF-ß signaling inactivation in CRC are modified by concurrent mutations in the tumor cell, and these concurrent mutations determine the ultimate biological effects of impaired TGF-ß signaling in the tumor. However, many of the mutations that cooperate with the deregulated TGF-ß signaling pathway in CRC remain unknown. Therefore, we sought to identify candidate driver genes that promote the formation of CRC in the setting of TGF-ß signaling inactivation. We performed a forward genetic screen in mice carrying conditionally inactivated alleles of the TGF-ß receptor, type II (Tgfbr2) using Sleeping Beauty (SB) transposon mediated mutagenesis. We used TAPDANCE and Gene-centric statistical methods to identify common insertion sites (CIS) and, thus, candidate tumor suppressor genes and oncogenes within the tumor genome. CIS analysis of multiple neoplasms from these mice identified many candidate Tgfbr2 cooperating genes and the Wnt/ß-catenin, Hippo and MAPK pathways as the most commonly affected pathways. Importantly, the majority of candidate genes were also found to be mutated in human CRC. The SB transposon system provides an unbiased method to identify Tgfbr2 cooperating genes in mouse CRC that are functionally relevant and that may provide further insight into the pathogenesis of human CRC.

Induction of KIAA1199/CEMIP is associated with colon cancer phenotype and poor patient survival.

Genes induced in colon cancer provide novel candidate biomarkers of tumor phenotype and aggressiveness. We originally identified KIAA1199 (now officially called CEMIP) as a transcript highly induced in colon cancer: initially designating the transcript as Colon Cancer Secreted Protein 1. We molecularly characterized CEMIP expression both at the mRNA and protein level and found it is a secreted protein induced an average of 54-fold in colon cancer. Knockout of CEMIPreduced the ability of human colon cancer cells to form xenograft tumors in athymic mice. Tumors that did grow had increased deposition of hyaluronan, linking CEMIP participation in hyaluronan degradation to the modulation of tumor phenotype. We find CEMIP mRNA overexpression correlates with poorer patient survival. In stage III only (n = 31) or in combined stage II plus stage III colon cancer cases (n = 73), 5-year overall survival was significantly better (p = 0.004 and p = 0.0003, respectively) among patients with low CEMIP expressing tumors than those with high CEMIP expressing tumors. These results demonstrate that CEMIP directly facilitates colon tumor growth, and high CEMIP expression correlates with poor outcome in stage III and in stages II+III combined cohorts. We present CEMIP as a candidate prognostic marker for colon cancer and a potential therapeutic target.

Combinatorial effects of multiple enhancer variants in linkage disequilibrium dictate levels of gene expression to confer susceptibility to common traits.

DNA variants (SNPs) that predispose to common traits often localize within noncoding regulatory elements such as enhancers. Moreover, loci identified by genome-wide association studies (GWAS) often contain multiple SNPs in linkage disequilibrium (LD), any of which may be causal. Thus, determining the effect of these multiple variant SNPs on target transcript levels has been a major challenge. Here, we provide evidence that for six common autoimmune disorders (rheumatoid arthritis, Crohn's disease, celiac disease, multiple sclerosis, lupus, and ulcerative colitis), the GWAS association arises from multiple polymorphisms in LD that map to clusters of enhancer elements active in the same cell type. This finding suggests a "multiple enhancer variant" hypothesis for common traits, where several variants in LD impact multiple enhancers and cooperatively affect gene expression. Using a novel method to delineate enhancer-gene interactions, we show that multiple enhancer variants within a given locus typically target the same gene. Using available data from HapMap and B lymphoblasts as a model system, we provide evidence at numerous loci that multiple enhancer variants cooperatively contribute to altered expression of their gene targets. The effects on target transcript levels tend to be modest and can be either gain- or loss-of-function. Additionally, the genes associated with multiple enhancer variants encode proteins that are often functionally related and enriched in common pathways. Overall, the multiple enhancer variant hypothesis offers a new paradigm by which noncoding variants can confer susceptibility to common traits.

Two Distinct Categories of Focal Deletions in Cancer Genomes.

One of the key questions about genomic alterations in cancer is whether they are functional in the sense of contributing to the selective advantage of tumor cells. The frequency with which an alteration occurs might reflect its ability to increase cancer cell growth, or alternatively, enhanced instability of a locus may increase the frequency with which it is found to be aberrant in tumors, regardless of oncogenic impact. Here we've addressed this on a genome-wide scale for cancer-associated focal deletions, which are known to pinpoint both tumor suppressor genes (tumor suppressors) and unstable loci. Based on DNA copy number analysis of over one-thousand human cancers representing ten different tumor types, we observed five loci with focal deletion frequencies above 5%, including the A2BP1 gene at 16p13.3 and the MACROD2 gene at 20p12.1. However, neither RNA expression nor functional studies support a tumor suppressor role for either gene. Further analyses suggest instead that these are sites of increased genomic instability and that they resemble common fragile sites (CFS). Genome-wide analysis revealed properties of CFS-like recurrent deletions that distinguish them from deletions affecting tumor suppressor genes, including their isolation at specific loci away from other genomic deletion sites, a considerably smaller deletion size, and dispersal throughout the affected locus rather than assembly at a common site of overlap. Additionally, CFS-like deletions have less impact on gene expression and are enriched in cell lines compared to primary tumors. We show that loci affected by CFS-like deletions are often distinct from known common fragile sites. Indeed, we find that each tumor tissue type has its own spectrum of CFS-like deletions, and that colon cancers have many more CFS-like deletions than other tumor types. We present simple rules that can pinpoint focal deletions that are not CFS-like and more likely to affect functional tumor suppressors.

H3K4me3 inversely correlates with DNA methylation at a large class of non-CpG-island-containing start sites.

BACKGROUND: In addition to mutations, epigenetic silencing of genes has been recognized as a fundamental mechanism that promotes human carcinogenesis. To date, characterization of epigenetic gene silencing has largely focused on genes in which silencing is mediated by hypermethylation of promoter-associated CpG islands, associated with loss of the H3K4me3 chromatin mark. Far less is known about promoters lacking CpG-islands or genes that are repressed by alternative mechanisms. METHODS: We performed integrative ChIP-chip, DNase-seq, and global gene expression analyses in colon cancer cells and normal colon mucosa to characterize chromatin features of both CpG-rich and CpG-poor promoters of genes that undergo silencing in colon cancer. RESULTS: Epigenetically repressed genes in colon cancer separate into two classes based on retention or loss of H3K4me3 at transcription start sites. Quantitatively, of transcriptionally repressed genes that lose H3K4me3 in colon cancer (K4-dependent genes), a large fraction actually lacks CpG islands. Nonetheless, similar to CpG-island containing genes, cytosines located near the start sites of K4-dependent genes become DNA hypermethylated, and repressed K4-dependent genes can be reactivated with 5-azacytidine. Moreover, we also show that when the H3K4me3 mark is retained, silencing of CpG island-associated genes can proceed through an alternative mechanism in which repressive chromatin marks are recruited. CONCLUSIONS: H3K4me3 equally protects from DNA methylation at both CpG-island and non-CpG island start sites in colon cancer. Moreover, the results suggest that CpG-rich genes repressed by loss of H3K4me3 and DNA methylation represent special instances of a more general epigenetic mechanism of gene silencing, one in which gene silencing is mediated by loss of H3K4me3 and methylation of non-CpG island promoter-associated cytosines.

Epigenomic enhancer profiling defines a signature of colon cancer.

Cancer is characterized by gene expression aberrations. Studies have largely focused on coding sequences and promoters, even though distal regulatory elements play a central role in controlling transcription patterns. We used the histone mark H3K4me1 to analyze gain and loss of enhancer activity genome-wide in primary colon cancer lines relative to normal colon crypts. We identified thousands of variant enhancer loci (VELs) that comprise a signature that is robustly predictive of the in vivo colon cancer transcriptome. Furthermore, VELs are enriched in haplotype blocks containing colon cancer genetic risk variants, implicating these genomic regions in colon cancer pathogenesis. We propose that reproducible changes in the epigenome at enhancer elements drive a specific transcriptional program to promote colon carcinogenesis.

BCL9 promotes tumor progression by conferring enhanced proliferative, metastatic, and angiogenic properties to cancer cells.

Several components of the Wnt signaling cascade have been shown to function either as tumor suppressor proteins or as oncogenes in multiple human cancers, underscoring the relevance of this pathway in oncogenesis and the need for further investigation of Wnt signaling components as potential targets for cancer therapy. Here, using expression profiling analysis as well as in vitro and in vivo functional studies, we show that the Wnt pathway component BCL9 is a novel oncogene that is aberrantly expressed in human multiple myeloma as well as colon carcinoma. We show that BCL9 enhances beta-catenin-mediated transcriptional activity regardless of the mutational status of the Wnt signaling components and increases cell proliferation, migration, invasion, and the metastatic potential of tumor cells by promoting loss of epithelial and gain of mesenchymal-like phenotype. Most importantly, BCL9 knockdown significantly increased the survival of xenograft mouse models of cancer by reducing tumor load, metastasis, and host angiogenesis through down-regulation of c-Myc, cyclin D1, CD44, and vascular endothelial growth factor expression by tumor cells. Together, these findings suggest that deregulation of BCL9 is an important contributing factor to tumor progression. The pleiotropic roles of BCL9 reported in this study underscore its value as a drug target for therapeutic intervention in several malignancies associated with aberrant Wnt signaling.

Discovery and scoring of protein interaction subnetworks discriminative of late stage human colon cancer.

We used a systems biology approach to identify and score protein interaction subnetworks whose activity patterns are discriminative of late stage human colorectal cancer (CRC) versus control in colonic tissue. We conducted two gel-based proteomics experiments to identify significantly changing proteins between normal and late stage tumor tissues obtained from an adequately sized cohort of human patients. A total of 67 proteins identified by these experiments was used to seed a search for protein-protein interaction subnetworks. A scoring scheme based on mutual information, calculated using gene expression data as a proxy for subnetwork activity, was developed to score the targets in the subnetworks. Based on this scoring, the subnetwork was pruned to identify the specific protein combinations that were significantly discriminative of late stage cancer versus control. These combinations could not be discovered using only proteomics data or by merely clustering the gene expression data. We then analyzed the resultant pruned subnetwork for biological relevance to human CRC. A number of the proteins in these smaller subnetworks have been associated with the progression (CSNK2A2, PLK1, and IGFBP3) or metastatic potential (PDGFRB) of CRC. Others have been recently identified as potential markers of CRC (IFITM1), and the role of others is largely unknown in this disease (CCT3, CCT5, CCT7, and GNA12). The functional interactions represented by these signatures provide new experimental hypotheses that merit follow-on validation for biological significance in this disease. Overall the method outlines a quantitative approach for integrating proteomics data, gene expression data, and the wealth of accumulated legacy experimental data to discover significant protein subnetworks specific to disease.

Mutational inactivation of TGFBR2 in microsatellite unstable colon cancer arises from the cooperation of genomic instability and the clonal outgrowth of transforming growth factor beta resistant cells.

The mutational inactivation of transforming growth factor beta receptor type II (TGFBR2) occurs in approximately 30% of colon cancers and promotes the formation of colon cancer by inhibiting the tumor suppressor activity of the TGFB signaling pathway. TGFBR2 mutations occur in >90% of microsatellite unstable (MSI) colon cancers and affect a polyadenine tract in exon 3 of TGFBR2, called BAT-RII, which is vulnerable to mutation in the setting of DNA mismatch repair (MMR) system deficiency. In light of the vulnerable nature of the BAT-RII tract in the setting of MMR inactivation and the favorable effects of TGFBR2 inactivation in colon cancer, analysis of TGFBR2 inactivation provides an opportunity to assess the roles of genomic instability vs. clonal selection in cells acquiring TGFBR2 BAT-RII tract mutations in MSI colon cancer formation. The contribution of genomic instability and/or clonal evolution to the mutational inactivation of TGBFR2 in MSI colon cancers has not been studied in a systematic way that would allow a determination of the relative contribution of these two mechanisms in the formation of MSI colon cancer. It has not been demonstrated whether the BAT-RII tract mutations are strictly a consequence of the BAT-RII region being hypermutable in the setting of MMR deficiency or if the mutations are rather a consequence of clonal selection pressure against the TGFB receptor. Through the use of defined cell line systems, we show that both genomic instability and clonal selection of TGFB resistant cells contribute to the high frequency of TGFBR2 mutations in MSI colon cancer.

Detection in fecal DNA of colon cancer-specific methylation of the nonexpressed vimentin gene.

BACKGROUND: Increased DNA methylation is an epigenetic alteration that is common in human cancers and is often associated with transcriptional silencing. Aberrantly methylated DNA has also been proposed as a potential tumor marker. However, genes such as vimentin, which are transcriptionally silent in normal epithelium, have not until now been considered as targets for cancer-associated aberrant methylation and for use as cancer markers. METHODS: We applied methylation-specific polymerase chain reaction to the vimentin gene, which is transcriptionally silent in normal colonocytes, and compared methylation of vimentin exon 1 in cancer tissues and in fecal DNA from colon cancer patients versus control samples from healthy subjects. RESULTS: Vimentin exon-1 sequences were unmethylated in 45 of 46 normal colon tissues. In contrast, vimentin exon-1 sequences were methylated in 83% (38 of 46) and 53% (57 of 107) of tumors from two independently collected groups of colon cancer patients. When evaluated as a marker for colon cancer detection in fecal DNA from another set of colon cancer patients, aberrant vimentin methylation was detected in fecal DNA from 43 of 94 patients, for a sensitivity of 46% (95% confidence interval [CI] = 35% to 56%). The sensitivity for detecting stage I and II cancers was 43% (26 of 60 case patients) (95% CI = 31% to 57%). Only 10% (20 of 198 case patients) of control fecal DNA samples from cancer-free individuals tested positive for vimentin methylation, for a specificity of 90% (95% CI = 85% to 94%). CONCLUSIONS: Aberrant methylation of exon-1 sequences within the nontranscribed vimentin gene is a novel molecular biomarker of colon cancer and can be successfully detected in fecal DNA to identify nearly half of individuals with colon cancer.

SLC5A8, a sodium transporter, is a tumor suppressor gene silenced by methylation in human colon aberrant crypt foci and cancers.

We identify a gene, SLC5A8, and show it is a candidate tumor suppressor gene whose silencing by aberrant methylation is a common and early event in human colon neoplasia. Aberrant DNA methylation has been implicated as a component of an epigenetic mechanism that silences genes in human cancers. Using restriction landmark genome scanning, we performed a global search to identify genes that would be aberrantly methylated at high frequency in human colon cancer. From among 1,231 genomic NotI sites assayed, site 3D41 was identified as methylated in 11 of 12 colon cancers profiled. Site 3D41 mapped to exon 1 of SLC5A8, a transcript that we assembled. In normal colon mucosa we found that SLC5A8 exon 1 is unmethylated and SLC5A8 transcript is expressed. In contrast, SLC5A8 exon 1 proved to be aberrantly methylated in 59% of primary colon cancers and 52% of colon cancer cell lines. SLC5A8 exon 1 methylated cells were uniformly silenced for SLC5A8 expression, but reactivated expression on treatment with a demethylating drug, 5-azacytidine. Transfection of SLC5A8 suppressed colony growth in each of three SLC5A8-deficient cell lines, but showed no suppressive effect in any of three SLC5A8-proficient cell lines. SLC5A8 exon 1 methylation is an early event, detectable in colon adenomas, and in even earlier microscopic colonic aberrant crypt foci. Structural homology and functional testing demonstrated that SLC5A8 is a member of the family of sodium solute symporters, which are now added as a class of candidate colon cancer suppressor genes.

HLTF gene silencing in human colon cancer.

Chromatin remodeling enzymes are increasingly implicated in a variety of important cellular functions. Various components of chromatin remodeling complexes, including several members of the SWI/SNF family, have been shown to be disrupted in cancer. In this study we identified as a target for gene inactivation in colon cancer the gene for helicase-like transcription factor (HLTF), a SWI/SNF family protein. Loss of HLTF expression accompanied by HLTF promoter methylation was noted in nine of 34 colon cancer cell lines. In these cell lines HLTF expression was restored by treatment with the demethylating agent 5-azacytidine. In further studies of primary colon cancer tissues, HLTF methylation was detected in 27 of 63 cases (43%). No methylation of HLTF was detected in breast or lung cancers, suggesting selection for HLTF methylation in colonic malignancies. Transfection of HLTF suppressed 75% of colony growth in each of three different HLTF-deficient cell lines, but showed no suppressive effect in any of three HLTF-proficient cell lines. These findings show that HLTF is a common target for methylation and epigenetic gene silencing in colon cancer and suggest HLTF is a candidate colon cancer suppressor gene.

Chromosomal autonomy of hMLH1 methylation in colon cancer.

Silencing of hMLH1 expression by aberrant hMLH1 promoter methylation accounts for the majority of sporadic colon cancers with microsatellite instability. We have previously shown hMLH1 silencing is biallelic and actively maintained. To study the mechanism of aberrant hMLH1 methylation, we assayed whether an hMLH1 methylated cell could transfer methylation and silencing to an exogenous hMLH1 promoter in somatic cell hybrids between hMLH1 methylated-silenced and hMLH1 unmethylated-expressing colon cancer cells. Conversely, we assayed whether these hybrids could reactivate expression of initially methylated and silenced hMLH1 alleles. Compellingly, within the hybrids each hMLH1 allele remained unchanged, retaining the expression status of its parental cell of origin. This chromosomal autonomy may not be simply determined by DNA methylation, as it is reasserted after experimentally forced demethylation of all hMLH1 alleles in the hybrids. Confirming findings included hMLH1 methylated cells being unable to methylate single transferred exogenous hMLH1 expressing chromosomes or transfected hMLH1 reporter constructs. hMLH1 silencing does not conform to either a dominant or recessive model, and is not determined by trans-acting factors differing between hMLH1 expressing or silenced genomes. We posit that hMLH1 methylation is dependent on and maintained by cis chromosomal marks, whose nature remains to be elucidated.