Passenger Mutations in More Than 2,500 Cancer Genomes: Overall Molecular Functional Impact and Consequences

Defining clinically useful biomarkers of immune checkpoint inhibitors in solid tumours

Although more than a decade has passed since the approval of immune checkpoint inhibitors (ICIs) for the treatment of melanoma and non-small-cell lung, breast and gastrointestinal cancers, many patients still show limited response. US Food and Drug Administration (FDA)-approved biomarkers include programmed cell death 1 ligand 1 (PDL1) expression, microsatellite status (that is, microsatellite instability-high (MSI-H)) and tumour mutational burden (TMB), but these have limited utility and/or lack standardized testing approaches for pan-cancer applications. Tissue-based analytes (such as tumour gene signatures, tumour antigen presentation or tumour microenvironment profiles) show a correlation with immune response, but equally, these demonstrate limited efficacy, as they represent a single time point and a single spatial assessment. Patient heterogeneity as well as inter- and intra-tumoural differences across different tissue sites and time points represent substantial challenges for static biomarkers. However, dynamic biomarkers such as longitudinal biopsies or novel, less-invasive markers such as blood-based biomarkers, radiomics and the gut microbiome show increasing potential for the dynamic identification of ICI response, and patient-tailored predictors identified through neoadjuvant trials or novel ex vivo tumour models can help to personalize treatment. In this Perspective, we critically assess the multiple new static, dynamic and patient-specific biomarkers, highlight the newest consortia and trial efforts, and provide recommendations for future clinical trials to make meaningful steps forwards in the field.

SAP30, an oncogenic driver of progression, poor survival, and drug resistance in neuroblastoma

Neuroblastoma is the most devastating extracranial solid malignancy in children. Despite an intense treatment regimen, the prognosis for high-risk neuroblastoma patients remains poor, with less than 40% survival. So far, MYCN amplification status is considered the most prognostic factor but corresponds to only ∼25% of neuroblastoma patients. Therefore, it is essential to identify a better prognosis and therapy response marker in neuroblastoma patients. We applied robust bioinformatic data mining tools, such as weighted gene co-expression network analysis, cisTarget, and single-cell regulatory network inference and clustering on two neuroblastoma patient datasets. We found Sin3A-associated protein 30 (SAP30), a driver transcription factor positively associated with high-risk, progression, stage 4, and poor survival in neuroblastoma patient cohorts. Tumors of high-risk neuroblastoma patients and relapse-specific patient-derived xenografts showed higher SAP30 levels. The advanced pharmacogenomic analysis and CRISPR-Cas9 screens indicated that SAP30 essentiality is associated with cisplatin resistance and further showed higher levels in cisplatin-resistant patient-derived xenograft tumor cell lines. Silencing of SAP30 induced cell death in vitro and led to a reduced tumor burden and size in vivo. Altogether, these results indicate that SAP30 is a better prognostic and cisplatin-resistance marker and thus a potential drug target in high-risk neuroblastoma.

The mitochondrial calcium uniporter: Balancing tumourigenic and anti-tumourigenic responses

Increased malignancy and poor treatability associated with solid tumour cancers have commonly been attributed to mitochondrial calcium (Ca2+) dysregulation. The mitochondrial Ca2+ uniporter complex (mtCU) is the predominant mode of Ca2+ uptake into the mitochondrial matrix. The main components of mtCU are the pore-forming mitochondrial Ca2+ uniporter (MCU) subunit, MCU dominant-negative beta (MCUb) subunit, essential MCU regulator (EMRE) and the gatekeeping mitochondrial Ca2+ uptake 1 and 2 (MICU1 and MICU2) proteins. In this review, we describe mtCU-mediated mitochondrial Ca2+ dysregulation in solid tumour cancer types, finding enhanced mtCU activity observed in colorectal cancer, breast cancer, oral squamous cell carcinoma, pancreatic cancer, hepatocellular carcinoma and embryonal rhabdomyosarcoma. By contrast, decreased mtCU activity is associated with melanoma, whereas the nature of mtCU dysregulation remains unclear in glioblastoma. Furthermore, we show that numerous polymorphisms associated with cancer may alter phosphorylation sites on the pore forming MCU and MCUb subunits, which cluster at interfaces with EMRE. We highlight downstream/upstream biomolecular modulators of MCU and MCUb that alter mtCU-mediated mitochondrial Ca2+ uptake and may be used as biomarkers or to aid in the development of novel cancer therapeutics. Additionally, we provide an overview of the current small molecule inhibitors of mtCU that interact with the Asp residue of the critical Asp-Ile-Met-Glu motif or through other allosteric regulatory mechanisms to block Ca2+ permeation. Finally, we describe the relationship between MCU- and MCUb-mediating microRNAs and mitochondrial Ca2+ uptake that should be considered in the discovery of new treatment approaches for cancer.

The war between the immune system and the tumor - using immune biomarkers as tracers

Nowadays, immunotherapy is one of the most promising anti-tumor therapeutic strategy. Specifically, immune-related targets can be used to predict the efficacy and side effects of immunotherapy and monitor the tumor immune response. In the past few decades, increasing numbers of novel immune biomarkers have been found to participate in certain links of the tumor immunity to contribute to the formation of immunosuppression and have entered clinical trials. Here, we systematically reviewed the oncogenesis and progression of cancer in the view of anti-tumor immunity, particularly in terms of tumor antigen expression (related to tumor immunogenicity) and tumor innate immunity to complement the cancer-immune cycle. From the perspective of integrated management of chronic cancer, we also appraised emerging factors affecting tumor immunity (including metabolic, microbial, and exercise-related markers). We finally summarized the clinical studies and applications based on immune biomarkers. Overall, immune biomarkers participate in promoting the development of more precise and individualized immunotherapy by predicting, monitoring, and regulating tumor immune response. Therefore, targeting immune biomarkers may lead to the development of innovative clinical applications.

Whole-genome mapping of APOBEC mutagenesis in metastatic urothelial carcinoma identifies driver hotspot mutations and a novel mutational signature

Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) enzymes mutate specific DNA sequences and hairpin-loop structures, challenging the distinction between passenger and driver hotspot mutations. Here, we characterized 115 whole genomes of metastatic urothelial carcinoma (mUC) to identify APOBEC mutagenic hotspot drivers. APOBEC-associated mutations were detected in 92% of mUCs and were equally distributed across the genome, while APOBEC hotspot mutations (ApoHMs) were enriched in open chromatin. Hairpin loops were frequent targets of didymi (twins in Greek), two hotspot mutations characterized by the APOBEC SBS2 signature, in conjunction with an uncharacterized mutational context (Ap[C>T]). Next, we developed a statistical framework that identified ApoHMs as drivers in coding and non-coding genomic regions of mUCs. Our results and statistical framework were validated in independent cohorts of 23 non-metastatic UCs and 3,744 samples of 17 metastatic cancers, identifying cancer-type-specific drivers. Our study highlights the role of APOBEC in cancer development and may contribute to developing novel targeted therapy options for APOBEC-driven cancers.

Germline structural variation globally impacts the cancer transcriptome including disease-relevant genes

