Absolute quantification of somatic DNA alterations in human cancer

Genomic profiling of mycosis fungoides identifies patients at high risk of disease progression

ABSTRACT

Mycosis fungoides (MF) is the most prevalent primary cutaneous T-cell lymphoma, with an indolent or aggressive course and poor survival. The pathogenesis of MF remains unclear, and prognostic factors in the early stages are not well established. Here, we characterized the most recurrent genomic alterations using whole-exome sequencing of 67 samples from 48 patients from Lille University Hospital (France), including 18 sequential samples drawn across stages of the malignancy. Genomic data were analyzed on the Broad Institute's Terra bioinformatics platform. We found that gain7q, gain10p15.1 (IL2RA and IL15RA), del10p11.22 (ZEB1), or mutations in JUNB and TET2 are associated with high-risk disease stages. Furthermore, gain7q, gain10p15.1 (IL2RA and IL15RA), del10p11.22 (ZEB1), and del6q16.3 (TNFAIP3) are coupled with shorter survival. Del6q16.3 (TNFAIP3) was a risk factor for progression in patients at low risk. By analyzing the clonal heterogeneity and the clonal evolution of the cohort, we defined different phylogenetic pathways of the disease with acquisition of JUNB, gain10p15.1 (IL2RA and IL15RA), or del12p13.1 (CDKN1B) at progression. These results establish the genomics and clonality of MF and identify potential patients at risk of progression, independent of their clinical stage.

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.

NOT gated T cells that selectively target EGFR and other widely expressed tumor antigens

Here, we show that a NOT gated cell therapy (Tmod) can exploit antigens such as epidermal growth factor receptor (EGFR) and human leukocyte antigen-E (HLA-E) which are widely expressed on cancer cells. Noncancerous cells-despite high expression of these antigens-are protected from cytotoxicity by the action of an inhibitory receptor ("blocker") via a mechanism that involves blocker modulation of CAR surface expression. The blocker is triggered by the product of a polymorphic HLA allele (e.g., HLA-A∗02) deleted in a significant subset of solid tumors via loss of heterozygosity. Moreover, Tmod constructs that target mouse homologs of EGFR or HLA-E for activation, and a mouse-equivalent of HLA-A∗02 for inhibition, protect mice from toxicity caused by the CAR alone. The blocker also controls graft vs. host response in allogeneic T cells in vitro, consistent with the use of Tmod cells for off-the-shelf therapy without additional gene-editing.

Microenvironment shapes small-cell lung cancer neuroendocrine states and presents therapeutic opportunities

Small-cell lung cancer (SCLC) is the most fatal form of lung cancer. Intratumoral heterogeneity, marked by neuroendocrine (NE) and non-neuroendocrine (non-NE) cell states, defines SCLC, but the cell-extrinsic drivers of SCLC plasticity are poorly understood. To map the landscape of SCLC tumor microenvironment (TME), we apply spatially resolved transcriptomics and quantitative mass spectrometry-based proteomics to metastatic SCLC tumors obtained via rapid autopsy. The phenotype and overall composition of non-malignant cells in the TME exhibit substantial variability, closely mirroring the tumor phenotype, suggesting TME-driven reprogramming of NE cell states. We identify cancer-associated fibroblasts (CAFs) as a crucial element of SCLC TME heterogeneity, contributing to immune exclusion, and predicting exceptionally poor prognosis. Our work provides a comprehensive map of SCLC tumor and TME ecosystems, emphasizing their pivotal role in SCLC's adaptable nature, opening possibilities for reprogramming the TME-tumor communications that shape SCLC tumor states.

CelFiE-ISH: a probabilistic model for multi-cell type deconvolution from single-molecule DNA methylation haplotypes

Deconvolution methods infer quantitative cell type estimates from bulk measurement of mixed samples including blood and tissue. DNA methylation sequencing measures multiple CpGs per read, but few existing deconvolution methods leverage this within-read information. We develop CelFiE-ISH, which extends an existing method (CelFiE) to use within-read haplotype information. CelFiE-ISH outperforms CelFiE and other existing methods, achieving 30% better accuracy and more sensitive detection of rare cell types. We also demonstrate the importance of marker selection and of tailoring markers for haplotype-aware methods. While here we use gold-standard short-read sequencing data, haplotype-aware methods will be well-suited for long-read sequencing.

Randomized Phase II Study Evaluating the Addition of Pembrolizumab to Radium-223 in Metastatic Castration-resistant Prostate Cancer

The checkpoint immunotherapeutic pembrolizumab induces responses in a small minority of patients with metastatic castration-resistant prostate cancer (mCRPC). Radium-223 (R223) may increase immunogenicity of bone metastases and increase pembrolizumab (P) activity. In a randomized phase II study, we assessed the effect of R223+P compared with R223 on tumor immune infiltration, safety, and clinical outcomes in patients with mCRPC. The primary endpoint was differences in CD4+ and CD8+ T-cell infiltrate in 8-week versus baseline bone metastasis biopsies; secondary endpoints were safety, radiographic progression-free survival (rPFS), and overall survival (OS). Of the 42 treated patients (29 R223+P, 13 R223), 18 R223+P and 8 R223 patients had evaluable paired tumor biopsies. Median fold-change of CD4+ T cells was -0.7 (range: -9.3 to 4.7) with R223+P and 0.1 (-11.1 to 3.7) with R223 (P = 0.66); for CD8+ T cells, median fold-change was -0.6 (-7.4 to 5.3) with R223+P and -1.3 (-3.1 to 4.8) with R223 (P = 0.66). Median rPFS and OS was 6.1 (95% confidence interval: 2.7-11.0) and 16.9 months [12.7-not reached (NR)], respectively, with R223+P and 5.7 (2.6-NR) and 16.0 (9.0-NR), respectively, with R223. Although R223+P was well tolerated with no unexpected toxicity, the combination did not improve efficacy. High-dimensional flow cytometry demonstrated minimal immune modulation with R223, whereas R223+P induced CTLA-4 expression on circulating CD4+ T cells. Clinical responders possessed lower circulating frequencies of Ki67+ T and myeloid cells at baseline and higher circulating frequencies of TIM-3+ T and myeloid cells by week 9. Although R223+P did not induce T-cell infiltration into the tumor microenvironment, exhaustion of induced peripheral T-cell immune responses may dampen the combination's clinical activity.

Adenoid cystic carcinoma of the Bartholin's gland is underpinned by MYB- and MYBL1- rearrangements

OBJECTIVE

Adenoid cystic carcinoma (AdCC) of the Bartholin's gland (AdCC-BG) is a very rare gynecologic vulvar malignancy. AdCC-BGs are slow-growing but locally aggressive and are associated with high recurrence rates. Here we sought to characterize the molecular underpinning of AdCC-BGs.

METHODS

AdCC-BGs (n = 6) were subjected to a combination of RNA-sequencing, targeted DNA-sequencing, reverse-transcription PCR, fluorescence in situ hybridization (FISH) and MYB immunohistochemistry (IHC). Clinicopathologic variables, somatic mutations, copy number alterations and chimeric transcripts were assessed.

RESULTS

All six AdCC-BGs were biphasic, composed of ductal and myoepithelial cells. Akin to salivary gland and breast AdCCs, three AdCC-BGs had the MYB::NFIB fusion gene with varying breakpoints, all of which were associated with MYB overexpression by IHC. Two AdCC-BGs were underpinned by MYBL1 fusion genes with different gene partners, including MYBL1::RAD51B and MYBL1::EWSR1 gene fusions, and showed MYB protein expression. Although the final AdCC-BG studied had MYB protein overexpression, no gene fusion was identified. AdCC-BGs harbored few additional somatic genetic alterations, and only few mutations in cancer-related genes were identified, including GNAQ, GNAS, KDM6A, AKT1 and BCL2, none of which were recurrent. Two AdCC-BGs, both with a MYB::NFIB fusion gene, developed metastatic disease.

CONCLUSIONS

