Emerging landscape of oncogenic signatures across human cancers. Nat Genet. 2013;45:1127-1133. [PMID: 24071851 PMCID: PMC4320046 DOI: 10.1038/ng.2762]

The genomic trajectory of ovarian high-grade serous carcinoma can be observed in STIC lesions

Ovarian high-grade serous carcinoma (HGSC) originates in the fallopian tube, with secretory cells carrying a TP53 mutation, known as p53 signatures, identified as potential precursors. p53 signatures evolve into serous tubal intraepithelial carcinoma (STIC) lesions, which in turn progress into invasive HGSC, which readily spreads to the ovary and disseminates around the peritoneal cavity. We recently investigated the genomic landscape of early- and late-stage HGSC and found higher ploidy in late-stage (median 3.1) than early-stage (median 2.0) samples. Here, to explore whether the high ploidy and possible whole-genome duplication (WGD) observed in late-stage disease were determined early in the evolution of HGSC, we analysed archival formalin-fixed paraffin-embedded (FFPE) samples from five HGSC patients. p53 signatures and STIC lesions were laser-capture microdissected and sequenced using shallow whole-genome sequencing (sWGS), while invasive ovarian/fallopian tube and metastatic carcinoma samples underwent macrodissection and were profiled using both sWGS and targeted next-generation sequencing. Results showed highly similar patterns of global copy number change between STIC lesions and invasive carcinoma samples within each patient. Ploidy changes were evident in STIC lesions, but not p53 signatures, and there was a strong correlation between ploidy in STIC lesions and invasive ovarian/fallopian tube and metastatic samples in each patient. The reconstruction of sample phylogeny for each patient from relative copy number indicated that high ploidy, when present, occurred early in the evolution of HGSC, which was further validated by copy number signatures in ovarian and metastatic tumours. These findings suggest that aberrant ploidy, suggestive of WGD, arises early in HGSC and is detected in STIC lesions, implying that the trajectory of HGSC may be determined at the earliest stages of tumour development. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Application of PARP inhibitors combined with immune checkpoint inhibitors in ovarian cancer

The advent of polyadenosine diphosphate ribose polymerase inhibitors (PARPi) has brought about significant changes in the field of ovarian cancer treatment. However, in 2022, Rucaparib, Olaparib, and Niraparib, had their marketing approval revoked for third-line and subsequent therapies due to an increased potential for adverse events. Consequently, the exploration of new treatment modalities remains imperative. Recently, the integration of PARPi with immune checkpoint inhibitors (ICIs) has emerged as a potential remedy option within the context of ovarian cancer. This article offers a comprehensive examination of the mechanisms and applications of PARPi and ICIs in the treatment of ovarian cancer. It synthesizes the existing evidence supporting their combined use and discusses key considerations that merit attention in ongoing development efforts.

Overcoming Resistance in Cancer Therapy: Computational Exploration of PIK3CA Mutations, Unveiling Novel Non-Toxic Inhibitors, and Molecular Insights Into Targeting PI3Kα

Phosphoinositide-3-kinases (PI3 K) are pivotal regulators of cell signaling implicated in various cancers. Particularly, mutations in the PIK3CA gene encoding the p110α catalytic subunit drive oncogenic signaling, making it an attractive therapeutic target. Our study conducted in silico exploration of 31 PIK3CA mutations across breast, endometrial, colon, and ovarian cancers, assessing their impacts on response to PI3Kα inhibitors and identifying potential non-toxic inhibitors and also elucidating their effects on protein stability and flexibility. Specifically, we observed significant alterations in the stability and flexibility of the PI3 K protein induced by these mutations. Through molecular docking analysis, we evaluated the binding interactions between the selected inhibitors and the PI3 K protein. The filtration of ligands involved calculating chemical descriptors, incorporating Veber and Lipinski rules, as well as IC50 values and toxicity predictions. This process reduced the initial dataset of 1394 ligands to 12 potential non-toxic inhibitors, and four reference inhibitors with significant biological activity in clinical trials were then chosen based on their physico-chemical properties. This analysis revealed Lig5's exceptional performance, exhibiting superior affinity and specificity compared to established reference inhibitors such as pictilisib. Lig5 formed robust binding interactions with the PI3 K protein, suggesting its potential as a highly effective therapeutic agent against PI3 K-driven cancers. Furthermore, molecular dynamics simulations provided valuable insights into Lig5's stability and its interactions with PI3 K over 100 ns. These simulations supported Lig5's potential as a versatile inhibitor capable of effectively targeting various mutational profiles of PI3 K, thereby mitigating issues related to resistance and toxicity commonly associated with current inhibitors.

Functional exploration of copy number alterations in a Drosophila model of triple-negative breast cancer

Accounting for 10-20% of breast cancer cases, triple-negative breast cancer (TNBC) is associated with a disproportionate number of breast cancer deaths. One challenge in studying TNBC is its genomic profile: with the exception of TP53 loss, most breast cancer tumors are characterized by a high number of copy number alterations (CNAs), making modeling the disease in whole animals challenging. We computationally analyzed 186 CNA regions previously identified in breast cancer tumors to rank genes within each region by likelihood of acting as a tumor driver. We then used a Drosophila p53-Myc TNBC model to identify 48 genes as functional drivers. To demonstrate the utility of this functional database, we established six 3-hit models; altering candidate genes led to increased aspects of transformation as well as resistance to the chemotherapeutic drug fluorouracil. Our work provides a functional database of CNA-associated TNBC drivers, and a template for an integrated computational/whole-animal approach to identify functional drivers of transformation and drug resistance within CNAs in other tumor types.

Assessing RNA-Seq Workflow Methodologies Using Shannon Entropy

RNA-seq faces persistent challenges due to the ongoing, expanding array of data processing workflows, none of which have yet achieved standardization to date. It is imperative to determine which method most effectively preserves biological facts. Here, we used Shannon entropy as a tool for depicting the biological status of a system. Thus, we assessed the measurement of Shannon entropy by several RNA-seq workflow approaches, such as DESeq2 and edgeR, but also by combining nine normalization methods with log2 fold change on paired samples of TCGA RNA-seq representing datasets of 515 patients and spanning 12 different cancer types with 5-year overall survival rates ranging from 20% to 98%. Our analysis revealed that TPM, RLE, and TMM normalization, coupled with a threshold of log2 fold change ≥1, for identifying differentially expressed genes, yielded the best results. We propose that Shannon entropy can serve as an objective metric for refining the optimization of RNA-seq workflows and mRNA sequencing technologies.

Selection forces underlying aneuploidy patterns in cancer

Aneuploidy, the presence of an aberrant number of chromosomes, has been associated with tumorigenesis for over a century. More recently, advances in karyotyping techniques have revealed its high prevalence in cancer: About 90% of solid tumors and 50-70% of hematopoietic cancers exhibit chromosome gains or losses. When analyzed at the level of specific chromosomes, there are strong patterns that are observed in cancer karyotypes both pan-cancer and for specific cancer types. These specific aneuploidy patterns correlate strongly with outcomes for tumor initiation, progression, metastasis formation, immune evasion and resistance to therapeutic treatment. Despite their prominence, understanding the basis underlying aneuploidy patterns in cancer has been challenging. Advances in genetic engineering and bioinformatic analyses now offer insights into the genetic determinants of aneuploidy pattern selection. Overall, there is substantial evidence that expression changes of particular genes can act as the positive selective forces for adaptation through aneuploidy. Recent findings suggest that multiple genes contribute to the selection of specific aneuploid chromosomes in cancer; however, further research is necessary to identify the most impactful driver genes. Determining the genetic basis and accompanying vulnerabilities of specific aneuploidy patterns is an essential step in selectively targeting these hallmarks of tumors.

Tumor-Agnostic Therapy-The Final Step Forward in the Cure for Human Neoplasms?

