An immunoevasive strategy through clinically-relevant pan-cancer genomic and transcriptomic alterations of JAK-STAT signaling components

BACKGROUND

Since its discovery almost three decades ago, the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway has paved the road for understanding inflammatory and immunity processes related to a wide range of human pathologies including cancer. Several studies have demonstrated the importance of JAK-STAT pathway components in regulating tumor initiation and metastatic progression, yet, the extent of how genetic alterations influence patient outcome is far from being understood.

METHODS

Focusing on 133 genes involved in JAK-STAT signaling, we investigated genomic, transcriptomic and clinical profiles of over 18,000 patients representing 21 diverse cancer types. We identified a core set of 28 putative gain- or loss-of-function JAK-STAT genes that correlated with survival outcomes using Cox proportional hazards regression and Kaplan-Meier analyses. Differential expression analyses between high- and low-expressing patient groups were performed to evaluate the consequences of JAK-STAT misexpression.

RESULTS

We found that copy number alterations underpinning transcriptional dysregulation of JAK-STAT pathway genes differ within and between cancer types. Integrated analyses uniting genomic and transcriptomic datasets revealed a core set of JAK-STAT pathway genes that correlated with survival outcomes in brain, renal, lung and endometrial cancers. High JAK-STAT scores were associated with increased mortality rates in brain and renal cancers, but not in lung and endometrial cancers where hyperactive JAK-STAT signaling is a positive prognostic factor. Patients with aberrant JAK-STAT signaling demonstrated pan-cancer molecular features associated with misexpression of genes in other oncogenic pathways (Wnt, MAPK, TGF-β, PPAR and VEGF). Brain and renal tumors with hyperactive JAK-STAT signaling had increased regulatory T cell gene (Treg) expression. A combined model uniting JAK-STAT and Tregs allowed further delineation of risk groups where patients with high JAK-STAT and Treg scores consistently performed the worst.

CONCLUSION

Providing a pan-cancer perspective of clinically-relevant JAK-STAT alterations, this study could serve as a framework for future research investigating anti-tumor immunity using combination therapy involving JAK-STAT and immune checkpoint inhibitors.

Defining the molecular evolution of extrauterine high grade serous carcinoma

OBJECTIVE

High grade serous carcinoma (HGSC) is the most common and most aggressive, subtype of epithelial ovarian cancer. It presents as advanced stage disease with poor prognosis. Recent pathological evidence strongly suggests HGSC arises from the fallopian tube via the precursor lesion; serous tubal intraepithelial carcinoma (STIC). However, further definition of the molecular evolution of HGSC has major implications for both clinical management and research. This study aims to more clearly define the molecular pathogenesis of HGSC.

METHODS

Six cases of HGSC were identified at the Northern Ireland Gynaecological Cancer Centre (NIGCC) that each contained ovarian HGSC (HGSC), omental HGSC (OMT), STIC, normal fallopian tube epithelium (FTE) and normal ovarian surface epithelium (OSE). The relevant formalin-fixed paraffin embedded (FFPE) tissue samples were retrieved from the pathology archive via the Northern Ireland Biobank following attaining ethical approval (NIB11:005). Full microarray-based gene expression profiling was performed on the cohort. The resulting data was analysed bioinformatically and the results were validated in a HGSC-specific in-vitro model.

RESULTS

The carcinogenesis of HGSC was investigated and showed the molecular profile of HGSC to be more closely related to normal FTE than OSE. STIC lesions also clustered closely with HGSC, indicating a common molecular origin.

CONCLUSION

This study provides strong evidence suggesting that extrauterine HGSC arises from the fimbria of the distal fallopian tube. Furthermore, several potential pathways were identified which could be targeted by novel therapies for HGSC. These findings have significant translational relevance for both primary prevention and clinical management of the disease.

Distinct molecular etiologies of male and female hepatocellular carcinoma

BACKGROUND

Sex-differences in cancer occurrence and mortality are evident across tumor types; men exhibit higher rates of incidence and often poorer responses to treatment. Targeted approaches to the treatment of tumors that account for these sex-differences require the characterization and understanding of the fundamental biological mechanisms that differentiate them. Hepatocellular Carcinoma (HCC) is the second leading cause of cancer death worldwide, with the incidence rapidly rising. HCC exhibits a male-bias in occurrence and mortality, but previous studies have failed to explore the sex-specific dysregulation of gene expression in HCC.

