DNA amplification is a ubiquitous mechanism of oncogene activation in lung and other cancers

Preclinical Models for Functional Precision Lung Cancer Research

Patient-centered precision oncology strives to deliver individualized cancer care. In lung cancer, preclinical models and technological innovations have become critical in advancing this approach. Preclinical models enable deeper insights into tumor biology and enhance the selection of appropriate systemic therapies across chemotherapy, targeted therapies, immunotherapies, antibody-drug conjugates, and emerging investigational treatments. While traditional human lung cancer cell lines offer a basic framework for cancer research, they often lack the tumor heterogeneity and intricate tumor-stromal interactions necessary to accurately predict patient-specific clinical outcomes. Patient-derived xenografts (PDXs), however, retain the original tumor's histopathology and genetic features, providing a more reliable model for predicting responses to systemic therapeutics, especially molecularly targeted therapies. For studying immunotherapies and antibody-drug conjugates, humanized PDX mouse models, syngeneic mouse models, and genetically engineered mouse models (GEMMs) are increasingly utilized. Despite their value, these in vivo models are costly, labor-intensive, and time-consuming. Recently, patient-derived lung cancer organoids (LCOs) have emerged as a promising in vitro tool for functional precision oncology studies. These LCOs demonstrate high success rates in growth and maintenance, accurately represent the histology and genomics of the original tumors and exhibit strong correlations with clinical treatment responses. Further supported by advancements in imaging, spatial and single-cell transcriptomics, proteomics, and artificial intelligence, these preclinical models are reshaping the landscape of drug development and functional precision lung cancer research. This integrated approach holds the potential to deliver increasingly accurate, personalized treatment strategies, ultimately enhancing patient outcomes in lung cancer.

Exploring the potential biological significance of KDELR family genes in lung adenocarcinoma

The Lys-Asp-Glu-Leu receptor (KDELR) family genes play critical roles in a variety of biological processes in different tumors. Our study aimed to provide a comprehensive analysis of the potential roles of KDELRs in lung adenocarcinoma (LUAD). Utilizing data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, as well as clinical samples, we conducted a series of analyses and validations using R software tools and various online resources. The results showed that KDELR family genes and proteins were highly expressed and associated with a poor prognosis of LUAD. Promoter hypomethylation and the competing endogenous RNA (ceRNA) network of PCAT6/hsa-miR-326/KDELR1 might be potential causes of aberrant KDELR1 overexpression in LUAD. Three key Transcription factors (TFs) (SPI1, EP300, and MAZ) and a TFs-miRNAs-KDELRs network (involving 11 TFs) might be involved in modulating KDELRs expression abnormalities. Gene Set Enrichment Analysis (GSEA) indicated enrichment of genes highly expressing KDELR1, KDELR2, and KDELR3 in MTORC1_SIGNALING, P53_PATHWAY, and ANGIOGENESIS. Negative correlations between KDELRs expression and CD8 + T cell infiltration, as well as CTLA-4 expression. Our multiple analyses suggested that the KDELRs are important signaling molecules in LUAD. These results provided novel insights for developing prognostic markers and novel therapies of LUAD.

Intestinal Tuft Cells: Morphology, Function, and Implications for Human Health

Tuft cells are a rare and morphologically distinct chemosensory cell type found throughout many organs, including the gastrointestinal tract. These cells were identified by their unique morphologies distinguished by large apical protrusions. Ultrastructural data have begun to describe the molecular underpinnings of their cytoskeletal features, and tuft cell-enriched cytoskeletal proteins have been identified, although the connection of tuft cell morphology to tuft cell functionality has not yet been established. Furthermore, tuft cells display variations in function and identity between and within tissues, leading to the delineation of distinct tuft cell populations. As a chemosensory cell type, they display receptors that are responsive to ligands specific for their environment. While many studies have demonstrated the tuft cell response to protists and helminths in the intestine, recent research has highlighted other roles of tuft cells as well as implicated tuft cells in other disease processes including inflammation, cancer, and viral infections. Here, we review the literature on the cytoskeletal structure of tuft cells. Additionally, we focus on new research discussing tuft cell lineage, ligand-receptor interactions, tuft cell tropism, and the role of tuft cells in intestinal disease. Finally, we discuss the implication of tuft cell-targeted therapies in human health and how the morphology of tuft cells may contribute to their functionality.

Assessment of MYC and TERT copy number variations in lung cancer using digital PCR

OBJECTIVE

Lung cancer is the second most frequent cancer type and the most common cause of cancer-related deaths worldwide. Alteration of gene copy numbers are associated with lung cancer and the determination of copy number variations (CNV) is appropriate for the discrimination between tumor and non-tumor tissue in lung cancer. As telomerase reverse transcriptase (TERT) and v-myc avian myelocytomatosis viral oncogene homolog (MYC) play a role in lung cancer the aims of this study were the verification of our recent results analyzing MYC CNV in tumor and non-tumor tissue of lung cancer patients using an independent study group and the assessment of TERT CNV as an additional marker.

RESULTS

TERT and MYC status was analyzed using digital PCR (dPCR) in tumor and adjacent non-tumor tissue samples of 114 lung cancer patients. The difference between tumor and non-tumor samples were statistically significant (p < 0.0001) for TERT and MYC. Using a predefined specificity of 99% a sensitivity of 41% and 51% was observed for TERT and MYC, respectively. For the combination of TERT and MYC the overall sensitivity increased to 60% at 99% specificity. We demonstrated that a combination of markers increases the performance in comparison to individual markers. Additionally, the determination of CNV using dPCR might be an appropriate tool in precision medicine.

lncENST Suppress the Warburg Effect Regulating the Tumor Progress by the Nkx2-5/ErbB2 Axis in Hepatocellular Carcinoma

The therapeutic efficacy of radiofrequency ablation (RFA) against liver cancer is often limited by proliferation and metastasis of residual tumor cells. These phenomena are closely associated with the Warburg effect, wherein ErbB2 is activated. While RFA inhibits the Warburg effect of residual tumor cells at the early stage, the specific mechanisms remain unclear. We explored the regulatory relationship between the long noncoding RNA ENST00000570843.1 (lncENST) and ErbB2 using lentiviral transfection of lncENST and ErbB2 overexpression/interference vectors in in vitro and in vivo models of hepatocellular carcinoma in the presence of sublethal heat at 50°C. ErbB2-mediated Warburg effect was suppressed by lncENST, as manifested by reduced glucose uptake and lactic acid production in SMMC-7721 cells. lncENST also increased tumor apoptosis and inhibited tumor progression in nude Balb/c mice for up to 28 days after RFA. Additionally, we predicted through bioinformatic analysis that the promoter of ErbB2 binds to the transcription factor Nkx2-5, resulting in a negative regulatory effect. This speculation was confirmed by chromatin immunoprecipitation of the Nkx2-5 protein and ErbB2, indicating that ErbB2 transcription was curbed by Nkx2-5. We propose that lncENST downplays the Warburg effect in residual tumor cells by downregulating ErbB2 via Nkx2-5 activation. This study is aimed at providing molecular targets that can prevent residual tumor cell proliferation after RFA, with clinical significance in hepatocellular carcinoma treatment.

Zoo: Selecting Transcriptomic and Methylomic Biomarkers by Ensembling Animal-Inspired Swarm Intelligence Feature Selection Algorithms

Biological omics data such as transcriptomes and methylomes have the inherent “large p small n” paradigm, i.e., the number of features is much larger than that of the samples. A feature selection (FS) algorithm selects a subset of the transcriptomic or methylomic biomarkers in order to build a better prediction model. The hidden patterns in the FS solution space make it challenging to achieve a feature subset with satisfying prediction performances. Swarm intelligence (SI) algorithms mimic the target searching behaviors of various animals and have demonstrated promising capabilities in selecting features with good machine learning performances. Our study revealed that different SI-based feature selection algorithms contributed complementary searching capabilities in the FS solution space, and their collaboration generated a better feature subset than the individual SI feature selection algorithms. Nine SI-based feature selection algorithms were integrated to vote for the selected features, which were further refined by the dynamic recursive feature elimination framework. In most cases, the proposed Zoo algorithm outperformed the existing feature selection algorithms on transcriptomics and methylomics datasets.

Lung Cancer Models Reveal Severe Acute Respiratory Syndrome Coronavirus 2-Induced Epithelial-to-Mesenchymal Transition Contributes to Coronavirus Disease 2019 Pathophysiology

INTRODUCTION

Coronavirus disease 2019 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which enters host cells through the cell surface proteins ACE2 and TMPRSS2.

METHODS

Using a variety of normal and malignant models and tissues from the aerodigestive and respiratory tracts, we investigated the expression and regulation of ACE2 and TMPRSS2.

RESULTS

