OncodriveCLUSTL: a sequence-based clustering method to identify cancer drivers

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

Copy number losses of oncogenes and gains of tumor suppressor genes generate common driver mutations

Cancer driver genes can undergo positive selection for various types of genetic alterations, including gain-of-function or loss-of-function mutations and copy number alterations (CNA). We investigated the landscape of different types of alterations affecting driver genes in 17,644 cancer exomes and genomes. We find that oncogenes may simultaneously exhibit signatures of positive selection and also negative selection in different gene segments, suggesting a method to identify additional tumor types where an oncogene is a driver or a vulnerability. Next, we characterize the landscape of CNA-dependent selection effects, revealing a general trend of increased positive selection on oncogene mutations not only upon CNA gains but also upon CNA deletions. Similarly, we observe a positive interaction between mutations and CNA gains in tumor suppressor genes. Thus, two-hit events involving point mutations and CNA are universally observed regardless of the type of CNA and may signal new therapeutic opportunities. An analysis with focus on the somatic CNA two-hit events can help identify additional driver genes relevant to a tumor type. By a global inference of point mutation and CNA selection signatures and interactions thereof across genes and tissues, we identify 9 evolutionary archetypes of driver genes, representing different mechanisms of (in)activation by genetic alterations.

Assessing the validity of driver gene identification tools for targeted genome sequencing data

Motivation

Most cancer driver gene identification tools have been developed for whole-exome sequencing data. Targeted sequencing is a popular alternative to whole-exome sequencing for large cancer studies due to its greater depth at a lower cost per tumor. Unlike whole-exome sequencing, targeted sequencing only enables mutation calling for a selected subset of genes. Whether existing driver gene identification tools remain valid in that context has not previously been studied.

Results

We evaluated the validity of seven popular driver gene identification tools when applied to targeted sequencing data. Based on whole-exome data of 14 different cancer types from TCGA, we constructed matching targeted datasets by keeping only the mutations overlapping with the pan-cancer MSK-IMPACT panel and, in the case of breast cancer, also the breast-cancer-specific B-CAST panel. We then compared the driver gene predictions obtained on whole-exome and targeted mutation data for each of the seven tools. Differences in how the tools model background mutation rates were the most important determinant of their validity on targeted sequencing data. Based on our results, we recommend OncodriveFML, OncodriveCLUSTL, 20/20+, dNdSCv, and ActiveDriver for driver gene identification in targeted sequencing data, whereas MutSigCV and DriverML are best avoided in that context.

Availability and implementation

Code for the analyses is available at https://github.com/SchmidtGroupNKI/TGSdrivergene_validity.

Analysis of 3760 hematologic malignancies reveals rare transcriptomic aberrations of driver genes

BACKGROUND

Rare oncogenic driver events, particularly affecting the expression or splicing of driver genes, are suspected to substantially contribute to the large heterogeneity of hematologic malignancies. However, their identification remains challenging.

METHODS

To address this issue, we generated the largest dataset to date of matched whole genome sequencing and total RNA sequencing of hematologic malignancies from 3760 patients spanning 24 disease entities. Taking advantage of our dataset size, we focused on discovering rare regulatory aberrations. Therefore, we called expression and splicing outliers using an extension of the workflow DROP (Detection of RNA Outliers Pipeline) and AbSplice, a variant effect predictor that identifies genetic variants causing aberrant splicing. We next trained a machine learning model integrating these results to prioritize new candidate disease-specific driver genes.

RESULTS

We found a median of seven expression outlier genes, two splicing outlier genes, and two rare splice-affecting variants per sample. Each category showed significant enrichment for already well-characterized driver genes, with odds ratios exceeding three among genes called in more than five samples. On held-out data, our integrative modeling significantly outperformed modeling based solely on genomic data and revealed promising novel candidate driver genes. Remarkably, we found a truncated form of the low density lipoprotein receptor LRP1B transcript to be aberrantly overexpressed in about half of hairy cell leukemia variant (HCL-V) samples and, to a lesser extent, in closely related B-cell neoplasms. This observation, which was confirmed in an independent cohort, suggests LRP1B as a novel marker for a HCL-V subclass and a yet unreported functional role of LRP1B within these rare entities.

CONCLUSIONS

Altogether, our census of expression and splicing outliers for 24 hematologic malignancy entities and the companion computational workflow constitute unique resources to deepen our understanding of rare oncogenic events in hematologic cancers.

Gastric cancer actionable genomic alterations across diverse populations worldwide and pharmacogenomics strategies based on precision oncology

