Numerous large-scale genomic studies of matched tumor-normal samples have established the somatic landscapes of most cancer types. However, the downstream analysis of data from somatic mutations entails a number of computational and statistical approaches, requiring usage of independent software and numerous tools. Here, we describe an R Bioconductor package, Maftools, which offers a multitude of analysis and visualization modules that are commonly used in cancer genomic studies, including driver gene identification, pathway, signature, enrichment, and association analyses. Maftools only requires somatic variants in Mutation Annotation Format (MAF) and is independent of larger alignment files. With the implementation of well-established statistical and computational methods, Maftools facilitates data-driven research and comparative analysis to discover novel results from publicly available data sets. In the present study, using three of the well-annotated cohorts from The Cancer Genome Atlas (TCGA), we describe the application of Maftools to reproduce known results. More importantly, we show that Maftools can also be used to uncover novel findings through integrative analysis.

Esophageal squamous cell carcinoma (ESCC) is among the most common malignancies, but little is known about its spatial intratumoral heterogeneity (ITH) and temporal clonal evolutionary processes. To address this, we performed multiregion whole-exome sequencing on 51 tumor regions from 13 ESCC cases and multiregion global methylation profiling for 3 of these 13 cases. We found an average of 35.8% heterogeneous somatic mutations with strong evidence of ITH. Half of the driver mutations located on the branches of tumor phylogenetic trees targeted oncogenes, including PIK3CA, NFE2L2 and MTOR, among others. By contrast, the majority of truncal and clonal driver mutations occurred in tumor-suppressor genes, including TP53, KMT2D and ZNF750, among others. Interestingly, phyloepigenetic trees robustly recapitulated the topological structures of the phylogenetic trees, indicating a possible relationship between genetic and epigenetic alterations. Our integrated investigations of spatial ITH and clonal evolution provide an important molecular foundation for enhanced understanding of tumorigenesis and progression in ESCC.

Key Points: MLL3 acts as tumor suppressor in FLT3-ITD AML. The existence of DNMT3A mutations in remission samples implies that the DNMT3A mutant clone can survive induction chemotherapy.

Long non-coding RNAs (lncRNAs) are expressed in tissue-specific pattern, but it is not clear how these are regulated. We aimed to identify squamous cell carcinoma (SCC)-specific lncRNAs and investigate mechanisms that control their expression and function. We studied expression patterns and functions of 4 SCC-specific lncRNAs. We obtained 113 esophageal SCC (ESCC) and matched non-tumor esophageal tissues from a hospital in Shantou City, China, and performed quantitative reverse transcription PCR assays to measure expression levels of LINC01503. We collected clinical data from patients and compared expression levels with survival times. LINC01503 was knocked down using small interfering RNAs and oligonucleotides in TE7, TE5 and KYSE510 cell lines and overexpressed in KYSE30 cells. Cells were analyzed by chromatin immunoprecipitation sequencing, luciferase reporter assays, colony formation, migration and invasion, and mass spectrometry analyses. Cells were injected into nude mice and growth of xenograft tumors was measured. LINC01503 interaction with proteins was studied using fluorescence in situ hybridization, RNA pulldown, and RNA immunoprecipitation analyses. We identified a lncRNA, LINC01503, which is regulated by a super-enhancer and is expressed at significantly higher levels in esophageal and head and neck SCCs than in non-tumor tissues. High levels in SCCs correlated with shorter survival times of patients. The transcription factor TP63 bound to the super-enhancer at the LINC01503 locus and activated its transcription. Expression of LINC01503 in ESCC cell lines increased their proliferation, colony formation, migration and invasion. Knockdown of LINC01503 in SCC cells reduced their proliferation, colony formation, migration and invasion, and growth of xenograft tumors in nude mice. Expression of LINC01503 in ESCC cell lines reduced ERK2 dephosphorylation by DUSP6, leading to activation of ERK signaling via MAPK. LINC01503 disrupted the interaction between EBP1 and the p85 subunit of PI3K, increasing AKT signaling. We identified a lncRNA, LINC01503, that is increased in SCC cells compared with non-tumor cells. Increased expression of LINC01503 promotes ESCC cell proliferation, migration, invasion, and growth of xenograft tumors. It might be developed as a biomarker of aggressive SCCs in patients.

Liposarcoma (LPS) is the most common type of soft tissue sarcoma accounting for 20% of all adult sarcomas. Due to absence of clinically effective treatment options in inoperable situations and resistance to chemotherapeutics, a critical need exists to identify novel therapeutic targets. We analyzed LPS genomic landscape using SNP arrays, whole exome sequencing and targeted exome sequencing to uncover the genomic information for development of specific anti-cancer targets. SNP array analysis indicated known amplified genes (MDM2, CDK4, HMGA2) and important novel genes (UAP1, MIR557, LAMA4, CPM, IGF2, ERBB3, IGF1R). Carboxypeptidase M (CPM), recurrently amplified gene in well-differentiated/de-differentiated LPS was noted as a putative oncogene involved in the EGFR pathway. Notable deletions were found at chromosome 1p (RUNX3, ARID1A), chromosome 11q (ATM, CHEK1) and chromosome 13q14.2 (MIR15A, MIR16-1). Significantly and recurrently mutated genes (false discovery rate < 0.05) included PLEC (27%), MXRA5 (21%), FAT3 (24%), NF1 (20%), MDC1 (10%), TP53 (7%) and CHEK2 (6%). Further, in vitro and in vivo functional studies provided evidence for the tumor suppressor role for Neurofibromin 1 (NF1) gene in different subtypes of LPS. Pathway analysis of recurrent mutations demonstrated signaling through MAPK, JAK-STAT, Wnt, ErbB, axon guidance, apoptosis, DNA damage repair and cell cycle pathways were involved in liposarcomagenesis. Interestingly, we also found mutational and copy number heterogeneity within a primary LPS tumor signifying the importance of multi-region sequencing for cancer-genome guided therapy. In summary, these findings provide insight into the genomic complexity of LPS and highlight potential druggable pathways for targeted therapeutic approach.