AdCC-BGs constitute a convergent phenotype, whereby activation of MYB or MYBL1 can be driven by the MYB::NFIB fusion gene or MYBL1 rearrangements. Our observations further support the notion that AdCCs, irrespective of organ site, constitute a genotypic-phenotypic correlation. Assessment of MYB or MYBL1 rearrangements may be used as an ancillary marker for the diagnosis of AdCC-BGs.

Neural-net-based cell deconvolution from DNA methylation reveals tumor microenvironment associated with cancer prognosis

DNA methylation is a pivotal epigenetic modification that defines cellular identity. While cell deconvolution utilizing this information is considered useful for clinical practice, current methods for deconvolution are limited in their accuracy and resolution. In this study, we collected DNA methylation data from 945 human samples derived from various tissues and tumor-infiltrating immune cells and trained a neural network model with them. The model, termed MEnet, predicted abundance of cell population together with the detailed immune cell status from bulk DNA methylation data, and showed consistency to those of flow cytometry and histochemistry. MEnet was superior to the existing methods in the accuracy, speed, and detectable cell diversity, and could be applicable for peripheral blood, tumors, cell-free DNA, and formalin-fixed paraffin-embedded sections. Furthermore, by applying MEnet to 72 intrahepatic cholangiocarcinoma samples, we identified immune cell profiles associated with cancer prognosis. We believe that cell deconvolution by MEnet has the potential for use in clinical settings.

Reconstructing phylogenetic trees from genome-wide somatic mutations in clonal samples

Phylogenetic trees are a powerful means to display the evolutionary history of species, pathogens and, more recently, individual cells of the human body. Whole-genome sequencing of laser capture microdissections or expanded stem cells has allowed the discovery of somatic mutations in clones, which can be used as natural barcodes to reconstruct the developmental history of individual cells. Here we describe Sequoia, our pipeline to reconstruct lineage trees from clones of normal cells. Candidate somatic mutations are called against the human reference genome and filtered to exclude germline mutations and artifactual variants. These filtered somatic mutations form the basis for phylogeny reconstruction using a maximum parsimony framework. Lastly, we use a maximum likelihood framework to explicitly map mutations to branches in the phylogenetic tree. The resulting phylogenies can then serve as a basis for many subsequent analyses, including investigating embryonic development, tissue dynamics in health and disease, and mutational signatures. Sequoia can be readily applied to any clonal somatic mutation dataset, including single-cell DNA sequencing datasets, using the commands and scripts provided. Moreover, Sequoia is highly flexible and can be easily customized. Typically, the runtime of the core script ranges from minutes to an hour for datasets with a moderate number (50,000-150,000) of variants. Competent bioinformatic skills, including in-depth knowledge of the R programming language, are required. A high-performance computing cluster (one that is capable of running mutation-calling algorithms and other aspects of the analysis at scale) is also required, especially if handling large datasets.

Stratification of glioma based on stemness scores in bulk and single-cell transcriptomes

BACKGROUND

Brain tumours are known to have a high mortality and morbidity rate due to their localised and frequent invasive growth. The concept that glioma resistance could originate from the dissimilarity in the vulnerability of clonogenic glial stem cells to chemotherapeutic drugs and radiation has driven the scientific community to reexamine the comprehension of glioma growth and strategies that target these cells or modify their stemness.

METHODS

Based on the enrichment scores of 12 stemness signatures, we identified glioma subtypes in both tumour bulks and single cells by clustering analysis. Furthermore, we comprehensively compared molecular and clinical features among the glioma subtypes.

RESULTS

Consistently, in seven different datasets, hierarchical clustering uncovered three subtypes of glioma, termed Stem-H, Stem-M, and Stem-L, with high, medium, and low stemness signatures, respectively. Stem-H and Stem-L exhibited the most unfavorable and favourable overall and disease-free survival, respectively. Stem-H showed the highest enrichment scores of the EMT, invasion, proliferation, differentiation, and metastasis processes signatures, while Stem-L displayed the lowest. Stem-H harboured a greater proportion of late-stage tumours compared to Stem-L. Moreover, Stem-H manifested higher tumour mutation burden, DNA damage repair and cell cycle activity, intratumour heterogeneity, and a more frequent incidence of TP53 and EGFR mutations than Stem-L. In contrast, Stem-L had higher O6-Methylguanine-DNA Methyltransferase (MGMT) methylation levels.

CONCLUSION

The classification of glioma based on stemness may offer new insights into the biology of the tumour, as well as more accurate clinical management of the disease.

A prognostic neural epigenetic signature in high-grade glioma

Neural-tumor interactions drive glioma growth as evidenced in preclinical models, but clinical validation is limited. We present an epigenetically defined neural signature of glioblastoma that independently predicts patients' survival. We use reference signatures of neural cells to deconvolve tumor DNA and classify samples into low- or high-neural tumors. High-neural glioblastomas exhibit hypomethylated CpG sites and upregulation of genes associated with synaptic integration. Single-cell transcriptomic analysis reveals a high abundance of malignant stemcell-like cells in high-neural glioblastoma, primarily of the neural lineage. These cells are further classified as neural-progenitor-cell-like, astrocyte-like and oligodendrocyte-progenitor-like, alongside oligodendrocytes and excitatory neurons. In line with these findings, high-neural glioblastoma cells engender neuron-to-glioma synapse formation in vitro and in vivo and show an unfavorable survival after xenografting. In patients, a high-neural signature is associated with decreased overall and progression-free survival. High-neural tumors also exhibit increased functional connectivity in magnetencephalography and resting-state magnet resonance imaging and can be detected via DNA analytes and brain-derived neurotrophic factor in patients' plasma. The prognostic importance of the neural signature was further validated in patients diagnosed with diffuse midline glioma. Our study presents an epigenetically defined malignant neural signature in high-grade gliomas that is prognostically relevant. High-neural gliomas likely require a maximized surgical resection approach for improved outcomes.

Ultrasensitive plasma-based monitoring of tumor burden using machine-learning-guided signal enrichment

In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole-genome sequencing (WGS). Here we now introduce MRD-EDGE, a machine-learning-guided WGS ctDNA single-nucleotide variant (SNV) and copy-number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by ~300× compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, vastly expanding its applicability within solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in TF in response to neoadjuvant immunotherapy in lung cancer and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables plasma-only (non-tumor-informed) disease monitoring in advanced melanoma and lung cancer, yielding clinically informative TF monitoring for patients on immune-checkpoint inhibition.

Genomic characterization and risk stratification of esophageal squamous dysplasia

Objectives

The majority of esophageal squamous dysplasia (ESD) patients progress slowly, while a subset of patients can undergo recurrence rapidly or progress to invasive cancer even after proper treatment. However, the molecular mechanisms underlying these clinical observations are still largely unknown.

Methods

By sequencing the genomic data of 160 clinical samples from 49 tumor-free ESD patients and 88 esophageal squamous cell carcinoma (ESCC) patients, we demonstrated lower somatic mutation and copy number alteration (CNA) burden in ESD compared with ESCC.

Results

Cross-species screening and functional assays identified ACSM5 as a novel driver gene for ESD progression. Furthermore, we revealed that miR-4292 promoted ESD progression and could serve as a non-invasive diagnostic marker for ESD.

Conclusions

These findings largely expanded our understanding of ESD genetics and tumorigenesis, which possessed promising significance for improving early diagnosis, reducing overtreatment, and identifying high-risk ESD patients.

Comprehensive proteogenomic characterization of rare kidney tumors

