Punctuated copy number evolution and clonal stasis in triple-negative breast cancer

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.

Portable platform for leukocyte extraction from blood using sheath-free microfluidic DLD

Leukocyte count is routinely performed for diagnostic purposes and is rapidly emerging as a significant biomarker for a wide array of diseases. Additionally, leukocytes have demonstrated considerable promise in novel cell-based immunotherapies. However, the direct retrieval of leukocytes from whole blood is a significant challenge due to their low abundance compared to erythrocytes. Here, we introduce a microfluidic-based platform that isolates and recovers leukocytes from diluted whole blood in a single step. Our platform utilizes a novel, sheathless method to initially sediment and focus blood cells into a dense stream while flowing through a tubing before entering the microfluidic device. A hexagonal-shaped structure, patterned at the device's inlet, directs all the blood cells against the channel's outer walls. The focused cells are then separated based on their size using the deterministic lateral displacement (DLD) microfluidic technique. We evaluated various parameters that could influence leukocyte separation, including different focusing structures (assessed both computationally and experimentally), the orientation of the tubing-chip interface, the effects of blood sample hematocrit (dilution), and flow rate. Our device demonstrated the ability to isolate leukocytes from diluted blood with a separation efficiency of 100%, a recovery rate of 76%, and a purity of 80%, while maintaining a cell viability of 98%. The device operates for over 30 min at a flow rate of 2 μL min-1. Furthermore, we developed a handheld pressure controller to drive fluid flow, enhancing the operability of our platform outside of central laboratories and enabling near-patient testing. Our platform can be integrated with downstream cell-based assays and analytical methods that require high leukocyte purity (80%), ranging from cell counting to diagnostics and cell culture applications.

NestedBD: Bayesian inference of phylogenetic trees from single-cell copy number profiles under a birth-death model

Copy number aberrations (CNAs) are ubiquitous in many types of cancer. Inferring CNAs from cancer genomic data could help shed light on the initiation, progression, and potential treatment of cancer. While such data have traditionally been available via "bulk sequencing," the more recently introduced techniques for single-cell DNA sequencing (scDNAseq) provide the type of data that makes CNA inference possible at the single-cell resolution. We introduce a new birth-death evolutionary model of CNAs and a Bayesian method, NestedBD, for the inference of evolutionary trees (topologies and branch lengths with relative mutation rates) from single-cell data. We evaluated NestedBD's performance using simulated data sets, benchmarking its accuracy against traditional phylogenetic tools as well as state-of-the-art methods. The results show that NestedBD infers more accurate topologies and branch lengths, and that the birth-death model can improve the accuracy of copy number estimation. And when applied to biological data sets, NestedBD infers plausible evolutionary histories of two colorectal cancer samples. NestedBD is available at https://github.com/Androstane/NestedBD .

High clonal diversity and spatial genetic admixture in early prostate cancer and surrounding normal tissue

Somatic copy number alterations (SCNAs) are pervasive in advanced human cancers, but their prevalence and spatial distribution in early-stage, localized tumors and their surrounding normal tissues are poorly characterized. Here, we perform multi-region, single-cell DNA sequencing to characterize the SCNA landscape across tumor-rich and normal tissue in two male patients with localized prostate cancer. We identify two distinct karyotypes: 'pseudo-diploid' cells harboring few SCNAs and highly aneuploid cells. Pseudo-diploid cells form numerous small-sized subclones ranging from highly spatially localized to broadly spread subclones. In contrast, aneuploid cells do not form subclones and are detected throughout the prostate, including normal tissue regions. Highly localized pseudo-diploid subclones are confined within tumor-rich regions and carry deletions in multiple tumor-suppressor genes. Our study reveals that SCNAs are widespread in normal and tumor regions across the prostate in localized prostate cancer patients and suggests that a subset of pseudo-diploid cells drive tumorigenesis in the aging prostate.

Machine-learning analysis reveals an important role for negative selection in shaping cancer aneuploidy landscapes

BACKGROUND

Aneuploidy, an abnormal number of chromosomes within a cell, is a hallmark of cancer. Patterns of aneuploidy differ across cancers, yet are similar in cancers affecting closely related tissues. The selection pressures underlying aneuploidy patterns are not fully understood, hindering our understanding of cancer development and progression.

RESULTS

Here, we apply interpretable machine learning methods to study tissue-selective aneuploidy patterns. We define 20 types of features corresponding to genomic attributes of chromosome-arms, normal tissues, primary tumors, and cancer cell lines (CCLs), and use them to model gains and losses of chromosome arms in 24 cancer types. To reveal the factors that shape the tissue-specific cancer aneuploidy landscapes, we interpret the machine learning models by estimating the relative contribution of each feature to the models. While confirming known drivers of positive selection, our quantitative analysis highlights the importance of negative selection for shaping aneuploidy landscapes. This is exemplified by tumor suppressor gene density being a better predictor of gain patterns than oncogene density, and vice versa for loss patterns. We also identify the importance of tissue-selective features and demonstrate them experimentally, revealing KLF5 as an important driver for chr13q gain in colon cancer. Further supporting an important role for negative selection in shaping the aneuploidy landscapes, we find compensation by paralogs to be among the top predictors of chromosome arm loss prevalence and demonstrate this relationship for one paralog interaction. Similar factors shape aneuploidy patterns in human CCLs, demonstrating their relevance for aneuploidy research.

CONCLUSIONS

Our quantitative, interpretable machine learning models improve the understanding of the genomic properties that shape cancer aneuploidy landscapes.

Tumor phylogeography reveals block-shaped spatial heterogeneity and the mode of evolution in Hepatocellular Carcinoma

