In silico saturation mutagenesis of cancer genes

Protein structural context of cancer mutations reveals molecular mechanisms and candidate driver genes

Advances in protein structure determination and modeling allow us to study the structural context of human genetic variants on an unprecedented scale. Here, we analyze millions of cancer-associated missense mutations based on their structural locations and predicted perturbative effects. By considering the collective properties of mutations at the level of individual proteins, we identify distinct patterns associated with tumor suppressors and oncogenes. Tumor suppressors are enriched in structurally damaging mutations, consistent with loss-of-function mechanisms, while oncogene mutations tend to be structurally mild, reflecting selection for gain-of-function driver mutations and against loss-of-function mutations. Although oncogenes are difficult to distinguish from genes with no role in cancer using only structural damage, we find that the three-dimensional clustering of mutations is highly predictive. These observations allow us to identify candidate driver genes and speculate about their molecular roles, which we expect will have general utility in the analysis of cancer sequencing data.

Worldwide analysis of actionable genomic alterations in lung cancer and targeted pharmacogenomic strategies

Based on data from the Global Cancer Statistics 2022, lung cancer stands as the most lethal cancer worldwide, with age-adjusted incidence and mortality rates of 23.6 and 16.9 per 100,000 people, respectively. Despite significant strides in precision oncology driven by large-scale international research consortia, there remains a critical need to deepen our understanding of the genomic landscape across diverse racial and ethnic groups. To address this challenge, we performed comprehensive in silico analyses and data mining to identify pathogenic variants in genes that drive lung cancer. We subsequently calculated the allele frequencies and assessed the deleteriousness of these oncogenic variants among populations such as African, Amish, Ashkenazi Jewish, East and South Asian, Finnish and non-Finnish European, Latino, and Middle Eastern. Our analysis examined 117,707 variants within 86 lung cancer-associated genes across 75,109 human genomes, uncovering 8042 variants that are known or predicted to be pathogenic. We prioritized variants based on their allele frequencies and deleterious scores, and identified those with potential significance for response to anti-cancer therapies through in silico drug simulations, current clinical pharmacogenomic guidelines, and ongoing late-stage clinical trials targeting lung cancer-driving proteins. In conclusion, it is crucial to unite global efforts to create public health policies that emphasize prevention strategies and ensure access to clinical trials, pharmacogenomic testing, and cancer research for these groups in developed nations.

CDMPred: a tool for predicting cancer driver missense mutations with high-quality passenger mutations

Most computational methods for predicting driver mutations have been trained using positive samples, while negative samples are typically derived from statistical methods or putative samples. The representativeness of these negative samples in capturing the diversity of passenger mutations remains to be determined. To tackle these issues, we curated a balanced dataset comprising driver mutations sourced from the COSMIC database and high-quality passenger mutations obtained from the Cancer Passenger Mutation database. Subsequently, we encoded the distinctive features of these mutations. Utilizing feature correlation analysis, we developed a cancer driver missense mutation predictor called CDMPred employing feature selection through the ensemble learning technique XGBoost. The proposed CDMPred method, utilizing the top 10 features and XGBoost, achieved an area under the receiver operating characteristic curve (AUC) value of 0.83 and 0.80 on the training and independent test sets, respectively. Furthermore, CDMPred demonstrated superior performance compared to existing state-of-the-art methods for cancer-specific and general diseases, as measured by AUC and area under the precision-recall curve. Including high-quality passenger mutations in the training data proves advantageous for CDMPred's prediction performance. We anticipate that CDMPred will be a valuable tool for predicting cancer driver mutations, furthering our understanding of personalized therapy.

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 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. In this study, we trained machine learning models for 12 of the most recurrent CH genes to identify their driver mutations. These models outperform 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. Significance: We developed and validated gene-specific machine learning models to identify CH driver mutations, showing their advantage with respect to expert-curated rules. These models can support the identification and clinical interpretation of CH mutations in newly sequenced individuals. See related commentary by Arends and Jaiswal, p. 1581.

Characterizing the Mutational Landscape of Diffuse Large B-Cell Lymphoma in a Prospective Cohort of Mexican Patients

Diffuse large B-cell lymphoma (DLBCL) is the most common B-cell malignancy worldwide. Molecular classifications have tried to improve cure rates. We prospectively examined and correlated the mutational landscape with the clinical features and outcomes of 185 Mexican patients (median age 59.3 years, 50% women) with newly diagnosed DLBCL. A customized panel of 79 genes was designed, based on previous international series. Most patients had ECOG performance status (PS) < 2 (69.2%), advanced-stage disease (72.4%), germinal-center phenotype (68.1%), and double-hit lymphomas (14.1%). One hundred and ten (59.5%) patients had at least one gene with driver mutations. The most common mutated genes were as follows: TP53, EZH2, CREBBP, NOTCH1, and KMT2D. The median follow-up was 42 months, and the 5-year relapse-free survival (RFS) and overall survival (OS) rates were 70% and 72%, respectively. In the multivariate analysis, both age > 50 years and ECOG PS > 2 were significantly associated with a worse OS. Our investigation did not reveal any discernible correlation between the presence of a specific mutation and survival. In conclusion, using a customized panel, we characterized the mutational landscape of a large cohort of Mexican DLBCL patients. These results need to be confirmed in further studies.

Genomic landscape of early-stage prostate adenocarcinoma in Mexican patients: an exploratory study

BACKGROUND

Health disparities have been highlighted among patient with prostate adenocarcinoma (PRAD) due to ethnicity. Mexican men present a more aggressive disease than other patients resulting in less favorable treatment outcome. We aimed to identify the mutational landscape which could help to reduce the health disparities among minority groups and generate the first genomics exploratory study of PRAD in Mexican patients.

METHODS

Paraffin-embedded formalin-fixed tumoral tissue from 20 Mexican patients with early-stage PRAD treated at The Instituto Nacional de Cancerología, Mexico City from 2017 to 2019 were analyzed. Tumoral DNA was prepared for whole exome sequencing, the resulting files were mapped against h19 using BWA-MEM. Strelka2 and Lancet packages were used to identify single nucleotide variants (SNV) and insertions or deletions. FACETS was used to determine somatic copy number alterations (SCNA). Cancer Genome Interpreter web interface was used to determine the clinical relevance of variants.

RESULTS