METHODS

Here, we characterize the sex-shared and sex-specific regulatory changes in HCC tumors in the TCGA LIHC cohort using combined and sex-stratified differential expression and eQTL analyses.

RESULTS

By using a sex-specific differential expression analysis of tumor and tumor-adjacent samples, we uncovered etiologically relevant genes and pathways differentiating male and female HCC. While both sexes exhibited activation of pathways related to apoptosis and cell cycle, males and females differed in the activation of several signaling pathways, with females showing PPAR pathway enrichment while males showed PI3K, PI3K/AKT, FGFR, EGFR, NGF, GF1R, Rap1, DAP12, and IL-2 signaling pathway enrichment. Using eQTL analyses, we discovered germline variants with differential effects on tumor gene expression between the sexes. 24.3% of the discovered eQTLs exhibit differential effects between the sexes, illustrating the substantial role of sex in modifying the effects of eQTLs in HCC. The genes that showed sex-specific dysregulation in tumors and those that harbored a sex-specific eQTL converge in clinically relevant pathways, suggesting that the molecular etiologies of male and female HCC are partially driven by differential genetic effects on gene expression.

CONCLUSIONS

Sex-stratified analyses detect sex-specific molecular etiologies of HCC. Overall, our results provide new insight into the role of inherited genetic regulation of transcription in modulating sex-differences in HCC etiology and provide a framework for future studies on sex-biased cancers.

Novel candidates in early-onset familial colorectal cancer

In 20-30% of patients suspected of a familial colorectal cancer (CRC) syndrome, no underlying genetic cause is detected. Recent advances in whole exome sequencing have generated evidence for new CRC-susceptibility genes including POLE, POLD1 and NTHL1¸ but many patients remain unexplained. Whole exome sequencing was performed on DNA from nine patients from five different families with familial clusters of CRC in which traditional genetic testing failed to yield a diagnosis. Variants were filtered by minor allele frequencies, followed by prioritization based on in silico prediction tools, and the presence in cancer susceptibility genes or genes in cancer-associated pathways. Effects of frameshift variants on protein structure were modeled using I-Tasser. One known pathogenic variant in POLD1 was detected (p.S478N), together with variants in 17 candidate genes not previously associated with CRC. Additional in silico analysis using SIFT, PROVEAN and PolyPhen on the 14 missense variants indicated a possible damaging effect in nine of 14 variants. Modeling of the insertions/deletions showed a damaging effect of two variants in NOTCH2 and CYP1B1. One family was explained by a mutation in a known familial CRC gene. In the remaining four families, the most promising candidates found are a frameshift NOTCH2 and a missense RAB25 variant. This study provides potential novel candidate variants in unexplained familial CRC patients, however, functional validation is imperative to confirm the role of these variants in CRC tumorigenesis. Additionally, while whole exome sequencing enables detection of variants throughout the exome, other causes explaining the familial phenotype such as multiple single nucleotide polymorphisms accumulating to a polygenic risk or epigenetic events, might be missed with this approach.

Defining the Boundary Between Myelodysplastic Syndromes and Myeloproliferative Neoplasms

In this article we provide a practical and comprehensive review of myeloid neoplasms with overlapping myelodysplastic (MDS) and myeloproliferative (MPN) features, with emphasis on recent updates in classification, particularly the utility of morphologic, cytogenetic, and molecular findings in better defining and classifying these disease entities. We provide the reader with a summary of the most recent developments and updates that have helped further our understanding of the genomic landscape, clinicopathologic features, and prognostic elements of myeloid neoplasms with MDS/MPN features.

Prognostic subclass of intrahepatic cholangiocarcinoma by integrative molecular-clinical analysis and potential targeted approach

BACKGROUND

Although molecular characterization of iCCA has been studied recently, integrative analysis of molecular and clinical characterization has not been fully established. If molecular features of iCCA can be predicted based on clinical findings, we can approach to distinguish targeted treatment. We analyzed RNA sequencing data annotated with clinicopathologic data to clarify molecular-specific clinical features and to evaluate potential therapies for molecular subtypes.

METHODS

We performed next-generation RNA sequencing of 30 surgically resected iCCA from Korean patients and the clinicopathologic features were analyzed. The RNA sequences from 32 iCCA resected from US patients were used for validation.

RESULTS