Solid tumors are complex ecosystems with heterogeneous 3D structures, but the spatial intra-tumor heterogeneity (sITH) at the macroscopic (i.e., whole tumor) level is under-explored. Using a phylogeographic approach, we sequence genomes and transcriptomes from 235 spatially informed sectors across 13 hepatocellular carcinomas (HCC), generating one of the largest datasets for studying sITH. We find that tumor heterogeneity in HCC segregates into spatially variegated blocks with large genotypic and phenotypic differences. By dissecting the transcriptomic heterogeneity, we discover that 30% of patients had a "spatially competing distribution" (SCD), where different spatial blocks have distinct transcriptomic subtypes co-existing within a tumor, capturing the critical transition period in disease progression. Interestingly, the tumor regions with more advanced transcriptomic subtypes (e.g., higher cell cycle) often take clonal dominance with a wider geographic range, rejecting neutral evolution for SCD patients. Extending the statistical tests for detecting natural selection to many non-SCD patients reveal varying levels of selective signal across different tumors, implying that many evolutionary forces including natural selection and geographic isolation can influence the overall pattern of sITH. Taken together, tumor phylogeography unravels a dynamic landscape of sITH, pinpointing important evolutionary and clinical consequences of spatial heterogeneity in cancer.

Punctuational evolution is pervasive in distal site metastatic colonization

The evolution of metastasis represents a lethal stage of cancer progression. Yet, the evolutionary kinetics of metastatic disease remain unresolved. Here, using single cell CRISPR-Cas9 lineage tracing data, we show that in metastatic disease, gradual molecular evolution is punctuated by episodes of rapid evolutionary change associated with lineage divergence. By measuring punctuational effects across the metastatic cascade, we show that punctuational effects contribute more to the molecular diversity at distal site metastases compared to the paired primary tumor, suggesting qualitatively different modes of evolution may drive primary and metastatic tumor progression. This is the first empirical evidence for distinct patterns of molecular evolution at early and late stages of metastasis and demonstrates the complex interplay of cell intrinsic and extrinsic factors that shape lethal cancer.

Exploiting the metabolic vulnerability of circulating tumour cells

Metastasis has a major part in the severity of disease and lethality of cancer. Circulating tumour cells (CTCs) represent a reservoir of metastatic precursors in circulation, most of which cannot survive due to hostile conditions in the bloodstream. Surviving cells colonise a secondary site based on a combination of physical, metabolic, and oxidative stress protection states required for that environment. Recent advances in CTC isolation methods and high-resolution 'omics technologies are revealing specific metabolic pathways that support this selection of CTCs. In this review, we discuss recent advances in our understanding of CTC biology and discoveries of adaptations in metabolic pathways during their selection. Understanding these traits and delineating mechanisms by which they confer acquired resistance or vulnerability in CTCs is crucial for developing successful prognostic and therapeutic strategies in cancer.

The reckoning of chromosomal instability: past, present, future

Quantitative measures of CIN are crucial to our understanding of its role in cancer. Technological advances have changed the way CIN is quantified, offering increased accuracy and insight. Here, we review measures of CIN through its rise as a field, discuss considerations for its measurement, and look forward to future quantification of CIN.

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

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

SIGNIFICANCE

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

A statistical learning method for simultaneous copy number estimation and subclone clustering with single-cell sequencing data

The availability of single-cell sequencing (SCS) enables us to assess intra-tumor heterogeneity and identify cellular subclones without the confounding effect of mixed cells. Copy number aberrations (CNAs) have been commonly used to identify subclones in SCS data using various clustering methods, as cells comprising a subpopulation are found to share a genetic profile. However, currently available methods may generate spurious results (e.g., falsely identified variants) in the procedure of CNA detection, thereby diminishing the accuracy of subclone identification within a large, complex cell population. In this study, we developed a subclone clustering method based on a fused lasso model, referred to as FLCNA, which can simultaneously detect CNAs in single-cell DNA sequencing (scDNA-seq) data. Spike-in simulations were conducted to evaluate the clustering and CNA detection performance of FLCNA, benchmarking it against existing copy number estimation methods (SCOPE, HMMcopy) in combination with commonly used clustering methods. Application of FLCNA to a scDNA-seq data set of breast cancer revealed different genomic variation patterns in neoadjuvant chemotherapy-treated samples and pretreated samples. We show that FLCNA is a practical and powerful method for subclone identification and CNA detection with scDNA-seq data.

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.

A Countable-Type Branching Process Model for the Tug-of-War Cancer Cell Dynamics

We consider a time-continuous Markov branching process of proliferating cells with a countable collection of types. Among-type transitions are inspired by the Tug-of-War process introduced by McFarland et al. (Proc Natl Acad Sci 111(42):15138-15143, 2014) as a mathematical model for competition of advantageous driver mutations and deleterious passenger mutations in cancer cells. We introduce a version of the model in which a driver mutation pushes the type of the cell L-units up, while a passenger mutation pulls it 1-unit down. The distribution of time to divisions depends on the type (fitness) of cell, which is an integer. The extinction probability given any initial cell type is strictly less than 1, which allows us to investigate the transition between types (type transition) in an infinitely long cell lineage of cells. The analysis leads to the result that under driver dominance, the type transition process escapes to infinity, while under passenger dominance, it leads to a limit distribution. Implications in cancer cell dynamics and population genetics are discussed.

