The life history of 21 breast cancers

Genetic landscape of interval and screen detected breast cancer

Interval breast cancers (IBCs) are cancers diagnosed between screening episodes. Understanding the biological differences between IBCs and screen-detected breast-cancers (SDBCs) has the potential to improve mammographic screening and patient management. We analysed and compared the genomic landscape of 288 IBCs and 473 SDBCs by whole genome sequencing of paired tumour-normal patient samples collected as part of the UK 100,000 Genomes Project. Compared to SDBCs, IBCs were more likely to be lobular, higher grade, and triple negative. A more aggressive clinical phenotype was reflected in IBCs displaying features of genomic instability including a higher mutation rate and number of chromosomal structural abnormalities, defective homologous recombination and TP53 mutations. We did not however, find evidence to indicate that IBCs are associated with a significantly different immune response. While IBCs do not represent a unique molecular class of invasive breast cancer they exhibit a more aggressive phenotype, which is likely to be a consequence of the timing of tumour initiation. This information is relevant both with respect to treatment as well as informing the screening interval for mammography.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

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

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

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

A glycolytic metabolite bypasses "two-hit" tumor suppression by BRCA2

Knudson's "two-hit" paradigm posits that carcinogenesis requires inactivation of both copies of an autosomal tumor suppressor gene. Here, we report that the glycolytic metabolite methylglyoxal (MGO) transiently bypasses Knudson's paradigm by inactivating the breast cancer suppressor protein BRCA2 to elicit a cancer-associated, mutational single-base substitution (SBS) signature in nonmalignant mammary cells or patient-derived organoids. Germline monoallelic BRCA2 mutations predispose to these changes. An analogous SBS signature, again without biallelic BRCA2 inactivation, accompanies MGO accumulation and DNA damage in Kras-driven, Brca2-mutant murine pancreatic cancers and human breast cancers. MGO triggers BRCA2 proteolysis, temporarily disabling BRCA2's tumor suppressive functions in DNA repair and replication, causing functional haploinsufficiency. Intermittent MGO exposure incites episodic SBS mutations without permanent BRCA2 inactivation. Thus, a metabolic mechanism wherein MGO-induced BRCA2 haploinsufficiency transiently bypasses Knudson's two-hit requirement could link glycolysis activation by oncogenes, metabolic disorders, or dietary challenges to mutational signatures implicated in cancer evolution.

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

Navigating Breast Cancer Oligometastasis and Oligoprogression: Current Landscape and Future Directions

PURPOSE

We examine the potential for curative approaches among metastatic breast cancer (MBC) patients by exploring the recent literature on local ablative therapies like surgery and stereotactic body radiation therapy (SBRT) in patients with oligometastatic (OM) breast cancer. We also cover therapies for MBC patients with oligoprogressive (OP) disease.

KEY FINDINGS

Surgery and SBRT have been studied for OM and OP breast cancer, mainly in retrospective or non-randomized trials. While many studies demonstrated favorable results, a cooperative study and single-institution trial found no support for surgery/SBRT in OM and OP cases, respectively.

CONCLUSION

While there is interest in applying local therapies to OM and OP breast cancer, the current randomized data does not back the routine use of surgery or SBRT, particularly when considering the potential for treatment-related toxicities. Future research should refine patient selection through advanced imaging and possibly explore these therapies specifically in patients with hormone receptor-positive or HER2-positive disease.

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.

Clonal expansion in normal tissues

Cancer originates from a single ancestral cell that acquires a driver mutation, which confers a growth or survival advantage, followed by the acquisition of additional driver mutations by descendant cells. Recently, it has become evident that somatic cell mutations accumulate in normal tissues with aging and exposure to environmental factors, such as alcohol, smoking, and UV rays, increases the mutation rate. Clones harboring driver mutations expand with age, leading to tissue remodeling. Lineage analysis of myeloproliferative neoplasms and der(1;16)-positive breast cancer revealed that driver mutations were acquired early in our lives and that the development of cancer takes decades, unveiling the previously unknown early process of cancer development. Evidence that clonal hematopoiesis affects various diseases, including nonneoplastic diseases, highlights the potential role of the identification and functional analysis of mutated clones in unraveling unknown pathologies. In this review, we summarize the recent updates on clonal expansion in normal tissues and the natural history of cancer revealed through lineage analysis of noncancerous and cancerous tissues.

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

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

Genomic evolution shapes prostate cancer disease type

The development of cancer is an evolutionary process involving the sequential acquisition of genetic alterations that disrupt normal biological processes, enabling tumor cells to rapidly proliferate and eventually invade and metastasize to other tissues. We investigated the genomic evolution of prostate cancer through the application of three separate classification methods, each designed to investigate a different aspect of tumor evolution. Integrating the results revealed the existence of two distinct types of prostate cancer that arise from divergent evolutionary trajectories, designated as the Canonical and Alternative evolutionary disease types. We therefore propose the evotype model for prostate cancer evolution wherein Alternative-evotype tumors diverge from those of the Canonical-evotype through the stochastic accumulation of genetic alterations associated with disruptions to androgen receptor DNA binding. Our model unifies many previous molecular observations, providing a powerful new framework to investigate prostate cancer disease progression.

Evolving copy number gains promote tumor expansion and bolster mutational diversification