Cancer accounted for 10 million deaths in 2020, nearly one in every six deaths annually. Despite advancements, the contemporary clinical management of human neoplasms faces a number of challenges. Surgical removal of tumor tissues is often not possible technically, while radiation and chemotherapy pose the risk of damaging healthy cells, tissues, and organs, presenting complex clinical challenges. These require a paradigm shift in developing new therapeutic modalities moving towards a more personalized and targeted approach. The tumor-agnostic philosophy, one of these new modalities, focuses on characteristic molecular signatures of transformed cells independently of their traditional histopathological classification. These include commonly occurring DNA aberrations in cancer cells, shared metabolic features of their homeostasis or immune evasion measures of the tumor tissues. The first dedicated, FDA-approved tumor-agnostic agent's profound progression-free survival of 78% in mismatch repair-deficient colorectal cancer paved the way for the accelerated FDA approvals of novel tumor-agnostic therapeutic compounds. Here, we review the historical background, current status, and future perspectives of this new era of clinical oncology.

Off-target autophagy inhibition by SHP2 allosteric inhibitors contributes to their antitumor activity in RAS-driven cancers

Aberrant activation of RAS/MAPK signaling is common in cancer, and efforts to inhibit pathway components have yielded drugs with promising clinical activities. Unfortunately, treatment-provoked adaptive resistance mechanisms inevitably develop, limiting their therapeutic potential. As a central node essential for receptor tyrosine kinase-mediated RAS activation, SHP2 has emerged as an attractive cancer target. Consequently, many SHP2 allosteric inhibitors are now in clinical testing. Here we discovered a previously unrecognized off-target effect associated with SHP2 allosteric inhibitors. We found that these inhibitors accumulate in the lysosome and block autophagic flux in an SHP2-independent manner. We showed that off-target autophagy inhibition by SHP2 allosteric inhibitors contributes to their antitumor activity. We also demonstrated that SHP2 allosteric inhibitors harboring this off-target activity not only suppress oncogenic RAS signaling but also overcome drug resistance such as MAPK rebound and protective autophagy in response to RAS/MAPK pathway blockage. Finally, we exemplified a therapeutic framework that harnesses both the on- and off-target activities of SHP2 allosteric inhibitors for improved treatment of mutant RAS-driven and drug-resistant malignancies such as pancreatic and colorectal cancers.

The yin and yang of chromosomal instability in prostate cancer

Metastatic prostate cancer remains an incurable lethal disease. Studies indicate that prostate cancer accumulates genomic changes during disease progression and displays the highest levels of chromosomal instability (CIN) across all types of metastatic tumours. CIN, which refers to ongoing chromosomal DNA gain or loss during mitosis, and derived aneuploidy, are known to be associated with increased tumour heterogeneity, metastasis and therapy resistance in many tumour types. Paradoxically, high CIN levels are also proposed to be detrimental to tumour cell survival, suggesting that cancer cells must develop adaptive mechanisms to ensure their survival. In the context of prostate cancer, studies indicate that CIN has a key role in disease progression and might also offer a therapeutic vulnerability that can be pharmacologically targeted. Thus, a comprehensive evaluation of the causes and consequences of CIN in prostate cancer, its contribution to aggressive advanced disease and a better understanding of the acquired CIN tolerance mechanisms can translate into new tumour classifications, biomarker development and therapeutic strategies.

Transcriptome analysis reveals PRKCA as a potential therapeutic target for overcoming cisplatin resistance in lung cancer through ferroptosis

Cisplatin-based chemotherapy is the current standard care for lung cancer patients; however, drug resistance frequently develops during treatment, thereby limiting therapeutic efficacy.The molecular mechanisms underlying cisplatin resistance remain elusive. In this study, we conducted an analysis of microarray data from the Gene Expression Omnibus (GEO) database under the accession numbers GSE21656, which encompassed expression profiling of cisplatin-resistant H460 (DDP-H460)and the parental cells (H460). Subsequently, we calculated the differentially expressed genes (DEGs) between DDP-H460 and H460. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of DEGs demonstrated significant impact on the Rap1, PI3K/AKT and MAPK signaling pathways. Moreover, protein and protein interaction (PPI) network analysis identified PRKCA, DET1, and UBE2N as hub genes that potentially contribute predominantly to cisplatin resistance. Ultimately, PRKCA was selected for validation due to its significant prognostic effect, which predicts unfavorable overall survival and disease-free survival in patients with lung cancer. Network analysis conducted on The Cancer Genome Atlas (TCGA) database revealed a strong gene-level correlation between PRKCA and TP53, CDKN2A, BYR2, TTN, KRAS, and PIK3CA; whereas at the protein level, it exhibited a high correlation with EGFR, Lck, Bcl2, and Syk. The in vitro experiments revealed that PRKCA was upregulated in the cisplatin-resistant A549 cells (DDP-A549), while knockdown of PRKCA increased DDP-A549 apoptosis upon cisplatin treatment. Moreover, we observed that PRKCA knockdown attenuated DDP-A549 proliferation, migration and invasion ability. Western blot analysis demonstrated that PRKCA knockdown downregulated phosphorylation of PI3K expression while upregulated the genes involved in ferroptosis signaling. In summary, our results elucidate the role of PRKCA in acquiring resistance to cisplatin and underscore its potential as a therapeutic target for cisplatin-resistant lung cancer.

Themis: advancing precision oncology through comprehensive molecular subtyping and optimization

Recent advances in tumor molecular subtyping have revolutionized precision oncology, offering novel avenues for patient-specific treatment strategies. However, a comprehensive and independent comparison of these subtyping methodologies remains unexplored. This study introduces 'Themis' (Tumor HEterogeneity analysis on Molecular subtypIng System), an evaluation platform that encapsulates a few representative tumor molecular subtyping methods, including Stemness, Anoikis, Metabolism, and pathway-based classifications, utilizing 38 test datasets curated from The Cancer Genome Atlas (TCGA) and significant studies. Our self-designed quantitative analysis uncovers the relative strengths, limitations, and applicability of each method in different clinical contexts. Crucially, Themis serves as a vital tool in identifying the most appropriate subtyping methods for specific clinical scenarios. It also guides fine-tuning existing subtyping methods to achieve more accurate phenotype-associated results. To demonstrate the practical utility, we apply Themis to a breast cancer dataset, showcasing its efficacy in selecting the most suitable subtyping methods for personalized medicine in various clinical scenarios. This study bridges a crucial gap in cancer research and lays a foundation for future advancements in individualized cancer therapy and patient management.

Adap-BDCM: Adaptive Bilinear Dynamic Cascade Model for Classification Tasks on CNV Datasets

Copy number variation (CNV) is an essential genetic driving factor of cancer formation and progression, making intelligent classification based on CNV feasible. However, there are a few challenges in the current machine learning and deep learning methods, such as the design of base classifier combination schemes in ensemble methods and the selection of layers of neural networks, which often result in low accuracy. Therefore, an adaptive bilinear dynamic cascade model (Adap-BDCM) is developed to further enhance the accuracy and applicability of these methods for intelligent classification on CNV datasets. In this model, a feature selection module is introduced to mitigate the interference of redundant information, and a bilinear model based on the gated attention mechanism is proposed to extract more beneficial deep fusion features. Furthermore, an adaptive base classifier selection scheme is designed to overcome the difficulty of manually designing base classifier combinations and enhance the applicability of the model. Lastly, a novel feature fusion scheme with an attribute recall submodule is constructed, effectively avoiding getting stuck in local solutions and missing some valuable information. Numerous experiments have demonstrated that our Adap-BDCM model exhibits optimal performance in cancer classification, stage prediction, and recurrence on CNV datasets. This study can assist physicians in making diagnoses faster and better.

Comprehensive profiling of cancer neoantigens from aberrant RNA splicing

BACKGROUND

Cancer neoantigens arise from protein-altering somatic mutations in tumor and rank among the most promising next-generation immuno-oncology agents when used in combination with immune checkpoint inhibitors. We previously developed a computational framework, REAL-neo, for identification, quality control, and prioritization of both class-I and class-II human leucocyte antigen (HLA)-presented neoantigens resulting from somatic single-nucleotide mutations, small insertions and deletions, and gene fusions. In this study, we developed a new module, SPLICE-neo, to identify neoantigens from aberrant RNA transcripts from two distinct sources: (1) DNA mutations within splice sites and (2) de novo RNA aberrant splicings.

METHODS