SILAC proteomics based on 3D cell spheroids unveils the role of RAC2 in regulating the crosstalk between triple-negative breast cancer cells and tumor-associated macrophages

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer, and is characterized by a high infiltration of tumor-associated macrophages (TAMs). TAMs contribute significantly to tumor progression by intricately interacting with tumor cells. Deeply investigating the interaction between TNBC cells and TAMs is of great importance for finding potential biomarkers and developing novel therapeutic strategies to further improve the clinical outcomes of TNBC patients. In this study, we confirmed the interplay using both 3D and 2D co-culture models. The stable-isotype labeling by amino acids in cell culture (SILAC)-based quantitative proteomics was conducted on 3D cell spheroids containing TNBC cells and macrophages to identify the potential candidate in regulating the crosstalk between TNBC and TAMs. Ras-related C3 botulinum toxin substrate 2 (RAC2) was identified as a potential molecule for further exploration, given its high expression in TNBC and positive correlation with M2 macrophage infiltration. The suppression of RAC2 inhibited TNBC cell proliferation, migration, invasion, and epithelial-mesenchymal transition (EMT) in vitro. Meanwhile, knocking down RAC2 in TNBC cells impaired macrophage recruitment and M2 polarization. Mechanistically, RAC2 exerted its roles in TNBC cells and TAMs by regulating the activation of P65 NF-κB and P38 MAPK, while TAMs further elevated RAC2 expression and P65 NF-κB activation by secreting soluble mediators including IL-10. These findings highlight the significance of RAC2 as a crucial molecule in the crosstalk between TNBC and TAMs, suggesting it could be a promising therapeutic target in TNBC.

Circulating tumor DNA-based copy-number profiles enable monitoring treatment effects during therapy in high-grade serous carcinoma

Circulating tumor DNA (ctDNA) analysis has emerged as a promising tool for detecting and profiling longitudinal genomics changes in cancer. While copy-number alterations (CNAs) play a major role in cancers, treatment effect monitoring using copy-number profiles has received limited attention as compared to mutations. A major reason for this is the insensitivity of CNA analysis for the real-life tumor-fraction ctDNA samples. We performed copy-number analysis on 152 plasma samples obtained from 29 patients with high-grade serous ovarian cancer (HGSC) using a sequencing panel targeting over 500 genes. Twenty-one patients had temporally matched tissue and plasma sample pairs, which enabled assessing concordance with tissues sequenced with the same panel or whole-genome sequencing and to evaluate sensitivity. Our approach could detect concordant CNA profiles in most plasma samples with as low as 5% tumor content and highly amplified regions in samples with ∼1% of tumor content. Longitudinal profiles showed changes in the CNA profiles in seven out of 11 patients with high tumor-content plasma samples at relapse. These changes included focal acquired or lost copy-numbers, even though most of the genome remained stable. Two patients displayed major copy-number profile changes during therapy. Our analysis revealed ctDNA-detectable subclonal selection resulting from both surgical operations and chemotherapy. Overall, longitudinal ctDNA data showed acquired and diminished CNAs at relapse when compared to pre-treatment samples. These results highlight the importance of genomic profiling during treatment as well as underline the usability of ctDNA.

Somatic Copy Number Alterations in Colorectal Cancer Lead to a Differentially Expressed ceRNA Network (ceRNet)

Colorectal cancer (CRC) represents the second deadliest malignancy worldwide. Around 75% of CRC patients exhibit high levels of chromosome instability that result in the accumulation of somatic copy number alterations. These alterations are associated with the amplification of oncogenes and deletion of tumor-ppressor genes and contribute to the tumoral phenotype in different malignancies. Even though this relationship is well known, much remains to be investigated regarding the effect of said alterations in long non-coding RNAs (lncRNAs) and, in turn, the impact these alterations have on the tumor phenotype. The present study aimed to evaluate the role of differentially expressed lncRNAs coded in regions with copy number alterations in colorectal cancer patient samples. We downloaded RNA-seq files of the Colorectal Adenocarcinoma Project from the The Cancer Genome Atlas (TCGA) repository (285 sequenced tumor tissues and 41 non-tumor tissues), evaluated differential expression, and mapped them over genome sequencing data with regions presenting copy number alterations. We obtained 78 differentially expressed (LFC > 1|< -1, padj < 0.05) lncRNAs, 410 miRNAs, and 5028 mRNAs and constructed a competing endogenous RNA (ceRNA) network, predicting significant lncRNA-miRNA-mRNA interactions. Said network consisted of 30 lncRNAs, 19 miRNAs, and 77 mRNAs. To understand the role that our ceRNA network played, we performed KEGG and GO analysis and found several oncogenic and anti-oncogenic processes enriched by the molecular players in our network. Finally, to evaluate the clinical relevance of the lncRNA expression, we performed survival analysis and found that C5orf64, HOTAIR, and RRN3P3 correlated with overall patient survival. Our results showed that lncRNAs coded in regions affected by SCNAs form a complex gene regulatory network in CCR.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.

SIGNIFICANCE

The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Application of single-cell sequencing to the research of tumor microenvironment

Single-cell sequencing is a technique for detecting and analyzing genomes, transcriptomes, and epigenomes at the single-cell level, which can detect cellular heterogeneity lost in conventional sequencing hybrid samples, and it has revolutionized our understanding of the genetic heterogeneity and complexity of tumor progression. Moreover, the tumor microenvironment (TME) plays a crucial role in the formation, development and response to treatment of tumors. The application of single-cell sequencing has ushered in a new age for the TME analysis, revealing not only the blueprint of the pan-cancer immune microenvironment, but also the heterogeneity and differentiation routes of immune cells, as well as predicting tumor prognosis. Thus, the combination of single-cell sequencing and the TME analysis provides a unique opportunity to unravel the molecular mechanisms underlying tumor development and progression. In this review, we summarize the recent advances in single-cell sequencing and the TME analysis, highlighting their potential applications in cancer research and clinical translation.

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.

Growth Dynamics of Ductal Carcinoma in Situ Recapitulate Normal Breast Development

