A scale-space method for detecting recurrent DNA copy number changes with analytical false discovery rate control

IDH1 mutation produces R-2-hydroxyglutarate (R-2HG) and induces mir-182-5p expression to regulate cell cycle and tumor formation in glioma

BACKGROUND

Mutations in isocitrate dehydrogenase 1 and 2 (IDH1 and IDH2), are present in most gliomas. IDH1 mutation is an important prognostic marker in glioma. However, its regulatory mechanism in glioma remains incompletely understood.

RESULTS

miR-182-5p expression was increased within IDH1-mutant glioma specimens according to TCGA, CGGA, and online dataset GSE119740, as well as collected clinical samples. (R)-2-hydroxyglutarate ((R)-2HG) treatment up-regulated the expression of miR-182-5p, enhanced glioma cell proliferation, and suppressed apoptosis; miR-182-5p inhibition partially eliminated the oncogenic effects of R-2HG upon glioma cells. By direct binding to Cyclin Dependent Kinase Inhibitor 2 C (CDKN2C) 3'UTR, miR-182-5p inhibited CDKN2C expression. Regarding cellular functions, CDKN2C knockdown promoted R-2HG-treated glioma cell viability, suppressed apoptosis, and relieved cell cycle arrest. Furthermore, CDKN2C knockdown partially attenuated the effects of miR-182-5p inhibition on cell phenotypes. Moreover, CDKN2C knockdown exerted opposite effects on cell cycle check point and apoptosis markers to those of miR-182-5p inhibition; also, CDKN2C knockdown partially attenuated the functions of miR-182-5p inhibition in cell cycle check point and apoptosis markers. The engineered CS-NPs (antagomir-182-5p) effectively encapsulated and delivered antagomir-182-5p, enhancing anti-tumor efficacy in vivo, indicating the therapeutic potential of CS-NPs(antagomir-182-5p) in targeting the miR-182-5p/CDKN2C axis against R-2HG-driven oncogenesis in mice models.

CONCLUSIONS

These insights highlight the potential of CS-NPs(antagomir-182-5p) to target the miR-182-5p/CDKN2C axis, offering a promising therapeutic avenue against R-2HG's oncogenic influence to glioma.

Large-scale Pan-cancer Cell Line Screening Identifies Actionable and Effective Drug Combinations

Oncology drug combinations can improve therapeutic responses and increase treatment options for patients. The number of possible combinations is vast and responses can be context-specific. Systematic screens can identify clinically relevant, actionable combinations in defined patient subtypes. We present data for 109 anticancer drug combinations from AstraZeneca's oncology small molecule portfolio screened in 755 pan-cancer cell lines. Combinations were screened in a 7 × 7 concentration matrix, with more than 4 million measurements of sensitivity, producing an exceptionally data-rich resource. We implement a new approach using combination Emax (viability effect) and highest single agent (HSA) to assess combination benefit. We designed a clinical translatability workflow to identify combinations with clearly defined patient populations, rationale for tolerability based on tumor type and combination-specific "emergent" biomarkers, and exposures relevant to clinical doses. We describe three actionable combinations in defined cancer types, confirmed in vitro and in vivo, with a focus on hematologic cancers and apoptotic targets.

SIGNIFICANCE

We present the largest cancer drug combination screen published to date with 7 × 7 concentration response matrices for 109 combinations in more than 750 cell lines, complemented by multi-omics predictors of response and identification of "emergent" combination biomarkers. We prioritize hits to optimize clinical translatability, and experimentally validate novel combination hypotheses. This article is featured in Selected Articles from This Issue, p. 695.

Comprehensive multi-omics analysis of the m7G in pan-cancer from the perspective of predictive, preventive, and personalized medicine

Background

The N7-methylguanosine modification (m7G) of the 5' cap structure in the mRNA plays a crucial role in gene expression. However, the relation between m7G and tumor immune remains unclear. Hence, we intended to perform a pan-cancer analysis of m7G which can help explore the underlying mechanism and contribute to predictive, preventive, and personalized medicine (PPPM / 3PM).