Patients were grouped into two subclasses on the basis of unsupervised clustering, which showed a difference in 5-year survival rates (48.5% vs 14.2%, p = 0.007) and similar survival outcome in the US samples. In subclass B (poor prognosis), both data sets were similar in higher carcinoembryonic antigen and cancer antigen 19-9 levels, underlying cholangitis, and bile duct-type pathology; in subclass A (better prognosis), there was more frequent viral hepatitis and cholangiolar-type pathology. On pathway analysis, subclass A had enriched liver-related signatures. Subclass B had enriched inflammation-related and TP53 pathways, with more frequent KRAS mutations. CCA cell lines with similar gene expression patterns of subclass A were sensitive to gemcitabine.

CONCLUSIONS

Two molecular subtypes of iCCA with distinct clinicopathological differences were identified. Knowledge of clinical and pathologic characteristics can predict molecular subtypes, and knowledge of different subtype signaling pathways may lead to more rational, targeted approaches to treatment.

Comprehensive validation of liquid-based cytology specimens for next-generation sequencing in cancer genome analysis

In addition to conventional cytology, liquid-based cytology (LBC) is also used for immunocytochemistry and gene analysis. However, an appropriate method to obtain high quality DNA for next-generation sequencing (NGS) using LBC specimens remains controversial. We determined the optimal conditions for fixation with an alcohol-based fixative for LBC and DNA extraction using cultured cancer cell lines and clinical specimens. The extracted DNA was processed for NGS after the DNA quality was confirmed based on the DNA concentration and degree of degradation. The optimal conditions for cultured cells to obtain high quality DNA were to fix the cells at a density of 6 × 10³ or 2 × 10⁴ cells/mL and to use the magnetic bead-based DNA extraction method. Even after storing the fixed cells for 90 days, DNA extracted using the above and other extraction kits, including membrane-based methods, did not undergo degradation. Furthermore, 5-year-old residual LBC samples demonstrated high DNA quality that was suitable for NGS. Furthermore, a cancer genome panel analysis was successfully performed with DNA extracted from cultured cells fixed at 6 × 10³ cells/mL for 90 days, and with DNA from residual LBC samples even after 1 year of storage. Residual LBC samples may be a useful source of DNA for clinical NGS to promote genome-based cancer medicine.

Genetic predisposition to colorectal cancer: syndromes, genes, classification of genetic variants and implications for precision medicine

This article reviews genes and syndromes associated with predisposition to colorectal cancer (CRC), with an overview of gene variant classification. We include updates on the application of preventive and therapeutic measures, focusing on the use of non-steroidal anti-inflammatory drugs (NSAIDs) and immunotherapy. Germline pathogenic variants in genes conferring high or moderate risk to cancer are detected in 6-10% of all CRCs and 20% of those diagnosed before age 50. CRC syndromes can be subdivided into nonpolyposis and polyposis entities, the most common of which are Lynch syndrome and familial adenomatous polyposis, respectively. In addition to known and novel genes associated with highly penetrant CRC risk, identification of pathogenic germline variants in genes associated with moderate-penetrance cancer risk and/or hereditary cancer syndromes not traditionally linked to CRC may have an impact on genetic testing, counseling, and surveillance. The use of multigene panels in genetic testing has exposed challenges in the classification of variants of uncertain significance. We provide an overview of the main classification systems and strategies for improving these. Finally, we highlight approaches for integrating chemoprevention in the care of individuals with genetic predisposition to CRC and use of targeted agents and immunotherapy for treatment of mismatch repair-deficient and hypermutant tumors. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Genomic Signature of Driver Genes Identified by Target Next-Generation Sequencing in Chinese Non-Small Cell Lung Cancer

BACKGROUND

Non-small cell lung cancer (NSCLC) is one of the most common human malignancies and the leading cause of cancer-related death. Over the past few decades, genomic alterations of cancer driver genes have been identified in NSCLC, and molecular testing and targeted therapies have become standard care for lung cancer patients. Here we studied the unique genomic profile of driver genes in Chinese patients with NSCLC by next-generation sequencing (NGS) assay.

MATERIALS AND METHODS

A total of 1,200 Chinese patients with NSCLC were enrolled in this study. The median age was 60 years (range: 26-89), and 83% cases were adenocarcinoma. NGS-based genomic profiling of major lung cancer-related genes was performed on formalin-fixed paraffin-embedded tumor samples and matched blood.

RESULTS

