SciClone: inferring clonal architecture and tracking the spatial and temporal patterns of tumor evolution

Integrative spatial and genomic analysis of tumor heterogeneity with Tumoroscope

Spatial and genomic heterogeneity of tumors are crucial factors influencing cancer progression, treatment, and survival. However, a technology for direct mapping the clones in the tumor tissue based on somatic point mutations is lacking. Here, we propose Tumoroscope, the first probabilistic model that accurately infers cancer clones and their localization in close to single-cell resolution by integrating pathological images, whole exome sequencing, and spatial transcriptomics data. In contrast to previous methods, Tumoroscope explicitly addresses the problem of deconvoluting the proportions of clones in spatial transcriptomics spots. Applied to a reference prostate cancer dataset and a newly generated breast cancer dataset, Tumoroscope reveals spatial patterns of clone colocalization and mutual exclusion in sub-areas of the tumor tissue. We further infer clone-specific gene expression levels and the most highly expressed genes for each clone. In summary, Tumoroscope enables an integrated study of the spatial, genomic, and phenotypic organization of tumors.

Whole-exome profiles of inflammatory breast cancer and pathological response to neoadjuvant chemotherapy

BACKGROUND

Neoadjuvant chemotherapy (NACT) became a standard treatment strategy for patients with inflammatory breast cancer (IBC) because of high disease aggressiveness. However, given the heterogeneity of IBC, no molecular feature reliably predicts the response to chemotherapy. Whole-exome sequencing (WES) of clinical tumor samples provides an opportunity to identify genomic alterations associated with chemosensitivity.

METHODS

We retrospectively applied WES to 44 untreated IBC primary tumor samples and matched normal DNA. The pathological response to NACT, assessed on operative specimen, distinguished the patients with versus without pathological complete response (pCR versus no-pCR respectively). We compared the mutational profiles, spectra and signatures, pathway mutations, copy number alterations (CNAs), HRD, and heterogeneity scores between pCR versus no-pCR patients.

RESULTS

The TMB, HRD, and mutational spectra were not different between the complete (N = 13) versus non-complete (N = 31) responders. The two most frequently mutated genes were TP53 and PIK3CA. They were more frequently mutated in the complete responders, but the difference was not significant. Only two genes, NLRP3 and SLC9B1, were significantly more frequently mutated in the complete responders (23% vs. 0%). By contrast, several biological pathways involved in protein translation, PI3K pathway, and signal transduction showed significantly higher mutation frequency in the patients with pCR. We observed a higher abundance of COSMIC signature 7 (due to ultraviolet light exposure) in tumors from complete responders. The comparison of CNAs of the 3808 genes included in the GISTIC regions between both patients' groups identified 234 genes as differentially altered. The CIN signatures were not differentially represented between the complete versus non-complete responders. Based on the H-index, the patients with heterogeneous tumors displayed a lower pCR rate (11%) than those with less heterogeneous tumors (35%).

CONCLUSIONS

This is the first study aiming at identifying correlations between the WES data of IBC samples and the achievement of pCR to NACT. Our results, obtained in this 44-sample series, suggest a few subtle genomic alterations associated with pathological response. Additional investigations are required in larger series.

Multiomic and clinical analysis of multiply recurrent meningiomas reveals risk factors, underlying biology, and insights into evolution

An important subset of meningiomas behaves aggressively and is characterized by multiple recurrences. We identify clinical, genetic, and epigenetic predictors of multiply recurrent meningiomas (MRMs) and evaluate the evolution of these meningiomas in patient-matched samples. On multivariable binomial logistic regression, MRMs were significantly associated with male sex (P = 0.012), subtotal resection (P = 0.001), higher number of meningiomas on presentation (P = 0.017), and histopathological sheeting (P = 0.002). Multiomic analysis of primary meningiomas revealed that MRMs have greater copy number losses (P = 0.0313) and increased DNA methylation (P = 0.0155). In meningioma cells with knockdown of EDNRB, a locus with greater promoter methylation and decreased gene expression in MRMs had increased proliferation (P < 0.0001). MRM recurrences were found to be similar to primaries but have a greater burden of copy number gains (P < 0.0001) and increased methylation (P = 0.0045). This clinical and multiomic investigation of MRMs harbors implications for the future development of biomarkers and therapeutic agents for these challenging tumors.

Parameter estimation from single patient, single time-point sequencing data of recurrent tumors

In this study, we develop consistent estimators for key parameters that govern the dynamics of tumor cell populations when subjected to pharmacological treatments. While these treatments often lead to an initial reduction in the abundance of drug-sensitive cells, a population of drug-resistant cells frequently emerges over time, resulting in cancer recurrence. Samples from recurrent tumors present as an invaluable data source that can offer crucial insights into the ability of cancer cells to adapt and withstand treatment interventions. To effectively utilize the data obtained from recurrent tumors, we derive several large number limit theorems, specifically focusing on the metrics that quantify the clonal diversity of cancer cell populations at the time of cancer recurrence. These theorems then serve as the foundation for constructing our estimators. A distinguishing feature of our approach is that our estimators only require a single time-point sequencing data from a single tumor, thereby enhancing the practicality of our approach and enabling the understanding of cancer recurrence at the individual level.

Molecular targets and strategies in the development of nucleic acid cancer vaccines: from shared to personalized antigens

Recent breakthroughs in cancer immunotherapies have emphasized the importance of harnessing the immune system for treating cancer. Vaccines, which have traditionally been used to promote protective immunity against pathogens, are now being explored as a method to target cancer neoantigens. Over the past few years, extensive preclinical research and more than a hundred clinical trials have been dedicated to investigating various approaches to neoantigen discovery and vaccine formulations, encouraging development of personalized medicine. Nucleic acids (DNA and mRNA) have become particularly promising platform for the development of these cancer immunotherapies. This shift towards nucleic acid-based personalized vaccines has been facilitated by advancements in molecular techniques for identifying neoantigens, antigen prediction methodologies, and the development of new vaccine platforms. Generating these personalized vaccines involves a comprehensive pipeline that includes sequencing of patient tumor samples, data analysis for antigen prediction, and tailored vaccine manufacturing. In this review, we will discuss the various shared and personalized antigens used for cancer vaccine development and introduce strategies for identifying neoantigens through the characterization of gene mutation, transcription, translation and post translational modifications associated with oncogenesis. In addition, we will focus on the most up-to-date nucleic acid vaccine platforms, discuss the limitations of cancer vaccines as well as provide potential solutions, and raise key clinical and technical considerations in vaccine development.

Aristolochic acid-related renal cell carcinoma exhibits a distinct tumor-immune microenvironment favoring response to immune checkpoint blockade

Immune checkpoint blockade (ICB) is currently the standard of care for metastatic renal cell carcinoma (RCC), but treatment responses remain unpredictable. Aristolochic acid (AA), a prevalent supplement additive in Taiwan, has been associated with RCC and induces signature mutations, although its effect on the tumor-immune microenvironment (TIME) is unclear. We aimed to investigate the immune profile of AA-positive RCCs and explore its potential role as a susceptible candidate for ICB. Tissue samples from 22 patients with clear cell RCC (ccRCC) were collected for whole-exome sequencing to determine the genetic features and AA mutational signature (the discovery cohort). The corresponding RNA was sent for NanoString PanCancer IO 360 gene expression analysis to explore the immunological features. The formalin-fixed, parafilm-embedded slides of ccRCCs were sent for multiplex immunohistochemistry/immunofluorescence stain using Vectra system to evaluate the TIME. Tissues from two patients with metastatic RCC demonstrating complete response to ICB were sent for studies to validate the findings (the index patients). The results showed that AA mutational signatures with high tumor mutational burden (TMB) were present in 31.81% of the tumors in the discovery cohort. Three distinct clusters were observed through NanoString analysis. Clusters 1 and 3 were composed mainly of AA-positive RCCs. Cluster 3 RCCs exhibited higher tumor inflammation signature scores and higher immune cell type scores. Vectra analysis revealed a higher percentage of CD15+ and BATF3+ cells in cluster 1, whereas the percentage of CD8+ cells was potentially higher in cluster 3. Strong AA mutational signatures were found in the tumors of two index patients, and both were grouped to cluster 3. In conclusion, AA may induce higher TMB and alter the immune microenvironment in RCCs, which makes the tumors more susceptible to ICB. © 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.

Exploring evolutionary trajectories in ovarian cancer patients by longitudinal analysis of ctDNA

OBJECTIVES

We analysed whether temporal heterogeneity of ctDNA encodes evolutionary patterns in ovarian cancer.

