Comprehensive analyses of partially methylated domains and differentially methylated regions in esophageal cancer reveal both cell-type- and cancer-specific epigenetic regulation

BACKGROUND

As one of the most common malignancies, esophageal cancer has two subtypes, squamous cell carcinoma and adenocarcinoma, arising from distinct cells-of-origin. Distinguishing cell-type-specific molecular features from cancer-specific characteristics is challenging.

RESULTS

We analyze whole-genome bisulfite sequencing data on 45 esophageal tumor and nonmalignant samples from both subtypes. We develop a novel sequence-aware method to identify large partially methylated domains (PMDs), revealing profound heterogeneity at both methylation level and genomic distribution of PMDs across tumor samples. We identify subtype-specific PMDs that are associated with repressive transcription, chromatin B compartments and high somatic mutation rate. While genomic locations of these PMDs are pre-established in normal cells, the degree of loss is significantly higher in tumors. We find that cell-type-specific deposition of H3K36me2 may underlie genomic distribution of PMDs. At a smaller genomic scale, both cell-type- and cancer-specific differentially methylated regions (DMRs) are identified for each subtype. Using binding motif analysis within these DMRs, we show that a cell-type-specific transcription factor HNF4A maintains the binding sites that it generates in normal cells, while establishing new binding sites cooperatively with novel partners such as FOSL1 in esophageal adenocarcinoma. Finally, leveraging pan-tissue single-cell and pan-cancer epigenomic datasets, we demonstrate that a substantial fraction of cell-type-specific PMDs and DMRs identified here in esophageal cancer are actually markers that co-occur in other cancers originating from related cell types.

CONCLUSIONS

These findings advance our understanding of DNA methylation dynamics at various genomic scales in normal and malignant states, providing novel mechanistic insights into cell-type- and cancer-specific epigenetic regulations.

DNA methylation and transcriptomic features are preserved throughout disease recurrence and chemoresistance in high grade serous ovarian cancers

Background

Little is known about the role of global DNA methylation in recurrence and chemoresistance of high grade serous ovarian cancer (HGSOC).

Methods

We performed whole genome bisulfite sequencing and transcriptome sequencing in 62 primary and recurrent tumors from 28 patients with stage III/IV HGSOC, of which 11 patients carried germline, pathogenic BRCA1 and/or BRCA2 mutations.

Results

Landscapes of genome-wide methylation (on average 24.2 million CpGs per tumor) and transcriptomes in primary and recurrent tumors showed extensive heterogeneity between patients but were highly preserved in tumors from the same patient. We identified significant differences in the burden of differentially methylated regions (DMRs) in tumors from BRCA1/2 compared to non-BRCA1/2 carriers (mean 659 DMRs and 388 DMRs in paired comparisons respectively). We identified overexpression of immune pathways in BRCA1/2 carriers compared to non-carriers, implicating an increased immune response in improved survival (P = 0.006) in these BRCA1/2 carriers.

Conclusion

These findings indicate methylome and gene expression programs established in the primary tumor are conserved throughout disease progression, even after extensive chemotherapy treatment, and that changes in methylation and gene expression are unlikely to serve as drivers for chemoresistance in HGSOC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02440-z.

ODP055 The methylation pattern of a unique Bilateral para-overian Adrenal Rest Tumor in a girl with Nicotinamide Nucleotide Transhydrogenase mutation

Background

Patients with NNT (Nicotinamide Nucleotide Transhydrogenase) gene mutations, a rare cause of glucocorticoid and mineralocorticoid deficiency require hormone replacement therapy. Adrenal Rest Tumor (ART) in females, reported so far only in noncompliant patients with congenital adrenal hyperplasia and elevated ACTH levels, is very rare (<20 cases world-wide). This study characterizes the pathophysiology, the molecular ontogeny and methylation analysis of a unique ART in a female with adrenal failure due to the G200S mutation in NNT Clinical presentation and Method: A 15-year-old girl, with homozygous G200S NNT-mutation followed for adrenal insufficiency reappeared to follow-up with severe virilization and elevated serum testosterone (28.3 nmol/l) and ACTH (> 1500 pmol/l). Pelvic MRI and Ultrasound demonstrated one sided paraovarian round tumor with pathological vascularization. Laparoscopic exploration revealed bilateral para ovarian mesosalpinx masses involving the serosa of the Fallopian tube, (3 and 1 cm in diameter). The testosterone level normalized within one day after surgical removal of those masses (0.2 nmol/l).

Results

Histopathology demonstrated a pattern of adrenal rest tissue with strong intracellular positive staining for adrenal markers such as SF-1, calretinin, MART1, inhibin and the pituitary marker ACTH. The staining for ovarian characteristic markers such as PAX 8 was negative. Studying mRNA extracted from the tissue by RT-PCR revealed the positive Gene expression of Cyp17a1, Cyp21a2 and Mc2r cDNA but not Pomc suggesting adrenal but not pituitary origin of the tissue. We further profiled the epigenomic profile of several adrenal rest tumors from both ovarian and testes origins using the Infinium Methylation EPIC array. We characterized the adrenal-specific features by comparing the tumors from the two different originating sites to published methylation array data on healthy adrenal tissue. We further investigated cancer-specific methylation changes to identify activated cancer pathways, and used the methylation arrays to identify somatic copy number alterations. Finally, we analyzed the likely developmental origin of these tumors by comparing to published methylation array data of developmental and adult reproductive tissues.

Conclusion

This study exemplifies severe virilization that resulted from a unique and rare type of ART in ovarian related tissue that was caused by incompliance to treatment in a patient with NNT gene mutation. The laparoscopic surgical findings indicate that imaging techniques may be insufficient in identification of such rest tumors and call for laparoscopy when clinical findings are suggestive. An early detection of this tumor could preserve fertility. Using histopathology markers cDNA studies and epigenomic profiling by methylation studies, our study shows for the first time that female ART originates from adrenal cells. The growth of a functional androgen producing "tumor" indicates that functional NNT protein is NOT required for androgen synthesis in contrast to glucocorticoids and that a zona reticularis similar tissue in ART is responsive to ACTH stimulation.

Presentation: No date and time listed

Detecting cell-of-origin and cancer-specific methylation features of cell-free DNA from Nanopore sequencing

The Oxford Nanopore (ONT) platform provides portable and rapid genome sequencing, and its ability to natively profile DNA methylation without complex sample processing is attractive for point-of-care real-time sequencing. We recently demonstrated ONT shallow whole-genome sequencing to detect copy number alterations (CNAs) from the circulating tumor DNA (ctDNA) of cancer patients. Here, we show that cell type and cancer-specific methylation changes can also be detected, as well as cancer-associated fragmentation signatures. This feasibility study suggests that ONT shallow WGS could be a powerful tool for liquid biopsy.

Using nucleosome levels, copy number alterations, and DNA methylation to profile disease state through liquid biopsy

e15014