Methods

The gene expression, genetic variation, clinical information, methylation, and digital pathological section from 33 cancer types were downloaded from the TCGA database. Immunohistochemistry (IHC) was used to validate the expression of the m7G regulator genes (m7RGs) hub-gene. The m7G score was calculated by single-sample gene-set enrichment analysis. The association of m7RGs with copy number variation, clinical features, immune-related genes, TMB, MSI, and tumor immune dysfunction and exclusion (TIDE) was comprehensively assessed. CellProfiler was used to extract pathological section characteristics. XGBoost and random forest were used to construct the m7G score prediction model. Single-cell transcriptome sequencing (scRNA-seq) was used to assess the activation state of the m7G in the tumor microenvironment.

Results

The m7RGs were highly expressed in tumors and most of the m7RGs are risk factors for prognosis. Moreover, the cellular pathway enrichment analysis suggested that m7G score was closely associated with invasion, cell cycle, DNA damage, and repair. In several cancers, m7G score was significantly negatively correlated with MSI and TMB and positively correlated with TIDE, suggesting an ICB marker potential. XGBoost-based pathomics model accurately predicts m7G scores with an area under the ROC curve (AUC) of 0.97. Analysis of scRNA-seq suggests that m7G differs significantly among cells of the tumor microenvironment. IHC confirmed high expression of EIF4E in breast cancer. The m7G prognostic model can accurately assess the prognosis of tumor patients with an AUC of 0.81, which was publicly hosted at https://pan-cancer-m7g.shinyapps.io/Panca-m7g/.

Conclusion

The current study explored for the first time the m7G in pan-cancer and identified m7G as an innovative marker in predicting clinical outcomes and immunotherapeutic efficacy, with the potential for deeper integration with PPPM. Combining m7G within the framework of PPPM will provide a unique opportunity for clinical intelligence and new approaches.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-022-00305-1.

Computational estimation of quality and clinical relevance of cancer cell lines

Immortal cancer cell lines (CCLs) are the most widely used system for investigating cancer biology and for the preclinical development of oncology therapies. Pharmacogenomic and genome‐wide editing screenings have facilitated the discovery of clinically relevant gene–drug interactions and novel therapeutic targets via large panels of extensively characterised CCLs. However, tailoring pharmacological strategies in a precision medicine context requires bridging the existing gaps between tumours and in vitro models. Indeed, intrinsic limitations of CCLs such as misidentification, the absence of tumour microenvironment and genetic drift have highlighted the need to identify the most faithful CCLs for each primary tumour while addressing their heterogeneity, with the development of new models where necessary. Here, we discuss the most significant limitations of CCLs in representing patient features, and we review computational methods aiming at systematically evaluating the suitability of CCLs as tumour proxies and identifying the best patient representative in vitro models. Additionally, we provide an overview of the applications of these methods to more complex models and discuss future machine‐learning‐based directions that could resolve some of the arising discrepancies.

CELLector: Genomics-Guided Selection of Cancer In Vitro Models

Selecting appropriate cancer models is a key prerequisite for maximizing translational potential and clinical relevance of in vitro oncology studies. We developed CELLector: an R package and R Shiny application allowing researchers to select the most relevant cancer cell lines in a patient-genomic-guided fashion. CELLector leverages tumor genomics to identify recurrent subtypes with associated genomic signatures. It then evaluates these signatures in cancer cell lines to prioritize their selection. This enables users to choose appropriate in vitro models for inclusion or exclusion in retrospective analyses and future studies. Moreover, this allows bridging outcomes from cancer cell line screens to precisely defined sub-cohorts of primary tumors. Here, we demonstrate the usefulness and applicability of CELLector, showing how it can aid prioritization of in vitro models for future development and unveil patient-derived multivariate prognostic and therapeutic markers. CELLector is freely available at https://ot-cellector.shinyapps.io/CELLector_App/ (code at https://github.com/francescojm/CELLector and https://github.com/francescojm/CELLector_App).