METHODS

Targeted sequencing of 275 cancer-associated genes was performed in a primary tumor biopsy and in ctDNA of six longitudinal plasma samples from 15 patients, using the Illumina platform.

RESULTS

While there was low overall concordance between the mutational spectrum of the primary tumor biopsies vs. ctDNA, TP53 variants were the most commonly shared somatic alterations. Up to three variant clusters were detected in each tumor biopsy, likely representing predominant clones of the primary tumor, most of them harbouring a TP53 variant. By tracing these clusters in ctDNA, we propose that liquid biopsy may allow to assess the contribution of ancestral clones of the tumor to relapsed abdominal masses, revealing two evolutionary patterns. In pattern#1, clusters detected in the primary tumor biopsy were likely relapse seeding clones, as they contributed a major share to ctDNA at relapse. In pattern#2, similar clusters were present in tumors and ctDNA; however, they were entirely cleared from liquid biopsy after chemotherapy and were undetectable at relapse. ctDNA private variants were present among both patterns, with some of them mirroring subclonal expansions after chemotherapy.

CONCLUSIONS

We demonstrate that tracing the temporal heterogeneity of ctDNA, even below exome scale resolution, deciphers evolutionary trajectories in ovarian cancer. Furthermore, we describe two evolutionary patterns that may help to identify relapse seeding clones for targeted therapy.

A Comparative Genomic Study of Conventional and Undifferentiated Melanoma

Undifferentiated melanoma, defined as melanoma that has lost all usual phenotypic and immunohistochemical characteristics of conventional melanoma, can pose significant diagnostic challenges. Molecular studies have advanced our understanding of undifferentiated melanoma by demonstrating that a subset of these tumors harbors known melanoma driver alterations in genes such as BRAF, NRAS, and NF1. However, there is a paucity of data describing genetic alterations that may distinguish undifferentiated melanoma from conventional melanoma. In this study, we directly compared the genomic profiles of undifferentiated melanoma to a cohort of conventional melanomas, including 14 undifferentiated melanoma cases (comprised of 2 primary cases, 2 cutaneous recurrences, and 10 metastases) and a cohort of 127 conventional melanomas including primary, recurrent, and metastatic cases. Targeted sequencing of 447 cancer-associated genes was performed, including identification of mutations and copy number alterations. NRAS was the most frequent melanoma driver in undifferentiated melanoma (8/14 cases, 57%), although notably, only 1 undifferentiated melanoma harbored an NRAS Q61R mutation. Compared with the conventional melanoma cohort, undifferentiated melanoma demonstrated statistically significant enrichment of pathogenic activating RAC1 mutations (6/14 total cases, 43%), including P29S (4/6 cases), P29L (1/6 cases), and D11E (1/6 cases). In addition to providing insight into the molecular pathogenesis of undifferentiated melanoma, these findings also suggest that RAS Q61R immunohistochemistry may have limited utility for its diagnosis. The presence of recurrent RAC1 mutations in undifferentiated melanoma is also notable as these alterations may contribute to mitogen-activated protein kinase pathway-targeted therapy resistance. Furthermore, the RAC1 alterations identified in this cohort have been shown to drive a melanocytic to mesenchymal switch in melanocytes, offering a possible explanation for the undifferentiated phenotype of these melanomas.

Reconstructing tumor clonal heterogeneity and evolutionary relationships based on tumor DNA sequencing data

The heterogeneity of tumor clones drives the selection and evolution of distinct tumor cell populations, resulting in an intricate and dynamic tumor evolution process. While tumor bulk DNA sequencing helps elucidate intratumor heterogeneity, challenges such as the misidentification of mutation multiplicity due to copy number variations and uncertainties in the reconstruction process hinder the accurate inference of tumor evolution. In this study, we introduce a novel approach, REconstructing Tumor Clonal Heterogeneity and Evolutionary Relationships (RETCHER), which characterizes more realistic cancer cell fractions by accurately identifying mutation multiplicity while considering uncertainty during the reconstruction process and the credibility and reasonableness of subclone clustering. This method comprehensively and accurately infers multiple forms of tumor clonal heterogeneity and phylogenetic relationships. RETCHER outperforms existing methods on simulated data and infers clearer subclone structures and evolutionary relationships in real multisample sequencing data from five tumor types. By precisely analysing the complex clonal heterogeneity within tumors, RETCHER provides a new approach to tumor evolution research and offers scientific evidence for developing precise and personalized treatment strategies. This approach is expected to play a significant role in tumor evolution research, clinical diagnosis, and treatment. RETCHER is available for free at https://github.com/zlsys3/RETCHER.

Comparative analysis of the mutational landscape and evolutionary patterns of pancreatic ductal adenocarcinoma metastases in the liver or peritoneum

Background

Pancreatic ductal adenocarcinoma (PDAC) often presents with liver or peritoneal metastases at diagnosis. Despite similar treatment approaches, patient outcomes vary between these metastatic sites. To improve targeted therapies for metastatic PDAC, a comprehensive analysis of the genetic profiles and evolutionary patterns at these distinct metastatic locations is essential.

Methods

We performed whole exome sequencing on 44 tissue samples from 27 PDAC patients, including primary tumours and matched liver or peritoneal metastases. We analysed somatic mutation profiles, signatures, and affected pathways for each group, and examined clonal evolution using subclonal architectures and phylogenetic trees.

Results

KRAS mutations remained the predominant driver alteration, with a prevalence of 89 % across all tumours. Notably, we observed site-specific differences in mutation frequencies, with KRAS alterations detected in 77.8 % (7/9) of peritoneal metastases and 87.5 % (7/8) of liver metastases. TP53 mutations exhibited a similar pattern, occurring in 55.6 % (5/9) of peritoneal and 37.5 % (3/8) of liver metastases. Intriguingly, we identified site-specific alterations in DNA repair pathway genes, including ATM and BRCA1, with distinct mutational profiles in liver versus peritoneal metastases. Furthermore, liver metastases demonstrated a significantly higher tumor mutational burden (TMB) compared to peritoneal metastases (median [IQR]: 2.14 [1.77-2.71] vs. 1.29 [1.21-1.69] mutations/Mb; P = 0.048).

Conclusions

In conclusion, metastasis of pancreatic cancer may be influenced by variables other than KRAS mutations, such as TP53. PDAC peritoneal and liver metastases may differ in potential therapeutic biomarkers. Further inquiry is needed on the biological mechanisms underlying metastasis and the treatment of diverse metastases.

CLEMENT: genomic decomposition and reconstruction of non-tumor subclones

Genome-level clonal decomposition of a single specimen has been widely studied; however, it is mostly limited to cancer research. In this study, we developed a new algorithm CLEMENT, which conducts accurate decomposition and reconstruction of multiple subclones in genome sequencing of non-tumor (normal) samples. CLEMENT employs the Expectation-Maximization (EM) algorithm with optimization strategies specific to non-tumor subclones, including false variant call identification, non-disparate clone fuzzy clustering, and clonal allele fraction confinement. In the simulation and in vitro cell line mixture data, CLEMENT outperformed current cancer decomposition algorithms in estimating the number of clones (root-mean-square-error = 0.58-0.78 versus 1.43-3.34) and in the variant-clone membership agreement (∼85.5% versus 70.1-76.7%). Additional testing on human multi-clonal normal tissue sequencing confirmed the accurate identification of subclones that originated from different cell types. Clone-level analysis, including mutational burden and signatures, provided a new understanding of normal-tissue composition. We expect that CLEMENT will serve as a crucial tool in the currently emerging field of non-tumor genome analysis.

Somatic mutation phasing and haplotype extension using linked-reads in multiple myeloma

Somatic mutation phasing informs our understanding of cancer-related events, like driver mutations. We generated linked-read whole genome sequencing data for 23 samples across disease stages from 14 multiple myeloma (MM) patients and systematically assigned somatic mutations to haplotypes using linked-reads. Here, we report the reconstructed cancer haplotypes and phase blocks from several MM samples and show how phase block length can be extended by integrating samples from the same individual. We also uncover phasing information in genes frequently mutated in MM, including DIS3, HIST1H1E, KRAS, NRAS, and TP53, phasing 79.4% of 20,705 high-confidence somatic mutations. In some cases, this enabled us to interpret clonal evolution models at higher resolution using pairs of phased somatic mutations. For example, our analysis of one patient suggested that two NRAS hotspot mutations occurred on the same haplotype but were independent events in different subclones. Given sufficient tumor purity and data quality, our framework illustrates how haplotype-aware analysis of somatic mutations in cancer can be beneficial for some cancer cases.