Germline variation and somatic alterations contribute to the molecular profile of cancers. We combine RNA with whole genome sequencing across 1,218 cancer patients to determine the extent germline structural variants (SVs) impact expression of nearby genes. For hundreds of genes, recurrent and common germline SV breakpoints within 100 kb associate with increased or decreased expression in tumors spanning various tissues of origin. A significant fraction of germline SV expression associations involves duplication of intergenic enhancers or 3' UTR disruption. Genes altered by both somatic and germline SVs include ATRX and CEBPA. Genes essential in cancer cell lines include BARD1 and IRS2. Genes with both expression and germline SV breakpoint patterns associated with patient survival include GCLM. Our results capture a class of phenotypic variation at work in the disease setting, including genes with cancer roles. Specific germline SVs represent potential cancer risk variants for genetic testing, including those involving genes with targeting implications.

Latent evolutionary signatures: a general framework for analysing music and cultural evolution

Cultural processes of change bear many resemblances to biological evolution. The underlying units of non-biological evolution have, however, remained elusive, especially in the domain of music. Here, we introduce a general framework to jointly identify underlying units and their associated evolutionary processes. We model musical styles and principles of organization in dimensions such as harmony and form as following an evolutionary process. Furthermore, we propose that such processes can be identified by extracting latent evolutionary signatures from musical corpora, analogously to identifying mutational signatures in genomics. These signatures provide a latent embedding for each song or musical piece. We develop a deep generative architecture for our model, which can be viewed as a type of variational autoencoder with an evolutionary prior constraining the latent space; specifically, the embeddings for each song are tied together via an energy-based prior, which encourages songs close in evolutionary space to share similar representations. As illustration, we analyse songs from the McGill Billboard dataset. We find frequent chord transitions and formal repetition schemes and identify latent evolutionary signatures related to these features. Finally, we show that the latent evolutionary representations learned by our model outperform non-evolutionary representations in such tasks as period and genre prediction.

Variant Allele Frequency Analysis of Circulating Tumor DNA as a Promising Tool in Assessing the Effectiveness of Treatment in Non-Small Cell Lung Carcinoma Patients

Despite the different possible paths of treatment, lung cancer remains one of the leading causes of death in oncological patients. New tools guiding the therapeutic process are under scientific investigation, and one of the promising indicators of the effectiveness of therapy in patients with NSCLC is variant allele frequency (VAF) analysis. VAF is a metric characterized as the measurement of the specific variant allele proportion within a genomic locus, and it can be determined using methods based on NGS or PCR. It can be assessed using not only tissue samples but also ctDNA (circulating tumor DNA) isolated from liquid biopsy. The non-invasive characteristic of liquid biopsy enables a more frequent collection of material and increases the potential of VAF analysis in monitoring therapy. Several studies have been performed on patients with NSCLC to evaluate the possibility of VAF usage. The research carried out so far demonstrates that the evaluation of VAF dynamics may be useful in monitoring tumor progression, remission, and recurrence during or after treatment. Moreover, the use of VAF analysis appears to be beneficial in making treatment decisions. However, several issues require better understanding and standardization before VAF testing can be implemented in clinical practice. In this review, we discuss the difficulties in the application of ctDNA VAF analysis in clinical routine, discussing the diagnostic and methodological challenges in VAF measurement in liquid biopsy. We highlight the possible applications of VAF-based measurements that are under consideration in clinical trials in the monitoring of personalized treatments for patients with NSCLC.

Aneuploid serves as a prognostic marker and favors immunosuppressive microenvironment in ovarian cancer

Ovarian cancer is the most lethal gynecologic neoplasm, and most patients experience recurrence and chemoresistance. Even the promising immunotherapy showed limited efficacy in ovarian cancer, probably due to the immunosuppressive microenvironment. However, the behind mechanisms of the immune exclusion or cold phenotype in ovarian cancer still remain to be explored. As a cancer dominated by copy number variations instead of mutations, ovarian cancer contains a high fraction of aneuploid, which might correlate with immune inhibition. Nevertheless, whether or how aneuploid affects ovarian cancer is still unclear. For exploring the role of aneuploid cancer cells and the potential ploidy-immune relationship, herein, the ploidy information was first comprehensively analyzed combining the karyotype data and copy number variation data obtained from Mitelman and cBioPortal databases, respectively. Ovarian cancer showed strong ploidy heterogeneity, with high fraction of aneuploid and recurrent arm-level and whole chromosome changes. Furthermore, clinical parameters were compared between the highly-aneuploid and the near-diploid ovarian cancers. Aneuploid indicated high grade, poor overall survival and poor disease-free survival in ovarian cancer. To understand the biofunction affected by aneuploid, the differentially expressed genes between the highly-aneuploid and the near-diploid groups were analyzed. Transcription data suggested that aneuploid cancer correlated with deregulated MHC expression, abnormal antigen presentation, and less infiltration of macrophages and activated T cells and higher level of T cell exclusion. Furthermore, the ploidy-MHC association was verified using the Human Protein Atlas database. All these data supported that aneuploid might be promising for cancer management and immune surveillance in ovarian cancer.

Chargaff's second parity rule lies at the origin of additive genetic interactions in quantitative traits to make omnigenic selection possible

Background

Francis Crick's central dogma provides a residue-by-residue mechanistic explanation of the flow of genetic information in living systems. However, this principle may not be sufficient for explaining how random mutations cause continuous variation of quantitative highly polygenic complex traits. Chargaff's second parity rule (CSPR), also referred to as intrastrand DNA symmetry, defined as near-exact equalities G ≈ C and A ≈ T within a single DNA strand, is a statistical property of cellular genomes. The phenomenon of intrastrand DNA symmetry was discovered more than 50 years ago; at present, it remains unclear what its biological role is, what the mechanisms are that force cellular genomes to comply strictly with CSPR, and why genomes of certain noncellular organisms have broken intrastrand DNA symmetry. The present work is aimed at studying a possible link between intrastrand DNA symmetry and the origin of genetic interactions in quantitative traits.

Methods

Computational analysis of single-nucleotide polymorphisms in human and mouse populations and of nucleotide composition biases at different codon positions in bacterial and human proteomes.

Results

The analysis of mutation spectra inferred from single-nucleotide polymorphisms observed in murine and human populations revealed near-exact equalities of numbers of reverse complementary mutations, indicating that random genetic variations obey CSPR. Furthermore, nucleotide compositions of coding sequences proved to be statistically interwoven via CSPR because pyrimidine bias at the 3rd codon position compensates purine bias at the 1st and 2nd positions.

Conclusions

According to Fisher's infinitesimal model, we propose that accumulation of reverse complementary mutations results in a continuous phenotypic variation due to small additive effects of statistically interwoven genetic variations. Therefore, additive genetic interactions can be inferred as a statistical entanglement of nucleotide compositions of separate genetic loci. CSPR challenges the neutral theory of molecular evolution-because all random mutations participate in variation of a trait-and provides an alternative solution to Haldane's dilemma by making a gene function diffuse. We propose that CSPR is symmetry of Fisher's infinitesimal model and that genetic information can be transferred in an implicit contactless manner.

Incorporating mutational heterogeneity to identify genes that are enriched for synonymous mutations in cancer

BACKGROUND

Synonymous mutations, which change the DNA sequence but not the encoded protein sequence, can affect protein structure and function, mRNA maturation, and mRNA half-lives. The possibility that synonymous mutations might be enriched in cancer has been explored in several recent studies. However, none of these studies control for all three types of mutational heterogeneity (patient, histology, and gene) that are known to affect the accurate identification of non-synonymous cancer-associated genes. Our goal is to adopt the current standard for non-synonymous mutations in an investigation of synonymous mutations.

RESULTS