Ductal carcinoma in situ (DCIS) and invasive breast cancer share many morphologic, proteomic, and genomic alterations. Yet in contrast to invasive cancer, many DCIS tumors do not progress and may remain indolent over decades. To better understand the heterogenous nature of this disease, we reconstructed the growth dynamics of 18 DCIS tumors based on the geo-spatial distribution of their somatic mutations. The somatic mutation topographies revealed that DCIS is multiclonal and consists of spatially discontinuous subclonal lesions. Here we show that this pattern of spread is consistent with a new 'Comet' model of DCIS tumorigenesis, whereby multiple subclones arise early and nucleate the buds of the growing tumor. The discontinuous, multiclonal growth of the Comet model is analogous to the branching morphogenesis of normal breast development that governs the rapid expansion of the mammary epithelium during puberty. The branching morphogenesis-like dynamics of the proposed Comet model diverges from the canonical model of clonal evolution, and better explains observed genomic spatial data. Importantly, the Comet model allows for the clinically relevant scenario of extensive DCIS spread, without being subjected to the selective pressures of subclone competition that promote the emergence of increasingly invasive phenotypes. As such, the normal cell movement inferred during DCIS growth provides a new explanation for the limited risk of progression in DCIS and adds biologic rationale for ongoing clinical efforts to reduce DCIS overtreatment.

Cracking the pattern of tumor evolution based on single-cell copy number alterations

Copy number alterations (CNAs) are a key characteristic of tumor development and progression. The accumulation of various CNAs during tumor development plays a critical role in driving tumor evolution. Heterogeneous clones driven by distinct CNAs have different selective advantages, leading to differential patterns of tumor evolution that are essential for developing effective cancer therapies. Recent advances in single-cell sequencing technology have enabled genome-wide copy number profiling of tumor cell populations at single-cell resolution. This has made it possible to explore the evolutionary patterns of CNAs and accurately discover the mechanisms of intra-tumor heterogeneity. Here, we propose a two-step statistical approach that distinguishes neutral, linear, branching and punctuated evolutionary patterns for a tumor cell population based on single-cell copy number profiles. We assessed our approach using a variety of simulated and real single-cell genomic and transcriptomic datasets, demonstrating its high accuracy and robustness in predicting tumor evolutionary patterns. We applied our approach to single-cell DNA sequencing data from 20 breast cancer patients and observed that punctuated evolution is the dominant evolutionary pattern in breast cancer. Similar conclusions were drawn when applying the approach to single-cell RNA sequencing data obtained from 132 various cancer patients. Moreover, we found that differential immune cell infiltration is associated with specific evolutionary patterns. The source code of our study is available at https://github.com/FangWang-SYSU/PTEM.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Cancer-inspired genomics mapper model for the generation of synthetic DNA sequences with desired genomics signatures

Genome data are crucial in modern medicine, offering significant potential for diagnosis and treatment. Thanks to technological advancements, many millions of healthy and diseased genomes have already been sequenced; however, obtaining the most suitable data for a specific study, and specifically for validation studies, remains challenging with respect to scale and access. Therefore, in silico genomics sequence generators have been proposed as a possible solution. However, the current generators produce inferior data using mostly shallow (stochastic) connections, detected with limited computational complexity in the training data. This means they do not take the appropriate biological relations and constraints, that originally caused the observed connections, into consideration. To address this issue, we propose cancer-inspired genomics mapper model (CGMM), that combines genetic algorithm (GA) and deep learning (DL) methods to tackle this challenge. CGMM mimics processes that generate genetic variations and mutations to transform readily available control genomes into genomes with the desired phenotypes. We demonstrate that CGMM can generate synthetic genomes of selected phenotypes such as ancestry and cancer that are indistinguishable from real genomes of such phenotypes, based on unsupervised clustering. Our results show that CGMM outperforms four current state-of-the-art genomics generators on two different tasks, suggesting that CGMM will be suitable for a wide range of purposes in genomic medicine, especially for much-needed validation studies.

Archival single-cell genomics reveals persistent subclones during DCIS progression

Ductal carcinoma in situ (DCIS) is a common precursor of invasive breast cancer. Our understanding of its genomic progression to recurrent disease remains poor, partly due to challenges associated with the genomic profiling of formalin-fixed paraffin-embedded (FFPE) materials. Here, we developed Arc-well, a high-throughput single-cell DNA-sequencing method that is compatible with FFPE materials. We validated our method by profiling 40,330 single cells from cell lines, a frozen tissue, and 27 FFPE samples from breast, lung, and prostate tumors stored for 3-31 years. Analysis of 10 patients with matched DCIS and cancers that recurred 2-16 years later show that many primary DCIS had already undergone whole-genome doubling and clonal diversification and that they shared genomic lineages with persistent subclones in the recurrences. Evolutionary analysis suggests that most DCIS cases in our cohort underwent an evolutionary bottleneck, and further identified chromosome aberrations in the persistent subclones that were associated with recurrence.

Characteristics of DNA macro-alterations in breast cancer with liver metastasis before treatment

BACKGROUND

Whole-genome doubling (WGD) has been observed in 30% of cancers, followed by a highly complex rearranged karyotype unfavourable to breast cancer's outcome. However, the macro-alterations that characterise liver metastasis in breast cancer(BC) are poorly understood. Here, we conducted a whole-genome sequencing analysis of liver metastases to explore the status and the time frame model of these macro-alterations in pre-treatment patients with metastatic breast cancer.

RESULTS