Patients were in an early clinical stage and had a mean age of 59.55 years (standard deviation [SD]: 7.1 years) with 90% of them having a Gleason Score of 7. Follow-up time was 48.50 months (SD: 32.77) with recurrences and progression in 30% and 15% of the patients, respectively. NUP98 (20%), CSMD3 (15%) and FAT1 (15%) were the genes most frequently affected by SNV; ARAF (75%) and ZNF419 (70%) were the most frequently affected by losses and gains SNCA's. One quarter of the patients had mutations useful as biomarkers for the use of PARP inhibitors, they comprise mutations in BRCA, RAD54L and ATM. SBS05, DBS03 and ID08 were the most common mutational signatures present in this cohort. No associations with recurrence or progression were identified.

CONCLUSIONS

This pilot study reveals the mutational landscape of early-stage prostate adenocarcinoma in Mexican men, providing a first approach to understand the mutational patterns and actionable mutations in early prostate cancer can inform personalized treatment approaches and reduce the underrepresentation in genomic cancer studies.

Unraveling druggable cancer-driving proteins and targeted drugs using artificial intelligence and multi-omics analyses

The druggable proteome refers to proteins that can bind to small molecules with appropriate chemical affinity, inducing a favorable clinical response. Predicting druggable proteins through screening and in silico modeling is imperative for drug design. To contribute to this field, we developed an accurate predictive classifier for druggable cancer-driving proteins using amino acid composition descriptors of protein sequences and 13 machine learning linear and non-linear classifiers. The optimal classifier was achieved with the support vector machine method, utilizing 200 tri-amino acid composition descriptors. The high performance of the model is evident from an area under the receiver operating characteristics (AUROC) of 0.975 ± 0.003 and an accuracy of 0.929 ± 0.006 (threefold cross-validation). The machine learning prediction model was enhanced with multi-omics approaches, including the target-disease evidence score, the shortest pathways to cancer hallmarks, structure-based ligandability assessment, unfavorable prognostic protein analysis, and the oncogenic variome. Additionally, we performed a drug repurposing analysis to identify drugs with the highest affinity capable of targeting the best predicted proteins. As a result, we identified 79 key druggable cancer-driving proteins with the highest ligandability, and 23 of them demonstrated unfavorable prognostic significance across 16 TCGA PanCancer types: CDKN2A, BCL10, ACVR1, CASP8, JAG1, TSC1, NBN, PREX2, PPP2R1A, DNM2, VAV1, ASXL1, TPR, HRAS, BUB1B, ATG7, MARK3, SETD2, CCNE1, MUTYH, CDKN2C, RB1, and SMARCA4. Moreover, we prioritized 11 clinically relevant drugs targeting these proteins. This strategy effectively predicts and prioritizes biomarkers, therapeutic targets, and drugs for in-depth studies in clinical trials. Scripts are available at https://github.com/muntisa/machine-learning-for-druggable-proteins .

Evolutionary dependency of cancer mutations in gene pairs inferred by nonsynonymous-synonymous mutation ratios

BACKGROUND

Determining the impact of somatic mutations requires understanding the functional relationship of genes acquiring mutations; however, it is largely unknown how mutations in functionally related genes influence each other.

METHODS

We employed non-synonymous-to-synonymous or dNdS ratios to evaluate the evolutionary dependency (ED) of gene pairs, assuming a mutation in one gene of a gene pair can affect the evolutionary fitness of mutations in its partner genes as mutation context. We employed PanCancer- and tumor type-specific mutational profiles to infer the ED of gene pairs and evaluated their biological relevance with respect to gene dependency and drug sensitivity.

RESULTS

We propose that dNdS ratios of gene pairs and their derived cdNS (context-dependent dNdS) scores as measure of ED distinguishing gene pairs either as synergistic (SYN) or antagonistic (ANT). Mutation contexts can induce substantial changes in the evolutionary fitness of mutations in the paired genes, e.g., IDH1 and IDH2 mutation contexts lead to substantial increase and decrease of dNdS ratios of ATRX indels and IDH1 missense mutations corresponding to SYN and ANT relationship with positive and negative cdNS scores, respectively. The impact of gene silencing or knock-outs on cell viability (genetic dependencies) often depends on ED, suggesting that ED can guide the selection of candidates for synthetic lethality such as TCF7L2-KRAS mutations. Using cell line-based drug sensitivity data, the effects of targeted agents on cell lines are often associated with mutations of genes exhibiting ED with the target genes, informing drug sensitizing or resistant mutations for targeted inhibitors, e.g., PRSS1 and CTCF mutations as resistant mutations to EGFR and BRAF inhibitors for lung adenocarcinomas and melanomas, respectively.

CONCLUSIONS

We propose that the ED of gene pairs evaluated by dNdS ratios can advance our understanding of the functional relationship of genes with potential biological and clinical implications.

Tumor Promoters and Opportunities for Molecular Cancer Prevention

Environmental carcinogens increase cancer incidence via both mutagenic and non-mutagenic mechanisms. There are over 500 known or suspected carcinogens classified by the International Agency for Research on Cancer. Sequencing of both cancerous and histologically non-cancerous tissue has been instrumental in improving our understanding of how environmental carcinogens cause cancer. Understanding how and defining which environmental or lifestyle exposures drive cancer will support cancer prevention. Recent research is revisiting the mechanisms of early tumorigenesis, paving the way for an era of molecular cancer prevention. Significance: Recent data have improved our understanding of how carcinogens cause cancer, which may reveal novel opportunities for molecular cancer prevention.

Saturation genome editing maps the functional spectrum of pathogenic VHL alleles

To maximize the impact of precision medicine approaches, it is critical to identify genetic variants underlying disease and to accurately quantify their functional effects. A gene exemplifying the challenge of variant interpretation is the von Hippel-Lindautumor suppressor (VHL). VHL encodes an E3 ubiquitin ligase that regulates the cellular response to hypoxia. Germline pathogenic variants in VHL predispose patients to tumors including clear cell renal cell carcinoma (ccRCC) and pheochromocytoma, and somatic VHL mutations are frequently observed in sporadic renal cancer. Here we optimize and apply saturation genome editing to assay nearly all possible single-nucleotide variants (SNVs) across VHL's coding sequence. To delineate mechanisms, we quantify mRNA dosage effects and compare functional effects in isogenic cell lines. Function scores for 2,268 VHL SNVs identify a core set of pathogenic alleles driving ccRCC with perfect accuracy, inform differential risk across tumor types and reveal new mechanisms by which variants impact function. These results have immediate utility for classifying VHL variants encountered clinically and illustrate how precise functional measurements can resolve pleiotropic and dosage-dependent genotype-phenotype relationships across complete genes.

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.