Introduction: Gastric cancer is one of the most prevalent types of cancer worldwide. The World Health Organization (WHO), the International Agency for Research on Cancer (IARC), and the Global Cancer Statistics (GLOBOCAN) reported an age standardized global incidence rate of 9.2 per 100,000 individuals for gastric cancer in 2022, with a mortality rate of 6.1. Despite considerable progress in precision oncology through the efforts of international consortia, understanding the genomic features and their influence on the effectiveness of anti-cancer treatments across diverse ethnic groups remains essential. Methods: Our study aimed to address this need by conducting integrated in silico analyses to identify actionable genomic alterations in gastric cancer driver genes, assess their impact using deleteriousness scores, and determine allele frequencies across nine global populations: European Finnish, European non-Finnish, Latino, East Asian, South Asian, African, Middle Eastern, Ashkenazi Jewish, and Amish. Furthermore, our goal was to prioritize targeted therapeutic strategies based on pharmacogenomics clinical guidelines, in silico drug prescriptions, and clinical trial data. Results: Our comprehensive analysis examined 275,634 variants within 60 gastric cancer driver genes from 730,947 exome sequences and 76,215 whole-genome sequences from unrelated individuals, identifying 13,542 annotated and predicted oncogenic variants. We prioritized the most prevalent and deleterious oncogenic variants for subsequent pharmacogenomics testing. Additionally, we discovered actionable genomic alterations in the ARID1A, ATM, BCOR, ERBB2, ERBB3, CDKN2A, KIT, PIK3CA, PTEN, NTRK3, TP53, and CDKN2A genes that could enhance the efficacy of anti-cancer therapies, as suggested by in silico drug prescription analyses, reviews of current pharmacogenomics clinical guidelines, and evaluations of phase III and IV clinical trials targeting gastric cancer driver proteins. Discussion: These findings underline the urgency of consolidating efforts to devise effective prevention measures, invest in genomic profiling for underrepresented populations, and ensure the inclusion of ethnic minorities in future clinical trials and cancer research in developed countries.

Uncovering the subtype-specific disease module and the development of drug response prediction models for glioma

The poor prognosis of glioma patients brought attention to the need for effective therapeutic approaches for precision therapy. Here, we deployed algorithms relying on network medicine and artificial intelligence to design the framework for subtype-specific target identification and drug response prediction in glioma. We identified the driver mutations that were differentially expressed in each subtype of lower-grade glioma and glioblastoma multiforme and were linked to cancer-specific processes. Driver mutations that were differentially expressed were also subjected to subtype-specific disease module identification. The drugs from the drug bank database were retrieved to target these disease modules. However, the efficacy of anticancer drugs depends on the molecular profile of the cancer and varies among cancer patients due to intratumor heterogeneity. Hence, we developed a deep-learning-based drug response prediction framework using the experimental drug screening data. Models for 30 drugs that can target the disease module were developed, where drug response measured by IC50 was considered a response and gene expression and mutation data were considered predictor variables. The model construction consists of three steps: feature selection, data integration, and classification. We observed the consistent performance of the models in training, test, and validation datasets. Drug responses were predicted for particular cell lines derived from distinct subtypes of gliomas. We found that subtypes of gliomas respond differently to the drug, highlighting the importance of subtype-specific drug response prediction. Therefore, the development of personalized therapy by integrating network medicine and a deep learning-based approach can lead to cancer-specific treatment and improved patient care.

Prediction of cancer driver genes and mutations: the potential of integrative computational frameworks

The vast amount of available sequencing data allows the scientific community to explore different genetic alterations that may drive cancer or favor cancer progression. Software developers have proposed a myriad of predictive tools, allowing researchers and clinicians to compare and prioritize driver genes and mutations and their relative pathogenicity. However, there is little consensus on the computational approach or a golden standard for comparison. Hence, benchmarking the different tools depends highly on the input data, indicating that overfitting is still a massive problem. One of the solutions is to limit the scope and usage of specific tools. However, such limitations force researchers to walk on a tightrope between creating and using high-quality tools for a specific purpose and describing the complex alterations driving cancer. While the knowledge of cancer development increases daily, many bioinformatic pipelines rely on single nucleotide variants or alterations in a vacuum without accounting for cellular compartments, mutational burden or disease progression. Even within bioinformatics and computational cancer biology, the research fields work in silos, risking overlooking potential synergies or breakthroughs. Here, we provide an overview of databases and datasets for building or testing predictive cancer driver tools. Furthermore, we introduce predictive tools for driver genes, driver mutations, and the impact of these based on structural analysis. Additionally, we suggest and recommend directions in the field to avoid silo-research, moving towards integrative frameworks.

Identification of Clonal Hematopoiesis Driver Mutations through In Silico Saturation Mutagenesis

Clonal hematopoiesis (CH) is a phenomenon of clonal expansion of hematopoietic stem cells driven by somatic mutations affecting certain genes. Recently, CH has been linked to the development of a number of hematologic malignancies, cardiovascular diseases and other conditions. Although the most frequently mutated CH driver genes have been identified, a systematic landscape of the mutations capable of initiating this phenomenon is still lacking. Here, we train high-quality machine-learning models for 12 of the most recurrent CH driver genes to identify their driver mutations. These models outperform an experimental base-editing approach and expert-curated rules based on prior knowledge of the function of these genes. Moreover, their application to identify CH driver mutations across almost half a million donors of the UK Biobank reproduces known associations between CH driver mutations and age, and the prevalence of several diseases and conditions. We thus propose that these models support the accurate identification of CH across healthy individuals.

Specific genetic aberrations of parathyroid in Chinese patients with tertiary hyperparathyroidism using whole-exome sequencing

Background

Tertiary hyperparathyroidism (THPT) is a peculiar subtype of hyperparathyroidism that usually develops from chronic kidney disease (CKD) and persists even after kidney transplantation. Unlike its precursor, secondary hyperparathyroidism (SHPT), THPT is characterized by uncontrolled high levels of calcium in the blood, which suggests the monoclonal or oligoclonal proliferation of parathyroid cells. However, the molecular abnormalities leading to THPT have not yet been fully understood.

Methods