Whole-genome sequencing was conducted in 11 paired primary tumours, lymph node metastasis, and liver metastasis fresh samples from four patients with late-stage breast cancer. We also chose five postoperative frozen specimens from patients with early-stage breast cancer before any treatment as control. Surprisingly, all four liver metastasis samples were classified as WGD + . However, the previous study reported that WGD happened in 30% of cancers and 2/5 in our early-stage samples. WGD was not observed in the two separate primary tumours and one lymph node metastasis of one patient with metastatic BC, but her liver metastasis showed an early burst of bi-allelic copy number gain. The phylogenetic tree proves her 4 tumour samples were the polyclonal origin and only one WGD + clone metastasis to the liver. Another 3 metastatic BC patients' primary tumour and lymph node metastasis experienced WGD as well as liver metastasis, and they all showed similar molecular time-frame of copy number(CN) gain across locations within the same patient. These patients' tumours were of monoclonal origin, and WGD happened in a founding clone before metastasis, explaining that all samples share the CN-gain time frame. After WGD, the genomes usually face instability to evolve other macro-alterations. For example, a greater quantity and variety of complex structural variations (SVs) were detected in WGD + samples. The breakpoints were enriched in the chr17: 39 Mb-40 Mb tile, which contained the HER2 gene, resulting in the formation of tyfonas, breakage-fusion-bridge cycles, and double minutes. These complex SVs may be involved in the evolutionary mechanisms of the dramatic increase of HER2 copy number.

CONCLUSION

Our work revealed that the WGD + clone might be a critical evolution step for liver metastasis and favoured following complex SV of breast cancer.

Interferon signaling promotes tolerance to chromosomal instability during metastatic evolution in renal cancer

Molecular routes to metastatic dissemination are critical determinants of aggressive cancers. Through in vivo CRISPR-Cas9 genome editing, we generated somatic mosaic genetically engineered models that faithfully recapitulate metastatic renal tumors. Disruption of 9p21 locus is an evolutionary driver to systemic disease through the rapid acquisition of complex karyotypes in cancer cells. Cross-species analysis revealed that recurrent patterns of copy number variations, including 21q loss and dysregulation of the interferon pathway, are major drivers of metastatic potential. In vitro and in vivo genomic engineering, leveraging loss-of-function studies, along with a model of partial trisomy of chromosome 21q, demonstrated a dosage-dependent effect of the interferon receptor genes cluster as an adaptive mechanism to deleterious chromosomal instability in metastatic progression. This work provides critical knowledge on drivers of renal cell carcinoma progression and defines the primary role of interferon signaling in constraining the propagation of aneuploid clones in cancer evolution.

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

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

CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples

Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq unravels T cell-related prognostic risk model and tumor immune microenvironment modulation in triple-negative breast cancer

BACKGROUND

Triple negative breast cancer (TNBC) is an aggressive and fatal malignancy. The current success of tumor immunotherapy has focused attention on intermediate T-cell subsets and the tumor microenvironment, which are essential for activation of the anti-tumor response. Therefore, both areas require further research to accelerate progress in developing tailored immunotherapeutic approaches for patients with TNBC.

METHODS

We obtained scRNA-seq data of TNBC from the GEO database. A multiplex strategy was used to analyze and identify the T-cell heterogeneity of TNBC. By combining the METABRIC and GEO databases, a prognostic risk model for T-cell marker genes was constructed and validated. In addition, the immune-infiltrating cells of TNBC was analyzed using CIBERSORT, and the association between the risk model and response to immunotherapy was investigated.

RESULTS

Based on scRNA-seq data, 25,932 cells were identified for multiple analyzes. T cells were studied with a focus on 2 subtypes, including CD8⁺ and CD4⁺. There were also communication relationships between T cells and multiple cell types. The results of the enrichment analysis showed that the T-cell marker genes were focused in pathways related to the immune system. In addition, OPTN, TMEM176A, PKM and HES1 deserve attention as prognostic markers in TNBC. The immune infiltration results showed that the high-risk group had significant immune cell infiltration and immunosuppression status.

CONCLUSION

This study provides a resource for understanding T-cell heterogeneity and the associated prognostic risk model for TNBC. The results show that the model helps predict prognosis and response to treatment in breast cancer.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

A statistical learning method for simultaneous copy number estimation and subclone clustering with single cell sequencing data

The availability of single cell sequencing (SCS) enables us to assess intra-tumor heterogeneity and identify cellular subclones without the confounding effect of mixed cells. Copy number aberrations (CNAs) have been commonly used to identify subclones in SCS data using various clustering methods, since cells comprising a subpopulation are found to share genetic profile. However, currently available methods may generate spurious results (e.g., falsely identified CNAs) in the procedure of CNA detection, hence diminishing the accuracy of subclone identification from a large complex cell population. In this study, we developed a CNA detection method based on a fused lasso model, referred to as FLCNA, which can simultaneously identify subclones in single cell DNA sequencing (scDNA-seq) data. Spike-in simulations were conducted to evaluate the clustering and CNA detection performance of FLCNA benchmarking to existing copy number estimation methods (SCOPE, HMMcopy) in combination with the existing and commonly used clustering methods. Interestingly, application of FLCNA to a real scDNA-seq dataset of breast cancer revealed remarkably different genomic variation patterns in neoadjuvant chemotherapy treated samples and pre-treated samples. We show that FLCNA is a practical and powerful method in subclone identification and CNA detection with scDNA-seq data.

Clonal origin and development of high hyperdiploidy in childhood acute lymphoblastic leukaemia

High hyperdiploid acute lymphoblastic leukemia (HeH ALL), one of the most common childhood malignancies, is driven by nonrandom aneuploidy (abnormal chromosome numbers) mainly comprising chromosomal gains. In this study, we investigate how aneuploidy in HeH ALL arises. Single cell whole genome sequencing of 2847 cells from nine primary cases and one normal bone marrow reveals that HeH ALL generally display low chromosomal heterogeneity, indicating that they are not characterized by chromosomal instability and showing that aneuploidy-driven malignancies are not necessarily chromosomally heterogeneous. Furthermore, most chromosomal gains are present in all leukemic cells, suggesting that they arose early during leukemogenesis. Copy number data from 577 primary cases reveals selective pressures that were used for in silico modeling of aneuploidy development. This shows that the aneuploidy in HeH ALL likely arises by an initial tripolar mitosis in a diploid cell followed by clonal evolution, in line with a punctuated evolution model.