Probability distribution of copy number alterations along the genome: an algorithm to distinguish different tumour profiles

Copy number alterations (CNAs) comprise deletions or amplifications of fragments of genomic material that are particularly common in cancer and play a major contribution in its development and progression. High resolution microarray-based genome-wide technologies have been widely used to detect CNAs, generating complex datasets that require further steps to allow for the determination of meaningful results. In this work, we propose a methodology to determine common regions of CNAs from these datasets, that in turn are used to infer the probability distribution of disease profiles in the population. This methodology was validated using simulated data and assessed using real data from Head and Neck Squamous Cell Carcinoma and Lung Adenocarcinoma, from the TCGA platform. Probability distribution profiles were produced allowing for the distinction between different phenotypic groups established within that cohort. This method may be used to distinguish between groups in the diseased population, within well-established degrees of confidence. The application of such methods may be of greater value in the clinical context both as a diagnostic or prognostic tool and, even as a useful way for helping to establish the most adequate treatment and care plans.

Identification of Intrinsic Drug Resistance and Its Biomarkers in High-Throughput Pharmacogenomic and CRISPR Screens

High-throughput drug screens in cancer cell lines test compounds at low concentrations, thereby enabling the identification of drug-sensitivity biomarkers, while resistance biomarkers remain underexplored. Dissecting meaningful drug responses at high concentrations is challenging due to cytotoxicity, i.e., off-target effects, thus limiting resistance biomarker discovery to frequently mutated cancer genes. To address this, we interrogate subpopulations carrying sensitivity biomarkers and consecutively investigate unexpectedly resistant (UNRES) cell lines for unique genetic alterations that may drive resistance. By analyzing the GDSC and CTRP datasets, we find 53 and 35 UNRES cases, respectively. For 24 and 28 of them, we highlight putative resistance biomarkers. We find clinically relevant cases such as EGFRT790M mutation in NCI-H1975 or PTEN loss in NCI-H1650 cells, in lung adenocarcinoma treated with EGFR inhibitors. Interrogating the underpinnings of drug resistance with publicly available CRISPR phenotypic assays assists in prioritizing resistance drivers, offering hypotheses for drug combinations.

Integrative analysis of genomic amplification-dependent expression and loss-of-function screen identifies ASAP1 as a driver gene in triple-negative breast cancer progression

The genetically heterogeneous triple-negative breast cancer (TNBC) continues to be an intractable disease, due to lack of effective targeted therapies. Gene amplification is a major event in tumorigenesis. Genes with amplification-dependent expression are being explored as therapeutic targets for cancer treatment. In this study, we have applied Analytical Multi-scale Identification of Recurring Events analysis and transcript quantification in the TNBC genome across 222 TNBC tumors and identified 138 candidate genes with positive correlation in copy number gain (CNG) and gene expression. siRNA-based loss-of-function screen of the candidate genes has validated EGFR, MYC, ASAP1, IRF2BP2, and CCT5 genes as drivers promoting proliferation in different TNBC cells. MYC, ASAP1, IRF2BP2, and CCT5 display frequent CNG and concurrent expression over 2173 breast cancer tumors (cBioPortal dataset). More frequently are MYC and ASAP1 amplified in TNBC tumors (>30%, n = 320). In particular, high expression of ASAP1, the ADP-ribosylation factor GTPase-activating protein, is significantly related to poor metastatic relapse-free survival of TNBC patients (n = 257, bc-GenExMiner). Furthermore, we have revealed that silencing of ASAP1 modulates numerous cytokine and apoptosis signaling components, such as IL1B, TRAF1, AIFM2, and MAP3K11 that are clinically relevant to survival outcomes of TNBC patients. ASAP1 has been reported to promote invasion and metastasis in various cancer cells. Our findings that ASAP1 is an amplification-dependent TNBC driver gene promoting TNBC cell proliferation, functioning upstream apoptosis components, and correlating to clinical outcomes of TNBC patients, support ASAP1 as a potential actionable target for TNBC treatment.