First, SPLICE-neo was used to profile all DNA splice-site mutations in 11,892 tumors from The Cancer Genome Atlas (TCGA) and identified 11 profiles of splicing donor or acceptor site gains or losses. Transcript isoforms resulting from the top seven most frequent profiles were computed using novel logic models. Second, SPLICE-neo identified de novo RNA splicing events using RNA sequencing reads mapped to novel exon junctions from either single, double, or multiple exon-skipping events. The aberrant transcripts from both sources were then ranked based on isoform expression levels and z-scores assuming that individual aberrant splicing events are rare. Finally, top-ranked novel isoforms were translated into protein, and the resulting neoepitopes were evaluated for neoantigen potential using REAL-neo. The top splicing neoantigen candidates binding to HLA-A*02:01 were validated using in vitro T2 binding assays.

RESULTS

We identified abundant splicing neoantigens in four representative TCGA cancers: BRCA, LUAD, LUSC, and LIHC. In addition to their substantial contribution to neoantigen load, several splicing neoantigens were potent tumor antigens with stronger bindings to HLA compared with the positive control of antigens from influenza virus.

CONCLUSIONS

SPLICE-neo is the first tool to comprehensively identify and prioritize splicing neoantigens from both DNA splice-site mutations and de novo RNA aberrant splicings. There are two major advances of SPLICE-neo. First, we developed novel logic models that assemble and prioritize full-length aberrant transcripts from DNA splice-site mutations. Second, SPLICE-neo can identify exon-skipping events involving more than two exons, which account for a quarter to one-third of all skipping events.

Evaluating the effectiveness of pre-operative diagnosis of ovarian cancer using minimally invasive liquid biopsies by combining serum human epididymis protein 4 and cell-free DNA in patients with an ovarian mass

OBJECTIVE

To assess the feasibility of scalable, objective, and minimally invasive liquid biopsy-derived biomarkers such as cell-free DNA copy number profiles, human epididymis protein 4 (HE4), and cancer antigen 125 (CA125) for pre-operative risk assessment of early-stage ovarian cancer in a clinically representative and diagnostically challenging population and to compare the performance of these biomarkers with the Risk of Malignancy Index (RMI).

METHODS

In this case-control study, we included 100 patients with an ovarian mass clinically suspected to be early-stage ovarian cancer. Of these 100 patients, 50 were confirmed to have a malignant mass (cases) and 50 had a benign mass (controls). Using WisecondorX, an algorithm used extensively in non-invasive prenatal testing, we calculated the benign-calibrated copy number profile abnormality score. This score represents how different a sample is from benign controls based on copy number profiles. We combined this score with HE4 serum concentration to separate cases and controls.

RESULTS

Combining the benign-calibrated copy number profile abnormality score with HE4, we obtained a model with a significantly higher sensitivity (42% vs 0%; p<0.002) at 99% specificity as compared with the RMI that is currently employed in clinical practice. Investigating performance in subgroups, we observed especially large differences in the advanced stage and non-high-grade serous ovarian cancer groups.

CONCLUSION

This study demonstrates that cell-free DNA can be successfully employed to perform pre-operative risk of malignancy assessment for ovarian masses; however, results warrant validation in a more extensive clinical study.

PRKDC-Mediated NHEJ May Play a Crucial Role in Aneuploidy of Chromosome 8-Driven Progression of Ovarian Cancer

High malignancy is a prominent characteristic of epithelial ovarian cancer (EOC), emphasizing the necessity for further elucidation of the potential mechanisms underlying cancer progression. Aneuploidy and copy number variation (CNV) partially contribute to the heightened malignancy observed in EOC; however, the precise features of aneuploidy and their underlying molecular patterns, as well as the relationship between CNV and aneuploidy in EOC, remain unclear. In this study, we employed single-cell sequencing data along with The Cancer Genome Atlas (TCGA) to investigate aneuploidy and CNV in EOC. The technique of fluorescence in situ hybridization (FISH) was employed using specific probes. The copy number variation within the genomic region of chromosome 8 (42754568-47889815) was assessed and utilized as a representative measure for the ploidy status of individual cells in chromosome 8. Differential expression analysis was performed between different subgroups based on chromosome 8 ploidy. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein-protein interaction (PPI), and hub-gene analyses were subsequently utilized to identify crucial genes involved. By classifying enriched tumor cells into distinct subtypes based on chromosome 8 ploidy combined with TCGA data integration, we identified key genes driving chromosome 8 aneuploidy in EOC, revealing that PRKDC gene involvement through the mediated non-homologous end-joining pathway may play a pivotal role in disease progression. Further validation through analysis of the GEO and TCGA database and survival assessment, considering both mRNA expression levels and CNV status of PRKDC, has confirmed its involvement in the progression of EOC. Further functional analysis revealed an upregulation of PRKDC in both ovarian EOC cells and tissues, with its expression showing a significant correlation with the extent of copy number variation (CNV) on chromosome 8. Taken together, CNV amplification and aneuploidy of chromosome 8 are important characteristics of EOC. PRKDC and the mediated NHEJ pathway may play a crucial role in driving aneuploidy on chromosome 8 during the progression of EOC.

Genetic interactions reveal distinct biological and therapeutic implications in breast cancer

Co-occurrence and mutual exclusivity of genomic alterations may reflect the existence of genetic interactions, potentially shaping distinct biological phenotypes and impacting therapeutic response in breast cancer. However, our understanding of them remains limited. Herein, we investigate a large-scale multi-omics cohort (n = 873) and a real-world clinical sequencing cohort (n = 4,405) including several clinical trials with detailed treatment outcomes and perform functional validation in patient-derived organoids, tumor fragments, and in vivo models. Through this comprehensive approach, we construct a network comprising co-alterations and mutually exclusive events and characterize their therapeutic potential and underlying biological basis. Notably, we identify associations between TP53mut-AURKAamp and endocrine therapy resistance, germline BRCA1mut-MYCamp and improved sensitivity to PARP inhibitors, and TP53mut-MYBamp and immunotherapy resistance. Furthermore, we reveal that precision treatment strategies informed by co-alterations hold promise to improve patient outcomes. Our study highlights the significance of genetic interactions in guiding genome-informed treatment decisions beyond single driver alterations.

Spatial transcriptomics reveals discrete tumour microenvironments and autocrine loops within ovarian cancer subclones

High-grade serous ovarian carcinoma (HGSOC) is genetically unstable and characterised by the presence of subclones with distinct genotypes. Intratumoural heterogeneity is linked to recurrence, chemotherapy resistance, and poor prognosis. Here, we use spatial transcriptomics to identify HGSOC subclones and study their association with infiltrating cell populations. Visium spatial transcriptomics reveals multiple tumour subclones with different copy number alterations present within individual tumour sections. These subclones differentially express various ligands and receptors and are predicted to differentially associate with different stromal and immune cell populations. In one sample, CosMx single molecule imaging reveals subclones differentially associating with immune cell populations, fibroblasts, and endothelial cells. Cell-to-cell communication analysis identifies subclone-specific signalling to stromal and immune cells and multiple subclone-specific autocrine loops. Our study highlights the high degree of subclonal heterogeneity in HGSOC and suggests that subclone-specific ligand and receptor expression patterns likely modulate how HGSOC cells interact with their local microenvironment.

CINner: modeling and simulation of chromosomal instability in cancer at single-cell resolution

Cancer development is characterized by chromosomal instability, manifesting in frequent occurrences of different genomic alteration mechanisms ranging in extent and impact. Mathematical modeling can help evaluate the role of each mutational process during tumor progression, however existing frameworks can only capture certain aspects of chromosomal instability (CIN). We present CINner, a mathematical framework for modeling genomic diversity and selection during tumor evolution. The main advantage of CINner is its flexibility to incorporate many genomic events that directly impact cellular fitness, from driver gene mutations to copy number alterations (CNAs), including focal amplifications and deletions, missegregations and whole-genome duplication (WGD). We apply CINner to find chromosome-arm selection parameters that drive tumorigenesis in the absence of WGD in chromosomally stable cancer types. We found that the selection parameters predict WGD prevalence among different chromosomally unstable tumors, hinting that the selective advantage of WGD cells hinges on their tolerance for aneuploidy and escape from nullisomy. Direct application of CINner to model the WGD proportion and fraction of genome altered (FGA) further uncovers the increase in CNA probabilities associated with WGD in each cancer type. CINner can also be utilized to study chromosomally stable cancer types, by applying a selection model based on driver gene mutations and focal amplifications or deletions. Finally, we used CINner to analyze the impact of CNA probabilities, chromosome selection parameters, tumor growth dynamics and population size on cancer fitness and heterogeneity. We expect that CINner will provide a powerful modeling tool for the oncology community to quantify the impact of newly uncovered genomic alteration mechanisms on shaping tumor progression and adaptation.