Somatic variant detection from multi-sampled genomic sequencing data of tumor specimens using the ith.Variant pipeline

A common technique for uncovering intra-tumor genomic heterogeneity (ITH) is variant detection. However, it can be challenging to reliably characterize ITH given uneven sample quality (e.g., depth of coverage, tumor purity, and subclonality). We describe a protocol for calling point mutations and copy number alterations using sequencing of multiple related clinical patient samples across diverse tissue, optimizing for sensitivity with specificity. The ith.Variant pipeline can be run on single- or multi-region whole-genome and whole-exome sequencing. For complete details on the use and execution of this protocol, please refer to Sun et al. (2017).1.

Deciphering breast cancer: from biology to the clinic

Breast cancer remains a leading cause of cancer-related mortality in women, reflecting profound disease heterogeneity, metastasis, and therapeutic resistance. Over the last decade, genomic and transcriptomic data have been integrated on an unprecedented scale and revealed distinct cancer subtypes, critical molecular drivers, clonal evolutionary trajectories, and prognostic signatures. Furthermore, multi-dimensional integration of high-resolution single-cell and spatial technologies has highlighted the importance of the entire breast cancer ecosystem and the presence of distinct cellular "neighborhoods." Clinically, a plethora of new targeted therapies has emerged, now being rapidly incorporated into routine care. Resistance to therapy, however, remains a crucial challenge for the field.

A review on genetic alterations in CNS metastases related to breast cancer treatment. Is there a role for liquid biopsies in CSF?

Acquired mutations or altered gene expression patterns in brain metastases (BM) and/or leptomeningeal metastases (LM) of breast cancer may play a role in therapy-resistance and offer new molecular targets and treatment options. Despite expanding knowledge of genetic alterations in breast cancer and their metastases, clinical applications for patients with central nervous system (CNS) metastases are currently limited. An emerging tool are DNA-techniques that may detect genetic alterations of the CNS metastases in the cerebrospinal fluid (CSF). In this review we discuss genetic studies in breast cancer and CNS metastases and the role of liquid biopsies in CSF.

Isogenic Cell Lines Derived from Specific Organ Metastases Exhibit Divergent Cytogenomic Aberrations

Aneuploidy, a deviation in chromosome numbers from the normal diploid set, is now recognized as a fundamental characteristic of all cancer types and is found in 70-90% of all solid tumors. The majority of aneuploidies are generated by chromosomal instability (CIN). CIN/aneuploidy is an independent prognostic marker of cancer survival and is a cause of drug resistance. Hence, ongoing research has been directed towards the development of therapeutics aimed at targeting CIN/aneuploidy. However, there are relatively limited reports on the evolution of CIN/aneuploidies within or across metastatic lesions. In this work, we built on our previous studies using a human xenograft model system of metastatic disease in mice that is based on isogenic cell lines derived from the primary tumor and specific metastatic organs (brain, liver, lung, and spine). As such, these studies were aimed at exploring distinctions and commonalities between the karyotypes; biological processes that have been implicated in CIN; single-nucleotide polymorphisms (SNPs); losses, gains, and amplifications of chromosomal regions; and gene mutation variants across these cell lines. Substantial amounts of inter- and intra-heterogeneity were found across karyotypes, along with distinctions between SNP frequencies across each chromosome of each metastatic cell line relative the primary tumor cell line. There were disconnects between chromosomal gains or amplifications and protein levels of the genes in those regions. However, commonalities across all cell lines provide opportunities to select biological processes as druggable targets that could have efficacy against the primary tumor, as well as metastases.

Spatiotemporal view of malignant histogenesis and macroevolution via formation of polyploid giant cancer cells

To understand how malignant tumors develop, we tracked cell membrane, nuclear membrane, spindle, and cell cycle dynamics in polyploid giant cancer cells (PGCCs) during the formation of high-grade serous carcinoma organoids using long-term time-lapse imaging. Single cells underwent traditional mitosis to generate tissue with uniform nuclear size, while others formed PGCCs via asymmetric mitosis, endoreplication, multipolar endomitosis, nuclear fusion, and karyokinesis without cytokinesis. PGCCs underwent restitution multipolar endomitosis, nuclear fragmentation, and micronuclei formation to increase nuclear contents and heterogeneity. At the cellular level, the development of PGCCs was associated with forming transient intracellular cells, termed fecundity cells. The fecundity cells can be decellularized to facilitate nuclear fusion and synchronized with other nuclei for subsequent nuclear replication. PGCCs can undergo several rounds of entosis to form complex tissue structures, termed fecundity structures. The formation of PGCCs via multiple modes of nuclear replication in the absence of cytokinesis leads to an increase in the nuclear-to-cytoplasmic (N/C) ratio and intracellular cell reproduction, which is remarkably similar to the mode of nuclear division during pre-embryogenesis. Our data support that PGCCs may represent a central regulator in malignant histogenesis, intratumoral heterogeneity, immune escape, and macroevolution via the de-repression of suppressed pre-embryogenic program in somatic cells.

Reconstructing mutational lineages in breast cancer by multi-patient-targeted single-cell DNA sequencing

Single-cell DNA sequencing (scDNA-seq) methods are powerful tools for profiling mutations in cancer cells; however, most genomic regions sequenced in single cells are non-informative. To overcome this issue, we developed a multi-patient-targeted (MPT) scDNA-seq method. MPT involves first performing bulk exome sequencing across a cohort of cancer patients to identify somatic mutations, which are then pooled together to develop a single custom targeted panel for high-throughput scDNA-seq using a microfluidics platform. We applied MPT to profile 330 mutations across 23,500 cells from 5 patients with triple negative-breast cancer (TNBC), which showed that 3 tumors were monoclonal and 2 tumors were polyclonal. From these data, we reconstructed mutational lineages and identified early mutational and copy-number events, including early TP53 mutations that occurred in all five patients. Collectively, our data suggest that MPT can overcome a major technical obstacle for studying tumor evolution using scDNA-seq by profiling information-rich mutation sites.