In this study, we analyzed DNA samples from hyperplastic parathyroid and corresponding blood cells of 11 patients with THPT using whole-exome sequencing (WES). We identified somatic single nucleotide variants (SNV) and insertions or deletions variants (INDEL) and performed driver mutation analysis, KEGG pathway, and GO functional enrichment analysis. To confirm the impact of selected driver mutated genes, we also tested their expression level in these samples using qRT-PCR.

Results

Following quality control and mutation filtering, we identified 17,401 mutations, comprising 6690 missense variants, 3078 frameshift variants, 2005 stop-gained variants, and 1630 synonymous variants. Copy number variants (CNV) analysis showed that chromosome 22 copy number deletion was frequently observed in 6 samples. Driver mutation analysis identified 179 statistically significant mutated genes, including recurrent missense mutations on TBX20, ATAD5, ZNF669, and NOX3 genes in 3 different patients. KEGG pathway analysis revealed two enriched pathways: non-homologous end-joining and cell cycle, with a sole gene, PRKDC, involved. GO analysis demonstrated significant enrichment of various cellular components and cytobiological processes associated with four genes, including GO items of positive regulation of developmental growth, protein ubiquitination, and positive regulation of the apoptotic process. Compared to blood samples, THPT samples exhibited lower expression levels of PRKDC, TBX20, ATAD5, and NOX3 genes. THPT samples with exon mutations had relatively lower expression levels of PRKDC, TBX20, and NOX3 genes compared to those without mutations, although the difference was not statistically significant.

Conclusion

This study provides a comprehensive landscape of the genetic characteristics of hyperplastic parathyroids in THPT, highlighting the involvement of multiple genes and pathways in the development and progression of this disease. The dominant mutations identified in our study depicted new insights into the pathogenesis and molecular characteristics of THPT.

PhenoDriver: interpretable framework for studying personalized phenotype-associated driver genes in breast cancer

Identifying personalized cancer driver genes and further revealing their oncogenic mechanisms is critical for understanding the mechanisms of cell transformation and aiding clinical diagnosis. Almost all existing methods primarily focus on identifying driver genes at the cohort or individual level but fail to further uncover their underlying oncogenic mechanisms. To fill this gap, we present an interpretable framework, PhenoDriver, to identify personalized cancer driver genes, elucidate their roles in cancer development and uncover the association between driver genes and clinical phenotypic alterations. By analyzing 988 breast cancer patients, we demonstrate the outstanding performance of PhenoDriver in identifying breast cancer driver genes at the cohort level compared to other state-of-the-art methods. Otherwise, our PhenoDriver can also effectively identify driver genes with both recurrent and rare mutations in individual patients. We further explore and reveal the oncogenic mechanisms of some known and unknown breast cancer driver genes (e.g. TP53, MAP3K1, HTT, etc.) identified by PhenoDriver, and construct their subnetworks for regulating clinical abnormal phenotypes. Notably, most of our findings are consistent with existing biological knowledge. Based on the personalized driver profiles, we discover two existing and one unreported breast cancer subtypes and uncover their molecular mechanisms. These results intensify our understanding for breast cancer mechanisms, guide therapeutic decisions and assist in the development of targeted anticancer therapies.

Whole-exome sequencing of Indian prostate cancer reveals a novel therapeutic target: POLQ

PURPOSE

Prostate cancer is the second most common cancer diagnosed worldwide and the third most common cancer among men in India. This study's objective was to characterise the mutational landscape of Indian prostate cancer using whole-exome sequencing to identify population-specific polymorphisms.

METHODS

Whole-exome sequencing was performed of 58 treatment-naive primary prostate tumors of Indian origin. Multiple computational and statistical analyses were used to profile the known common mutations, other deleterious mutations, driver genes, prognostic biomarkers, and gene signatures unique to each clinical parameter. Cox analysis was performed to validate survival-associated genes. McNemar test identified genes significant to recurrence and receiver-operating characteristic (ROC) analysis was conducted to determine its accuracy. OncodriveCLUSTL algorithm was used to deduce driver genes. The druggable target identified was modeled with its known inhibitor using Autodock.

RESULTS

TP53 was the most commonly mutated gene in our cohort. Three novel deleterious variants unique to the Indian prostate cancer subtype were identified: POLQ, FTHL17, and OR8G1. COX regression analysis identified ACSM5, a mitochondrial gene responsible for survival. CYLC1 gene, which encodes for sperm head cytoskeletal protein, was identified as an unfavorable prognostic biomarker indicative of recurrence. The novel POLQ mutant, also identified as a driver gene, was evaluated as the druggable target in this study. POLQ, a DNA repair enzyme implicated in various cancer types, is overexpressed and is associated with a poor prognosis. The mutant POLQ was subjected to structural analysis and modeled with its known inhibitor novobiocin resulting in decreased binding efficiency necessitating the development of a better drug.

CONCLUSION

In this pilot study, the molecular profiling using multiple computational and statistical analyses revealed distinct polymorphisms in the Indian prostate cancer cohort. The mutational signatures identified provide a valuable resource for prognostic stratification and targeted treatment strategies for Indian prostate cancer patients. The DNA repair enzyme, POLQ, was identified as the druggable target in this study.

Multi-omics integration analysis of GPCRs in pan-cancer to uncover inter-omics relationships and potential driver genes

