Hotspots of aberrant enhancer activity punctuate the colorectal cancer epigenome

DKK1-SE recruits AP1 to activate the target gene DKK1 thereby promoting pancreatic cancer progression

Super-enhancers are a class of DNA cis-regulatory elements that can regulate cell identity, cell fate, stem cell pluripotency, and even tumorigenesis. Increasing evidence shows that epigenetic modifications play an important role in the pathogenesis of various types of cancer. However, the current research is far from enough to reveal the complex mechanism behind it. This study found a super-enhancer enriched with abnormally active histone modifications in pancreatic ductal adenocarcinoma (PDAC), called DKK1-super-enhancer (DKK1-SE). The major active component of DKK1-SE is component enhancer e1. Mechanistically, AP1 induces chromatin remodeling in component enhancer e1 and activates the transcriptional activity of DKK1. Moreover, DKK1 was closely related to the malignant clinical features of PDAC. Deletion or knockdown of DKK1-SE significantly inhibited the proliferation, colony formation, motility, migration, and invasion of PDAC cells in vitro, and these phenomena were partly mitigated upon rescuing DKK1 expression. In vivo, DKK1-SE deficiency not only inhibited tumor proliferation but also reduced the complexity of the tumor microenvironment. This study identifies that DKK1-SE drives DKK1 expression by recruiting AP1 transcription factors, exerting oncogenic effects in PDAC, and enhancing the complexity of the tumor microenvironment.

Endogenous retroviruses mediate transcriptional rewiring in response to oncogenic signaling in colorectal cancer

Cancer cells exhibit rewired transcriptional regulatory networks that promote tumor growth and survival. However, the mechanisms underlying the formation of these pathological networks remain poorly understood. Through a pan-cancer epigenomic analysis, we found that primate-specific endogenous retroviruses (ERVs) are a rich source of enhancers displaying cancer-specific activity. In colorectal cancer and other epithelial tumors, oncogenic MAPK/AP1 signaling drives the activation of enhancers derived from the primate-specific ERV family LTR10. Functional studies in colorectal cancer cells revealed that LTR10 elements regulate tumor-specific expression of multiple genes associated with tumorigenesis, such as ATG12 and XRCC4. Within the human population, individual LTR10 elements exhibit germline and somatic structural variation resulting from a highly mutable internal tandem repeat region, which affects AP1 binding activity. Our findings reveal that ERV-derived enhancers contribute to transcriptional dysregulation in response to oncogenic signaling and shape the evolution of cancer-specific regulatory networks.

The epigenetic landscape in intestinal stem cells and its deregulation in colorectal cancer

Epigenetic mechanisms play a pivotal role in controlling gene expression and cellular plasticity in both normal physiology and pathophysiological conditions. These mechanisms are particularly important in the regulation of stem cell self-renewal and differentiation, both in embryonic development and within adult tissues. A prime example of this finely tuned epigenetic control is observed in the gastrointestinal lining, where the small intestine undergoes renewal approximately every 3-5 days. How various epigenetic mechanisms modulate chromatin functions in intestinal stem cells (ISCs) is currently an active area of research. In this review, we discuss the main epigenetic mechanisms that control ISC differentiation under normal homeostasis. Furthermore, we explore the dysregulation of these mechanisms in the context of colorectal cancer (CRC) development. By outlining the main epigenetic mechanisms contributing to CRC, we highlight the recent therapeutics development and future directions for colorectal cancer research.

Comprehensive multimodal and multiomic profiling reveals epigenetic and transcriptional reprogramming in lung tumors

Epigenomic mechanisms are critically involved in mediation of genetic and environmental factors that underlie cancer development. Histone modifications represent highly informative epigenomic marks that reveal activation and repression of gene activities and dysregulation of transcriptional control due to tumorigenesis. Here, we present a comprehensive epigenomic and transcriptomic mapping of 18 tumor and 20 non-neoplastic tissues from non-small cell lung adenocarcinoma patients. Our profiling covers 5 histone marks including activating (H3K4me3, H3K4me1, and H3K27ac) and repressive (H3K27me3 and H3K9me3) marks and the transcriptome using only 20 mg of tissue per sample, enabled by low-input omic technologies. Using advanced integrative bioinformatic analysis, we uncovered cancer-driving signaling cascade networks, changes in 3D genome modularity, and differential expression and functionalities of transcription factors and noncoding RNAs. Many of these identified genes and regulatory molecules showed no significant change in their expression or a single epigenomic modality, emphasizing the power of integrative multimodal and multiomic analysis using patient samples.

Genome-wide interaction study of dietary intake of fibre, fruits, and vegetables with risk of colorectal cancer

BACKGROUND

Consumption of fibre, fruits and vegetables have been linked with lower colorectal cancer (CRC) risk. A genome-wide gene-environment (G × E) analysis was performed to test whether genetic variants modify these associations.

METHODS

A pooled sample of 45 studies including up to 69,734 participants (cases: 29,896; controls: 39,838) of European ancestry were included. To identify G × E interactions, we used the traditional 1--degree-of-freedom (DF) G × E test and to improve power a 2-step procedure and a 3DF joint test that investigates the association between a genetic variant and dietary exposure, CRC risk and G × E interaction simultaneously.

FINDINGS

The 3-DF joint test revealed two significant loci with p-value <5 × 10-8. Rs4730274 close to the SLC26A3 gene showed an association with fibre (p-value: 2.4 × 10-3) and G × fibre interaction with CRC (OR per quartile of fibre increase = 0.87, 0.80, and 0.75 for CC, TC, and TT genotype, respectively; G × E p-value: 1.8 × 10-7). Rs1620977 in the NEGR1 gene showed an association with fruit intake (p-value: 1.0 × 10-8) and G × fruit interaction with CRC (OR per quartile of fruit increase = 0.75, 0.65, and 0.56 for AA, AG, and GG genotype, respectively; G × E -p-value: 0.029).

INTERPRETATION

We identified 2 loci associated with fibre and fruit intake that also modify the association of these dietary factors with CRC risk. Potential mechanisms include chronic inflammatory intestinal disorders, and gut function. However, further studies are needed for mechanistic validation and replication of findings.

FUNDING

National Institutes of Health, National Cancer Institute. Full funding details for the individual consortia are provided in acknowledgments.

Two genome-wide interaction loci modify the association of nonsteroidal anti-inflammatory drugs with colorectal cancer

Regular, long-term aspirin use may act synergistically with genetic variants, particularly those in mechanistically relevant pathways, to confer a protective effect on colorectal cancer (CRC) risk. We leveraged pooled data from 52 clinical trial, cohort, and case-control studies that included 30,806 CRC cases and 41,861 controls of European ancestry to conduct a genome-wide interaction scan between regular aspirin/nonsteroidal anti-inflammatory drug (NSAID) use and imputed genetic variants. After adjusting for multiple comparisons, we identified statistically significant interactions between regular aspirin/NSAID use and variants in 6q24.1 (top hit rs72833769), which has evidence of influencing expression of TBC1D7 (a subunit of the TSC1-TSC2 complex, a key regulator of MTOR activity), and variants in 5p13.1 (top hit rs350047), which is associated with expression of PTGER4 (codes a cell surface receptor directly involved in the mode of action of aspirin). Genetic variants with functional impact may modulate the chemopreventive effect of regular aspirin use, and our study identifies putative previously unidentified targets for additional mechanistic interrogation.

Fusobacterium nucleatum infection modulates the transcriptome and epigenome of HCT116 colorectal cancer cells in an oxygen-dependent manner

Fusobacterium nucleatum, a gram-negative oral bacterium, has been consistently validated as a strong contributor to the progression of several types of cancer, including colorectal (CRC) and pancreatic cancer. While previous in vitro studies have shown that intracellular F. nucleatum enhances malignant phenotypes such as cell migration, the dependence of this regulation on features of the tumor microenvironment (TME) such as oxygen levels are wholly uncharacterized. Here we examine the influence of hypoxia in facilitating F. nucleatum invasion and its effects on host responses focusing on changes in the global epigenome and transcriptome. Using a multiomic approach, we analyze epigenomic alterations of H3K27ac and global transcriptomic alterations sustained within a hypoxia and normoxia conditioned CRC cell line HCT116 at 24 h following initial infection with F. nucleatum. Our findings reveal that intracellular F. nucleatum activates signaling pathways and biological processes in host cells similar to those induced upon hypoxia conditioning in the absence of infection. Furthermore, we show that a hypoxic TME favors F. nucleatum invasion and persistence and therefore infection under hypoxia may amplify malignant transformation by exacerbating the effects induced by hypoxia alone. These results motivate future studies to investigate host-microbe interactions in tumor tissue relevant conditions that more accurately define parameters for targeted cancer therapies.