Breast cancer heterogeneity and its implication in personalized precision therapy

Breast cancer heterogeneity determines cancer progression, treatment effects, and prognosis. However, the precise mechanism for this heterogeneity remains unknown owing to its complexity. Here, we summarize the origins of breast cancer heterogeneity and its influence on disease progression, recurrence, and therapeutic resistance. We review the possible mechanisms of heterogeneity and the research methods used to analyze it. We also highlight the importance of cell interactions for the origins of breast cancer heterogeneity, which can be further categorized into cooperative and competitive interactions. Finally, we provide new insights into precise individual treatments based on heterogeneity.

scONE-seq: A single-cell multi-omics method enables simultaneous dissection of phenotype and genotype heterogeneity from frozen tumors

Single-cell multi-omics can provide a unique perspective on tumor cellular heterogeneity. Most previous single-cell whole-genome RNA sequencing (scWGS-RNA-seq) methods demonstrate utility with intact cells from fresh samples. Among them, many are not applicable to frozen samples that cannot produce intact single-cell suspensions. We have developed scONE-seq, a versatile scWGS-RNA-seq method that amplifies single-cell DNA and RNA without separating them from each other and hence is compatible with frozen biobanked samples. We benchmarked scONE-seq against existing methods using fresh and frozen samples to demonstrate its performance in various aspects. We identified a unique transcriptionally normal-like tumor clone by analyzing a 2-year frozen astrocytoma sample, demonstrating that performing single-cell multi-omics interrogation on biobanked tissue by scONE-seq could enable previously unidentified discoveries in tumor biology.

MEDICC2: whole-genome doubling aware copy-number phylogenies for cancer evolution

Aneuploidy, chromosomal instability, somatic copy-number alterations, and whole-genome doubling (WGD) play key roles in cancer evolution and provide information for the complex task of phylogenetic inference. We present MEDICC2, a method for inferring evolutionary trees and WGD using haplotype-specific somatic copy-number alterations from single-cell or bulk data. MEDICC2 eschews simplifications such as the infinite sites assumption, allowing multiple mutations and parallel evolution, and does not treat adjacent loci as independent, allowing overlapping copy-number events. Using simulations and multiple data types from 2780 tumors, we use MEDICC2 to demonstrate accurate inference of phylogenies, clonal and subclonal WGD, and ancestral copy-number states.

Genomic Profiling Identifies Putative Pathogenic Alterations in NSCLC Brain Metastases

Introduction

Brain metastases (BM) severely affect the prognosis and quality of life of patients with NSCLC. Recently, molecularly targeted agents were found to have promising activity against BM in patients with NSCLC whose primary tumors carry "druggable" mutations. Nevertheless, it remains critical to identify specific pathogenic alterations that drive NSCLC-BM and that can provide novel and more effective therapeutic targets.

Methods

To identify potentially targetable pathogenic alterations in NSCLC-BM, we profiled somatic copy number alterations (SCNAs) in 51 matched pairs of primary NSCLC and BM samples from 33 patients with lung adenocarcinoma and 18 patients with lung squamous cell carcinoma. In addition, we performed multiregion copy number profiling on 15 BM samples and whole-exome sequencing on 40 of 51 NSCLC-BM pairs.

Results

BM consistently had a higher burden of SCNAs compared with the matched primary tumors, and SCNAs were typically homogeneously distributed within BM, suggesting BM do not undergo extensive evolution once formed. By comparing focal SCNAs in matched NSCLC-BM pairs, we identified putative BM-driving alterations affecting multiple cancer genes, including several potentially targetable alterations in genes such as CDK12, DDR2, ERBB2, and NTRK1, which we validated in an independent cohort of 84 BM samples. Finally, we identified putative pathogenic alterations in multiple cancer genes, including genes involved in epigenome editing and 3D genome organization, such as EP300, CTCF, and STAG2, which we validated by targeted sequencing of an independent cohort of 115 BM samples.

Conclusions

Our study represents the most comprehensive genomic characterization of NSCLC-BM available to date, paving the way to functional studies aimed at assessing the potential of the identified pathogenic alterations as clinical biomarkers and targets.

Chromosomal Instability, Selection and Competition: Factors That Shape the Level of Karyotype Intra-Tumor Heterogeneity

Simple Summary

Each cancer consists of billions of cells. These cells are far from identical; hence, the population of cells that constitute a tumor is heterogeneous. A salient property that varies between cells in a tumor is their karyotype, the number and configuration of the chromosomes. The level of karyotype heterogeneity can be used to predict the survival of a patient. In this review, we describe the processes that shape the level of karyotype heterogeneity in a cancer.

Abstract

Intra-tumor heterogeneity (ITH) is a pan-cancer predictor of survival, with high ITH being correlated to a dismal prognosis. The level of ITH is, hence, a clinically relevant characteristic of a malignancy. ITH of karyotypes is driven by chromosomal instability (CIN). However, not all new karyotypes generated by CIN are viable or competitive, which limits the amount of ITH. Here, we review the cellular processes and ecological properties that determine karyotype ITH. We propose a framework to understand karyotype ITH, in which cells with new karyotypes emerge through CIN, are selected by cell intrinsic and cell extrinsic selective pressures, and propagate through a cancer in competition with other malignant cells. We further discuss how CIN modulates the cell phenotype and immune microenvironment, and the implications this has for the subsequent selection of karyotypes. Together, we aim to provide a comprehensive overview of the biological processes that shape the level of karyotype heterogeneity.

Patterns of metastases progression- The linear parallel ratio

Background