We find that ACE2 expression is restricted to a select population of epithelial cells. Notably, infection with SARS-CoV-2 in cancer cell lines, bronchial organoids, and patient nasal epithelium induces metabolic and transcriptional changes consistent with epithelial-to-mesenchymal transition (EMT), including up-regulation of ZEB1 and AXL, resulting in an increased EMT score. In addition, a transcriptional loss of genes associated with tight junction function occurs with SARS-CoV-2 infection. The SARS-CoV-2 receptor, ACE2, is repressed by EMT through the transforming growth factor-β, ZEB1 overexpression, and onset of EGFR tyrosine kinase inhibitor resistance. This suggests a novel model of SARS-CoV-2 pathogenesis in which infected cells shift toward an increasingly mesenchymal state, associated with a loss of tight junction components with acute respiratory distress syndrome-protective effects. AXL inhibition and ZEB1 reduction, as with bemcentinib, offer a potential strategy to reverse this effect.

CONCLUSIONS

These observations highlight the use of aerodigestive and, especially, lung cancer model systems in exploring the pathogenesis of SARS-CoV-2 and other respiratory viruses and offer important insights into the potential mechanisms underlying the morbidity and mortality of coronavirus disease 2019 in healthy patients and patients with cancer alike.

Copy number variation in triple negative breast cancer samples associated with lymph node metastasis

Triple negative breast cancer (TNBC) is a highly metastatic and aggressive subtype of breast cancer and cases presenting with lymph node involvement have worse outcomes. This study aimed to determine the regions of copy number variation (CNV) associated with lymph node metastasis in TNBC patients. CNV analyses were performed in a study cohort of 23 invasive ductal carcinomas (IDCs), 12 lymph node metastases (LNmets), and 7 normal adjacent tissues (NATs); as well as in an independent cohort containing 70 TNBC IDCs and the same 7 NATs. CNV-associated genes were analyzed using GO-enrichment and Pathway analysis. The prognostic role for genes showing CNV-based changes in messenger RNA expression was determined using the Kaplan-Meier plotter database. For the IDCs, there were a number of variations that were common in both the study and independent cohorts in the amplified regions of 1q, 8q, 19 (p and q), 2p, 5p and the deleted regions in 8p followed by 5q, and 19p. The most frequently amplified regions in the LNmets of the study cohort were 4q28.3, 2p, 3q24, 1q21.2, 10p, 12p11.1, 8q, 20p11.22-20p11.21, 21q22.13, 6p22.1 and the most frequently deleted regions were in 1p36.23, 4q21.1 and 5q. A total of 686 (441 amplified and 245 deleted) genes were associated with LNmets. The LNmet-associated genes were highly enriched for “regulation of complement activation,” “regulation of protein activation cascade,” “regulation of humoral immune response,” “oxytocin signalling pathway,” and “TRAIL binding” pathways. Moreover, 6/686 LNmet-associated genes showed CNV-based changes in their mRNA expression of which, high expression of ASPM and KIF14 was significantly associated with worse relapse-free survival. This study has identified several CNV regions in TNBC that could play a major role in metastasis to the lymph node.

Extracellular signal-regulated kinase mediates chromatin rewiring and lineage transformation in lung cancer

Lineage transformation between lung cancer subtypes is a poorly understood phenomenon associated with resistance to treatment and poor patient outcomes. Here, we aimed to model this transition to define underlying biological mechanisms and identify potential avenues for therapeutic intervention. Small cell lung cancer (SCLC) is neuroendocrine in identity and, in contrast to non-SCLC (NSCLC), rarely contains mutations that drive the MAPK pathway. Likewise, NSCLCs that transform to SCLC concomitantly with development of therapy resistance downregulate MAPK signaling, suggesting an inverse relationship between pathway activation and lineage state. To test this, we activated MAPK in SCLC through conditional expression of mutant KRAS or EGFR, which revealed suppression of the neuroendocrine differentiation program via ERK. We found that ERK induces the expression of ETS factors that mediate transformation into a NSCLC-like state. ATAC-seq demonstrated ERK-driven changes in chromatin accessibility at putative regulatory regions and global chromatin rewiring at neuroendocrine and ETS transcriptional targets. Further, ERK-mediated induction of ETS factors as well as suppression of neuroendocrine differentiation were dependent on histone acetyltransferase activities of CBP/p300. Overall, we describe how the ERK-CBP/p300-ETS axis promotes a lineage shift between neuroendocrine and non-neuroendocrine lung cancer phenotypes and provide rationale for the disruption of this program during transformation-driven resistance to targeted therapy.

Identification of a targetable KRAS-mutant epithelial population in non-small cell lung cancer

Lung cancer is the leading cause of cancer deaths. Tumor heterogeneity, which hampers development of targeted therapies, was herein deconvoluted via single cell RNA sequencing in aggressive human adenocarcinomas (carrying Kras-mutations) and comparable murine model. We identified a tumor-specific, mutant-KRAS-associated subpopulation which is conserved in both human and murine lung cancer. We previously reported a key role for the oncogene BMI-1 in adenocarcinomas. We therefore investigated the effects of in vivo PTC596 treatment, which affects BMI-1 activity, in our murine model. Post-treatment, MRI analysis showed decreased tumor size, while single cell transcriptomics concomitantly detected near complete ablation of the mutant-KRAS-associated subpopulation, signifying the presence of a pharmacologically targetable, tumor-associated subpopulation. Our findings therefore hold promise for the development of a targeted therapy for KRAS-mutant adenocarcinomas.

Association of Body Mass Index With Somatic Mutations in Breast Cancer

Background

The relationship between body mass index (BMI) and the prognosis or treatment response in patients with breast cancer has been demonstrated in previous studies, but the somatic mutation profiles in breast cancer patients with different BMIs have not been explored.

Methods

In the present study, the somatic mutation profiles in 421 female breast cancer patients who were stratified into three subgroups based on BMI (normal weight, overweight/obese, and underweight) were investigated. Capture-based targeted sequencing was performed using a panel comprising 520 cancer-related genes.

Results

A total of 3547 mutations were detected in 390 genes. In breast cancer patients with different BMI statuses, the tumors exhibited high mutation frequency and burden. TP53 was the most common gene in the three groups, followed by PIK3CA, ERBB2, and CDK12. Meanwhile, the mutation hotspots in TP53 and PIK3CA were the same in the three BMI groups. More JAK1 mutations were identified in underweight patients than those in normal patients. Except for JAK1, differentially mutated genes in postmenopausal patients were completely different from those in premenopausal patients. The distribution of mutation types was significantly different among BMI groups in the postmenopausal group. Underweight patients in the postmenopausal group harbored more TP53 mutations, more amplifications, and more mutations in genes involved in the WNT signaling pathway.

Conclusions

Our next-generation sequencing (NGS)-based gene panel analysis revealed the gene expression profiles of breast cancer patients with different BMI statuses. Although genes with high mutation frequency and burden were found in different BMI groups, some subtle differences could not be ignored. JAK1 mutations might play a vital role in the progression of breast cancer in underweight patients, and this needs further analysis. Postmenopausal underweight patients with breast cancer have more aggressive characteristics, such as TP53 mutations, more amplifications, and more mutations in genes involved in the WNT signaling pathway. This study provides new evidence for understanding the characteristics of breast cancer patients with different BMIs.

A dynamic recursive feature elimination framework (dRFE) to further refine a set of OMIC biomarkers

MOTIVATION

A feature selection algorithm may select the subset of features with the best associations with the class labels. The recursive feature elimination (RFE) is a heuristic feature screening framework and has been widely used to select the biological OMIC biomarkers. This study proposed a dynamic recursive feature elimination (dRFE) framework with more flexible feature elimination operations. The proposed dRFE was comprehensively compared with 11 existing feature selection algorithms and five classifiers on the eight difficult transcriptome datasets from a previous study, the ten newly collected transcriptome datasets and the five methylome datasets.

RESULTS

The experimental data suggested that the regular RFE framework did not perform well, and dRFE outperformed the existing feature selection algorithms in most cases. The dRFE-detected features achieved Acc=1.0000 for the two methylome datasets GSE53045 and GSE66695. The best prediction accuracies of the dRFE-detected features were 0.9259, 0.9424, and 0.8601 for the other three methylome datasets GSE74845, GSE103186, and GSE80970, respectively. Four transcriptome datasets received Acc=1.0000 using the dRFE-detected features, and the prediction accuracies for the other six newly collected transcriptome datasets were between 0.6301 and 0.9917.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Lung cancer models reveal SARS-CoV-2-induced EMT contributes to COVID-19 pathophysiology