Here, we create an algorithm, MutSigCVsyn, an adaptation of MutSigCV, to identify cancer-associated genes that are enriched for synonymous mutations based on a non-coding background model that takes into account the mutational heterogeneity across these levels. Using MutSigCVsyn, we first analyzed 2572 cancer whole-genome samples from the Pan-cancer Analysis of Whole Genomes (PCAWG) to identify non-synonymous cancer drivers as a quality control. Indicative of the algorithm accuracy we find that 58.6% of these candidate genes were also found in Cancer Census Gene (CGC) list, and 66.2% were found within the PCAWG cancer driver list. We then applied it to identify 30 putative cancer-associated genes that are enriched for synonymous mutations within the same samples. One of the promising gene candidates is the B cell lymphoma 2 (BCL-2) gene. BCL-2 regulates apoptosis by antagonizing the action of proapoptotic BCL-2 family member proteins. The synonymous mutations in BCL2 are enriched in its anti-apoptotic domain and likely play a role in cancer cell proliferation.

CONCLUSION

Our study introduces MutSigCVsyn, an algorithm that accounts for mutational heterogeneity at patient, histology, and gene levels, to identify cancer-associated genes that are enriched for synonymous mutations using whole genome sequencing data. We identified 30 putative candidate genes that will benefit from future experimental studies on the role of synonymous mutations in cancer biology.

Heterogeneity and individualized treatment of microenvironment in glioblastoma (Review)

The heterogeneity of glioblastoma can suppress immune cell function and lead to immune evasion, which presents a challenge in developing effective molecular therapies for tumor cells. However, the study of tumor immune heterogeneity holds great potential for clinical immunotherapy. Liquid biopsy is a useful tool for accurately monitoring dynamic changes in tumor immune heterogeneity and the tumor microenvironment. This paper explores the heterogeneity of glioblastoma and the immune microenvironment, providing a therapeutic basis for individualized treatment. Using liquid biopsy technology as a new diagnostic method, innovative treatment strategies may be implemented for patients with glioblastoma to improve their outcomes.

Clustering malignant cell states using universally variable genes

Single-cell RNA sequencing (scRNA-seq) has revealed important insights into the heterogeneity of malignant cells. However, sample-specific genomic alterations often confound such analysis, resulting in patient-specific clusters that are difficult to interpret. Here, we present a novel approach to address the issue. By normalizing gene expression variances to identify universally variable genes (UVGs), we were able to reduce the formation of sample-specific clusters and identify underlying molecular hallmarks in malignant cells. In contrast to highly variable genes vulnerable to a specific sample bias, UVGs led to better detection of clusters corresponding to distinct malignant cell states. Our results demonstrate the utility of this approach for analyzing scRNA-seq data and suggest avenues for further exploration of malignant cell heterogeneity.

VIPpred: a novel model for predicting variant impact on phosphorylation events driving carcinogenesis

Disrupted protein phosphorylation due to genetic variation is a widespread phenomenon that triggers oncogenic transformation of healthy cells. However, few relevant phosphorylation disruption events have been verified due to limited biological experimental methods. Because of the lack of reliable benchmark datasets, current bioinformatics methods primarily use sequence-based traits to study variant impact on phosphorylation (VIP). Here, we increased the number of experimentally supported VIP events from less than 30 to 740 by manually curating and reanalyzing multi-omics data from 916 patients provided by the Clinical Proteomic Tumor Analysis Consortium. To predict VIP events in cancer cells, we developed VIPpred, a machine learning method characterized by multidimensional features that exhibits robust performance across different cancer types. Our method provided a pan-cancer landscape of VIP events, which are enriched in cancer-related pathways and cancer driver genes. We found that variant-induced increases in phosphorylation events tend to inhibit the protein degradation of oncogenes and promote tumor suppressor protein degradation. Our work provides new insights into phosphorylation-related cancer biology as well as novel avenues for precision therapy.

Global analysis of suppressor mutations that rescue human genetic defects

BACKGROUND

Genetic suppression occurs when the deleterious effects of a primary "query" mutation, such as a disease-causing mutation, are rescued by a suppressor mutation elsewhere in the genome.

METHODS

To capture existing knowledge on suppression relationships between human genes, we examined 2,400 published papers for potential interactions identified through either genetic modification of cultured human cells or through association studies in patients.

RESULTS

The resulting network encompassed 476 unique suppression interactions covering a wide spectrum of diseases and biological functions. The interactions frequently linked genes that operate in the same biological process. Suppressors were strongly enriched for genes with a role in stress response or signaling, suggesting that deleterious mutations can often be buffered by modulating signaling cascades or immune responses. Suppressor mutations tended to be deleterious when they occurred in absence of the query mutation, in apparent contrast with their protective role in the presence of the query. We formulated and quantified mechanisms of genetic suppression that could explain 71% of interactions and provided mechanistic insight into disease pathology. Finally, we used these observations to predict suppressor genes in the human genome.

CONCLUSIONS

The global suppression network allowed us to define principles of genetic suppression that were conserved across diseases, model systems, and species. The emerging frequency of suppression interactions among human genes and range of underlying mechanisms, together with the prevalence of suppression in model organisms, suggest that compensatory mutations may exist for most genetic diseases.

Precision medicine in the era of multi-omics: can the data tsunami guide rational treatment decision?

Precision medicine for cancer is rapidly moving to an approach that integrates multiple dimensions of the biology in order to model mechanisms of cancer progression in each patient. The discovery of multiple drivers per tumor challenges medical decision that faces several treatment options. Drug sensitivity depends on the actionability of the target, its clonal or subclonal origin and coexisting genomic alterations. Sequencing has revealed a large diversity of drivers emerging at treatment failure, which are potential targets for clinical trials or drug repurposing. To effectively prioritize therapies, it is essential to rank genomic alterations based on their proven actionability. Moving beyond primary drivers, the future of precision medicine necessitates acknowledging the intricate spatial and temporal heterogeneity inherent in cancer. The advent of abundant complex biological data will make artificial intelligence algorithms indispensable for thorough analysis. Here, we will discuss the advancements brought by the use of high-throughput genomics, the advantages and limitations of precision medicine studies and future perspectives in this field.

Global impact of somatic structural variation on the cancer proteome

Both proteome and transcriptome data can help assess the relevance of non-coding somatic mutations in cancer. Here, we combine mass spectrometry-based proteomics data with whole genome sequencing data across 1307 human tumors spanning various tissues to determine the extent somatic structural variant (SV) breakpoint patterns impact protein expression of nearby genes. We find that about 25% of the hundreds of genes with SV-associated cis-regulatory alterations at the mRNA level are similarly associated at the protein level. SVs associated with enhancer hijacking, retrotransposon translocation, altered DNA methylation, or fusion transcripts are implicated in protein over-expression. SVs combined with altered protein levels considerably extend the numbers of patients with tumors somatically altered for critical pathways. We catalog both SV breakpoint patterns involving patient survival and genes with nearby SV breakpoints associated with increased cell dependency in cancer cell lines. Pan-cancer proteogenomics identifies targetable non-coding alterations, by virtue of the associated deregulated genes.

SNPeffect 5.0: large-scale structural phenotyping of protein coding variants extracted from next-generation sequencing data using AlphaFold models

BACKGROUND