Integrated multiomic profiling of breast cancer in the Chinese population reveals patient stratification and therapeutic vulnerabilities

Molecular profiling guides precision treatment of breast cancer; however, Asian patients are underrepresented in publicly available large-scale studies. We established a comprehensive multiomics cohort of 773 Chinese patients with breast cancer and systematically analyzed their genomic, transcriptomic, proteomic, metabolomic, radiomic and digital pathology characteristics. Here we show that compared to breast cancers in white individuals, Asian individuals had more targetable AKT1 mutations. Integrated analysis revealed a higher proportion of HER2-enriched subtype and correspondingly more frequent ERBB2 amplification and higher HER2 protein abundance in the Chinese HR+HER2+ cohort, stressing anti-HER2 therapy for these individuals. Furthermore, comprehensive metabolomic and proteomic analyses revealed ferroptosis as a potential therapeutic target for basal-like tumors. The integration of clinical, transcriptomic, metabolomic, radiomic and pathological features allowed for efficient stratification of patients into groups with varying recurrence risks. Our study provides a public resource and new insights into the biology and ancestry specificity of breast cancer in the Asian population, offering potential for further precision treatment approaches.

Evaluation of AJCC staging system and proposal of a novel stage grouping system in retroperitoneal liposarcoma: the Fudan Zhongshan experience

Background

Overall survival (OS) varies significantly among individuals with heterogeneous retroperitoneal liposarcoma (RPLS), even among those with the same clinical stage. Improved staging of RPLS is a critical unmet need, given the disappointing results of external validations of the 8th American Joint Committee on Cancer (AJCC) TNM staging system.

Methods

The cohort study included 220 consecutive patients who underwent surgical resection for primary RPLS at the largest sarcoma centre of Fudan University in China from September 2009 to August 2021, combined with 277 adult patients with RPLS in the SEER database from 1975 to 2020. Data analysis was performed from December 2021 to December 2022. Patients were retrospectively restaged according to the 8th and 7th editions of the TNM staging system as well as the new TNM (nTNM) staging system. The primary endpoint was overall survival (OS). Comparative analysis of postoperative survival was performed using the Kaplan-Meier method, and differences between subgroups were tested using the log-rank test. The OS prediction nomogram was generated based on baseline variables and tumour characteristics. Harrell's consistency index (C-index), area under the curve (AUC) of receiver operating characteristic curves (ROC), and calibration curves were used to evaluate the performance of the nomogram.

Results

A total of 497 patients were enrolled in the study, including 282 (56.7%) male patients. The median follow-up was 51 months (interquartile range, IQR, 23-83), and the OS rates at 1, 3, and 5 years were 87.9%, 75.3%, and 64.9%, respectively. According to the staging distribution of the AJCC 7th edition, 6 patients were stage IA (1.2%), 189 patients were stage IB (38%), 12 patients were stage IIA (2.4%), 150 patients were stage IIB (30.1%), 131 patients were stage III (26.3%), and 9 patients were stage IV (1.8%). With the 8th edition staging, this distribution changed: 6 patients (1.2%) were stage IA, 189 patients (38%) were stage IB, 12 patients (2.4%) were stage II, 24 patients (4.8%) were stage IIIA, 257 patients (51.7%) were stage IIIB, and 9 patients (1.8%) were stage IV. 182 patients (36.6%) were reclassified according to the nTNM staging system with the new T stage classification. The C-index and log-rank score improved after implementation of nTNM implementation. The nTNM system was associated with improved identification of high-risk patients compared with the AJCC 7th and 8th TNM. The FNCLCC stage proved to be highly prognostic with significant intergroup differences in OS. The calibration curve shows a high degree of agreement between the actual OS rate and the nomogram estimated OS rate.

Conclusion

Compared with 8th AJCC TNM, 7th AJCC TNM staging system showed a more homogeneous staging distribution and a slight improvement in the prognostic accuracy of RPLS. The revised T-stage and nTNM systems showed better risk stratification performance. The FNCLCC stage was found to have high prognostic value, further emphasising histological grade is the least negligible prognostic factor in predicting patient survival. The constructed nomogram model enables individualized prognostic analysis and helps to develop risk-adapted therapy for RPLS patients.

eEF1A2 promotes PTEN-GSK3β-SCF complex-dependent degradation of Aurora kinase A and is inactivated in breast cancer

The translation elongation factor eEF1A promotes protein synthesis. Its methylation by METTL13 increases its activity, supporting tumor growth. However, in some cancers, a high abundance of eEF1A isoforms is associated with a good prognosis. Here, we found that eEF1A2 exhibited oncogenic or tumor-suppressor functions depending on its interaction with METTL13 or the phosphatase PTEN, respectively. METTL13 and PTEN competed for interaction with eEF1A2 in the same structural domain. PTEN-bound eEF1A2 promoted the ubiquitination and degradation of the mitosis-promoting Aurora kinase A in the S and G2 phases of the cell cycle. eEF1A2 bridged the interactions between the SKP1-CUL1-FBXW7 (SCF) ubiquitin ligase complex, the kinase GSK3β, and Aurora-A, thereby facilitating the phosphorylation of Aurora-A in a degron site that was recognized by FBXW7. Genetic ablation of Eef1a2 or Pten in mice resulted in a greater abundance of Aurora-A and increased cell cycling in mammary tumors, which was corroborated in breast cancer tissues from patients. Reactivating this pathway using fimepinostat, which relieves inhibitory signaling directed at PTEN and increases FBXW7 expression, combined with inhibiting Aurora-A with alisertib, suppressed breast cancer cell proliferation in culture and tumor growth in vivo. The findings demonstrate a therapeutically exploitable, tumor-suppressive role for eEF1A2 in breast cancer.

Mapping crossover events of mouse meiotic recombination by restriction fragment ligation-based Refresh-seq

Single-cell whole-genome sequencing methods have undergone great improvements over the past decade. However, allele dropout, which means the inability to detect both alleles simultaneously in an individual diploid cell, largely restricts the application of these methods particularly for medical applications. Here, we develop a new single-cell whole-genome sequencing method based on third-generation sequencing (TGS) platform named Refresh-seq (restriction fragment ligation-based genome amplification and TGS). It is based on restriction endonuclease cutting and ligation strategy in which two alleles in an individual cell can be cut into equal fragments and tend to be amplified simultaneously. As a new single-cell long-read genome sequencing method, Refresh-seq features much lower allele dropout rate compared with SMOOTH-seq. Furthermore, we apply Refresh-seq to 688 sperm cells and 272 female haploid cells (secondary polar bodies and parthenogenetic oocytes) from F1 hybrid mice. We acquire high-resolution genetic map of mouse meiosis recombination at low sequencing depth and reveal the sexual dimorphism in meiotic crossovers. We also phase the structure variations (deletions and insertions) in sperm cells and female haploid cells with high precision. Refresh-seq shows great performance in screening aneuploid sperm cells and oocytes due to the low allele dropout rate and has great potential for medical applications such as preimplantation genetic diagnosis.

scAbsolute: measuring single-cell ploidy and replication status

Cancer cells often exhibit DNA copy number aberrations and can vary widely in their ploidy. Correct estimation of the ploidy of single-cell genomes is paramount for downstream analysis. Based only on single-cell DNA sequencing information, scAbsolute achieves accurate and unbiased measurement of single-cell ploidy and replication status, including whole-genome duplications. We demonstrate scAbsolute's capabilities using experimental cell multiplets, a FUCCI cell cycle expression system, and a benchmark against state-of-the-art methods. scAbsolute provides a robust foundation for single-cell DNA sequencing analysis across different technologies and has the potential to enable improvements in a number of downstream analyses.