Non-clear cell renal cell carcinomas (non-ccRCCs) encompass diverse malignant and benign tumors. Refinement of differential diagnosis biomarkers, markers for early prognosis of aggressive disease, and therapeutic targets to complement immunotherapy are current clinical needs. Multi-omics analyses of 48 non-ccRCCs compared with 103 ccRCCs reveal proteogenomic, phosphorylation, glycosylation, and metabolic aberrations in RCC subtypes. RCCs with high genome instability display overexpression of IGF2BP3 and PYCR1. Integration of single-cell and bulk transcriptome data predicts diverse cell-of-origin and clarifies RCC subtype-specific proteogenomic signatures. Expression of biomarkers MAPRE3, ADGRF5, and GPNMB differentiates renal oncocytoma from chromophobe RCC, and PIGR and SOSTDC1 distinguish papillary RCC from MTSCC. This study expands our knowledge of proteogenomic signatures, biomarkers, and potential therapeutic targets in non-ccRCC.

Copy number alterations: a catastrophic orchestration of the breast cancer genome

Breast cancer (BCa) is a prevalent malignancy that predominantly affects women around the world. Somatic copy number alterations (CNAs) are tumor-specific amplifications or deletions of DNA segments that often drive BCa development and therapy resistance. Hence, the complex patterns of CNAs complement BCa classification systems. In addition, understanding the precise contributions of CNAs is essential for tailoring personalized treatment approaches. This review highlights how tumor evolution drives the acquisition of CNAs, which in turn shape the genomic landscapes of BCas. It also discusses advanced methodologies for identifying recurrent CNAs, studying CNAs in BCa and their clinical impact.

Molecular alterations in claudin 18 suppressed and non-suppressed gastric adenocarcinomas to guide targeted therapies

BACKGROUND

Gastric adenocarcinoma represents an aggressive type of cancer and an important cause of cancer mortality. Progress in gastric cancer therapeutics has resulted from a better understanding of the molecular pathogenesis of the disease and introduction of targeted therapies, but most gastric cancer patients still rely on non-targeted chemotherapy as the mainstay of treatment for advanced disease.

METHODS

An analysis of publicly available series from The Cancer Genome Atlas (TCGA) gastric cancer cohort was undertaken to delineate the clinical and genomic landscape of gastric cancers with suppressed expression of claudin 18 compared with cancers with non-suppressed claudin 18. Claudin 18 suppressed cancers were defined as having an mRNA expression z-score relative to normal samples (log RNA Seq V2) of less than -1. Claudin 18 non-suppressed cancers were defined as having an mRNA expression z-score relative to normal samples (log RNA Seq V2) above 0.5.

RESULTS

Gastric cancers with claudin 18 mRNA suppression represented 7.7% of the gastric adenocarcinomas of TCGA cohort, while non-suppressed cancers represented 46.6% of the cases. The two groups did not differ in clinical and genomic characteristics, such as mean age, histology, grade, and stage. The mutation landscape of claudin 18 suppressed cases included high mutation rates of TP53, of genes of the WNT/β-catenin pathway and of ubiquitin ligase FBXW7. Moreover, a subset of both claudin 18 suppressed and non-suppressed cancers displayed mutations in Mismatch Repair (MMR) associated genes or a high tumor mutation burden (TMB). At the mRNA expression level, claudin 18 suppressed gastric cancers showed up-regulation of EMT core transcription factor Snail 2 and down-regulation of genes of HLA cluster. The survival of gastric cancer patients with claudin 18 mRNA suppression was not significantly different compared with patients with non-suppressed claudin 18.

CONCLUSION

Sub-sets of gastric cancers with claudin 18 mRNA suppression displayed characteristics of potential therapeutic interest, such as mutations in WNT and PI3K pathways and MMR defects. These may guide the development of alternative targeted therapies, in this sub-set of gastric cancers which are not candidates for claudin 18 targeting therapies.

Classification of multiple primary lung cancer in patients with multifocal lung cancer: assessment of a machine learning approach using multidimensional genomic data

Background

Multiple primary lung cancer (MPLC) is an increasingly well-known clinical phenomenon. However, its molecular characterizations are poorly understood, and still lacks of effective method to distinguish it from intrapulmonary metastasis (IM). Herein, we propose an identification model based on molecular multidimensional analysis in order to accurately optimize treatment.

Methods

A total of 112 Chinese lung cancers harboring at least two tumors (n = 270) were enrolled. We retrospectively selected 74 patients with 121 tumor pairs and randomly divided the tumor pairs into a training cohort and a test cohort in a 7:3 ratio. A novel model was established in training cohort, optimized for MPLC identification using comprehensive genomic profiling analyzed by a broad panel with 808 cancer-related genes, and evaluated in the test cohort and a prospective validation cohort of 38 patients with 112 tumors.

Results

We found differences in molecular characterizations between the two diseases and rigorously selected the characterizations to build an identification model. We evaluated the performance of the classifier using the test cohort data and observed an 89.5% percent agreement (PA) for MPLC and a 100.0% percent agreement for IM. The model showed an excellent area under the curve (AUC) of 0.947 and a 91.3% overall accuracy. Similarly, the assay achieved a considerable performance in the independent validation set with an AUC of 0.938 and an MPLC predictive value of 100%. More importantly, the MPLC predictive value of the classification achieved 100% in both the test set and validation cohort. Compared to our previous mutation-based method, the classifier showed better κ consistencies with clinical classification among all 112 patients (0.84 vs. 0.65, p <.01).

Conclusion

These data provide novel evidence of MPLC-specific genomic characteristics and demonstrate that our one-step molecular classifier can accurately classify multifocal lung tumors as MPLC or IM, which suggested that broad panel NGS may be a useful tool for assisting with differential diagnoses.

Assay Validation of Cell-Free DNA Shallow Whole-Genome Sequencing to Determine Tumor Fraction in Advanced Cancers

Blood-based liquid biopsy is increasingly used in clinical care of patients with cancer, and fraction of tumor-derived DNA in circulation (tumor fraction; TFx) has demonstrated clinical validity across multiple cancer types. To determine TFx, shallow whole-genome sequencing of cell-free DNA (cfDNA) can be performed from a single blood sample, using an established computational pipeline (ichorCNA), without prior knowledge of tumor mutations, in a highly cost-effective manner. We describe assay validation of this approach to facilitate broad clinical application, including evaluation of assay sensitivity, precision, repeatability, reproducibility, pre-analytic factors, and DNA quality/quantity. Sensitivity to detect TFx of 3% (lower limit of detection) was 97.2% to 100% at 1× and 0.1× mean sequencing depth, respectively. Precision was demonstrated on distinct sequencing instruments (HiSeqX and NovaSeq) with no observable differences. The assay achieved prespecified 95% agreement of TFx across replicates of the same specimen (repeatability) and duplicate samples in different batches (reproducibility). Comparison of samples collected in EDTA and Streck tubes from single venipuncture in 23 patients demonstrated that EDTA or Streck tubes were comparable if processed within 8 hours. On the basis of a range of DNA inputs (1 to 50 ng), 20 ng cfDNA is the preferred input, with 5 ng minimum acceptable. Overall, this shallow whole-genome sequencing of cfDNA and ichorCNA approach offers sensitive, precise, and reproducible quantitation of TFx, facilitating assay application in clinical cancer care.

Identification of CNKSR1 as a biomarker for "cold" tumor microenvironment in lung adenocarcinoma: An integrative analysis based on a novel workflow

Background

Therapies targeting PD1/PD-L1 pathway have revolutionized the treatment of lung cancer. However, anti-PD1/PD-L1 therapies have proven beneficial for only a select group of lung adenocarcinoma (LUAD) patients and generally do not work for immuno-cold tumors characterized by a lack of immune cell infiltration. Identifying novel biomarkers is vital to broad therapeutic options for LUAD patients with no response to anti-PD1/PD-L1 immunotherapies.

Methods

Our study has developed a novel strategy to identify a promising biomarker that addresses the limitations of anti-PD1/PD-L1 immunotherapy in treating immunological cold tumors. We exacted LUAD RNA-seq data from the Cancer Genome Atlas database (TCGA). Using several machine learning methods, we identified the candidate biomarker. Based on the expression level of PD-L1 and the identified biomarker, samples were categorized into four groups. We further used ESTIMATE, ssGSEA, and CIBERSORT algorithms to calculate the immune infiltration level of each group. The results were validated in three independent bulk datasets and one scRNA-seq dataset. Immunohistochemistry (IHC) assessments were performed in clinical samples to further evaluate the coexpression of CNKSR1 and PD-L1, and to compare CD8 + T cell infiltration among groups.