The mutagenic forces shaping the genomic landscape of lung cancer in never smokers

Lung cancer in never smokers (LCINS) accounts for up to 25% of all lung cancers and has been associated with exposure to secondhand tobacco smoke and air pollution in observational studies. Here, we evaluate the mutagenic exposures in LCINS by examining deep whole-genome sequencing data from a large international cohort of 871 treatment-naïve LCINS recruited from 28 geographical locations within the Sherlock-Lung study. KRAS mutations were 3.8-fold more common in adenocarcinomas of never smokers from North America and Europe, while a 1.6-fold higher prevalence of EGFR and TP53 mutations was observed in adenocarcinomas from East Asia. Signature SBS40a, with unknown cause, was found in most samples and accounted for the largest proportion of single base substitutions in adenocarcinomas, being enriched in EGFR-mutated cases. Conversely, the aristolochic acid signature SBS22a was almost exclusively observed in patients from Taipei. Even though LCINS exposed to secondhand smoke had an 8.3% higher mutational burden and 5.4% shorter telomeres, passive smoking was not associated with driver mutations in cancer driver genes or the activities of individual mutational signatures. In contrast, patients from regions with high levels of air pollution were more likely to have TP53 mutations while exhibiting shorter telomeres and an increase in most types of somatic mutations, including a 3.9-fold elevation of signature SBS4 (q-value=3.1 × 10-5), previously linked mainly to tobacco smoking, and a 76% increase of clock-like signature SBS5 (q-value=5.0 × 10-5). A positive dose-response effect was observed with air pollution levels, which correlated with both a decrease in telomere length and an elevation in somatic mutations, notably attributed to signatures SBS4 and SBS5. Our results elucidate the diversity of mutational processes shaping the genomic landscape of lung cancer in never smokers.

DNA lesion bypass and the stochastic dynamics of transcription-coupled repair

DNA base damage is a major source of oncogenic mutations and disruption to gene expression. The stalling of RNA polymerase II (RNAP) at sites of DNA damage and the subsequent triggering of repair processes have major roles in shaping the genome-wide distribution of mutations, clearing barriers to transcription, and minimizing the production of miscoded gene products. Despite its importance for genetic integrity, key mechanistic features of this transcription-coupled repair (TCR) process are controversial or unknown. Here, we exploited a well-powered in vivo mammalian model system to explore the mechanistic properties and parameters of TCR for alkylation damage at fine spatial resolution and with discrimination of the damaged DNA strand. For rigorous interpretation, a generalizable mathematical model of DNA damage and TCR was developed. Fitting experimental data to the model and simulation revealed that RNA polymerases frequently bypass lesions without triggering repair, indicating that small alkylation adducts are unlikely to be an efficient barrier to gene expression. Following a burst of damage, the efficiency of transcription-coupled repair gradually decays through gene bodies with implications for the occurrence and accurate inference of driver mutations in cancer. The reinitation of transcription from the repair site is not a general feature of transcription-coupled repair, and the observed data is consistent with reinitiation never taking place. Collectively, these results reveal how the directional but stochastic activity of TCR shapes the distribution of mutations following DNA damage.

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.

How CRISPR Is Revolutionizing the Generation of New Models for Cancer Research

Cancers arise through acquisition of mutations in genes that regulate core biological processes like cell proliferation and cell death. Decades of cancer research have led to the identification of genes and mutations causally involved in disease development and evolution, yet defining their precise function across different cancer types and how they influence therapy responses has been challenging. Mouse models have helped define the in vivo function of cancer-associated alterations, and genome-editing approaches using CRISPR have dramatically accelerated the pace at which these models are developed and studied. Here, we highlight how CRISPR technologies have impacted the development and use of mouse models for cancer research and discuss the many ways in which these rapidly evolving platforms will continue to transform our understanding of this disease.

A novel bioinformatic approach reveals cooperation between Cancer/Testis genes in basal-like breast tumors

Breast cancer is the most prevalent type of cancer in women worldwide. Within breast tumors, the basal-like subtype has the worst prognosis, prompting the need for new tools to understand, detect, and treat these tumors. Certain germline-restricted genes show aberrant expression in tumors and are known as Cancer/Testis genes; their misexpression has diagnostic and therapeutic applications. Here we designed a new bioinformatic approach to examine Cancer/Testis gene misexpression in breast tumors. We identify several new markers in Luminal and HER-2 positive tumors, some of which predict response to chemotherapy. We then use machine learning to identify the two Cancer/Testis genes most associated with basal-like breast tumors: HORMAD1 and CT83. We show that these genes are expressed by tumor cells and not by the microenvironment, and that they are not expressed by normal breast progenitors; in other words, their activation occurs de novo. We find these genes are epigenetically repressed by DNA methylation, and that their activation upon DNA demethylation is irreversible, providing a memory of past epigenetic disturbances. Simultaneous expression of both genes in breast cells in vitro has a synergistic effect that increases stemness and activates a transcriptional profile also observed in double-positive tumors. Therefore, we reveal a functional cooperation between Cancer/Testis genes in basal breast tumors; these findings have consequences for the understanding, diagnosis, and therapy of the breast tumors with the worst outcomes.

Discovery of Novel Potential Prognostic Markers and Targeted Therapy to Overcome Chemotherapy Resistance in an Advanced-Stage Wilms Tumor

Wilms tumor (WT), the most prevalent type of renal cancer in children, exhibits overall survival rates exceeding 90%. However, chemotherapy resistance, which occurs in approximately 10% of WT cases, is a major challenge for the treatment of WT, particularly for advanced-stage patients. In this study, we aimed to discover potential mutation markers and drug targets associated with chemotherapy resistance in advanced-stage WT. We performed exome sequencing to detect somatic mutations and molecular targets in 43 WT samples, comprising 26 advanced-stage WTs, of which 7 cases were chemotherapy-resistant. Our analysis revealed four genes (ALPK2, C16orf96, PRKDC, and SVIL) that correlated with chemotherapy resistance and reduced disease-free survival in advanced-stage WT. Additionally, we identified driver mutations in 55 genes within the chemotherapy-resistant group, including 14 druggable cancer driver genes. Based on the mutation profiles of the resistant WT samples, we propose potential therapeutic strategies involving platinum-based agents, PARP inhibitors, and antibiotic/antineoplastic agents. Our findings provide insights into the genetic landscape of WT and offer potential avenues for targeted treatment, particularly for patients with chemotherapy resistance.