COVID-19 is an infectious disease caused by SARS-CoV-2, which enters host cells via the cell surface proteins ACE2 and TMPRSS2. Using a variety of normal and malignant models and tissues from the aerodigestive and respiratory tracts, we investigated the expression and regulation of ACE2 and TMPRSS2. We find that ACE2 expression is restricted to a select population of highly epithelial cells. Notably, infection with SARS-CoV-2 in cancer cell lines, bronchial organoids, and patient nasal epithelium, induces metabolic and transcriptional changes consistent with epithelial to mesenchymal transition (EMT), including upregulation of ZEB1 and AXL, resulting in an increased EMT score. Additionally, a transcriptional loss of genes associated with tight junction function occurs with SARS-CoV-2 infection. The SARS-CoV-2 receptor, ACE2, is repressed by EMT via TGFbeta, ZEB1 overexpression and onset of EGFR TKI inhibitor resistance. This suggests a novel model of SARS-CoV-2 pathogenesis in which infected cells shift toward an increasingly mesenchymal state, associated with a loss of tight junction components with acute respiratory distress syndrome-protective effects. AXL-inhibition and ZEB1-reduction, as with bemcentinib, offers a potential strategy to reverse this effect. These observations highlight the utility of aerodigestive and, especially, lung cancer model systems in exploring the pathogenesis of SARS-CoV-2 and other respiratory viruses, and offer important insights into the potential mechanisms underlying the morbidity and mortality of COVID-19 in healthy patients and cancer patients alike.

Digital PCR for the Analysis of MYC Copy Number Variation in Lung Cancer

Background

MYC (v-myc avian myelocytomatosis viral oncogene homolog) is one of the most frequently amplified genes in lung tumors. For the analysis of gene copy number variations, dPCR (digital PCR) is an appropriate tool. The aim of our study was the assessment of dPCR for the detection of MYC copy number variations (CNV) in lung tissue considering clinicopathological parameters. Material and Methods. MYC status was analyzed with dPCR as well as qPCR (quantitative PCR) using gDNA (genomic DNA) from tumor and adjacent nontumor tissue samples of lung cancer patients. The performance of MYC was estimated based on the AUC (area under curve).

Results

The results of the MYC amplification correlated significantly between dPCR and qPCR (r _S = 0.81, P < 0.0001). The MYC copy number revealed by dPCR showed statistically significant differences between tumor and adjacent nontumor tissues. For discrimination, a sensitivity of 43% and a specificity of 99% were calculated, representing 55 true-positive and one false-positive tests. No statistically significant differences could be observed for age, sex, and smoking status or the clinicopathological parameters (histological subtype, grade, and stage).

Conclusion

The results of the study show that dPCR is an accurate and reliable method for the determination of MYC copy numbers. The application is characterized by high specificity and moderate sensitivity. MYC amplification is a common event in lung cancer patients, and it is indicated that the determination of the MYC status might be useful in clinical diagnostics.

Supramolecular Polysaccharide Nanotheranostics that Inhibit Cancer Cells Growth and Monitor Targeted Therapy Response

Targeted anticancer therapies directed against specific molecular drivers of tumors are emerging as effective treatment strategies, however, monitoring their response is still challenging. Current clinical imaging techniques that measure either morphological or metabolic changes in the tumor are not always indicative of clinical outcome due to delayed or variable responses. Here, dual-stage polysaccharide-based supramolecular nanotheranostics (SPN) were designed that enable co-delivery of kinase inhibitor and an activatable imaging probe.

Methods: The SPNs were prepared by supramolecular assembly of two components, polysaccharide construct conjugated with kinase inhibitor-function activatable probe and kinase inhibitor- β-cyclodextrin conjugate. Physiochemical characterization of SPNs including size, stability, zeta potential and pH-responsiveness of the assembly was performed. The efficacy of SPNs in inducing cancer cell death by inhibition of kinase signaling and imaging the response was evaluated in murine BRAFV^600E melanoma (D4M) and triple-negative breast cancer (4T1) cell lines. Finally, the in vivo efficacy was investigated in D4M melanoma tumor model.

Results: The polysaccharide-constructs along with kinase inhibitor- β-cyclodextrin conjugates self-assemble to produce SPNs of around 200 nm in diameter and were stable for over a week under physiologically relevant conditions. The SPNs exhibited enhanced cytotoxic effect and significant inhibition of kinase signaling as compared to the free inhibitor. In vitro imaging studies confirmed their enzyme-activatable therapy response tracking abilities both in cancer cells and tumor spheroids. Furthermore, SPN treated mice exhibited better tumor growth inhibition as compared to the control groups and therapy response could be imaged at both early (24-48h) and later time points.

Conclusion: These findings demonstrate that the supramolecular polysaccharide nanotheranostics can not only inhibit kinase signaling pathway in aggressive tumor, but also monitor targeted therapy response early.

Targeting the cyclin-dependent kinase 5 in metastatic melanoma

The cyclin-dependent kinase 5 (CDK5), originally described as a neuronal-specific kinase, is also frequently activated in human cancers. Using conditional CDK5 knockout mice and a mouse model of highly metastatic melanoma, we found that CDK5 is dispensable for the growth of primary tumors. However, we observed that ablation of CDK5 completely abrogated the metastasis, revealing that CDK5 is essential for the metastatic spread. In mouse and human melanoma cells CDK5 promotes cell invasiveness by directly phosphorylating an intermediate filament protein, vimentin, thereby inhibiting assembly of vimentin filaments. Chemical inhibition of CDK5 blocks the metastatic spread of patient-derived melanomas in patient-derived xenograft (PDX) mouse models. Hence, inhibition of CDK5 might represent a very potent therapeutic strategy to impede the metastatic dissemination of malignant cells.

Identification of PGRMC1 as a Candidate Oncogene for Head and Neck Cancers and Its Involvement in Metabolic Activities

Progesterone Receptor Membrane Component 1 (PGRMC1/Sigma-2 receptor) is located on chromosome Xq21 and encodes a haem-containing protein that interacts with epidermal growth factor receptor (EGFR) and cytochromes P450, with function in tumor proliferation and chemoresistance. Although the over-expression of PGRMC1 reported in many different types of human cancers, systematic analysis of its oncogenic role of PGRMC1 has not been performed for any cancer. In this work, we analyzed the transcriptomics, genomics, and clinical data of 498 head-neck squamous cell carcinoma (HNSC) samples from the public-accessible database, The Cancer Genome Atlas (TCGA). The Cox regression was performed to calculate the hazard ratio (HR) of PGRMC1 expression as a prognosis feature for overall survival (OS). Our results demonstrated that PGRMC1 expression served as a predictor for worse OS (HR = 1.95, P = 0.0005) in head-neck squamous cell carcinoma. And the over-expression of PGRMC1 was strongly correlated with various metabolic process activity as well as cancer metastasis and cell proliferation features in human head-neck squamous cell carcinoma patient's cohort. Besides, the over-expression and unfavorable prognosis value of PGRMC1 were also observed in many other cancer types. This study provides insights into the potential oncogenic functional significance of PGRMC1 in cancer research.

A kinase of many talents: non-neuronal functions of CDK5 in development and disease

The cyclin-dependent kinase 5 (CDK5) represents an unusual member of the family of cyclin-dependent kinases, which is activated upon binding to non-cyclin p35 and p39 proteins. The role of CDK5 in the nervous system has been very well established. In addition, there is growing evidence that CDK5 is also active in non-neuronal tissues, where it has been postulated to affect a variety of functions such as the immune response, angiogenesis, myogenesis, melanogenesis and regulation of insulin levels. Moreover, high levels of CDK5 have been observed in different tumour types, and CDK5 was proposed to play various roles in the tumorigenic process. In this review, we discuss these various CDK5 functions in normal physiology and disease, and highlight the therapeutic potential of targeting CDK5.

Interpreting pathways to discover cancer driver genes with Moonlight

Cancer driver gene alterations influence cancer development, occurring in oncogenes, tumor suppressors, and dual role genes. Discovering dual role cancer genes is difficult because of their elusive context-dependent behavior. We define oncogenic mediators as genes controlling biological processes. With them, we classify cancer driver genes, unveiling their roles in cancer mechanisms. To this end, we present Moonlight, a tool that incorporates multiple -omics data to identify critical cancer driver genes. With Moonlight, we analyze 8000+ tumor samples from 18 cancer types, discovering 3310 oncogenic mediators, 151 having dual roles. By incorporating additional data (amplification, mutation, DNA methylation, chromatin accessibility), we reveal 1000+ cancer driver genes, corroborating known molecular mechanisms. Additionally, we confirm critical cancer driver genes by analysing cell-line datasets. We discover inactivation of tumor suppressors in intron regions and that tissue type and subtype indicate dual role status. These findings help explain tumor heterogeneity and could guide therapeutic decisions.

Epithelial tumor suppressor ELF3 is a lineage-specific amplified oncogene in lung adenocarcinoma