Next-generation sequencing technologies yield large numbers of genetic alterations, of which a subset are missense variants that alter an amino acid in the protein product. These variants can have a potentially destabilizing effect leading to an increased risk of misfolding and aggregation. Multiple software tools exist to predict the effect of single-nucleotide variants on proteins, however, a pipeline integrating these tools while starting from an NGS data output list of variants is lacking.

RESULTS

The previous version SNPeffect 4.0 (De Baets in Nucleic Acids Res 40(D1):D935-D939, 2011) provided an online database containing pre-calculated variant effects and low-throughput custom variant analysis. Here, we built an automated and parallelized pipeline that analyzes the impact of missense variants on the aggregation propensity and structural stability of proteins starting from the Variant Call Format as input. The pipeline incorporates the AlphaFold Protein Structure Database to achieve high coverage for structural stability analyses using the FoldX force field. The effect on aggregation-propensity is analyzed using the established predictors TANGO and WALTZ. The pipeline focuses solely on the human proteome and can be used to analyze proteome stability/damage in a given sample based on sequencing results.

CONCLUSION

We provide a bioinformatics pipeline that allows structural phenotyping from sequencing data using established stability and aggregation predictors including FoldX, TANGO, and WALTZ; and structural proteome coverage provided by the AlphaFold database. The pipeline and installation guide are freely available for academic users on https://github.com/vibbits/snpeffect and requires a computer cluster.

LARP4A and LARP4B in cancer: The new kids on the block

Recent developments have mounted a stunning body of evidence underlying the importance of RNA binding proteins (RBPs) in cancer research. In this minireview we focus on LARP4A and LARP4B, two paralogs belonging to the superfamily of La-related proteins, and provide a critical overview of current research, including their roles in cancer pathogenesis and cell proliferation, migration, cell cycle and apoptosis. We highlight current controversies surrounding LARP4A and LARP4B and conclude that their complex roles in tumorigenesis are cell-, tissue- and context-dependent, warning that caution must be exercised before categorising either protein as an oncoprotein or tumour-suppressor. We also reveal that LARP4A and LARP4B have often been confused with one another, adding uncertainty in delineating their functions. We suggest that further functional and mechanistic studies of LARP4 proteins present significant challenges for future investigations to recognise the vital contributions of these RBPs in cancer research.

Tumour mutations in long noncoding RNAs enhance cell fitness

Long noncoding RNAs (lncRNAs) are linked to cancer via pathogenic changes in their expression levels. Yet, it remains unclear whether lncRNAs can also impact tumour cell fitness via function-altering somatic "driver" mutations. To search for such driver-lncRNAs, we here perform a genome-wide analysis of fitness-altering single nucleotide variants (SNVs) across a cohort of 2583 primary and 3527 metastatic tumours. The resulting 54 mutated and positively-selected lncRNAs are significantly enriched for previously-reported cancer genes and a range of clinical and genomic features. A number of these lncRNAs promote tumour cell proliferation when overexpressed in in vitro models. Our results also highlight a dense SNV hotspot in the widely-studied NEAT1 oncogene. To directly evaluate the functional significance of NEAT1 SNVs, we use in cellulo mutagenesis to introduce tumour-like mutations in the gene and observe a significant and reproducible increase in cell fitness, both in vitro and in a mouse model. Mechanistic studies reveal that SNVs remodel the NEAT1 ribonucleoprotein and boost subnuclear paraspeckles. In summary, this work demonstrates the utility of driver analysis for mapping cancer-promoting lncRNAs, and provides experimental evidence that somatic mutations can act through lncRNAs to enhance pathological cancer cell fitness.

Unified views on variant impact across many diseases

Genomic studies of human disorders are often performed by distinct research communities (i.e., focused on rare diseases, common diseases, or cancer). Despite underlying differences in the mechanistic origin of different disease categories, these studies share the goal of identifying causal genomic events that are critical for the clinical manifestation of the disease phenotype. Moreover, these studies face common challenges, including understanding the complex genetic architecture of the disease, deciphering the impact of variants on multiple scales, and interpreting noncoding mutations. Here, we highlight these challenges in depth and argue that properly addressing them will require a more unified vocabulary and approach across disease communities. Toward this goal, we present a unified perspective on relating variant impact to various genomic disorders.

Unraveling the Drivers of Tumorigenesis in the Context of Evolution: Theoretical Models and Bioinformatics Tools

Cancer originates from somatic cells that have accumulated mutations. These mutations alter the phenotype of the cells, allowing them to escape homeostatic regulation that maintains normal cell numbers. The emergence of malignancies is an evolutionary process in which the random accumulation of somatic mutations and sequential selection of dominant clones cause cancer cells to proliferate. The development of technologies such as high-throughput sequencing has provided a powerful means to measure subclonal evolutionary dynamics across space and time. Here, we review the patterns that may be observed in cancer evolution and the methods available for quantifying the evolutionary dynamics of cancer. An improved understanding of the evolutionary trajectories of cancer will enable us to explore the molecular mechanism of tumorigenesis and to design tailored treatment strategies.

The significance of co-mutations in EGFR-mutated non-small cell lung cancer: Optimizing the efficacy of targeted therapies?

Non-small cell lung cancer (NSCLC) is the most common cause of cancer death worldwide. In non-squamous NSCLC, the identification of oncogenic drivers and the development of target-specific molecules led to remarkable progress in therapeutic strategies and overall survival over the last decade. Nevertheless, responses are limited by systematically acquired mechanisms of resistance early on after starting a targeted therapy. Moreover, mounting evidence has demonstrated that each oncogenic-driven cluster is actually heterogeneous in terms of molecular features, clinical behaviour, and sensitivity to targeted therapy. In this review, we aimed to examine the prognostic and predictive significance of oncogene-driven co-mutations, focusing mainly on EGFR and TP53. A narrative review was performed by searching MEDLINE databases for English articles published over the last decade (from January 2012 until November 2022). The bibliographies of key references were manually reviewed to select those eligible for the topic. The genetic landscape of EGFR-mutated NSCLC is more complicated than what is known so far. In particular, the occurrence of TP53 co-mutations stratify patients carrying EGFR mutations in terms of treatment response. The study provides a deeper understanding of the mechanisms underlying the variability of the genetic landscape of EGFR-mutated NSCLC and summarizes notably the clinical importance of TP53 co-mutations for an open avenue to more properly addressing the clinical decision-making in the near future.

Cell Adaptive Fitness and Cancer Evolutionary Dynamics

Genome instability of cancer cells translates into increased entropy and lower information processing capacity, leading to metabolic reprograming toward higher energy states, presumed to be aligned with a cancer growth imperative. Dubbed as the cell adaptive fitness, the proposition postulates that the coupling between cell signaling and metabolism constrains cancer evolutionary dynamics along trajectories privileged by the maintenance of metabolic sufficiency for survival. In particular, the conjecture postulates that clonal expansion becomes restricted when genetic alterations induce a sufficiently high level of disorder, that is, high entropy, in the regulatory signaling network, abrogating as a result the ability of cancer cells to successfully replicate, leading to a stage of clonal stagnation. The proposition is analyzed in the context of an in-silico model of tumor evolutionary dynamics to illustrate how cell-inherent adaptive fitness may predictably constrain clonal evolution of tumors, which would have significant implications for the design of adaptive cancer therapies.

Estimation of Neutral Mutation Rates and Quantification of Somatic Variant Selection Using cancereffectsizeR