Results

After comprehensive analyses, CNKSR1 was identified as a novel promising biomarker for immuno-cold LUAD. CNKSR1 mRNA expression levels exhibited a negative correlation with both PD-L1 mRNA expression and the extent of immune cell infiltration in LUAD. Besides, in contrast to the significant association between the expression of PD-L1 and the majority of other well-established or widely studied immune checkpoint molecules, a mutually exclusive expression pattern is observed between CNKSR1 and these molecules. The aforementioned results were consistent in validation datasets. The prognostic model built based on the CNKSR1 coexpression module also showed robust predictive performance. Additionally, IHC assessments have confirmed that the coexpression of CNKSR1 and PD-L1 is rare in LUAD samples. Notably, LUADs in the high-CNKSR1 group, characterized by high CNKSR1 but low PD- L1 expression, demonstrated reduced infiltration of CD8+ T cells.

Conclusions

In summary, CNKSR1 emerges as a promising biomarker for immune-cold LUADs, and the study into CNKSR1 modulating T-cell infiltration may lead to the identification of compensatory molecules to enhance the effectiveness of current immunotherapy for LUAD.

Long-Term Follow-Up Defines the Population That Benefits from Early Interception in a High-Risk Smoldering Multiple Myeloma Clinical Trial Using the Combination of Ixazomib, Lenalidomide, and Dexamethasone

Background

Early therapeutic intervention in high-risk SMM (HR-SMM) has demonstrated benefit in previous studies of lenalidomide with or without dexamethasone. Triplets and quadruplet studies have been examined in this same population. However, to date, none of these studies examined the impact of depth of response on long-term outcomes of participants treated with lenalidomide-based therapy, and whether the use of the 20/2/20 model or the addition of genomic alterations can further define the population that would benefit the most from early therapeutic intervention. Here, we present the results of the phase II study of the combination of ixazomib, lenalidomide, and dexamethasone in patients with HR-SMM with long-term follow-up and baseline single-cell tumor and immune sequencing that help refine the population to be treated for early intervention studies.

Methods

This is a phase II trial of ixazomib, lenalidomide, and dexamethasone (IRD) in HR-SMM. Patients received 9 cycles of induction therapy with ixazomib 4mg on days 1, 8, and 15; lenalidomide 25mg on days 1-21; and dexamethasone 40mg on days 1, 8, 15, and 22. The induction phase was followed by maintenance with ixazomib 4mg on days 1, 8, and 15; and lenalidomide 15mg d1-21 for 15 cycles for 24 months of treatment. The primary endpoint was progression-free survival after 2 years of therapy. Secondary endpoints included depth of response, biochemical progression, and correlative studies included single-cell RNA sequencing and/or whole-genome sequencing of the tumor and single-cell sequencing of immune cells at baseline.

Results

Fifty-five patients, with a median age of 64, were enrolled in the study. The overall response rate was 93%, with 31% of patients achieving a complete response and 45% achieving a very good partial response or better. The most common grade 3 or greater treatment-related hematologic toxicities were neutropenia (16 patients; 29%), leukopenia (10 patients; 18%), lymphocytopenia (8 patients; 15%), and thrombocytopenia (4 patients; 7%). Non-hematologic grade 3 or greater toxicities included hypophosphatemia (7 patients; 13%), rash (5 patients; 9%), and hypokalemia (4 patients; 7%). After a median follow-up of 50 months, the median progression-free survival (PFS) was 48.6 months (95% CI: 39.9 - not reached; NR) and median overall survival has not been reached. Patients achieving VGPR or better had a significantly better progression-free survival (p<0.001) compared to those who did not achieve VGPR (median PFS 58.2 months vs. 31.3 months). Biochemical progression preceded or was concurrent with the development of SLiM-CRAB criteria in eight patients during follow-up, indicating that biochemical progression is a meaningful endpoint that correlates with the development of end-organ damage. High-risk 20/2/20 participants had the worst PFS compared to low- and intermediate-risk participants. The use of whole genome or single-cell sequencing of tumor cells identified high-risk aberrations that were not identified by FISH alone and aided in the identification of participants at risk of progression. scRNA-seq analysis revealed a positive correlation between MHC class I expression and response to proteasome inhibition and at the same time a decreased proportion of GZMB+ T cells within the clonally expanded CD8+ T cell population correlated with suboptimal response.

Conclusions

Ixazomib, lenalidomide and dexamethasone in HR-SMM demonstrates significant clinical activity with an overall favorable safety profile. Achievement of VGPR or greater led to significant improvement in time to progression, suggesting that achieving deep response is beneficial in HR-SMM. Biochemical progression correlates with end-organ damage. Patients with high-risk FISH and lack of deep response had poor outcomes. ClinicalTrials.gov identifier: ( NCT02916771 ).

Multicellular ecotypes shape progression of lung adenocarcinoma from ground-glass opacity toward advanced stages

Lung adenocarcinoma is a type of cancer that exhibits a wide range of clinical radiological manifestations, from ground-glass opacity (GGO) to pure solid nodules, which vary greatly in terms of their biological characteristics. Our current understanding of this heterogeneity is limited. To address this gap, we analyze 58 lung adenocarcinoma patients via machine learning, single-cell RNA sequencing (scRNA-seq), and whole-exome sequencing, and we identify six lung multicellular ecotypes (LMEs) correlating with distinct radiological patterns and cancer cell states. Notably, GGO-associated neoantigens in early-stage cancers are recognized by CD8+ T cells, indicating an immune-active environment, while solid nodules feature an immune-suppressive LME with exhausted CD8+ T cells, driven by specific stromal cells such as CTHCR1+ fibroblasts. This study also highlights EGFR(L858R) neoantigens in GGO samples, suggesting potential CD8+ T cell activation. Our findings offer valuable insights into lung adenocarcinoma heterogeneity, suggesting avenues for targeted therapies in early-stage disease.

Differential whole-genome doubling based signatures for improvement on clinical outcomes and drug response in patients with breast cancer

Whole genome doublings (WGD), a hallmark of human cancer, is pervasive in breast cancer patients. However, the molecular mechanism of the complete impact of WGD on survival and treatment response in breast cancer remains unclear. To address this, we performed a comprehensive and systematic analysis of WGD, aiming to identify distinct genetic alterations linked to WGD and highlight its improvement on clinical outcomes and treatment response for breast cancer. A linear regression model along with weighted gene co-expression network analysis (WGCNA) was applied on The Cancer Genome Atlas (TCGA) dataset to identify critical genes related to WGD. Further Cox regression models with random selection were used to optimize the most useful prognostic markers in the TCGA dataset. The clinical implication of the risk model was further assessed through prognostic impact evaluation, tumor stratification, functional analysis, genomic feature difference analysis, drug response analysis, and multiple independent datasets for validation. Our findings revealed a high aneuploidy burden, chromosomal instability (CIN), copy number variation (CNV), and mutation burden in breast tumors exhibiting WGD events. Moreover, 247 key genes associated with WGD were identified from the distinct genomic patterns in the TCGA dataset. A risk model consisting of 22 genes was optimized from the key genes. High-risk breast cancer patients were more prone to WGD and exhibited greater genomic diversity compared to low-risk patients. Some oncogenic signaling pathways were enriched in the high-risk group, while primary immune deficiency pathways were enriched in the low-risk group. We also identified a risk gene, ANLN (anillin), which displayed a strong positive correlation with two crucial WGD genes, KIF18A and CCNE2. Tumors with high expression of ANLN were more prone to WGD events and displayed worse clinical survival outcomes. Furthermore, the expression levels of these risk genes were significantly associated with the sensitivities of BRCA cell lines to multiple drugs, providing valuable insights for targeted therapies. These findings will be helpful for further improvement on clinical outcomes and contribution to drug development in breast cancer.