G protein-coupled receptors (GPCRs) are the largest drug target family. Unfortunately, applications of GPCRs in cancer therapy are scarce due to very limited knowledge regarding their correlations with cancers. Multi-omics data enables systematic investigations of GPCRs, yet their effective integration remains a challenge due to the complexity of the data. Here, we adopt two types of integration strategies, multi-staged and meta-dimensional approaches, to fully characterize somatic mutations, somatic copy number alterations (SCNAs), DNA methylations, and mRNA expressions of GPCRs in 33 cancers. Results from the multi-staged integration reveal that GPCR mutations cannot well predict expression dysregulation. The correlations between expressions and SCNAs are primarily positive, while correlations of the methylations with expressions and SCNAs are bimodal with negative correlations predominating. Based on these correlations, 32 and 144 potential cancer-related GPCRs driven by aberrant SCNA and methylation are identified, respectively. In addition, the meta-dimensional integration analysis is carried out by using deep learning models, which predict more than one hundred GPCRs as potential oncogenes. When comparing results between the two integration strategies, 165 cancer-related GPCRs are common in both, suggesting that they should be prioritized in future studies. However, 172 GPCRs emerge in only one, indicating that the two integration strategies should be considered concurrently to complement the information missed by the other such that obtain a more comprehensive understanding. Finally, correlation analysis further reveals that GPCRs, in particular for the class A and adhesion receptors, are generally immune-related. In a whole, the work is for the first time to reveal the associations between different omics layers and highlight the necessity of combing the two strategies in identifying cancer-related GPCRs.

Pan-cancer clinical impact of latent drivers from double mutations

Here, we discover potential 'latent driver' mutations in cancer genomes. Latent drivers have low frequencies and minor observable translational potential. As such, to date they have escaped identification. Their discovery is important, since when paired in cis, latent driver mutations can drive cancer. Our comprehensive statistical analysis of the pan-cancer mutation profiles of ~60,000 tumor sequences from the TCGA and AACR-GENIE cohorts identifies significantly co-occurring potential latent drivers. We observe 155 same gene double mutations of which 140 individual components are cataloged as latent drivers. Evaluation of cell lines and patient-derived xenograft response data to drug treatment indicate that in certain genes double mutations may have a prominent role in increasing oncogenic activity, hence obtaining a better drug response, as in PIK3CA. Taken together, our comprehensive analyses indicate that same-gene double mutations are exceedingly rare phenomena but are a signature for some cancer types, e.g., breast, and lung cancers. The relative rarity of doublets can be explained by the likelihood of strong signals resulting in oncogene-induced senescence, and by doublets consisting of non-identical single residue components populating the background mutational load, thus not identified.

Immunogenomic analysis of human brain metastases reveals diverse immune landscapes across genetically distinct tumors

Brain metastases (BrMs) are the most common form of brain tumors in adults and frequently originate from lung and breast primary cancers. BrMs are associated with high mortality, emphasizing the need for more effective therapies. Genetic profiling of primary tumors is increasingly used as part of the effort to guide targeted therapies against BrMs, and immune-based strategies for the treatment of metastatic cancer are gaining momentum. However, the tumor immune microenvironment (TIME) of BrM is extremely heterogeneous, and whether specific genetic profiles are associated with distinct immune states remains unknown. Here, we perform an extensive characterization of the immunogenomic landscape of human BrMs by combining whole-exome/whole-genome sequencing, RNA sequencing of immune cell populations, flow cytometry, immunofluorescence staining, and tissue imaging analyses. This revealed unique TIME phenotypes in genetically distinct lung- and breast-BrMs, thereby enabling the development of personalized immunotherapies tailored by the genetic makeup of the tumors.

Comparative Genomics Provides Etiologic and Biological Insight into Melanoma Subtypes

Melanoma is a cancer of melanocytes, with multiple subtypes based on body site location. Cutaneous melanoma is associated with skin exposed to ultraviolet radiation; uveal melanoma occurs in the eyes; mucosal melanoma occurs in internal mucous membranes; and acral melanoma occurs on the palms, soles, and nail beds. Here, we present the largest whole-genome sequencing study of melanoma to date, with 570 tumors profiled, as well as methylation and RNA sequencing for subsets of tumors. Uveal melanoma is genomically distinct from other melanoma subtypes, harboring the lowest tumor mutation burden and with significantly mutated genes in the G-protein signaling pathway. Most cutaneous, acral, and mucosal melanomas share alterations in components of the MAPK, PI3K, p53, p16, and telomere pathways. However, the mechanism by which these pathways are activated or inactivated varies between melanoma subtypes. Additionally, we identify potential novel germline predisposition genes for some of the less common melanoma subtypes.

SIGNIFICANCE

This is the largest whole-genome analysis of melanoma to date, comprehensively comparing the genomics of the four major melanoma subtypes. This study highlights both similarities and differences between the subtypes, providing insights into the etiology and biology of melanoma. This article is highlighted in the In This Issue feature, p. 2711.

Discovering the drivers of clonal hematopoiesis

Mutations in genes that confer a selective advantage to hematopoietic stem cells (HSCs) drive clonal hematopoiesis (CH). While some CH drivers have been identified, the compendium of all genes able to drive CH upon mutations in HSCs remains incomplete. Exploiting signals of positive selection in blood somatic mutations may be an effective way to identify CH driver genes, analogously to cancer. Using the tumor sample in blood/tumor pairs as reference, we identify blood somatic mutations across more than 12,000 donors from two large cancer genomics cohorts. The application of IntOGen, a driver discovery pipeline, to both cohorts, and more than 24,000 targeted sequenced samples yields a list of close to 70 genes with signals of positive selection in CH, available at http://www.intogen.org/ch. This approach recovers known CH genes, and discovers other candidates.