Super enhancer lncRNAs: a novel hallmark in cancer

Super enhancers (SEs) consist of clusters of enhancers, harboring an unusually high density of transcription factors, mediator coactivators and epigenetic modifications. SEs play a crucial role in the maintenance of cancer cell identity and promoting oncogenic transcription. Super enhancer lncRNAs (SE-lncRNAs) refer to either transcript from SEs locus or interact with SEs, whose transcriptional activity is highly dependent on SEs. Moreover, these SE-lncRNAs can interact with their associated enhancer regions in cis and modulate the expression of oncogenes or key signal pathways in cancers. Inhibition of SEs would be a promising therapy for cancer. In this review, we summarize the research of SE-lncRNAs in different kinds of cancers so far and decode the mechanism of SE-lncRNAs in carcinogenesis to provide novel ideas for the cancer therapy.

Targeting super-enhancer activity for colorectal cancer therapy

In addition to genetic variants and copy number alterations, epigenetic deregulation of oncogenes and tumor suppressors is a major contributor in cancer development and propagation. Regulatory elements for gene transcription regulation can be found in promoters which are located in the vicinity of transcription start sites but also at a distance, in enhancer sites, brought to interact with proximal sites when occupied by enhancer protein complexes. These sites provide most of the specific regulatory sequences recognized by transcription factors. A sub-set of enhancers characterized by a longer structure and stronger activity, called super-enhancers, are critical for the expression of specific genes, usually associated with individual cell type identity and function. Super-enhancers show deregulation in cancer, which may have profound repercussions for cancer cell survival and response to therapy. Dysfunction of super-enhancers may result from multiple mechanisms that include changes in their sequence, alterations in the topological neighborhoods where they belong, and alterations in the proteins that mediate their function, such as transcription factors and epigenetic modifiers. These can become potential targets for therapeutic interventions. Genes that are targets of super-enhancers are cell and cancer type specific and could also be of interest for therapeutic targeting. In colorectal cancer, a super-enhancer regulated and over-expressed oncogene is MYC, under the influence of the WNT/β-catenin pathway. Identification and targeting of additional oncogenes regulated by super-enhancers in colorectal cancer may pave the way for combination therapies targeting the super-enhancer machinery and signal transduction pathways that regulate the specific transcription factors operative on them.

TSCRE: a comprehensive database for tumor-specific cis-regulatory elements