HetRCNA: A Novel Method to Identify Recurrent Copy Number Alternations from Heterogeneous Tumor Samples Based on Matrix Decomposition Framework

A common strategy to discovering cancer associated copy number aberrations (CNAs) from a cohort of cancer samples is to detect recurrent CNAs (RCNAs). Although the previous methods can successfully identify communal RCNAs shared by nearly all tumor samples, detecting subgroup-specific RCNAs and their related subgroup samples from cancer samples with heterogeneity is still invalid for these existing approaches. In this paper, we introduce a novel integrated method called HetRCNA, which can identify statistically significant subgroup-specific RCNAs and their related subgroup samples. Based on matrix decomposition framework with weight constraint, HetRCNA can successfully measure the subgroup samples by coefficients of left vectors with weight constraint and subgroup-specific RCNAs by coefficients of the right vectors and significance test. When we evaluate HetRCNA on simulated dataset, the results show that HetRCNA gives the best performances among the competing methods and is robust to the noise factors of the simulated data. When HetRCNA is applied on a real breast cancer dataset, our approach successfully identifies a bunch of RCNA regions and the result is highly correlated with the results of the other two investigated approaches. Notably, the genomic regions identified by HetRCNA harbor many breast cancer related genes reported by previous researches.

Gene networks in cancer are biased by aneuploidies and sample impurities

Gene regulatory network inference is a standard technique for obtaining structured regulatory information from, for instance, gene expression measurements. Methods performing this task have been extensively evaluated on synthetic, and to a lesser extent real data sets. In contrast to these test evaluations, applications to gene expression data of human cancers are often limited by fewer samples and more potential regulatory links, and are biased by copy number aberrations as well as cell mixtures and sample impurities. Here, we take networks inferred from TCGA cohorts as an example to show that (1) transcription factor annotations are essential to obtain reliable networks, and (2) even for state of the art methods, we expect that between 20 and 80% of edges are caused by copy number changes and cell mixtures rather than transcription factor regulation.

Comparative oncogenomics identifies combinations of driver genes and drug targets in BRCA1-mutated breast cancer

BRCA1-mutated breast cancer is primarily driven by DNA copy-number alterations (CNAs) containing large numbers of candidate driver genes. Validation of these candidates requires novel approaches for high-throughput in vivo perturbation of gene function. Here we develop genetically engineered mouse models (GEMMs) of BRCA1-deficient breast cancer that permit rapid introduction of putative drivers by either retargeting of GEMM-derived embryonic stem cells, lentivirus-mediated somatic overexpression or in situ CRISPR/Cas9-mediated gene disruption. We use these approaches to validate Myc, Met, Pten and Rb1 as bona fide drivers in BRCA1-associated mammary tumorigenesis. Iterative mouse modeling and comparative oncogenomics analysis show that MYC-overexpression strongly reshapes the CNA landscape of BRCA1-deficient mammary tumors and identify MCL1 as a collaborating driver in these tumors. Moreover, MCL1 inhibition potentiates the in vivo efficacy of PARP inhibition (PARPi), underscoring the therapeutic potential of this combination for treatment of BRCA1-mutated cancer patients with poor response to PARPi monotherapy.

Analysis of high-throughput biological data using their rank values

High-throughput biological technologies are routinely used to generate gene expression profiling or cytogenetics data. To achieve high performance, methods available in the literature become more specialized and often require high computational resources. Here, we propose a new versatile method based on the data-ordering rank values. We use linear algebra, the Perron-Frobenius theorem and also extend a method presented earlier for searching differentially expressed genes for the detection of recurrent copy number aberration. A result derived from the proposed method is a one-sample Student's t-test based on rank values. The proposed method is to our knowledge the only that applies to gene expression profiling and to cytogenetics data sets. This new method is fast, deterministic, and requires a low computational load. Probabilities are associated with genes to allow a statistically significant subset selection in the data set. Stability scores are also introduced as quality parameters. The performance and comparative analyses were carried out using real data sets. The proposed method can be accessed through an R package available from the CRAN (Comprehensive R Archive Network) website: https://cran.r-project.org/web/packages/fcros .