Molecular patterns of resistance to immune checkpoint blockade in melanoma

Immune checkpoint blockade (ICB) has improved outcome for patients with metastatic melanoma but not all benefit from treatment. Several immune- and tumor intrinsic features are associated with clinical response at baseline. However, we need to further understand the molecular changes occurring during development of ICB resistance. Here, we collect biopsies from a cohort of 44 patients with melanoma after progression on anti-CTLA4 or anti-PD1 monotherapy. Genetic alterations of antigen presentation and interferon gamma signaling pathways are observed in approximately 25% of ICB resistant cases. Anti-CTLA4 resistant lesions have a sustained immune response, including immune-regulatory features, as suggested by multiplex spatial and T cell receptor (TCR) clonality analyses. One anti-PD1 resistant lesion harbors a distinct immune cell niche, however, anti-PD1 resistant tumors are generally immune poor with non-expanded TCR clones. Such immune poor microenvironments are associated with melanoma cells having a de-differentiated phenotype lacking expression of MHC-I molecules. In addition, anti-PD1 resistant tumors have reduced fractions of PD1+ CD8+ T cells as compared to ICB naïve metastases. Collectively, these data show the complexity of ICB resistance and highlight differences between anti-CTLA4 and anti-PD1 resistance that may underlie differential clinical outcomes of therapy sequence and combination.

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.

Leveraging new methods for comprehensive characterization of mitochondrial DNA in esophageal squamous cell carcinoma

BACKGROUND

Mitochondria play essential roles in tumorigenesis; however, little is known about the contribution of mitochondrial DNA (mtDNA) to esophageal squamous cell carcinoma (ESCC). Whole-genome sequencing (WGS) is by far the most efficient technology to fully characterize the molecular features of mtDNA; however, due to the high redundancy and heterogeneity of mtDNA in regular WGS data, methods for mtDNA analysis are far from satisfactory.

METHODS

Here, we developed a likelihood-based method dMTLV to identify low-heteroplasmic mtDNA variants. In addition, we described fNUMT, which can simultaneously detect non-reference nuclear sequences of mitochondrial origin (non-ref NUMTs) and their derived artifacts. Using these new methods, we explored the contribution of mtDNA to ESCC utilizing the multi-omics data of 663 paired tumor-normal samples.

RESULTS

dMTLV outperformed the existing methods in sensitivity without sacrificing specificity. The verification using Nanopore long-read sequencing data showed that fNUMT has superior specificity and more accurate breakpoint identification than the current methods. Leveraging the new method, we identified a significant association between the ESCC overall survival and the ratio of mtDNA copy number of paired tumor-normal samples, which could be potentially explained by the differential expression of genes enriched in pathways related to metabolism, DNA damage repair, and cell cycle checkpoint. Additionally, we observed that the expression of CBWD1 was downregulated by the non-ref NUMTs inserted into its intron region, which might provide precursor conditions for the tumor cells to adapt to a hypoxic environment. Moreover, we identified a strong positive relationship between the number of mtDNA truncating mutations and the contribution of signatures linked to tumorigenesis and treatment response.

CONCLUSIONS

Our new frameworks promote the characterization of mtDNA features, which enables the elucidation of the landscapes and roles of mtDNA in ESCC essential for extending the current understanding of ESCC etiology. dMTLV and fNUMT are freely available from https://github.com/sunnyzxh/dMTLV and https://github.com/sunnyzxh/fNUMT , respectively.

The correlation between clinical outcomes and genomic analysis with high risk factors for the progression of osteosarcoma

Osteosarcoma (OS) is a rare but aggressive malignancy. Despite previous reports, molecular characterization of this disease is not well understood, and little is known regarding OS in Chinese patients. Herein, we analyzed the genomic signatures of 73 Chinese OS cases. TP53, NCOR1, LRP1B, ATRX, RB1, and TFE3 were the most frequently mutated gene in our OS cohort. In addition, the genomic analysis of Western OS patients was performed. Notably, there were remarkable disparities in mutational landscape, base substitution pattern, and tumor mutational burden between the Chinese and Western OS cohorts. Specific molecular mechanisms, including DNA damage repair (DDR) gene mutations, copy number variation (CNV) presence, aneuploidy, and intratumoral heterogeneity, were associated with disease progression. Additionally, 30.1% of OS patients carried clinically actionable alterations, which were mainly enriched in PI3K, MAPK, DDR, and RTK signaling pathways. A specific molecular subtype incorporating DDR alterations and CNVs was significantly correlated with distant metastasis-free survival and event-free survival, and this correlation was observed in all subgroups with different characteristics. These findings comprehensively elucidated the genomic profile and revealed novel prognostic factors in OS, which would contribute to understanding this disease and promoting precision medicine of this population.

Identifying cancer subtypes based on embryonic and hematopoietic stem cell signatures in pan-cancer

PURPOSE

Cancer cells with stem cell-like properties may contribute to cancer development and therapy resistance. The advancement of multi-omics technology has sparked interest in exploring cancer stemness from a multi-omics perspective. However, there is a limited number of studies that have attempted to subtype cancer by combining different types of stem cell signatures.

METHODS

In this study, 10,323 cancer specimens from 33 TCGA cancer types were clustered based on the enrichment scores of six stemness gene sets, representing two types of stem cell backgrounds: embryonic stem cells (ESCs) and hematopoietic stem cells (HSCs).

RESULTS

We identified four subtypes of pan-cancer, termed StC1, StC2, StC3 and StC4, which displayed distinct molecular and clinical features, including stemness, genome integrity, intratumor heterogeneity, methylation levels, tumor microenvironment, tumor progression, responses to chemotherapy and immunotherapy, and survival prognosis. Importantly, this subtyping method for pan-cancer is reproducible at the protein level.

CONCLUSION

Our findings indicate that the ESC signature is an adverse prognostic factor in cancer, while the HSC signature and ratio of HSC/ESC signatures are positive prognostic factors. The subtyping of cancer based on ESC and HSC signatures may provide insights into cancer biology and clinical implications of cancer.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

CDK12 Loss Promotes Prostate Cancer Development While Exposing Vulnerabilities to Paralog-Based Synthetic Lethality

Biallelic loss of cyclin-dependent kinase 12 (CDK12) defines a unique molecular subtype of metastatic castration-resistant prostate cancer (mCRPC). It remains unclear, however, whether CDK12 loss per se is sufficient to drive prostate cancer development-either alone, or in the context of other genetic alterations-and whether CDK12-mutant tumors exhibit sensitivity to specific pharmacotherapies. Here, we demonstrate that tissue-specific Cdk12 ablation is sufficient to induce preneoplastic lesions and robust T cell infiltration in the mouse prostate. Allograft-based CRISPR screening demonstrated that Cdk12 loss is positively associated with Trp53 inactivation but negatively associated with Pten inactivation-akin to what is observed in human mCRPC. Consistent with this, ablation of Cdk12 in prostate organoids with concurrent Trp53 loss promotes their proliferation and ability to form tumors in mice, while Cdk12 knockout in the Pten-null prostate cancer mouse model abrogates tumor growth. Bigenic Cdk12 and Trp53 loss allografts represent a new syngeneic model for the study of androgen receptor (AR)-positive, luminal prostate cancer. Notably, Cdk12/Trp53 loss prostate tumors are sensitive to immune checkpoint blockade. Cdk12-null organoids (either with or without Trp53 co-ablation) and patient-derived xenografts from tumors with CDK12 inactivation are highly sensitive to inhibition or degradation of its paralog kinase, CDK13. Together, these data identify CDK12 as a bona fide tumor suppressor gene with impact on tumor progression and lends support to paralog-based synthetic lethality as a promising strategy for treating CDK12-mutant mCRPC.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

Exemestane plus everolimus and palbociclib in metastatic breast cancer: clinical response and genomic/transcriptomic determinants of resistance in a phase I/II trial