TMED2 promotes glioma tumorigenesis by being involved in EGFR recycling transport

Aberrant epidermal growth factor receptor (EGFR) signaling is the core signaling commonly activated in glioma. The transmembrane emp24 protein transport domain protein 2 (TMED2) interacts with cargo proteins involved in protein sorting and transport between endoplasmic reticulum (ER) and Golgi apparatus. In this study, we found the correlation between TMED2 with glioma progression and EGFR signaling through database analysis. Moreover, we demonstrated that TMED2 is essential for glioma cell proliferation, migration, and invasion at the cellular levels, as well as tumor formation in mouse models, underscoring its significance in the pathobiology of gliomas. Mechanistically, TMED2 was found to enhance EGFR-AKT signaling by facilitating EGFR recycling, thereby providing the initial evidence of TMED2's involvement in the membrane protein recycling process. In summary, our findings shed light on the roles and underlying mechanisms of TMED2 in the regulation of glioma tumorigenesis and EGFR signaling, suggesting that targeting TMED2 could emerge as a promising therapeutic strategy for gliomas and other tumors associated with aberrant EGFR signaling.

Tracing back primed resistance in cancer via sister cells

Exploring non-genetic evolution of cell states during cancer treatments has become attainable by recent advances in lineage-tracing methods. However, transcriptional changes that drive cells into resistant fates may be subtle, necessitating high resolution analysis. Here, we present ReSisTrace that uses shared transcriptomic features of sister cells to predict the states priming treatment resistance. Applying ReSisTrace in ovarian cancer cells perturbed with olaparib, carboplatin or natural killer (NK) cells reveals pre-resistant phenotypes defined by proteostatic and mRNA surveillance features, reflecting traits enriched in the upcoming subclonal selection. Furthermore, we show that DNA repair deficiency renders cells susceptible to both DNA damaging agents and NK killing in a context-dependent manner. Finally, we leverage the obtained pre-resistance profiles to predict and validate small molecules driving cells to sensitive states prior to treatment. In summary, ReSisTrace resolves pre-existing transcriptional features of treatment vulnerability, facilitating both molecular patient stratification and discovery of synergistic pre-sensitizing therapies.

Genome-wide quantification of copy-number aberration impact on gene expression in ovarian high-grade serous carcinoma

Copy-number alterations (CNAs) are a hallmark of cancer and can regulate cancer cell states via altered gene expression values. Herein, we have developed a copy-number impact (CNI) analysis method that quantifies the degree to which a gene expression value is impacted by CNAs and leveraged this analysis at the pathway level. Our results show that a high CNA is not necessarily reflected at the gene expression level, and our method is capable of detecting genes and pathways whose activity is strongly influenced by CNAs. Furthermore, the CNI analysis enables unbiased categorization of CNA categories, such as deletions and amplifications. We identified six CNI-driven pathways associated with poor treatment response in ovarian high-grade serous carcinoma (HGSC), which we found to be the most CNA-driven cancer across 14 cancer types. The key driver in most of these pathways was amplified wild-type KRAS, which we validated functionally using CRISPR modulation. Our results suggest that wild-type KRAS amplification is a driver of chemotherapy resistance in HGSC and may serve as a potential treatment target.

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.

The practical utility of AI-assisted molecular profiling in the diagnosis and management of cancer of unknown primary: an updated review

Cancer of unknown primary (CUP) presents a complex diagnostic challenge, characterized by metastatic tumors of unknown tissue origin and a dismal prognosis. This review delves into the emerging significance of artificial intelligence (AI) and machine learning (ML) in transforming the landscape of CUP diagnosis, classification, and treatment. ML approaches, trained on extensive molecular profiling data, have shown promise in accurately predicting tissue of origin. Genomic profiling, encompassing driver mutations and copy number variations, plays a pivotal role in CUP diagnosis by providing insights into tumor type-specific oncogenic alterations. Mutational signatures (MS), reflecting somatic mutation patterns, offer further insights into CUP diagnosis. Known MS with established etiology, such as ultraviolet (UV) light-induced DNA damage and tobacco exposure, have been identified in cases of dedifferentiated/transdifferentiated melanoma and carcinoma. Deep learning models that integrate gene expression data and DNA methylation patterns offer insights into tissue lineage and tumor classification. In digital pathology, machine learning algorithms analyze whole-slide images to aid in CUP classification. Finally, precision oncology, guided by molecular profiling, offers targeted therapies independent of primary tissue identification. Clinical trials assigning CUP patients to molecularly guided therapies, including targetable alterations and tumor mutation burden as an immunotherapy biomarker, have resulted in improved overall survival in a subset of patients. In conclusion, AI- and ML-driven approaches are revolutionizing CUP management by enhancing diagnostic accuracy. Precision oncology utilizing enhanced molecular profiling facilitates the identification of targeted therapies that transcend the need to identify the tissue of origin, ultimately improving patient outcomes.

Aneuploidy and complex genomic rearrangements in cancer evolution

Mutational processes that alter large genomic regions occur frequently in developing tumors. They range from simple copy number gains and losses to the shattering and reassembly of entire chromosomes. These catastrophic events, such as chromothripsis, chromoplexy and the formation of extrachromosomal DNA, affect the expression of many genes and therefore have a substantial effect on the fitness of the cells in which they arise. In this review, we cover large genomic alterations, the mechanisms that cause them and their effect on tumor development and evolution.

From Clustering to Cluster Explanations via Neural Networks

A recent trend in machine learning has been to enrich learned models with the ability to explain their own predictions. The emerging field of explainable AI (XAI) has so far mainly focused on supervised learning, in particular, deep neural network classifiers. In many practical problems, however, the label information is not given and the goal is instead to discover the underlying structure of the data, for example, its clusters. While powerful methods exist for extracting the cluster structure in data, they typically do not answer the question why a certain data point has been assigned to a given cluster. We propose a new framework that can, for the first time, explain cluster assignments in terms of input features in an efficient and reliable manner. It is based on the novel insight that clustering models can be rewritten as neural networks-or "neuralized." Cluster predictions of the obtained networks can then be quickly and accurately attributed to the input features. Several showcases demonstrate the ability of our method to assess the quality of learned clusters and to extract novel insights from the analyzed data and representations.

Somatic mouse models of gastric cancer reveal genotype-specific features of metastatic disease

Metastatic gastric carcinoma is a highly lethal cancer that responds poorly to conventional and molecularly targeted therapies. Despite its clinical relevance, the mechanisms underlying the behavior and therapeutic response of this disease are poorly understood owing, in part, to a paucity of tractable models. Here we developed methods to somatically introduce different oncogenic lesions directly into the murine gastric epithelium. Genotypic configurations observed in patients produced metastatic gastric cancers that recapitulated the histological, molecular and clinical features of all nonviral molecular subtypes of the human disease. Applying this platform to both wild-type and immunodeficient mice revealed previously unappreciated links between the genotype, organotropism and immune surveillance of metastatic cells, which produced distinct patterns of metastasis that were mirrored in patients. Our results establish a highly portable platform for generating autochthonous cancer models with flexible genotypes and host backgrounds, which can unravel mechanisms of gastric tumorigenesis or test new therapeutic concepts.

Intramedullary Spinal Cord Tumors: Whole-Genome Sequencing to Assist Management and Prognosis

Intramedullary spinal cord tumors (IMSCTs) harbor unique genetic mutations which may play a role in prognostication and management. To this end, we present the largest cohort of IMSCTs with genetic characterization in the literature from our multi-site institutional registry. A total of 93 IMSCT patient records were reviewed from the years 1999 to 2020. Out of these, 61 complied with all inclusion criteria, 14 of these patients had undergone genetic studies with 8 undergoing whole-genomic sequencing. Univariate analyses were used to assess any factors associated with progression-free survival (PFS) using the Cox proportional hazards model. Firth's penalized likelihood approach was used to account for the low event rates. Fisher's exact test was performed to compare whole-genome analyses and specific gene mutations with progression. PFS (months) was given as a hazard ratio. Only the absence of copy neutral loss of heterozygosity (LOH) was shown to be significant (0.05, p = 0.008). Additionally, higher risk of recurrence/progression was associated with LOH (p = 0.0179). Our results suggest LOH as a genetic predictor of shorter progression-free survival, particularly within ependymoma and glioblastoma tumor types. Further genomic research with larger multi-institutional datasets should focus on these mutations as possible prognostic factors.

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.