Cis-regulatory elements (CREs) and super cis-regulatory elements (SCREs) are non-coding DNA regions which influence the transcription of nearby genes and play critical roles in development. Dysregulated CRE and SCRE activities have been reported to alter the expression of oncogenes and tumor suppressors, thereby regulating cancer hallmarks. To address the strong need for a comprehensive catalogue of dysregulated CREs and SCREs in human cancers, we present TSCRE (http://tscre.zsqylab.com/), an open resource providing tumor-specific and cell type-specific CREs and SCREs derived from the re-analysis of publicly available histone modification profiles. Currently, TSCRE contains 1 864 941 dysregulated CREs and 68 253 dysregulated SCREs identified from 1366 human patient samples spanning 17 different cancer types and 9 histone marks. Over 95% of these elements have been validated in public resources. TSCRE offers comprehensive annotations for each element, including associated genes, expression patterns, clinical prognosis, somatic mutations, transcript factor binding sites, cancer-type specificity, and drug response. Additionally, TSCRE integrates pathway and transcript factor enrichment analyses for each study, enabling in-depth functional and mechanistic investigations. Furthermore, TSCRE provides an interactive interface for users to explore any CRE and SCRE of interest. We believe TSCRE will be a highly valuable platform for the community to discover candidate cancer biomarkers.

Enhancer Recognition: A Transformer Encoder-Based Method with WGAN-GP for Data Augmentation

Enhancers are located upstream or downstream of key deoxyribonucleic acid (DNA) sequences in genes and can adjust the transcription activity of neighboring genes. Identifying enhancers and determining their functions are important for understanding gene regulatory networks and expression regulatory mechanisms. However, traditional enhancer recognition relies on manual feature engineering, which is time-consuming and labor-intensive, making it difficult to perform large-scale recognition analysis. In addition, if the original dataset is too small, there is a risk of overfitting. In recent years, emerging methods, such as deep learning, have provided new insights for enhancing identification. However, these methods also present certain challenges. Deep learning models typically require a large amount of high-quality data, and data acquisition demands considerable time and resources. To address these challenges, in this paper, we propose a data-augmentation method based on generative adversarial networks to solve the problem of small datasets. Moreover, we used regularization methods such as weight decay to improve the generalizability of the model and alleviate overfitting. The Transformer encoder was used as the main component to capture the complex relationships and dependencies in enhancer sequences. The encoding layer was designed based on the principle of k-mers to preserve more information from the original DNA sequence. Compared with existing methods, the proposed approach made significant progress in enhancing the accuracy and strength of enhancer identification and prediction, demonstrating the effectiveness of the proposed method. This paper provides valuable insights for enhancer analysis and is of great significance for understanding gene regulatory mechanisms and studying disease correlations.

Cooperation of MLL1 and Jun in controlling H3K4me3 on enhancers in colorectal cancer

BACKGROUND

Enhancer dysregulation is one of the important features for cancer cells. Enhancers enriched with H3K4me3 have been implicated to play important roles in cancer. However, their detailed features and regulatory mechanisms have not been well characterized.

RESULTS

Here, we profile the landscape of H3K4me3-enriched enhancers (m3Es) in 43 pairs of colorectal cancer (CRC) samples. M3Es are widely distributed in CRC and averagely possess around 10% of total active enhancers. We identify 1322 gain variant m3Es and 367 lost variant m3Es in CRC. The target genes of the gain m3Es are enriched in immune response pathways. We experimentally prove that repression of CBX8 and RPS6KA5 m3Es inhibits target gene expression in CRC. Furthermore, we find histone methyltransferase MLL1 is responsible for depositing H3K4me3 on the identified Vm3Es. We demonstrate that the transcription factor AP1/JUN interacts with MLL1 and regulates m3E activity. Application of a small chemical inhibitor for MLL1 activity, OICR-9429, represses target gene expression of the identified Vm3Es, enhances anti-tumor immunity and inhibits CRC growth in an animal model.

CONCLUSIONS

Taken together, our study illustrates the genome-wide landscape and the regulatory mechanisms of m3Es in CRC, and reveals potential novel strategies for cancer treatment.

Genome-wide interaction analysis of folate for colorectal cancer risk

BACKGROUND

Epidemiological and experimental evidence suggests that higher folate intake is associated with decreased colorectal cancer (CRC) risk; however, the mechanisms underlying this relationship are not fully understood. Genetic variation that may have a direct or indirect impact on folate metabolism can provide insights into folate's role in CRC.

OBJECTIVES

Our aim was to perform a genome-wide interaction analysis to identify genetic variants that may modify the association of folate on CRC risk.

METHODS

We applied traditional case-control logistic regression, joint 3-degree of freedom, and a 2-step weighted hypothesis approach to test the interactions of common variants (allele frequency >1%) across the genome and dietary folate, folic acid supplement use, and total folate in relation to risk of CRC in 30,550 cases and 42,336 controls from 51 studies from 3 genetic consortia (CCFR, CORECT, GECCO).

RESULTS

Inverse associations of dietary, total folate, and folic acid supplement with CRC were found (odds ratio [OR]: 0.93; 95% confidence interval [CI]: 0.90, 0.96; and 0.91; 95% CI: 0.89, 0.94 per quartile higher intake, and 0.82 (95% CI: 0.78, 0.88) for users compared with nonusers, respectively). Interactions (P-interaction < 5×10-8) of folic acid supplement and variants in the 3p25.2 locus (in the region of Synapsin II [SYN2]/tissue inhibitor of metalloproteinase 4 [TIMP4]) were found using traditional interaction analysis, with variant rs150924902 (located upstream to SYN2) showing the strongest interaction. In stratified analyses by rs150924902 genotypes, folate supplementation was associated with decreased CRC risk among those carrying the TT genotype (OR: 0.82; 95% CI: 0.79, 0.86) but increased CRC risk among those carrying the TA genotype (OR: 1.63; 95% CI: 1.29, 2.05), suggesting a qualitative interaction (P-interaction = 1.4×10-8). No interactions were observed for dietary and total folate.

CONCLUSIONS

Variation in 3p25.2 locus may modify the association of folate supplement with CRC risk. Experimental studies and studies incorporating other relevant omics data are warranted to validate this finding.

Genome-wide study of genetic polymorphisms predictive for outcome from first-line oxaliplatin-based chemotherapy in colorectal cancer patients

We conducted the first large genome-wide association study to identify novel genetic variants that predict better (or poorer) prognosis in colorectal cancer patients receiving standard first-line oxaliplatin-based chemotherapy vs chemotherapy without oxaliplatin. We used data from two phase III trials, NCCTG N0147 and NCCTG N9741 and a population-based patient cohort, DACHS. Multivariable Cox proportional hazards models were employed, including an interaction term between each SNP and type of treatment for overall survival (OS) and progression-free survival. The analysis was performed for studies individually, and the results were combined using fixed-effect meta-analyses separately for resected stage III colon cancer (3098 patients from NCCTG N0147 and 549 patients from DACHS) and mCRC (505 patients from NCCTG N9741 and 437 patients from DACHS). We further performed gene-based analysis as well as in silico bioinformatics analysis for CRC-relevant functional genomic annotation of identified loci. In stage III colon cancer patients, a locus on chr22 (rs11912167) was associated with significantly poorer OS after oxaliplatin-based chemotherapy vs chemotherapy without oxaliplatin (Pinteraction  < 5 × 10-8 ). For mCRC patients, three loci on chr1 (rs1234556), chr12 (rs11052270) and chr15 (rs11858406) were found to be associated with differential OS (P < 5 × 10-7 ). The locus on chr1 located in the intronic region of RCSD1 was replicated in an independent cohort of 586 mCRC patients from ALGB/SWOG 80405 (Pinteraction  = .04). The GWA gene-based analysis yielded for RCSD1 the most significant association with differential OS in mCRC (P = 6.6 × 10-6 ). With further investigation into its biological mechanisms, this finding could potentially be used to individualize first-line treatment and improve clinical outcomes.

Refined read-out: The hUHRF1 Tandem-Tudor domain prefers binding to histone H3 tails containing K4me1 in the context of H3K9me2/3

UHRF1 is an essential chromatin protein required for DNA methylation maintenance, mammalian development, and gene regulation. We investigated the Tandem-Tudor domain (TTD) of human UHRF1 that is known to bind H3K9me2/3 histones and is a major driver of UHRF1 localization in cells. We verified binding to H3K9me2/3 but unexpectedly discovered stronger binding to H3 peptides and mononucleosomes containing K9me2/3 with additional K4me1. We investigated the combined binding of TTD to H3K4me1-K9me2/3 versus H3K9me2/3 alone, engineered mutants with specific and differential changes of binding, and discovered a novel read-out mechanism for H3K4me1 in an H3K9me2/3 context that is based on the interaction of R207 with the H3K4me1 methyl group and on counting the H-bond capacity of H3K4. Individual TTD mutants showed up to a 10,000-fold preference for the double-modified peptides, suggesting that after a conformational change, WT TTD could exhibit similar effects. The frequent appearance of H3K4me1-K9me2 regions in human chromatin demonstrated in our TTD chromatin pull-down and ChIP-western blot data suggests that it has specific biological roles. Chromatin pull-down of TTD from HepG2 cells and full-length murine UHRF1 ChIP-seq data correlate with H3K4me1 profiles indicating that the H3K4me1-K9me2/3 interaction of TTD influences chromatin binding of full-length UHRF1. We demonstrate the H3K4me1-K9me2/3 specific binding of UHRF1-TTD to enhancers and promoters of cell-type-specific genes at the flanks of cell-type-specific transcription factor binding sites, and provided evidence supporting an H3K4me1-K9me2/3 dependent and TTD mediated downregulation of these genes by UHRF1. All these findings illustrate the important physiological function of UHRF1-TTD binding to H3K4me1-K9me2/3 double marks in a cellular context.

A bioinformatics analysis of potential cellular communication networks in non-alcoholic steatohepatitis and colorectal adenoma using scRNA-seq and bulk-seq

Background

Non-alcoholic fatty liver disease (NAFLD) is the global most common chronic liver disease. Non-alcoholic steatohepatitis (NASH), an inflammatory subtype of NAFLD, has been shown to significantly increase the risk of colorectal adenoma (CRA). Therefore, from the perspective of bioinformatics analysis, the potential mechanisms of NASH/NAFLD-CRA can be explored.

Methods

In this study, we screened the differentially expressed genes (DEGs) and core effect pathways between NASH and CRA by analyzing the single-cell data of CRA patients and the high-throughput sequencing data (GSE37364 and GSE89632) in the online database. We screened therapeutic targets and biomarkers through gene function classification, pathway enrichment analysis, and protein-protein interaction network analysis. In terms of single cell data, we screened the core effect pathway and specific signal pathway of cell communication through cell annotation and cell communication analyses. The purpose of the study was to find potential biomarkers, therapeutic targets, and related effect pathways of NASH-CRA.

Results

NASH-CRA comorbidities were concentrated in inflammatory regulation-related pathways, and the core genes of disease progression included IL1B, FOSL1, EGR1, MYC, PTGS2, and FOS. The results suggested the key pathway of NASH-CRA might be the WNT pathway. The main cell signal communication pathways included WNT2B - (FZD6 + LRP5) and WNT2B - (FZD6 + LRP6). The send-receive process occurred in embryonic stem cells.

Conclusions

The core genes of NASH-CRA (FOS, EGR1, MYC, PTGS2, FOSL1, and IL1B) may participate in inflammation and immune responses through up-regulation in the process of disease occurrence, interfering with the pathophysiological process of CRA and NASH. NASH-CRA produces cell signal communication in the WNT pathway sent by WNT2B and received by FZD6, LRP5, and LRP6 in embryonic stem cells. These findings may help formulate early diagnosis and treatment strategies for CRA in NAFLD/NASH patients, and further explore corresponding prognostic markers and potential approaches. The significance of scRNA-seq in exploring tumor heterogeneity lies in promoting our understanding of the expression program of tumor related genes in tumor development patterns. However, the biggest challenge is that this analysis may miss out on some biologically significant gene expression programs.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses

BACKGROUND

Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis.

METHODS

We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test).

RESULTS

Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09).

DISCUSSION

These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Super-enhancer landscape rewiring in cancer: The epigenetic control at distal sites

Super-enhancers evolve as elements at the top of the hierarchical control of gene expression. They are important end-gatherers of signaling pathways that control stemness, differentiation or adaptive responses. Many epigenetic regulations focus on these regions, and not surprisingly, during the process of tumorigenesis, various alterations can account for their dysfunction. Super-enhancers are emerging as key drivers of the aberrant gene expression landscape that sustain the aggressiveness of cancer cells. In this review, we will describe and discuss about the structure of super-enhancers, their epigenetic regulation, and the major changes affecting their functionality in cancer.

Genome-wide Interaction Study with Smoking for Colorectal Cancer Risk Identifies Novel Genetic Loci Related to Tumor Suppression, Inflammation, and Immune Response

BACKGROUND

Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer.

METHODS

A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia.

RESULTS

Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10-8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20-1.30] compared with the other genotypes (OR <1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10-8) and 8q24.23 (rs7005722, P = 2.88 × 10-8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09-1.16) compared with the other genotypes (OR <1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07-1.28) compared with the other genotypes (OR <1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33).

CONCLUSIONS

Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response.

IMPACT

These findings can guide potential prevention treatments.

Epigenetic effects of high-fat diet on intestinal tumorigenesis in C57BL/6J-Apc Min/+ mice

Aim