Linear and parallel are the two leading models of metastatic progression. In this study we propose a simple way to differentiate between them. While the linear model predicts accumulation of genetic and epigenetic alterations within the primary tumor by founder cells before spreading as waves of metastases, the parallel model suggests preclinical distribution of less advanced disseminated tumor cells with independent selection and expansion at the ectopic sites. Due to identical clonal origin and time of dispatching, linear metastases are expected to have comparable diameters in any specific organ while parallel metastases are expected to appear in variable sizes.

Methods and findings

Retrospective revision of chest CT of oncological patients with lung metastases was performed. Metastasis number and largest diameters were recorded. The sum number of metastases with a similar diameter (c) and those without (i) was counted and the linear/parallel ratio (LPR) was calculated for each patient using the formula (∑c-∑i)/(∑c+∑i). A LPR ratio of 1 implies pure linear progression pattern and -1 pure parallel. 12,887 metastases were measured in 503 patients with nine malignancy types. The median LPR of the entire group was 0.71 (IQR 0.14–0.93). In carcinomas of the pancreas, prostate, and thyroid the median LPR was 1. Median LPRs were 0.91, 0.65, 0.60, 0.58, 0.50 and 0.43 in renal cell carcinomas, melanomas, colorectal, breast, bladder, and sarcomas, respectively.

Conclusions

Metastatic spread of thyroid, pancreas, and prostate tumors is almost exclusively by a linear route. The spread of kidney, melanoma, colorectal, breast, bladder and sarcoma is both linear and parallel with increasing dominance of the parallel route in this order. These findings can explain and predict the clinical and genomic features of these tumors and can potentially be used for evaluation of metastatic origin in the individual patient.

The role of aging in cancer

Many cancers show an increase in incidence with age, and age is the biggest single risk factor for many cancers. However, the molecular basis of this relationship is poorly understood. Through a collection of review articles, our thematic issue discusses the link between aging and cancer in aspects including somatic mutations, proteostasis, mitochondria, metabolism, senescence, epigenetic regulation, immune regulation, DNA damage, and telomere function.

Modes of Selection in Tumors as Reflected by Two Mathematical Models and Site Frequency Spectra

The tug-of-war model was developed in a series of papers of McFarland and co-authors to account for existence of mutually counteracting rare advantageous driver mutations and more frequent slightly deleterious passenger mutations in cancer. In its original version, it was a state-dependent branching process. Because of its formulation, the tug-of-war model is of importance for tackling the problem as to whether evolution of cancerous tumors is "Darwinian" or "non-Darwinian." We define two Time-Continuous Markov Chain versions of the model, including identical mutation processes but adopting different drift and selection components. In Model A, drift and selection process preserves expected fitness whereas in Model B it leads to non-decreasing expected fitness. We investigate these properties using mathematical analysis and extensive simulations, which detect the effect of the so-called drift barrier in Model B but not in Model A. These effects are reflected in different structure of clone genealogies in the two models. Our work is related to the past theoretical work in the field of evolutionary genetics, concerning the interplay among mutation, drift and selection, in absence of recombination (asexual reproduction), where epistasis plays a major role. Finally, we use the statistics of mutation frequencies known as the Site Frequency Spectra (SFS), to compare the variant frequencies in DNA of sequenced HER2+ breast cancers, to those based on Model A and B simulations. The tumor-based SFS are better reproduced by Model A, pointing out a possible selection pattern of HER2+ tumor evolution. To put our models in context, we carried out an exploratory study of how publicly accessible data from breast, prostate, skin and ovarian cancers fit a range of models found in the literature.

Untangling the web of intratumour heterogeneity

Intratumour heterogeneity (ITH) is a hallmark of cancer that drives tumour evolution and disease progression. Technological and computational advances have enabled us to assess ITH at unprecedented depths, yet this accumulating knowledge has not had a substantial clinical impact. This is in part due to a limited understanding of the functional relevance of ITH and the inadequacy of preclinical experimental models to reproduce it. Here, we discuss progress made in these areas and illuminate future directions.

J-SPACE: a Julia package for the simulation of spatial models of cancer evolution and of sequencing experiments

Background

The combined effects of biological variability and measurement-related errors on cancer sequencing data remain largely unexplored. However, the spatio-temporal simulation of multi-cellular systems provides a powerful instrument to address this issue. In particular, efficient algorithmic frameworks are needed to overcome the harsh trade-off between scalability and expressivity, so to allow one to simulate both realistic cancer evolution scenarios and the related sequencing experiments, which can then be used to benchmark downstream bioinformatics methods.

Result

We introduce a Julia package for SPAtial Cancer Evolution (J-SPACE), which allows one to model and simulate a broad set of experimental scenarios, phenomenological rules and sequencing settings.Specifically, J-SPACE simulates the spatial dynamics of cells as a continuous-time multi-type birth-death stochastic process on a arbitrary graph, employing different rules of interaction and an optimised Gillespie algorithm. The evolutionary dynamics of genomic alterations (single-nucleotide variants and indels) is simulated either under the Infinite Sites Assumption or several different substitution models, including one based on mutational signatures. After mimicking the spatial sampling of tumour cells, J-SPACE returns the related phylogenetic model, and allows one to generate synthetic reads from several Next-Generation Sequencing (NGS) platforms, via the ART read simulator. The results are finally returned in standard FASTA, FASTQ, SAM, ALN and Newick file formats.

Conclusion

J-SPACE is designed to efficiently simulate the heterogeneous behaviour of a large number of cancer cells and produces a rich set of outputs. Our framework is useful to investigate the emergent spatial dynamics of cancer subpopulations, as well as to assess the impact of incomplete sampling and of experiment-specific errors. Importantly, the output of J-SPACE is designed to allow the performance assessment of downstream bioinformatics pipelines processing NGS data. J-SPACE is freely available at: https://github.com/BIMIB-DISCo/J-Space.jl.