Human Fallopian Tube-Derived Organoids with TP53 and RAD51D Mutations Recapitulate an Early Stage High-Grade Serous Ovarian Cancer Phenotype In Vitro

RAD51D mutations have been implicated in the transformation of normal fallopian tube epithelial (FTE) cells into high-grade serous ovarian cancer (HGSOC), one of the most prevalent and aggressive gynecologic malignancies. Currently, no suitable model exists to elucidate the role of RAD51D in disease initiation and progression. Here, we established organoids from primary human FTE and introduced TP53 as well as RAD51D knockdown to enable the exploration of their mutational impact on FTE lesion generation. We observed that TP53 deletion rescued the adverse effects of RAD51D deletion on the proliferation, stemness, senescence, and apoptosis of FTE organoids. RAD51D deletion impaired the homologous recombination (HR) function and induced G2/M phase arrest, whereas concurrent TP53 deletion mitigated G0/G1 phase arrest and boosted DNA replication when combined with RAD51D mutation. The co-deletion of TP53 and RAD51D downregulated cilia assembly, development, and motility, but upregulated multiple HGSOC-associated pathways, including the IL-17 signaling pathway. IL-17A treatment significantly improved cell viability. TP53 and RAD51D co-deleted organoids exhibited heightened sensitivity to platinum, poly-ADP ribose polymerase inhibitors (PARPi), and cell cycle-related medication. In summary, our research highlighted the use of FTE organoids with RAD51D mutations as an invaluable in vitro platform for the early detection of carcinogenesis, mechanistic exploration, and drug screening.

KNSTRN Is a Prognostic Biomarker That Is Correlated with Immune Infiltration in Breast Cancer and Promotes Cell Cycle and Proliferation

Kinetochore-localized astrin/SPAG5-binding protein (KNSTRN) promotes the progression of bladder cancer and lung adenocarcinoma. However, its expression and biological function in breast cancer remain largely unknown. Therefore, this study aimed to analyze KNSTRN expression, prognoses, correlation with immune infiltration, expression-associated genes, and regulated signaling pathways to characterize its role in regulating the cell cycle using both bioinformatics and in vitro functional experiments. Analyses of The Cancer Genome Atlas, Gene Expression Omnibus, TIMER, and The Human Protein Atlas databases revealed a significant upregulation of KNSTRN transcript and protein levels in breast cancer. Kaplan-Meier survival analyses demonstrated a significant association between high expression of KNSTRN and poor overall survival, relapse-free survival, post-progression survival, and distant metastases-free survival in patients with breast cancer. Furthermore, multivariate Cox regression analyses confirmed that KNSTRN is an independent prognostic factor for breast cancer. Immune infiltration analysis indicated a positive correlation between KNSTRN expression and T regulatory cell infiltration while showing a negative correlation with Tgd and natural killer cell infiltration. Gene set enrichment analysis along with single-cell transcriptome data analysis suggested that KNSTRN promoted cell cycle progression by regulating the expression of key cell cycle proteins. The overexpression and silencing of KNSTRN in vitro, respectively, promoted and inhibited the proliferation of breast cancer cells. The overexpression of KNSTRN enhanced the expression of key cell cycle regulators, including CDK4, CDK6, and cyclin D3, thereby accelerating the G1/S phase transition and leading to aberrant proliferation of breast cancer cells. In conclusion, our study demonstrates that KNSTRN functions as an oncogene in breast cancer by regulating immune response, promoting G1/S transition, and facilitating breast cancer cell proliferation. Moreover, KNSTRN has potential as a molecular biomarker for diagnostic and prognostic prediction in breast cancer.

Biology-driven therapy advances in high-grade serous ovarian cancer

Following a period of slow progress, the completion of genome sequencing and the paradigm shift relative to the cell of origin for high grade serous ovarian cancer (HGSOC) led to a new perspective on the biology and therapeutic solutions for this deadly cancer. Experimental models were revisited to address old questions, and improved tools were generated. Additional pathways emerging as drivers of ovarian tumorigenesis and key dependencies for therapeutic targeting, in particular, VEGF-driven angiogenesis and homologous recombination deficiency, were discovered. Molecular profiling of histological subtypes of ovarian cancer defined distinct genetic events for each entity, enabling the first attempts toward personalized treatment. Armed with this knowledge, HGSOC treatment was revised to include new agents. Among them, PARP inhibitors (PARPis) were shown to induce unprecedented improvement in clinical benefit for selected subsets of patients. Research on mechanisms of resistance to PARPis is beginning to discover vulnerabilities and point to new treatment possibilities. This Review highlights these advances, the remaining challenges, and unsolved problems in the field.

Exploring the basis of heterogeneity of cancer aggressiveness among the mutated POLE variants

Germline pathogenic variants in the exonuclease domain of the replicative DNA polymerase Pol ε encoded by the POLE gene, predispose essentially to colorectal and endometrial tumors by inducing an ultramutator phenotype. It is still unclear whether all the POLE alterations influence similar strength tumorigenesis, immune microenvironment, and treatment response. In this review, we summarize the current understanding of the mechanisms and consequences of POLE mutations in human malignancies; we highlight the heterogeneity of mutation rate and cancer aggressiveness among POLE variants, propose some mechanistic basis underlining such heterogeneity, and discuss novel considerations for the choice and efficacy of therapies of POLE tumors.

Adaptive Control of Tumor Growth

Cancer treatment optimizations select the most optimum combinations of drugs, sequencing schedules, and appropriate doses that would limit toxicity and yield an improved patient quality of life. However, these optimizations often lack an adequate consideration of cancer's near-infinite potential for evolutionary adaptation to therapeutic interventions. Adapting cancer therapy based on monitored tumor burden and clonal composition is an intuitively sound approach to the treatment of cancer as an inherently complex and adaptive system. The adaptation would be driven by clinical outcome setpoints embodying the aims to thwart therapeutic resistance and maintain a long-term management of the disease or even a cure. However, given the nonlinear, stochastic dynamics of tumor response to therapeutic interventions, adaptive therapeutic strategies may at least need a one-step-ahead prediction of tumor burden to maintain their control over tumor growth dynamics. The article explores the feasibility of adaptive cancer treatment driven by tumor state feedback assuming cell adaptive fitness to be the underlying source of phenotypic plasticity and pathway entropy as a biomarker of tumor growth trajectory. The exploration is undertaken using deterministic and stochastic models of tumor growth dynamics.

Genomic alterations driving precancerous to cancerous lesions in esophageal cancer development

Esophageal squamous cell carcinoma (ESCC) develops through a series of increasingly abnormal precancerous lesions. Previous studies have revealed the striking differences between normal esophageal epithelium and ESCC in copy number alterations (CNAs) and mutations in genes driving clonal expansion. However, due to limited data on early precancerous lesions, the timing of these transitions and which among them are prerequisites for malignant transformation remained unclear. Here, we analyze 1,275 micro-biopsies from normal esophagus, early and late precancerous lesions, and esophageal cancers to decipher the genomic alterations at each stage. We show that the frequency of TP53 biallelic inactivation increases dramatically in early precancerous lesion stage while CNAs and APOBEC mutagenesis substantially increase at late stages. TP53 biallelic loss is the prerequisite for the development of CNAs of genes in cell cycle, DNA repair, and apoptosis pathways, suggesting it might be one of the earliest steps initiating malignant transformation.

Nuclear Glycoprotein A Repetitions Predominant (GARP) Is a Common Trait of Glioblastoma Stem-like Cells and Correlates with Poor Survival in Glioblastoma Patients