Obesity and obesogenic diets might partly accelerate cancer development through epigenetic mechanisms. To determine these early effects, we investigated the impact of three days of a high-fat diet on epigenomic and transcriptomic changes in Apc ^Min/+ murine intestinal epithelia.

Method

ChIP-Seq and RNA-Seq were performed on small intestinal epithelia of WT and Apc ^Min/+ male mice fed high-fat diet (HFD) or low-fat diet (LFD) for three days to identify genomic regions associated with differential H3K27ac levels as a marker of variant enhancer loci (VELs) as well as differentially expressed genes (DEGs).

Results

Regarding epigenetic and transcriptomic changes, diet type (LFD vs. HFD) showed a significant impact, and genotype (WT vs.Apc ^Min/+) showed a small impact. Compared to LFD, HFD resulted in 1306 gained VELs, 230 lost VELs, 133 upregulated genes, and 127 downregulated genes in WT mice, with 1056 gained VELs, 371 lost VELs, 222 upregulated genes, and 182 downregulated genes in Apc ^Min/+ mice. Compared to the WT genotype, the Apc ^Min/+ genotype resulted in zero changed VELs for either diet type group, 21 DEGs for LFD, and 48 DEGs for HFD. Most gained VELs, and upregulated genes were associated with lipid metabolic processes. Gained VELs were also associated with Wnt signaling. Downregulated genes were associated with antigen presentation and processing.

Conclusion

Three days of HFD-induced epigenomic and transcriptomic changes involving metabolic and immunologic pathways that may promote tumor growth in the genetically predisposed murine intestine without affecting key cancer signaling pathways.

Regulatory enhancer profiling of mesenchymal-type gastric cancer reveals subtype-specific epigenomic landscapes and targetable vulnerabilities

OBJECTIVE

Gastric cancer (GC) comprises multiple molecular subtypes. Recent studies have highlighted mesenchymal-subtype GC (Mes-GC) as a clinically aggressive subtype with few treatment options. Combining multiple studies, we derived and applied a consensus Mes-GC classifier to define the Mes-GC enhancer landscape revealing disease vulnerabilities.

DESIGN

Transcriptomic profiles of ~1000 primary GCs and cell lines were analysed to derive a consensus Mes-GC classifier. Clinical and genomic associations were performed across >1200 patients with GC. Genome-wide epigenomic profiles (H3K27ac, H3K4me1 and assay for transposase-accessible chromatin with sequencing (ATAC-seq)) of 49 primary GCs and GC cell lines were generated to identify Mes-GC-specific enhancer landscapes. Upstream regulators and downstream targets of Mes-GC enhancers were interrogated using chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing, CRISPR/Cas9 editing, functional assays and pharmacological inhibition.

RESULTS

We identified and validated a 993-gene cancer-cell intrinsic Mes-GC classifier applicable to retrospective cohorts or prospective single samples. Multicohort analysis of Mes-GCs confirmed associations with poor patient survival, therapy resistance and few targetable genomic alterations. Analysis of enhancer profiles revealed a distinctive Mes-GC epigenomic landscape, with TEAD1 as a master regulator of Mes-GC enhancers and Mes-GCs exhibiting preferential sensitivity to TEAD1 pharmacological inhibition. Analysis of Mes-GC super-enhancers also highlighted NUAK1 kinase as a downstream target, with synergistic effects observed between NUAK1 inhibition and cisplatin treatment.

CONCLUSION

Our results establish a consensus Mes-GC classifier applicable to multiple transcriptomic scenarios. Mes-GCs exhibit a distinct epigenomic landscape, and TEAD1 inhibition and combinatorial NUAK1 inhibition/cisplatin may represent potential targetable options.

A distal super-enhancer activates oncogenic ETS2 via recruiting MECOM in inflammatory bowel disease and colorectal cancer

Abnormal activities of distal cis-regulatory elements (CREs) contribute to the initiation and progression of cancer. Gain of super-enhancer (SE), a highly active distal CRE, is essential for the activation of key oncogenes in various cancers. However, the mechanism of action for most tumor-specific SEs still largely remains elusive. Here, we report that a candidate oncogene ETS2 was activated by a distal SE in inflammatory bowel disease (IBD) and colorectal cancer (CRC). The SE physically interacted with the ETS2 promoter and was required for the transcription activation of ETS2. Strikingly, the ETS2-SE activity was dramatically upregulated in both IBD and CRC tissues when compared to normal colon controls and was strongly correlated with the level of ETS2 expression. The tumor-specific activation of ETS2-SE was further validated by increased enhancer RNA transcription from this region in CRC. Intriguingly, a known IBD-risk SNP resides in the ETS2-SE and the genetic variant modulated the level of ETS2 expression through affecting the binding of an oncogenic transcription factor MECOM. Silencing of MECOM induced significant downregulation of ETS2 in CRC cells, and the level of MECOM and ETS2 correlated well with each other in CRC and IBD samples. Functionally, MECOM and ETS2 were both required for maintaining the colony-formation and sphere-formation capacities of CRC cells and MECOM was crucial for promoting migration. Taken together, we uncovered a novel disease-specific SE that distantly drives oncogenic ETS2 expression in IBD and CRC and delineated a mechanistic link between non-coding genetic variation and epigenetic regulation of gene transcription.

Systematic discovery and functional dissection of enhancers needed for cancer cell fitness and proliferation

A scarcity of functionally validated enhancers in the human genome presents a significant hurdle to understanding how these cis-regulatory elements contribute to human diseases. We carry out highly multiplexed CRISPR-based perturbation and sequencing to identify enhancers required for cell proliferation and fitness in 10 human cancer cell lines. Our results suggest that the cell fitness enhancers, unlike their target genes, display high cell-type specificity of chromatin features. They typically adopt a modular structure, comprised of activating elements enriched for motifs of oncogenic transcription factors, surrounded by repressive elements enriched for motifs recognized by transcription factors with tumor suppressor functions. We further identify cell fitness enhancers that are selectively accessible in clinical tumor samples, and the levels of chromatin accessibility are associated with patient survival. These results reveal functional enhancers across multiple cancer cell lines, characterize their context-dependent chromatin organization, and yield insights into altered transcription programs in cancer cells.

A local tumor microenvironment acquired super-enhancer induces an oncogenic driver in colorectal carcinoma

Tumors exhibit enhancer reprogramming compared to normal tissue. The etiology is largely attributed to cell-intrinsic genomic alterations. Here, using freshly resected primary CRC tumors and patient-matched adjacent normal colon, we find divergent epigenetic landscapes between CRC tumors and cell lines. Intriguingly, this phenomenon extends to highly recurrent aberrant super-enhancers gained in CRC over normal. We find one such super-enhancer activated in epithelial cancer cells due to surrounding inflammation in the tumor microenvironment. We restore this super-enhancer and its expressed gene, PDZK1IP1, following treatment with cytokines or xenotransplantation into nude mice, thus demonstrating cell-extrinsic etiology. We demonstrate mechanistically that PDZK1IP1 enhances the reductive capacity CRC cancer cells via the pentose phosphate pathway. We show this activation enables efficient growth under oxidative conditions, challenging the previous notion that PDZK1IP1 acts as a tumor suppressor in CRC. Collectively, these observations highlight the significance of epigenomic profiling on primary specimens.

Multi-Omics Analysis of Molecular Characteristics and Carcinogenic Effect of NFE2L3 in Pan-Cancer

NFE2L3, also known as NFE2L3, is a nuclear transcription factor associated with the pathogenesis and progression of human tumors. To systematically and comprehensively investigate the role of NFE2L3 in tumors, a pan-cancer analysis was performed using multi-omics data, including gene expression analysis, diagnostic and prognostic analysis, epigenetic methylation analysis, gene alteration analysis, immune feature analysis, functional enrichment analysis, and tumor cell functional status analysis. Furthermore, the molecular mechanism of NFE2L3 in liver hepatocellular carcinoma (LIHC) was explored. The relationship between NFE2L3 expression and survival prognosis of patients with LIHC was analyzed and a nomogram prediction model was constructed. Our study showed that NFE2L3 expression was upregulated in most cancers, suggesting that NFE2L3 may play an important role in promoting cancer progression. NFE2L3 expression is closely related to DNA methylation, genetic alteration, immune signature, and tumor cell functional status in pan-cancers. Furthermore, NFE2L3 was demonstrated to be an independent risk factor for LIHC, and the nomogram model based on NFE2L3 expression had good prediction efficiency for the overall survival of patients with LIHC. In summary, our study indicated that NFE2L3 may be an important molecular biomarker for the diagnosis and prognosis of pan-cancer. NFE2L3 is expected to be a potential molecular target for the treatment of tumors.