The timing and fitness effect of somatic copy number alterations (SCNA) in cancer evolution remains poorly understood. Here we present a framework to determine the timing of a clonal SCNA that encompasses multiple gains. This involves calculating the proportion of time from its last gain to the onset of population expansion (lead time) as well as the proportion of time prior to its first gain (initiation time). Our method capitalizes on the observation that a genomic segment, while in a specific copy number (CN) state, accumulates point mutations proportionally to its CN. Analyzing 184 whole genome sequenced samples from 75 patients across five tumor types, we commonly observe late gains following early initiating events, occurring just before the clonal expansion relevant to the sampling. These include gains acquired after genome doubling in more than 60% of cases. Notably, mathematical modeling suggests that late clonal gains may contain final-expansion drivers. Lastly, SCNAs bolster mutational diversification between subpopulations, exacerbating the circle of proliferation and increasing heterogeneity.

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.

Deep whole-genome analysis of 494 hepatocellular carcinomas

Over half of hepatocellular carcinoma (HCC) cases diagnosed worldwide are in China1-3. However, whole-genome analysis of hepatitis B virus (HBV)-associated HCC in Chinese individuals is limited4-8, with current analyses of HCC mainly from non-HBV-enriched populations9,10. Here we initiated the Chinese Liver Cancer Atlas (CLCA) project and performed deep whole-genome sequencing (average depth, 120×) of 494 HCC tumours. We identified 6 coding and 28 non-coding previously undescribed driver candidates. Five previously undescribed mutational signatures were found, including aristolochic-acid-associated indel and doublet base signatures, and a single-base-substitution signature that we termed SBS_H8. Pentanucleotide context analysis and experimental validation confirmed that SBS_H8 was distinct to the aristolochic-acid-associated SBS22. Notably, HBV integrations could take the form of extrachromosomal circular DNA, resulting in elevated copy numbers and gene expression. Our high-depth data also enabled us to characterize subclonal clustered alterations, including chromothripsis, chromoplexy and kataegis, suggesting that these catastrophic events could also occur in late stages of hepatocarcinogenesis. Pathway analysis of all classes of alterations further linked non-coding mutations to dysregulation of liver metabolism. Finally, we performed in vitro and in vivo assays to show that fibrinogen alpha chain (FGA), determined as both a candidate coding and non-coding driver, regulates HCC progression and metastasis. Our CLCA study depicts a detailed genomic landscape and evolutionary history of HCC in Chinese individuals, providing important clinical implications.

Prediction of plasma ctDNA fraction and prognostic implications of liquid biopsy in advanced prostate cancer