Differential impact of intratumor heterogeneity (ITH) on survival outcomes in early-stage lung squamous and adenocarcinoma based on tumor mutational burden (TMB)

Background

Molecular biomarkers are reshaping patient stratification and treatment decisions, yet their precise use and best implementation remain uncertain. Intratumor heterogeneity (ITH), an area of increasing research interest with prognostic value across various conditions, lacks defined clinical relevance in certain non-small cell lung cancer (NSCLC) subtypes. Exploring the relationship between ITH and tumor mutational burden (TMB) is crucial, as their interplay might reveal distinct patient subgroups. This study evaluates how the ITH-TMB dynamic affects prognosis across the two main histological subtypes of NSCLC, squamous cell and adenocarcinoma, with a specific focus on early-stage cases to address their highly unmet clinical needs.

Methods

We stratify a cohort of 741 early-stage NSCLC patients from The Cancer Genome Atlas (TCGA) based on ITH and TMB and evaluate differences in clinical outcomes. Additionally, we compare driver mutations and the tumor microenvironment (TME) between high and low ITH groups.

Results

In lung squamous cell carcinoma (LUSC), high ITH predicts an extended progression-free survival (PFS) (median: 21 vs. 14 months, P=0.01), while in lung adenocarcinoma (LUAD), high ITH predicts a reduced PFS (median: 15 vs. 20 months, P=0.04). This relationship is driven by the low TMB subset of patients. Additionally, we found that CD8 T cells were enriched in better-performing subgroups, regardless of histologic subtype or ITH status.

Conclusions

There are significant differences in clinical outcomes, driver mutations, and the TME between high and low ITH groups among early-stage NSCLC patients. These differences may have treatment implications, necessitating further validation in other NSCLC datasets.

A Novel Affordable and Reliable Framework for Accurate Detection and Comprehensive Analysis of Somatic Mutations in Cancer

Accurate detection and analysis of somatic variants in cancer involve multiple third-party tools with complex dependencies and configurations, leading to laborious, error-prone, and time-consuming data conversions. This approach lacks accuracy, reproducibility, and portability, limiting clinical application. Musta was developed to address these issues as an end-to-end pipeline for detecting, classifying, and interpreting cancer mutations. Musta is based on a Python command-line tool designed to manage tumor-normal samples for precise somatic mutation analysis. The core is a Snakemake-based workflow that covers all key cancer genomics steps, including variant calling, mutational signature deconvolution, variant annotation, driver gene detection, pathway analysis, and tumor heterogeneity estimation. Musta is easy to install on any system via Docker, with a Makefile handling installation, configuration, and execution, allowing for full or partial pipeline runs. Musta has been validated at the CRS4-NGS Core facility and tested on large datasets from The Cancer Genome Atlas and the Beijing Institute of Genomics. Musta has proven robust and flexible for somatic variant analysis in cancer. It is user-friendly, requiring no specialized programming skills, and enables data processing with a single command line. Its reproducibility ensures consistent results across users following the same protocol.

Evolutionary history of adenomas to colorectal cancer in FAP families

Objective

Familial adenomatous polyposis (FAP) is a genetic syndrome characterized by multiple polyps at various evolutionary stages, which, if left untreated, inevitably progress to colorectal cancer (CRC). In this study, we present a comprehensive analysis of the evolutionary history of FAP-CRC from precancerous adenoma to carcinoma.

Design

Tissues were collected from gastrointestinal endoscopy or surgical resection. Exome sequencing was performed on multiple regions of adenocarcinoma (n = 8), villous adenoma (n = 10), tubular adenoma (n = 9) and blood samples were obtained from 9 patients belonging to 7 Chinese FAP families. Phylogenetic trees were reconstructed, and evolutionary analysis was conducted to reveal the temporal sequence of events leading to CRC.

Results

Inherited germline mutation sites in APC gene were identified in FAP01 (p.S1281*, COSM19212), FAP03 (p.S384Tfs*19), FAP04 (p.E1538*, COSM6041693), FAP05 (p.Q1062*, COSM3696862), and FAP07-FAP09 (p.V677Sfs*3). Notably, p.V677Sfs*3 mutation was recognized as a novel germline mutation in APC, supported by evidence of genotype-phenotype correlation in pedigree analysis. Adenomas exhibited lower mutational rates than FAP-CRC and displayed recurrent alterations in well-known chromosomal instability (CIN) genes (APC, RAS, SMAD4 and TP53) and DNA damage repair genes (SUZ12, KMT2C, BCLAF1, RUNX1, and ARID1B), suggesting the presence of genomic instability. Furthermore, a progressive increase in the HRD score (a measure of "genomic scars") was observed from tubular adenomas to villous adenomas and ultimately to carcinomas. TP53 emerged as the primary driver gene for adenoma-carcinoma transition, with driver mutations consistently appearing simultaneously rather than sequentially acquired from adenomas to carcinomas. Clonal evolution demonstrated that liver metastases can originate from the same cancer-primed cell present in a primary cancerous lesion.

Conclusion

We identified a novel pathogenic variant in APC, namely, p.V677Sfs*3. The process of carcinogenesis in FAP-CRC supports the classical cancerization model, where an initial APC mutation leads to the activation of the WNT signaling pathway and CIN. Subsequently, additional mutations occur in other putative CIN genes (e.g., DNA repair, chromatin remodeling), ultimately leading to the development of microsatellite stable (MSS) tumors. Our study provides a comprehensive understanding of the genomic landscapes that underlie the transition from adenoma to carcinoma.

Integrating Multisector Molecular Characterization into Personalized Peptide Vaccine Design for Patients with Newly Diagnosed Glioblastoma

PURPOSE

Outcomes for patients with glioblastoma (GBM) remain poor despite multimodality treatment with surgery, radiation, and chemotherapy. There are few immunotherapy options due to the lack of tumor immunogenicity. Several clinical trials have reported promising results with cancer vaccines. To date, studies have used data from a single tumor site to identify targetable antigens, but this approach limits the antigen pool and is antithetical to the heterogeneity of GBM. We have implemented multisector sequencing to increase the pool of neoantigens across the GBM genomic landscape that can be incorporated into personalized peptide vaccines called NeoVax.

PATIENTS AND METHODS

In this study, we report the findings of four patients enrolled onto the NeoVax clinical trial (NCT0342209).

RESULTS

Immune reactivity to NeoVax neoantigens was assessed in peripheral blood mononuclear cells pre- and post-NeoVax for patients 1 to 3 using IFNγ-ELISPOT assay. A statistically significant increase in IFNγ producing T cells at the post-NeoVax time point for several neoantigens was observed. Furthermore, a post-NeoVax tumor biopsy was obtained from patient 3 and, upon evaluation, revealed evidence of infiltrating, clonally expanded T cells.

CONCLUSIONS

Collectively, our findings suggest that NeoVax stimulated the expansion of neoantigen-specific effector T cells and provide encouraging results to aid in the development of future neoantigen vaccine-based clinical trials in patients with GBM. Herein, we demonstrate the feasibility of incorporating multisector sampling in cancer vaccine design and provide information on the clinical applicability of clonality, distribution, and immunogenicity of the neoantigen landscape in patients with GBM.

Resolving tumor evolution: a phylogenetic approach

The evolutionary dynamics of cancer, characterized by its profound heterogeneity, demand sophisticated tools for a holistic understanding. This review delves into tumor phylogenetics, an essential approach bridging evolutionary biology with oncology, offering unparalleled insights into cancer's evolutionary trajectory. We provide an overview of the workflow, encompassing study design, data acquisition, and phylogeny reconstruction. Notably, the integration of diverse data sets emerges as a transformative step, enhancing the depth and breadth of evolutionary insights. With this integrated perspective, tumor phylogenetics stands poised to redefine our understanding of cancer evolution and influence therapeutic strategies.

Enforcing Temporal Consistency in Migration History Inference

In addition to undergoing evolution, members of biological populations may also migrate between locations. Examples include the spread of tumor cells from the primary tumor to distant metastases or the spread of pathogens from one host to another. One may represent migration histories by assigning a location label to each vertex of a given phylogenetic tree such that an edge connecting vertices with distinct locations represents a migration. Some biological populations undergo comigration, a phenomenon where multiple taxa from distinct lineages simultaneously comigrate from one location to another. In this work, we show that a previous problem statement for inferring migration histories that are parsimonious in terms of migrations and comigrations may lead to temporally inconsistent solutions. To remedy this deficiency, we introduce precise definitions of temporal consistency of comigrations in a phylogenetic tree, leading to three successive problems. First, we formulate the temporally consistent comigration problem to check if a set of comigrations is temporally consistent and provide a linear time algorithm for solving this problem. Second, we formulate the parsimonious consistent comigrations (PCC) problem, which aims to find comigrations given a location labeling of a phylogenetic tree. We show that PCC is NP-hard. Third, we formulate the parsimonious consistent comigration history (PCCH) problem, which infers the migration history given a phylogenetic tree and locations of its extant vertices only. We show that PCCH is NP-hard as well. On the positive side, we propose integer linear programming models to solve the PCC and PCCH problems. We demonstrate our algorithms on simulated and real data.