Beyond GWAS of Colorectal Cancer: Evidence of Interaction with Alcohol Consumption and Putative Causal Variant for the 10q24.2 Region

BACKGROUND

Currently known associations between common genetic variants and colorectal cancer explain less than half of its heritability of 25%. As alcohol consumption has a J-shape association with colorectal cancer risk, nondrinking and heavy drinking are both risk factors for colorectal cancer.

METHODS

Individual-level data was pooled from the Colon Cancer Family Registry, Colorectal Transdisciplinary Study, and Genetics and Epidemiology of Colorectal Cancer Consortium to compare nondrinkers (≤1 g/day) and heavy drinkers (>28 g/day) with light-to-moderate drinkers (1-28 g/day) in GxE analyses. To improve power, we implemented joint 2df and 3df tests and a novel two-step method that modifies the weighted hypothesis testing framework. We prioritized putative causal variants by predicting allelic effects using support vector machine models.

RESULTS

For nondrinking as compared with light-to-moderate drinking, the hybrid two-step approach identified 13 significant SNPs with pairwise r2 > 0.9 in the 10q24.2/COX15 region. When stratified by alcohol intake, the A allele of lead SNP rs2300985 has a dose-response increase in risk of colorectal cancer as compared with the G allele in light-to-moderate drinkers [OR for GA genotype = 1.11; 95% confidence interval (CI), 1.06-1.17; OR for AA genotype = 1.22; 95% CI, 1.14-1.31], but not in nondrinkers or heavy drinkers. Among the correlated candidate SNPs in the 10q24.2/COX15 region, rs1318920 was predicted to disrupt an HNF4 transcription factor binding motif.

CONCLUSIONS

Our study suggests that the association with colorectal cancer in 10q24.2/COX15 observed in genome-wide association study is strongest in nondrinkers. We also identified rs1318920 as the putative causal regulatory variant for the region.

IMPACT

The study identifies multifaceted evidence of a possible functional effect for rs1318920.

Chromatin state dynamics confers specific therapeutic strategies in enhancer subtypes of colorectal cancer

OBJECTIVE

Enhancer aberrations are beginning to emerge as a key epigenetic feature of colorectal cancers (CRC), however, a comprehensive knowledge of chromatin state patterns in tumour progression, heterogeneity of these patterns and imparted therapeutic opportunities remain poorly described.

DESIGN

We performed comprehensive epigenomic characterisation by mapping 222 chromatin profiles from 69 samples (33 colorectal adenocarcinomas, 4 adenomas, 21 matched normal tissues and 11 colon cancer cell lines) for six histone modification marks: H3K4me3 for Pol II-bound and CpG-rich promoters, H3K4me1 for poised enhancers, H3K27ac for enhancers and transcriptionally active promoters, H3K79me2 for transcribed regions, H3K27me3 for polycomb repressed regions and H3K9me3 for heterochromatin.

RESULTS

We demonstrate that H3K27ac-marked active enhancer state could distinguish between different stages of CRC progression. By epigenomic editing, we present evidence that gains of tumour-specific enhancers for crucial oncogenes, such as ASCL2 and FZD10, was required for excessive proliferation. Consistently, combination of MEK plus bromodomain inhibition was found to have synergistic effects in CRC patient-derived xenograft models. Probing intertumour heterogeneity, we identified four distinct enhancer subtypes (EPIgenome-based Classification, EpiC), three of which correlate well with previously defined transcriptomic subtypes (consensus molecular subtypes, CMSs). Importantly, CMS2 can be divided into two EpiC subgroups with significant survival differences. Leveraging such correlation, we devised a combinatorial therapeutic strategy of enhancer-blocking bromodomain inhibitors with pathway-specific inhibitors (PARPi, EGFRi, TGFβi, mTORi and SRCi) for EpiC groups.

CONCLUSION

Our data suggest that the dynamics of active enhancer underlies CRC progression and the patient-specific enhancer patterns can be leveraged for precision combination therapy.

Nutritional Epigenetics in Cancer

Alterations in the epigenome are well known to affect cancer development and progression. Epigenetics is highly influenced by the environment, including diet, which is a source of metabolic substrates that influence the synthesis of cofactors or substrates for chromatin and RNA modifying enzymes. In addition, plants are a common source of bioactives that can directly modify the activity of these enzymes. Here, we review and discuss the impact of diet on epigenetic mechanisms, including chromatin and RNA regulation, and its potential implications for cancer prevention and treatment.

Recurrent Superenhancer of the Oncogene POU5F1B in Colorectal Cancers

Superenhancer usages in single cancer form such as colorectal cancer (CRC) may provide novel efficient targeting candidates. It is unclear whether CRC contains recurrent superenhancers that confer a predisposition to malignancy. We investigated the superenhancer profile of CRC cell line HCT116 and compared it to that of a healthy sigmoid colon. We found that HCT116 had lost most of the normal colon superenhancer activities but gained a new set of tumor-favoring superenhancers that facilitate tumor proliferation, growth signalling, and hypoxia resistance. Inhibiting the superenhancers by JQ-1 treatment had significantly decreased the colony formation capability of HCT116. Then, by comparing the superenhancer genes and robust CRC upregulated genes, we identified a superenhancer associated with a common CRC upregulated oncogene, POU5f1B. POU5f1B overexpression is related to the worse outcome in CRCs. Via performing ChIP-PCR in 35 clinical samples and investigating CRC anti-H3K27ac ChiP-seq public dataset consisting of 36 samples, we further identified that the superenhancer of oncogene POU5F1B is recurrently activated in CRCs, taking 62 and 72 per cent, respectively. Moreover, JQ-1 treatment successfully inhibited the POU5F1B expression in 5 out of 6 POU5F1B superenhancer-positive samples. Therefore, we concluded that the superenhancer activation of POU5F1B contributes partially to its high expression in CRCs, in addition to the well-known gene amplification aetiology.

A β-Catenin-TCF-Sensitive Locus Control Region Mediates GUCY2C Ligand Loss in Colorectal Cancer

BACKGROUND & AIMS

Sporadic colorectal cancers arise from initiating mutations in APC, producing oncogenic β-catenin/TCF-dependent transcriptional reprogramming. Similarly, the tumor suppressor axis regulated by the intestinal epithelial receptor GUCY2C is among the earliest pathways silenced in tumorigenesis. Retention of the receptor, but loss of its paracrine ligands, guanylin and uroguanylin, is an evolutionarily conserved feature of colorectal tumors, arising in the earliest dysplastic lesions. Here, we examined a mechanism of GUCY2C ligand transcriptional silencing by β-catenin/TCF signaling.

METHODS

We performed RNA sequencing analysis of 4 unique conditional human colon cancer cell models of β-catenin/TCF signaling to map the core Wnt-transcriptional program. We then performed a comparative analysis of orthogonal approaches, including luciferase reporters, chromatin immunoprecipitation sequencing, CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats) knockout, and CRISPR epigenome editing, which were cross-validated with human tissue chromatin immunoprecipitation sequencing datasets, to identify functional gene enhancers mediating GUCY2C ligand loss.

RESULTS

RNA sequencing analyses reveal the GUCY2C hormones as 2 of the most sensitive targets of β-catenin/TCF signaling, reflecting transcriptional repression. The GUCY2C hormones share an insulated genomic locus containing a novel locus control region upstream of the guanylin promoter that mediates the coordinated silencing of both genes. Targeting this region with CRISPR epigenome editing reconstituted GUCY2C ligand expression, overcoming gene inactivation by mutant β-catenin/TCF signaling.

CONCLUSIONS

These studies reveal DNA elements regulating corepression of GUCY2C ligand transcription by β-catenin/TCF signaling, reflecting a novel pathophysiological step in tumorigenesis. They offer unique genomic strategies that could reestablish hormone expression in the context of canonical oncogenic mutations to reconstitute the GUCY2C axis and oppose transformation.

Genome-wide profiling in colorectal cancer identifies PHF19 and TBC1D16 as oncogenic super enhancers