No consensus strategies exist for prognosticating metastatic castration-resistant prostate cancer (mCRPC). Circulating tumor DNA fraction (ctDNA%) is increasingly reported by commercial and laboratory tests but its utility for risk stratification is unclear. Here, we intersect ctDNA%, treatment outcomes, and clinical characteristics across 738 plasma samples from 491 male mCRPC patients from two randomized multicentre phase II trials and a prospective province-wide blood biobanking program. ctDNA% correlates with serum and radiographic metrics of disease burden and is highest in patients with liver metastases. ctDNA% strongly predicts overall survival, progression-free survival, and treatment response independent of therapeutic context and outperformed established prognostic clinical factors. Recognizing that ctDNA-based biomarker genotyping is limited by low ctDNA% in some patients, we leverage the relationship between clinical prognostic factors and ctDNA% to develop a clinically-interpretable machine-learning tool that predicts whether a patient has sufficient ctDNA% for informative ctDNA genotyping (available online: https://www.ctDNA.org ). Our results affirm ctDNA% as an actionable tool for patient risk stratification and provide a practical framework for optimized biomarker testing.

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.

Breast Tumor Metastasis and Its Microenvironment: It Takes Both Seed and Soil to Grow a Tumor and Target It for Treatment

Metastasis remains a major challenge in treating breast cancer. Breast tumors metastasize to organ-specific locations such as the brain, lungs, and bone, but why some organs are favored over others remains unclear. Breast tumors also show heterogeneity, plasticity, and distinct microenvironments. This contributes to treatment failure and relapse. The interaction of breast cancer cells with their metastatic microenvironment has led to the concept that primary breast cancer cells act as seeds, whereas the metastatic tissue microenvironment (TME) is the soil. Improving our understanding of this interaction could lead to better treatment strategies for metastatic breast cancer. Targeted treatments for different subtypes of breast cancers have improved overall patient survival, even with metastasis. However, these targeted treatments are based upon the biology of the primary tumor and often these patients' relapse, after therapy, with metastatic tumors. The advent of immunotherapy allowed the immune system to target metastatic tumors. Unfortunately, immunotherapy has not been as effective in metastatic breast cancer relative to other cancers with metastases, such as melanoma. This review will describe the heterogeneic nature of breast cancer cells and their microenvironments. The distinct properties of metastatic breast cancer cells and their microenvironments that allow interactions, especially in bone and brain metastasis, will also be described. Finally, we will review immunotherapy approaches to treat metastatic breast tumors and discuss future therapeutic approaches to improve treatments for metastatic breast cancer.

Ferroptosis-Based Therapeutic Strategies toward Precision Medicine for Cancer

Ferroptosis is a type of iron-dependent programmed cell death characterized by the dysregulation of iron metabolism and the accumulation of lipid peroxides. This nonapoptotic mode of cell death is implicated in various physiological and pathological processes. Recent findings have underscored its potential as an innovative strategy for cancer treatment, particularly against recalcitrant malignancies that are resistant to conventional therapies. This article focuses on ferroptosis-based therapeutic strategies for precision cancer treatment, covering the molecular mechanisms of ferroptosis, four major types of ferroptosis inducers and their inhibitory effects on diverse carcinomas, the detection of ferroptosis by fluorescent probes, and their implementation in image-guided therapy. These state-of-the-art tactics have manifested enhanced selectivity and efficacy against malignant carcinomas. Given that the administration of ferroptosis in cancer therapy is still at a burgeoning stage, some major challenges and future perspectives are discussed for the clinical translation of ferroptosis into precision cancer treatment.

Intra-prostatic tumour evolution, steps in metastatic spread and histogenomic associations revealed by integration of multi-region whole-genome sequencing with histopathological features

BACKGROUND

Extension of prostate cancer beyond the primary site by local invasion or nodal metastasis is associated with poor prognosis. Despite significant research on tumour evolution in prostate cancer metastasis, the emergence and evolution of cancer clones at this early stage of expansion and spread are poorly understood. We aimed to delineate the routes of evolution and cancer spread within the prostate and to seminal vesicles and lymph nodes, linking these to histological features that are used in diagnostic risk stratification.

METHODS

We performed whole-genome sequencing on 42 prostate cancer samples from the prostate, seminal vesicles and lymph nodes of five treatment-naive patients with locally advanced disease. We spatially mapped the clonal composition of cancer across the prostate and the routes of spread of cancer cells within the prostate and to seminal vesicles and lymph nodes in each individual by analysing a total of > 19,000 copy number corrected single nucleotide variants.

RESULTS

In each patient, we identified sample locations corresponding to the earliest part of the malignancy. In patient 10, we mapped the spread of cancer from the apex of the prostate to the seminal vesicles and identified specific genomic changes associated with the transformation of adenocarcinoma to amphicrine morphology during this spread. Furthermore, we show that the lymph node metastases in this patient arose from specific cancer clones found at the base of the prostate and the seminal vesicles. In patient 15, we observed increased mutational burden, altered mutational signatures and histological changes associated with whole genome duplication. In all patients in whom histological heterogeneity was observed (4/5), we found that the distinct morphologies were located on separate branches of their respective evolutionary trees.

CONCLUSIONS

Our results link histological transformation with specific genomic alterations and phylogenetic branching. These findings have implications for diagnosis and risk stratification, in addition to providing a rationale for further studies to characterise the genetic changes causally linked to morphological transformation. Our study demonstrates the value of integrating multi-region sequencing with histopathological data to understand tumour evolution and identify mechanisms of prostate cancer spread.

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.

Quantifying the seed sensitivity of cancer subclonal reconstruction algorithms

Background

Intra-tumoural heterogeneity complicates cancer prognosis and impairs treatment success. One of the ways subclonal reconstruction (SRC) quantifies intra-tumoural heterogeneity is by estimating the number of subclones present in bulk DNA sequencing data. SRC algorithms are probabilistic and need to be initialized by a random seed. However, the seeds used in bioinformatics algorithms are rarely reported in the literature. Thus, the impact of the initializing seed on SRC solutions has not been studied. To address this gap, we generated a set of ten random seeds to systematically benchmark the seed sensitivity of three probabilistic SRC algorithms: PyClone-VI, DPClust, and PhyloWGS.

Results

We characterized the seed sensitivity of three algorithms across fourteen whole-genome sequences of head and neck squamous cell carcinoma and nine SRC pipelines, each composed of a single nucleotide variant caller, a copy number aberration caller and an SRC algorithm. This led to a total of 1470 subclonal reconstructions, including 1260 single-region and 210 multi-region reconstructions. The number of subclones estimated per patient vary across SRC pipelines, but all three SRC algorithms show substantial seed sensitivity: subclone estimates vary across different seeds for the same set of input using the same SRC algorithm. No seed consistently estimated the mode number of subclones across all patients for any SRC algorithm.

Conclusions

These findings highlight the variability in quantifying intra-tumoural heterogeneity introduced by the seed sensitivity of probabilistic SRC algorithms. We recommend that authors, reviewers and editors adopt guidelines to both report and randomize seed choices. It may also be valuable to consider seed-sensitivity in the benchmarking of newly developed SRC algorithms. These findings may be of interest in other areas of bioinformatics where seeded probabilistic algorithms are used and suggest consideration of formal seed reporting standards to enhance reproducibility.

Spatial architectures of somatic mutations in normal prostate, benign prostatic hyperplasia and coexisting prostate cancer

This study aimed to identify somatic mutations in nontumor cells (NSMs) in normal prostate and benign prostatic hyperplasia (BPH) and to determine their relatedness to prostate cancer (PCA). From 22 PCA patients, two prostates were sampled for 3-dimensional mapping (50 normal, 46 BPH and 1 PCA samples), and 20 prostates were trio-sampled (two normal or BPH samples and one PCA sample) and analyzed by whole-genome sequencing. Normal and BPH tissues harbored several driver NSMs and copy number alterations (CNAs), including in FOXA1, but the variations exhibited low incidence, rare recurrence, and rare overlap with PCAs. CNAs, structural variants, and mutation signatures were similar between normal and BPH samples, while BPHs harbored a higher mutation burden, shorter telomere length, larger clone size, and more private NSMs than normal prostates. We identified peripheral-zonal dominance and right-side asymmetry in NSMs, but the asymmetry was heterogeneous between samples. In one normal prostate, private oncogenic RAS-signaling NSMs were detected, suggesting convergence in clonal maintenance. Early embryonic mutations exhibited two distinct distributions, characterized as layered and mixed patterns. Our study identified that the BPH genome differed from the normal prostate genome but was still closer to the normal genome than to the PCA genome, suggesting that BPH might be more related to aging or environmental stress than to tumorigenic processes.

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.

The molecular map of CLL and Richter's syndrome

Clonal expansion of B-cells, from the early stages of monoclonal B-cell lymphocytosis through to chronic lymphocytic leukemia (CLL), and then in some cases to Richter's syndrome (RS) provides a comprehensive model of cancer evolution, notable for the marked morphological transformation and distinct clinical phenotypes. High-throughput sequencing of large cohorts of patients and single-cell studies have generated a molecular map of CLL and more recently, of RS, yielding fundamental insights into these diseases and of clonal evolution. A selection of CLL driver genes have been functionally interrogated to yield novel insights into the biology of CLL. Such findings have the potential to impact patient care through risk stratification, treatment selection and drug discovery. However, this molecular map remains incomplete, with extant questions concerning the origin of the B-cell clone, the role of the TME, inter- and intra-compartmental heterogeneity and of therapeutic resistance mechanisms. Through the application of multi-modal single-cell technologies across tissues, disease states and clinical contexts, these questions can now be addressed with the answers holding great promise of generating translatable knowledge to improve patient care.

CONIPHER: a computational framework for scalable phylogenetic reconstruction with error correction

Intratumor heterogeneity provides the fuel for the evolution and selection of subclonal tumor cell populations. However, accurate inference of tumor subclonal architecture and reconstruction of tumor evolutionary histories from bulk DNA sequencing data remains challenging. Frequently, sequencing and alignment artifacts are not fully filtered out from cancer somatic mutations, and errors in the identification of copy number alterations or complex evolutionary events (e.g., mutation losses) affect the estimated cellular prevalence of mutations. Together, such errors propagate into the analysis of mutation clustering and phylogenetic reconstruction. In this Protocol, we present a new computational framework, CONIPHER (COrrecting Noise In PHylogenetic Evaluation and Reconstruction), that accurately infers subclonal structure and phylogenetic relationships from multisample tumor sequencing, accounting for both copy number alterations and mutation errors. CONIPHER has been used to reconstruct subclonal architecture and tumor phylogeny from multisample tumors with high-depth whole-exome sequencing from the TRACERx421 dataset, as well as matched primary-metastatic cases. CONIPHER outperforms similar methods on simulated datasets, and in particular scales to a large number of tumor samples and clones, while completing in under 1.5 h on average. CONIPHER enables automated phylogenetic analysis that can be effectively applied to large sequencing datasets generated with different technologies. CONIPHER can be run with a basic knowledge of bioinformatics and R and bash scripting languages.

Reconstructing cell lineage trees with genomic barcoding: approaches and applications

In multicellular organisms, developmental history of cell divisions and functional annotation of terminal cells can be organized into a cell lineage tree (CLT). The reconstruction of the CLT has long been a major goal in developmental biology and other related fields. Recent technological advancements, especially those in editable genomic barcodes and single-cell high-throughput sequencing, have sparked a new wave of experimental methods for reconstructing CLTs. Here we review the existing experimental approaches to the reconstruction of CLT, which are broadly categorized as either image-based or DNA barcode-based methods. In addition, we present a summary of the related literature based on the biological insight provided by the obtained CLTs. Moreover, we discuss the challenges that will arise as more and better CLT data become available in the near future. Genomic barcoding-based CLT reconstructions and analyses, due to their wide applicability and high scalability, offer the potential for novel biological discoveries, especially those related to general and systemic properties of the developmental process.

Genomic and transcriptomic characteristics of 12 novel primary cell lines derived from three patients with cholangiocarcinoma

Cholangiocarcinoma (CCA) is an aggressive bile duct malignancy with poor prognosis. To improve our understanding of the biological characteristics of CCA and develop effective therapies, appropriate preclinical models are required. Here, we established and characterized 12 novel patient-derived primary cancer cell (PDPC) models using multi-region sampling. At the genomic level of PDPCs, we observed not only commonly mutated genes, such as TP53, JAK3, and KMT2C, consistent with the reports in CCA, but also specific mutation patterns in each cell line. In addition, specific expression patterns with distinct biological functions and pathways involved were also observed in the PDPCs at the transcriptomic level. Furthermore, the drug-sensitivity results revealed that the PDPCs exhibited different responses to the six commonly used compounds. Our findings indicate that the established PDPCs can serve as novel in vitro reliable models to provide a crucial molecular basis for improving the understanding of tumorigenesis and its treatment.

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.

Druggable growth dependencies and tumor evolution analysis in patient-derived organoids of neuroendocrine neoplasms from multiple body sites

Neuroendocrine neoplasms (NENs) comprise well-differentiated neuroendocrine tumors (NETs) and poorly differentiated neuroendocrine carcinomas (NECs). Treatment options for patients with NENs are limited, in part due to lack of accurate models. We establish patient-derived tumor organoids (PDTOs) from pulmonary NETs and derive PDTOs from an understudied subtype of NEC, large cell neuroendocrine carcinoma (LCNEC), arising from multiple body sites. PDTOs maintain the gene expression patterns, intra-tumoral heterogeneity, and evolutionary processes of parental tumors. Through hypothesis-driven drug sensitivity analyses, we identify ASCL1 as a potential biomarker for response of LCNEC to treatment with BCL-2 inhibitors. Additionally, we discover a dependency on EGF in pulmonary NET PDTOs. Consistent with these findings, we find that, in an independent cohort, approximately 50% of pulmonary NETs express EGFR. This study identifies an actionable vulnerability for a subset of pulmonary NETs, emphasizing the utility of these PDTO models.

Whole genome sequencing refines stratification and therapy of patients with clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer, but a comprehensive description of its genomic landscape is lacking. We report the whole genome sequencing of 778 ccRCC patients enrolled in the 100,000 Genomes Project, providing the most detailed somatic mutational landscape to date. We identify new driver genes, which as well as emphasising the major role of epigenetic regulation in ccRCC highlight additional biological pathways extending opportunities for drug repurposing. Genomic characterisation identified patients with divergent clinical outcome; higher number of structural copy number alterations associated with poorer prognosis, whereas VHL mutations were independently associated with a better prognosis. The twin observations that higher T-cell infiltration is associated with better outcome and that genetically predicted immune evasion is not common supports the rationale for immunotherapy. These findings should inform personalised surveillance and treatment strategies for ccRCC patients.

Completing a genomic characterisation of microscopic tumour samples with copy number

BACKGROUND

Genomic insights in settings where tumour sample sizes are limited to just hundreds or even tens of cells hold great clinical potential, but also present significant technical challenges. We previously developed the DigiPico sequencing platform to accurately identify somatic mutations from such samples.

RESULTS

Here, we complete this genomic characterisation with copy number. We present a novel protocol, PicoCNV, to call allele-specific somatic copy number alterations from picogram quantities of tumour DNA. We find that PicoCNV provides exactly accurate copy number in 84% of the genome for even the smallest samples, and demonstrate its clinical potential in maintenance therapy.

CONCLUSIONS

PicoCNV complements our existing platform, allowing for accurate and comprehensive genomic characterisations of cancers in settings where only microscopic samples are available.

Circulating tumor DNA as liquid biopsy in lung cancer: Biological characteristics and clinical integration

Lung cancer maintains high morbidity and mortality rate globally despite significant advancements in diagnosis and treatment in the era of precision medicine. Pathological analysis of tumor tissue, the current gold standard for lung cancer diagnosis, is intrusive and intrinsically confined to evaluating the limited amount of tissues that could be physically extracted. However, tissue biopsy has several limitations, including the invasiveness of the procedure and difficulty in obtaining samples for patients at advanced stages., there Additionally,has been no major breakthrough in tumor biomarkers with high specificity and sensitivity, particularly for early-stage lung cancer. Liquid biopsy has been considered a feasible auxiliary tool for tearly dianosis, evaluating treatment responses and monitoring prognosis of lung cancer. Circulating tumor DNA (ctDNA), an ideal biomarker of liquid biopsy, has emerged as one of the most reliable tools for monitoring tumor processes at molecular levels. Herein, this review focuses on tumor heterogeneity to elucidate the superiority of liquid biopsy and retrospectively discussdeciphersolution. We systematically elaborate ctDNA biological characteristics, introduce methods for ctDNA detection, and discuss the current role of plasma ctDNA in lung cancer management. Finally, we summarize the drawbacks of ctDNA analysis and highlight its potential clinical application in lung cancer.

Phenotypic noise and plasticity in cancer evolution

Non-genetic alterations can produce changes in a cell's phenotype. In cancer, these phenomena can influence a cell's fitness by conferring access to heritable, beneficial phenotypes. Herein, we argue that current discussions of 'phenotypic plasticity' in cancer evolution ignore a salient feature of the original definition: namely, that it occurs in response to an environmental change. We suggest 'phenotypic noise' be used to distinguish non-genetic changes in phenotype that occur independently from the environment. We discuss the conceptual and methodological techniques used to identify these phenomena during cancer evolution. We propose that the distinction will guide efforts to define mechanisms of phenotype change, accelerate translational work to manipulate phenotypes through treatment, and, ultimately, improve patient outcomes.

Multi-scale characterisation of homologous recombination deficiency in breast cancer

BACKGROUND

Homologous recombination is a robust, broadly error-free mechanism of double-strand break repair, and deficiencies lead to PARP inhibitor sensitivity. Patients displaying homologous recombination deficiency can be identified using 'mutational signatures'. However, these patterns are difficult to reliably infer from exome sequencing. Additionally, as mutational signatures are a historical record of mutagenic processes, this limits their utility in describing the current status of a tumour.

METHODS

We apply two methods for characterising homologous recombination deficiency in breast cancer to explore the features and heterogeneity associated with this phenotype. We develop a likelihood-based method which leverages small insertions and deletions for high-confidence classification of homologous recombination deficiency for exome-sequenced breast cancers. We then use multinomial elastic net regression modelling to develop a transcriptional signature of heterogeneous homologous recombination deficiency. This signature is then applied to single-cell RNA-sequenced breast cancer cohorts enabling analysis of homologous recombination deficiency heterogeneity and differential patterns of tumour microenvironment interactivity.

RESULTS

We demonstrate that the inclusion of indel events, even at low levels, improves homologous recombination deficiency classification. Whilst BRCA-positive homologous recombination deficient samples display strong similarities to those harbouring BRCA1/2 defects, they appear to deviate in microenvironmental features such as hypoxic signalling. We then present a 228-gene transcriptional signature which simultaneously characterises homologous recombination deficiency and BRCA1/2-defect status, and is associated with PARP inhibitor response. Finally, we show that this signature is applicable to single-cell transcriptomics data and predict that these cells present a distinct milieu of interactions with their microenvironment compared to their homologous recombination proficient counterparts, typified by a decreased cancer cell response to TNFα signalling.

CONCLUSIONS

We apply multi-scale approaches to characterise homologous recombination deficiency in breast cancer through the development of mutational and transcriptional signatures. We demonstrate how indels can improve homologous recombination deficiency classification in exome-sequenced breast cancers. Additionally, we demonstrate the heterogeneity of homologous recombination deficiency, especially in relation to BRCA1/2-defect status, and show that indications of this feature can be captured at a single-cell level, enabling further investigations into interactions between DNA repair deficient cells and their tumour microenvironment.

Evolutionary determinants of curability in cancer

The emergence of drug-resistant cells, most of which have a mutated TP53 gene, prevents curative treatment in most advanced and common metastatic cancers of adults. Yet, a few, rarer malignancies, all of which are TP53 wild type, have high cure rates. In this Perspective, we discuss how common features of curable cancers offer insights into the evolutionary and developmental determinants of drug resistance. Acquired loss of TP53 protein function is the most common genetic change in cancer. This probably reflects positive selection in the context of strong ecosystem pressures including microenvironmental hypoxia. Loss of TP53's functions results in multiple fitness benefits and enhanced evolvability of cancer cells. TP53-null cells survive apoptosis, and tolerate potent oncogenic signalling, DNA damage and genetic instability. In addition, critically, they provide an expanded pool of self-renewing, or stem, cells, the primary units of evolutionary selection in cancer, making subsequent adaptation to therapeutic challenge by drug resistance highly probable. The exceptional malignancies that are curable, including the common genetic subtype of childhood acute lymphoblastic leukaemia and testicular seminoma, differ from the common adult cancers in originating prenatally from embryonic or fetal cells that are developmentally primed for TP53-dependent apoptosis. Plus, they have other genetic and phenotypic features that enable dissemination without exposure to selective pressures for TP53 loss, retaining their intrinsic drug hypersensitivity.

Clonal tracking in cancer and metastasis

The eradication of many cancers has proven challenging due to the presence of functionally and genetically heterogeneous clones maintained by rare cancer stem cells (CSCs), which contribute to disease progression, treatment refractoriness, and late relapse. The characterization of functional CSC activity has necessitated the development of modern clonal tracking strategies. This review describes viral-based and CRISPR-Cas9-based cellular barcoding, lineage tracing, and imaging-based approaches. DNA-based cellular barcoding technology is emerging as a powerful and robust strategy that has been widely applied to in vitro and in vivo model systems, including patient-derived xenograft models. This review also highlights the potential of these methods for use in the clinical and drug discovery contexts and discusses the important insights gained from such approaches.

A pan-tissue survey of mosaic chromosomal alterations in 948 individuals

Genetic mutations accumulate in an organism's body throughout its lifetime. While somatic single-nucleotide variants have been well characterized in the human body, the patterns and consequences of large chromosomal alterations in normal tissues remain largely unknown. Here, we present a pan-tissue survey of mosaic chromosomal alterations (mCAs) in 948 healthy individuals from the Genotype-Tissue Expression project, augmenting RNA-based allelic imbalance estimation with haplotype phasing. We found that approximately a quarter of the individuals carry a clonally-expanded mCA in at least one tissue, with incidence strongly correlated with age. The prevalence and genome-wide patterns of mCAs vary considerably across tissue types, suggesting tissue-specific mutagenic exposure and selection pressures. The mCA landscapes in normal adrenal and pituitary glands resemble those in tumors arising from these tissues, whereas the same is not true for the esophagus and skin. Together, our findings show a widespread age-dependent emergence of mCAs across normal human tissues with intricate connections to tumorigenesis.

Human APOBEC3B promotes tumor development in vivo including signature mutations and metastases

The antiviral DNA cytosine deaminase APOBEC3B has been implicated as a source of mutation in many cancers. However, despite years of work, a causal relationship has yet to be established in vivo. Here, we report a murine model that expresses tumor-like levels of human APOBEC3B. Animals expressing full-body APOBEC3B appear to develop normally. However, adult males manifest infertility, and older animals of both sexes show accelerated rates of carcinogenesis, visual and molecular tumor heterogeneity, and metastasis. Both primary and metastatic tumors exhibit increased frequencies of C-to-T mutations in TC dinucleotide motifs consistent with the established biochemical activity of APOBEC3B. Enrichment for APOBEC3B-attributable single base substitution mutations also associates with elevated levels of insertion-deletion mutations and structural variations. APOBEC3B catalytic activity is required for all of these phenotypes. Together, these studies provide a cause-and-effect demonstration that human APOBEC3B is capable of driving both tumor initiation and evolution in vivo.

Cancer origin tracing and timing in two high-risk prostate cancers using multisample whole genome analysis: prospects for personalized medicine

BACKGROUND

Prostate cancer (PrCa) genomic heterogeneity causes resistance to therapies such as androgen deprivation. Such heterogeneity can be deciphered in the context of evolutionary principles, but current clinical trials do not include evolution as an essential feature. Whether or not analysis of genomic data in an evolutionary context in primary prostate cancer can provide unique added value in the research and clinical domains remains an open question.

METHODS

We used novel processing techniques to obtain whole genome data together with 3D anatomic and histomorphologic analysis in two men (GP5 and GP12) with high-risk PrCa undergoing radical prostatectomy. A total of 22 whole genome-sequenced sites (16 primary cancer foci and 6 lymph node metastatic) were analyzed using evolutionary reconstruction tools and spatio-evolutionary models. Probability models were used to trace spatial and chronological origins of the primary tumor and metastases, chart their genetic drivers, and distinguish metastatic and non-metastatic subclones.

RESULTS

In patient GP5, CDK12 inactivation was among the first mutations, leading to a PrCa tandem duplicator phenotype and initiating the cancer around age 50, followed by rapid cancer evolution after age 57, and metastasis around age 59, 5 years prior to prostatectomy. In patient GP12, accelerated cancer progression was detected after age 54, and metastasis occurred around age 56, 3 years prior to prostatectomy. Multiple metastasis-originating events were identified in each patient and tracked anatomically. Metastasis from prostate to lymph nodes occurred strictly ipsilaterally in all 12 detected events. In this pilot, metastatic subclone content analysis appears to substantially enhance the identification of key drivers. Evolutionary analysis' potential impact on therapy selection appears positive in these pilot cases.

CONCLUSIONS

PrCa evolutionary analysis allows tracking of anatomic site of origin, timing of cancer origin and spread, and distinction of metastatic-capable from non-metastatic subclones. This enables better identification of actionable targets for therapy. If extended to larger cohorts, it appears likely that similar analyses could add substantial biological insight and clinically relevant value.

Mutational topography reflects clinical neuroblastoma heterogeneity

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although clinical risk-defining genomic alterations exist in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography and mutational signatures derived from all variant classes, we identified co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is associated with differences in the mutational processes driving these scenarios, linking risk-defining pathognomonic variants to distinct molecular processes. Whereas high-risk MYCN-amplified neuroblastomas were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas were marked by footprints of chromosome mis-segregation and TOP1 mutational activity. Furthermore, analysis of subclonal mutations uncovered differential activity of these processes through neuroblastoma evolution. Thus, clinical heterogeneity of neuroblastoma patients can be linked to differences in the mutational processes that are active in their tumors.

Characterization of Cancer Evolution Landscape Based on Accurate Detection of Somatic Mutations in Single Tumor Cells

Accurate detection of somatic mutations in single tumor cells is greatly desired as it allows us to quantify the single-cell mutation burden and construct the mutation-based phylogenetic tree. Here we developed scNanoSeq chemistry and profiled 842 single cells from 21 human breast cancer samples. The majority of the mutation-based phylogenetic trees comprise a characteristic stem evolution followed by the clonal sweep. We observed the subtype-dependent lengths in the stem evolution. To explain this phenomenon, we propose that the differences are related to different reprogramming required for different subtypes of breast cancer. Furthermore, we reason that the time that the tumor-initiating cell took to acquire the critical clonal-sweep-initiating mutation by random chance set the time limit for the reprogramming process. We refer to this model as a reprogramming and critical mutation co-timing (RCMC) subtype model. Next, in the sweeping clone, we observed that tumor cells undergo a branched evolution with rapidly decreasing selection. In the most recent clades, effectively neutral evolution has been reached, resulting in a substantially large number of mutational heterogeneities. Integrative analysis with 522-713X ultra-deep bulk whole genome sequencing (WGS) further validated this evolution mode. Mutation-based phylogenetic trees also allow us to identify the early branched cells in a few samples, whose phylogenetic trees support the gradual evolution of copy number variations (CNVs). Overall, the development of scNanoSeq allows us to unveil novel insights into breast cancer evolution.

Acute myeloid leukemia: from NGS, through scRNA-seq, to CAR-T. dissect cancer heterogeneity and tailor the treatment

Acute myeloid leukemia (AML) is a malignant blood cancer with marked cellular heterogeneity due to altered maturation and differentiation of myeloid blasts, the possible causes of which are transcriptional or epigenetic alterations, impaired apoptosis, and excessive cell proliferation. This neoplasm has a high rate of resistance to anticancer therapies and thus a high risk of relapse and mortality because of both the biological diversity of the patient and intratumoral heterogeneity due to the acquisition of new somatic changes. For more than 40 years, the old gold standard "one size fits all" treatment approach included intensive chemotherapy treatment with anthracyclines and cytarabine.The manuscript first traces the evolution of the understanding of the pathology from the 1970s to the present. The enormous strides made in its categorization prove to be crucial for risk stratification, enabling an increasingly personalized diagnosis and treatment approach.Subsequently, we highlight how, over the past 15 years, technological advances enabling single cell RNA sequencing and T-cell modification based on the genomic tools are affecting the classification and treatment of AML. At the dawn of the new millennium, the advent of high-throughput next-generation sequencing technologies has enabled the profiling of patients evidencing different facets of the same disease, stratifying risk, and identifying new possible therapeutic targets that have subsequently been validated. Currently, the possibility of investigating tumor heterogeneity at the single cell level, profiling the tumor at the time of diagnosis or after treatments exist. This would allow the identification of underrepresented cellular subclones or clones resistant to therapeutic approaches and thus responsible for post-treatment relapse that would otherwise be difficult to detect with bulk investigations on the tumor biopsy. Single-cell investigation will then allow even greater personalization of therapy to the genetic and transcriptional profile of the tumor, saving valuable time and dangerous side effects. The era of personalized medicine will take a huge step forward through the disclosure of each individual piece of the complex puzzle that is cancer pathology, to implement a "tailored" therapeutic approach based also on engineered CAR-T cells.

Genomic signatures of past and present chromosomal instability in Barrett's esophagus and early esophageal adenocarcinoma

The progression of precancerous lesions to malignancy is often accompanied by increasing complexity of chromosomal alterations but how these alterations arise is poorly understood. Here we perform haplotype-specific analysis of chromosomal copy-number evolution in the progression of Barrett's esophagus (BE) to esophageal adenocarcinoma (EAC) on multiregional whole-genome sequencing data of BE with dysplasia and microscopic EAC foci. We identify distinct patterns of copy-number evolution indicating multigenerational chromosomal instability that is initiated by cell division errors but propagated only after p53 loss. While abnormal mitosis, including whole-genome duplication, underlies chromosomal copy-number changes, segmental alterations display signatures of successive breakage-fusion-bridge cycles and chromothripsis of unstable dicentric chromosomes. Our analysis elucidates how multigenerational chromosomal instability generates copy-number variation in BE cells, precipitates complex alterations including DNA amplifications, and promotes their independent clonal expansion and transformation. In particular, we suggest sloping copy-number variation as a signature of ongoing chromosomal instability that precedes copy-number complexity. These findings suggest copy-number heterogeneity in advanced cancers originates from chromosomal instability in precancerous cells and such instability may be identified from the presence of sloping copy-number variation in bulk sequencing data.

Non-homogeneous Poisson and renewal processes as spatial models for cancer mutation

Advances in sequencing technology assisted biologists in revealing signatures of DNA cancer mutation process and in demonstrating the mutagenesis behind. However, most of these signatures proposed by majority of work focus only on the type and frequency of mutations, without considering spatial information which is non-negligible in exploring mechanisms of mutation occurrence, e.g., Kataegis. Statistical characterization of the distance between consecutive mutations can give us relative spatial information; however, it ignores location information which is as important as distance information. In this work, we assume that DNA cancer mutations are location-dependent and that integrating the two variables, location and inter-distance, is beneficial to study DNA cancer mutation processes more accurately. Particularly, instead of following a specific distribution over the whole DNA sequence, we found out that the distribution of distance between successive mutations alternates between exponential and power-law distributions. Apart from this, the parameters of either of the two distributions vary with DNA locations. The cancers with kataegis phenomenon, a specific mutation pattern caused by abnormal activity of APOBEC protein family, are more likely to be accompanied by higher parameter values of distance distribution, implying higher occurrence rate of mutation. Therefore, we propose non-homogeneous Poisson and Renewal processes to spatially model DNA cancer mutations and to describe mutation patterns quantitatively and more accurately through a statistical perspective.

Refphase: Multi-sample phasing reveals haplotype-specific copy number heterogeneity

Most computational methods that infer somatic copy number alterations (SCNAs) from bulk sequencing of DNA analyse tumour samples individually. However, the sequencing of multiple tumour samples from a patient's disease is an increasingly common practice. We introduce Refphase, an algorithm that leverages this multi-sampling approach to infer haplotype-specific copy numbers through multi-sample phasing. We demonstrate Refphase's ability to infer haplotype-specific SCNAs and characterise their intra-tumour heterogeneity, to uncover previously undetected allelic imbalance in low purity samples, and to identify parallel evolution in the context of whole genome doubling in a pan-cancer cohort of 336 samples from 99 tumours.

Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole-genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole-genome sequencing (WGS) but still generate hundreds of artifacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolbox (PTATO), to accurately detect single base substitutions, insertions-deletions (indels), and structural variants in PTA-based WGS data. PTATO includes a machine learning approach and filtering based on recurrence to distinguish PTA artifacts from true mutations with high sensitivity (up to 90%), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that hematopoietic stem cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS, have normal somatic single base substitution burdens but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.

Genetic inhibitors of APOBEC3B-induced mutagenesis

The cytidine deaminases APOBEC3A (A3A) and APOBEC3B (A3B) are prominent mutators of human cancer genomes. However, tumor-specific genetic modulators of APOBEC-induced mutagenesis are poorly defined. Here, we used a screen to identify 61 gene deletions that increase A3B-induced mutations in yeast. We also determined whether each deletion was epistatic with Ung1 loss, which indicated whether the encoded factors participate in the homologous recombination (HR)-dependent bypass of A3B/Ung1-dependent abasic sites or suppress A3B-catalyzed deamination by protecting against aberrant formation of single-stranded DNA (ssDNA). We found that the mutation spectra of A3B-induced mutations revealed genotype-specific patterns of strand-specific ssDNA formation and nucleotide incorporation across APOBEC-induced lesions. Combining these three metrics, we were able to establish a multifactorial signature of APOBEC-induced mutations specific to (1) failure to remove H3K56 acetylation, (2) defective CTF18-RFC complex function, and (3) defective HR-mediated bypass of APOBEC-induced lesions. We extended these results by analyzing mutation data for human tumors and found BRCA1/2-deficient breast cancers display three- to fourfold more APOBEC-induced mutations. Mirroring our results in yeast, Rev1-mediated C-to-G substitutions are mainly responsible for increased APOBEC-signature mutations in BRCA1/2-deficient tumors, and these mutations associate with lagging strand synthesis during replication. These results identify important factors that influence DNA replication dynamics and likely the abundance of APOBEC-induced mutation during tumor progression. They also highlight a novel role for BRCA1/2 during HR-dependent lesion bypass of APOBEC-induced lesions during cancer cell replication.

Validation of liquid biopsy for ESR1-mutation analysis in hormone-sensitive breast cancer: a pooled meta-analysis

Several retrospective and prospective studies have shown that genomic alterations in Estrogen-receptor one (ESR1) can be characterized not only in tissue samples but also by sequencing circulating tumor DNA (ctDNA) in liquid biopsy. Therefore, liquid biopsy is a potential noninvasive surrogate for tissue biopsy. This meta-analysis was designed to compare the prevalence of ESR 1 mutation detected with liquid biopsy and tissue biopsy. A pooled meta-analysis of studies published between 1 January 2007 and 1 March 2021 was conducted regarding the methodologies used for ESR1 mutation analysis. Liquid biopsy is a valid, inexpensive, and attractive noninvasive alternative to tumor biopsies for the identification of ESR1 mutations. Liquid biopsy for ESR 1 analysis would facilitate regular testing, allowing monitoring of the sensitivity to ET and guiding treatment strategies.

Early clonal extinction in glioblastoma progression revealed by genetic barcoding

Glioblastoma progression in its early stages remains poorly understood. Here, we transfer PDGFB and genetic barcodes in mouse brain to initiate gliomagenesis and enable direct tracing of glioblastoma evolution from its earliest possible stage. Unexpectedly, we observe a high incidence of clonal extinction events and progressive divergence in clonal sizes, even after the acquisition of malignant phenotype. Computational modeling suggests these dynamics result from clonal-based cell-cell competition. Through bulk and single-cell transcriptome analyses, coupled with lineage tracing, we reveal that Myc transcriptional targets have the strongest correlation with clonal size imbalances. Moreover, we show that the downregulation of Myc expression is sufficient to drive competitive dynamics in intracranially transplanted gliomas. Our findings provide insights into glioblastoma evolution that are inaccessible using conventional retrospective approaches, highlighting the potential of combining clonal tracing and transcriptomic analyses in this field.