Glioblastoma (GB) is notoriously resistant to therapy. GB genesis and progression are driven by glioblastoma stem-like cells (GSCs). One goal for improving treatment efficacy and patient outcomes is targeting GSCs. Currently, there are no universal markers for GSCs. Glycoprotein A repetitions predominant (GARP), an anti-inflammatory protein expressed by activated regulatory T cells, was identified as a possible marker for GSCs. This study evaluated GARP for the detection of human GSCs utilizing a multidimensional experimental design that replicated several features of GB: (1) intratumoral heterogeneity, (2) cellular hierarchy (GSCs with varied degrees of self-renewal and differentiation), and (3) longitudinal GSC evolution during GB recurrence (GSCs from patient-matched newly diagnosed and recurrent GB). Our results indicate that GARP is expressed by GSCs across various cellular states and disease stages. GSCs with an increased GARP expression had reduced self-renewal but no alterations in proliferative capacity or differentiation commitment. Rather, GARP correlated inversely with the expression of GFAP and PDGFR-α, markers of astrocyte or oligodendrocyte differentiation. GARP had an abnormal nuclear localization (GARPNU+) in GSCs and was negatively associated with patient survival. The uniformity of GARP/GARPNU+ expression across different types of GSCs suggests a potential use of GARP as a marker to identify GSCs.

Clinical implication of genetic intratumor heterogeneity for targeted therapy in head and neck cancer

BACKGROUND

Genomic profiling is increasingly used both in therapeutic decision-making and as inclusion criteria for trials testing targeted therapies. However, the mutational landscape may vary across different areas of a tumor and intratumor heterogeneity will challenge treatments or clinical decisions based on single tumor biopsies. The purpose of this study was to assess the clinical relevance of genetic intratumor heterogeneity in head and neck squamous cell carcinomas (HNSCC) using the ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT).

MATERIALS AND METHODS

This prospective study included 33 whole tumor specimens from 28 patients with primary or recurrent HNSCC referred for surgery. Three tumor blocks were selected from central, semi-peripheral, and peripheral positions, mimicking biopsies in three different locations. Genetic analysis of somatic copy number alterations (SCNAs) was performed on the three biopsies using Oncoscan, focusing on 45 preselected HNSCC genes of interest. Clinical relevance was assessed using the ESCAT score to investigate whether and how treatment decisions would change based on the three biopsies from the same tumor.

RESULTS

The SCNAs identified among 45 preselected genes within the three tumor biopsies derived from the same tumor revealed distinct variations. The detected discrepancies could potentially influence treatment approaches or clinical decisions in 36% of the patients if only one tumor biopsy was used. Recurrent tumors exhibited significantly higher variation in SCNAs than primary tumors (p = .024). No significant correlation between tumor size and heterogeneity (p = .7) was observed.

CONCLUSION

In 36% of patients diagnosed with HNSCC, clinically significant intratumor heterogeneity was observed which may have implications for patient management. This finding substantiates the need for future studies that specifically investigate the clinical implications associated with intratumor heterogeneity.

Classification of tumor types using XGBoost machine learning model: a vector space transformation of genomic alterations

BACKGROUND

Machine learning (ML) represents a powerful tool to capture relationships between molecular alterations and cancer types and to extract biological information. Here, we developed a plain ML model aimed at distinguishing cancer types based on genetic lesions, providing an additional tool to improve cancer diagnosis, particularly for tumors of unknown origin.

METHODS

TCGA data from 9,927 samples spanning 32 different cancer types were downloaded from cBioportal. A vector space model type data transformation technique was designed to build consistently homogeneous new datasets containing, as predictive features, calls for somatic point mutations and copy number variations at chromosome arm-level, thus allowing the use of the XGBoost classifier models. Considering the imbalance in the dataset, due to large difference in the number of cases for each tumor, two preprocessing strategies were considered: i) setting a percentage cut-off threshold to remove less represented cancer types, ii) dividing cancer types into different groups based on biological criteria and training a specific XGBoost model for each of them. The performance of all trained models was mainly assessed by the out-of-sample balanced accuracy (BACC) and the AUC scores.

RESULTS

The XGBoost classifier achieved the best performance (BACC 77%; AUC 97%) on a dataset containing the 10 most represented tumor types. Moreover, dividing the 18 most represented cancers into three different groups (endocrine-related carcinomas, other carcinomas and other cancers),such analysis models achieved 78%, 71% and 86% BACC, respectively, with AUC scores greater than 96%. In addition, the model capable of linking each group to a specific cancer type reached 81% BACC and 94% AUC. Overall, the diagnostic potential of our model was comparable/higher with respect to others already described in literature and based on similar molecular data and ML approaches.

CONCLUSIONS

A boosted ML approach able to accurately discriminate different cancer types was developed. The methodology builds datasets simpler and more interpretable than the original data, while keeping enough information to accurately train standard ML models without resorting to sophisticated Deep Learning architectures. In combination with histopathological examinations, this approach could improve cancer diagnosis by using specific DNA alterations, processed by a replicable and easy-to-use automated technology. The study encourages new investigations which could further increase the classifier's performance, for example by considering more features and dividing tumors into their main molecular subtypes.

The Pathology according to p53 Pathway

BACKGROUND

Observations play a pivotal role in the progress of science, including in pathology. The cause of a disease such as cancer is analyzed by breaking it down into smaller organs, tissues, cells, and molecules. The current standard cancer diagnostic procedure, microscopic observation, relies on preserved morphological characteristics. In contrast, molecular analyses explore oncogenic pathway activation that leads to genetic mutations and aberrant protein expression. Such molecular analyses could potentially identify therapeutic targets and has gained considerable attention in clinical oncology.

SUMMARY

This review summarizes the cardinal biomarkers of the p53 pathway, p53, p16, and mouse double minute 2 (MDM2), in the context of traditional surgical pathology and emerging genomic oncology. The p53 pathway, which is dysregulated in more than a half of all cancers, can be applied in several diagnostic settings. A four-classification model of immunophenotype for p53 pathway gene status, tumor types with a high frequency of abnormalities for each p53 pathway gene, and a minimal p53 pathway immunohistochemical panel is also described.

KEY MESSAGES

Immunohistochemistry of oncogenic signals should be interpreted according to molecular findings based on genomic oncology, in addition to the microscopic findings of diagnostic pathology.

Oncogenic pathways refine a new perspective on the classification of hepatocellular carcinoma

BACKGROUND

Genetic alterations in oncogenic pathways are critical for cancer initiation, development, and treatment resistance. However, studies are limited regarding pathways correlated with prognosis, sorafenib, and transcatheter arterial chemoembolization (TACE) in hepatocellular carcinoma (HCC).

METHODS

In this study, 1928 patients from 11 independent datasets and a clinical in-house cohort were screened to explore the relationships among canonical pathway alterations in HCC patients. The molecular mechanisms, biological functions, immune landscape, and clinical outcomes among three heterogeneous phenotypes were further explored.

RESULTS

We charted the detailed landscape of pathway alterations in the TCGA-LIHC cohort, screened three pivotal pathways (p53, PI3K, and WNT), identified co-occurrence patterns and mutual exclusively, and stratified patients into three altered-pathway dominant phenotypes (ADPs). P53|PI3K ADP characterized by genomic instability (e.g., highest TMB, FGA, FGG, and FGL) indicated an unfavorable prognosis. While, patients in WNT ADP suggested a median prognosis, enhanced immune activation, and sensitivity to PD-L1 therapy. Remarkably, sorafenib and TACE exhibited efficacy for patients in WNT ADP and low frequent alteration phenotype (LFP). Additionally, ADP could work independently of common clinical traits (e.g., AJCC stage) and previous molecular classifications (e.g., iCluster, serum biomarkers).

CONCLUSIONS

ADP provides a new perspective for identifying patients at high risk of recurrence and could optimize precision treatment to improve the clinical outcomes in HCC.

Permission to pass: on the role of p53 as a gatekeeper for aneuploidy