Approximately 73.9% of patients with NSCLC harbored at least one actionable alteration recommended by the National Comprehensive Cancer Network guideline, including epidermal growth factor receptor (EGFR), ALK, ERBB2, MET, BRAF, RET, and ROS1. Twenty-seven patients (2.2%) harbored inherited germline mutations of cancer susceptibility genes. The frequencies of EGFR genomic alterations (both mutations and amplification) and ALK rearrangement were identified as 50.1% and 7.8% in Chinese NSCLC populations, respectively, and significantly higher than the Western population. Fifty-six distinct uncommon EGFR mutations other than L858R, exon19del, exon20ins, or T790M were identified in 18.9% of patients with EGFR-mutant NSCLC. About 7.4% of patients harbored both sensitizing and uncommon mutations, and 11.6% of patients harbored only uncommon EGFR mutations. The uncommon EGFR mutations more frequently combined with the genomic alterations of ALK, CDKN2A, NTRK3, TSC2, and KRAS. In patients <40 years of age, the ALK-positive percentage was up to 28.2%. Moreover, 3.2% of ALK-positive patients harbored multi ALK rearrangements, and seven new partner genes were identified.

CONCLUSION

More unique features of cancer driver genes in Chinese NSCLC were identified by next-generation sequencing. These findings highlighted that NGS technology is more feasible and necessary than other molecular testing methods, and suggested that the special strategies are needed for drug development and targeted therapy for Chinese patients with NSCLC.

IMPLICATIONS FOR PRACTICE

Molecular targeted therapy is now the standard first-line treatment for patients with advanced non-small cell lung cancer (NSCLC). Samples of 1,200 Chinese patients with NSCLC were analyzed through next-generation sequencing to characterize the unique feature of uncommon EGFR mutations and ALK fusion. The results showed that 7.4% of EGFR-mutant patients harbored both sensitizing and uncommon mutations and 11.6% harbored only uncommon mutations. Uncommon EGFR mutations more frequently combined with the genomic alterations of ALK, CDKN2A, NTRK3, TSC2, and KRAS. ALK fusion was more common in younger patients, and the frequency decreased monotonically with age. 3.2% of ALK-positive patients harbored multi ALK rearrangement, and seven new partner genes were identified.

MicroRNAs and their role in environmental chemical carcinogenesis

MicroRNAs (miRNAs) are a class of small, noncoding RNA species that play crucial roles across many biological processes and in the pathogenesis of major diseases, including cancer. Recent studies suggest that the expression of miRNA is altered by certain environmental chemicals, including metals, organic pollutants, cigarette smoke, pesticides and carcinogenic drugs. In addition, extensive studies have indicated the existence and importance of miRNA in different cancers, suggesting that cancer-related miRNAs could serve as potential markers for chemically induced cancers. The altered expression of miRNA was considered to be a vital pathogenic role in xenobiotic-induced cancer development. However, the significance of miRNA in the etiology of cancer and the exact mechanisms by which environmental factors alter miRNA expression remain relatively unexplored. Hence, understanding the interaction of miRNAs with environmental chemicals will provide important information on mechanisms underlying the pathogenesis of chemically induced cancers, and effectively diagnose and treat human cancers resulting from chronic or acute carcinogen exposure. This study presents the current evidence that the miRNA deregulation induced by various chemical carcinogens, different cancers caused by environmental carcinogens and the potentially related genes in the onset or progression of cancer. For each carcinogen, the specifically expressed miRNA may be considered as the early biomarkers of the cancer process. In this review, we also summarize various target genes of the altered miRNA, oncogenes or anti-oncogenes, and the existing evidence regarding the gene regulation mechanisms of cancer caused by environmentally induced miRNA alteration. The future perspective of miRNA may become attractive targets for the diagnosis and treatment of carcinogen-induced cancer.

Genetics for the Women's Health Trainee: A Five-Module Curriculum

Introduction

Genetics is ubiquitous in OB-GYN. However, data suggest that trainees feel underprepared to counsel patients about genetic testing, the nuances of which are becoming increasingly complicated. We sought to develop and implement a genetics curriculum for OB-GYN residents.

Methods

This five-module (screening for fetal aneuploidy, prenatal diagnostic testing, prenatal carrier screening, pedigrees, and cancer genetics), interactive, case-based curriculum is linked to Council on Resident Education in Obstetrics and Gynecology objectives and can stand alone or work as part of an ultrasound or obstetrics rotation. Each module, containing objectives, assigned readings, and cases with answers, is used in a small-group format and can be completed in 20-30 minutes prior to the start of a clinical day. Modules were implemented at two academic centers with first-year OB-GYN residents. Qualitative real-time feedback and summative quantitative feedback from OB-GYN residents were obtained.