Identifying the genetic drivers of clonal haematopoiesis (CH) has been challenging due to their low frequencies and a lack of adequate tools. Here, the authors use a reverse calling to detect blood somatic mutations and the IntOGen pipeline to identify CH drivers in large cancer genomics data sets based on signals of positive selection.

Novel Calcium-Binding Ablating Mutations Induce Constitutive RET Activity and Drive Tumorigenesis

Distinguishing oncogenic mutations from variants of unknown significance (VUS) is critical for precision cancer medicine. Here, computational modeling of 71,756 RET variants for positive selection together with functional assays of 110 representative variants identified a three-dimensional cluster of VUSs carried by multiple human cancers that cause amino acid substitutions in the calmodulin-like motif (CaLM) of RET. Molecular dynamics simulations indicated that CaLM mutations decrease interactions between Ca2+ and its surrounding residues and induce conformational distortion of the RET cysteine-rich domain containing the CaLM. RET-CaLM mutations caused ligand-independent constitutive activation of RET kinase by homodimerization mediated by illegitimate disulfide bond formation. RET-CaLM mutants possessed oncogenic and tumorigenic activities that could be suppressed by tyrosine kinase inhibitors targeting RET. This study identifies calcium-binding ablating mutations as a novel type of oncogenic mutation of RET and indicates that in silico-driven annotation of VUSs of druggable oncogenes is a promising strategy to identify targetable driver mutations.

SIGNIFICANCE

Comprehensive proteogenomic and in silico analyses of a vast number of VUSs identify a novel set of oncogenic and druggable mutations in the well-characterized RET oncogene.

Hepatocellular Carcinoma in Mongolia Delineates Unique Molecular Traits and a Mutational Signature Associated with Environmental Agents

PURPOSE

Mongolia has the world's highest incidence of hepatocellular carcinoma (HCC), with ∼100 cases/100,000 inhabitants, although the reasons for this have not been thoroughly delineated.

EXPERIMENTAL DESIGN

We performed a molecular characterization of Mongolian (n = 192) compared with Western (n = 187) HCCs by RNA sequencing and whole-exome sequencing to unveil distinct genomic and transcriptomic features associated with environmental factors in this population.

RESULTS

Mongolian patients were younger, with higher female prevalence, and with predominantly HBV-HDV coinfection etiology. Mongolian HCCs presented significantly higher rates of protein-coding mutations (121 vs. 70 mutations per tumor in Western), and in specific driver HCC genes (i.e., APOB and TSC2). Four mutational signatures characterized Mongolian samples, one of which was novel (SBS Mongolia) and present in 25% of Mongolian HCC cases. This signature showed a distinct substitution profile with a high proportion of T>G substitutions and was significantly associated with a signature of exposure to the environmental agent dimethyl sulfate (71%), a 2A carcinogenic associated with coal combustion. Transcriptomic-based analysis delineated three molecular clusters, two not present in Western HCC; one with a highly inflamed profile and the other significantly associated with younger female patients.

CONCLUSIONS

Mongolian HCC has unique molecular traits with a high mutational burden and a novel mutational signature associated with genotoxic environmental factors present in this country.

Identification of Candidate Genes in Breast Cancer Induced by Estrogen Plus Progestogens Using Bioinformatic Analysis

Menopausal hormone therapy (MHT) was widely used to treat menopause-related symptoms in menopausal women. However, MHT therapies were controversial with the increased risk of breast cancer because of different estrogen and progestogen combinations, and the molecular basis behind this phenomenon is currently not understood. To address this issue, we identified differentially expressed genes (DEGs) between the estrogen plus progestogens treatment (EPT) and estrogen treatment (ET) using the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) data. As a result, a total of 96 upregulated DEGs were first identified. Seven DEGs related to the cell cycle (CCNE2, CDCA5, RAD51, TCF19, KNTC1, MCM10, and NEIL3) were validated by RT-qPCR. Specifically, these seven DEGs were increased in EPT compared to ET (p < 0.05) and had higher expression levels in breast cancer than adjacent normal tissues (p < 0.05). Next, we found that estrogen receptor (ER)-positive breast cancer patients with a higher CNNE2 expression have a shorter overall survival time (p < 0.05), while this effect was not observed in the other six DEGs (p > 0.05). Interestingly, the molecular docking results showed that CCNE2 might bind to 17β-estradiol (−6.791 kcal/mol), progesterone (−6.847 kcal/mol), and medroxyprogesterone acetate (−6.314 kcal/mol) with a relatively strong binding affinity, respectively. Importantly, CNNE2 protein level could be upregulated with EPT and attenuated by estrogen receptor antagonist, acolbifene and had interactions with cancer driver genes (AKT1 and KRAS) and high mutation frequency gene (TP53 and PTEN) in breast cancer patients. In conclusion, the current study showed that CCNE2, CDCA5, RAD51, TCF19, KNTC1, MCM10, and NEIL3 might contribute to EPT-related tumorigenesis in breast cancer, with CCNE2 might be a sensitive risk indicator of breast cancer risk in women using MHT.

Improvements in Quality Control and Library Preparation for Targeted Sequencing Allowed Detection of Potentially Pathogenic Alterations in Circulating Cell-Free DNA Derived from Plasma of Brain Tumor Patients