Somatic nucleotide mutations can contribute to cancer cell survival, proliferation, and pathogenesis. Although research has focused on identifying which mutations are "drivers" versus "passengers," quantifying the proliferative effects of specific variants within clinically relevant contexts could reveal novel aspects of cancer biology. To enable researchers to estimate these cancer effects, we developed cancereffectsizeR, an R package that organizes somatic variant data, facilitates mutational signature analysis, calculates site-specific mutation rates, and tests models of selection. Built-in models support effect estimation from single nucleotides to genes. Users can also estimate epistatic effects between paired sets of variants, or design and test custom models. The utility of cancer effect was validated by showing in a pan-cancer dataset that somatic variants classified as likely pathogenic or pathogenic in ClinVar exhibit substantially higher effects than most other variants. Indeed, cancer effect was a better predictor of pathogenic status than variant prevalence or functional impact scores. In addition, the application of this approach toward pairwise epistasis in lung adenocarcinoma showed that driver mutations in BRAF, EGFR, or KRAS typically reduce selection for alterations in the other two genes. Companion reference data packages support analyses using the hg19 or hg38 human genome builds, and a reference data builder enables use with any species or custom genome build with available genomic and transcriptomic data. A reference manual, tutorial, and public source code repository are available at https://townsend-lab-yale.github.io/cancereffectsizeR. Comprehensive estimation of cancer effects of somatic mutations can provide insights into oncogenic trajectories, with implications for cancer prognosis and treatment.

SIGNIFICANCE

An R package provides streamlined, customizable estimation of underlying nucleotide mutation rates and of the oncogenic and epistatic effects of mutations in cancer cohorts.

Privacy-preserving cancer type prediction with homomorphic encryption

Cancer genomics tailors diagnosis and treatment based on an individual's genetic information and is the crux of precision medicine. However, analysis and maintenance of high volume of genetic mutation data to build a machine learning (ML) model to predict the cancer type is a computationally expensive task and is often outsourced to powerful cloud servers, raising critical privacy concerns for patients' data. Homomorphic encryption (HE) enables computation on encrypted data, thus, providing cryptographic guarantees to protect privacy. But restrictive overheads of encrypted computation deter its usage. In this work, we explore the challenges of privacy preserving cancer type prediction using a dataset consisting of more than 2 million genetic mutations from 2713 patients for several cancer types by building a highly accurate ML model and then implementing its privacy preserving version in HE. Our solution for cancer type inference encodes somatic mutations based on their impact on the cancer genomes into the feature space and then uses statistical tests for feature selection. We propose a fast matrix multiplication algorithm for HE-based model. Our final model achieves 0.98 micro-average area under curve improving accuracy from 70.08 to 83.61% , being 550 times faster than the standard matrix multiplication-based privacy-preserving models. Our tool can be found at https://github.com/momalab/octal-candet .

Shared Cancer Dataset Analysis Identifies and Predicts the Quantitative Effects of Pan-Cancer Somatic Driver Variants

Driver mutations endow tumors with selective advantages and produce an array of pathogenic effects. Determining the function of somatic variants is important for understanding cancer biology and identifying optimal therapies. Here, we compiled a shared dataset from several cancer genomic databases. Two measures were applied to 535 cancer genes based on observed and expected frequencies of driver variants as derived from cancer-specific rates of somatic mutagenesis. The first measure comprised a binary classifier based on a binomial test; the second was tumor variant amplitude (TVA), a continuous measure representing the selective advantage of individual variants. TVA outperformed all other computational tools in terms of its correlation with experimentally derived functional scores of cancer mutations. TVA also highly correlated with drug response, overall survival, and other clinical implications in relevant cancer genes. This study demonstrates how a selective advantage measure based on a large cancer dataset significantly impacts our understanding of the spectral effect of driver variants in cancer. The impact of this information will increase as cancer treatment becomes more precise and personalized to tumor-specific mutations.

SIGNIFICANCE

A new selective advantage estimation assists in oncogenic driver identification and relative effect measurements, enabling better prognostication, therapy selection, and prioritization.

Cancer fitness genes: emerging therapeutic targets for metastasis

Development of cancer therapeutics has traditionally focused on targeting driver oncogenes. Such an approach is limited by toxicity to normal tissues and treatment resistance. A class of 'cancer fitness genes' with crucial roles in metastasis have been identified. Elevated or altered activities of these genes do not directly cause cancer; instead, they relieve the stresses that tumor cells encounter and help them adapt to a changing microenvironment, thus facilitating tumor progression and metastasis. Importantly, as normal cells do not experience high levels of stress under physiological conditions, targeting cancer fitness genes is less likely to cause toxicity to noncancerous tissues. Here, we summarize the key features and function of cancer fitness genes and discuss their therapeutic potential.

Evolutionary Analysis of TCGA Data Using Over- and Under- Mutated Genes Identify Key Molecular Pathways and Cellular Functions in Lung Cancer Subtypes

We identify critical conserved and mutated genes through a theoretical model linking a gene’s fitness contribution to its observed mutational frequency in a clinical cohort. “Passenger” gene mutations do not alter fitness and have mutational frequencies determined by gene size and the mutation rate. Driver mutations, which increase fitness (and proliferation), are observed more frequently than expected. Non-synonymous mutations in essential genes reduce fitness and are eliminated by natural selection resulting in lower prevalence than expected. We apply this “evolutionary triage” principle to TCGA data from EGFR-mutant, KRAS-mutant, and NEK (non-EGFR/KRAS) lung adenocarcinomas. We find frequent overlap of evolutionarily selected non-synonymous gene mutations among the subtypes suggesting enrichment for adaptations to common local tissue selection forces. Overlap of conserved genes in the LUAD subtypes is rare suggesting negative evolutionary selection is strongly dependent on initiating mutational events during carcinogenesis. Highly expressed genes are more likely to be conserved and significant changes in expression (>20% increased/decreased) are common in genes with evolutionarily selected mutations but not in conserved genes. EGFR-mut cancers have fewer average mutations (89) than KRAS-mut (228) and NEK (313). Subtype-specific variation in conserved and mutated genes identify critical molecular components in cell signaling, extracellular matrix remodeling, and membrane transporters. These findings demonstrate subtype-specific patterns of co-adaptations between the defining driver mutation and somatically conserved genes as well as novel insights into epigenetic versus genetic contributions to cancer evolution.

Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

Loss of Epigenetic Regulation Disrupts Lineage Integrity, Induces Aberrant Alveogenesis, and Promotes Breast Cancer

Systematically investigating the scores of genes mutated in cancer and discerning disease drivers from inconsequential bystanders is a prerequisite for precision medicine but remains challenging. Here, we developed a somatic CRISPR/Cas9 mutagenesis screen to study 215 recurrent "long-tail" breast cancer genes, which revealed epigenetic regulation as a major tumor-suppressive mechanism. We report that components of the BAP1 and COMPASS-like complexes, including KMT2C/D, KDM6A, BAP1, and ASXL1/2 ("EpiDrivers"), cooperate with PIK3CAH1047R to transform mouse and human breast epithelial cells. Mechanistically, we find that activation of PIK3CAH1047R and concomitant EpiDriver loss triggered an alveolar-like lineage conversion of basal mammary epithelial cells and accelerated formation of luminal-like tumors, suggesting a basal origin for luminal tumors. EpiDriver mutations are found in ∼39% of human breast cancers, and ∼50% of ductal carcinoma in situ express casein, suggesting that lineage infidelity and alveogenic mimicry may significantly contribute to early steps of breast cancer etiology.

SIGNIFICANCE