Background: Nucleosomes are the repeating unit of chromatin that contain important signals for proper genomic function and transcriptional activity. Upon cell death, chromatin is fragmented, and nucleosomes are released into circulation and thus, detectable in liquid biopsy. Cell death naturally occurs over time and increases in a variety of diseases resulting in elevated levels of circulating nucleosomes and the genetic and epigenetic signatures they carry. These signatures contain important information about the cells from which they were derived, as well as the reason for them being in circulation. Methods: We measured circulating nucleosome levels, using Volition’s Nu.Q H3.1 sandwich ELISA across a healthy cohort spanning 40-85 years of age as well as patients whom had been diagnosed with Non-Hodgkin's Lymphoma (NHL), some of which were untreated and others that had undergone treatment. Additionally, to demonstrate that the correlation between nucleosome (Nu.Q) levels and copy number alterations is specific to cancer, we assessed the genomic integrity of an NHL sample with a Nu.Q level of 744 ng/ml and a sepsis sample with Nu.Q level of 3,281 ng/ml using low coverage Oxford Nanopore whole genome sequencing to identify genomic amplifications and deletions. We further used cell type specific methylation patterns to identify the cell of origin of cfDNA in each sample. Results: Importantly, we found no significant correlation between nucleosome level and age in our healthy cohort (R2= 0.0004; 40-49: N = 10, mean = 39 ng/ml; 50-59: N = 10, mean = 18 ng/ml; 60-69: N = 10, mean = 35 ng/ml; 70-79: N = 15, mean = 34 ng/ml; 80-85: N = 5, mean = 32 ng/ml). We did however find elevated nucleosome levels in untreated NHL patients (N = 4, mean 646 ng/ml) and that nucleosome levels in NHL patients that had undergone treatment were reduced (N = 29, mean = 180 ng/ml). Using copy number changes, we found that the fraction of cfDNA estimated to be derived from the tumor in the NHL sample was 54% and were able to detect amplifications in chromosomes 10, 11 and 18 and deletions in chromosomes 1, 8, 17, 18 and 23, whereas we found no copy number alterations present in the cfDNA sample derived from a sepsis patient and a tumor fraction of 0. Furthermore, using cell type specific methylation patterns to deconvolute the data we found that a large fraction of the cfDNA in the NHL sample was derived from B-cells consistent with the NHL sample being of B-cell origin whereas we determined that the cell of origin of the cfDNA in the sepsis sample was Neutrophils. Conclusions: This data shows that nucleosome levels, as measured by Nu.Q, can be used in conjunction with shallow sequencing and DNA methylation profiles to define tumor fraction and cell of origin in liquid biopsy.

Activation of bivalent factor DLX5 cooperates with master regulator TP63 to promote squamous cell carcinoma

To reconstruct systematically hyperactive transcription factor (TF)-dependent transcription networks in squamous cell carcinomas (SCCs), a computational method (ELMER) was applied to 1293 pan-SCC patient samples, and 44 hyperactive SCC TFs were identified. As a top candidate, DLX5 exhibits a notable bifurcate re-configuration of its bivalent promoter in cancer. Specifically, DLX5 maintains a bivalent state in normal tissues; its promoter is hypermethylation, leading to DLX5 transcriptional silencing in esophageal adenocarcinoma (EAC). In stark contrast, DLX5 promoter gains active histone marks and becomes transcriptionally activated in ESCC, which is directly mediated by SOX2. Functionally, silencing of DLX5 substantially inhibits SCC viability both in vitro and in vivo. Mechanistically, DLX5 cooperates with TP63 in regulating ∼2000 enhancers and promoters, which converge on activating cancer-promoting pathways. Together, our data establish a novel and strong SCC-promoting factor and elucidate a new epigenomic mechanism - bifurcate chromatin re-configuration - during cancer development.

Expanding the MECP2 network using comparative genomics reveals potential therapeutic targets for Rett syndrome

Inactivating mutations in the Methyl-CpG Binding Protein 2 (MECP2) gene are the main cause of Rett syndrome (RTT). Despite extensive research into MECP2 function, no treatments for RTT are currently available. Here, we used an evolutionary genomics approach to construct an unbiased MECP2 gene network, using 1028 eukaryotic genomes to prioritize proteins with strong co-evolutionary signatures with MECP2. Focusing on proteins targeted by FDA-approved drugs led to three promising targets, two of which were previously linked to MECP2 function (IRAK, KEAP1) and one that was not (EPOR). The drugs targeting these three proteins (Pacritinib, DMF, and EPO) were able to rescue different phenotypes of MECP2 inactivation in cultured human neural cell types, and appeared to converge on Nuclear Factor Kappa B (NF-κB) signaling in inflammation. This study highlights the potential of comparative genomics to accelerate drug discovery, and yields potential new avenues for the treatment of RTT.

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.

Whole-genome characterization of lung adenocarcinomas lacking the RTK/RAS/RAF pathway

RTK/RAS/RAF pathway alterations (RPAs) are a hallmark of lung adenocarcinoma (LUAD). In this study, we use whole-genome sequencing (WGS) of 85 cases found to be RPA(-) by previous studies from The Cancer Genome Atlas (TCGA) to characterize the minority of LUADs lacking apparent alterations in this pathway. We show that WGS analysis uncovers RPA(+) in 28 (33%) of the 85 samples. Among the remaining 57 cases, we observe focal deletions targeting the promoter or transcription start site of STK11 (n = 7) or KEAP1 (n = 3), and promoter mutations associated with the increased expression of ILF2 (n = 6). We also identify complex structural variations associated with high-level copy number amplifications. Moreover, an enrichment of focal deletions is found in TP53 mutant cases. Our results indicate that RPA(-) cases demonstrate tumor suppressor deletions and genome instability, but lack unique or recurrent genetic lesions compensating for the lack of RPAs. Larger WGS studies of RPA(-) cases are required to understand this important LUAD subset.

Ovarian Cancer Risk Variants Are Enriched in Histotype-Specific Enhancers and Disrupt Transcription Factor Binding Sites

Quantifying the functional effects of complex disease risk variants can provide insights into mechanisms underlying disease biology. Genome-wide association studies have identified 39 regions associated with risk of epithelial ovarian cancer (EOC). The vast majority of these variants lie in the non-coding genome, where they likely function through interaction with gene regulatory elements. In this study we first estimated the heritability explained by known common low penetrance risk alleles for EOC. The narrow sense heritability (hg2) of EOC overall and high-grade serous ovarian cancer (HGSOCs) were estimated to be 5%-6%. Partitioned SNP heritability across broad functional categories indicated a significant contribution of regulatory elements to EOC heritability. We collated epigenomic profiling data for 77 cell and tissue types from Roadmap Epigenomics and ENCODE, and from H3K27Ac ChIP-seq data generated in 26 ovarian cancer and precursor-related cell and tissue types. We identified significant enrichment of risk single-nucleotide polymorphisms (SNPs) in active regulatory elements marked by H3K27Ac in HGSOCs. To further investigate how risk SNPs in active regulatory elements influence predisposition to ovarian cancer, we used motifbreakR to predict the disruption of transcription factor binding sites. We identified 469 candidate causal risk variants in H3K27Ac peaks that are predicted to significantly break transcription factor (TF) motifs. The most frequently broken motif was REST (p value = 0.0028), which has been reported as both a tumor suppressor and an oncogene. Overall, these systematic functional annotations with epigenomic data improve interpretation of EOC risk variants and shed light on likely cells of origin.

Before and After: Comparison of Legacy and Harmonized TCGA Genomic Data Commons' Data