The landscape of cyclin-dependent kinase 4/6 inhibitor (CDK4/6i) resistance is still being elucidated and the optimal subsequent therapy to overcome resistance remains uncertain. Here we present the final results of a phase Ib/IIa, open-label trial (NCT02871791) of exemestane plus everolimus and palbociclib for CDK4/6i-resistant metastatic breast cancer. The primary objective of phase Ib was to evaluate safety and tolerability and determine the maximum tolerated dose/recommended phase II dose (100 mg palbociclib, 5 mg everolimus, 25 mg exemestane). The primary objective of phase IIa was to determine the clinical benefit rate (18.8%, n = 6/32), which did not meet the predefined endpoint (65%). Secondary objectives included pharmacokinetic profiling (phase Ib), objective response rate, disease control rate, duration of response, and progression free survival (phase IIa), and correlative multi-omics analysis to investigate biomarkers of resistance to CDK4/6i. All participants were female. Multi-omics data from the phase IIa patients (n = 24 tumor/17 blood biopsy exomes; n = 27 tumor transcriptomes) showed potential mechanisms of resistance (convergent evolution of HER2 activation, BRAFV600E), identified joint genomic/transcriptomic resistance features (ESR1 mutations, high estrogen receptor pathway activity, and a Luminal A/B subtype; ERBB2/BRAF mutations, high RTK/MAPK pathway activity, and a HER2-E subtype), and provided hypothesis-generating results suggesting that mTOR pathway activation correlates with response to the trial's therapy. Our results illustrate how genome and transcriptome sequencing may help better identify patients likely to respond to CDK4/6i therapies.

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.

Clinical value of FAT1 mutations to indicate the immune response in colorectal cancer patients

Immunotherapy is currently approved for CRC whose tumors have high MSI-H. To find additional biomarkers for immunotherapy in CRC, targeted sequencing was performed on tumor tissues from a discovery cohort of 161 CRC patients. Validation cohorts from the cBioPortal were also used for survival and tumor cell infiltration analyses. The FAT1-mutated CRC group often co-occurred with MSI events and displayed a higher tumor mutational burden compared to the FAT1 wild-type CRC. Overall survival was higher in patients with FAT1 mutations than in patients with wild type FAT1. The altered PI3K-AKT pathway and immune pathways were enriched in the FAT1-mutated CRC. A higher infiltration rate of immune cells including CD4+ T cells, CD8+ T cells, macrophages M1 and regulatory T cells were also observed in the colorectal tumors with FAT1 mutation compared to tumors with wild type FAT1. The results showed that CRC patients with FAT1 mutations exhibited an immunotherapy-favorable profile.

Whole-genome doubling is a double-edged sword: the heterogeneous role of whole-genome doubling in various cancer types

Whole-genome doubling (WGD), characterized by the duplication of an entire set of chromosomes, is commonly observed in various tumors, occurring in approximately 30-40% of patients with different cancer types. The effect of WGD on tumorigenesis varies depending on the context, either promoting or suppressing tumor progression. Recent advances in genomic technologies and large-scale clinical investigations have led to the identification of the complex patterns of genomic alterations underlying WGD and their functional consequences on tumorigenesis progression and prognosis. Our comprehensive review aims to summarize the causes and effects of WGD on tumorigenesis, highlighting its dualistic influence on cancer cells. We then introduce recent findings on WGD-associated molecular signatures and genetic aberrations and a novel subtype related to WGD. Finally, we discuss the clinical implications of WGD in cancer subtype classification and future therapeutic interventions. Overall, a comprehensive understanding of WGD in cancer biology is crucial to unraveling its complex role in tumorigenesis and identifying novel therapeutic strategies. [BMB Reports 2024; 57(3): 125-134].

Evolutionary trajectories of small cell lung cancer under therapy

The evolutionary processes that underlie the marked sensitivity of small cell lung cancer (SCLC) to chemotherapy and rapid relapse are unknown1-3. Here we determined tumour phylogenies at diagnosis and throughout chemotherapy and immunotherapy by multiregion sequencing of 160 tumours from 65 patients. Treatment-naive SCLC exhibited clonal homogeneity at distinct tumour sites, whereas first-line platinum-based chemotherapy led to a burst in genomic intratumour heterogeneity and spatial clonal diversity. We observed branched evolution and a shift to ancestral clones underlying tumour relapse. Effective radio- or immunotherapy induced a re-expansion of founder clones with acquired genomic damage from first-line chemotherapy. Whereas TP53 and RB1 alterations were exclusively part of the common ancestor, MYC family amplifications were frequently not constituents of the founder clone. At relapse, emerging subclonal mutations affected key genes associated with SCLC biology, and tumours harbouring clonal CREBBP/EP300 alterations underwent genome duplications. Gene-damaging TP53 alterations and co-alterations of TP53 missense mutations with TP73, CREBBP/EP300 or FMN2 were significantly associated with shorter disease relapse following chemotherapy. In summary, we uncover key processes of the genomic evolution of SCLC under therapy, identify the common ancestor as the source of clonal diversity at relapse and show central genomic patterns associated with sensitivity and resistance to chemotherapy.

Detecting Somatic Mutations Without Matched Normal Samples Using Long Reads

DNA sequencing of tumours to identify somatic mutations has become a critical tool to guide the type of treatment given to cancer patients. The gold standard for mutation calling is comparing sequencing data from the tumour to a matched normal sample to avoid mis-classifying inherited SNPs as mutations. This procedure works extremely well, but in certain situations only a tumour sample is available. While approaches have been developed to find mutations without a matched normal, they have limited accuracy or require specific types of input data (e.g. ultra-deep sequencing). Here we explore the application of single molecule long read sequencing to calling somatic mutations without matched normal samples. We develop a simple theoretical framework to show how haplotype phasing is an important source of information for determining whether a variant is a somatic mutation. We then use simulations to assess the range of experimental parameters (tumour purity, sequencing depth) where this approach is effective. These ideas are developed into a prototype somatic mutation caller, smrest, and its use is demonstrated on two highly mutated cancer cell lines. Finally, we argue that this approach has potential to measure clinically important biomarkers that are based on the genome-wide distribution of mutations: tumour mutation burden and mutation signatures.

Pre-existing subclones determine radioresistance in rectal cancer organoids

More than half of all patients with cancer receive radiation therapy, but resistance is commonly observed. Currently, it is unknown whether resistance to radiation therapy is acquired or inherently present. Here, we employed organoids derived from rectal cancer and single-cell whole-genome sequencing to investigate the long-term evolution of subclones in response to radiation. Comparing single-cell whole-genome karyotypes between in-vitro-unirradiated and -irradiated organoids revealed three patterns of subclonal evolution: (1) subclonal persistence, (2) subclonal extinction, and (3) subclonal expansion. Organoids in which subclonal shifts occurred (i.e., expansion or extinction) became more resistant to radiation. Although radioresistant subclones did not share recurrent copy-number alterations that could explain their radioresistance, resistance was associated with reduced chromosomal instability, an association that was also observed in 529 human cancer cell lines. These data suggest that resistance to radiation is inherently present and associated with reduced chromosomal instability.

Converging and evolving immuno-genomic routes toward immune escape in breast cancer

The interactions between tumor and immune cells along the course of breast cancer progression remain largely unknown. Here, we extensively characterize multiple sequential and parallel multiregion tumor and blood specimens of an index patient and a cohort of metastatic triple-negative breast cancers. We demonstrate that a continuous increase in tumor genomic heterogeneity and distinct molecular clocks correlated with resistance to treatment, eventually allowing tumors to escape from immune control. TCR repertoire loses diversity over time, leading to convergent evolution as breast cancer progresses. Although mixed populations of effector memory and cytotoxic single T cells coexist in the peripheral blood, defects in the antigen presentation machinery coupled with subdued T cell recruitment into metastases are observed, indicating a potent immune avoidance microenvironment not compatible with an effective antitumor response in lethal metastatic disease. Our results demonstrate that the immune responses against cancer are not static, but rather follow dynamic processes that match cancer genomic progression, illustrating the complex nature of tumor and immune cell interactions.