Establishment of Human Pituitary Neuroendocrine Tumor Derived Organoid and Its Pilot Application for Drug Screening

CONTEXT

Precision medicine for pituitary neuroendocrine tumors (PitNETs) is limited by the lack of reliable research models.

OBJECTIVE

To generate patient-derived organoids (PDOs), which could serve as a platform for personalized drug screening for PitNET patients.

DESIGN

From July 2019 to May 2022, a total of 32 human PitNET specimens were collected for the establishment of organoids with an optimized culture protocol.

SETTING

This study was conducted at Sun Yat-Sen University Cancer Center.

PATIENTS

PitNET patients who were pathologically confirmed were enrolled in this study.

INTERVENTIONS

Histological staining and whole-exome sequencing were utilized to confirm the pathologic and genomic features of PDOs. A drug response assay on PDOs was also performed.

MAIN OUTCOME MEASURES

PDOs retained key genetic and morphological features of their parental tumors.

RESULTS

PDOs were successfully established from various types of PitNET samples with an overall success rate of 87.5%. Clinical nonfunctioning PitNETs-derived organoids (22/23, 95.7%) showed a higher likelihood of successful generation compared to those from functioning PitNETs (6/9, 66.7%). Preservation of cellular structure, subtype-specific neuroendocrine profiles, mutational features, and tumor microenvironment heterogeneity from parental tumors was observed. A distinctive response profile in drug tests was observed among the organoids from patients with different subtypes of PitNETs. With the validation of key characteristics from parental tumors in histological, genomic, and microenvironment heterogeneity consistency assays, we demonstrated the predictive value of the PDOs in testing individual drugs.

CONCLUSION

The established PDOs, retaining typical features of parental tumors, indicate a translational significance in innovating personalized treatment for refractory PitNETs.

Mutational landscape of inflammatory breast cancer

BACKGROUND

Inflammatory breast cancer (IBC) is the most pro-metastatic form of BC. Better understanding of its enigmatic pathophysiology is crucial. We report here the largest whole-exome sequencing (WES) study of clinical IBC samples.

METHODS

We retrospectively applied WES to 54 untreated IBC primary tumor samples and matched normal DNA. The comparator samples were 102 stage-matched non-IBC samples from TCGA. We compared the somatic mutational profiles, spectra and signatures, copy number alterations (CNAs), HRD and heterogeneity scores, and frequencies of actionable genomic alterations (AGAs) between IBCs and non-IBCs. The comparisons were adjusted for the molecular subtypes.

RESULTS

The number of somatic mutations, TMB, and mutational spectra were not different between IBCs and non-IBCs, and no gene was differentially mutated or showed differential frequency of CNAs. Among the COSMIC signatures, only the age-related signature was more frequent in non-IBCs than in IBCs. We also identified in IBCs two new mutational signatures not associated with any environmental exposure, one of them having been previously related to HIF pathway activation. Overall, the HRD score was not different between both groups, but was higher in TN IBCs than TN non-IBCs. IBCs were less frequently classified as heterogeneous according to heterogeneity H-index than non-IBCs (21% vs 33%), and clonal mutations were more frequent and subclonal mutations less frequent in IBCs. More than 50% of patients with IBC harbored at least one high-level of evidence (LOE) AGA (OncoKB LOE 1-2, ESCAT LOE I-II), similarly to patients with non-IBC.

CONCLUSIONS

We provide the largest mutational landscape of IBC. Only a few subtle differences were identified with non-IBCs. The most clinically relevant one was the higher HRD score in TN IBCs than in TN non-IBCs, whereas the most intriguing one was the smaller intratumor heterogeneity of IBCs.

Treatment with the apoptosis inhibitor Asunercept reduces clone sizes in patients with lower risk Myelodysplastic Neoplasms

In low-risk Myelodysplastic Neoplasms (MDS), increased activity of apoptosis-promoting factors such as tumor necrosis factor (TNFα) and pro-apoptotic Fas ligand (CD95L) have been described as possible pathomechanisms leading to impaired erythropoiesis. Asunercept (APG101) is a novel therapeutic fusion protein blocking CD95, which has previously shown partial efficacy in reducing transfusion requirement in a clinical phase I trial for low-risk MDS patients (NCT01736436; 2012-11-26). In the current study we aimed to evaluate the effect of Asunercept therapy on the clonal bone marrow composition to identify potential biomarkers to predict response. Bone marrow samples of n = 12 low-risk MDS patients from the above referenced clinical trial were analyzed by serial deep whole exome sequencing in a total of n = 58 time points. We could distinguish a mean of 3.5 molecularly defined subclones per patient (range 2-6). We observed a molecular response defined as reductions of dominant clone sizes by a variant allele frequency (VAF) decrease of at least 10% (mean 20%, range: 10.5-39.2%) in dependency of Asunercept treatment in 9 of 12 (75%) patients. Most of this decline in clonal populations was observed after completion of 12 weeks treatment. Particularly early and pronounced reductions of clone sizes were found in subclones driven by mutations in genes involved in regulation of methylation (n = 1 DNMT3A, n = 1 IDH2, n = 1 TET2). Our results suggest that APG101 could be efficacious in reducing clone sizes of mutated hematopoietic cells in the bone marrow of Myelodysplastic Neoplasms, which warrants further investigation.

Characterization of tumor evolution by functional clonality and phylogenetics in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a molecularly heterogeneous solid malignancy, and its fitness may be shaped by how its tumor cells evolve. However, ability to monitor tumor cell evolution is hampered by the presence of numerous passenger mutations that do not provide any biological consequences. Here we develop a strategy to determine the tumor clonality of three independent HCC cohorts of 524 patients with diverse etiologies and race/ethnicity by utilizing somatic mutations in cancer driver genes. We identify two main types of tumor evolution, i.e., linear, and non-linear models where non-linear type could be further divided into classes, which we call shallow branching and deep branching. We find that linear evolving HCC is less aggressive than other types. GTF2IRD2B mutations are enriched in HCC with linear evolution, while TP53 mutations are the most frequent genetic alterations in HCC with non-linear models. Furthermore, we observe significant B cell enrichment in linear trees compared to non-linear trees suggesting the need for further research to uncover potential variations in immune cell types within genomically determined phylogeny types. These results hint at the possibility that tumor cells and their microenvironment may collectively influence the tumor evolution process.

Genomic characterization of peritoneal lavage cytology-positive gastric cancer

Objective

Positive peritoneal lavege cytology (CY1) gastric cancer is featured by dismal prognosis, with high risks of peritoneal metastasis. However, there is a lack of evidence on pathogenic mechanism and signature of CY1 and there is a continuous debate on CY1 therapy. Therefore, exploring the mechanism of CY1 is crucial for treatment strategies and targets for CY1 gastric cancer.

Methods

In order to figure out specific driver genes and marker genes of CY1 gastric cancer, and ultimately offer clues for potential marker and risk assessment of CY1, 17 cytology-positive gastric cancer patients and 31 matched cytology-negative gastric cancer patients were enrolled in this study. The enrollment criteria were based on the results of diagnostic laparoscopy staging and cytology inspection of exfoliated cells. Whole exome sequencing was then performed on tumor samples to evaluate genomic characterization of cytology-positive gastric cancer.

Results

Least absolute shrinkage and selection operator (LASSO) algorithm identified 43 cytology-positive marker genes, while MutSigCV identified 42 cytology-positive specific driver genes. CD3G and CDKL2 were both driver and marker genes of CY1. Regarding mutational signatures, driver gene mutation and tumor subclone architecture, no significant differences were observed between CY1 and negative peritoneal lavege cytology (CY0).

Conclusions

There might not be distinct differences between CY1 and CY0, and CY1 might represent the progression of CY0 gastric cancer rather than constituting an independent subtype. This genomic analysis will thus provide key molecular insights into CY1, which may have a direct effect on treatment recommendations for CY1 and CY0 patients, and provides opportunities for genome-guided clinical trials and drug development.