Infrequently mutated genes comprise most of the mutational burden in breast tumors but are poorly understood. In vivo CRISPR screening identified functional tumor suppressors that converged on epigenetic regulation. Loss of epigenetic regulators accelerated tumorigenesis and revealed lineage infidelity and aberrant expression of alveogenesis genes as potential early events in tumorigenesis. This article is highlighted in the In This Issue feature, p. 2711.

LOF variants identifying candidate genes of laterality defects patients with congenital heart disease

Defects in laterality pattern can result in abnormal positioning of the internal organs during the early stages of embryogenesis, as manifested in heterotaxy syndrome and situs inversus, while laterality defects account for 3~7% of all congenital heart defects (CHDs). However, the pathogenic mechanism underlying most laterality defects remains unknown. In this study, we recruited 70 laterality defect patients with CHDs to identify candidate disease genes by exome sequencing. We then evaluated rare, loss-of-function (LOF) variants, identifying candidates by referring to previous literature. We chose TRIP11, DNHD1, CFAP74, and EGR4 as candidates from 776 LOF variants that met the initial screening criteria. After the variants-to-gene mapping, we performed function research on these candidate genes. The expression patterns and functions of these four candidate genes were studied by whole-mount in situ hybridization, gene knockdown, and gene rescue methods in zebrafish models. Among the four genes, trip11, dnhd1, and cfap74 morphant zebrafish displayed abnormalities in both cardiac looping and expression patterns of early signaling molecules, suggesting that these genes play important roles in the establishment of laterality patterns. Furthermore, we performed immunostaining and high-speed cilia video microscopy to investigate Kupffer's vesicle organogenesis and ciliogenesis of morphant zebrafish. Impairments of Kupffer's vesicle organogenesis or ciliogenesis were found in trip11, dnhd1, and cfap74 morphant zebrafish, which revealed the possible pathogenic mechanism of their LOF variants in laterality defects. These results highlight the importance of rare, LOF variants in identifying disease-related genes and identifying new roles for TRIP11, DNHD1, and CFAP74 in left-right patterning. Additionally, these findings are consistent with the complex genetics of laterality defects.

Neoplasia-associated Chromosome Translocations Resulting in Gene Truncation

Chromosomal translocations in cancer as well as benign neoplasias typically lead to the formation of fusion genes. Such genes may encode chimeric proteins when two protein-coding regions fuse in-frame, or they may result in deregulation of genes via promoter swapping or translocation of the gene into the vicinity of a highly active regulatory element. A less studied consequence of chromosomal translocations is the fusion of two breakpoint genes resulting in an out-of-frame chimera. The breaks then occur in one or both protein-coding regions forming a stop codon in the chimeric transcript shortly after the fusion point. Though the latter genetic events and mechanisms at first awoke little research interest, careful investigations have established them as neither rare nor inconsequential. In the present work, we review and discuss the truncation of genes in neoplastic cells resulting from chromosomal rearrangements, especially from seemingly balanced translocations.

MACHETE identifies interferon-encompassing chromosome 9p21.3 deletions as mediators of immune evasion and metastasis

The most prominent homozygous deletions in cancer affect chromosome 9p21.3 and eliminate CDKN2A/B tumor suppressors, disabling a cell-intrinsic barrier to tumorigenesis. Half of 9p21.3 deletions, however, also encompass a type I interferon (IFN) gene cluster; the consequences of this co-deletion remain unexplored. To functionally dissect 9p21.3 and other large genomic deletions, we developed a flexible deletion engineering strategy, MACHETE (molecular alteration of chromosomes with engineered tandem elements). Applying MACHETE to a syngeneic mouse model of pancreatic cancer, we found that co-deletion of the IFN cluster promoted immune evasion, metastasis and immunotherapy resistance. Mechanistically, IFN co-deletion disrupted type I IFN signaling in the tumor microenvironment, leading to marked changes in infiltrating immune cells and escape from CD8+ T-cell surveillance, effects largely driven by the poorly understood interferon epsilon. These results reveal a chromosomal deletion that disables both cell-intrinsic and cell-extrinsic tumor suppression and provide a framework for interrogating large deletions in cancer and beyond.

Deep Learning Approaches for Detection of Breast Adenocarcinoma Causing Carcinogenic Mutations

Genes are composed of DNA and each gene has a specific sequence. Recombination or replication within the gene base ends in a permanent change in the nucleotide collection in a DNA called mutation and some mutations can lead to cancer. Breast adenocarcinoma starts in secretary cells. Breast adenocarcinoma is the most common of all cancers that occur in women. According to a survey within the United States of America, there are more than 282,000 breast adenocarcinoma patients registered each 12 months, and most of them are women. Recognition of cancer in its early stages saves many lives. A proposed framework is developed for the early detection of breast adenocarcinoma using an ensemble learning technique with multiple deep learning algorithms, specifically: Long Short-Term Memory (LSTM), Gated Recurrent Units (GRU), and Bi-directional LSTM. There are 99 types of driver genes involved in breast adenocarcinoma. This study uses a dataset of 4127 samples including men and women taken from more than 12 cohorts of cancer detection institutes. The dataset encompasses a total of 6170 mutations that occur in 99 genes. On these gene sequences, different algorithms are applied for feature extraction. Three types of testing techniques including independent set testing, self-consistency testing, and a 10-fold cross-validation test is applied to validate and test the learning approaches. Subsequently, multiple deep learning approaches such as LSTM, GRU, and bi-directional LSTM algorithms are applied. Several evaluation metrics are enumerated for the validation of results including accuracy, sensitivity, specificity, Mathew’s correlation coefficient, area under the curve, training loss, precision, recall, F1 score, and Cohen’s kappa while the values obtained are 99.57, 99.50, 99.63, 0.99, 1.0, 0.2027, 99.57, 99.57, 99.57, and 99.14 respectively.

Chromatin accessibility of primary human cancers ties regional mutational processes and signatures with tissues of origin

Somatic mutations in cancer genomes are associated with DNA replication timing (RT) and chromatin accessibility (CA), however these observations are based on normal tissues and cell lines while primary cancer epigenomes remain uncharacterised. Here we use machine learning to model megabase-scale mutation burden in 2,500 whole cancer genomes and 17 cancer types via a compendium of 900 CA and RT profiles covering primary cancers, normal tissues, and cell lines. CA profiles of primary cancers, rather than those of normal tissues, are most predictive of regional mutagenesis in most cancer types. Feature prioritisation shows that the epigenomes of matching cancer types and organ systems are often the strongest predictors of regional mutation burden, highlighting disease-specific associations of mutational processes. The genomic distributions of mutational signatures are also shaped by the epigenomes of matched cancer and tissue types, with SBS5/40, carcinogenic and unknown signatures most accurately predicted by our models. In contrast, fewer associations of RT and regional mutagenesis are found. Lastly, the models highlight genomic regions with overrepresented mutations that dramatically exceed epigenome-derived expectations and show a pan-cancer convergence to genes and pathways involved in development and oncogenesis, indicating the potential of this approach for coding and non-coding driver discovery. The association of regional mutational processes with the epigenomes of primary cancers suggests that the landscape of passenger mutations is predominantly shaped by the epigenomes of cancer cells after oncogenic transformation.

Small extracellular vesicle DNA-mediated horizontal gene transfer as a driving force for tumor evolution: Facts and riddles