Gene function in cancer is often cell type-specific. The epithelial cell-specific transcription factor ELF3 is a documented tumor suppressor in many epithelial tumors yet displays oncogenic properties in others. Here, we show that ELF3 is an oncogene in the adenocarcinoma subtype of lung cancer (LUAD), providing genetic, functional, and clinical evidence of subtype specificity. We discover a region of focal amplification at chromosome 1q32.1 encompassing the ELF3 locus in LUAD which is absent in the squamous subtype. Gene dosage and promoter hypomethylation affect the locus in up to 80% of LUAD analyzed. ELF3 expression was required for tumor growth and a pan-cancer expression network analysis supports its subtype and tissue specificity. We further show that ELF3 displays strong prognostic value in LUAD but not LUSC. We conclude that, contrary to many other tumors of epithelial origin, ELF3 is an oncogene and putative therapeutic target in LUAD.

Tissue context can dictate why a gene can have seemingly opposing functions in different settings. ELF3 is tumor suppressive in many cancers of epithelial origin but in lung cancer, the authors describe an oncogenic role in the adenocarcinoma histology of non-small cell lung cancer.

Integrative Genomic Analyses Identifies GGA2 as a Cooperative Driver of EGFR-Mediated Lung Tumorigenesis

INTRODUCTION

Targeted therapies for lung adenocarcinoma (LUAD) have improved patient outcomes; however, drug resistance remains a major problem. One strategy to achieve durable response is to develop combination-based therapies that target both mutated oncogenes and key modifiers of oncogene-driven tumorigenesis. This is based on the premise that mutated oncogenes, although necessary, are not sufficient for malignant transformation. We aimed to uncover genetic alterations that cooperate with mutant EGFR during LUAD development.

METHODS

We performed integrative genomic analyses, combining copy number, gene expression and mutational information for over 500 LUAD tumors. Co-immunoprecipitation and Western blot analysis were performed in LUAD cell lines to confirm candidate interactions while RNA interference and gene overexpression were used for in vitro and in vivo functional assessment.

RESULTS

We identified frequent amplifications/deletions of chromosomal regions affecting the activity of genes specifically in the context of EGFR mutation, including amplification of the mutant EGFR allele and deletion of dual specificity phosphatase 4 (DUSP4), which have both previously been reported. In addition, we identified the novel amplification of a segment of chromosome arm 16p in mutant-EGFR tumors corresponding to increased expression of Golgi Associated, Gamma Adaptin Ear Containing, ARF Binding Protein 2 (GGA2), which functions in protein trafficking and sorting. We found that GGA2 interacts with EGFR, increases EGFR protein levels and modifies EGFR degradation after ligand stimulation. Furthermore, we show that overexpression of GGA2 enhances EGFR mediated transformation while GGA2 knockdown reduces the colony and tumor forming ability of EGFR mutant LUAD.

CONCLUSIONS

These data suggest that overexpression of GGA2 in LUAD tumors results in the accumulation of EGFR protein and increased EGFR signaling, which helps drive tumor progression. Thus, GGA2 plays a cooperative role with EGFR during LUAD development and is a potential therapeutic target for combination-based strategies in LUAD.

The oncogenic neurotrophin receptor tropomyosin-related kinase variant, TrkAIII

Oncogenes derived from the neurotrophin receptor tropomyosin-related kinase TrkA act as drivers in sub-populations of a wide-range of human cancers. This, combined with a recent report that both adult and childhood cancers driven by novel oncogenic TrkA chimeric-fusions exhibit profound, long-lived therapeutic responses to the Trk inhibitor Larotrectinib, highlights the need to improve clinical detection of TrkA oncogene-driven cancers in order to maximise this novel therapeutic potential. Cancers potentially driven by TrkA oncogenes include a proportion of paediatric neuroblastomas (NBs) that express the alternative TrkA splice variant TrkAIII, which exhibits exon 6, 7 and 9 skipping and oncogenic-activity that depends upon deletion of the extracellular D4 Ig-like domain. In contrast to fully spliced TrkA, which exhibits tumour suppressor activity in NB and associates with good prognosis, TrkAIII associates with advanced stage metastatic disease, post therapeutic relapse and worse prognosis, induces malignant transformation of NIH-3T3 cells and exhibits oncogenic activity in NB models. TrkAIII induction in NB cells is stress-regulated by conditions that mimic hypoxia or perturbate the ER with potential to change TrkA tumour-suppressing signals into oncogenic TrkAIII signals within the stressful tumour microenvironment. In contrast to cell surface TrkA, TrkAIII re-localises to intracellular pre-Golgi membranes, centrosomes and mitochondria, within which it exhibits spontaneous ligand-independent activation, triggering a variety of mechanisms that promote tumorigenicity and malignant behaviour, which impact the majority of cancer hallmarks. In this review, we present updates on TrkAIII detection and association with human malignancies, the multiple ways TrkAIII exerts oncogenic activity and potential therapeutic approaches for TrkAIII expressing cancers, with particular reference to NB.

The Emerging Role of Cdk5 in Cancer

Cdk5 is an atypical cyclin-dependent kinase that is well characterized for its role in the central nervous system rather than in the cell cycle. However Cdk5 has been recently implicated in the development and progression of a variety of cancers including breast, lung, colon, pancreatic, melanoma, thyroid and brain tumors. This broad pro-tumorigenic role makes Cdk5 a promising drug target for the development of new cancer therapies. Here we review the contribution of Cdk5 to molecular mechanisms that confer upon tumors the ability to grow, proliferate and disseminate to secondary organs, as well as resistance to chemotherapies. We subsequently discuss existing and new strategies for targeting Cdk5 and its downstream mechanisms as anti-cancer treatments.

eTumorType, An Algorithm of Discriminating Cancer Types for Circulating Tumor Cells or Cell-free DNAs in Blood

With the technology development on detecting circulating tumor cells (CTCs) and cell-free DNAs (cfDNAs) in blood, serum, and plasma, non-invasive diagnosis of cancer becomes promising. A few studies reported good correlations between signals from tumor tissues and CTCs or cfDNAs, making it possible to detect cancers using CTCs and cfDNAs. However, the detection cannot tell which cancer types the person has. To meet these challenges, we developed an algorithm, eTumorType, to identify cancer types based on copy number variations (CNVs) of the cancer founding clone. eTumorType integrates cancer hallmark concepts and a few computational techniques such as stochastic gradient boosting, voting, centroid, and leading patterns. eTumorType has been trained and validated on a large dataset including 18 common cancer types and 5327 tumor samples. eTumorType produced high accuracies (0.86-0.96) and high recall rates (0.79-0.92) for predicting colon, brain, prostate, and kidney cancers. In addition, relatively high accuracies (0.78-0.92) and recall rates (0.58-0.95) have also been achieved for predicting ovarian, breast luminal, lung, endometrial, stomach, head and neck, leukemia, and skin cancers. These results suggest that eTumorType could be used for non-invasive diagnosis to determine cancer types based on CNVs of CTCs and cfDNAs.

Global analysis of chromosome 1 genes among patients with lung adenocarcinoma, squamous carcinoma, large-cell carcinoma, small-cell carcinoma, or non-cancer

The present study aimed at investigating genetic variations, specific signal pathways, or biological processes of chromosome 1 genes between subtypes and stages of lung cancer and prediction of selected targeting genes for patient survival rate. About 537 patients with lung adenocarcinoma (ADC), 140 with lung squamous carcinoma (SCC), 9 with lung large-cell carcinoma (LCC), 56 with small-cell lung cancer (SCLC), and 590 without caner were integrated from 16 databases and analyzed in the present study. Three (ASPM, CDC20, KIAA1799) or 28 genes significantly up- or down-expressed in four subtypes of lung cancer. The activated cell division and down-regulated immune responses were identified in patients with lung cancer. Keratinocyte development associated genes S100 and SPRR families dominantly up-expressed in SCC and AKT3 and NRAS in SCLC. Subtype-specific genes of ADC, SCC, LCC, or SCLC were also identified. C1orf106, CAPN8, CDC20, COL11A1, CRABP2, and NBPF9 up-expressed at four stages of ADC. Fifty six related with keratinocytes or potassium channels up-expressed in three stages of SCC. CDC20, IL10, ECM1, GABPB2, CRABP2, and COL11A1 significantly predicted the poor overall survival of ADC patients and S100A2 and TIMM17A in SCC patients. Our data indicate that a number of altered chromosome 1 genes have the subtype and stage specificities of lung cancer and can be considered as diagnostic and prognosis biomarkers.

CDK5 functions as a tumor promoter in human colorectal cancer via modulating the ERK5-AP-1 axis