Colorectal cancer is one of the most common cancers in the world. Although genomic mutations and single nucleotide polymorphisms have been extensively studied, the epigenomic status in colorectal cancer patient tissues remains elusive. Here, together with genomic and transcriptomic analysis, we use ChIP-Seq to profile active enhancers at the genome wide level in colorectal cancer paired patient tissues (tumor and adjacent tissues from the same patients). In total, we sequence 73 pairs of colorectal cancer tissues and generate 147 H3K27ac ChIP-Seq, 144 RNA-Seq, 147 whole genome sequencing and 86 H3K4me3 ChIP-Seq samples. Our analysis identifies 5590 gain and 1100 lost variant enhancer loci in colorectal cancer, and 334 gain and 121 lost variant super enhancer loci. Multiple key transcription factors in colorectal cancer are predicted with motif analysis and core regulatory circuitry analysis. Further experiments verify the function of the super enhancers governing PHF19 and TBC1D16 in regulating colorectal cancer tumorigenesis, and KLF3 is identified as an oncogenic transcription factor in colorectal cancer. Taken together, our work provides an important epigenomic resource and functional factors for epigenetic studies in colorectal cancer.

Hepatocyte nuclear factor HNF1A is a potential regulator in shaping the super-enhancer landscape in colorectal cancer liver metastasis

Super-enhancers (SEs) play essential roles in colorectal cancer (CRC) progression. However, how the SE landscape is orchestrated by transcriptional regulators and evolves is not clear. Using de novo motif analysis, we show that the hepatocyte nuclear factor 1 (HNF1)-binding motif is enriched in SEs in cell lines derived from liver metastases, but not in those from primary tumors. This finding was further validated by extending the method to pancreatic cancer and a pair of isogenic CRC lines. Next, we revealed HNF1-alpha (HNF1A) was majorly expressed and upregulated in CRC liver metastatic cell lines. Clinically, HNF1A was remarkably upregulated in synchronous liver metastases as compared to localized tumors. Collectively, our study implicates HNF1A as a key regulator in shaping the SE landscape in CRC liver metastasis.

Role of Enhancers in Development and Diseases

Enhancers are cis-regulatory elements containing short DNA sequences that serve as binding sites for pioneer/regulatory transcription factors, thus orchestrating the regulation of genes critical for lineage determination. The activity of enhancer elements is believed to be determined by transcription factor binding, thus determining the cell state identity during development. Precise spatio-temporal control of the transcriptome during lineage specification requires the coordinated binding of lineage-specific transcription factors to enhancers. Thus, enhancers are the primary determinants of cell identity. Numerous studies have explored the role and mechanism of enhancers during development and disease, and various basic questions related to the functions and mechanisms of enhancers have not yet been fully answered. In this review, we discuss the recently published literature regarding the roles of enhancers, which are critical for various biological processes governing development. Furthermore, we also highlight that altered enhancer landscapes provide an essential context to understand the etiologies and mechanisms behind numerous complex human diseases, providing new avenues for effective enhancer-based therapeutic interventions.

Integrated Analysis of Whole Genome and Epigenome Data Using Machine Learning Technology: Toward the Establishment of Precision Oncology

With the completion of the International Human Genome Project, we have entered what is known as the post-genome era, and efforts to apply genomic information to medicine have become more active. In particular, with the announcement of the Precision Medicine Initiative by U.S. President Barack Obama in his State of the Union address at the beginning of 2015, "precision medicine," which aims to divide patients and potential patients into subgroups with respect to disease susceptibility, has become the focus of worldwide attention. The field of oncology is also actively adopting the precision oncology approach, which is based on molecular profiling, such as genomic information, to select the appropriate treatment. However, the current precision oncology is dominated by a method called targeted-gene panel (TGP), which uses next-generation sequencing (NGS) to analyze a limited number of specific cancer-related genes and suggest optimal treatments, but this method causes the problem that the number of patients who benefit from it is limited. In order to steadily develop precision oncology, it is necessary to integrate and analyze more detailed omics data, such as whole genome data and epigenome data. On the other hand, with the advancement of analysis technologies such as NGS, the amount of data obtained by omics analysis has become enormous, and artificial intelligence (AI) technologies, mainly machine learning (ML) technologies, are being actively used to make more efficient and accurate predictions. In this review, we will focus on whole genome sequencing (WGS) analysis and epigenome analysis, introduce the latest results of omics analysis using ML technologies for the development of precision oncology, and discuss the future prospects.

A pan-cancer analysis of CpG Island gene regulation reveals extensive plasticity within Polycomb target genes

CpG Island promoter genes make up more than half of human genes, and a subset regulated by Polycomb-Repressive Complex 2 (PRC2+-CGI) become DNA hypermethylated and silenced in cancer. Here, we perform a systematic analysis of CGI genes across TCGA cancer types, finding that PRC2+-CGI genes are frequently prone to transcriptional upregulation as well. These upregulated PRC2+-CGI genes control important pathways such as Epithelial-Mesenchymal Transition (EMT) and TNFα-associated inflammatory response, and have greater cancer-type specificity than other CGI genes. Using publicly available chromatin datasets and genetic perturbations, we show that transcription factor binding sites (TFBSs) within distal enhancers underlie transcriptional activation of PRC2+-CGI genes, coinciding with loss of the PRC2-associated mark H3K27me3 at the linked promoter. In contrast, PRC2-free CGI genes are predominantly regulated by promoter TFBSs which are common to most cancer types. Surprisingly, a large subset of PRC2+-CGI genes that are upregulated in one cancer type are also hypermethylated/silenced in at least one other cancer type, underscoring the high degree of regulatory plasticity of these genes, likely derived from their complex regulatory control during normal development.

Epigenomic landscape of human colorectal cancer unveils an aberrant core of pan-cancer enhancers orchestrated by YAP/TAZ

Cancer is characterized by pervasive epigenetic alterations with enhancer dysfunction orchestrating the aberrant cancer transcriptional programs and transcriptional dependencies. Here, we epigenetically characterize human colorectal cancer (CRC) using de novo chromatin state discovery on a library of different patient-derived organoids. By exploring this resource, we unveil a tumor-specific deregulated enhancerome that is cancer cell-intrinsic and independent of interpatient heterogeneity. We show that the transcriptional coactivators YAP/TAZ act as key regulators of the conserved CRC gained enhancers. The same YAP/TAZ-bound enhancers display active chromatin profiles across diverse human tumors, highlighting a pan-cancer epigenetic rewiring which at single-cell level distinguishes malignant from normal cell populations. YAP/TAZ inhibition in established tumor organoids causes extensive cell death unveiling their essential role in tumor maintenance. This work indicates a common layer of YAP/TAZ-fueled enhancer reprogramming that is key for the cancer cell state and can be exploited for the development of improved therapeutic avenues.

The role of epigenetic deregulation in colorectal cancer (CRC) is not fully understood yet. Here the authors use patient-derived organoids, epigenomics and single-cell RNA-seq to reveal that YAP/TAZ are key regulators that bind to active enhancers in CRC and promote tumour survival.

iEnhancer-GAN: A Deep Learning Framework in Combination with Word Embedding and Sequence Generative Adversarial Net to Identify Enhancers and Their Strength

As critical components of DNA, enhancers can efficiently and specifically manipulate the spatial and temporal regulation of gene transcription. Malfunction or dysregulation of enhancers is implicated in a slew of human pathology. Therefore, identifying enhancers and their strength may provide insights into the molecular mechanisms of gene transcription and facilitate the discovery of candidate drug targets. In this paper, a new enhancer and its strength predictor, iEnhancer-GAN, is proposed based on a deep learning framework in combination with the word embedding and sequence generative adversarial net (Seq-GAN). Considering the relatively small training dataset, the Seq-GAN is designed to generate artificial sequences. Given that each functional element in DNA sequences is analogous to a “word” in linguistics, the word segmentation methods are proposed to divide DNA sequences into “words”, and the skip-gram model is employed to transform the “words” into digital vectors. In view of the powerful ability to extract high-level abstraction features, a convolutional neural network (CNN) architecture is constructed to perform the identification tasks, and the word vectors of DNA sequences are vertically concatenated to form the embedding matrices as the input of the CNN. Experimental results demonstrate the effectiveness of the Seq-GAN to expand the training dataset, the possibility of applying word segmentation methods to extract “words” from DNA sequences, the feasibility of implementing the skip-gram model to encode DNA sequences, and the powerful prediction ability of the CNN. Compared with other state-of-the-art methods on the training dataset and independent test dataset, the proposed method achieves a significantly improved overall performance. It is anticipated that the proposed method has a certain promotion effect on enhancer related fields.