Malignant gliomas are the most frequent primary brain tumors in adults. They are genetically heterogenous and invariably recur due to incomplete surgery and therapy resistance. Circulating tumor DNA (ctDNA) is a component of circulating cell-free DNA (ccfDNA) and represents genetic material that originates from the primary tumor or metastasis. Brain tumors are frequently located in the eloquent brain regions, which makes biopsy difficult or impossible due to severe postoperative complications. The analysis of ccfDNA from a patient's blood presents a plausible and noninvasive alternative. In this study, freshly frozen tumors and corresponding blood samples were collected from 84 brain tumor patients and analyzed by targeted next-generation sequencing (NGS). The cohort included 80 glioma patients, 2 metastatic cancer patients, and 2 primary CNS lymphoma (PCNSL) patients. We compared the pattern of genetic alterations in the tumor DNA (tDNA) with that of ccfDNA. The implemented technical improvements in quality control and library preparation allowed for the detection of ctDNA in 8 out of 84 patients, including 5 out of 80 glioma patients. In 32 out of 84 patients, we found potentially pathogenic genetic alterations in ccfDNA that were not detectable in tDNA. While sequencing ccfDNA from plasma has a low efficacy as a diagnostic tool for glioma patients, we concluded that further improvements in sample processing and library preparation can make liquid biopsy a valuable diagnostic tool for glioma patients.

Novel Driver Strength Index highlights important cancer genes in TCGA PanCanAtlas patients

Background

Cancer driver genes are usually ranked by mutation frequency, which does not necessarily reflect their driver strength. We hypothesize that driver strength is higher for genes preferentially mutated in patients with few driver mutations overall, because these few mutations should be strong enough to initiate cancer.

Methods

We propose formulas for the Driver Strength Index (DSI) and the Normalized Driver Strength Index (NDSI), the latter independent of gene mutation frequency. We validate them using TCGA PanCanAtlas datasets, established driver prediction algorithms and custom computational pipelines integrating SNA, CNA and aneuploidy driver contributions at the patient-level resolution.

Results

DSI and especially NDSI provide substantially different gene rankings compared to the frequency approach. E.g., NDSI prioritized members of specific protein families, including G proteins GNAQ, GNA11 and GNAS, isocitrate dehydrogenases IDH1 and IDH2, and fibroblast growth factor receptors FGFR2 and FGFR3. KEGG analysis shows that top NDSI-ranked genes comprise EGFR/FGFR2/GNAQ/GNA11-NRAS/HRAS/KRAS-BRAF pathway, AKT1-MTOR pathway, and TCEB1-VHL-HIF1A pathway.

Conclusion

Our indices are able to select for driver gene attributes not selected by frequency sorting, potentially for driver strength. Genes and pathways prioritized are likely the strongest contributors to cancer initiation and progression and should become future therapeutic targets.

Detection of oncogenic and clinically actionable mutations in cancer genomes critically depends on variant calling tools

Motivation

The analysis of cancer genomes provides fundamental information about its etiology, the processes driving cell transformation or potential treatments. While researchers and clinicians are often only interested in the identification of oncogenic mutations, actionable variants or mutational signatures, the first crucial step in the analysis of any tumor genome is the identification of somatic variants in cancer cells (i.e. those that have been acquired during their evolution). For that purpose, a wide range of computational tools have been developed in recent years to detect somatic mutations in sequencing data from tumor samples. While there have been some efforts to benchmark somatic variant calling tools and strategies, the extent to which variant calling decisions impact the results of downstream analyses of tumor genomes remains unknown.

Results

Here, we quantify the impact of variant calling decisions by comparing the results obtained in three important analyses of cancer genomics data (identification of cancer driver genes, quantification of mutational signatures and detection of clinically actionable variants) when changing the somatic variant caller (MuSE, MuTect2, SomaticSniper and VarScan2) or the strategy to combine them (Consensus of two, Consensus of three and Union) across all 33 cancer types from The Cancer Genome Atlas. Our results show that variant calling decisions have a significant impact on these analyses, creating important differences that could even impact treatment decisions for some patients. Moreover, the Consensus of three calling strategy to combine the output of multiple variant calling tools, a very widely used strategy by the research community, can lead to the loss of some cancer driver genes and actionable mutations. Overall, our results highlight the limitations of widespread practices within the cancer genomics community and point to important differences in critical analyses of tumor sequencing data depending on variant calling, affecting even the identification of clinically actionable variants.

Availability and implementation

Code is available at https://github.com/carlosgarciaprieto/VariantCallingClinicalBenchmark.

Supplementary information

Supplementary data are available at Bioinformatics online.

Integrative Genomic Analyses of 1,145 Patient Samples Reveal New Biomarkers in Esophageal Squamous Cell Carcinoma

Due to the lack of effective diagnostic markers and therapeutic targets, esophageal squamous cell carcinoma (ESCC) shows a poor 5 years survival rate of less than 30%. To explore the potential therapeutic targets of ESCC, we integrated and reanalyzed the mutation data of WGS (whole genome sequencing) or WES (whole exome sequencing) from a total of 1,145 samples in 7 large ESCC cohorts, including 270 ESCC gene expression data. Two new mutation signatures and 20 driver genes were identified in our study. Among them, AP3S1, MUC16, and RPS15 were reported for the first time. We also discovered that the KMT2D was associated with the multiple clinical characteristics of ESCC, and KMT2D knockdown cells showed enhanced cell migration and cell invasion. Furthermore, a few neoantigens were shared between ESCC patients. For ESCC, compared to TMB, neoantigen might be treated as a better immunotherapy biomarker. Our research expands the understanding of ESCC mutations and helps the identification of ESCC biomarkers, especially for immunotherapy biomarkers.