Abnormal expression of cyclin-dependent kinase 5 (CDK5) has been found in several human cancers, whereas the role of CDK5 in the malignant development of colorectal cancer (CRC) has not been well characterized. Here we investigated the role of CDK5 in CRC and found that its expression was much higher in CRC tissues than that in normal tissues with a higher expression level of CDK5 closely correlating to advanced American Joint Committee on Cancer (AJCC) stage, poor differentiation, increased tumor size and poor prognosis of CRC. Biological function experiments showed that CDK5 regulated CRC cell proliferation and metastasis ability. Whole-genome microarray analysis, co-immunoprecipitation, in vitro kinase assay, western blotting, luciferase reporter assays and electrophoretic mobility shift assay (EMSA) showed that CDK5 could directly phosphorylate ERK5 at threonine (Thr) 732 and finally modulate the oncogenic ERK5–AP-1 axis. Further researches showed that CDK5–ERK5–AP-1 axis could promote progression of CRC carcinogenesis and had a significant correlation in human CRC samples. In summary, this study revealed the functional and mechanistic links between CDK5 and the oncogenic ERK5–AP-1 signaling pathway in the pathogenesis of CRC. These findings suggest that CDK5 has an important role in CRC development and may serve as a potential therapeutic target for CRC.

Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing

BACKGROUND

B chromosomes are dispensable and variable karyotypic elements found in some species of animals, plants and fungi. They often originate from duplications and translocations of host genomic regions or result from hybridization. In most species, little is known about their DNA content. Here we perform high-throughput sequencing and analysis of B chromosomes of roe deer and brocket deer, the only representatives of Cetartiodactyla known to have B chromosomes.

RESULTS

In this study we developed an approach to identify genomic regions present on chromosomes by high-throughput sequencing of DNA generated from flow-sorted chromosomes using degenerate-oligonucleotide-primed PCR. Application of this method on small cattle autosomes revealed a previously described KIT gene region translocation associated with colour sidedness. Implementing this approach to B chromosomes from two cervid species, Siberian roe deer (Capreolus pygargus) and grey brocket deer (Mazama gouazoubira), revealed dramatically different genetic content: roe deer B chromosomes consisted of two duplicated genomic regions (a total of 1.42-1.98 Mbp) involving three genes, while grey brocket deer B chromosomes contained 26 duplicated regions (a total of 8.28-9.31 Mbp) with 34 complete and 21 partial genes, including KIT and RET protooncogenes, previously found on supernumerary chromosomes in canids. Sequence variation analysis of roe deer B chromosomes revealed a high frequency of mutations and increased heterozygosity due to either amplification within B chromosomes or divergence between different Bs. In contrast, grey brocket deer B chromosomes were found to be more homogeneous and resembled autosomes in patterns of sequence variation. Similar tendencies were observed in repetitive DNA composition.

CONCLUSIONS

Our data demonstrate independent origins of B chromosomes in the grey brocket and roe deer. We hypothesize that the B chromosomes of these two cervid species represent different stages of B chromosome sequences evolution: probably nascent and similar to autosomal copies in brocket deer, highly derived in roe deer. Based on the presence of the same orthologous protooncogenes in canids and brocket deer Bs we argue that genomic regions involved in B chromosome formation are not random. In addition, our approach is also applicable to the characterization of other evolutionary and clinical rearrangements.

Conservation of protein abundance patterns reveals the regulatory architecture of the EGFR-MAPK pathway

Various genetic mutations associated with cancer are known to alter cell signaling, but it is not clear whether they dysregulate signaling pathways by altering the abundance of pathway proteins. Using a combination of RNA sequencing and ultrasensitive targeted proteomics, we defined the primary components-16 core proteins and 10 feedback regulators-of the epidermal growth factor receptor (EGFR)-mitogen-activated protein kinase (MAPK) pathway in normal human mammary epithelial cells and then quantified their absolute abundance across a panel of normal and breast cancer cell lines as well as fibroblasts. We found that core pathway proteins were present at very similar concentrations across all cell types, with a variance similar to that of proteins previously shown to display conserved abundances across species. In contrast, EGFR and transcriptionally controlled feedback regulators were present at highly variable concentrations. The absolute abundance of most core proteins was between 50,000 and 70,000 copies per cell, but the adaptors SOS1, SOS2, and GAB1 were found at far lower amounts (2000 to 5000 copies per cell). MAPK signaling showed saturation in all cells between 3000 and 10,000 occupied EGFRs, consistent with the idea that adaptors limit signaling. Our results suggest that the relative stoichiometry of core MAPK pathway proteins is very similar across different cell types, with cell-specific differences mostly restricted to variable amounts of feedback regulators and receptors. The low abundance of adaptors relative to EGFR could be responsible for previous observations that only a fraction of total cell surface EGFR is capable of rapid endocytosis, high-affinity binding, and mitogenic signaling.

The well-accepted notion that gene amplification contributes to increased expression still remains, after all these years, a reasonable but unproven assumption

“Gene amplification causes overexpression” is a longstanding and well-accepted concept in cancer genetics. However, raking the whole literature, we find only statistical analyses showing a positive correlation between gene copy number and expression level, but do not find convincing experimental corroboration for this notion, for most of the amplified oncogenes in cancers. Since an association does not need to be an actual causal relation, in our opinion, this widespread notion still remains a reasonable but unproven assumption awaiting experimental verification.

An immunohistochemical study of cyclin-dependent kinase 5 (CDK5) expression in non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC): a possible prognostic biomarker

Background

Cyclin-dependent kinase 5 (CDK5) is an atypical CDK which plays a vital role in several cancers via regulating migration and motility of cancer cells. However, the clinicopathological impact and function of CDK5 in lung cancer remain poorly understood. The present study was aimed at exploring expression and clinicopathological significance of CDK5 in lung cancer.

Methods

There were 395 samples of lung tissue including 365 lung tumors (339 non-small cell lung cancers and 26 small cell lung cancers) and 30 samples of normal lung. CDK5 expression was detected by immunohistochemistry on lung tissue microarrays.

Results

Over expression was detected in lung cancer compared with normal lung tissues (P = 0.001). Furthermore, area under curve (AUC) of receiver operating characteristic (ROC) of CDK5 was 0.685 (95 % CI 0.564~0.751, P = 0.004). In lung cancer, we also discovered close correlations between CDK5 and pathological grading (r = 0.310, P < 0.001), TNM stage (r = 0.155, P = 0.003), and lymph node metastasis (r = 0.279, P < 0.001) by using Spearman analysis. In two subgroups of non-small cell lung cancer (NSCLC) and small cell lung cancer (SCLC), the expression of CDK5 was also higher than that of normal lung tissue, respectively (P = 0.001 and P = 0.004). Moreover, in NSCLCs, Spearman analysis revealed that expression of CDK5 was correlated with TNM stages (r = 0.129, P = 0.017), lymph node metastasis (r = 0.365, P < 0.001), and pathological grading (r = 0.307, P < 0.001), respectively. The significant correlation was also found between CDK5 expression and TNM stages (r = 0.415, P = 0.049) and lymphatic metastasis (r = 0.469, P = 0.024) in SCLCs.

Conclusions

The results of this present study suggest that the CDK5 expression is associated with several clinicopathological factors linked with poorer prognosis.

Patterns of Chromosomal Aberrations in Solid Tumors

Chromosomal abnormalities are a defining feature of solid tumors. Such cytogenetic alterations are mainly classified into structural chromosomal aberrations and copy number alterations, giving rise to aneuploid karyotypes. The increasing detection of these genetic changes allowed the description of specific tumor entities and the associated patterns of gene expression. In fact, tumor-specific landscapes of gross genomic copy number changes, including aneuploidies of entire chromosome arms and chromosomes result in a global deregulation of the transcriptome of cancer cells. Furthermore, the molecular characterization of cytogenetic abnormalities has provided insights into the mechanisms of tumorigenesis and has, in a few instances, led to the clinical implementation of effective diagnostic and prognostic tools, as well as treatment strategies that target a specific genetic abnormality.

TrkA is amplified in malignant melanoma patients and induces an anti-proliferative response in cell lines

BACKGROUND

The nerve growth factor (NGF) receptor tyrosine-kinase TrkA is a well-known determinant of the melanocytic lineage, through modulation of the MAPK and AKT cascades. While TrkA gene is frequently rearranged in cancers, its involvement in malignant melanoma (MM) development is still unclear.

METHODS

We analyzed a dataset of primary cutaneous MM (n = 31) by array comparative genomic hybridization (aCGH), to identify genomic amplifications associated with tumor progression. The analysis was validated by genomic quantitative PCR (qPCR) on an extended set of cases (n = 64) and the results were correlated with the clinical outcome. To investigate TrkA molecular pathways and cellular function, we generated inducible activation of the NGF-TrkA signaling in human MM cell lines.

RESULTS