We present a systematic analysis of the effects of synchronizing a large-scale, deeply characterized, multi-omic dataset to the current human reference genome, using updated software, pipelines, and annotations. For each of 5 molecular data platforms in The Cancer Genome Atlas (TCGA)-mRNA and miRNA expression, single nucleotide variants, DNA methylation and copy number alterations-comprehensive sample, gene, and probe-level studies were performed, towards quantifying the degree of similarity between the 'legacy' GRCh37 (hg19) TCGA data and its GRCh38 (hg38) version as 'harmonized' by the Genomic Data Commons. We offer gene lists to elucidate differences that remained after controlling for confounders, and strategies to mitigate their impact on biological interpretation. Our results demonstrate that the hg19 and hg38 TCGA datasets are very highly concordant, promote informed use of either legacy or harmonized omics data, and provide a rubric that encourages similar comparisons as new data emerge and reference data evolve.

ELMER v.2: an R/Bioconductor package to reconstruct gene regulatory networks from DNA methylation and transcriptome profiles

Motivation

DNA methylation has been used to identify functional changes at transcriptional enhancers and other cis-regulatory modules (CRMs) in tumors and other disease tissues. Our R/Bioconductor package ELMER (Enhancer Linking by Methylation/Expression Relationships) provides a systematic approach that reconstructs altered gene regulatory networks (GRNs) by combining enhancer methylation and gene expression data derived from the same sample set.

Results

We present a completely revised version 2 of ELMER that provides numerous new features including an optional web-based interface and a new Supervised Analysis mode to use pre-defined sample groupings. We show that Supervised mode significantly increases statistical power and identifies additional GRNs and associated Master Regulators, such as SOX11 and KLF5 in Basal-like breast cancer.

Availability and implementation

ELMER v.2 is available as an R/Bioconductor package at http://bioconductor.org/packages/ELMER/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Non-coding somatic mutations converge on the PAX8 pathway in ovarian cancer

The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineage-specific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

Lineage-Specific Epigenomic and Genomic Activation of Oncogene HNF4A Promotes Gastrointestinal Adenocarcinomas

Gastrointestinal adenocarcinomas (GIAC) of the tubular gastrointestinal (GI) tract including esophagus, stomach, colon, and rectum comprise most GI cancers and share a spectrum of genomic features. However, the unified epigenomic changes specific to GIAC are poorly characterized. Using 907 GIAC samples from The Cancer Genome Atlas, we applied mathematical algorithms to large-scale DNA methylome and transcriptome profiles to reconstruct transcription factor (TF) networks and identify a list of functionally hyperactive master regulator (MR) TF shared across different GIAC. The top candidate HNF4A exhibited prominent genomic and epigenomic activation in a GIAC-specific manner. A complex interplay between the HNF4A promoter and three distal enhancer elements was coordinated by GIAC-specific MRTF including ELF3, GATA4, GATA6, and KLF5. HNF4A also self-regulated its own promoter and enhancers. Functionally, HNF4A promoted cancer proliferation and survival by transcriptional activation of many downstream targets, including HNF1A and factors of interleukin signaling, in a lineage-specific manner. Overall, our study provides new insights into the GIAC-specific gene regulatory networks and identifies potential therapeutic strategies against these common cancers. SIGNIFICANCE: These findings show that GIAC-specific master regulatory transcription factors control HNF4A via three distal enhancers to promote GIAC cell proliferation and survival. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/13/2722/F1.large.jpg.

Master transcription factors form interconnected circuitry and orchestrate transcriptional networks in oesophageal adenocarcinoma

OBJECTIVE

While oesophageal squamous cell carcinoma remains infrequent in Western populations, the incidence of oesophageal adenocarcinoma (EAC) has increased sixfold to eightfold over the past four decades. We aimed to characterise oesophageal cancer-specific and subtypes-specific gene regulation patterns and their upstream transcription factors (TFs).  DESIGN: To identify regulatory elements, we profiled fresh-frozen oesophageal normal samples, tumours and cell lines with chromatin immunoprecipitation sequencing (ChIP-Seq). Mathematical modelling was performed to establish (super)-enhancers landscapes and interconnected transcriptional circuitry formed by master TFs. Coregulation and cooperation between master TFs were investigated by ChIP-Seq, circularised chromosome conformation capture sequencing and luciferase assay. Biological functions of candidate factors were evaluated both in vitro and in vivo.

RESULTS

We found widespread and pervasive alterations of the (super)-enhancer reservoir in both subtypes of oesophageal cancer, leading to transcriptional activation of a myriad of novel oncogenes and signalling pathways, some of which may be exploited pharmacologically (eg, leukemia inhibitory factor (LIF) pathway). Focusing on EAC, we bioinformatically reconstructed and functionally validated an interconnected circuitry formed by four master TFs-ELF3, KLF5, GATA6 and EHF-which promoted each other's expression by interacting with each super-enhancer. Downstream, these master TFs occupied almost all EAC super-enhancers and cooperatively orchestrated EAC transcriptome. Each TF within the transcriptional circuitry was highly and specifically expressed in EAC and functionally promoted EAC cell proliferation and survival.

CONCLUSIONS

By establishing cancer-specific and subtype-specific features of the EAC epigenome, our findings promise to transform understanding of the transcriptional dysregulation and addiction of EAC, while providing molecular clues to develop novel therapeutic modalities against this malignancy.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Abstract AP11: IDENTIFYING FUNCTIONAL NONCODING SOMATIC MUTATIONS IN OVARIAN CANCER

Epithelial ovarian cancer (OC) is a heterogeneous disease that stratifies into different histologic subtypes including high-grade serous (HGSOC), clear cell (CCOC), endometrioid (EnOC) and mucinous (MOC) ovarian cancer. HGSOC is the most common, and each histotype is characterized by largely distinct germline genetics, somatic alterations and clinical biomarkers. Recently, whole genome sequencing (WGS) studies have catalogued genome-wide somatic variation for most OC histotypes. These data demonstrate that OC histotypes harbor thousands of noncoding somatic mutations and our next major challenge is to distinguish the few important noncoding somatic drivers from the thousands of passenger mutations.

We hypothesized that driver noncoding somatic mutations impact disease development and progression through altering the sequence of regulatory elements (REs), such as enhancers and promoters, resulting in perturbed expression of genes involved in neoplastic development of OC. To systematically address this hypothesis, we established genome-wide H3K27ac epigenomic profiles, annotating active REs for the different ovarian cancer histotypes using chromatin immunoprecipitation sequencing (ChIP-seq) in 20 fresh frozen primary OC tissue samples—five tumors for each major histotype. In parallel, we performed transcriptional profiling using RNA sequencing (RNA-seq). Together, these two datasets enabled us to evaluate epigenetic alterations and the transcriptome.

We identified histotype-specific active REs, and common active REs across all histotypes. First, we used the RNA-seq data to assess the effect of the histotype-specific REs in gene expression, and to find target genes of cis-REs and novel histotype-specific biomarkers. Next, we integrated these unique profiles with WGS data from 232 OCs (169 HGSOCs, 35 CCOCs and 28 EnOCs). The number of somatic single nucleotide variants per sample range from 481 to 40,764 (mean=7199, sd = 5751). Of these, 9.4% were noncoding and overlapped active REs in OC. Using a Poisson binomial distribution, we tested the significance of the observed number of mutated samples for any given active RE. We identified several frequently mutated REs in HGSOC, CCOC and EnOC, including promoters (POLR3E and CHCHD6 for HGSOC, CA5BP1 and C16orf87 for EnOC, DKC1 and SLIT3 for CCOC+EnOC), enhancers and super enhancers (KLF6 for HGSOC). We developed a novel approach to identifying target genes associated with a collection of frequently mutated enhancers by correlating enhancer activity (H3K27ac ChIP-Seq signal) vs. gene expression (RNA-seq signal), identifying HOXD cluster genes and C19orf44 as targets of mutated REs in OC.