Aneuploidy-the karyotype state in which the number of chromosomes deviates from a multiple of the haploid chromosome set-is common in cancer, where it is thought to facilitate tumor initiation and progression. However, it is poorly tolerated in healthy cells: during development and tissue homeostasis, aneuploid cells are efficiently cleared from the population. It is still largely unknown how cancer cells become, and adapt to being, aneuploid. P53, the gatekeeper of the genome, has been proposed to guard against aneuploidy. Aneuploidy in cancer genomes strongly correlates with mutations in TP53, and p53 is thought to prevent the propagation of aneuploid cells. Whether p53 also participates in preventing the mistakes in cell division that lead to aneuploidy is still under debate. In this review, we summarize the current understanding of the role of p53 in protecting cells from aneuploidy, and we explore the consequences of functional p53 loss for the propagation of aneuploidy in cancer.

Performance of computational algorithms to deconvolve heterogeneous bulk ovarian tumor tissue depends on experimental factors

BACKGROUND

Single-cell gene expression profiling provides unique opportunities to understand tumor heterogeneity and the tumor microenvironment. Because of cost and feasibility, profiling bulk tumors remains the primary population-scale analytical strategy. Many algorithms can deconvolve these tumors using single-cell profiles to infer their composition. While experimental choices do not change the true underlying composition of the tumor, they can affect the measurements produced by the assay.

RESULTS

We generated a dataset of high-grade serous ovarian tumors with paired expression profiles from using multiple strategies to examine the extent to which experimental factors impact the results of downstream tumor deconvolution methods. We find that pooling samples for single-cell sequencing and subsequent demultiplexing has a minimal effect. We identify dissociation-induced differences that affect cell composition, leading to changes that may compromise the assumptions underlying some deconvolution algorithms. We also observe differences across mRNA enrichment methods that introduce additional discrepancies between the two data types. We also find that experimental factors change cell composition estimates and that the impact differs by method.

CONCLUSIONS

Previous benchmarks of deconvolution methods have largely ignored experimental factors. We find that methods vary in their robustness to experimental factors. We provide recommendations for methods developers seeking to produce the next generation of deconvolution approaches and for scientists designing experiments using deconvolution to study tumor heterogeneity.

Molecular landscape and functional characterization of centrosome amplification in ovarian cancer

High-grade serous ovarian carcinoma (HGSOC) is characterised by poor outcome and extreme chromosome instability (CIN). Therapies targeting centrosome amplification (CA), a key mediator of chromosome missegregation, may have significant clinical utility in HGSOC. However, the prevalence of CA in HGSOC, its relationship to genomic biomarkers of CIN and its potential impact on therapeutic response have not been defined. Using high-throughput multi-regional microscopy on 287 clinical HGSOC tissues and 73 cell lines models, here we show that CA through centriole overduplication is a highly recurrent and heterogeneous feature of HGSOC and strongly associated with CIN and genome subclonality. Cell-based studies showed that high-prevalence CA is phenocopied in ovarian cancer cell lines, and that high CA is associated with increased multi-treatment resistance; most notably to paclitaxel, the commonest treatment used in HGSOC. CA in HGSOC may therefore present a potential driver of tumour evolution and a powerful biomarker for response to standard-of-care treatment.

KEAP1/NRF2 Mutations in Stem Cells Define an Aggressive Subset of Head and Neck Cancer Patients Who Have a Poor Prognosis, Lung Metastasis, and Therapeutic Failure

Mutations in Keap1/Nrf2 in head and neck cancer result in abnormal cell growth. Progenitor cells, bulk tumor cells, and head and neck cancer stem cells (HN-CSCs) may all harbor these mutations. Nevertheless, whether Keap1/Nrf2 mutations in HN-CSCs have an impact on clinical outcomes is unknown. Cancerous HN-CSCs and benign stem cells were obtained from freshly resected head and neck cancer patients (n = 50) via flow cytometry cell sorting and tested for Keap1/Nrf2 mutations. The existence of Keap1/Nrf2 mutations in HN-CSCs, as well as their correlations with tumor mutations, pathologic tumor stage, tumor histologic grades, lung metastasis, treatment outcomes, and the patient's age and conditions, are assessed at the last follow-up visit. Thirteen tumors were found to have Keap1/Nrf2 mutations in their HN-CSCs. More than half of the lung metastases and disease progression occurred in HN-CSCs with mutations. Patients whose tumors carried Keap1/Nrf2 mutations in their HN-CSCs had significantly shorter progression-free survival, overall survival, and time of treatment failure than their non-HN-CSC counterparts. These associations were partly driven by HN-CSCs, in which Keap1/Nrf2 mutations were overrepresented in fast progressors and associated with an increased risk of disease progression. Our findings suggest that molecular genotyping of HN-CSCs may facilitate personalized treatment strategies and assist in identifying patients who are likely to benefit from chemotherapy.

HCK Promotes High-Grade Serous Ovarian Cancer Tumorigenesis through CD44 and NOTCH3 Signaling

High-grade serous ovarian cancer (HGSOC) is a highly aggressive and lethal subtype of ovarian cancer. While most patients initially respond to standard-of-care treatment, the majority will eventually relapse and succumb to their disease. Despite significant advances in our understanding of this disease, the mechanisms that govern the distinctions between HGSOC with good and poor prognosis remain unclear. In this study, we implemented a proteogenomic approach to analyze gene expression, proteomic and phosphoproteomic profiles of HGSOC tumor samples to identify molecular pathways that distinguish HGSOC tumors relative to clinical outcome. Our analyses identify significant upregulation of hematopoietic cell kinase (HCK) expression and signaling in poor prognostic HGSOC patient samples. Analyses of independent gene expression datasets and IHC of patient samples confirmed increased HCK signaling in tumors relative to normal fallopian or ovarian samples and demonstrated aberrant expression in tumor epithelial cells. Consistent with the association between HCK expression and tumor aggressiveness in patient samples, in vitro phenotypic studies showed that HCK can, in part, promote cell proliferation, colony formation, and invasive capacity of cell lines. Mechanistically, HCK mediates these phenotypes, partly through CD44 and NOTCH3-dependent signaling, and inhibiting CD44 or NOTCH3 activity, either genetically or through gamma-secretase inhibitors, can revert HCK-driven phenotypes.

IMPLICATIONS

Collectively, these studies establish that HCK acts as an oncogenic driver of HGSOC through aberrant activation of CD44 and NOTCH3 signaling and identifies this network as a potential therapeutic opportunity in a subset of patients with aggressive and recurrent HGSOC.

Predicting the likelihood of a BRCA1/2 pathogenic variant being somatic by testing only tumour DNA in non-mucinous high-grade epithelial ovarian cancer

AIMS

Clinical guidelines recommend testing both germline and tumour DNA for BRCA1/2 pathogenic variants (PVs) in non-mucinous high-grade epithelial ovarian cancer (NMEOC). In this study, we show that some tumour BRCA1/2 PVs are highly likely to be somatic based on certain clinical and variant characteristics, meaning it may not be necessary to test all NMEOC cases for germline BRCA1/2 PVs.

METHODS

An observational study that included all tumour BRCA1/2 PVs detected in cases of NMEOC in the Northwest of England between July 2017 and February 2022. All tumour BRCA1/2 PVs were compared with PVs recorded in a prospectively gathered pan-cancer germline BRCA1/2 (gBRCA) testing database for the same geographical region (gBRCA1 PVs=910 and gBRCA2 PVs=922). Tumour BRCA1/2 PVs were categorised as common (≥1%), uncommon (<1%) or absent from the germline database.

RESULTS

One hundred and thirteen tumour BRCA1/2 PVs were detected in 111 NMEOC cases. There were 69 germline and 44 somatic variants. The mean age at diagnosis for gBRCA and somatic BRCA1/2 (sBRCA) PVs was 56.9 and 68.5 years, respectively (Student's t-test p<0.0001). All sBRCA PVs were detected in non-familial cases. All tumour BRCA1/2 PVs with a variant allele frequency (VAF) <35% in non-familial cases were somatic variants. Eighty-one per cent of germline-tumour BRCA1/2 PVs were present (common=31, uncommon=25) in the gBRCA testing database, while 89% of somatic-tumour BRCA1/2 PVs were absent (n=39).

CONCLUSIONS

We predict the likelihood of a tumour BRCA1/2 PV being somatic is 99.8% in non-familial cases of NMEOC diagnosed aged ≥75, where the VAF is ≤30% and there is no regional germline commonality.