We identified amplification of 1q23.1, where the TrkA locus resides, as a candidate hotspot implicated in the progression of MM. Across 40 amplicons detected, segmental amplification of 1q23.1 showed the strongest association with tumor thickness. By validation of the analysis, TrkA gene amplification emerged as a frequent event in primary melanomas (50 % of patients), and correlated with worse clinical outcome. However, experiments in cell lines revealed that induction of the NGF-TrkA signaling produced a phenotype of dramatic suppression of cell proliferation through inhibition of cell division and pronounced intracellular vacuolization, in a way straightly dependent on NGF activation of TrkA. These events were triggered via MAPK activity but not via AKT, and involved p21(cip1) protein increase, compatibly with a mechanism of oncogene-induced growth arrest.

CONCLUSIONS

Taken together, our findings point to TrkA as a candidate oncogene in MM and support a model in which the NGF-TrkA-MAPK pathway may mediate a trade-off between neoplastic transformation and adaptive anti-proliferative response.

Targeting cyclin dependent kinase 5 in hepatocellular carcinoma--A novel therapeutic approach

BACKGROUND & AIMS

For a long time cyclin dependent kinase 5 (Cdk5) was thought to be exclusively important in neuronal cells. However, increasing evidence recently suggests a function of Cdk5 in cancer progression. In this study, we examined the role of Cdk5 and its therapeutic accessibility in hepatocellular carcinoma (HCC), a highly chemoresistant cancer with poor prognosis and paramount clinical importance in order to develop novel targeted therapies for systemic treatment.

METHODS

Expression and activity of Cdk5 was analyzed in a human HCC tissue microarray, human patient samples and HCC cell lines. To characterize Cdk5 functions and signaling pathways in HCC, we applied genetic downregulation and pharmacologic inhibition in various approaches including cell based assays and mouse xenograft models.

RESULTS

Expression and activity of Cdk5 was increased in human HCC tissues as compared to normal liver tissues. Functional ablation of Cdk5 significantly decreased HCC cell proliferation and clonogenic survival. Moreover, genetic and pharmacological inhibition of Cdk5 showed in vivo efficacy in HCC xenograft mouse models. Investigating the mechanisms behind these functional effects revealed that Cdk5 is most active in the nucleus of cells in G2/M phase. Cdk5 regulates DNA damage response by phosphorylating ataxia telangiectasia mutated (ATM) kinase and thereby influencing its downstream cascade. Consequently, combination of Cdk5 inhibition with DNA-damage-inducing chemotherapeutics synergistically inhibited HCC tumor progression in vitro and in vivo.

CONCLUSIONS

In summary, we introduce Cdk5 as a novel drugable target for HCC treatment and suggest the combination of Cdk5 inhibition and DNA damaging agents as a novel therapeutic approach.

MYC and human telomerase gene (TERC) copy number gain in early-stage non-small cell lung cancer

OBJECTIVES

We investigated the frequency of MYC and TERC increased gene copy number (GCN) in early-stage non-small cell lung cancer (NSCLC) and evaluated the correlation of these genomic imbalances with clinicopathologic parameters and outcome.

MATERIALS AND METHODS

Tumor tissues were obtained from 113 resected NSCLCs. MYC and TERC GCNs were tested by fluorescence in situ hybridization (FISH) according to the University of Colorado Cancer Center (UCCC) criteria and based on the receiver operating characteristic (ROC) classification.

RESULTS

When UCCC criteria were applied, 41 (36%) cases for MYC and 41 (36%) cases for TERC were considered FISH-positive. MYC and TERC concurrent FISH-positive was observed in 12 cases (11%): 2 (17%) cases with gene amplification and 10 (83%) with high polysomy. By using the ROC analysis, high MYC (mean ≥ 2.83 copies/cell) and TERC (mean ≥ 2.65 copies/cell) GCNs were observed in 60 (53.1%) cases and 58 (51.3%) cases, respectively. High TERC GCN was associated with squamous cell carcinoma (SCC) histology (P=0.001). In univariate analysis, increased MYC GCN was associated with shorter overall survival (P=0.032 [UCCC criteria] or P=0.02 [ROC classification]), whereas high TERC GCN showed no association. In multivariate analysis including stage and age, high MYC GCN remained significantly associated with worse overall survival using both the UCCC criteria (P=0.02) and the ROC classification (P=0.008).

CONCLUSIONS

Our results confirm MYC as frequently amplified in early-stage NSCLC and increased MYC GCN as a strong predictor of worse survival. Increased TERC GCN does not have prognostic impact but has strong association with squamous histology.

DNA copy number analysis of metastatic urothelial carcinoma with comparison to primary tumors

Background

To date, there have been no reports characterizing the genome-wide somatic DNA chromosomal copy-number alteration landscape in metastatic urothelial carcinoma. We sought to characterize the DNA copy-number profile in a cohort of metastatic samples and compare them to a cohort of primary urothelial carcinoma samples in order to identify changes that are associated with progression from primary to metastatic disease.

Methods

Using molecular inversion probe array analysis we compared genome-wide chromosomal copy-number alterations between 30 metastatic and 29 primary UC samples. Whole transcriptome RNA-Seq analysis was also performed in primary and matched metastatic samples which was available for 9 patients.

Results

Based on a focused analysis of 32 genes in which alterations may be clinically actionable, there were significantly more amplifications/deletions in metastases (8.6% vs 4.5%, p < 0.001). In particular, there was a higher frequency of E2F3 amplification in metastases (30% vs 7%, p = 0.046). Paired primary and metastatic tissue was available for 11 patients and 3 of these had amplifications of potential clinical relevance in metastases that were not in the primary tumor including ERBB2, CDK4, CCND1, E2F3, and AKT1. The transcriptional activity of these amplifications was supported by RNA expression data.

Conclusions

The discordance in alterations between primary and metastatic tissue may be of clinical relevance in the era of genomically directed precision cancer medicine.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1192-2) contains supplementary material, which is available to authorized users.

EZH2 inhibition sensitizes BRG1 and EGFR mutant lung tumours to TopoII inhibitors

Non-small-cell lung cancer is the leading cause of cancer-related death worldwide. Chemotherapies such as the topoisomerase II (TopoII) inhibitor etoposide effectively reduce disease in a minority of patients with this cancer; therefore, alternative drug targets, including epigenetic enzymes, are under consideration for therapeutic intervention. A promising potential epigenetic target is the methyltransferase EZH2, which in the context of the polycomb repressive complex 2 (PRC2) is well known to tri-methylate histone H3 at lysine 27 (H3K27me3) and elicit gene silencing. Here we demonstrate that EZH2 inhibition has differential effects on the TopoII inhibitor response of non-small-cell lung cancers in vitro and in vivo. EGFR and BRG1 mutations are genetic biomarkers that predict enhanced sensitivity to TopoII inhibitor in response to EZH2 inhibition. BRG1 loss-of-function mutant tumours respond to EZH2 inhibition with increased S phase, anaphase bridging, apoptosis and TopoII inhibitor sensitivity. Conversely, EGFR and BRG1 wild-type tumours upregulate BRG1 in response to EZH2 inhibition and ultimately become more resistant to TopoII inhibitor. EGFR gain-of-function mutant tumours are also sensitive to dual EZH2 inhibition and TopoII inhibitor, because of genetic antagonism between EGFR and BRG1. These findings suggest an opportunity for precision medicine in the genetically complex disease of non-small-cell lung cancer.

Identifying core gene modules in glioblastoma based on multilayer factor-mediated dysfunctional regulatory networks through integrating multi-dimensional genomic data

The driver genetic aberrations collectively regulate core cellular processes underlying cancer development. However, identifying the modules of driver genetic alterations and characterizing their functional mechanisms are still major challenges for cancer studies. Here, we developed an integrative multi-omics method CMDD to identify the driver modules and their affecting dysregulated genes through characterizing genetic alteration-induced dysregulated networks. Applied to glioblastoma (GBM), the CMDD identified a core gene module of 17 genes, including seven known GBM drivers, and their dysregulated genes. The module showed significant association with shorter survival of GBM. When classifying driver genes in the module into two gene sets according to their genetic alteration patterns, we found that one gene set directly participated in the glioma pathway, while the other indirectly regulated the glioma pathway, mostly, via their dysregulated genes. Both of the two gene sets were significant contributors to survival and helpful for classifying GBM subtypes, suggesting their critical roles in GBM pathogenesis. Also, by applying the CMDD to other six cancers, we identified some novel core modules associated with overall survival of patients. Together, these results demonstrate integrative multi-omics data can identify driver modules and uncover their dysregulated genes, which is useful for interpreting cancer genome.

Loss of the Notch effector RBPJ promotes tumorigenesis

Kulic et al. show that RBPJ, a transcriptional repressor of Notch, is frequently deleted in human cancers and can function as a tumor suppressor. Loss of RBPJ acts to derepress target gene promoters, allowing Notch-independent activation by alternate transcription factors that promote tumor growth.