In conclusion, we have used a multi-omics method to identify functional, driver noncoding somatic mutations for ovarian cancer based on their interaction with disease-specific regulatory elements. This represents a powerful way to distinguish important noncoding somatic drivers from a much larger number of passenger mutations that accumulate during tumor development.

Citation Format: Rosario I. Corona, Ji-Heui Seo, Dennis J. Hazelett, Xianzhi Lin, Paulette Y. Mhawech-Fauceglia, Jenny Lester, Sohrab P. Shah, David G. Huntsman, Beth Y. Karlan, Alexander Gusev, Benjamin P. Berman, Matthew L. Freedman, Simon A. Gayther, Kate Lawrenson. IDENTIFYING FUNCTIONAL NONCODING SOMATIC MUTATIONS IN OVARIAN CANCER [abstract]. In: Proceedings of the 12th Biennial Ovarian Cancer Research Symposium; Sep 13-15, 2018; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2019;25(22 Suppl):Abstract nr AP11.

GENAVi: a shiny web application for gene expression normalization, analysis and visualization

BACKGROUND

The development of next generation sequencing (NGS) methods led to a rapid rise in the generation of large genomic datasets, but the development of user-friendly tools to analyze and visualize these datasets has not developed at the same pace. This presents a two-fold challenge to biologists; the expertise to select an appropriate data analysis pipeline, and the need for bioinformatics or programming skills to apply this pipeline. The development of graphical user interface (GUI) applications hosted on web-based servers such as Shiny can make complex workflows accessible across operating systems and internet browsers to those without programming knowledge.

RESULTS

We have developed GENAVi (Gene Expression Normalization Analysis and Visualization) to provide a user-friendly interface for normalization and differential expression analysis (DEA) of human or mouse feature count level RNA-Seq data. GENAVi is a GUI based tool that combines Bioconductor packages in a format for scientists without bioinformatics expertise. We provide a panel of 20 cell lines commonly used for the study of breast and ovarian cancer within GENAVi as a foundation for users to bring their own data to the application. Users can visualize expression across samples, cluster samples based on gene expression or correlation, calculate and plot the results of principal components analysis, perform DEA and gene set enrichment and produce plots for each of these analyses. To allow scalability for large datasets we have provided local install via three methods. We improve on available tools by offering a range of normalization methods and a simple to use interface that provides clear and complete session reporting and for reproducible analysis.

CONCLUSION

The development of tools using a GUI makes them practical and accessible to scientists without bioinformatics expertise, or access to a data analyst with relevant skills. While several GUI based tools are currently available for RNA-Seq analysis we improve on these existing tools. This user-friendly application provides a convenient platform for the normalization, analysis and visualization of gene expression data for scientists without bioinformatics expertise.

RNA-Binding Protein ZFP36L1 Suppresses Hypoxia and Cell-Cycle Signaling

ZFP36L1 is a tandem zinc-finger RNA-binding protein that recognizes conserved adenylate-uridylate-rich elements (ARE) located in 3'untranslated regions (UTR) to mediate mRNA decay. We hypothesized that ZFP36L1 is a negative regulator of a posttranscriptional hub involved in mRNA half-life regulation of cancer-related transcripts. Analysis of in silico data revealed that ZFP36L1 was significantly mutated, epigenetically silenced, and downregulated in a variety of cancers. Forced expression of ZFP36L1 in cancer cells markedly reduced cell proliferation in vitro and in vivo, whereas silencing of ZFP36L1 enhanced tumor cell growth. To identify direct downstream targets of ZFP36L1, systematic screening using RNA pull-down of wild-type and mutant ZFP36L1 as well as whole transcriptome sequencing of bladder cancer cells {plus minus} tet-on ZFP36L1 was performed. A network of 1,410 genes was identified as potential direct targets of ZFP36L1. These targets included a number of key oncogenic transcripts such as HIF1A, CCND1, and E2F1. ZFP36L1 specifically bound to the 3'UTRs of these targets for mRNA degradation, thus suppressing their expression. Dual luciferase reporter assays and RNA electrophoretic mobility shift assays showed that wild-type, but not zinc-finger mutant ZFP36L1, bound to HIF1A 3'UTR and mediated HIF1A mRNA degradation, leading to reduced expression of HIF1A and its downstream targets. Collectively, our findings reveal an indispensable role of ZFP36L1 as a posttranscriptional safeguard against aberrant hypoxic signaling and abnormal cell-cycle progression. SIGNIFICANCE: RNA-binding protein ZFP36L1 functions as a tumor suppressor by regulating the mRNA stability of a number of mRNAs involved in hypoxia and cell-cycle signaling.

Super-enhancer-associated MEIS1 promotes transcriptional dysregulation in Ewing sarcoma in co-operation with EWS-FLI1

As the second most common malignant bone tumor in children and adolescents, Ewing sarcoma is initiated and exacerbated by a chimeric oncoprotein, most commonly, EWS-FLI1. In this study, we apply epigenomic analysis to characterize the transcription dysregulation in this cancer, focusing on the investigation of super-enhancer and its associated transcriptional regulatory mechanisms. We demonstrate that super-enhancer-associated transcripts are significantly enriched in EWS-FLI1 target genes, contribute to the aberrant transcriptional network of the disease, and mediate the exceptional sensitivity of Ewing sarcoma to transcriptional inhibition. Through integrative analysis, we identify MEIS1 as a super-enhancer-driven oncogene, which co-operates with EWS-FLI1 in transcriptional regulation, and plays a key pro-survival role in Ewing sarcoma. Moreover, APCDD1, another super-enhancer-associated gene, acting as a downstream target of both MEIS1 and EWS-FLI1, is also characterized as a novel tumor-promoting factor in this malignancy. These data delineate super-enhancer-mediated transcriptional deregulation in Ewing sarcoma, and uncover numerous candidate oncogenes which can be exploited for further understanding of the molecular pathogenesis for this disease.

Mapping Citrullinated Sites in Multiple Organs of Mice Using Hypercitrullinated Library

Protein citrullination (or deimination), an irreversible post-translational modification, has been implicated in several physiological and pathological processes, including gene expression regulation, apoptosis, rheumatoid arthritis, and Alzheimer's disease. Several research studies have been carried out on citrullination under many conditions. However, until now, challenges in sample preparation and data analysis have made it difficult to confidently identify a citrullinated protein and assign the citrullinated site. To overcome these limitations, we generated a mouse hyper-citrullinated spectral library and set up coordinates to confidently identify and validate citrullinated sites. Using this workflow, we detect a four-fold increase in citrullinated proteome coverage across six mouse organs compared with the current state-of-the art techniques. Our data reveal that the subcellular distribution of citrullinated proteins is tissue-type-dependent and that citrullinated targets are involved in fundamental physiological processes, including the metabolic process. These data represent the first report of a hyper-citrullinated library for the mouse and serve as a central resource for exploring the role of citrullination in this organism.

The chromatin accessibility landscape of primary human cancers