The basis of the conventional gene-centric view on tumor evolution is that vertically inherited mutations largely define the properties of tumor cells. In recent years, however, accumulating evidence shows that both the tumor cells and their microenvironment may acquire external, non-vertically inherited genetic properties via horizontal gene transfer (HGT), particularly through small extracellular vesicles (sEVs). Many phases of sEV-mediated HGT have been described, such as DNA packaging into small vesicles, their release, uptake by recipient cells, and incorporation of sEV-DNA into the recipient genome to modify the phenotype and properties of cells. Recent techniques in sEV separation, genome sequencing and editing, as well as the identification of new secretion mechanisms, shed light on a number of additional details of this phenomenon. Here, we discuss the key features of this form of gene transfer and make an attempt to draw relevant conclusions on the contribution of HGT to tumor evolution.

Pan-cancer analyses of synonymous mutations based on tissue-specific codon optimality

•

Developed tissue-specific codon optimality in 29 human tissues.

•

Applied these to analyze synonymous mutations in ∼10,000 tumor and normal samples.

•

Synonymous mutations frequently increase optimal codons in most cancer types.

•

Synonymous mutations frequently increase optimal codons cell cycle-related genes.

•

Frequency of optimal codon gain relates to proliferation, DDR deficiency, and survival.

Codon optimality has been demonstrated to be an important determinant of mRNA stability and expression levels in multiple model organisms and human cell lines. However, tissue-specific codon optimality has not been developed to investigate how codon optimality is usually perturbed by somatic synonymous mutations in human cancers. Here, we determined tissue-specific codon optimality in 29 human tissues based on mRNA expression data from the Genotype-Tissue Expression project. We found that optimal codons were associated with differentiation, whereas non-optimal codons were correlated with proliferation. Furthermore, codons biased toward differentiation displayed greater tissue specificity in codon optimality, and the tissue specificity of codon optimality was primarily present in amino acids with high degeneracy of the genetic code. By applying tissue-specific codon optimality to somatic synonymous mutations in 8532 tumor samples across 24 cancer types and to those in 416 normal cells across six human tissues, we found that synonymous mutations frequently increased optimal codons in tumor cells and cancer-related genes (e.g., genes involved in cell cycle). Furthermore, an elevated frequency of optimal codon gain was found to promote tumor cell proliferation in three cancer types characterized by DNA damage repair deficiency and could act as a prognostic biomarker for patients with triple-negative breast cancer. In summary, this study profiled tissue-specific codon optimality in human tissues, revealed alterations in codon optimality caused by synonymous mutations in human cancers, and highlighted the non-negligible role of optimal codon gain in tumorigenesis and therapeutics.

SYSMut: decoding the functional significance of rare somatic mutations in cancer

Current tailored-therapy efforts in cancer are largely focused on a small number of highly recurrently mutated driver genes but therapeutic targeting of these oncogenes remains challenging. However, the vast number of genes mutated infrequently across cancers has received less attention, in part, due to a lack of understanding of their biological significance. We present SYSMut, an extendable systems biology platform that can robustly infer the biologic consequences of somatic mutations by integrating routine multiomics profiles in primary tumors. We establish SYSMut's improved performance vis-à-vis state-of-the-art driver gene identification methodologies by recapitulating the functional impact of known driver genes, while additionally identifying novel functionally impactful mutated genes across 29 cancers. Subsequent application of SYSMut on low-frequency gene mutations in head and neck squamous cell (HNSC) cancers, followed by molecular and pharmacogenetic validation, revealed the lipidogenic network as a novel therapeutic vulnerability in aggressive HNSC cancers. SYSMut is thus a robust scalable framework that enables the discovery of new targetable avenues in cancer.

Pan-Cancer Analysis for Immune Cell Infiltration and Mutational Signatures Using Non-Negative Canonical Correlation Analysis

Mutational signatures indicate the mutational processes and substitution patterns in cancer cell genomes. However, the functional consequences of mutational signatures remain unclear, and there have been no comprehensive systematic studies to examine the relationships between the mutational signatures and the immune cell infiltration. Here, the relationship between mutational signatures and immune cell infiltration using non-negative canonical correlation analysis based on 8927 patients across 25 tumor types was investigated. By inspecting mutational signatures with the maximal coefficients determined by the non-negative canonical correlation analysis, the study identified mutational signatures related to immune cell infiltration composed of tumor microenvironments. The analysis was validated by showing that the genes associated with the identified mutational signatures were linked to overall survival by a Kaplan–Meier curve and a log-rank test and were mainly related to immunity by gene set enrichment analysis. These results will help expand our knowledge of tumor biology and recognize the functional roles and associations of immune systems with mutational signatures.

Inferring Potential Cancer Driving Synonymous Variants

Synonymous single nucleotide variants (sSNVs) are often considered functionally silent, but a few cases of cancer-causing sSNVs have been reported. From available databases, we collected four categories of sSNVs: germline, somatic in normal tissues, somatic in cancerous tissues, and putative cancer drivers. We found that screening sSNVs for recurrence among patients, conservation of the affected genomic position, and synVep prediction (synVep is a machine learning-based sSNV effect predictor) recovers cancer driver variants (termed proposed drivers) and previously unknown putative cancer genes. Of the 2.9 million somatic sSNVs found in the COSMIC database, we identified 2111 proposed cancer driver sSNVs. Of these, 326 sSNVs could be further tagged for possible RNA splicing effects, RNA structural changes, and affected RBP motifs. This list of proposed cancer driver sSNVs provides computational guidance in prioritizing the experimental evaluation of synonymous mutations found in cancers. Furthermore, our list of novel potential cancer genes, galvanized by synonymous mutations, may highlight yet unexplored cancer mechanisms.

Heterogeneity of the tumor immune microenvironment and its clinical relevance

During the course of tumorigenesis and subsequent metastasis, malignant cells gradually diversify and become more heterogeneous. Consequently, the tumor mass might be infiltrated by diverse immune-related components, including the cytokine/chemokine environment, cytotoxic activity, or immunosuppressive elements. This immunological heterogeneity is universally presented spatially or varies temporally along with tumor evolution or therapeutic intervention across almost all solid tumors. The heterogeneity of anti-tumor immunity shows a profound association with the progression of disease and responsiveness to treatment, particularly in the realm of immunotherapy. Therefore, an accurate understanding of tumor immunological heterogeneity is essential for the development of effective therapies. Facilitated by multi-regional and -omics sequencing, single cell sequencing, and longitudinal liquid biopsy approaches, recent studies have demonstrated the potential to investigate the complexity of immunological heterogeneity of the tumors and its clinical relevance in immunotherapy. Here, we aimed to review the mechanism underlying the heterogeneity of the immune microenvironment. We also explored how clinical assessments of tumor heterogeneity might facilitate the development of more effective personalized therapies.

SWAN pathway-network identification of common aneuploidy-based oncogenic drivers

Haploinsufficiency drives Darwinian evolution. Siblings, while alike in many aspects, differ due to monoallelic differences inherited from each parent. In cancer, solid tumors exhibit aneuploid genetics resulting in hundreds to thousands of monoallelic gene-level copy-number alterations (CNAs) in each tumor. Aneuploidy patterns are heterogeneous, posing a challenge to identify drivers in this high-noise genetic environment. Here, we developed Shifted Weighted Annotation Network (SWAN) analysis to assess biology impacted by cumulative monoallelic changes. SWAN enables an integrated pathway-network analysis of CNAs, RNA expression, and mutations via a simple web platform. SWAN is optimized to best prioritize known and novel tumor suppressors and oncogenes, thereby identifying drivers and potential druggable vulnerabilities within cancer CNAs. Protein homeostasis, phospholipid dephosphorylation, and ion transport pathways are commonly suppressed. An atlas of CNA pathways altered in each cancer type is released. These CNA network shifts highlight new, attractive targets to exploit in solid tumors.