Aberrant Notch activity is oncogenic in several malignancies, but it is unclear how expression or function of downstream elements in the Notch pathway affects tumor growth. Transcriptional regulation by Notch is dependent on interaction with the DNA-binding transcriptional repressor, RBPJ, and consequent derepression or activation of associated gene promoters. We show here that RBPJ is frequently depleted in human tumors. Depletion of RBPJ in human cancer cell lines xenografted into immunodeficient mice resulted in activation of canonical Notch target genes, and accelerated tumor growth secondary to reduced cell death. Global analysis of activated regions of the genome, as defined by differential acetylation of histone H4 (H4ac), revealed that the cell death pathway was significantly dysregulated in RBPJ-depleted tumors. Analysis of transcription factor binding data identified several transcriptional activators that bind promoters with differential H4ac in RBPJ-depleted cells. Functional studies demonstrated that NF-κB and MYC were essential for survival of RBPJ-depleted cells. Thus, loss of RBPJ derepresses target gene promoters, allowing Notch-independent activation by alternate transcription factors that promote tumorigenesis.

Molecular characterization of immortalized normal and dysplastic oral cell lines

BACKGROUND

Cell lines have been developed for modeling cancer and cancer progression. The molecular background of these cell lines is often unknown to those using them to model disease behaviors. As molecular alterations are the ultimate drivers of cell phenotypes, having an understanding of the molecular make-up of these systems is critical for understanding the disease biology modeled.

METHODS

Six immortalized normal, one immortalized dysplasia, one self-immortalized dysplasia, and two primary normal cell lines derived from oral tissues were analyzed for DNA copy number changes and changes in both mRNA and miRNA expression using SMRT-v.2 genome-wide tiling comparative genomic hybridization arrays, Agilent Whole Genome 4x44k expression arrays, and Exiqon V2.M-RT-PCR microRNA Human panels.

RESULTS

DNA copy number alterations were detected in both normal and dysplastic immortalized cell lines-as well as in the single non-immortalized dysplastic cell line. These lines were found to have changes in expression of genes related to cell cycle control as well as alterations in miRNAs that are deregulated in clinical oral squamous cell carcinoma tissues. Immortal lines-whether normal or dysplastic-had increased disruption in expression relative to primary lines. All data are available as a public resource.

CONCLUSIONS

Molecular profiling experiments have identified DNA, mRNA, and miRNA alterations for a panel of normal and dysplastic oral tissue cell lines. These data are a valuable resource to those modeling diseases of the oral mucosa, and give insight into the selection of model cell lines and the interpretation of data from those lines.

Proteomic analysis of the human cyclin-dependent kinase family reveals a novel CDK5 complex involved in cell growth and migration

Cyclin-dependent kinases (CDKs) are the catalytic subunits of a family of mammalian heterodimeric serine/threonine kinases that play critical roles in the control of cell-cycle progression, transcription, and neuronal functions. However, the functions, substrates, and regulation of many CDKs are poorly understood. To systematically investigate these features of CDKs, we conducted a proteomic analysis of the CDK family and identified their associated protein complexes in two different cell lines using a modified SAINT (Significance Analysis of INTeractome) method. The mass spectrometry data were deposited to ProteomeXchange with identifier PXD000593 and DOI 10.6019/PXD000593. We identified 753 high-confidence candidate interaction proteins (HCIPs) in HEK293T cells and 352 HCIPs in MCF10A cells. We subsequently focused on a neuron-specific CDK, CDK5, and uncovered two novel CDK5-binding partners, KIAA0528 and fibroblast growth factor (acidic) intracellular binding protein (FIBP), in non-neuronal cells. We showed that these three proteins form a stable complex, with KIAA0528 and FIBP being required for the assembly and stability of the complex. Furthermore, CDK5-, KIAA0528-, or FIBP-depleted breast cancer cells displayed impaired proliferation and decreased migration, suggesting that this complex is required for cell growth and migration in non-neural cells. Our study uncovers new aspects of CDK functions, which provide direction for further investigation of these critical protein kinases.

Recurring DNA copy number gain at chromosome 9p13 plays a role in the activation of multiple candidate oncogenes in progressing oral premalignant lesions

Genomic alteration at chromosome 9p has been previously reported as a frequent and critical event in oral premalignancy. While this alteration is typically reported as a loss driven by selection for CDKN2A deactivation (at 9p21.3), we detect a recurrent DNA copy number gain of ∼2.49 Mbp at chromosome 9p13 in oral premalignant lesions (OPLs) that later progressed to invasive lesions. This recurrent alteration event has been validated using fluorescence in situ hybridization in an independent set of OPLs. Analysis of publicly available gene expression datasets aided in identifying three oncogene candidates that may have driven selection for DNA copy number increases in this region (VCP, DCTN3, and STOML2). We performed in vitro silencing and activation experiments for each of these genes in oral cancer cell lines and found that each gene is independently capable of upregulating proliferation and anchorage-independent growth. We next analyzed the activity of each of these genes in biopsies of varying histological grades that were obtained from a diseased oral tissue field in a single patient, finding further molecular evidence of parallel activation of VCP, DCTN3, and STOML2 during progression from normal healthy tissue to invasive oral carcinoma. Our results support the conclusion that DNA gain at 9p13 is important to the earliest stages of oral tumorigenesis and that this alteration event likely contributes to the activation of multiple oncogene candidates capable of governing oral cancer phenotypes.

Co-expression network analysis identifies Spleen Tyrosine Kinase (SYK) as a candidate oncogenic driver in a subset of small-cell lung cancer

Background

Oncogenic mechanisms in small-cell lung cancer remain poorly understood leaving this tumor with the worst prognosis among all lung cancers. Unlike other cancer types, sequencing genomic approaches have been of limited success in small-cell lung cancer, i.e., no mutated oncogenes with potential driver characteristics have emerged, as it is the case for activating mutations of epidermal growth factor receptor in non-small-cell lung cancer. Differential gene expression analysis has also produced SCLC signatures with limited application, since they are generally not robust across datasets. Nonetheless, additional genomic approaches are warranted, due to the increasing availability of suitable small-cell lung cancer datasets. Gene co-expression network approaches are a recent and promising avenue, since they have been successful in identifying gene modules that drive phenotypic traits in several biological systems, including other cancer types.

Results

We derived an SCLC-specific classifier from weighted gene co-expression network analysis (WGCNA) of a lung cancer dataset. The classifier, termed SCLC-specific hub network (SSHN), robustly separates SCLC from other lung cancer types across multiple datasets and multiple platforms, including RNA-seq and shotgun proteomics. The classifier was also conserved in SCLC cell lines. SSHN is enriched for co-expressed signaling network hubs strongly associated with the SCLC phenotype. Twenty of these hubs are actionable kinases with oncogenic potential, among which spleen tyrosine kinase (SYK) exhibits one of the highest overall statistical associations to SCLC. In patient tissue microarrays and cell lines, SCLC can be separated into SYK-positive and -negative. SYK siRNA decreases proliferation rate and increases cell death of SYK-positive SCLC cell lines, suggesting a role for SYK as an oncogenic driver in a subset of SCLC.

Conclusions

SCLC treatment has thus far been limited to chemotherapy and radiation. Our WGCNA analysis identifies SYK both as a candidate biomarker to stratify SCLC patients and as a potential therapeutic target. In summary, WGCNA represents an alternative strategy to large scale sequencing for the identification of potential oncogenic drivers, based on a systems view of signaling networks. This strategy is especially useful in cancer types where no actionable mutations have emerged.

YEATS4 is a novel oncogene amplified in non-small cell lung cancer that regulates the p53 pathway

Genetic analyses of lung cancer have helped found new treatments in this disease. We conducted an integrative analysis of gene expression and copy number in 261 non-small cell lung cancers (NSCLC) relative to matched normal tissues to define novel candidate oncogenes, identifying 12q13-15 and more specifically the YEATS4 gene as amplified and overexpressed in ~20% of the NSCLC cases examined. Overexpression of YEATS4 abrogated senescence in human bronchial epithelial cells. Conversely, RNAi-mediated attenuation of YEATS4 in human lung cancer cells reduced their proliferation and tumor growth, impairing colony formation and inducing cellular senescence. These effects were associated with increased levels of p21WAF1 and p53 and cleavage of PARP, implicating YEATS4 as a negative regulator of the p21-p53 pathway. We also found that YEATS4 expression affected cellular responses to cisplastin, with increased levels associated with resistance and decreased levels with sensitivity. Taken together, our findings reveal YEATS4 as a candidate oncogene amplified in NSCLC, and a novel mechanism contributing to NSCLC pathogenesis.

Systematic identification of molecular subtype-selective vulnerabilities in non-small-cell lung cancer