Machine learning-based extrachromosomal DNA identification in large-scale cohorts reveals its clinical implications in cancer

The clinical implications of extrachromosomal DNA (ecDNA) in cancer therapy remain largely elusive. Here, we present a comprehensive analysis of ecDNA amplification spectra and their association with clinical and molecular features in multiple cohorts comprising over 13,000 pan-cancer patients. Using our developed computational framework, GCAP, and validating it with multifaceted approaches, we reveal a consistent pan-cancer pattern of mutual exclusivity between ecDNA amplification and microsatellite instability (MSI). In addition, we establish the role of ecDNA amplification as a risk factor and refine genomic subtypes in a cohort from 1015 colorectal cancer patients. Importantly, our investigation incorporates data from four clinical trials focused on anti-PD-1 immunotherapy, demonstrating the pivotal role of ecDNA amplification as a biomarker for guiding checkpoint blockade immunotherapy in gastrointestinal cancer. This finding represents clinical evidence linking ecDNA amplification to the effectiveness of immunotherapeutic interventions. Overall, our study provides a proof-of-concept of identifying ecDNA amplification from cancer whole-exome sequencing (WES) data, highlighting the potential of ecDNA amplification as a valuable biomarker for facilitating personalized cancer treatment.

Genomic and epigenomic basis of breast invasive lobular carcinomas lacking CDH1 genetic alterations

CDH1 (E-cadherin) bi-allelic inactivation is the hallmark alteration of breast invasive lobular carcinoma (ILC), resulting in its discohesive phenotype. A subset of ILCs, however, lack CDH1 genetic/epigenetic inactivation, and their genetic underpinning is unknown. Through clinical targeted sequencing data reanalysis of 364 primary ILCs, we identified 25 ILCs lacking CDH1 bi-allelic genetic alterations. CDH1 promoter methylation was frequent (63%) in these cases. Targeted sequencing reanalysis revealed 3 ILCs harboring AXIN2 deleterious fusions (n = 2) or loss-of-function mutation (n = 1). Whole-genome sequencing of 3 cases lacking bi-allelic CDH1 genetic/epigenetic inactivation confirmed the AXIN2 mutation and no other cell-cell adhesion genetic alterations but revealed a new CTNND1 (p120) deleterious fusion. AXIN2 knock-out in MCF7 cells resulted in lobular-like features, including increased cellular migration and resistance to anoikis. Taken together, ILCs lacking CDH1 genetic/epigenetic alterations are driven by inactivating alterations in other cell adhesion genes (CTNND1 or AXIN2), endorsing a convergent phenotype in ILC.

Allosteric PI3Kα Inhibition Overcomes On-target Resistance to Orthosteric Inhibitors Mediated by Secondary PIK3CA Mutations

PIK3CA mutations occur in ∼8% of cancers, including ∼40% of HR-positive breast cancers, where the PI3K-alpha (PI3Kα)-selective inhibitor alpelisib is FDA approved in combination with fulvestrant. Although prior studies have identified resistance mechanisms, such as PTEN loss, clinically acquired resistance to PI3Kα inhibitors remains poorly understood. Through serial liquid biopsies and rapid autopsies in 39 patients with advanced breast cancer developing acquired resistance to PI3Kα inhibitors, we observe that 50% of patients acquire genomic alterations within the PI3K pathway, including PTEN loss and activating AKT1 mutations. Notably, although secondary PIK3CA mutations were previously reported to increase sensitivity to PI3Kα inhibitors, we identified emergent secondary resistance mutations in PIK3CA that alter the inhibitor binding pocket. Some mutations had differential effects on PI3Kα-selective versus pan-PI3K inhibitors, but resistance induced by all mutations could be overcome by the novel allosteric pan-mutant-selective PI3Kα-inhibitor RLY-2608. Together, these findings provide insights to guide strategies to overcome resistance in PIK3CA-mutated cancers.

SIGNIFICANCE

In one of the largest patient cohorts analyzed to date, this study defines the clinical landscape of acquired resistance to PI3Kα inhibitors. Genomic alterations within the PI3K pathway represent a major mode of resistance and identify a novel class of secondary PIK3CA resistance mutations that can be overcome by an allosteric PI3Kα inhibitor. See related commentary by Gong and Vanhaesebroeck, p. 204 . See related article by Varkaris et al., p. 240 . This article is featured in Selected Articles from This Issue, p. 201.

Rare subclonal sequencing of breast cancers indicates putative metastatic driver mutations are predominately acquired after dissemination

BACKGROUND

Evolutionary models of breast cancer progression differ on the extent to which metastatic potential is pre-encoded within primary tumors. Although metastatic recurrences often harbor putative driver mutations that are not detected in their antecedent primary tumor using standard sequencing technologies, whether these mutations were acquired before or after dissemination remains unclear.

METHODS

To ascertain whether putative metastatic driver mutations initially deemed specific to the metastasis by whole exome sequencing were, in actuality, present within rare ancestral subclones of the primary tumors from which they arose, we employed error-controlled ultra-deep sequencing (UDS-UMI) coupled with FFPE artifact mitigation by uracil-DNA glycosylase (UDG) to assess the presence of 132 "metastasis-specific" mutations within antecedent primary tumors from 21 patients. Maximum mutation detection sensitivity was ~1% of primary tumor cells. A conceptual framework was developed to estimate relative likelihoods of alternative models of mutation acquisition.

RESULTS

The ancestral primary tumor subclone responsible for seeding the metastasis was identified in 29% of patients, implicating several putative drivers in metastatic seeding including LRP5 A65V and PEAK1 K140Q. Despite this, 93% of metastasis-specific mutations in putative metastatic driver genes remained undetected within primary tumors, as did 96% of metastasis-specific mutations in known breast cancer drivers, including ERRB2 V777L, ESR1 D538G, and AKT1 D323H. Strikingly, even in those cases in which the rare ancestral subclone was identified, 87% of metastasis-specific putative driver mutations remained undetected. Modeling indicated that the sequential acquisition of multiple metastasis-specific driver or passenger mutations within the same rare subclonal lineage of the primary tumor was highly improbable.

CONCLUSIONS

Our results strongly suggest that metastatic driver mutations are sequentially acquired and selected within the same clonal lineage both before, but more commonly after, dissemination from the primary tumor, and that these mutations are biologically consequential. Despite inherent limitations in sampling archival primary tumors, our findings indicate that tumor cells in most patients continue to undergo clinically relevant genomic evolution after their dissemination from the primary tumor. This provides further evidence that metastatic recurrence is a multi-step, mutation-driven process that extends beyond primary tumor dissemination and underscores the importance of longitudinal tumor assessment to help guide clinical decisions.

A novel pipeline for prioritizing cancer type-specific therapeutic vulnerabilities using DepMap identifies PAK2 as a target in head and neck squamous cell carcinomas

There is limited guidance on exploiting the genome-wide loss-of-function CRISPR screens in cancer Dependency Map (DepMap) to identify new targets for individual cancer types. This study integrated multiple tools to filter these data in order to seek new therapeutic targets specific to head and neck squamous cell carcinoma (HNSCC). The resulting pipeline prioritized 143 targetable dependencies that represented both well-studied targets and emerging target classes like mitochondrial carriers and RNA-binding proteins. In total, 14 targets had clinical inhibitors used for other cancers or nonmalignant diseases that hold near-term potential to repurpose for HNSCC therapy. Comparing inhibitor response data that were publicly available for 13 prioritized targets between the cell lines with high vs. low dependency on each target uncovered novel therapeutic potential for the PAK2 serine/threonine kinase. PAK2 gene dependency was found to be associated with wild-type p53, low PAK2 mRNA, and diploid status of the 3q amplicon containing PAK2. These findings establish a generalizable pipeline to prioritize clinically relevant targets for individual cancer types using DepMap. Its application to HNSCC highlights novel relevance for PAK2 inhibition and identifies biomarkers of PAK2 inhibitor response.