APOBEC shapes tumor evolution and age at onset of lung cancer in smokers

APOBEC enzymes are part of the innate immunity and are responsible for restricting viruses and retroelements by deaminating cytosine residues1,2. Most solid tumors harbor different levels of somatic mutations attributed to the off-target activities of APOBEC3A (A3A) and/or APOBEC3B (A3B)3-6. However, how APOBEC3A/B enzymes shape the tumor evolution in the presence of exogenous mutagenic processes is largely unknown. Here, by combining deep whole-genome sequencing with multi-omics profiling of 309 lung cancers from smokers with detailed tobacco smoking information, we identify two subtypes defined by low (LAS) and high (HAS) APOBEC mutagenesis. LAS are enriched for A3B-like mutagenesis and KRAS mutations, whereas HAS for A3A-like mutagenesis and TP53 mutations. Unlike APOBEC3A, APOBEC3B expression is strongly associated with an upregulation of the base excision repair pathway. Hypermutation by unrepaired A3A and tobacco smoking mutagenesis combined with TP53-induced genomic instability can trigger senescence7, apoptosis8, and cell regeneration9, as indicated by high expression of pulmonary healing signaling pathway, stemness markers and distal cell-of-origin in HAS. The expected association of tobacco smoking variables (e.g., time to first cigarette) with genomic/epigenomic changes are not observed in HAS, a plausible consequence of frequent cell senescence or apoptosis. HAS have more neoantigens, slower clonal expansion, and older age at onset compared to LAS, particularly in heavy smokers, consistent with high proportions of newly generated, unmutated cells and frequent immuno-editing. These findings show how heterogeneity in mutational burden across co-occurring mutational processes and cell types contributes to tumor development, with important clinical implications.

New pan-ALK inhibitor-resistant EML4::ALK mutations detected by liquid biopsy in lung cancer patients

ALK and ROS1 fusions are effectively targeted by tyrosine kinase inhibitors (TKIs), however patients inevitably relapse after an initial response, often due to kinase domain mutations. We investigated circulating DNA from TKI-relapsed NSCLC patients by deep-sequencing. New EML4::ALK substitutions, L1198R, C1237Y and L1196P, were identified in the plasma of NSCLC ALK patients and characterized in a Ba/F3 cell model. Variants C1237Y and L1196P demonstrated pan-inhibitor resistance across 5 clinical and 2 investigational TKIs.

A comprehensive clinically informed map of dependencies in cancer cells and framework for target prioritization

Genetic screens in cancer cell lines inform gene function and drug discovery. More comprehensive screen datasets with multi-omics data are needed to enhance opportunities to functionally map genetic vulnerabilities. Here, we construct a second-generation map of cancer dependencies by annotating 930 cancer cell lines with multi-omic data and analyze relationships between molecular markers and cancer dependencies derived from CRISPR-Cas9 screens. We identify dependency-associated gene expression markers beyond driver genes, and observe many gene addiction relationships driven by gain of function rather than synthetic lethal effects. By combining clinically informed dependency-marker associations with protein-protein interaction networks, we identify 370 anti-cancer priority targets for 27 cancer types, many of which have network-based evidence of a functional link with a marker in a cancer type. Mapping these targets to sequenced tumor cohorts identifies tractable targets in different cancer types. This target prioritization map enhances understanding of gene dependencies and identifies candidate anti-cancer targets for drug development.

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.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

A Learning Program for Treatment Recommendations by Molecular Tumor Boards and Artificial Intelligence

Importance

Substantial heterogeneity exists in treatment recommendations across molecular tumor boards (MTBs), especially for biomarkers with low evidence levels; therefore, the learning program is essential.

Objective

To determine whether a learning program sharing treatment recommendations for biomarkers with low evidence levels contributes to the standardization of MTBs and to investigate the efficacy of an artificial intelligence (AI)-based annotation system.

Design, Setting, and Participants

This prospective quality improvement study used 50 simulated cases to assess concordance of treatment recommendations between a central committee and participants. Forty-seven participants applied from April 7 to May 13, 2021. Fifty simulated cases were randomly divided into prelearning and postlearning evaluation groups to assess similar concordance based on previous investigations. Participants included MTBs at hub hospitals, treating physicians at core hospitals, and AI systems. Each participant made treatment recommendations for each prelearning case from registration to June 30, 2021; participated in the learning program on July 18, 2021; and made treatment recommendations for each postlearning case from August 3 to September 30, 2021. Data were analyzed from September 2 to December 10, 2021.

Exposures

The learning program shared the methodology of making appropriate treatment recommendations, especially for biomarkers with low evidence levels.

Main Outcomes and Measures

The primary end point was the proportion of MTBs that met prespecified accreditation criteria for postlearning evaluations (approximately 90% concordance with high evidence levels and approximately 40% with low evidence levels). Key secondary end points were chronological enhancements in the concordance of treatment recommendations on postlearning evaluations from prelearning evaluations. Concordance of treatment recommendations by an AI system was an exploratory end point.

Results

Of the 47 participants who applied, 42 were eligible. The accreditation rate of the MTBs was 55.6% (95% CI, 35.3%-74.5%; P < .001). Concordance in MTBs increased from 58.7% (95% CI, 52.8%-64.4%) to 67.9% (95% CI, 61.0%-74.1%) (odds ratio, 1.40 [95% CI, 1.06-1.86]; P = .02). In postlearning evaluations, the concordance of treatment recommendations by the AI system was significantly higher than that of MTBs (88.0% [95% CI, 68.7%-96.1%]; P = .03).

Conclusions and Relevance

The findings of this quality improvement study suggest that use of a learning program improved the concordance of treatment recommendations provided by MTBs to central ones. Treatment recommendations made by an AI system showed higher concordance than that for MTBs, indicating the potential clinical utility of the AI system.

Hotspot propensity across mutational processes