Sequential antigen loss and branching evolution in lymphoma after CD19- and CD20-targeted T-cell-redirecting therapy

ABSTRACT

CD19 chimeric antigen receptor (CAR) T cells and CD20 targeting T-cell-engaging bispecific antibodies (bispecs) have been approved in B-cell non-Hodgkin lymphoma lately, heralding a new clinical setting in which patients are treated with both approaches, sequentially. The aim of our study was to investigate the selective pressure of CD19- and CD20-directed therapy on the clonal architecture in lymphoma. Using a broad analytical pipeline on 28 longitudinally collected specimen from 7 patients, we identified truncating mutations in the gene encoding CD20 conferring antigen loss in 80% of patients relapsing from CD20 bispecs. Pronounced T-cell exhaustion was identified in cases with progressive disease and retained CD20 expression. We also confirmed CD19 loss after CAR T-cell therapy and reported the case of sequential CD19 and CD20 loss. We observed branching evolution with re-emergence of CD20+ subclones at later time points and spatial heterogeneity for CD20 expression in response to targeted therapy. Our results highlight immunotherapy as not only an evolutionary bottleneck selecting for antigen loss variants but also complex evolutionary pathways underlying disease progression from these novel therapies.

Natural History of Germline BRCA1 Mutated and BRCA Wild-type Triple-negative Breast Cancer

We report a deep next-generation sequencing analysis of 13 sequentially obtained tumor samples, eight sequentially obtained circulating tumor DNA (ctDNA) samples and three germline DNA samples over the life history of 3 patients with triple-negative breast cancer (TNBC), 2 of whom had germline pathogenic BRCA1 mutation, to unravel tumor evolution. Tumor tissue from all timepoints and germline DNA was subjected to whole-exome sequencing (WES), custom amplicon deep sequencing (30,000X) of a WES-derived somatic mutation panel, and SNP arrays for copy-number variation (CNV), while whole transcriptome sequencing (RNA-seq) was performed only on somatic tumor.There was enrichment of homologous recombination deficiency signature in all tumors and widespread CNV, which remained largely stable over time. Somatic tumor mutation numbers varied between patients and within each patient (range: 70-216, one outlier). There was minimal mutational overlap between patients with TP53 being the sole commonly mutated gene, but there was substantial overlap in sequential samples in each patient. Each patient's tumor contained a founding ("stem") clone at diagnosis, which persisted over time, from which all other clones ("subclone") were derived ("branching evolution"), which contained mutations in well-characterized cancer-related genes like PDGFRB, ARID2, TP53 (Patient_02), TP53, BRAF, BRIP1, CSF3R (Patient_04), and TP53, APC, EZH2 (Patient_07). Including stem and subclones, tumors from all patients were polyclonal at diagnosis and during disease progression. ctDNA recapitulated most tissue-derived stem clonal and subclonal mutations while detecting some additional subclonal mutations. RNA-seq revealed a stable basal-like pattern, with most highly expressed variants belonging to stem clone.

SIGNIFICANCE

In germline BRCA1 mutated and BRCA wild-type patients, TNBC shows a branching evolutionary pattern of mutations with a single founding clone, are polyclonal throughout their disease course, and have widespread copy-number aberrations. This evolutionary pattern may be associated with treatment resistance or sensitivity and could be therapeutically exploited.

A Bayesian framework to study tumor subclone-specific expression by combining bulk DNA and single-cell RNA sequencing data

Genetic and gene expression heterogeneity is an essential hallmark of many tumors, allowing the cancer to evolve and to develop resistance to treatment. Currently, the most commonly used data types for studying such heterogeneity are bulk tumor/normal whole-genome or whole-exome sequencing (WGS, WES); and single-cell RNA sequencing (scRNA-seq), respectively. However, tools are currently lacking to link genomic tumor subclonality with transcriptomic heterogeneity by integrating genomic and single-cell transcriptomic data collected from the same tumor. To address this gap, we developed scBayes, a Bayesian probabilistic framework that uses tumor subclonal structure inferred from bulk DNA sequencing data to determine the subclonal identity of cells from single-cell gene expression (scRNA-seq) measurements. Grouping together cells representing the same genetically defined tumor subclones allows comparison of gene expression across different subclones, or investigation of gene expression changes within the same subclone across time (i.e., progression, treatment response, or relapse) or space (i.e., at multiple metastatic sites and organs). We used simulated data sets, in silico synthetic data sets, as well as biological data sets generated from cancer samples to extensively characterize and validate the performance of our method, as well as to show improvements over existing methods. We show the validity and utility of our approach by applying it to published data sets and recapitulating the findings, as well as arriving at novel insights into cancer subclonal expression behavior in our own data sets. We further show that our method is applicable to a wide range of single-cell sequencing technologies including single-cell DNA sequencing as well as Smart-seq and 10x Genomics scRNA-seq protocols.

Computational validation of clonal and subclonal copy number alterations from bulk tumor sequencing using CNAqc

Copy number alterations (CNAs) are among the most important genetic events in cancer, but their detection from sequencing data is challenging because of unknown sample purity, tumor ploidy, and general intra-tumor heterogeneity. Here, we present CNAqc, an evolution-inspired method to perform the computational validation of clonal and subclonal CNAs detected from bulk DNA sequencing. CNAqc is validated using single-cell data and simulations, is applied to over 4000 TCGA and PCAWG samples, and is incorporated into the validation process for the clinically accredited bioinformatics pipeline at Genomics England. CNAqc is designed to support automated quality control procedures for tumor somatic data validation.

Computational immunogenomic approaches to predict response to cancer immunotherapies

Cancer immunogenomics is an emerging field that bridges genomics and immunology. The establishment of large-scale genomic collaborative efforts along with the development of new single-cell transcriptomic techniques and multi-omics approaches have enabled characterization of the mutational and transcriptional profiles of many cancer types and helped to identify clinically actionable alterations as well as predictive and prognostic biomarkers. Researchers have developed computational approaches and machine learning algorithms to accurately obtain clinically useful information from genomic and transcriptomic sequencing data from bulk tissue or single cells and explore tumours and their microenvironment. The rapid growth in sequencing and computational approaches has resulted in the unmet need to understand their true potential and limitations in enabling improvements in the management of patients with cancer who are receiving immunotherapies. In this Review, we describe the computational approaches currently available to analyse bulk tissue and single-cell sequencing data from cancer, stromal and immune cells, as well as how best to select the most appropriate tool to address various clinical questions and, ultimately, improve patient outcomes.

Temporal genomic analysis of melanoma rejection identifies regulators of tumor immune evasion

Decreased intra-tumor heterogeneity (ITH) correlates with increased patient survival and immunotherapy response. However, even highly homogenous tumors may display variability in their aggressiveness, and how immunologic-factors impinge on their aggressiveness remains understudied. Here we studied the mechanisms responsible for the immune-escape of murine tumors with low ITH. We compared the temporal growth of homogeneous, genetically-similar single-cell clones that are rejected vs. those that are not-rejected after transplantation in-vivo using single-cell RNA sequencing and immunophenotyping. Non-rejected clones showed high infiltration of tumor-associated-macrophages (TAMs), lower T-cell infiltration, and increased T-cell exhaustion compared to rejected clones. Comparative analysis of rejection-associated gene expression programs, combined with in-vivo CRISPR knockout screens of candidate mediators, identified Mif (macrophage migration inhibitory factor) as a regulator of immune rejection. Mif knockout led to smaller tumors and reversed non-rejection-associated immune composition, particularly, leading to the reduction of immunosuppressive macrophage infiltration. Finally, we validated these results in melanoma patient data.

A novel subclonal rearrangement of the STRN3::PDGFRB gene in de novo acute myeloid leukemia with NPM1 mutation and its leukemogenic effects