Results

Twenty-one OB-GYN residents completed the curriculum, which was well received by trainees and program directors. All residents (100%) felt the curriculum increased knowledge of prenatal genetics and felt more comfortable counseling patients after completion. Seventy-three percent enjoyed the discussion/case-based format; associated articles were found helpful by 100% of trainees. Facilitators enjoyed teaching the curriculum and felt learner knowledge improved dramatically.

Discussion

These low-cost modules were easy to implement and resulted in increased knowledge and confidence in prenatal and cancer genetics. Designed to stand alone and take as little as 20 minutes, the modules provide a helpful adjunct to a women's health rotation or didactic curriculum.

Loss at 16q22.1 identified as a risk factor for intrahepatic recurrence in hepatocellular carcinoma and screening of differentially expressed genes

Copy number alteration (CNA) of chromosome 16, a frequent genetic event in tumors including hepatocellular carcinoma (HCC), has been associated with HCC etiology of hepatitis B virus (HBV) and with clinical outcomes in multiple types of cancer. This study identified CNAs in chromosome 16 in relation to intrahepatic recurrence of HCC in a population with high HBV prevalence, and further screened for differentially expressed genes in recurrence-related CNAs. Array comparative genomic hybridization and expression arrays were used to detect CNAs and gene expression differences, respectively. The associations between CNAs and intrahepatic recurrence were analyzed on 66 patients, follow-up period of 3-73 months. One hundred and nine cases were further evaluated regarding the differentially expressed genes. Losses at 16q and 16p were detected in 62.1% and 51.5% of the 66 cases, respectively. The most recurrent CNAs (with frequency >20%) were losses at 16p13.3-13.2, 16p13.11, 16q11.2-22.1, 16q22.1, 16q22.2-24.2 and 16q24.2. Of the CNAs, 16q22.1 loss was significantly associated with unfavorable intrahepatic recurrence-free survival (P = 0.025). Multivariate Cox analysis identified 16q22.1 loss as an independent risk factor for intrahepatic recurrence (HR = 2.32, 95% CI = 1.26-4.27). A panel of 21 genes, including TRADD, PSMB10, THAP11, CTCF and ESRP2, were significantly downregulated in HCCs with 16q22.1 loss compared to those without the loss. These results suggest that loss at 16q22.1 was associated with increased risk for intrahepatic recurrence of HCC, at least in the HBV-prevalence population. Multiple downregulated genes correlated with the loss were screened.

Transcriptional Regulation in Prostate Cancer

Prostate cancer development involves corruption of the normal prostate transcriptional network, following deregulated expression or mutation of key transcription factors. Here, we provide an overview of the transcription factors that are important in normal prostate homeostasis (NKX3-1, p63, androgen receptor [AR]), primary prostate cancer (ETS family members, c-MYC), castration-resistant prostate cancer (AR, FOXA1), and AR-independent castration-resistant neuroendocrine prostate cancer (RB1, p53, N-MYC). We use functional (in vitro and in vivo) as well as clinical data to discuss evidence that unveils their roles in the initiation and progression of prostate cancer, with an emphasis on results of chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq).

Creutzfeldt astrocytes may be seen in IDH-wildtype glioblastoma and retain expression of DNA repair and chromatin binding proteins