Comprehensive patient-level classification and quantification of driver events in TCGA PanCanAtlas cohorts

There is a growing need to develop novel therapeutics for targeted treatment of cancer. The prerequisite to success is the knowledge about which types of molecular alterations are predominantly driving tumorigenesis. To shed light onto this subject, we have utilized the largest database of human cancer mutations–TCGA PanCanAtlas, multiple established algorithms for cancer driver prediction (2020plus, CHASMplus, CompositeDriver, dNdScv, DriverNet, HotMAPS, OncodriveCLUSTL, OncodriveFML) and developed four novel computational pipelines: SNADRIF (Single Nucleotide Alteration DRIver Finder), GECNAV (Gene Expression-based Copy Number Alteration Validator), ANDRIF (ANeuploidy DRIver Finder) and PALDRIC (PAtient-Level DRIver Classifier). A unified workflow integrating all these pipelines, algorithms and datasets at cohort and patient levels was created. We have found that there are on average 12 driver events per tumour, of which 0.6 are single nucleotide alterations (SNAs) in oncogenes, 1.5 are amplifications of oncogenes, 1.2 are SNAs in tumour suppressors, 2.1 are deletions of tumour suppressors, 1.5 are driver chromosome losses, 1 is a driver chromosome gain, 2 are driver chromosome arm losses, and 1.5 are driver chromosome arm gains. The average number of driver events per tumour increases with age (from 7 to 15) and cancer stage (from 10 to 15) and varies strongly between cancer types (from 1 to 24). Patients with 1 and 7 driver events per tumour are the most frequent, and there are very few patients with more than 40 events. In tumours having only one driver event, this event is most often an SNA in an oncogene. However, with increasing number of driver events per tumour, the contribution of SNAs decreases, whereas the contribution of copy-number alterations and aneuploidy events increases.

By analysing genomic and transcriptomic data from 10000 cancer patients through our custom-built computational pipelines and previously established third-party algorithms, we have found that half of all driver events in a patient’s tumour appear to be gains and losses of chromosomal arms and whole chromosomes. We therefore suggest that future therapeutics development efforts should be focused on targeting aneuploidy. We have also found that approximately a third of driver events in a patient are whole gene amplifications and deletions. Thus, therapies aimed at copy-number alterations also appear very promising. On the other hand, drugs aiming at point mutations are predicted to be less successful, as these alterations are responsible for just a couple of drivers per tumour. One notable exception are patients having only one driver event in their tumours, as this event is almost always a single nucleotide alteration in an oncogene.

Multidimensional Mutational Profiling of the Indian HNSCC Sub-Population Provides IRAK1, a Novel Driver Gene and Potential Druggable Target

Head and neck squamous cell carcinomas (HNSCC) include heterogeneous group of tumors, classified according to their anatomical site. It is the sixth most prevalent cancer globally. Among South Asian countries, India accounts for 40% of HNC malignancies with significant morbidity and mortality. In the present study, we have performed exome sequencing and analysis of 51 Head and Neck squamous cell carcinoma samples. Besides known mutations in the oncogenes and tumour suppressors, we have identified novel gene signatures differentiating buccal, alveolar, and tongue cancers. Around 50% of the patients showed mutation in tumour suppressor genes TP53 and TP63. Apart from the known mutations, we report novel mutations in the genes AKT1, SPECC1, and LRP1B, which are linked with tumour progression and patient survival. A highly curated process was developed to identify survival signatures. 36 survival-related genes were identified based on the correlation of functional impact of variants identified using exome-seq with gene expression from transcriptome data (GEPIA database) and survival. An independent LASSO regression analysis was also performed. Survival signatures common to both the methods led to identification of 4 dead and 3 alive gene signatures, the accuracy of which was confirmed by performing a ROC analysis (AUC=0.79 and 0.91, respectively). Also, machine learning-based driver gene prediction tool resulted in the identification of IRAK1 as the driver (p-value = 9.7 e-08) and also as an actionable mutation. Modelling of the IRAK1 mutation showed a decrease in its binding to known IRAK1 inhibitors.

Genomic Analysis Reveals Heterogeneity Between Lesions in Synchronous Primary Right-Sided and Left-Sided Colon Cancer

Background: The synchronous primary right-sided and left-sided colon cancer (sRL-CC) is a peculiar subtype of colorectal cancer. However, the genomic landscape of sRL-CC remains elusive. Methods: Twenty-eight paired tumor samples and their corresponding normal mucosa samples from 14 patients were collected from the Second Affiliated Hospital of Harbin Medical University from 2011 to 2018. The clinical-pathological data were obtained, and whole-exome sequencing was performed based on formalin-fixed and paraffin-embedded samples of these patients, and then, comprehensive bioinformatic analyses were conducted. Results: Both the lesions of sRL-CC presented dissimilar histological grade and differentiation. Based on sequencing data, few overlapping SNV signatures, onco-driver gene mutations, and SMGs were identified. Moreover, the paired lesions harbored a different distribution of copy number variants (CNVs) and loss of heterozygosity. The clonal architecture analysis demonstrated the polyclonal origin of sRL-CC and inter-cancerous heterogeneity between two lesions. Conclusion: Our work provides evidence that lesions of sRL-CC share few overlapping mutational signatures and CNVs, and may originate from different clones.