Chromosome translocations in the 5q31-33 region are associated with a range of hematologic malignancies, some of which involve the platelet-derived growth factor receptor beta (PDGFRB) gene. We report a case of acute myeloid leukemia (AML) with a mutation in the NPM1 gene (NPM1-mut AML) and a subclonal gene rearrangement involving the PDGFRB gene. We identified a novel fusion gene, STRN3::PDGFRB, resulting from t(5;14) (q32;q12) chromosomal rearrangement. Sequential FISH confirmed that ~15% of leukemic cells carried the PDGFRB gene rearrangement, which suggests that STRN3::PDGFRB is a previously unreported fusion gene in a subclone. Reverse transcription PCR (RT-PCR) and Sanger sequencing confirmed that the fusion gene consisted of STRN3 exon 7 fused to PDGFRB exon 11, resulting in a chimeric protein containing the coiled-coil domain of striatin-3 and the transmembrane and intracellular tyrosine kinase domains of the PDGFRB. The new protein exhibited distinct cytoplasmic localization and had leukemogenic effects, as demonstrated by its ability to transform Ba/F3 cells to growth factor independence and cause a fatal myelodysplastic/myeloproliferative neoplasm (MDS/MPN)-like disease in mice, which then transformant to T-cell lymphoblastic lymphoma in secondary recipients. Ba/F3 cells expressing STRN3::PDGFRB or ETV6::PDGFRB were sensitive to tyrosine kinase inhibitors (TKIs) and selinexor, but in vitro experiments showed that the combination of imatinib and selinexor had a marked synergistic effect, although only the imatinib alone group could prolong the survival of T-cell blast transformation recipient mice. Our findings demonstrate the leukemogenic effects of the novel fusion gene and provide insights into the clone evolution of AML, which can be influenced by therapy selection. Furthermore, our results provide insight into the potential therapeutic options for patients with this type of mutation, as well as the need for careful consideration of treatment selection to prevent undesirable side effects.

Clinical-genomic determinants of immune checkpoint blockade response in head and neck squamous cell carcinoma

BACKGROUNDRecurrent and/or metastatic (R/M) head and neck squamous cell carcinoma (HNSCC) is generally an incurable disease, with patients experiencing median survival of under 10 months and significant morbidity. While immune checkpoint blockade (ICB) drugs are effective in approximately 20% of patients, the remaining experience limited clinical benefit and are exposed to potential adverse effects and financial costs. Clinically approved biomarkers, such as tumor mutational burden (TMB), have a modest predictive value in HNSCC.METHODSWe analyzed clinical and genomic features, generated using whole-exome sequencing, in 133 ICB-treated patients with R/M HNSCC, of whom 69 had virus-associated and 64 had non-virus-associated tumors.RESULTSHierarchical clustering of genomic data revealed 6 molecular subtypes characterized by a wide range of objective response rates and survival after ICB therapy. The prognostic importance of these 6 subtypes was validated in an external cohort. A random forest-based predictive model, using several clinical and genomic features, predicted progression-free survival (PFS), overall survival (OS), and response with greater accuracy than did a model based on TMB alone. Recursive partitioning analysis identified 3 features (systemic inflammatory response index, TMB, and smoking signature) that classified patients into risk groups with accurate discrimination of PFS and OS.CONCLUSIONThese findings shed light on the immunogenomic characteristics of HNSCC tumors that drive differential responses to ICB and identify a clinical-genomic classifier that outperformed the current clinically approved biomarker of TMB. This validated predictive tool may help with clinical risk stratification in patients with R/M HNSCC for whom ICB is being considered.FUNDINGFundación Alfonso Martín Escudero, NIH R01 DE027738, US Department of Defense CA210784, The Geoffrey Beene Cancer Research Center, The MSKCC Population Science Research Program, the Jayme Flowers Fund, the Sebastian Nativo Fund, and the NIH/NCI Cancer Center Support Grant P30 CA008748.

Advances in novel strategies for isolation, characterization, and analysis of CTCs and ctDNA

Over the past decade, the detection and analysis of liquid biopsy biomarkers such as circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA) have advanced significantly. They have received recognition for their clinical usefulness in detecting cancer at an early stage, monitoring disease, and evaluating treatment response. The emergence of liquid biopsy has been a helpful development, as it offers a minimally invasive, rapid, real-time monitoring, and possible alternative to traditional tissue biopsies. In resource-limited settings, the ideal platform for liquid biopsy should not only extract more CTCs or ctDNA from a minimal sample volume but also accurately represent the molecular heterogeneity of the patient's disease. This review covers novel strategies and advancements in CTC and ctDNA-based liquid biopsy platforms, including microfluidic applications and comprehensive analysis of molecular complexity. We discuss these systems' operational principles and performance efficiencies, as well as future opportunities and challenges for their implementation in clinical settings. In addition, we emphasize the importance of integrated platforms that incorporate machine learning and artificial intelligence in accurate liquid biopsy detection systems, which can greatly improve cancer management and enable precision diagnostics.

Proteomic profiling of gastric cancer with peritoneal metastasis identifies a protein signature associated with immune microenvironment and patient outcome

BACKGROUND

Peritoneal metastasis (PM) frequently occurs in patients with gastric cancer (GC) and is a major cause of mortality. Risk stratification for PM can optimize decision making in GC treatment.

METHODS

A total of 25 GC patients (13 with synchronous, 6 with metachronous PM and 6 PM-free) were included in this study. Quantitative proteomics by high-depth tandem mass tags labeling and whole-exome sequencing were conducted in primary GC and PM samples. Proteomic signature and prognostic model were established by machine learning algorithms in PM and PM-free GC, then validated in two external cohorts. Tumor-infiltrating immune cells in GC were analyzed by CIBERSORT.

RESULTS

Heterogeneity between paired primary and PM samples was observed at both genomic and proteomic levels. Compared to primary GC, proteome of PM samples was enriched in RNA binding and extracellular exosomes. 641 differently expressed proteins (DEPs) between primary GC of PM group and PM-free group were screened, which were enriched in extracellular exosome and cell adhesion pathways. Subsequently, a ten-protein signature was derived based on DEPs by machine learning. This signature was significantly associated with patient prognosis in internal cohort and two external proteomic datasets of diffuse and mixed type GC. Tumor-infiltrating immune cell analysis showed that the signature was associated with immune microenvironment of GC.

CONCLUSIONS

We characterized proteomic features that were informative for PM progression of GC. A protein signature associated with immune microenvironment and patient outcome was derived, and it could guide risk stratification and individualized treatment.

Single-cell methylation sequencing data reveal succinct metastatic migration histories and tumor progression models

Recent studies exploring the impact of methylation in tumor evolution suggest that although the methylation status of many of the CpG sites are preserved across distinct lineages, others are altered as the cancer progresses. Because changes in methylation status of a CpG site may be retained in mitosis, they could be used to infer the progression history of a tumor via single-cell lineage tree reconstruction. In this work, we introduce the first principled distance-based computational method, Sgootr, for inferring a tumor's single-cell methylation lineage tree and for jointly identifying lineage-informative CpG sites that harbor changes in methylation status that are retained along the lineage. We apply Sgootr on single-cell bisulfite-treated whole-genome sequencing data of multiregionally sampled tumor cells from nine metastatic colorectal cancer patients, as well as multiregionally sampled single-cell reduced-representation bisulfite sequencing data from a glioblastoma patient. We show that the tumor lineages constructed reveal a simple model underlying tumor progression and metastatic seeding. A comparison of Sgootr against alternative approaches shows that Sgootr can construct lineage trees with fewer migration events and with more in concordance with the sequential-progression model of tumor evolution, with a running time a fraction of that used in prior studies. Lineage-informative CpG sites identified by Sgootr are in inter-CpG island (CGI) regions, as opposed to intra-CGIs, which have been the main regions of interest in genomic methylation-related analyses.

Precise microdissection of gastric mixed adeno-neuroendocrine carcinoma dissects its genomic landscape and evolutionary clonal origins

Gastric mixed adenoneuroendocrine carcinoma (MANEC) is a clinically aggressive and heterogeneous tumor composed of adenocarcinoma (ACA) and neuroendocrine carcinoma (NEC). The genomic properties and evolutionary clonal origins of MANEC remain unclear. We conduct whole-exome and multiregional sequencing on 101 samples from 33 patients to elucidate their evolutionary paths. We identify four significantly mutated genes, TP53, RB1, APC, and CTNNB1. MANEC resembles chromosomal instability stomach adenocarcinoma in that whole-genome doubling in MANEC is predominant and occurs earlier than most copy-number losses. All tumors are of monoclonal origin, and NEC components show more aggressive genomic properties than their ACA counterparts. The phylogenetic trees show two tumor divergence patterns, including sequential and parallel divergence. Furthermore, ACA-to-NEC rather than NEC-to-ACA transition is confirmed by immunohistochemistry on 6 biomarkers in ACA- and NEC-dominant regions. These results provide insights into the clonal origin and tumor differentiation of MANEC.

Novel Genetic Subtypes of Urothelial Carcinoma With Differential Outcomes on Immune Checkpoint Blockade

PURPOSE

Immune checkpoint blockade (ICB) therapy has significantly improved clinical outcomes in bladder cancer. Identification of correlates of benefit is critical to select appropriate therapy for individual patients.

METHODS

To reveal genetic variables associated with benefit from ICB, we performed whole-exome sequencing on tumor specimens from 88 patients with advanced bladder cancer treated with ICB.