RUBIC identifies driver genes by detecting recurrent DNA copy number breaks

The frequent recurrence of copy number aberrations across tumour samples is a reliable hallmark of certain cancer driver genes. However, state-of-the-art algorithms for detecting recurrent aberrations fail to detect several known drivers. In this study, we propose RUBIC, an approach that detects recurrent copy number breaks, rather than recurrently amplified or deleted regions. This change of perspective allows for a simplified approach as recursive peak splitting procedures and repeated re-estimation of the background model are avoided. Furthermore, we control the false discovery rate on the level of called regions, rather than at the probe level, as in competing algorithms. We benchmark RUBIC against GISTIC2 (a state-of-the-art approach) and RAIG (a recently proposed approach) on simulated copy number data and on three SNP6 and NGS copy number data sets from TCGA. We show that RUBIC calls more focal recurrent regions and identifies a much larger fraction of known cancer genes.

Analysis of cancer genome sequencing data has been used to predict genes associated with the pathogenesis of cancer. Here, the authors propose a new algorithm entitled RUBIC that predicts breaks in DNA as opposed to previously published methods that predict amplifications and deletions of DNA.

A Landscape of Pharmacogenomic Interactions in Cancer

Systematic studies of cancer genomes have provided unprecedented insights into the molecular nature of cancer. Using this information to guide the development and application of therapies in the clinic is challenging. Here, we report how cancer-driven alterations identified in 11,289 tumors from 29 tissues (integrating somatic mutations, copy number alterations, DNA methylation, and gene expression) can be mapped onto 1,001 molecularly annotated human cancer cell lines and correlated with sensitivity to 265 drugs. We find that cell lines faithfully recapitulate oncogenic alterations identified in tumors, find that many of these associate with drug sensitivity/resistance, and highlight the importance of tissue lineage in mediating drug response. Logic-based modeling uncovers combinations of alterations that sensitize to drugs, while machine learning demonstrates the relative importance of different data types in predicting drug response. Our analysis and datasets are rich resources to link genotypes with cellular phenotypes and to identify therapeutic options for selected cancer sub-populations.

•

We integrate heterogeneous molecular data of 11,289 tumors and 1,001 cell lines

•

We measure the response of 1,001 cancer cell lines to 265 anti-cancer drugs

•

We uncover numerous oncogenic aberrations that sensitize to an anti-cancer drug

•

Our study forms a resource to identify therapeutic options for cancer sub-populations

Discovering Recurrent Copy Number Aberrations in Complex Patterns via Non-Negative Sparse Singular Value Decomposition

Recurrent copy number aberrations (RCNAs) in multiple cancer samples are strongly associated with tumorigenesis, and RCNA discovery is helpful to cancer research and treatment. Despite the emergence of numerous RCNA discovering methods, most of them are unable to detect RCNAs in complex patterns that are influenced by complicating factors including aberration in partial samples, co-existing of gains and losses and normal-like tumor samples. Here, we propose a novel computational method, called non-negative sparse singular value decomposition (NN-SSVD), to address the RCNA discovering problem in complex patterns. In NN-SSVD, the measurement of RCNA is based on the aberration frequency in a part of samples rather than all samples, which can circumvent the complexity of different RCNA patterns. We evaluate NN-SSVD on synthetic dataset by comparison on detection scores and Receiver Operating Characteristics curves, and the results show that NN-SSVD outperforms existing methods in RCNA discovery and demonstrate more robustness to RCNA complicating factors. Applying our approach on a breast cancer dataset, we successfully identify a number of genomic regions that are strongly correlated with previous studies, which harbor a bunch of known breast cancer associated genes.

Identification of cancer-driver genes in focal genomic alterations from whole genome sequencing data