A New Insight for the Identification of Oncogenic Variants in Breast and Prostate Cancers in Diverse Human Populations, With a Focus on Latinos

Background: Breast cancer (BRCA) and prostate cancer (PRCA) are the most commonly diagnosed cancer types in Latin American women and men, respectively. Although in recent years large-scale efforts from international consortia have focused on improving precision oncology, a better understanding of genomic features of BRCA and PRCA in developing regions and racial/ethnic minority populations is still required.

Methods: To fill in this gap, we performed integrated in silico analyses to elucidate oncogenic variants from BRCA and PRCA driver genes; to calculate their deleteriousness scores and allele frequencies from seven human populations worldwide, including Latinos; and to propose the most effective therapeutic strategies based on precision oncology.

Results: We analyzed 339,100 variants belonging to 99 BRCA and 82 PRCA driver genes and identified 18,512 and 15,648 known/predicted oncogenic variants, respectively. Regarding known oncogenic variants, we prioritized the most frequent and deleterious variants of BRCA (n = 230) and PRCA (n = 167) from Latino, African, Ashkenazi Jewish, East Asian, South Asian, European Finnish, and European non-Finnish populations, to incorporate them into pharmacogenomics testing. Lastly, we identified which oncogenic variants may shape the response to anti-cancer therapies, detailing the current status of pharmacogenomics guidelines and clinical trials involved in BRCA and PRCA cancer driver proteins.

Conclusion: It is imperative to unify efforts where developing countries might invest in obtaining databases of genomic profiles of their populations, and developed countries might incorporate racial/ethnic minority populations in future clinical trials and cancer researches with the overall objective of fomenting pharmacogenomics in clinical practice and public health policies.

Improving existing analysis pipeline to identify and analyze cancer driver genes using multi-omics data

The cumulative of genes carrying mutations is vital for the establishment and development of cancer. However, this driver gene exploring research line has selected and used types of tools and models of analysis unsystematically and discretely. Also, the previous studies may have neglected low-frequency drivers and seldom predicted subgroup specificities of identified driver genes. In this study, we presented an improved driver gene identification and analysis pipeline that comprises the four most widely focused analyses for driver genes: enrichment analysis, clinical feature association with expression profiles of identified driver genes as well as with their functional modules, and patient stratification by existing advanced computational tools integrating multi-omics data. The improved pipeline's general usability was demonstrated straightforwardly for breast cancer, validated by some independent databases. Accordingly, 31 validated driver genes, including four novel ones, were discovered. Subsequently, we detected cancer-related significantly enriched gene ontology terms and pathways, probable drug targets, two co-expressed modules associated significantly with several clinical features, such as number of positive lymph nodes, Nottingham prognostic index, and tumor stage, and two biologically distinct groups of BRCA patients. Data and source code of the case study can be downloaded at https://github.com/hauldhut/drivergene.

A compendium of mutational cancer driver genes

A fundamental goal in cancer research is to understand the mechanisms of cell transformation. This is key to developing more efficient cancer detection methods and therapeutic approaches. One milestone towards this objective is the identification of all the genes with mutations capable of driving tumours. Since the 1970s, the list of cancer genes has been growing steadily. Because cancer driver genes are under positive selection in tumorigenesis, their observed patterns of somatic mutations across tumours in a cohort deviate from those expected from neutral mutagenesis. These deviations, which constitute signals of positive selection, may be detected by carefully designed bioinformatics methods, which have become the state of the art in the identification of driver genes. A systematic approach combining several of these signals could lead to a compendium of mutational cancer genes. In this Review, we present the Integrative OncoGenomics (IntOGen) pipeline, an implementation of such an approach to obtain the compendium of mutational cancer drivers. Its application to somatic mutations of more than 28,000 tumours of 66 cancer types reveals 568 cancer genes and points towards their mechanisms of tumorigenesis. The application of this approach to the ever-growing datasets of somatic tumour mutations will support the continuous refinement of our knowledge of the genetic basis of cancer.

HotSpotAnnotations-a database for hotspot mutations and annotations in cancer

Hotspots, recurrently mutated DNA positions in cancer, are thought to be oncogenic drivers because random chance is unlikely and the knowledge of clear examples of oncogenic hotspots in genes like BRAF, IDH1, KRAS and NRAS among many other genes. Hotspots are attractive because provide opportunities for biomedical research and novel treatments. Nevertheless, recent evidence, such as DNA hairpins for APOBEC3A, suggests that a considerable fraction of hotspots seem to be passengers rather than drivers. To document hotspots, the database HotSpotsAnnotations is proposed. For this, a statistical model was implemented to detect putative hotspots, which was applied to TCGA cancer datasets covering 33 cancer types, 10 182 patients and 3 175 929 mutations. Then, genes and hotspots were annotated by two published methods (APOBEC3A hairpins and dN/dS ratio) that may inform and warn researchers about possible false functional hotspots. Moreover, manual annotation from users can be added and shared. From the 23 198 detected as possible hotspots, 4435 were selected after false discovery rate correction and minimum mutation count. From these, 305 were annotated as likely for APOBEC3A whereas 442 were annotated as unlikely. To date, this is the first database dedicated to annotating hotspots for possible false functional hotspots.