Astrocytes with multiple micronuclei ("Creutzfeldt cells") in a brain biopsy are classically associated with demyelinating disease. However, glioblastoma may also have prominent Creutzfeldt astrocytes, along with granular mitoses. Therefore, Creutzfeldt cells may raise the diagnostic dilemma of high-grade glioma vs tumefactive demyelination. While cases of glioblastoma (GBM) with Creutzfeldt astrocytes have been reported, their clinicopathologic spectrum and genetic features are not understood. Studies have proposed that micronuclei in Creutzfeldt cells are a consequence of DNA damage, or may be susceptible to DNA damage and chromothripsis, but their biology in the context of glioblastoma remains unclear. Based on a challenging index case of GBM with mild hypercellularity, Creutzfeldt astrocytes, and granular mitoses on biopsy, we searched our archives for additional cases with similar histopathologic features. We identified 13 cases, reviewed their clinico-radiologic and pathologic features, and examined them for recurrent genetic alterations via NGS (9 cases) and for evidence of DNA damage by immunohistochemistry for DNA repair and chromatin remodeling proteins. We found that Creutzfeldt cell-rich GBMs were IDH-wildtype with no recurring genetic alterations. To test our hypothesis that micronuclei demonstrate loss of DNA repair or chromatin remodeling proteins, we examined the expression of various proteins (MDM2, p53, MLH1, MSH2, PMS2, MSH6, ATRX, INI1, SATB2, Ki67, pHH3) in Creutzfeldt cell rich-GBM. There was intact expression of DNA repair and chromatin remodeling proteins, with accumulation of p53 and reduced MDM2 expression within micronuclei. In contrast, granular mitoses showed pHH3 expression, confirming these cells are undergoing mitotic division, with no accumulation of p53 and reduced expression of DNA repair proteins. Our results emphasize that Creutzfeldt cells are part of the morphologic spectrum of IDH-wildtype glioblastoma. We did not find a role for DNA damage in the generation of Creutzfeldt cells, as both DNA repair and chromatin remodeling protein expression was retained in these cells.

The NETest: The Clinical Utility of Multigene Blood Analysis in the Diagnosis and Management of Neuroendocrine Tumors

The neuroendocrine neoplasms test (NETest) is a multianalyte liquid biopsy that measures neuroendocrine tumor gene expression in blood. This unique signature precisely defines the biological activity of an individual tumor in real time. The assay meets the 3 critical requirements of an optimal biomarker: diagnostic accuracy, prognostic value, and predictive therapeutic assessment. NETest performance metrics are sensitivity and specificity and in head-to-head comparison are 4-fold to 10-fold more accurate than chromogranin A. NETest accurately identifies completeness of surgery and response to somatostatin analogs. Clinical registry data demonstrate significant clinical utility in watch/wait programs.

An automated microfluidic gene-editing platform for deciphering cancer genes

Gene-editing techniques such as RNA-guided endonuclease systems are becoming increasingly popular for phenotypic screening. Such screens are normally conducted in arrayed or pooled formats. There has been considerable interest in recent years to find new technological methods for conducting these gene-editing assays. We report here the first digital microfluidic method that can automate arrayed gene-editing in mammalian cells. Specifically, this method was useful in culturing lung cancer cells for up to six days, as well as implementing automated gene transfection and knockout procedures. In addition, a standardized imaging pipeline to analyse fluorescently labelled cells was also designed and implemented during these procedures. A gene editing assay for interrogating the MAPK/ERK pathway was performed to show the utility of our platform and to determine the effects of knocking out the RAF1 gene in lung cancer cells. In addition to gene knockout, we also treated the cells with an inhibitor, Sorafenib Tosylate, to determine the effects of enzymatic inhibition. The combination of enzymatic inhibition and guide targeting on device resulted in lower drug concentrations for achieving half-inhibitory effects (IC50) compared to cells treated only with the inhibitor, confirming that lung cancer cells are being successfully edited on the device. We propose that this system will be useful for other types of gene-editing assays and applications related to personalized medicine.

The Discovery of Mutated Driver Pathways in Cancer: Models and Algorithms

The pathogenesis of cancer in human is still poorly understood. With the rapid development of high-throughput sequencing technologies, huge volumes of cancer genomics data have been generated. Deciphering that data poses great opportunities and challenges to computational biologists. One of such key challenges is to distinguish driver mutations, genes as well as pathways from passenger ones. Mutual exclusivity of gene mutations (each patient has no more than one mutation in the gene set) has been observed in various cancer types and thus has been used as an important property of a driver gene set or pathway. In this article, we aim to review the recent development of computational models and algorithms for discovering driver pathways or modules in cancer with the focus on mutual exclusivity-based ones.

A proteomic landscape of diffuse-type gastric cancer

The diffuse-type gastric cancer (DGC) is a subtype of gastric cancer with the worst prognosis and few treatment options. Here we present a dataset from 84 DGC patients, composed of a proteome of 11,340 gene products and mutation information of 274 cancer driver genes covering paired tumor and nearby tissue. DGC can be classified into three subtypes (PX1-3) based on the altered proteome alone. PX1 and PX2 exhibit dysregulation in the cell cycle and PX2 features an additional EMT process; PX3 is enriched in immune response proteins, has the worst survival, and is insensitive to chemotherapy. Data analysis revealed four major vulnerabilities in DGC that may be targeted for treatment, and allowed the nomination of potential immunotherapy targets for DGC patients, particularly for those in PX3. This dataset provides a rich resource for information and knowledge mining toward altered signaling pathways in DGC and demonstrates the benefit of proteomic analysis in cancer molecular subtyping.