Nonlobular Invasive Breast Carcinomas with Biallelic Pathogenic CDH1 Somatic Alterations: A Histologic, Immunophenotypic, and Genomic Characterization

CDH1 encodes for E-cadherin, and its loss of function is the hallmark of invasive lobular carcinoma (ILC). Albeit vanishingly rare, biallelic CDH1 alterations may be found in nonlobular breast carcinomas (NL-BCs). We sought to determine the clinicopathologic characteristics and repertoire of genetic alterations of NL-BCs harboring CDH1 biallelic genetic alterations. Analysis of 5842 breast cancers (BCs) subjected to clinical tumor-normal sequencing with an FDA-cleared multigene panel was conducted to identify BCs with biallelic CDH1 pathogenic/likely pathogenic somatic mutations lacking lobular features. The genomic profiles of NL-BCs with CDH1 biallelic genetic alterations were compared with those of ILCs and invasive ductal carcinomas (IDCs), matched by clinicopathologic characteristics. Of the 896 CDH1-altered BCs, 889 samples were excluded based on the diagnosis of invasive mixed ductal/lobular carcinoma or ILC or the detection of monoallelic CDH1 alterations. Only 7 of the 5842 (0.11%) BCs harbored biallelic CDH1 alterations and lacked lobular features. Of these, 4/7 (57%) cases were ER-positive/HER2-negative, 1/7 (14%) was ER-positive/HER2-positive, and 2/7 (29%) were ER-negative/HER2-negative. In total, 5/7 (71%) were of Nottingham grade 2, and 2/7 (29%) were of grade 3. The NL-BCs with CDH1 biallelic genetic alterations included a mucinous carcinoma (n = 1), IDCs with focal nested growth (n = 2), IDC with solid papillary (n = 1) or apocrine (n = 2) features, and an IDC of no special type (NST; n = 1). E-cadherin expression, as detected by immunohistochemistry, was absent (3/5) or aberrant (discontinuous membranous/cytoplasmic/granular; 2/5). However, NL-BCs with CDH1 biallelic genetic alterations displayed recurrent genetic alterations, including TP53, PIK3CA (57%, 4/7; each), FGFR1, and NCOR1 (28%, 2/7, each) alterations. Compared with CDH1 wild-type IDC-NSTs, NL-BCs less frequently harbored GATA3 mutations (0% vs 47%, P = .03), but no significant differences were detected when compared with matched ILCs. Therefore, NL-BCs with CDH1 biallelic genetic alterations are vanishingly rare, predominantly comprise IDCs with special histologic features, and have genomic features akin to luminal B ER-positive BCs.

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.

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.

Investigating USP42 Mutation as Underlying Cause of Familial Non-Medullary Thyroid Carcinoma

In a family with Familial Non-Medullary Thyroid Carcinoma (FNMTC), our investigation using Whole-Exome Sequencing (WES) uncovered a novel germline USP42 mutation [p.(Gly486Arg)]. USP42 is known for regulating p53, cell cycle arrest, and apoptosis, and for being reported as overexpressed in breast and gastric cancer patients. Recently, a USP13 missense mutation was described in FNMTC, suggesting a potential involvement in thyroid cancer. Aiming to explore the USP42 mutation as an underlying cause of FNMTC, our team validated the mutation in blood and tissue samples from the family. Using immunohistochemistry, the expression of USP42, Caspase-3, and p53 was assessed. The USP42 gene was silenced in human thyroid Nthy-Ori 3-1 cells using siRNAs. Subsequently, expression, viability, and morphological assays were conducted. p53, Cyclin D1, p21, and p27 proteins were evaluated by Western blot. USP42 protein was confirmed in all family members and was found to be overexpressed in tumor samples, along with an increased expression of p53 and cleaved Caspase-3. siRNA-mediated USP42 downregulation in Nthy-Ori 3-1 cells resulted in reduced cell viability, morphological changes, and modifications in cell cycle-related proteins. Our results suggest a pivotal role of USP42 mutation in thyroid cell biology, and this finding indicates that USP42 may serve as a new putative target in FNMTC.

labelSeg: segment annotation for tumor copy number alteration profiles

Somatic copy number alterations (SCNAs) are a predominant type of oncogenomic alterations that affect a large proportion of the genome in the majority of cancer samples. Current technologies allow high-throughput measurement of such copy number aberrations, generating results consisting of frequently large sets of SCNA segments. However, the automated annotation and integration of such data are particularly challenging because the measured signals reflect biased, relative copy number ratios. In this study, we introduce labelSeg, an algorithm designed for rapid and accurate annotation of CNA segments, with the aim of enhancing the interpretation of tumor SCNA profiles. Leveraging density-based clustering and exploiting the length-amplitude relationships of SCNA, our algorithm proficiently identifies distinct relative copy number states from individual segment profiles. Its compatibility with most CNA measurement platforms makes it suitable for large-scale integrative data analysis. We confirmed its performance on both simulated and sample-derived data from The Cancer Genome Atlas reference dataset, and we demonstrated its utility in integrating heterogeneous segment profiles from different data sources and measurement platforms. Our comparative and integrative analysis revealed common SCNA patterns in cancer and protein-coding genes with a strong correlation between SCNA and messenger RNA expression, promoting the investigation into the role of SCNA in cancer development.

Temporal change of DNA methylation subclasses between matched newly diagnosed and recurrent glioblastoma

The longitudinal transition of phenotypes is pivotal in glioblastoma treatment resistance and DNA methylation emerged as an important tool for classifying glioblastoma phenotypes. We aimed to characterize DNA methylation subclass heterogeneity during progression and assess its clinical impact. Matched tissues from 47 glioblastoma patients were subjected to DNA methylation profiling, including CpG-site alterations, tissue and serum deconvolution, mass spectrometry, and immunoassay. Effects of clinical characteristics on temporal changes and outcomes were studied. Among 47 patients, 8 (17.0%) had non-matching classifications at recurrence. In the remaining 39 cases, 28.2% showed dominant DNA methylation subclass transitions, with 72.7% being a mesenchymal subclass. In general, glioblastomas with a subclass transition showed upregulated metabolic processes. Newly diagnosed glioblastomas with mesenchymal transition displayed increased stem cell-like states and decreased immune components at diagnosis and exhibited elevated immune signatures and cytokine levels in serum. In contrast, tissue of recurrent glioblastomas with mesenchymal transition showed increased immune components but decreased stem cell-like states. Survival analyses revealed comparable outcomes for patients with and without subclass transitions. This study demonstrates a temporal heterogeneity of DNA methylation subclasses in 28.2% of glioblastomas, not impacting patient survival. Changes in cell state composition associated with subclass transition may be crucial for recurrent glioblastoma targeted therapies.

T proliferating cells derived autophagy signature associated with prognosis and immunotherapy resistance in a pan-cancer analysis

Despite autophagy modulating tumor immunity in the tumor microenvironment (TME), the immunotherapeutic efficacy and potential mechanism of autophagy signature was not explicit. We manually curated an autophagy gene set and defined a pan-cancer autophagy signature by comparing malignant tissues and normal tissues in The Cancer Genome Atlas (TCGA) cohort. The pan-cancer autophagy signature was derived from T proliferating cells as demonstrated in multiple single-cell RNA sequencing (scRNA-seq) datasets. The pan-cancer autophagy signature could influence the cell-cell interactions in the TME and predict the responsiveness of immune checkpoint inhibitors (ICIs) in the metastatic renal cell carcinoma, non-small cell lung cancer, bladder cancer, and melanoma cohorts. Metabolism inactivation accompanied with dysregulation of autophagy was investigated with transcriptomic and proteomic data. The immunotherapeutic predictive role and mechanism regulation of the autophagy signature was validated in an in-house cohort. Our study provides valuable insights into the mechanisms of ICI resistance.