We present the genome-wide chromatin accessibility profiles of 410 tumor samples spanning 23 cancer types from The Cancer Genome Atlas (TCGA). We identify 562,709 transposase-accessible DNA elements that substantially extend the compendium of known cis-regulatory elements. Integration of ATAC-seq (the assay for transposase-accessible chromatin using sequencing) with TCGA multi-omic data identifies a large number of putative distal enhancers that distinguish molecular subtypes of cancers, uncovers specific driving transcription factors via protein-DNA footprints, and nominates long-range gene-regulatory interactions in cancer. These data reveal genetic risk loci of cancer predisposition as active DNA regulatory elements in cancer, identify gene-regulatory interactions underlying cancer immune evasion, and pinpoint noncoding mutations that drive enhancer activation and may affect patient survival. These results suggest a systematic approach to understanding the noncoding genome in cancer to advance diagnosis and therapy.

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

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

Identification of distinct mutational patterns and new driver genes in oesophageal squamous cell carcinomas and adenocarcinomas

OBJECTIVES

Oesophageal squamous cell carcinoma (OSCC) and adenocarcinoma (OAC) are distinct cancers in terms of a number of clinical and epidemiological characteristics, complicating the design of clinical trials and biomarker developments. We analysed 1048 oesophageal tumour-germline pairs from both subtypes, to characterise their genomic features, and biological and clinical significance.

DESIGN

Previously exome-sequenced samples were re-analysed to identify significantly mutated genes (SMGs) and mutational signatures. The biological functions of novel SMGs were investigated using cell line and xenograft models. We further performed whole-genome bisulfite sequencing and chromatin immunoprecipitation (ChIP)-seq to characterise epigenetic alterations.

RESULTS

OSCC and OAC displayed nearly mutually exclusive sets of driver genes, indicating that they follow independent developmental paths. The combined sample size allowed the statistical identification of a number of novel subtype-specific SMGs, mutational signatures and prognostic biomarkers. Particularly, we identified a novel mutational signature similar to Catalogue Of Somatic Mutations In Cancer (COSMIC)signature 16, which has prognostic value in OSCC. Two newly discovered SMGs, CUL3 and ZFP36L2, were validated as important tumour-suppressors specific to the OSCC subtype. We further identified their additional loss-of-function mechanisms. CUL3 was homozygously deleted specifically in OSCC and other squamous cell cancers (SCCs). Notably, ZFP36L2 is associated with super-enhancer in healthy oesophageal mucosa; DNA hypermethylation in its super-enhancer reduced active histone markers in squamous cancer cells, suggesting an epigenetic inactivation of a super-enhancer-associated SCC suppressor.

CONCLUSIONS

These data comprehensively contrast differences between OSCC and OAC at both genomic and epigenomic levels, and reveal novel molecular features for further delineating the pathophysiological mechanisms and treatment strategies for these cancers.

Coronary artery calcium quantification using contrast-enhanced dual-energy computed tomography scans in comparison with unenhanced single-energy scans

Extracting coronary artery calcium (CAC) scores from contrast-enhanced computed tomography (CT) images using dual-energy (DE) based material decomposition has been shown feasible, mainly through patient studies. However, the quantitative performance of such DE-based CAC scores, particularly per stenosis, is underexamined due to lack of reference standard and repeated scans. In this work we conducted a comprehensive quantitative comparative analysis of CAC scores obtained with DE and compare to conventional unenhanced single-energy (SE) CT scans through phantom studies. Synthetic vessels filled with iodinated blood mimicking material and containing calcium stenoses of different sizes and densities were scanned with a third generation dual-source CT scanner in a chest phantom using a DE coronary CT angiography protocol with three exposures/CTDIvol: auto-mAs/8 mGy (automatic exposure), 160 mAs/20 mGy and 260 mAs/34 mGy and 10 repeats. As a control, a set of vessel phantoms without iodine was scanned using a standard SE CAC score protocol (3 mGy). Calcium volume, mass and Agatston scores were estimated for each stenosis. For DE dataset, image-based three-material decomposition was applied to remove iodine before scoring. Performance of DE-based calcium scores were analyzed on a per-stenosis level and compared to SE-based scores. There was excellent correlation between the DE- and SE-based scores (correlation coefficient r: 0.92-0.98). Percent bias for the calcium volume and mass scores varied as a function of stenosis size and density for both modalities. Precision (coefficient of variation) improved with larger and denser stenoses for both DE- and SE-based calcium scores. DE-based scores (20 mGy and 34 mGy) provided comparable per-stenosis precision to SE-based (3 mGy). Our findings suggest that on a per-stenosis level, DE-based CAC scores from contrast-enhanced CT images can achieve comparable quantification performance to conventional SE-based scores. However, DE-based CAC scoring required more dose compared with SE for high per-stenosis precision so some caution is necessary with clinical DE-based CAC scoring.

Abstract 395: Identifying the functional drivers of noncoding somatic mutations in ovarian cancer

Epithelial ovarian cancer (OC) is a heterogeneous disease that has been stratified into different histologic subtypes: high-grade serous (HGSOC), clear cell (CCOC), endometrioid (EnOC) and mucinous (MOC). HGSOC is the most common, but every histotype is characterized by largely distinct germline genetics, somatic alterations and clinical biomarkers. Recently, whole-genome sequencing (WGS) studies have catalogued genome-wide somatic variation for most OC histotypes. These data demonstrate that OC histotypes harbor thousands of noncoding somatic mutations and our next major challenge is to distinguish the few important noncoding somatic drivers from the thousands of passenger mutations. We hypothesized that driver noncoding somatic mutations impact disease development and progression through altering the sequence of regulatory elements (REs), such as enhancers and promoters, eventually resulting in perturbation of the expression of target genes. To systematically address this hypothesis, we established genome-wide H3K27ac epigenomic profiles, annotating active REs for the different ovarian cancer histotypes using chromatin immunoprecipitation sequencing (ChIP-seq) in 20 fresh frozen primary OC tissue samples—five tumors for each major histotype. In parallel, we performed transcriptional profiling using RNA sequencing (RNA-seq). Together, these two datasets enabled us to evaluate epigenetic alterations and the transcriptome. We identified histotype-specific active REs, and common active REs across all histotypes. We used the RNA-seq data to assess the effect of the histotype-specific REs in gene expression, and to find target genes of cis-REs and novel histotype-specific biomarkers. Next, we integrated these unique profiles with WGS data from 232 OCs (169 HGSOCs, 35 CCOCs and 28 EnOCs). The number of somatic single-nucleotide variants per sample range from 481 to 40,764 (mean=7199, sd = 5751). Of these, 9.4% were noncoding and overlapped active REs in OC. Using a Poisson binomial distribution, we tested the significance of the observed number of mutated samples for any given active RE. Importantly, the method adjusted for interpatient mutation rate heterogeneity. Using the Benjamini-Hochberg procedure for p-value correction, we identified several significantly mutated active REs, including the promoters of POLR3E and WDR74. In conclusion, we have used an integrative method to identify functional, driver noncoding somatic mutations for ovarian cancer based on their interaction with disease-specific regulatory elements, and their putative target genes. This represents a powerful way to distinguish important noncoding somatic drivers from a much larger number of passenger mutations that accumulate during tumor development. The approach can be applied to other cancer types if sufficient depth and quality of genetic, transcriptomic and epigenomic datasets are available.