The 2010-2020 National Natural Science Foundation of China: Radiation Therapy

OBJECTIVE

This study explains the application number and funding rate of oncology projects undertaken by the National Natural Science Foundation of China (NSFC), with focus on tumor radiotherapy-related research over the past 11 years.

METHODS

A stratified analysis was carried out on the application and funding status of tumor radiotherapy studies in different NSFC project categories, different research areas, and different tumor types. Research areas that required specific focus, such as immunology-related radiotherapy, multimodality imaging and radiomics, and post-radiotherapy organ injury, were separately analyzed.

RESULTS

The status and development trends of various related research fields were studied, and the research results were presented with the support of the NSFC, in order to provide reference for future applications and funding allocations.

CONCLUSION

The number of applications for funding increases every year. Although the total number of funded projects has also increased every year, the funding rate has decreased year by year. Projects on radiotherapy and immunization have been at the forefront in recent years, and the funding rate for these projects increases yearly.

Influence of T Cell-Mediated Immune Surveillance on Somatic Mutation Occurrences in Melanoma

Background

Neoantigens are presented on the cancer cell surface by peptide-restricted human leukocyte antigen (HLA) proteins and can subsequently activate cognate T cells. It has been hypothesized that the observed somatic mutations in tumors are shaped by immunosurveillance.

Methods

We investigated all somatic mutations identified in The Cancer Genome Atlas (TCGA) Skin Cutaneous Melanoma (SKCM) samples. By applying a computational algorithm, we calculated the binding affinity of the resulting neo-peptides and their corresponding wild-type peptides with the major histocompatibility complex (MHC) Class I complex. We then examined the relationship between binding affinity alterations and mutation frequency.

Results

Our results show that neoantigens derived from recurrent mutations tend to have lower binding affinities with the MHC Class I complex compared to peptides from non-recurrent mutations. Tumor samples harboring recurrent SKCM mutations exhibited lower immune infiltration levels, indicating a relatively colder immune microenvironment.

Conclusions

These results suggested that the occurrences of somatic mutations in melanoma have been shaped by immunosurveillance. Mutations that lead to neoantigens with high MHC class I binding affinity are more likely to be eliminated and thus are less likely to be present in tumors.

Incidental germline findings during molecular profiling of tumor tissues for precision oncology: molecular survey and methodological obstacles

Background

A fraction of patients referred for complex molecular profiling of biopsied tumors may harbor germline variants in genes associated with the development of hereditary cancer syndromes (HCS). Neither the bioinformatic analysis nor the reporting of such incidental germline findings are standardized.

Methods

Data from Next-Generation Sequencing (NGS) of biopsied tumor samples referred for complex molecular profiling were analyzed for germline variants in HCS-associated genes. Analysis of variant origin was performed employing bioinformatic algorithms followed by manual curation. When possible, the origin of the variant was validated by Sanger sequencing of the sample of normal tissue. The variants’ pathogenicity was assessed according to ACMG/AMP.

Results

Tumors were sampled from 183 patients (Males: 75 [41.0%]; Females: 108 [59.0%]; mean [SD] age, 57.7 [13.3] years) and analysed by targeted NGS. The most common tumor types were colorectal (19%), pancreatic (13%), and lung cancer (10%). A total of 56 sequence variants in genes associated with HCS were detected in 40 patients. Of them, 17 variants found in 14 patients were predicted to be of germline origin, with 6 variants interpreted as pathogenic (PV) or likely pathogenic (LPV), and 9 as variants of uncertain significance (VUS). For the 41 out of 42 (97%) missense variants in HCS-associated genes, the results of computational prediction of variant origin were concordant with that of experimental examination. We estimate that Sanger sequencing of a sample of normal tissue would be required for ~ 1–7% of the total assessed cases with PV or LPV, when necessity to follow with genetic counselling referral in ~ 2–15% of total assessed cases (PV, LPV or VUS found in HCS genes).

Conclusion

Incidental findings of pathogenic germline variants are common in data from cancer patients referred for complex molecular profiling. We propose an algorithm for the management of patients with newly detected variants in genes associated with HCS.

Dynamic Cancer Cell Heterogeneity: Diagnostic and Therapeutic Implications

Though heterogeneity of cancers is recognized and has been much discussed in recent years, the concept often remains overlooked in different routine examinations. Indeed, in clinical or biological articles, reviews, and textbooks, cancers and cancer cells are generally presented as evolving distinct entities rather than as an independent heterogeneous cooperative cell population with its self-oriented biology. There are, therefore, conceptual gaps which can mislead the interpretations/diagnostic and therapeutic approaches. In this short review, we wish to summarize and discuss various aspects of this dynamic evolving heterogeneity and its biological, pathological, clinical, diagnostic, and therapeutic implications, using thyroid carcinoma as an illustrative example.

Rapid testing of candidate oncogenes and tumour suppressor genes in signal transduction and neoplastic transformation

The COSMIC database (version 94) lists 576 genes in the Cancer Gene Census which have a defined function as drivers of malignancy (oncogenes) or as tumour suppressors (Tier 1). In addition, there are 147 genes with similar functions, but which are less well characterised (Tier 2). Furthermore, next-generation sequencing projects in the context of precision oncology activities are constantly discovering new ones. Since cancer genes differ from their wild-type precursors in numerous molecular and biochemical properties and exert significant differential effects on downstream processes, simple assays that can uncover oncogenic or anti-oncogenic functionality are desirable and may precede more sophisticated analyses. We describe simple functional assays for PTPN11 (protein-tyrosine phosphatase, non-receptor-type 11)/SHP2 mutants, which are typically found in RASopathies and exhibit potential oncogenic activity. We have also designed a functional test for lysyl oxidase (LOX), a prototypical class II tumour suppressor gene whose loss of function may contribute to neoplastic transformation by RAS oncogenes. Moreover, we applied this test to analyse three co-regulated, RAS-responsive genes for transformation-suppressive activity. The integration of these tests into systems biology studies will contribute to a better understanding of cellular networks in cancer.

Rearrangement-mediated cis-regulatory alterations in advanced patient tumors reveal interactions with therapy

The global impact of somatic structural variants (SVs) on gene regulation in advanced tumors with complex treatment histories has been mostly uncharacterized. Here, using whole-genome and RNA sequencing from 570 recurrent or metastatic tumors, we report the altered expression of hundreds of genes in association with nearby SV breakpoints, including oncogenes and G-protein-coupled receptor-related genes such as PLEKHG2. A significant fraction of genes with SV-expression associations correlate with worse patient survival in primary and advanced cancers, including SRD5A1. In many instances, SV-expression associations involve retrotransposons being translocated near genes. High overall SV burden is associated with treatment with DNA alkylating agents or taxanes and altered expression of metabolism-associated genes. SV-expression associations within tumors from topoisomerase I inhibitor-treated patients include chromatin-related genes. Within anthracycline-treated tumors, SV breakpoints near chromosome 1p genes include PDE4B. Patient treatment and history can help understand the widespread SV-mediated cis-regulatory alterations found in cancer.