The sparsity of mutations observed across tumours hinders our ability to study mutation rate variability at nucleotide resolution. To circumvent this, here we investigated the propensity of mutational processes to form mutational hotspots as a readout of their mutation rate variability at single base resolution. Mutational signatures 1 and 17 have the highest hotspot propensity (5-78 times higher than other processes). After accounting for trinucleotide mutational probabilities, sequence composition and mutational heterogeneity at 10 Kbp, most (94-95%) signature 17 hotspots remain unexplained, suggesting a significant role of local genomic features. For signature 1, the inclusion of genome-wide distribution of methylated CpG sites into models can explain most (80-100%) of the hotspot propensity. There is an increased hotspot propensity of signature 1 in normal tissues and de novo germline mutations. We demonstrate that hotspot propensity is a useful readout to assess the accuracy of mutation rate models at nucleotide resolution. This new approach and the findings derived from it open up new avenues for a range of somatic and germline studies investigating and modelling mutagenesis.

A novel patient-derived meningioma spheroid model as a tool to study and treat epithelial-to-mesenchymal transition (EMT) in meningiomas

Meningiomas are the most common intracranial brain tumours. These tumours are heterogeneous and encompass a wide spectrum of clinical aggressivity. Treatment options are limited to surgery and radiotherapy and have a risk of post-operative morbidities and radiation neurotoxicity, reflecting the need for new therapies. Three-dimensional (3D) patient-derived cell culture models have been shown to closely recapitulate in vivo tumour biology, including microenvironmental interactions and have emerged as a robust tool for drug development. Here, we established a novel easy-to-use 3D patient-derived meningioma spheroid model using a scaffold-free approach. Patient-derived meningioma spheroids were characterised and compared to patient tissues and traditional monolayer cultures by histology, genomics, and transcriptomics studies. Patient-derived meningioma spheroids closely recapitulated morphological and molecular features of matched patient tissues, including patient histology, genomic alterations, and components of the immune microenvironment, such as a CD68 + and CD163 + positive macrophage cell population. Comprehensive transcriptomic profiling revealed an increase in epithelial-to-mesenchymal transition (EMT) in meningioma spheroids compared to traditional monolayer cultures, confirming this model as a tool to elucidate EMT in meningioma. Therefore, as proof of concept study, we developed a treatment strategy to target EMT in meningioma. We found that combination therapy using the MER tyrosine kinase (MERTK) inhibitor UNC2025 and the histone deacetylase (HDAC) inhibitor Trichostatin A (TSA) effectively decreased meningioma spheroid viability and proliferation. Furthermore, we demonstrated this combination therapy significantly increased the expression of the epithelial marker E-cadherin and had a repressive effect on WHO grade 2-derived spheroid invasion, which is suggestive of a partial reversal of EMT in meningioma spheroids.

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.

Neurodevelopmental disorders and cancer networks share pathways, but differ in mechanisms, signaling strength, and outcome

Epidemiological studies suggest that individuals with neurodevelopmental disorders (NDDs) are more prone to develop certain types of cancer. Notably, however, the case statistics can be impacted by late discovery of cancer in individuals afflicted with NDDs, such as intellectual disorders, autism, and schizophrenia, which may bias the numbers. As to NDD-associated mutations, in most cases, they are germline while cancer mutations are sporadic, emerging during life. However, somatic mosaicism can spur NDDs, and cancer-related mutations can be germline. NDDs and cancer share proteins, pathways, and mutations. Here we ask (i) exactly which features they share, and (ii) how, despite their commonalities, they differ in clinical outcomes. To tackle these questions, we employed a statistical framework followed by network analysis. Our thorough exploration of the mutations, reconstructed disease-specific networks, pathways, and transcriptome levels and profiles of autism spectrum disorder (ASD) and cancers, point to signaling strength as the key factor: strong signaling promotes cell proliferation in cancer, and weaker (moderate) signaling impacts differentiation in ASD. Thus, we suggest that signaling strength, not activating mutations, can decide clinical outcome.

AQcalc: A web server that identifies weak molecular interactions in protein structures