Epigenetic alterations in the gastrointestinal tract: Current and emerging use for biomarkers of cancer

Colorectal cancer is a leading cause of cancer related deaths worldwide. One of the hallmarks of cancer and a fundamental trait of virtually all gastrointestinal cancers is genomic and epigenomic DNA alterations. Cancer cells acquire genetic and epigenetic alterations that drive the initiation and progression of the cancers by altering the molecular and cell biological process of the cells. These alterations, as well as other host and microenvironment factors, ultimately mediate the initiation and progression of cancers, including colorectal cancer. Epigenetic alterations, which include changes affecting DNA methylation, histone modifications, chromatin structure, and noncoding RNA expression, have emerged as a major class of molecular alteration in colon polyps and colorectal cancer. The classes of epigenetic alterations, their status in colorectal polyps and cancer, their effects on neoplasm biology, and their application to clinical care will be discussed.

Multi-omics characterization of WNT pathway reactivation to ameliorate BET inhibitor resistance in liver cancer cells

The Bromodomain and Extra-terminal domain (BET) proteins are promising targets in treating cancers. Although BET inhibitors have been in clinical trials, they are limited by lacking of suitable biomarkers to indicate drug responses in different cancers. Here we identify DHRS2, ETV4 and NOTUM as potential biomarkers to indicate drug resistance in liver cancer cells of a recently discovered BET inhibitor, Hjp-6-171. Furthermore, we confirm that reactivation of WNT pathway, the target of NOTUM, contributes to the drug sensitivity restoration in Hjp-6-171 resistant cells. Specially, combinations of Hjp-6-171 and a GSK3β inhibitor CHIR-98014 show remarkable therapeutic effects in vitro and in vivo. Integrating RNA-seq and ChIP-seq data, we reveal the expression signature of β-catenin regulated genes is contrary in sensitive cells to that in resistant cells. We propose WNT signaling molecules such as β-catenin and ETV4 to be candidate biomarkers to indicate BET inhibitor responses in liver cancer patients.

Genetic Effects on Transcriptome Profiles in Colon Epithelium Provide Functional Insights for Genetic Risk Loci

BACKGROUND & AIMS

The association of genetic variation with tissue-specific gene expression and alternative splicing guides functional characterization of complex trait-associated loci and may suggest novel genes implicated in disease. Here, our aims were as follows: (1) to generate reference profiles of colon mucosa gene expression and alternative splicing and compare them across colon subsites (ascending, transverse, and descending), (2) to identify expression and splicing quantitative trait loci (QTLs), (3) to find traits for which identified QTLs contribute to single-nucleotide polymorphism (SNP)-based heritability, (4) to propose candidate effector genes, and (5) to provide a web-based visualization resource.

METHODS

We collected colonic mucosal biopsy specimens from 485 healthy adults and performed bulk RNA sequencing. We performed genome-wide SNP genotyping from blood leukocytes. Statistical approaches and bioinformatics software were used for QTL identification and downstream analyses.

RESULTS

We provided a complete quantification of gene expression and alternative splicing across colon subsites and described their differences. We identified thousands of expression and splicing QTLs and defined their enrichment at genome-wide regulatory regions. We found that part of the SNP-based heritability of diseases affecting colon tissue, such as colorectal cancer and inflammatory bowel disease, but also of diseases affecting other tissues, such as psychiatric conditions, can be explained by the identified QTLs. We provided candidate effector genes for multiple phenotypes. Finally, we provided the Colon Transcriptome Explorer web application.

CONCLUSIONS

We provide a large characterization of gene expression and splicing across colon subsites. Our findings provide greater etiologic insight into complex traits and diseases influenced by transcriptomic changes in colon tissue.

SNPs associated with colorectal cancer at 15q13.3 affect risk enhancers that modulate GREM1 gene expression

Several genome wide association studies of colorectal cancer (CRC) have identified single nucleotide polymorphisms (SNPs) on chromosome 15q13.3 associated with CRC risk. To identify functional variant(s) underlying this association, we investigated SNPs in linkage disequilibrium with the risk‐associated SNP rs4779584 that overlapped regulatory regions/enhancer elements characterized in colon‐related tissues and cells. We identified several SNP‐containing regulatory regions that exhibited enhancer activity in vitro, including one SNP (rs1406389) that correlated with allele‐specific effects on enhancer activity. Deletion of either this enhancer or another enhancer that had previously been reported in this region correlated with decreased expression of GREM1 following CRISPR/Cas9 genome editing. That GREM1 is one target of these enhancers was further supported by an expression quantitative trait loci correlation between rs1406389 and GREM1 expression in the transverse but not sigmoid colon in the Genotype‐Tissue Expression dataset. Taken together, we conclude that the 15q13.3 region contains at least two functional variants that map to distinct enhancers and impact CRC risk through modulation of GREM1 expression.

Super enhancers-Functional cores under the 3D genome

Complex biochemical reactions take place in the nucleus all the time. Transcription machines must follow the rules. The chromatin state, especially the three-dimensional structure of the genome, plays an important role in gene regulation and expression. The super enhancers are important for defining cell identity in mammalian developmental processes and human diseases. It has been shown that the major components of transcriptional activation complexes are recruited by super enhancer to form phase-separated condensates. We summarize the current knowledge about super enhancer in the 3D genome. Furthermore, a new related transcriptional regulation model from super enhancer is outlined to explain its role in the mammalian cell progress.

Epigenetic Alterations in the Gastrointestinal Tract: Current and Emerging Use for Biomarkers of Cancer

Colorectal cancer, liver cancer, stomach cancer, pancreatic cancer, and esophageal cancer are leading causes of cancer-related deaths worldwide. A fundamental trait of virtually all gastrointestinal cancers is genomic and epigenomic DNA alterations. Cancer cells acquire genetic and epigenetic alterations that drive the initiation and progression of the cancers by altering the molecular and cell biological processes of the cells. These alterations, as well as other host and microenvironment factors, ultimately mediate the clinical behavior of the precancers and cancers and can be used as biomarkers for cancer risk determination, early detection of cancer and precancer, determination of the prognosis of cancer and prediction of the response to therapy. Epigenetic alterations have emerged as one of most robust classes of biomarkers and are the basis for a growing number of clinical tests for cancer screening and surveillance.

An integrated multi-omics approach to identify regulatory mechanisms in cancer metastatic processes

BACKGROUND

Metastatic progress is the primary cause of death in most cancers, yet the regulatory dynamics driving the cellular changes necessary for metastasis remain poorly understood. Multi-omics approaches hold great promise for addressing this challenge; however, current analysis tools have limited capabilities to systematically integrate transcriptomic, epigenomic, and cistromic information to accurately define the regulatory networks critical for metastasis.

RESULTS

To address this limitation, we use a purposefully generated cellular model of colon cancer invasiveness to generate multi-omics data, including expression, accessibility, and selected histone modification profiles, for increasing levels of invasiveness. We then adopt a rigorous probabilistic framework for joint inference from the resulting heterogeneous data, along with transcription factor binding profiles. Our approach uses probabilistic graphical models to leverage the functional information provided by specific epigenomic changes, models the influence of multiple transcription factors simultaneously, and automatically learns the activating or repressive roles of cis-regulatory events. Global analysis of these relationships reveals key transcription factors driving invasiveness, as well as their likely target genes. Disrupting the expression of one of the highly ranked transcription factors JunD, an AP-1 complex protein, confirms functional relevance to colon cancer cell migration and invasion. Transcriptomic profiling confirms key regulatory targets of JunD, and a gene signature derived from the model demonstrates strong prognostic potential in TCGA colorectal cancer data.

CONCLUSIONS

Our work sheds new light into the complex molecular processes driving colon cancer metastasis and presents a statistically sound integrative approach to analyze multi-omics profiles of a dynamic biological process.