Citation Format: Rosario I. Corona, Ji-Heui Seo, Dennis J. Hazelett, Xianzhi Lin, Paulette Y. Mhawech-Fauceglia, Jenny Lester, Sohrab Shah, David G. Huntsman, Beth Y. Karlan, Benjamin P. Berman, Matthew L. Freedman, Simon A. Gayther, Kate Lawrenson. Identifying the functional drivers of noncoding somatic mutations in ovarian cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 395.

TCGAbiolinksGUI: A graphical user interface to analyze cancer molecular and clinical data

The GDC (Genomic Data Commons) data portal provides users with data from cancer genomics studies. Recently, we developed the R/Bioconductor TCGAbiolinks package, which allows users to search, download and prepare cancer genomics data for integrative data analysis. The use of this package requires users to have advanced knowledge of R thus limiting the number of users. To overcome this obstacle and improve the accessibility of the package by a wider range of users, we developed a graphical user interface (GUI) using Shiny available through the package TCGAbiolinksGUI. The TCGAbiolinksGUI package is freely available within the Bioconductor project at http://bioconductor.org/packages/TCGAbiolinksGUI/. Links to the GitHub repository, a demo version of the tool, a docker image and PDF/video tutorials are available from the TCGAbiolinksGUI site.

DNA methylation loss in late-replicating domains is linked to mitotic cell division

DNA methylation loss occurs frequently in cancer genomes, primarily within lamina-associated, late-replicating regions termed Partially Methylated Domains (PMDs). We profiled 39 diverse primary tumors and 8 matched adjacent tissues using Whole-Genome Bisulfite Sequencing (WGBS), and analyzed them alongside 343 additional human and 206 mouse WGBS datasets. We identified a local CpG sequence context associated with preferential hypomethylation in PMDs. Analysis of CpGs in this context (“Solo-WCGWs”) revealed previously undetected PMD hypomethylation in almost all healthy tissue types. PMD hypomethylation increased with age, beginning during fetal development, and appeared to track the accumulation of cell divisions. In cancer, PMD hypomethylation depth correlated with somatic mutation density and cell-cycle gene expression, consistent with its reflection of mitotic history, and suggesting its application as a mitotic clock. We propose that late replication leads to lifelong progressive methylation loss, which acts as a biomarker for cellular aging and which may contribute to oncogenesis.

Alpha-oxoglutarate inhibits the proliferation of immortalized normal bladder epithelial cells via an epigenetic switch involving ARID1A

Interstitial cystitis (IC) is a chronic urinary tract disease that is characterized by unpleasant sensations, such as persistent pelvic pain, in the absence of infection or other identifiable causes. We previously performed comprehensive metabolomics profiling of urine samples from IC patients using nuclear magnetic resonance and gas-chromatography/mass spectrometry and found that urinary α-oxoglutarate (α-OG), was significantly elevated. α-OG, a tricarboxylic acid (TCA) cycle intermediate, reportedly functions to suppress the proliferation of immortalized normal human bladder epithelial cells. Here, we identified AT-rich interactive domain 1 A (ARID1A), a key chromatin remodeler, as being hypomethylated and upregulated by α-OG treatment. This was done through EPIC DNA methylation profiling and subsequent biochemical approaches, including quantitative RT-PCR and western blot analyses. Furthermore, we found that α-OG almost completely suppresses ten-eleven translocation (TET) activity, but does not affect DNA methyltransferase (DNMT) activity. Altogether, our studies reveal the potential role of α-OG in epigenetic remodeling through its effects on ARID1A and TET expression in the bladder. This may provide a new possible therapeutic strategy in treating IC.

Quantitative assessment of nonsolid pulmonary nodule volume with computed tomography in a phantom study

Background

To assess the volumetric measurement of small (≤1 cm) nonsolid nodules with computed tomography (CT), focusing on the interaction of state of the art iterative reconstruction (IR) methods and dose with nodule densities, sizes, and shapes.

Methods

Twelve synthetic nodules [5 and 10 mm in diameter, densities of -800, -630 and -10 Hounsfield units (HU), spherical and spiculated shapes] were scanned within an anthropomorphic phantom. Dose [computed tomography scan dose index (CTDI_vol)] ranged from standard (4.1 mGy) to below screening levels (0.3 mGy). Data was reconstructed using filtered back-projection and two state-of-the-art IR methods (adaptive and model-based). Measurements were extracted with a previously validated matched filter-based estimator. Analysis of accuracy and precision was based on evaluation of percent bias (PB) and the repeatability coefficient (RC) respectively.

Results

Density had the most important effect on measurement error followed by the interaction of density with nodule size. The nonsolid -630 HU nodules had high accuracy and precision at levels comparable to solid (-10 HU) nonsolid, regardless of reconstruction method and with CTDI_vol as low as 0.6 mGy. PB was <5% and <11% for the 10- and 5-mm in nominal diameter -630 HU nodules respectively, and RC was <5% and <12% for the same nodules. For nonsolid -800 HU nodules, PB increased to <11% and <30% for the 10- and 5-mm nodules respectively, whereas RC increased slightly overall but varied widely across dose and reconstruction algorithms for the 5-mm nodules. Model-based IR improved measurement accuracy for the 5-mm, low-density (-800, -630 HU) nodules. For other nodules the effect of reconstruction method was small. Dose did not affect volumetric accuracy and only affected slightly the precision of 5-mm nonsolid nodules.

Conclusions

Reasonable values of both accuracy and precision were achieved for volumetric measurements of all 10-mm nonsolid nodules, and for the 5-mm nodules with -630 HU or higher density, when derived from scans acquired with below screening dose levels as low as 0.6 mGy and regardless of reconstruction algorithm.

Exosomal MicroRNA Transfer Into Macrophages Mediates Cellular Postconditioning

BACKGROUND

Cardiosphere-derived cells (CDCs) confer cardioprotection in acute myocardial infarction by distinctive macrophage (Mϕ) polarization. Here we demonstrate that CDC-secreted exosomes (CDC_exo) recapitulate the cardioprotective effects of CDC therapy known as cellular postconditioning.

METHODS

Rats and pigs underwent myocardial infarction induced by ischemia/reperfusion before intracoronary infusion of CDC_exo, inert fibroblast exosomes (Fb_exo; control), or vehicle. Two days later, infarct size was quantified. Macrophages were isolated from cardiac tissue or bone marrow for downstream analyses. RNA sequencing was used to determine exosome content and alterations in gene expression profiles in Mϕ.

RESULTS

Administration of CDC_exo but not Fb_exo after reperfusion reduces infarct size in rat and pig models of myocardial infarction. Furthermore, CDC_exo reduce the number of CD68+ Mϕ within infarcted tissue and modify the polarization state of Mϕ so as to mimic that induced by CDCs. CDC_exo are enriched in several miRNAs (including miR-146a, miR-181b, and miR-126) relative to Fb_exo. Reverse pathway analysis of whole-transcriptome data from CDC_exo-primed Mϕ implicated miR-181b as a significant (P=1.3x10^-21) candidate mediator of CDC-induced Mϕ polarization, and PKCδ (protein kinase C δ) as a downstream target. Otherwise inert Fb_exo loaded selectively with miR-181b alter Mϕ phenotype and confer cardioprotective efficacy in a rat model of myocardial infarction. Adoptive transfer of PKCδ-suppressed Mϕ recapitulates cardioprotection.