DNA copy number alterations (CNAs) are the main genomic events that occur during the initiation and development of cancer. Distinguishing driver aberrant regions from passenger regions, which might contain candidate target genes for cancer therapies, is an important issue. Several methods for identifying cancer-driver genes from multiple cancer patients have been developed for single nucleotide polymorphism (SNP) arrays. However, for NGS data, methods for the SNP array cannot be directly applied because of different characteristics of NGS such as higher resolutions of data without predefined probes and incorrectly mapped reads to reference genomes. In this study, we developed a wavelet-based method for identification of focal genomic alterations for sequencing data (WIFA-Seq). We applied WIFA-Seq to whole genome sequencing data from glioblastoma multiforme, ovarian serous cystadenocarcinoma and lung adenocarcinoma, and identified focal genomic alterations, which contain candidate cancer-related genes as well as previously known cancer-driver genes.

Integration of genomic, transcriptomic and proteomic data identifies two biologically distinct subtypes of invasive lobular breast cancer

Invasive lobular carcinoma (ILC) is the second most frequently occurring histological breast cancer subtype after invasive ductal carcinoma (IDC), accounting for around 10% of all breast cancers. The molecular processes that drive the development of ILC are still largely unknown. We have performed a comprehensive genomic, transcriptomic and proteomic analysis of a large ILC patient cohort and present here an integrated molecular portrait of ILC. Mutations in CDH1 and in the PI3K pathway are the most frequent molecular alterations in ILC. We identified two main subtypes of ILCs: (i) an immune related subtype with mRNA up-regulation of PD-L1, PD-1 and CTLA-4 and greater sensitivity to DNA-damaging agents in representative cell line models; (ii) a hormone related subtype, associated with Epithelial to Mesenchymal Transition (EMT), and gain of chromosomes 1q and 8q and loss of chromosome 11q. Using the somatic mutation rate and eIF4B protein level, we identified three groups with different clinical outcomes, including a group with extremely good prognosis. We provide a comprehensive overview of the molecular alterations driving ILC and have explored links with therapy response. This molecular characterization may help to tailor treatment of ILC through the application of specific targeted, chemo- and/or immune-therapies.

Computational characterisation of cancer molecular profiles derived using next generation sequencing

Our current understanding of cancer genetics is grounded on the principle that cancer arises from a clone that has accumulated the requisite somatically acquired genetic aberrations, leading to the malignant transformation. It also results in aberrent of gene and protein expression. Next generation sequencing (NGS) or deep sequencing platforms are being used to create large catalogues of changes in copy numbers, mutations, structural variations, gene fusions, gene expression, and other types of information for cancer patients. However, inferring different types of biological changes from raw reads generated using the sequencing experiments is algorithmically and computationally challenging. In this article, we outline common steps for the quality control and processing of NGS data. We highlight the importance of accurate and application-specific alignment of these reads and the methodological steps and challenges in obtaining different types of information. We comment on the importance of integrating these data and building infrastructure to analyse it. We also provide exhaustive lists of available software to obtain information and point the readers to articles comparing software for deeper insight in specialised areas. We hope that the article will guide readers in choosing the right tools for analysing oncogenomic datasets.

Identifying driver mutations in sequenced cancer genomes: computational approaches to enable precision medicine

High-throughput DNA sequencing is revolutionizing the study of cancer and enabling the measurement of the somatic mutations that drive cancer development. However, the resulting sequencing datasets are large and complex, obscuring the clinically important mutations in a background of errors, noise, and random mutations. Here, we review computational approaches to identify somatic mutations in cancer genome sequences and to distinguish the driver mutations that are responsible for cancer from random, passenger mutations. First, we describe approaches to detect somatic mutations from high-throughput DNA sequencing data, particularly for tumor samples that comprise heterogeneous populations of cells. Next, we review computational approaches that aim to predict driver mutations according to their frequency of occurrence in a cohort of samples, or according to their predicted functional impact on protein sequence or structure. Finally, we review techniques to identify recurrent combinations of somatic mutations, including approaches that examine mutations in known pathways or protein-interaction networks, as well as de novo approaches that identify combinations of mutations according to statistical patterns of mutual exclusivity. These techniques, coupled with advances in high-throughput DNA sequencing, are enabling precision medicine approaches to the diagnosis and treatment of cancer.