Ethnic specific differences in endocrine neoplasms: The role of susceptibility genes

Ethnic disparity in disease incidence, prevalence, and outcome has been documented for a variety of diseases and cancers. Dr. LaSalle D. Leffall was one of the first to note that genetic susceptibility is one important aspect that needed to be studied to better understand cancer disparity. In this article, we cover disparity in endocrine neoplasm presentation and new, emerging genetic data that may explain this disparity.

Promoter-bound METTL3 maintains myeloid leukaemia by m6A-dependent translation control

N6-methyladenosine (m6A) is an abundant internal RNA modification in both coding and non-coding RNAs that is catalysed by the METTL3-METTL14 methyltransferase complex. However, the specific role of these enzymes in cancer is still largely unknown. Here we define a pathway that is specific for METTL3 and is implicated in the maintenance of a leukaemic state. We identify METTL3 as an essential gene for growth of acute myeloid leukaemia cells in two distinct genetic screens. Downregulation of METTL3 results in cell cycle arrest, differentiation of leukaemic cells and failure to establish leukaemia in immunodeficient mice. We show that METTL3, independently of METTL14, associates with chromatin and localizes to the transcriptional start sites of active genes. The vast majority of these genes have the CAATT-box binding protein CEBPZ present at the transcriptional start site, and this is required for recruitment of METTL3 to chromatin. Promoter-bound METTL3 induces m6A modification within the coding region of the associated mRNA transcript, and enhances its translation by relieving ribosome stalling. We show that genes regulated by METTL3 in this way are necessary for acute myeloid leukaemia. Together, these data define METTL3 as a regulator of a chromatin-based pathway that is necessary for maintenance of the leukaemic state and identify this enzyme as a potential therapeutic target for acute myeloid leukaemia.

Automated deconvolution of structured mixtures from heterogeneous tumor genomic data

With increasing appreciation for the extent and importance of intratumor heterogeneity, much attention in cancer research has focused on profiling heterogeneity on a single patient level. Although true single-cell genomic technologies are rapidly improving, they remain too noisy and costly at present for population-level studies. Bulk sequencing remains the standard for population-scale tumor genomics, creating a need for computational tools to separate contributions of multiple tumor clones and assorted stromal and infiltrating cell populations to pooled genomic data. All such methods are limited to coarse approximations of only a few cell subpopulations, however. In prior work, we demonstrated the feasibility of improving cell type deconvolution by taking advantage of substructure in genomic mixtures via a strategy called simplicial complex unmixing. We improve on past work by introducing enhancements to automate learning of substructured genomic mixtures, with specific emphasis on genome-wide copy number variation (CNV) data, as well as the ability to process quantitative RNA expression data, and heterogeneous combinations of RNA and CNV data. We introduce methods for dimensionality estimation to better decompose mixture model substructure; fuzzy clustering to better identify substructure in sparse, noisy data; and automated model inference methods for other key model parameters. We further demonstrate their effectiveness in identifying mixture substructure in true breast cancer CNV data from the Cancer Genome Atlas (TCGA). Source code is available at https://github.com/tedroman/WSCUnmix.

Cancer Subtype Discovery Based on Integrative Model of Multigenomic Data

One major goal of large-scale cancer omics study is to understand molecular mechanisms of cancer and find new biomedical targets. To deal with the high-dimensional multidimensional cancer omics data (DNA methylation, mRNA expression, etc.), which can be used to discover new insight on identifying cancer subtypes, clustering methods are usually used to find an effective low-dimensional subspace of the original data and then cluster cancer samples in the reduced subspace. However, due to data-type diversity and big data volume, few methods can integrate these data and map them into an effective low-dimensional subspace. In this paper, we develop a dimension-reduction and data-integration method for indentifying cancer subtypes, named Scluster. First, Scluster, respectively, projects the different original data into the principal subspaces by an adaptive sparse reduced-rank regression method. Then, a fused patient-by-patient network is obtained for these subgroups through a scaled exponential similarity kernel method. Finally, candidate cancer subtypes are identified using spectral clustering method. We demonstrate the efficiency of our Scluster method using three cancers by jointly analyzing mRNA expression, miRNA expression, and DNA methylation data. The evaluation results and analyses show that Scluster is effective for predicting survival and identifies novel cancer subtypes of large-scale multi-omics data.