Context-specific molecular vulnerabilities that arise during tumor evolution represent an attractive intervention target class. However, the frequency and diversity of somatic lesions detected among lung tumors can confound efforts to identify these targets. To confront this challenge, we have applied parallel screening of chemical and genetic perturbations within a panel of molecularly annotated NSCLC lines to identify intervention opportunities tightly linked to molecular response indicators predictive of target sensitivity. Anchoring this analysis on a matched tumor/normal cell model from a lung adenocarcinoma patient identified three distinct target/response-indicator pairings that are represented with significant frequencies (6%-16%) in the patient population. These include NLRP3 mutation/inflammasome activation-dependent FLIP addiction, co-occurring KRAS and LKB1 mutation-driven COPI addiction, and selective sensitivity to a synthetic indolotriazine that is specified by a seven-gene expression signature. Target efficacies were validated in vivo, and mechanism-of-action studies informed generalizable principles underpinning cancer cell biology.

Elevated expression of BIRC6 protein in non-small-cell lung cancers is associated with cancer recurrence and chemoresistance

INTRODUCTION

Non-small-cell lung cancer (NSCLC) is an aggressive, highly chemoresistant disease. Reliable prognostic assays and more effective treatments are critically required. BIRC6 (baculoviral inhibitors of apoptosis proteins repeat-containing 6) protein is a member of the inhibitors of apoptosis protein family thought to play an important role in the progression or chemoresistance of many cancers. In this study, we investigated whether BIRC6 expression can be used as a prognostic marker or potential therapeutic target for NSCLC.

METHODS

In a retrospective analysis, BIRC6 protein expression was determined for 78 resected primary NSCLCs and nine benign lung tissues. Twenty-nine chemoresistant or chemosensitive subrenal capsule NSCLC tissue xenografts were assessed for BIRC6 expression, using immunohistochemistry, and 13 of them for BIRC6 gene copy number, using array comparative genomic hybridization analysis. The effect of small interfering RNA-induced BIRC6 knockdown on the growth of human NSCLC cell cultures and apoptosis (in combination with cisplatin) was investigated.

RESULTS

Elevated BIRC6 protein expression in NSCLC tissues was associated with poor 3-year relapse-free patient survival, lymph node involvement, and advanced pathological tumor, node, metastasis stage. In patient-derived lung squamous cell carcinoma xenografts, chemoresistance was associated with elevated BIRC6 expression and increased gene copy number. Small interfering RNA-induced BIRC6 down-regulation inhibited growth of the NSCLC cells and sensitized the cells to cisplatin.

CONCLUSIONS

BIRC6 may play an important role in the malignant progression and chemoresistance of NSCLC. Elevated BIRC6 protein expression may serve as a predictive marker for chemoresistance of NSCLCs and a poor prognostic factor for NSCLC patients. Down-regulation of the BIRC6 gene as a therapeutic approach may be effective, especially in combination with conventional chemotherapeutics.

Identification of TRPC6 as a possible candidate target gene within an amplicon at 11q21-q22.2 for migratory capacity in head and neck squamous cell carcinomas

BACKGROUND

Cytogenetic and gene expression analyses in head and neck squamous cell carcinomas (HNSCC) have allowed identification of genomic aberrations that may contribute to cancer pathophysiology. Nevertheless, the molecular consequences of numerous genetic alterations still remain unclear.

METHODS

To identify novel genes implicated in HNSCC pathogenesis, we analyzed the genomic alterations present in five HNSCC-derived cell lines by array CGH, and compared high level focal gene amplifications with gene expression levels to identify genes whose expression is directly impacted by these genetic events. Next, we knocked down TRPC6, one of the most highly amplified and over-expressed genes, to characterize the biological roles of TRPC6 in carcinogenesis. Finally, real time PCR was performed to determine TRPC6 gene dosage and mRNA levels in normal mucosa and human HNSCC tissues.

RESULTS

The data showed that the HNSCC-derived cell lines carry most of the recurrent genomic abnormalities previously described in primary tumors. High-level genomic amplifications were found at four chromosomal sites (11q21-q22.2, 18p11.31-p11.21, 19p13.2-p13.13, and 21q11) with associated gene expression changes in selective candidate genes suggesting that they may play an important role in the malignant behavior of HNSCC. One of the most dramatic alterations of gene transcription involved the TRPC6 gene (located at 11q21-q22.2) which has been recently implicated in tumour invasiveness. siRNA-induced knockdown of TRPC6 expression in HNSCC-derived cells dramatically inhibited HNSCC-cell invasion but did not significantly alter cell proliferation. Importantly, amplification and concomitant overexpression of TRPC6 was also found in HNSCC tumour samples.

CONCLUSIONS

Altogether, these data show that TRPC6 is likely to be a target for 11q21-22.2 amplification that confers enhanced invasive behavior to HNSCC cells. Therefore, TRPC6 may be a promising therapeutic target in the treatment of HNSCC.

Efficacy of two fluorescence in situ hybridization (FISH) probes for diagnosing malignant pleural effusions

It is difficult to differentiate tumor cells in pleural fluid from reactive benign mesothelium. Fluorescence in situ hybridization (FISH) can increase diagnostic accuracy. Two hundred pleural fluid samples were analyzed by using FISH probes for chromosomes 11 and 17. Histological analysis was used to diagnose cancer. Clinical, radiological, and histological data were used to exclude malignancy. Eighty-two pleural effusion samples had positive cytology, 51 were benign, and 67 were atypical, but inconclusive. The 82 positive cases were confirmed to be malignant. Among the 51 negative cytology cases, videothoracoscopy-guided pleural biopsy revealed malignancy in three; aneuploid cells were detected by FISH in all cases. In 43 of the 67 cases with inconclusive cytology, malignancy was confirmed based on histology and fluorescence in situ hybridization. One case of parapneumonic effusion with no evidence of cancer during clinical follow-up had a suspicious cytology and positive fluorescence in situ hybridization result. The remaining 23 cases had no histological, radiological, clinical, or genetic evidence of malignancy. This study demonstrated that cytogenetic analysis of fresh pleural fluid samples using only two FISH probes is a valuable ancillary method for the identification of malignant pleural effusion, particularly in cases in which oncotic cytology is inconclusive.

Identifying master regulators of cancer and their downstream targets by integrating genomic and epigenomic features

Vast amounts of molecular data characterizing the genome, epigenome and transcriptome are becoming available for a variety of cancers. The current challenge is to integrate these diverse layers of molecular biology information to create a more comprehensive view of key biological processes underlying cancer. We developed a biocomputational algorithm that integrates copy number, DNA methylation, and gene expression data to study master regulators of cancer and identify their targets. Our algorithm starts by generating a list of candidate driver genes based on the rationale that genes that are driven by multiple genomic events in a subset of samples are unlikely to be randomly deregulated. We then select the master regulators from the candidate driver and identify their targets by inferring the underlying regulatory network of gene expression. We applied our biocomputational algorithm to identify master regulators and their targets in glioblastoma multiforme (GBM) and serous ovarian cancer. Our results suggest that the expression of candidate drivers is more likely to be influenced by copy number variations than DNA methylation. Next, we selected the master regulators and identified their downstream targets using module networks analysis. As a proof-of-concept, we show that the GBM and ovarian cancer module networks recapitulate known processes in these cancers. In addition, we identify master regulators that have not been previously reported and suggest their likely role. In summary, focusing on genes whose expression can be explained by their genomic and epigenomic aberrations is a promising strategy to identify master regulators of cancer.

Cliques for the identification of gene signatures for colorectal cancer across population

Background

Colorectal cancer (CRC) is one of the most commonly diagnosed cancers worldwide. Studies have correlated risk of CRC development with dietary habits and environmental conditions. Gene signatures for any disease can identify the key biological processes, which is especially useful in studying cancer development. Such processes can be used to evaluate potential drug targets. Though recognition of CRC gene-signatures across populations is crucial to better understanding potential novel treatment options for CRC, it remains a challenging task.

Results

We developed a topological and biological feature-based network approach for identifying the gene signatures across populations. In this work, we propose a novel approach of using cliques to understand the variability within population. Cliques are more conserved and co-expressed, therefore allowing identification and comparison of cliques across a population which can help researchers study gene variations. Our study was based on four publicly available expression datasets belonging to four different populations across the world. We identified cliques of various sizes (0 to 7) across the four population networks. Cliques of size seven were further analyzed across populations for their commonality and uniqueness. Forty-nine common cliques of size seven were identified. These cliques were further analyzed based on their connectivity profiles. We found associations between the cliques and their connectivity profiles across networks. With these clique connectivity profiles (CCPs), we were able to identify the divergence among the populations, important biological processes (cell cycle, signal transduction, and cell differentiation), and related gene pathways. Therefore the genes identified in these cliques and their connectivity profiles can be defined as the gene-signatures across populations. In this work we demonstrate the power and effectiveness of cliques to study CRC across populations.

Conclusions

We developed a new approach where cliques and their connectivity profiles helped elucidate the variation and similarity in CRC gene profiles across four populations with unique dietary habits.