Weak molecular interactions play an important role in protein structure and function. Computational tools that identify weak molecular interactions are, therefore, valuable for the study of proteins. Here, we present AQcalc, a web server (https://aqcalcbiocomputing.com/) that can be used to identify anion-quadrupole (AQ) interactions, which are weak interactions involving aromatic residue (Trp, Tyr, and Phe) ring edges and anions (Asp, Glu, and phosphate ion) both within proteins and at their interfaces (protein-protein, protein-nucleic acids, and protein-lipid bilayer). AQcalc identifies AQ interactions as well as clusters involving AQ, cation-π, and salt bridges, among others. Utilizing AQcalc we analyzed weak interactions in protein models, even in the absence of experimental structures, to understand the contributions of weak interactions to deleterious structural changes, including those associated with oncogenic and germline disease variants. We identified several deleterious variants with disrupted AQ interactions (comparable in frequency to cation-π disruptions). Amyloid fibrils utilize AQ to bury anions at frequencies that far exceed those observed for globular proteins. AQ interactions were detected three and five times more frequently than the hydrogen-bonded AQ (HBAQ) in fibril structures and protein-lipid bilayer interfaces, respectively. By contrast, AQ and HBAQ interactions were detected with similar frequencies in globular proteins. Collectively, these findings suggest AQcalc will be effective in facilitating fine structural analysis. As other web utilities designed to identify protein residue interaction networks do not report AQ interactions, wide use of AQcalc will enrich our understanding of residue interaction networks and facilitate hypothesis testing by identifying and experimentally characterizing these comparably weak but important interactions.

Pan-cancer proteogenomics connects oncogenic drivers to functional states

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Integrated multi-omics analysis reveals the molecular interplay between circadian clocks and cancer pathogenesis

Circadian rhythms (CRs) are fundamental biological processes that significantly impact human well-being. Disruption of these rhythms can trigger insufficient neurocognitive development, insomnia, mental disorders, cardiovascular diseases, metabolic dysfunctions, and cancer. The field of chronobiology has increased our understanding of how rhythm disturbances contribute to cancer pathogenesis, and how circadian timing influences the efficacy of cancer treatments. As the circadian clock steadily gains recognition as an emerging factor in tumorigenesis, a thorough and comprehensive multi-omics analysis of CR genes/proteins has never been performed. To shed light on this, we performed, for the first time, an integrated data analysis encompassing genomic/transcriptomic alterations across 32 cancer types (n = 10,918 tumors) taken from the PanCancer Atlas, unfavorable prognostic protein analysis, protein-protein interactomics, and shortest distance score pathways to cancer hallmark phenotypes. This data mining strategy allowed us to unravel 31 essential CR-related proteins involved in the signaling crossroad between circadian rhythms and cancer. In the context of drugging the clock, we identified pharmacogenomic clinical annotations and drugs currently in late phase clinical trials that could be considered as potential cancer therapeutic strategies. These findings highlight the diverse roles of CR-related genes/proteins in the realm of cancer research and therapy.

De novo identification of complex traits associated with asthma

Introduction

Asthma is a heterogeneous inflammatory disease often associated with other complex phenotypes. Identifying asthma-associated diseases and uncovering the molecular mechanisms mediating their interaction can help detangle the heterogeneity of asthma. Network analysis is a powerful approach for untangling such inter-disease relationships.

Methods

Here, we integrated information on physical contacts between common single nucleotide polymorphisms (SNPs) and gene expression with expression quantitative trait loci (eQTL) data from the lung and whole blood to construct two tissue-specific spatial gene regulatory networks (GRN). We then located the asthma GRN (level 0) within each tissue-specific GRN by identifying the genes that are functionally affected by asthma-associated spatial eQTLs. Curated protein interaction partners were subsequently identified up to four edges or levels away from the asthma GRN. The eQTLs spatially regulating genes on levels 0-4 were queried against the GWAS Catalog to identify the traits enriched (hypergeometric test; FDR ≤ 0.05) in each level.

Results

We identified 80 and 82 traits significantly enriched in the lung and blood GRNs, respectively. All identified traits were previously reported to be comorbid or associated (positively or negatively) with asthma (e.g., depressive symptoms and lung cancer), except 8 traits whose association with asthma is yet to be confirmed (e.g., reticulocyte count). Our analysis additionally pinpoints the variants and genes that link asthma to the identified asthma-associated traits, a subset of which was replicated in a comorbidity analysis using health records of 26,781 asthma patients in New Zealand.

Discussion

Our discovery approach identifies enriched traits in the regulatory space proximal to asthma, in the tissue of interest, without a priori selection of the interacting traits. The predictions it makes expand our understanding of possible shared molecular interactions and therapeutic targets for asthma, where no cure is currently available.

Highlights from the 1st European cancer dependency map symposium and workshop

The systematic identification of tumour vulnerabilities through perturbational experiments on cancer models, including genome editing and drug screens, is playing a crucial role in combating cancer. This collective effort is known as the Cancer Dependency Map (DepMap). The 1st European Cancer Dependency Map Symposium (EuroDepMap), held in Milan last May, featured talks, a roundtable discussion, and a poster session, showcasing the latest discoveries and future challenges related to the DepMap. The symposium aimed to facilitate interactions among participants across Europe, encourage idea exchange with leading experts, and present their work and future projects. Importantly, it sparked discussions on future endeavours, such as screening more complex cancer models and accounting for tumour evolution.

Actionability classification of variants of unknown significance correlates with functional effect

Genomically-informed therapy requires consideration of the functional impact of genomic alterations on protein expression and/or function. However, a substantial number of variants are of unknown significance (VUS). The MD Anderson Precision Oncology Decision Support (PODS) team developed an actionability classification scheme that categorizes VUS as either "Unknown" or "Potentially" actionable based on their location within functional domains and/or proximity to known oncogenic variants. We then compared PODS VUS actionability classification with results from a functional genomics platform consisting of mutant generation and cell viability assays. 106 (24%) of 438 VUS in 20 actionable genes were classified as oncogenic in functional assays. Variants categorized by PODS as Potentially actionable (N = 204) were more likely to be oncogenic than those categorized as Unknown (N = 230) (37% vs 13%, p = 4.08e-09). Our results demonstrate that rule-based actionability classification of VUS can identify patients more likely to have actionable variants for consideration with genomically-matched therapy.

Impact of risk factors on early cancer evolution

Recent identification of oncogenic cells within healthy tissues and the prevalence of indolent cancers found incidentally at autopsies reveal a greater complexity in tumor initiation than previously appreciated. The human body contains roughly 40 trillion cells of 200 different types that are organized within a complex three-dimensional matrix, necessitating exquisite mechanisms to restrain aberrant outgrowth of malignant cells that have the capacity to kill the host. Understanding how this defense is overcome to trigger tumorigenesis and why cancer is so extraordinarily rare at the cellular level is vital to future prevention therapies. In this review, we discuss how early initiated cells are protected from further tumorigenesis and the non-mutagenic pathways by which cancer risk factors promote tumor growth. By nature, the absence of permanent genomic alterations potentially renders these tumor-promoting mechanisms clinically targetable. Finally, we consider existing strategies for early cancer interception with perspectives on the next steps for molecular cancer prevention.

Challenging gold standard hematology diagnostics through the introduction of whole genome sequencing and artificial intelligence

The diagnosis of hematological malignancies is rather complex and requires the application of a plethora of different assays, techniques and methodologies. Some of the methods, like cytomorphology, have been in use for decades, while other methods, such as next-generation sequencing or even whole genome sequencing (WGS), are relatively new. The application of the methods and the evaluation of the results require distinct skills and knowledge and place different demands on the practitioner. However, even with experienced hematologists, diagnostic ambiguity remains a regular occurrence and the comprehensive analysis of high-dimensional WGS data soon exceeds any human's capacity. Hence, in order to reduce inter-observer variability and to improve the timeliness and accuracy of diagnoses, machine learning based approaches have been developed to assist in the decision making process. Moreover, to achieve the goal of precision oncology, comprehensive genomic profiling is increasingly being incorporated into routine standard of care.

Lung adenocarcinoma promotion by air pollutants

A complete understanding of how exposure to environmental substances promotes cancer formation is lacking. More than 70 years ago, tumorigenesis was proposed to occur in a two-step process: an initiating step that induces mutations in healthy cells, followed by a promoter step that triggers cancer development1. Here we propose that environmental particulate matter measuring ≤2.5 μm (PM2.5), known to be associated with lung cancer risk, promotes lung cancer by acting on cells that harbour pre-existing oncogenic mutations in healthy lung tissue. Focusing on EGFR-driven lung cancer, which is more common in never-smokers or light smokers, we found a significant association between PM2.5 levels and the incidence of lung cancer for 32,957 EGFR-driven lung cancer cases in four within-country cohorts. Functional mouse models revealed that air pollutants cause an influx of macrophages into the lung and release of interleukin-1β. This process results in a progenitor-like cell state within EGFR mutant lung alveolar type II epithelial cells that fuels tumorigenesis. Ultradeep mutational profiling of histologically normal lung tissue from 295 individuals across 3 clinical cohorts revealed oncogenic EGFR and KRAS driver mutations in 18% and 53% of healthy tissue samples, respectively. These findings collectively support a tumour-promoting role for  PM2.5 air pollutants  and provide impetus for public health policy initiatives to address air pollution to reduce disease burden.

Computational Methods Summarizing Mutational Patterns in Cancer: Promise and Limitations for Clinical Applications

Since the rise of next-generation sequencing technologies, the catalogue of mutations in cancer has been continuously expanding. To address the complexity of the cancer-genomic landscape and extract meaningful insights, numerous computational approaches have been developed over the last two decades. In this review, we survey the current leading computational methods to derive intricate mutational patterns in the context of clinical relevance. We begin with mutation signatures, explaining first how mutation signatures were developed and then examining the utility of studies using mutation signatures to correlate environmental effects on the cancer genome. Next, we examine current clinical research that employs mutation signatures and discuss the potential use cases and challenges of mutation signatures in clinical decision-making. We then examine computational studies developing tools to investigate complex patterns of mutations beyond the context of mutational signatures. We survey methods to identify cancer-driver genes, from single-driver studies to pathway and network analyses. In addition, we review methods inferring complex combinations of mutations for clinical tasks and using mutations integrated with multi-omics data to better predict cancer phenotypes. We examine the use of these tools for either discovery or prediction, including prediction of tumor origin, treatment outcomes, prognosis, and cancer typing. We further discuss the main limitations preventing widespread clinical integration of computational tools for the diagnosis and treatment of cancer. We end by proposing solutions to address these challenges using recent advances in machine learning.

Artificial intelligence in hematological diagnostics: Game changer or gadget?

The future of clinical diagnosis and treatment of hematologic diseases will inevitably involve the integration of artificial intelligence (AI)-based systems into routine practice to support the hematologists' decision making. Several studies have shown that AI-based models can already be used to automatically differentiate cells, reliably detect malignant cell populations, support chromosome banding analysis, and interpret clinical variants, contributing to early disease detection and prognosis. However, even the best tool can become useless if it is misapplied or the results are misinterpreted. Therefore, in order to comprehensively judge and correctly apply newly developed AI-based systems, the hematologist must have a basic understanding of the general concepts of machine learning. In this review, we provide the hematologist with a comprehensive overview of various machine learning techniques, their current implementations and approaches in different diagnostic subfields (e.g., cytogenetics, molecular genetics), and the limitations and unresolved challenges of the systems.

Transmicron: accurate prediction of insertion probabilities improves detection of cancer driver genes from transposon mutagenesis screens

Transposon screens are powerful in vivo assays used to identify loci driving carcinogenesis. These loci are identified as Common Insertion Sites (CISs), i.e. regions with more transposon insertions than expected by chance. However, the identification of CISs is affected by biases in the insertion behaviour of transposon systems. Here, we introduce Transmicron, a novel method that differs from previous methods by (i) modelling neutral insertion rates based on chromatin accessibility, transcriptional activity and sequence context and (ii) estimating oncogenic selection for each genomic region using Poisson regression to model insertion counts while controlling for neutral insertion rates. To assess the benefits of our approach, we generated a dataset applying two different transposon systems under comparable conditions. Benchmarking for enrichment of known cancer genes showed improved performance of Transmicron against state-of-the-art methods. Modelling neutral insertion rates allowed for better control of false positives and stronger agreement of the results between transposon systems. Moreover, using Poisson regression to consider intra-sample and inter-sample information proved beneficial in small and moderately-sized datasets. Transmicron is open-source and freely available. Overall, this study contributes to the understanding of transposon biology and introduces a novel approach to use this knowledge for discovering cancer driver genes.

Base editing screens map mutations affecting interferon-γ signaling in cancer

Interferon-γ (IFN-γ) signaling mediates host responses to infection, inflammation and anti-tumor immunity. Mutations in the IFN-γ signaling pathway cause immunological disorders, hematological malignancies, and resistance to immune checkpoint blockade (ICB) in cancer; however, the function of most clinically observed variants remains unknown. Here, we systematically investigate the genetic determinants of IFN-γ response in colorectal cancer cells using CRISPR-Cas9 screens and base editing mutagenesis. Deep mutagenesis of JAK1 with cytidine and adenine base editors, combined with pathway-wide screens, reveal loss-of-function and gain-of-function mutations, including causal variants in hematological malignancies and mutations detected in patients refractory to ICB. We functionally validate variants of uncertain significance in primary tumor organoids, where engineering missense mutations in JAK1 enhanced or reduced sensitivity to autologous tumor-reactive T cells. We identify more than 300 predicted missense mutations altering IFN-γ pathway activity, generating a valuable resource for interpreting gene variant function.

An evolution-based machine learning to identify cancer type-specific driver mutations

Identifying cancer type-specific driver mutations is crucial for illuminating distinct pathologic mechanisms across various tumors and providing opportunities of patient-specific treatment. However, although many computational methods were developed to predict driver mutations in a type-specific manner, the methods still have room to improve. Here, we devise a novel feature based on sequence co-evolution analysis to identify cancer type-specific driver mutations and construct a machine learning (ML) model with state-of-the-art performance. Specifically, relying on 28 000 tumor samples across 66 cancer types, our ML framework outperformed current leading methods of detecting cancer driver mutations. Interestingly, the cancer mutations identified by sequence co-evolution feature are frequently observed in interfaces mediating tissue-specific protein-protein interactions that are known to associate with shaping tissue-specific oncogenesis. Moreover, we provide pre-calculated potential oncogenicity on available human proteins with prediction scores of all possible residue alterations through user-friendly website (http://sbi.postech.ac.kr/w/cancerCE). This work will facilitate the identification of cancer type-specific driver mutations in newly sequenced tumor samples.

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.

Interrogating the Human Diplome: Computational Methods, Emerging Applications, and Challenges

Human DNA sequencing protocols have revolutionized human biology, biomedical science, and clinical practice, but still have very important limitations. One limitation is that most protocols do not separate or assemble (i.e., "phase") the nucleotide content of each of the maternally and paternally derived chromosomal homologs making up the 22 autosomal pairs and the chromosomal pair making up the pseudo-autosomal region of the sex chromosomes. This has led to a dearth of studies and a consequent underappreciation of many phenomena of fundamental importance to basic and clinical genomic science. We discuss a few protocols for obtaining phase information as well as their limitations, including those that could be used in tumor phasing settings. We then describe a number of biological and clinical phenomena that require phase information. These include phenomena that require precise knowledge of the nucleotide sequence in a chromosomal segment from germline or somatic cells, such as DNA binding events, and insight into unique cis vs. trans-acting functionally impactful variant combinations-for example, variants implicated in a phenotype governed by compound heterozygosity. In addition, we also comment on the need for reliable and consensus-based diploid-context computational workflows for variant identification as well as the need for laboratory-based functional verification strategies for validating cis vs. trans effects of variant combinations. We also briefly describe available resources, example studies, as well as areas of further research, and ultimately argue that the science behind the study of human diploidy, referred to as "diplomics," which will be enabled by nucleotide-level resolution of phased genomes, is a logical next step in the analysis of human genome biology.

Reducing uncertainty in genetic testing with Saturation Genome Editing

Accurate interpretation of human genetic data is critical for optimizing outcomes in the era of genomic medicine. Powerful methods for testing genetic variants for functional effects are allowing researchers to characterize thousands of variants across disease genes. Here, we review experimental tools enabling highly scalable assays of variants, focusing specifically on Saturation Genome Editing (SGE). We discuss examples of how this technique is being implemented for variant testing at scale and describe how SGE data for BRCA1 have been clinically validated and used to aid variant interpretation. The initial success at predicting variant pathogenicity with SGE has spurred efforts to expand this and related techniques to many more genes.

Genome-wide mapping of somatic mutation rates uncovers drivers of cancer

Identification of cancer driver mutations that confer a proliferative advantage is central to understanding cancer; however, searches have often been limited to protein-coding sequences and specific non-coding elements (for example, promoters) because of the challenge of modeling the highly variable somatic mutation rates observed across tumor genomes. Here we present Dig, a method to search for driver elements and mutations anywhere in the genome. We use deep neural networks to map cancer-specific mutation rates genome-wide at kilobase-scale resolution. These estimates are then refined to search for evidence of driver mutations under positive selection throughout the genome by comparing observed to expected mutation counts. We mapped mutation rates for 37 cancer types and applied these maps to identify putative drivers within intronic cryptic splice regions, 5' untranslated regions and infrequently mutated genes. Our high-resolution mutation rate maps, available for web-based exploration, are a resource to enable driver discovery genome-wide.

Introducing AI to the molecular tumor board: one direction toward the establishment of precision medicine using large-scale cancer clinical and biological information

Since U.S. President Barack Obama announced the Precision Medicine Initiative in his New Year's State of the Union address in 2015, the establishment of a precision medicine system has been emphasized worldwide, particularly in the field of oncology. With the advent of next-generation sequencers specifically, genome analysis technology has made remarkable progress, and there are active efforts to apply genome information to diagnosis and treatment. Generally, in the process of feeding back the results of next-generation sequencing analysis to patients, a molecular tumor board (MTB), consisting of experts in clinical oncology, genetic medicine, etc., is established to discuss the results. On the other hand, an MTB currently involves a large amount of work, with humans searching through vast databases and literature, selecting the best drug candidates, and manually confirming the status of available clinical trials. In addition, as personalized medicine advances, the burden on MTB members is expected to increase in the future. Under these circumstances, introducing cutting-edge artificial intelligence (AI) technology and information and communication technology to MTBs while reducing the burden on MTB members and building a platform that enables more accurate and personalized medical care would be of great benefit to patients. In this review, we introduced the latest status of elemental technologies that have potential for AI utilization in MTB, and discussed issues that may arise in the future as we progress with AI implementation.

Random forest classifier improving phenylketonuria screening performance in two Chinese populations

Phenylketonuria (PKU) is a genetic disorder with amino acid metabolic defect, which does great harms to the development of newborns and children. Early diagnosis and treatment can effectively prevent the disease progression. Here we developed a PKU screening model using random forest classifier (RFC) to improve PKU screening performance with excellent sensitivity, false positive rate (FPR) and positive predictive value (PPV) in all the validation dataset and two testing Chinese populations. RFC represented outstanding advantages comparing several different classification models based on machine learning and the traditional logistic regression model. RFC is promising to be applied to neonatal PKU screening.

Somatic Mutational Landscape in Mexican Patients: CDH1 Mutations and chr20q13.33 Amplifications Are Associated with Diffuse-Type Gastric Adenocarcinoma

The Hispanic population, compared with other ethnic groups, presents a more aggressive gastric cancer phenotype with higher frequency of diffuse-type gastric adenocarcinoma (GA); this could be related to the mutational landscape of GA in these patients. Using whole-exome sequencing, we sought to present the mutational landscape of GA from 50 Mexican patients who were treated at The Instituto Nacional de Cancerología from 2019 to 2020. We performed a comprehensive statistical analysis to explore the relationship of the genomic variants and clinical data such as tumor histology and presence of signet-ring cell, H. pylori, and EBV. We describe a potentially different mutational landscape between diffuse and intestinal GA in Mexican patients. Patients with intestinal-type GA tended to present a higher frequency of NOTCH1 mutations, copy number gains in cytobands 13.14, 10q23.33, and 12q25.1, and copy number losses in cytobands 7p12, 14q24.2, and 11q13.1; whereas patients with diffuse-type GA tended to present a high frequency of CDH1 mutations and CNV gains in cytobands 20q13.33 and 22q11.21. This is the first description of a mutational landscape of GA in Mexican patients to better understand tumorigenesis in Hispanic patients and lay the groundwork for discovering potential biomarkers and therapeutic targets.

Cancer-driving mutations are enriched in genic regions intolerant to germline variation

Large reference datasets of protein-coding variation in human populations have allowed us to determine which genes and genic subregions are intolerant to germline genetic variation. There is also a growing number of genes implicated in severe Mendelian diseases that overlap with genes implicated in cancer. We hypothesized that cancer-driving mutations might be enriched in genic subregions that are depleted of germline variation relative to somatic variation. We introduce a new metric, OncMTR (oncology missense tolerance ratio), which uses 125,748 exomes in the Genome Aggregation Database (gnomAD) to identify these genic subregions. We demonstrate that OncMTR can significantly predict driver mutations implicated in hematologic malignancies. Divergent OncMTR regions were enriched for cancer-relevant protein domains, and overlaying OncMTR scores on protein structures identified functionally important protein residues. Last, we performed a rare variant, gene-based collapsing analysis on an independent set of 394,694 exomes from the UK Biobank and find that OncMTR markedly improves genetic signals for hematologic malignancies.