A composite network of conserved and tissue specific gene interactions reveals possible genetic interactions in glioma

Differential co-expression network analyses have recently become an important step in the investigation of cellular differentiation and dysfunctional gene-regulation in cell and tissue disease-states. The resulting networks have been analyzed to identify and understand pathways associated with disorders, or to infer molecular interactions. However, existing methods for differential co-expression network analysis are unable to distinguish between various forms of differential co-expression. To close this gap, here we define the three different kinds (conserved, specific, and differentiated) of differential co-expression and present a systematic framework, CSD, for differential co-expression network analysis that incorporates these interactions on an equal footing. In addition, our method includes a subsampling strategy to estimate the variance of co-expressions. Our framework is applicable to a wide variety of cases, such as the study of differential co-expression networks between healthy and disease states, before and after treatments, or between species. Applying the CSD approach to a published gene-expression data set of cerebral cortex and basal ganglia samples from healthy individuals, we find that the resulting CSD network is enriched in genes associated with cognitive function, signaling pathways involving compounds with well-known roles in the central nervous system, as well as certain neurological diseases. From the CSD analysis, we identify a set of prominent hubs of differential co-expression, whose neighborhood contains a substantial number of genes associated with glioblastoma. The resulting gene-sets identified by our CSD analysis also contain many genes that so far have not been recognized as having a role in glioblastoma, but are good candidates for further studies. CSD may thus aid in hypothesis-generation for functional disease-associations.

Two different, mutually exclusively distributed, TP53 mutations in ovarian and peritoneal tumor tissues of a serous ovarian cancer patient: indicative for tumor origin?

High-grade serous ovarian cancer (HGSOC) is characterized by a TP53 mutation rate of up to 96.7% and associated with a more aggressive tumor biology. The origin of HGSOC is thought to arise either from fallopian tube secretory cells or the ovarian surface epithelium/inclusion cysts, the former with more evidence. Peritoneal tumor spread is heterogeneous, either excessive in the peritoneum (with miliary appearance) or more confined to the ovaries with only few (bigger and exophytically growing) peritoneal implants. Using RNA sequencing and DNA digital droplet polymerase chain reaction (PCR), we identified two different functional TP53 mutations in one HGSOC patient: one exclusively in the ovarian tumor mass and the other exclusively in ascites tumor cells, peritoneal tumor masses, and a lymph node metastasis. In blood, both mutations could be detected, the one from the peritoneal tumors with much higher frequency, presumably because of the higher tumor load. We conclude that this mutually exclusive distribution of two different TP53 mutations in different tumor tissues indicates the development of two independent carcinomas in the peritoneal cavity, probably one originating from a precancerous lesion in the fallopian tube and the other from the ovaries. In addition, in the patient's ascites CD45 and EpCAM, double-positive cells were found-proliferating but testing negative for the above-mentioned TP53 mutations. This mutually exclusive distribution of two TP53 mutations is probably further evidence that HGSOC can originate either from the fallopian tube or (more seldom) the ovaries, the former more prone for excessive peritoneal tumor spread.

Transplantation of engineered organoids enables rapid generation of metastatic mouse models of colorectal cancer

Colorectal cancer (CRC) is a leading cause of death in the developed world, yet facile preclinical models that mimic the natural stages of CRC progression are lacking. Through the orthotopic engraftment of colon organoids we describe a broadly usable immunocompetent CRC model that recapitulates the entire adenoma-adenocarcinoma-metastasis axis in vivo. The engraftment procedure takes less than 5 minutes, shows efficient tumor engraftment in two-thirds of mice, and can be achieved using organoids derived from genetically engineered mouse models (GEMMs), wild-type organoids engineered ex vivo, or from patient-derived human CRC organoids. In this model, we describe the genotype and time-dependent progression of CRCs from adenocarcinoma (6 weeks), to local disseminated disease (11-12 weeks), and spontaneous metastasis (>20 weeks). Further, we use the system to show that loss of dysregulated Wnt signaling is critical for the progression of disseminated CRCs. Thus, our approach provides a fast and flexible means to produce tailored CRC mouse models for genetic studies and pre-clinical investigation.