RESULTS

We identified several genetic factors that correlated with progression-free and overall survival after ICB therapy including ARID1A mutation, tumor mutational burden, intratumoral heterogeneity, the ratio of nonsynonymous to synonymous mutations in the immunopeptidome (immune dN/dS), and tumor cell purity. In addition, we noted that neutrophil-to-lymphocyte ratio and smoking history were negatively associated with overall survival. These genetic characteristics define four molecular subtypes demonstrating differential sensitivity to ICB. We validated the association of these four subtypes with clinical benefit from ICB in an independent cohort (IMvigor210). Finally, we showed that these molecular subtypes also correlate with outcome, although with distinct relationships, among patients not treated with ICB using The Cancer Genome Atlas (TCGA) bladder cancer cohort. Using parallel RNA sequencing data, the subtypes were also shown to correlate with immune infiltration and inflammation, respectively, in the IMvigor210 and TCGA cohorts.

CONCLUSION

Together, our study defines molecular subgroups of bladder cancer that influence benefit from ICB.

Inferring evolutionary trajectories from cross-sectional transcriptomic data to mirror lung adenocarcinoma progression

Lung adenocarcinoma (LUAD) is a deadly tumor with dynamic evolutionary process. Although much endeavors have been made in identifying the temporal patterns of cancer progression, it remains challenging to infer and interpret the molecular alterations associated with cancer development and progression. To this end, we developed a computational approach to infer the progression trajectory based on cross-sectional transcriptomic data. Analysis of the LUAD data using our approach revealed a linear trajectory with three different branches for malignant progression, and the results showed consistency in three independent cohorts. We used the progression model to elucidate the potential molecular events in LUAD progression. Further analysis showed that overexpression of BUB1B, BUB1 and BUB3 promoted tumor cell proliferation and metastases by disturbing the spindle assembly checkpoint (SAC) in the mitosis. Aberrant mitotic spindle checkpoint signaling appeared to be one of the key factors promoting LUAD progression. We found the inferred cancer trajectory allows to identify LUAD susceptibility genetic variations using genome-wide association analysis. This result shows the opportunity for combining analysis of candidate genetic factors with disease progression. Furthermore, the trajectory showed clear evident mutation accumulation and clonal expansion along with the LUAD progression. Understanding how tumors evolve and identifying mutated genes will help guide cancer management. We investigated the clonal architectures and identified distinct clones and subclones in different LUAD branches. Validation of the model in multiple independent data sets and correlation analysis with clinical results demonstrate that our method is effective and unbiased.

An unusual ectopic thymoma clonal evolution analysis: A case report

Thymomas and thymic carcinomas are rare and primary tumors of the mediastinum which is derived from the thymic epithelium. Thymomas are the most common primary anterior mediastinal tumor, while ectopic thymomas are rarer. Mutational profiles of ectopic thymomas may help expand our understanding of the occurrence and treatment options of these tumors. In this report, we sought to elucidate the mutational profiles of two ectopic thymoma nodules to gain deeper understanding of the molecular genetic information of this rare tumor and to provide guidance treatment options. We presented a case of 62-year-old male patient with a postoperative pathological diagnosis of type A mediastinal thymoma and ectopic pulmonary thymoma. After mediastinal lesion resection and thoracoscopic lung wedge resection, the mediastinal thymoma was completely removed, and the patient recovered from the surgery and no recurrence was found by examination until now. Whole exome sequencing was performed on both mediastinal thymoma and ectopic pulmonary thymoma tissue samples of the patient and clonal evolution analysis were further conducted to analyze the genetic characteristics. We identified eight gene mutations that were co-mutated in both lesions. Consistent with a previous exome sequencing analysis of thymic epithelial tumor, HRAS was also observed in both mediastinal lesion and lung lesion tissues. We also evaluated the intratumor heterogeneity of non-silent mutations. The results showed that the mediastinal lesion tissue has higher degree of heterogeneity and the lung lesion tissue has relatively low amount of variant heterogeneity in the detected variants. Through pathology and genomics sequencing detection, we initially revealed the genetic differences between mediastinal thymoma and ectopic thymoma, and clonal evolution analysis showed that these two lesions originated from multi-ancestral regions.

Inhibition of lysyl oxidases synergizes with 5-azacytidine to restore erythropoiesis in myelodysplastic and myeloid malignancies

Limited response rates and frequent relapses during standard of care with hypomethylating agents in myelodysplastic neoplasms (MN) require urgent improvement of this treatment indication. Here, by combining 5-azacytidine (5-AZA) with the pan-lysyl oxidase inhibitor PXS-5505, we demonstrate superior restoration of erythroid differentiation in hematopoietic stem and progenitor cells (HSPCs) of MN patients in 20/31 cases (65%) versus 9/31 cases (29%) treated with 5-AZA alone. This effect requires direct contact of HSPCs with bone marrow stroma components and is dependent on integrin signaling. We further confirm these results in vivo using a bone marrow niche-dependent MN xenograft model in female NSG mice, in which we additionally demonstrate an enforced reduction of dominant clones as well as significant attenuation of disease expansion and normalization of spleen sizes. Overall, these results lay out a strong pre-clinical rationale for efficacy of combination treatment of 5-AZA with PXS-5505 especially for anemic MN.

Reconstructing Clonal Evolution-A Systematic Evaluation of Current Bioinformatics Approaches

The accurate reconstruction of clonal evolution, including the identification of newly developing, highly aggressive subclones, is essential for the application of precision medicine in cancer treatment. Reconstruction, aiming for correct variant clustering and clonal evolution tree reconstruction, is commonly performed by tedious manual work. While there is a plethora of tools to automatically generate reconstruction, their reliability, especially reasons for unreliability, are not systematically assessed. We developed clevRsim-an approach to simulate clonal evolution data, including single-nucleotide variants as well as (overlapping) copy number variants. From this, we generated 88 data sets and performed a systematic evaluation of the tools for the reconstruction of clonal evolution. The results indicate a major negative influence of a high number of clones on both clustering and tree reconstruction. Low coverage as well as an extreme number of time points usually leads to poor clustering results. An underlying branched independent evolution hampers correct tree reconstruction. A further major decline in performance could be observed for large deletions and duplications overlapping single-nucleotide variants. In summary, to explore the full potential of reconstructing clonal evolution, improved algorithms that can properly handle the identified limitations are greatly needed.

Effective personalized neoantigen vaccine plus anti-PD-1 in a PD-1 blockade-resistant lung cancer patient

Background: Although significant progress has been made in immune checkpoint inhibitor (ICI) treatment of advanced squamous cell carcinoma (SqCC), most patients still experience acquired drug resistance. Methods: We used a dendritic cell-based neoantigen vaccine combined with ICIs to treat advanced SqCC in a PD-1 blockade-resistant patient. Results: The follow-up of this patient after 12 months revealed significant tumor regression. We also identified a new JAK1 ICI-resistant mutation that could become a potential universal neoantigen target for tumor vaccines. Conclusion: Individualized management of advanced SqCC through a combined neoantigen vaccine and ICI administration could yield beneficial clinical outcomes. Vaccines targeting anti-PD-1-resistant JAK1 mutations might be of particular benefit to a specific group of solid tumor patients.

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

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

Challenges in neoantigen-directed therapeutics

A fundamental prerequisite for the efficacy of cancer immunotherapy is the presence of functional, antigen-specific T cells within the tumor. Neoantigen-directed therapy is a promising strategy that aims at targeting the host's immune response against tumor-specific antigens, thereby eradicating cancer cells. Initial forays have been made in clinical environments utilizing vaccines and adoptive cell therapy; however, many challenges lie ahead. We provide an in-depth overview of the current state of the field with an emphasis on in silico neoantigen discovery and the clinical aspects that need to be addressed to unlock the full potential of this therapy.

A Random Forest Genomic Classifier for Tumor Agnostic Prediction of Response to Anti-PD1 Immunotherapy