CONCLUSIONS

Our data support the hypothesis that exosomal transfer of miR-181b from CDCs into Mϕ reduces PKCδ transcript levels and underlies the cardioprotective effects of CDCs administered after reperfusion.

Genomic and Epigenomic Heterogeneity of Hepatocellular Carcinoma

Understanding the intratumoral heterogeneity of hepatocellular carcinoma is instructive for developing personalized therapy and identifying molecular biomarkers. Here we applied whole-exome sequencing to 69 samples from 11 patients to resolve the genetic architecture of subclonal diversification. Spatial genomic diversity was found in all 11 hepatocellular carcinoma cases, with 29% of driver mutations being heterogeneous, including TERT, ARID1A, NOTCH2, and STAG2. Similar with other cancer types, TP53 mutations were always shared between all tumor regions, that is, located on the "trunk" of the evolutionary tree. In addition, we found that variants within several drug targets such as KIT, SYK, and PIK3CA were mutated in a fully clonal manner, indicating their therapeutic potentials for hepatocellular carcinoma. Temporal dissection of mutational signatures suggested that mutagenic processes associated with exposure to aristolochic acid and aflatoxin might play a more important role in early, as opposed to late, stages of hepatocellular carcinoma development. Moreover, we observed extensive intratumoral epigenetic heterogeneity in hepatocellular carcinoma based on multiple independent analytical methods and showed that intratumoral methylation heterogeneity might play important roles in the biology of hepatocellular carcinoma cells. Our results also demonstrated prominent heterogeneity of intratumoral methylation even in a stable hepatocellular carcinoma genome. Together, these findings highlight widespread intratumoral heterogeneity at both the genomic and epigenomic levels in hepatocellular carcinoma and provide an important molecular foundation for better understanding the pathogenesis of this malignancy. Cancer Res; 77(9); 2255-65. ©2017 AACR.

FOXC1 Activates Smoothened-Independent Hedgehog Signaling in Basal-like Breast Cancer

The mesoderm- and epithelial-mesenchymal transition-associated transcription factor FOXC1 is specifically overexpressed in basal-like breast cancer (BLBC), but its biochemical function is not understood. Here, we demonstrate that FOXC1 controls cancer stem cell (CSC) properties enriched in BLBC cells via activation of Smoothened (SMO)-independent Hedgehog (Hh) signaling. This non-canonical activation of Hh is specifically mediated by Gli2. Furthermore, we show that the N-terminal domain of FOXC1 (aa 1-68) binds directly to an internal region (aa 898-1168) of Gli2, enhancing the DNA-binding and transcription-activating capacity of Gli2. FOXC1 expression correlates with that of Gli2 and its targets in human breast cancers. Moreover, FOXC1 overexpression reduces sensitivity to anti-Hedgehog (Hh) inhibitors in BLBC cells and xenograft tumors. Together, these findings reveal FOXC1-mediated non-canonical Hh signaling that determines the BLBC stem-like phenotype and anti-Hh sensitivity, supporting inhibition of FOXC1 pathways as potential approaches for improving BLBC treatment.

Demystifying the secret mission of enhancers: linking distal regulatory elements to target genes

Enhancers are short regulatory sequences bound by sequence-specific transcription factors and play a major role in the spatiotemporal specificity of gene expression patterns in development and disease. While it is now possible to identify enhancer regions genomewide in both cultured cells and primary tissues using epigenomic approaches, it has been more challenging to develop methods to understand the function of individual enhancers because enhancers are located far from the gene(s) that they regulate. However, it is essential to identify target genes of enhancers not only so that we can understand the role of enhancers in disease but also because this information will assist in the development of future therapeutic options. After reviewing models of enhancer function, we discuss recent methods for identifying target genes of enhancers. First, we describe chromatin structure-based approaches for directly mapping interactions between enhancers and promoters. Second, we describe the use of correlation-based approaches to link enhancer state with the activity of nearby promoters and/or gene expression. Third, we describe how to test the function of specific enhancers experimentally by perturbing enhancer–target relationships using high-throughput reporter assays and genome editing. Finally, we conclude by discussing as yet unanswered questions concerning how enhancers function, how target genes can be identified, and how to distinguish direct from indirect changes in gene expression mediated by individual enhancers.

Abstract 2210: DNA methylation to identify non-coding gene regulatory regions in tumor samples

Altered control of transcriptional regulation via mutations to transcription factors and chromatin modifiers is a hallmark of cancer. Many popular methods for studying transcriptional regulation, such as ChIP-seq, DNase-I hypersensitivity, and Chromatin Conformation Capture, can be challenging to apply to clinical samples. DNA methylation has long been the chromatin biomarker of choice in disease studies because of its stability under clinical collection conditions. We have used DNA methylation profiles of thousands of tumors within the TCGA to reconstruct cancer-specific transcription factor regulatory networks, including non-coding transcription factor binding sites in enhancer regions. We use methyl-expression QTL (M-&gt;E QTL) analysis to associate changes at individual binding sites with their gene targets, often hundreds of thousands of base pairs away. This new approach allows for the first time the identification of in vivo enhancer changes in thousands of individuals, and provides a technology platform to identify key regulatory events in cancer progression.

Note: This abstract was not presented at the meeting.

Citation Format: Benjamin P. Berman. DNA methylation to identify non-coding gene regulatory regions in tumor samples. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2210. doi:10.1158/1538-7445.AM2015-2210

Gender-specific postnatal demethylation and establishment of epigenetic memory

DNA methylation patterns are set up in a relatively fixed programmed manner during normal embryonic development and are then stably maintained. Using genome-wide analysis, we discovered a postnatal pathway involving gender-specific demethylation that occurs exclusively in the male liver. This demodification is programmed to take place at tissue-specific enhancer sequences, and our data show that the methylation state at these loci is associated with and appears to play a role in the transcriptional regulation of nearby genes. This process is mediated by the secretion of testosterone at the time of sexual maturity, but the resulting methylation profile is stable and therefore can serve as an epigenetic memory even in the absence of this inducer. These findings add a new dimension to our understanding of the role of DNA methylation in vivo and provide the foundations for deciphering how environment can impact on the epigenetic regulation of genes in general.

The role of DNA methylation in directing the functional organization of the cancer epigenome

The holistic role of DNA methylation in the organization of the cancer epigenome is not well understood. Here we perform a comprehensive, high-resolution analysis of chromatin structure to compare the landscapes of HCT116 colon cancer cells and a DNA methylation-deficient derivative. The NOMe-seq accessibility assay unexpectedly revealed symmetrical and transcription-independent nucleosomal phasing across active, poised, and inactive genomic elements. DNA methylation abolished this phasing primarily at enhancers and CpG island (CGI) promoters, with little effect on insulators and non-CGI promoters. Abolishment of DNA methylation led to the context-specific reestablishment of the poised and active states of normal colon cells, which were marked in methylation-deficient cells by distinct H3K27 modifications and the presence of either well-phased nucleosomes or nucleosome-depleted regions, respectively. At higher-order genomic scales, we found that long, H3K9me3-marked domains had lower accessibility, consistent with a more compact chromatin structure. Taken together, our results demonstrate the nuanced and context-dependent role of DNA methylation in the functional, multiscale organization of cancer epigenomes.

Global loss of DNA methylation uncovers intronic enhancers in genes showing expression changes

