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.

PathwayPCA: an R/Bioconductor Package for Pathway Based Integrative Analysis of Multi-Omics Data.

The authors present pathwayPCA, an R/Bioconductor package for integrative pathway analysis that utilizes modern statistical methodology, including supervised and adaptive, elastic-net, sparse principal component analysis. pathwayPCA can be applied to continuous, binary, and survival outcomes in studies with multiple covariates and/or interaction effects. It outperforms several alternative methods at identifying disease-associated pathways in integrative analysis using both simulated and real datasets. In addition, several case studies are provided to illustrate pathwayPCA analysis with gene selection, estimating, and visualizing sample-specific pathway activities, identifying sex-specific pathway effects in kidney cancer, and building integrative models for predicting patient prognosis. pathwayPCA is an open-source R package, freely available through the Bioconductor repository. pathwayPCA is expected to be a useful tool for empowering the wider scientific community to analyze and interpret the wealth of available proteomics data, along with other types of molecular data recently made available by Clinical Proteomic Tumor Analysis Consortium and other large consortiums.

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.

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.