Tumor mutational burden (TMB), a surrogate for tumor neoepitope burden, is used as a pan-tumor biomarker to identify patients who may benefit from anti-program cell death 1 (PD1) immunotherapy, but it is an imperfect biomarker. Multiple additional genomic characteristics are associated with anti-PD1 responses, but the combined predictive value of these features and the added informativeness of each respective feature remains unknown. We evaluated whether machine learning (ML) approaches using proposed determinants of anti-PD1 response derived from whole exome sequencing (WES) could improve prediction of anti-PD1 responders over TMB alone. Random forest classifiers were trained on publicly available anti-PD1 data (n = 104), and subsequently tested on an independent anti-PD1 cohort (n = 69). Both the training and test datasets included a range of cancer types such as non-small cell lung cancer (NSCLC), head and neck squamous cell carcinoma (HNSCC), melanoma, and smaller numbers of patients from other tumor types. Features used include summaries such as TMB and number of frameshift mutations, as well as more gene-level features such as counts of mutations associated with immune checkpoint response and resistance. Both ML algorithms demonstrated area under the receiver-operator curves (AUC) that exceeded TMB alone (AUC 0.63 "human-guided," 0.64 "cluster," and 0.58 TMB alone). Mutations within oncogenes disproportionately modulate anti-PD1 responses relative to their overall contribution to tumor neoepitope burden. The use of a ML algorithm evaluating multiple proposed genomic determinants of anti-PD1 responses modestly improves performance over TMB alone, highlighting the need to integrate other biomarkers to further improve model performance.

Characterization of genomic alterations and neoantigens and analysis of immune infiltration identified therapeutic and prognostic biomarkers in adenocarcinoma at the gastroesophageal junction

Objectives

Adenocarcinoma at the gastroesophageal junction (ACGEJ) refers to a malignant tumor that occurs at the esophagogastric junction. Despite some progress in targeted therapies for HER2, FGFR2, EGFR, MET, Claudin 18.2 and immune checkpoints in ACGEJ tumors, the 5-year survival rate of patients remains poor. Thus, it is urgent to explore genomic alterations and neoantigen characteristics of tumors and identify CD8+ T-cell infiltration-associated genes to find potential therapeutic targets and develop a risk model to predict ACGEJ patients’ overall survival (OS).

Methods

Whole-exome sequencing (WES) was performed on 55 paired samples from Chinese ACGEJ patients. Somatic mutations and copy number variations were detected by Strelka2 and FACETS, respectively. SigProfiler and SciClone were employed to decipher the mutation signature and clonal structure of each sample, respectively. Neoantigens were predicted using the MuPeXI pipeline. RNA sequencing (RNA-seq) data of ACGEJ samples from our previous studies and The Cancer Genome Atlas (TCGA) were used to identify genes significantly associated with CD8+ T-cell infiltration by weighted gene coexpression network analysis (WGCNA). To construct a risk model, we conducted LASSO and univariate and multivariate Cox regression analyses.

Results

Recurrent MAP2K7, RNF43 and RHOA mutations were found in ACGEJ tumors. The COSMIC signature SBS17 was associated with ACGEJ progression. CCNE1 and VEGFA were identified as putative CNV driver genes. PI3KCA and TP53 mutations conferred selective advantages to cancer cells. The Chinese ACGEJ patient neoantigen landscape was revealed for the first time, and 58 potential neoantigens common to TSNAdb and IEDB were identified. Compared with Siewert type II samples, Siewert type III samples had significant enrichment of the SBS17 signature, a lower TNFRSF14 copy number, a higher proportion of samples with complex clonal architecture and a higher neoantigen load. We identified 10 important CD8+ T-cell infiltration-related Hub genes (CCL5, CD2, CST7, GVINP1, GZMK, IL2RB, IKZF3, PLA2G2D, P2RY10 and ZAP70) as potential therapeutic targets from the RNA-seq data. Seven CD8+ T-cell infiltration-related genes (ADAM28, ASPH, CAMK2N1, F2R, STAP1, TP53INP2, ZC3H3) were selected to construct a prognostic model. Patients classified as high risk based on this model had significantly worse OS than low-risk patients, which was replicated in the TCGA-ACGEJ cohort.

Conclusions

This study provides new neoantigen-based immunotherapeutic targets for ACGEJ treatment and effective disease prognosis biomarkers.

Pharmacogenomic Profiling of Pediatric Acute Myeloid Leukemia to Identify Therapeutic Vulnerabilities and Inform Functional Precision Medicine

Despite the expanding portfolio of targeted therapies for adults with acute myeloid leukemia (AML), direct implementation in children is challenging due to inherent differences in underlying genetics. Here we established the pharmacologic profile of pediatric AML by screening myeloblast sensitivity to approved and investigational agents, revealing candidates of immediate clinical relevance. Drug responses ex vivo correlated with patient characteristics, exhibited age-specific alterations, and concorded with activities in xenograft models. Integration with genomic data uncovered new gene-drug associations, suggesting actionable therapeutic vulnerabilities. Transcriptome profiling further identified gene-expression signatures associated with on- and off-target drug responses. We also demonstrated the feasibility of drug screening-guided treatment for children with high-risk AML, with two evaluable cases achieving remission. Collectively, this study offers a high-dimensional gene-drug clinical data set that could be leveraged to research the unique biology of pediatric AML and sets the stage for realizing functional precision medicine for the clinical management of the disease.

SIGNIFICANCE

We conducted integrated drug and genomic profiling of patient biopsies to build the functional genomic landscape of pediatric AML. Age-specific differences in drug response and new gene-drug interactions were identified. The feasibility of functional precision medicine-guided management of children with high-risk AML was successfully demonstrated in two evaluable clinical cases. This article is highlighted in the In This Issue feature, p. 476.

Netie: inferring the evolution of neoantigen-T cell interactions in tumors

Neoantigens are the key targets of antitumor immune responses from cytotoxic T cells and play a critical role in affecting tumor progressions and immunotherapy treatment responses. However, little is known about how the interaction between neoantigens and T cells ultimately affects the evolution of cancerous masses. Here, we develop a hierarchical Bayesian model, named neoantigen-T cell interaction estimation (netie) to infer the history of neoantigen-CD8⁺ T cell interactions in tumors. Netie was systematically validated and applied to examine the molecular patterns of 3,219 tumors, compiled from a panel of 18 cancer types. We showed that tumors with an increase in immune selection pressure over time are associated with T cells that have an activation-related expression signature. We also identified a subset of exhausted cytotoxic T cells postimmunotherapy associated with tumor clones that newly arise after treatment. These analyses demonstrate how netie enables the interrogation of the relationship between individual neoantigen repertoires and the tumor molecular profiles. We found that a T cell inflammation gene expression profile (TIGEP) is more predictive of patient outcomes in the tumors with an increase in immune pressure over time, which reveals a curious synergy between T cells and neoantigen distributions. Overall, we provide a new tool that is capable of revealing the imprints left by neoantigens during each tumor's developmental process and of predicting how tumors will progress under further pressure of the host's immune system.

Bioinformatics roadmap for therapy selection in cancer genomics

Tumour heterogeneity is one of the main characteristics of cancer and can be categorised into inter- or intratumour heterogeneity. This heterogeneity has been revealed as one of the key causes of treatment failure and relapse. Precision oncology is an emerging field that seeks to design tailored treatments for each cancer patient according to epidemiological, clinical and omics data. This discipline relies on bioinformatics tools designed to compute scores to prioritise available drugs, with the aim of helping clinicians in treatment selection. In this review, we describe the current approaches for therapy selection depending on which type of tumour heterogeneity is being targeted and the available next-generation sequencing data. We cover intertumour heterogeneity studies and individual treatment selection using genomics variants, expression data or multi-omics strategies. We also describe intratumour dissection through clonal inference and single-cell transcriptomics, in each case providing bioinformatics tools for tailored treatment selection. Finally, we discuss how these therapy selection workflows could be integrated into the clinical practice.

Big data in basic and translational cancer research

Historically, the primary focus of cancer research has been molecular and clinical studies of a few essential pathways and genes. Recent years have seen the rapid accumulation of large-scale cancer omics data catalysed by breakthroughs in high-throughput technologies. This fast data growth has given rise to an evolving concept of 'big data' in cancer, whose analysis demands large computational resources and can potentially bring novel insights into essential questions. Indeed, the combination of big data, bioinformatics and artificial intelligence has led to notable advances in our basic understanding of cancer biology and to translational advancements. Further advances will require a concerted effort among data scientists, clinicians, biologists and policymakers. Here, we review the current state of the art and future challenges for harnessing big data to advance cancer research and treatment.

Reconstructing cancer phylogenies using Pairtree, a clone tree reconstruction algorithm

Pairtree is a clone tree reconstruction algorithm that uses somatic point mutations to build clone trees describing the evolutionary history of individual cancers. Using the Pairtree software package, we describe steps to preprocess somatic mutation data, cluster mutations into subclones, search for clone trees, and visualize clone trees. Pairtree builds clone trees using up to 100 samples from a single cancer with at least 30 subclonal populations. For complete details on the use and execution of this protocol, please refer to Wintersinger et al. (2022).