BACKGROUND

Gene expression is epigenetically regulated by a combination of histone modifications and methylation of CpG dinucleotides in promoters. In normal cells, CpG-rich promoters are typically unmethylated, marked with histone modifications such as H3K4me3, and are highly active. During neoplastic transformation, CpG dinucleotides of CG-rich promoters become aberrantly methylated, corresponding with the removal of active histone modifications and transcriptional silencing. Outside of promoter regions, distal enhancers play a major role in the cell type-specific regulation of gene expression. Enhancers, which function by bringing activating complexes to promoters through chromosomal looping, are also modulated by a combination of DNA methylation and histone modifications.

RESULTS

Here we use HCT116 colorectal cancer cells with and without mutations in DNA methyltransferases, the latter of which results in a 95% reduction in global DNA methylation levels. These cells are used to study the relationship between DNA methylation, histone modifications, and gene expression. We find that the loss of DNA methylation is not sufficient to reactivate most of the silenced promoters. In contrast, the removal of DNA methylation results in the activation of a large number of enhancer regions as determined by the acquisition of active histone marks.

CONCLUSIONS

Although the transcriptome is largely unaffected by the loss of DNA methylation, we identify two distinct mechanisms resulting in the upregulation of distinct sets of genes. One is a direct result of DNA methylation loss at a set of promoter regions and the other is due to the presence of new intragenic enhancers.

Identification and characterization of functional risk variants for colorectal cancer mapping to chromosome 11q23.1

Genome-wide association studies of colorectal cancer (CRC) have identified a number of common variants associated with modest risk, including rs3802842 at chromosome 11q23.1. Several genes map to this region but rs3802842 does not map to any known transcribed or regulatory sequences. We reasoned, therefore, that rs3802842 is not the functional single-nucleotide polymorphism (SNP), but is in linkage disequilibrium (LD) with a functional SNP(s). We performed ChIP-seq for histone modifications in SW480 and HCT-116 CRC cells, and incorporated ChIP-seq and DNase I hypersensitivity data available through ENCODE within a 137-kb genomic region containing rs3802842 on 11q23.1. We identified SNP rs10891246 in LD with rs3802842 that mapped within a bidirectional promoter region of genes C11orf92 and C11orf93. Following mutagenesis to the risk allele, the promoter demonstrated lower levels of reporter gene expression. A second SNP rs7130173 was identified in LD with rs3802842 that mapped to a candidate enhancer region, which showed strong unidirectional activity in both HCT-116 and SW480 CRC cells. The risk allele of rs7130173 demonstrated reduced enhancer activity compared with the common allele, and reduced nuclear protein binding affinity in electromobility shift assays compared with the common allele suggesting differential transcription factor (TF) binding. SNPs rs10891246 and rs7130173 are on the same haplotype, and expression quantitative trait loci (eQTL) analyses of neighboring genes implicate C11orf53, C11orf92 and C11orf93 as candidate target genes. These data imply that rs10891246 and rs7130173 are functional SNPs mapping to 11q23.1 and that C11orf53, C11orf92 and C11orf93 represent novel candidate target genes involved in CRC etiology.

Integrated transcriptomic and epigenomic analysis of primary human lung epithelial cell differentiation

Elucidation of the epigenetic basis for cell-type specific gene regulation is key to gaining a full understanding of how the distinct phenotypes of differentiated cells are achieved and maintained. Here we examined how epigenetic changes are integrated with transcriptional activation to determine cell phenotype during differentiation. We performed epigenomic profiling in conjunction with transcriptomic profiling using in vitro differentiation of human primary alveolar epithelial cells (AEC). This model recapitulates an in vivo process in which AEC transition from one differentiated cell type to another during regeneration following lung injury. Interrogation of histone marks over time revealed enrichment of specific transcription factor binding motifs within regions of changing chromatin structure. Cross-referencing of these motifs with pathways showing transcriptional changes revealed known regulatory pathways of distal alveolar differentiation, such as the WNT and transforming growth factor beta (TGFB) pathways, and putative novel regulators of adult AEC differentiation including hepatocyte nuclear factor 4 alpha (HNF4A), and the retinoid X receptor (RXR) signaling pathways. Inhibition of the RXR pathway confirmed its functional relevance for alveolar differentiation. Our incorporation of epigenetic data allowed specific identification of transcription factors that are potential direct upstream regulators of the differentiation process, demonstrating the power of this approach. Integration of epigenomic data with transcriptomic profiling has broad application for the identification of regulatory pathways in other models of differentiation.

Understanding the role of epigenetic control of gene expression is critical to the full description of biological processes, such as development and regeneration. Herein we utilize the differentiation of cells from the distal lung to gain insight into the correlation between the epigenetic landscape, molecular signaling events, and eventual changes in transcription and phenotype. We found that by integrating epigenetic profiling with whole genome transcriptomic data we were able to determine which molecular signaling events were activated and repressed during adult alveolar epithelial cell differentiation, and we identified epigenetic changes that contributed to these changes. Furthermore, we validated the role of one of these predicted but not previously identified pathways, retinoid X receptor signaling, in this process.

FunciSNP: an R/bioconductor tool integrating functional non-coding data sets with genetic association studies to identify candidate regulatory SNPs

Single nucleotide polymorphisms (SNPs) are increasingly used to tag genetic loci associated with phenotypes such as risk of complex diseases. Technically, this is done genome-wide without prior restriction or knowledge of biological feasibility in scans referred to as genome-wide association studies (GWAS). Depending on the linkage disequilibrium (LD) structure at a particular locus, such tagSNPs may be surrogates for many thousands of other SNPs, and it is difficult to distinguish those that may play a functional role in the phenotype from those simply genetically linked. Because a large proportion of tagSNPs have been identified within non-coding regions of the genome, distinguishing functional from non-functional SNPs has been an even greater challenge. A strategy was recently proposed that prioritizes surrogate SNPs based on non-coding chromatin and epigenomic mapping techniques that have become feasible with the advent of massively parallel sequencing. Here, we introduce an R/Bioconductor software package that enables the identification of candidate functional SNPs by integrating information from tagSNP locations, lists of linked SNPs from the 1000 genomes project and locations of chromatin features which may have functional significance.

AVAILABILITY

FunciSNP is available from Bioconductor (bioconductor.org).

Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules

DNA methylation and nucleosome positioning work together to generate chromatin structures that regulate gene expression. Nucleosomes are typically mapped using nuclease digestion requiring significant amounts of material and varying enzyme concentrations. We have developed a method (NOMe-seq) that uses a GpC methyltransferase (M.CviPI) and next generation sequencing to generate a high resolution footprint of nucleosome positioning genome-wide using less than 1 million cells while retaining endogenous DNA methylation information from the same DNA strand. Using a novel bioinformatics pipeline, we show a striking anti-correlation between nucleosome occupancy and DNA methylation at CTCF regions that is not present at promoters. We further show that the extent of nucleosome depletion at promoters is directly correlated to expression level and can accommodate multiple nucleosomes and provide genome-wide evidence that expressed non-CpG island promoters are nucleosome-depleted. Importantly, NOMe-seq obtains DNA methylation and nucleosome positioning information from the same DNA molecule, giving the first genome-wide DNA methylation and nucleosome positioning correlation at the single molecule, and thus, single cell level, that can be used to monitor disease progression and response to therapy.