Potential Mechanism of Immune Evasion Associated with the Master Regulator ASCL2 in Microsatellite Stability in Colorectal Cancer

Colorectal cancer (CRC) has two major subtypes, microsatellite instability (MSI) and microsatellite stability (MSS) based on the genomic instability. In this study, using computational programs, we identified 9 master transcription factors (TFs) based on epigenomic profiling in MSS CRC samples. Notably, unbiased gene set enrichment analysis (GSEA) showed that several master TFs were strongly associated with immune-related functions in TCGA MSS CRC tissues, such as interferon gamma (IFN-γ) and interferon alpha (IFN-α) responses. Focusing to the top candidate, ASCL2, we found that CD8⁺ T cell infiltration was low in ASCL2 overexpressed MSS CRC samples. Compared with other gastrointestinal (GI) cancers (gastric cancer, MSI CRC, and esophageal cancer), ASCL2 is specifically upregulated in MSS CRC. Moreover, we identified 28 candidate genes in IFN-γ and IFN-α response pathways which were negatively correlated with ASCL2. Together, these results link transcriptional dysregulation with the immune evasion in MSS CRC, which may advance the understanding of immune resistance and contribute to developing novel treatments of MSS CRC.

Transcriptionally active enhancers in human cancer cells

The growth of human cancer cells is driven by aberrant enhancer and gene transcription activity. Here, we use transient transcriptome sequencing (TT-seq) to map thousands of transcriptionally active putative enhancers in fourteen human cancer cell lines covering seven types of cancer. These enhancers were associated with cell type-specific gene expression, enriched for genetic variants that predispose to cancer, and included functionally verified enhancers. Enhancer-promoter (E-P) pairing by correlation of transcription activity revealed ~ 40,000 putative E-P pairs, which were depleted for housekeeping genes and enriched for transcription factors, cancer-associated genes, and 3D conformational proximity. The cell type specificity and transcription activity of target genes increased with the number of paired putative enhancers. Our results represent a rich resource for future studies of gene regulation by enhancers and their role in driving cancerous cell growth.

Common variants in signaling transcription-factor-binding sites drive phenotypic variability in red blood cell traits

Genome-wide association studies identify genomic variants associated with human traits and diseases. Most trait-associated variants are located within cell-type-specific enhancers, but the molecular mechanisms governing phenotypic variation are less well understood. Here, we show that many enhancer variants associated with red blood cell (RBC) traits map to enhancers that are co-bound by lineage-specific master transcription factors (MTFs) and signaling transcription factors (STFs) responsive to extracellular signals. The majority of enhancer variants reside on STF and not MTF motifs, perturbing DNA binding by various STFs (BMP/TGF-β-directed SMADs or WNT-induced TCFs) and affecting target gene expression. Analyses of engineered human blood cells and expression quantitative trait loci verify that disrupted STF binding leads to altered gene expression. Our results propose that the majority of the RBC-trait-associated variants that reside on transcription-factor-binding sequences fall in STF target sequences, suggesting that the phenotypic variation of RBC traits could stem from altered responsiveness to extracellular stimuli.

The Multifaceted Output of c-Jun Biological Activity: Focus at the Junction of CD8 T Cell Activation and Exhaustion

c-Jun is a major component of the dimeric transcription factor activator protein-1 (AP-1), a paradigm for transcriptional response to extracellular signaling, whose components are basic-Leucine Zipper (bZIP) transcription factors of the Jun, Fos, activating transcription factor (ATF), ATF-like (BATF) and Jun dimerization protein 2 (JDP2) gene families. Extracellular signals regulate c-Jun/AP-1 activity at multiple levels, including transcriptional and posttranscriptional regulation of c-Jun expression and transactivity, in turn, establishing the magnitude and the duration of c-Jun/AP-1 activation. Another important level of c-Jun/AP-1 regulation is due to the capability of Jun family members to bind DNA as a heterodimer with every other member of the AP-1 family, and to interact with other classes of transcription factors, thereby acquiring the potential to integrate diverse extrinsic and intrinsic signals into combinatorial regulation of gene expression. Here, we review how these features of c-Jun/AP-1 regulation underlie the multifaceted output of c-Jun biological activity, eliciting quite distinct cellular responses, such as neoplastic transformation, differentiation and apoptosis, in different cell types. In particular, we focus on the current understanding of the role of c-Jun/AP-1 in the response of CD8 T cells to acute infection and cancer. We highlight the transcriptional and epigenetic regulatory mechanisms through which c-Jun/AP-1 participates in the productive immune response of CD8 T cells, and how its downregulation may contribute to the dysfunctional state of tumor infiltrating CD8 T cells. Additionally, we discuss recent insights pointing at c-Jun as a suitable target for immunotherapy-based combination approaches to reinvigorate anti-tumor immune functions.

Genetic Variant Associated With Survival of Patients With Stage II-III Colon Cancer

BACKGROUND & AIMS

Many genetic variants have been associated with colorectal cancer risk, although few have been associated with survival times of patients. Identification of genetic variants associated with survival times might improve our understanding of disease progression and aid in outcome prediction. We performed a genome-wide association study to identify variants associated with colon cancer survival time.

METHODS

We performed a post hoc analysis of data from NCCTG N0147 (Alliance), a randomized phase 3 trial of patients with resected stage III colon cancer, and from NSABP C-08 (NRG), a phase 3 trial that compared therapy regimens for patients with resected stage II or III colon cancer. Genotype analyses were performed on DNA from blood samples from 4974 patients. We used Cox proportional hazards regression to evaluate the association of each single nucleotide polymorphism with times of overall survival and disease-free survival, adjusting for age at diagnosis, sex, treatment group, and principal components of genetic ancestry. We performed the analysis for studies N0147 and C-08 separately, and results were combined in a fixed-effects meta-analysis.

RESULTS

A locus on chromosome 7p15.2 was significantly associated with overall survival time (P ≤ 5x10-08). The most significant variant at this locus, rs76766811 (P = 1.6x10-08), is common among African Americans (minor allele frequency, approximately 18%) but rare in European Americans (minor allele frequency <0.1%). Within strata of self-reported ancestry, this variant was associated with times of overall survival and disease-free survival in only African Americans (hazard ratio for overall survival, 2.82; 95% CI, 1.88-4.23; P = 5.0x10-07 and hazard ratio for disease-free survival, 2.27; 95% CI, 1.62-3.18; P = 1.8x10-06).

CONCLUSIONS

In an analysis of data from 2 trials of patients with stage II or III colon cancer, we identified rs76766811 as a potential prognostic variant in African American patients. This finding should be confirmed in additional study populations. ClinicalTrials.gov Identifiers: NCT00096278 (NSABP C-08) and NCT00079274 (NCCTG N0147).

Aberrant super-enhancer landscape reveals core transcriptional regulatory circuitry in lung adenocarcinoma

Lung adenocarcinoma (LUAD) relies on dysregulated gene expression to sustain its infinite growth and progression. Emerging evidence indicates that aberrant transcriptional program results from core transcriptional regulatory circuitry (CRC) which is driven by super-enhancers (SEs). In this study, by integrating profiles of H3K27Ac chromatin immunoprecipitation sequencing (ChIP-seq) from normal adult lung and LUAD cell lines, we revealed that widespread alterations of the super-enhancer were presence during lung carcinogenesis. With SE-based modeling of regulatory circuits and assessments of transcription factor (TF) dependencies, we reconstructed an interconnected transcriptional regulation network formed by three master TFs, including ELF3, EHF, and TGIF1, all of which promoted each other’s expression that confirmed by ChIP-qPCR and western blot. Loss-of function assay revealed that each of them is essential for LUAD cells survival, invasion and metastasis. Meanwhile, the rescue assay also illustrated the transacting transcriptional regulatory circuitry. In addition, the mRNA levels of ELF3, EHF, and TGIF1 were differentially expressed in LUAD tumors and peritumoral tissue. IHC of serial sections revealed that high expressions of CRC (ELF3/EHF/TGIF1-High) were closely associated with high proliferative activity in tumor tissue and poor prognosis on patients with LUAD. Finally, we used small molecular inhibitors to perturb the transcriptional circuitry, also exhibited a prominent anti-cancer effect in vitro. Our findings reveal the mechanism of the transcriptional dysregulation and addiction of LUAD.