Single-cell DNA sequencing reveals a late-dissemination model in metastatic colorectal cancer

Inference of chromosome selection parameters and missegregation rate in cancer from DNA-sequencing data

Aneuploidy is frequently observed in cancers and has been linked to poor patient outcome. Analysis of aneuploidy in DNA-sequencing (DNA-seq) data necessitates untangling the effects of the Copy Number Aberration (CNA) occurrence rates and the selection coefficients that act upon the resulting karyotypes. We introduce a parameter inference algorithm that takes advantage of both bulk and single-cell DNA-seq cohorts. The method is based on Approximate Bayesian Computation (ABC) and utilizes CINner, our recently introduced simulation algorithm of chromosomal instability in cancer. We examine three groups of statistics to summarize the data in the ABC routine: (A) Copy Number-based measures, (B) phylogeny tip statistics, and (C) phylogeny balance indices. Using these statistics, our method can recover both the CNA probabilities and selection parameters from ground truth data, and performs well even for data cohorts of relatively small sizes. We find that only statistics in groups A and C are well-suited for identifying CNA probabilities, and only group A carries the signals for estimating selection parameters. Moreover, the low number of CNA events at large scale compared to cell counts in single-cell samples means that statistics in group B cannot be estimated accurately using phylogeny reconstruction algorithms at the chromosome level. As data from both bulk and single-cell DNA-sequencing techniques becomes increasingly available, our inference framework promises to facilitate the analysis of distinct cancer types, differentiation between selection and neutral drift, and prediction of cancer clonal dynamics.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Cancer stem cells: advances in knowledge and implications for cancer therapy

Cancer stem cells (CSCs), a small subset of cells in tumors that are characterized by self-renewal and continuous proliferation, lead to tumorigenesis, metastasis, and maintain tumor heterogeneity. Cancer continues to be a significant global disease burden. In the past, surgery, radiotherapy, and chemotherapy were the main cancer treatments. The technology of cancer treatments continues to develop and advance, and the emergence of targeted therapy, and immunotherapy provides more options for patients to a certain extent. However, the limitations of efficacy and treatment resistance are still inevitable. Our review begins with a brief introduction of the historical discoveries, original hypotheses, and pathways that regulate CSCs, such as WNT/β-Catenin, hedgehog, Notch, NF-κB, JAK/STAT, TGF-β, PI3K/AKT, PPAR pathway, and their crosstalk. We focus on the role of CSCs in various therapeutic outcomes and resistance, including how the treatments affect the content of CSCs and the alteration of related molecules, CSCs-mediated therapeutic resistance, and the clinical value of targeting CSCs in patients with refractory, progressed or advanced tumors. In summary, CSCs affect therapeutic efficacy, and the treatment method of targeting CSCs is still difficult to determine. Clarifying regulatory mechanisms and targeting biomarkers of CSCs is currently the mainstream idea.

A war on many fronts: cross disciplinary approaches for novel cancer treatment strategies

Cancer is a disease characterized by uncontrolled cellular growth where cancer cells take advantage of surrounding cellular populations to obtain resources and promote invasion. Carcinomas are the most common type of cancer accounting for almost 90% of cancer cases. One of the major subtypes of carcinomas are adenocarcinomas, which originate from glandular cells that line certain internal organs. Cancers such as breast, prostate, lung, pancreas, colon, esophageal, kidney are often adenocarcinomas. Current treatment strategies include surgery, chemotherapy, radiation, targeted therapy, and more recently immunotherapy. However, patients with adenocarcinomas often develop resistance or recur after the first line of treatment. Understanding how networks of tumor cells interact with each other and the tumor microenvironment is crucial to avoid recurrence, resistance, and high-dose therapy toxicities. In this review, we explore how mathematical modeling tools from different disciplines can aid in the development of effective and personalized cancer treatment strategies. Here, we describe how concepts from the disciplines of ecology and evolution, economics, and control engineering have been applied to mathematically model cancer dynamics and enhance treatment strategies.

Heterogeneity and molecular landscape of melanoma: implications for targeted therapy

Uveal cancer (UM) offers a complex molecular landscape characterized by substantial heterogeneity, both on the genetic and epigenetic levels. This heterogeneity plays a critical position in shaping the behavior and response to therapy for this uncommon ocular malignancy. Targeted treatments with gene-specific therapeutic molecules may prove useful in overcoming radiation resistance, however, the diverse molecular makeups of UM call for a patient-specific approach in therapy procedures. We need to understand the intricate molecular landscape of UM to develop targeted treatments customized to each patient's specific genetic mutations. One of the promising approaches is using liquid biopsies, such as circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA), for detecting and monitoring the disease at the early stages. These non-invasive methods can help us identify the most effective treatment strategies for each patient. Single-cellular is a brand-new analysis platform that gives treasured insights into diagnosis, prognosis, and remedy. The incorporation of this data with known clinical and genomics information will give a better understanding of the complicated molecular mechanisms that UM diseases exploit. In this review, we focused on the heterogeneity and molecular panorama of UM, and to achieve this goal, the authors conducted an exhaustive literature evaluation spanning 1998 to 2023, using keywords like "uveal melanoma, "heterogeneity". "Targeted therapies"," "CTCs," and "single-cellular analysis".

Molecular perspectives on systemic priming and concomitant immunity in colorectal carcinoma

The progression of metastasis, a complex systemic disease, is facilitated by interactions between tumor cells and their isolated microenvironments. Over the past few decades, researchers have investigated the metastatic spread of cancer extensively, identifying multiple stages in the process, such as intravasation, extravasation, tumor latency, and the development of micrometastasis and macrometastasis. The premetastatic niche is established in target organs by the accumulation of aberrant immune cells and extracellular matrix proteins. The "seed and soil" idea, which has become widely known and accepted, is being used to this day to guide cancer studies. Changes in the local and systemic immune systems have a major impact on whether an infection spreads or not. The belief that the immune response may play a role in slowing tumor growth and may be beneficial against the metastatic disease underpins the responsiveness shown in the immunological landscape of metastasis. Various hypotheses on the phylogenesis of metastases have been proposed in the past. The primary tumor's secreting factors shape the intratumoral microenvironment and the immune landscape, allowing this progress to be made. Therefore, it is evident that among disseminated tumor cells, there are distinct phenotypes that either carry budding for metastasis or have the ability to obtain this potential or in systemic priming through contact with substantial metastatic niches that have implications for medicinal chemistry. Concurrent immunity signals that the main tumor induces an immune response that may not be strong enough to eradicate the tumor. Immunotherapy's success with some cancer patients shows that it is possible to effectively destroy even advanced-stage tumors by modifying the microenvironment and tumor-immune cell interactions. This review focuses on the metastasome in colorectal carcinoma and the therapeutic implications of site-specific metastasis, systemic priming, tumor spread, and the relationship between the immune system and metastasis.

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.

Inferring Markov Chains to Describe Convergent Tumor Evolution With CIMICE

The field of tumor phylogenetics focuses on studying the differences within cancer cell populations. Many efforts are done within the scientific community to build cancer progression models trying to understand the heterogeneity of such diseases. These models are highly dependent on the kind of data used for their construction, therefore, as the experimental technologies evolve, it is of major importance to exploit their peculiarities. In this work we describe a cancer progression model based on Single Cell DNA Sequencing data. When constructing the model, we focus on tailoring the formalism on the specificity of the data. We operate by defining a minimal set of assumptions needed to reconstruct a flexible DAG structured model, capable of identifying progression beyond the limitation of the infinite site assumption. Our proposal is conservative in the sense that we aim to neither discard nor infer knowledge which is not represented in the data. We provide simulations and analytical results to show the features of our model, test it on real data, show how it can be integrated with other approaches to cope with input noise. Moreover, our framework can be exploited to produce simulated data that follows our theoretical assumptions. Finally, we provide an open source R implementation of our approach, called CIMICE, that is publicly available on BioConductor.

Integrating SNVs and CNAs on a phylogenetic tree from single-cell DNA sequencing data

Single-cell DNA sequencing enables the construction of evolutionary trees that can reveal how tumors gain mutations and grow. Different whole-genome amplification procedures render genomic materials of different characteristics, often suitable for the detection of either single-nucleotide variation or copy number aberration, but not ideally for both. Consequently, this hinders the inference of a comprehensive phylogenetic tree and limits opportunities to investigate the interplay of SNVs and CNAs. Existing methods such as SCARLET and COMPASS require that the SNVs and CNAs are detected from the same sets of cells, which is technically challenging. Here we present a novel computational tool, SCsnvcna, that places SNVs on a tree inferred from CNA signals, whereas the sets of cells rendering the SNVs and CNAs are independent, offering a more practical solution in terms of the technical challenges. SCsnvcna is a Bayesian probabilistic model using both the genotype constraints on the tree and the cellular prevalence to search the optimal solution. Comprehensive simulations and comparison with seven state-of-the-art methods show that SCsnvcna is robust and accurate in a variety of circumstances. Particularly, SCsnvcna most frequently produces the lowest error rates, with ability to scale to a wide range of numerical values for leaf nodes in the tree, SNVs, and SNV cells. The application of SCsnvcna to two published colorectal cancer data sets shows highly consistent placement of SNV cells and SNVs with the original study while also supporting a refined placement of ATP7B, illustrating SCsnvcna's value in analyzing complex multitumor samples.

ConDoR: tumor phylogeny inference with a copy-number constrained mutation loss model

A tumor contains a diverse collection of somatic mutations that reflect its past evolutionary history and that range in scale from single nucleotide variants (SNVs) to large-scale copy-number aberrations (CNAs). However, no current single-cell DNA sequencing (scDNA-seq) technology produces accurate measurements of both SNVs and CNAs, complicating the inference of tumor phylogenies. We introduce a new evolutionary model, the constrained k-Dollo model, that uses SNVs as phylogenetic markers but constrains losses of SNVs according to clusters of cells. We derive an algorithm, ConDoR, that infers phylogenies from targeted scDNA-seq data using this model. We demonstrate the advantages of ConDoR on simulated and real scDNA-seq data.

Assessing the performance of methods for cell clustering from single-cell DNA sequencing data

BACKGROUND

Many cancer genomes have been known to contain more than one subclone inside one tumor, the phenomenon of which is called intra-tumor heterogeneity (ITH). Characterizing ITH is essential in designing treatment plans, prognosis as well as the study of cancer progression. Single-cell DNA sequencing (scDNAseq) has been proven effective in deciphering ITH. Cells corresponding to each subclone are supposed to carry a unique set of mutations such as single nucleotide variations (SNV). While there have been many studies on the cancer evolutionary tree reconstruction, not many have been proposed that simply characterize the subclonality without tree reconstruction. While tree reconstruction is important in the study of cancer evolutionary history, typically they are computationally expensive in terms of running time and memory consumption due to the huge search space of the tree structure. On the other hand, subclonality characterization of single cells can be converted into a cell clustering problem, the dimension of which is much smaller, and the turnaround time is much shorter. Despite the existence of a few state-of-the-art cell clustering computational tools for scDNAseq, there lacks a comprehensive and objective comparison under different settings.

RESULTS

In this paper, we evaluated six state-of-the-art cell clustering tools-SCG, BnpC, SCClone, RobustClone, SCITE and SBMClone-on simulated data sets given a variety of parameter settings and a real data set. We designed a simulator specifically for cell clustering, and compared these methods' performances in terms of their clustering accuracy, specificity and sensitivity and running time. For SBMClone, we specifically designed an ultra-low coverage large data set to evaluate its performance in the face of an extremely high missing rate.

CONCLUSION

From the benchmark study, we conclude that BnpC and SCG's clustering accuracy are the highest and comparable to each other. However, BnpC is more advantageous in terms of running time when cell number is high (> 1500). It also has a higher clustering accuracy than SCG when cluster number is high (> 16). SCClone's accuracy in estimating the number of clusters is the highest. RobustClone and SCITE's clustering accuracy are the lowest for all experiments. SCITE tends to over-estimate the cluster number and has a low specificity, whereas RobustClone tends to under-estimate the cluster number and has a much lower sensitivity than other methods. SBMClone produced reasonably good clustering (V-measure > 0.9) when coverage is > = 0.03 and thus is highly recommended for ultra-low coverage large scDNAseq data sets.

[Advances of spatial omics in the individualized diagnosis and treatment of head and neck cancer]

Spatialomics is another research hotspot of biotechnology after single-cell sequencing technology, which can make up for the defect that single-cell sequencing technology can not obtain cell spatial distribution information. Spatialomics mainly studies the relative position of cells in tissue samples to reveal the effect of cell spatial distribution on diseases. In recent years, spatialomics has made new progress in the pathogenesis, target exploration, drug development and many other aspects of head and neck tumors. This paper summarizes the latest progress of spatialomics in the diagnosis and treatment of head and neck cancer.

Barcoded multiple displacement amplification for high coverage sequencing in spatial genomics

Determining mutational landscapes in a spatial context is essential for understanding genetically heterogeneous cell microniches. Current approaches, such as Multiple Displacement Amplification (MDA), offer high genome coverage but limited multiplexing, which hinders large-scale spatial genomic studies. Here, we introduce barcoded MDA (bMDA), a technique that achieves high-coverage genomic analysis of low-input DNA while enhancing the multiplexing capabilities. By incorporating cell barcodes during MDA, bMDA streamlines library preparation in one pot, thereby overcoming a key bottleneck in spatial genomics. We apply bMDA to the integrative spatial analysis of triple-negative breast cancer tissues by examining copy number alterations, single nucleotide variations, structural variations, and kataegis signatures for each spatial microniche. This enables the assessment of subclonal evolutionary relationships within a spatial context. Therefore, bMDA has emerged as a scalable technology with the potential to advance the field of spatial genomics significantly.

Improving cellular phylogenies through the integrated use of mutation order and optimality principles

The study of tumor evolution is being revolutionalized by single-cell sequencing technologies that survey the somatic variation of cancer cells. In these endeavors, reliable inference of the evolutionary relationship of single cells is a key step. However, single-cell sequences contain many errors and missing bases, which necessitate advancing standard molecular phylogenetics approaches for applications in analyzing these datasets. We have developed a computational approach that integratively applies standard phylogenetic optimality principles and patterns of co-occurrence of sequence variations to produce more expansive and accurate cellular phylogenies from single-cell sequence datasets. We found the new approach to also perform well for CRISPR/Cas9 genome editing datasets, suggesting that it can be useful for various applications. We apply the new approach to some empirical datasets to showcase its use for reconstructing recurrent mutations and mutational reversals as well as for phylodynamics analysis to infer metastatic cell migrations between tumors.

A weighted distance-based approach for deriving consensus tumor evolutionary trees

MOTIVATION

The acquisition of somatic mutations by a tumor can be modeled by a type of evolutionary tree. However, it is impossible to observe this tree directly. Instead, numerous algorithms have been developed to infer such a tree from different types of sequencing data. But such methods can produce conflicting trees for the same patient, making it desirable to have approaches that can combine several such tumor trees into a consensus or summary tree. We introduce The Weighted m-Tumor Tree Consensus Problem (W-m-TTCP) to find a consensus tree among multiple plausible tumor evolutionary histories, each assigned a confidence weight, given a specific distance measure between tumor trees. We present an algorithm called TuELiP that is based on integer linear programming which solves the W-m-TTCP, and unlike other existing consensus methods, allows the input trees to be weighted differently.

RESULTS

On simulated data we show that TuELiP outperforms two existing methods at correctly identifying the true underlying tree used to create the simulations. We also show that the incorporation of weights can lead to more accurate tree inference. On a Triple-Negative Breast Cancer dataset, we show that including confidence weights can have important impacts on the consensus tree identified.

AVAILABILITY

An implementation of TuELiP and simulated datasets are available at https://bitbucket.org/oesperlab/consensus-ilp/src/main/.

Insights into the metastatic cascade through research autopsies

Metastasis is a complex process and the leading cause of cancer-related death globally. Recent studies have demonstrated that genomic sequencing data from paired primary and metastatic tumours can be used to trace the evolutionary origins of cells responsible for metastasis. This approach has yielded new insights into the genomic alterations that engender metastatic potential, and the mechanisms by which cancer spreads. Given that the reliability of these approaches is contingent upon how representative the samples are of primary and metastatic tumour heterogeneity, we review insights from studies that have reconstructed the evolution of metastasis within the context of their cohorts and designs. We discuss the role of research autopsies in achieving the comprehensive sampling necessary to advance the current understanding of metastasis.

Single-cell technologies in multiple myeloma: new insights into disease pathogenesis and translational implications

Multiple myeloma (MM) is a hematological malignancy characterized by clonal proliferation of plasma cells. Although therapeutic advances have been made to improve clinical outcomes and to prolong patients' survival in the past two decades, MM remains largely incurable. Single-cell sequencing (SCS) is a powerful method to dissect the cellular and molecular landscape at single-cell resolution, instead of providing averaged results. The application of single-cell technologies promises to address outstanding questions in myeloma biology and has revolutionized our understanding of the inter- and intra-tumor heterogeneity, tumor microenvironment, and mechanisms of therapeutic resistance in MM. In this review, we summarize the recently developed SCS methodologies and latest MM research progress achieved by single-cell profiling, including information regarding the cancer and immune cell landscapes, tumor heterogeneities, underlying mechanisms and biomarkers associated with therapeutic response and resistance. We also discuss future directions of applying transformative SCS approaches with contribution to clinical translation.

Testing for Phylogenetic Signal in Single-Cell RNA-Seq Data

Phylogenetic methods are emerging as a useful tool to understand cancer evolutionary dynamics, including tumor structure, heterogeneity, and progression. Most currently used approaches utilize either bulk whole genome sequencing or single-cell DNA sequencing and are based on calling copy number alterations and single nucleotide variants (SNVs). Single-cell RNA sequencing (scRNA-seq) is commonly applied to explore differential gene expression of cancer cells throughout tumor progression. The method exacerbates the single-cell sequencing problem of low yield per cell with uneven expression levels. This accounts for low and uneven sequencing coverage and makes SNV detection and phylogenetic analysis challenging. In this article, we demonstrate for the first time that scRNA-seq data contain sufficient evolutionary signal and can also be utilized in phylogenetic analyses. We explore and compare results of such analyses based on both expression levels and SNVs called from scRNA-seq data. Both techniques are shown to be useful for reconstructing phylogenetic relationships between cells, reflecting the clonal composition of a tumor. Both standardized expression values and SNVs appear to be equally capable of reconstructing a similar pattern of phylogenetic relationship. This pattern is stable even when phylogenetic uncertainty is taken in account. Our results open up a new direction of somatic phylogenetics based on scRNA-seq data. Further research is required to refine and improve these approaches to capture the full picture of somatic evolutionary dynamics in cancer.

State-dependent evolutionary models reveal modes of solid tumour growth

Spatial properties of tumour growth have profound implications for cancer progression, therapeutic resistance and metastasis. Yet, how spatial position governs tumour cell division remains difficult to evaluate in clinical tumours. Here, we demonstrate that faster division on the tumour periphery leaves characteristic genetic patterns, which become evident when a phylogenetic tree is reconstructed from spatially sampled cells. Namely, rapidly dividing peripheral lineages branch more extensively and acquire more mutations than slower-dividing centre lineages. We develop a Bayesian state-dependent evolutionary phylodynamic model (SDevo) that quantifies these patterns to infer the differential division rates between peripheral and central cells. We demonstrate that this approach accurately infers spatially varying birth rates of simulated tumours across a range of growth conditions and sampling strategies. We then show that SDevo outperforms state-of-the-art, non-cancer multi-state phylodynamic methods that ignore differential sequence evolution. Finally, we apply SDevo to single-time-point, multi-region sequencing data from clinical hepatocellular carcinomas and find evidence of a three- to six-times-higher division rate on the tumour edge. With the increasing availability of high-resolution, multi-region sequencing, we anticipate that SDevo will be useful in interrogating spatial growth restrictions and could be extended to model non-spatial factors that influence tumour progression.

Orchestration of Collective Migration and Metastasis by Tumor Cell Clusters

Metastatic dissemination has lethal consequences for cancer patients. Accruing evidence supports the hypothesis that tumor cells can migrate and metastasize as clusters of cells while maintaining contacts with one another. Collective metastasis enables tumor cells to colonize secondary sites more efficiently, resist cell death, and evade the immune system. On the other hand, tumor cell clusters face unique challenges for dissemination particularly during systemic dissemination. Here, we review recent progress toward understanding how tumor cell clusters overcome these disadvantages as well as mechanisms they utilize to gain advantages throughout the metastatic process. We consider useful models for studying collective metastasis and reflect on how the study of collective metastasis suggests new opportunities for eradicating and preventing metastatic disease.

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.

ConDoR: Tumor phylogeny inference with a copy-number constrained mutation loss model

Tumors consist of subpopulations of cells that harbor distinct collections of somatic mutations. These mutations range in scale from single nucleotide variants (SNVs) to large-scale copy-number aberrations (CNAs). While many approaches infer tumor phylogenies using SNVs as phylogenetic markers, CNAs that overlap SNVs may lead to erroneous phylogenetic inference. Specifically, an SNV may be lost in a cell due to a deletion of the genomic segment containing the SNV. Unfortunately, no current single-cell DNA sequencing (scDNA-seq) technology produces accurate measurements of both SNVs and CNAs. For instance, recent targeted scDNA-seq technologies, such as Mission Bio Tapestri, measure SNVs with high fidelity in individual cells, but yield much less reliable measurements of CNAs. We introduce a new evolutionary model, the constrained k-Dollo model, that uses SNVs as phylogenetic markers and partial information about CNAs in the form of clustering of cells with similar copy-number profiles. This copy-number clustering constrains where loss of SNVs can occur in the phylogeny. We develop ConDoR (Constrained Dollo Reconstruction), an algorithm to infer tumor phylogenies from targeted scDNA-seq data using the constrained k-Dollo model. We show that ConDoR outperforms existing methods on simulated data. We use ConDoR to analyze a new multi-region targeted scDNA-seq dataset of 2153 cells from a pancreatic ductal adenocarcinoma (PDAC) tumor and produce a more plausible phylogeny compared to existing methods that conforms to histological results for the tumor from a previous study. We also analyze a metastatic colorectal cancer dataset, deriving a more parsimonious phylogeny than previously published analyses and with a simpler monoclonal origin of metastasis compared to the original study.

Code availability

Software is available at https://github.com/raphael-group/constrained-Dollo.

A phylogenetic approach to inferring the order in which mutations arise during cancer progression

Although the role of evolutionary process in cancer progression is widely accepted, increasing attention is being given to the evolutionary mechanisms that can lead to differences in clinical outcome. Recent studies suggest that the temporal order in which somatic mutations accumulate during cancer progression is important. Single-cell sequencing (SCS) provides a unique opportunity to examine the effect that the mutation order has on cancer progression and treatment effect. However, the error rates associated with single-cell sequencing are known to be high, which greatly complicates the task. We propose a novel method for inferring the order in which somatic mutations arise within an individual tumor using noisy data from single-cell sequencing. Our method incorporates models at two levels in that the evolutionary process of somatic mutation within the tumor is modeled along with the technical errors that arise from the single-cell sequencing data collection process. Through analyses of simulations across a wide range of realistic scenarios, we show that our method substantially outperforms existing approaches for identifying mutation order. Most importantly, our method provides a unique means to capture and quantify the uncertainty in the inferred mutation order along a given phylogeny. We illustrate our method by analyzing data from colorectal and prostate cancer patients, in which our method strengthens previously reported mutation orders. Our work is an important step towards producing meaningful prediction of mutation order with high accuracy and measuring the uncertainty of predicted mutation order in cancer patients, with the potential to lead to new insights about the evolutionary trajectories of cancer.

Intracellular and extracellular factors of colorectal cancer liver metastasis: a pivotal perplex to be fully elucidated

Metastasis is the leading cause of death in colorectal cancer (CRC) patients, and the liver is the most common site of metastasis. Tumor cell metastasis can be thought of as an invasion-metastasis cascade and metastatic organotropism is thought to be a process that relies on the intrinsic properties of tumor cells and their interactions with molecules and cells in the microenvironment. Many studies have provided new insights into the molecular mechanism and contributing factors involved in CRC liver metastasis for a better understanding of the organ-specific metastasis process. The purpose of this review is to summarize the theories that explain CRC liver metastasis at multiple molecular dimensions (including genetic and non-genetic factors), as well as the main factors that cause CRC liver metastasis. Many findings suggest that metastasis may occur earlier than expected and with specific organ-anchoring property. The emergence of potential metastatic clones, the timing of dissemination, and the distinct routes of metastasis have been explained by genomic studies. The main force of CRC liver metastasis is also thought to be epigenetic alterations and dynamic phenotypic traits. Furthermore, we review key extrinsic factors that influence CRC cell metastasis and liver tropisms, such as pre-niches, tumor stromal cells, adhesion molecules, and immune/inflammatory responses in the tumor microenvironment. In addition, biomarkers associated with early diagnosis, prognosis, and recurrence of liver metastasis from CRC are summarized to enlighten potential clinical practice, including some markers that can be used as therapeutic targets to provide new perspectives for the treatment strategies of CRC liver metastasis.

Molecular pathogenesis, targeted therapies, and future perspectives for gastric cancer

Gastric cancer is a major source of global cancer mortality with limited treatment options and poor patient survival. As our molecular understanding of gastric cancer improves, we are now beginning to recognize that these cancers are a heterogeneous group of diseases with incredibly unique pathogeneses and active oncogenic pathways. It is this molecular diversity and oftentimes lack of common oncogenic driver mutations that bestow the poor treatment responses that oncologists often face when treating gastric cancer. In this review, we will examine the treatments for gastric cancer including up-to-date molecularly targeted therapies and immunotherapies. We will then review the molecular subtypes of gastric cancer to highlight the diversity seen in this disease. We will then shift our discussion to basic science and gastric cancer mouse models as tools to study gastric cancer molecular heterogeneity. Furthermore, we will elaborate on a molecular process termed paligenosis and the cyclical hit model as key events during gastric cancer initiation that impart nondividing mature differentiated cells the ability to re-enter the cell cycle and accumulate disparate genomic mutations during years of chronic inflammation and injury. As our basic science understanding of gastric cancer advances, so too must our translational and clinical efforts. We will end with a discussion regarding single-cell molecular analyses and cancer organoid technologies as future translational avenues to advance our understanding of gastric cancer heterogeneity and to design precision-based gastric cancer treatments. Elucidation of interpatient and intratumor heterogeneity is the only way to advance future cancer prevention, diagnoses and treatment.

Genetic and microenvironmental intra-tumor heterogeneity impacts colorectal cancer evolution and metastatic development

Colorectal cancer (CRC) is a highly diverse disease, where different genomic instability pathways shape genetic clonal diversity and tumor microenvironment. Although intra-tumor heterogeneity has been characterized in primary tumors, its origin and consequences in CRC outcome is not fully understood. Therefore, we assessed intra- and inter-tumor heterogeneity of a prospective cohort of 136 CRC samples. We demonstrate that CRC diversity is forged by asynchronous forms of molecular alterations, where mutational and chromosomal instability collectively boost CRC genetic and microenvironment intra-tumor heterogeneity. We were able to depict predictor signatures of cancer-related genes that can foresee heterogeneity levels across the different tumor consensus molecular subtypes (CMS) and primary tumor location. Finally, we show that high genetic and microenvironment heterogeneity are associated with lower metastatic potential, whereas late-emerging copy number variations favor metastasis development and polyclonal seeding. This study provides an exhaustive portrait of the interplay between genetic and microenvironment intra-tumor heterogeneity across CMS subtypes, depicting molecular events with predictive value of CRC progression and metastasis development.

Accurate tumor clonal structures require single-cell analysis

Tumor clonal structure is closely related to future progression, which has been mainly investigated as mutation abundance clustering in bulk samples. With relatively limited studies at single-cell resolution, a systematic comparison of the two approaches is still lacking. Here, using bulk and single-cell mutational data from the liver and colorectal cancers, we checked whether co-mutations determined by single-cell analysis had corresponding bulk variant allele frequency (VAF) peaks. While bulk analysis suggested the absence of subclonal peaks and, possibly, neutral evolution in some cases, the single-cell analysis identified coexisting subclones. The overlaps of bulk VAF ranges for co-mutations from different subclones made it difficult to separate them. Complex subclonal structures and dynamic evolution could be hidden under the seemingly clonal neutral pattern at the bulk level, suggesting single-cell analysis is necessary to avoid underestimation of tumor heterogeneity.

Intestinal cellular heterogeneity and disease development revealed by single-cell technology

The intestinal epithelium is responsible for food digestion and nutrient absorption and plays a critical role in hormone secretion, microorganism defense, and immune response. These functions depend on the integral single-layered intestinal epithelium, which shows diversified cell constitution and rapid self-renewal and presents powerful regeneration plasticity after injury. Derailment of homeostasis of the intestine epithelium leads to the development of diseases, most commonly including enteritis and colorectal cancer. Therefore, it is important to understand the cellular characterization of the intestinal epithelium at the molecular level and the mechanisms underlying its homeostatic maintenance. Single-cell technologies allow us to gain molecular insights at the single-cell level. In this review, we summarize the single-cell RNA sequencing applications to understand intestinal cell characteristics, spatiotemporal evolution, and intestinal disease development.

Putative anoikis resistant subpopulations are enriched in lymph node metastases and indicate adverse prognosis in colorectal carcinoma

Anoikis refers to apoptosis induced by the loss of contact with the extracellular matrix. Anoikis resistance is essential for metastasis. We have recently shown that it is possible to quantitatively evaluate putative anoikis resistant (AR) subpopulations in colorectal carcinoma (CRC). Abundance of these multi-cell structures is an independent marker of adverse prognosis. Here, we have quantified putative AR subpopulations in lymph node (LN) metastases of CRC and evaluated their prognostic value and relationship with the characteristics of primary tumors. A case series included 137 unselected CRC patients, 54 with LN metastases. Areal densities (structures/mm²) of putative AR structures in primary tumors had been analyzed previously and now were determined from all nodal metastases (n = 183). Areal density of putative AR structures was higher in LN metastases than in primary tumors. Variation of the areal density within different LN metastases of a single patient was lower than between metastases of different patients. Abundance of putative AR structures in LN metastases was associated with shorter cancer specific survival (p = 0.013), and this association was independent of T and N stages. Abundance of putative AR structures in primary tumors and LN metastases had a cumulative adverse effect on prognosis. Enrichment of putative AR subpopulations in LN metastases suggest that in metastasis formation, there is a selection favoring cells capable of forming these structures. Higher intra-case constancy relative to inter-case variation suggests that such selection is stable in metastasis development. Our findings indirectly support the biological validity of our concept of putative AR structures.

Single-cell genomic and transcriptomic landscapes of primary and metastatic colorectal cancer tumors

BACKGROUND

Colorectal cancer (CRC) ranks as the second-leading cause of cancer-related death worldwide with metastases being the main cause of cancer-related death. Here, we investigated the genomic and transcriptomic alterations in matching adjacent normal tissues, primary tumors, and metastatic tumors of CRC patients.

METHODS

We performed whole genome sequencing (WGS), multi-region whole exome sequencing (WES), simultaneous single-cell RNA-Seq, and single-cell targeted cDNA Sanger sequencing on matching adjacent normal tissues, primary tumors, and metastatic tumors from 12 metastatic colorectal cancer patients (n=84 for genomes, n=81 for exomes, n=9120 for single cells). Patient-derived tumor organoids were used to estimate the anti-tumor effects of a PPAR inhibitor, and self-renewal and differentiation ability of stem cell-like tumor cells.

RESULTS

We found that the PPAR signaling pathway was prevalently and aberrantly activated in CRC tumors. Blocking of PPAR pathway both suppressed the growth and promoted the apoptosis of CRC organoids in vitro, indicating that aberrant activation of the PPAR signaling pathway plays a critical role in CRC tumorigenesis. Using matched samples from the same patient, distinct origins of the metastasized tumors between lymph node and liver were revealed, which was further verified by both copy number variation and mitochondrial mutation profiles at single-cell resolution. By combining single-cell RNA-Seq and single-cell point mutation identification by targeted cDNA Sanger sequencing, we revealed important phenotypic differences between cancer cells with and without critical point mutations (KRAS and TP53) in the same patient in vivo at single-cell resolution.

CONCLUSIONS

Our data provides deep insights into how driver mutations interfere with the transcriptomic state of cancer cells in vivo at a single-cell resolution. Our findings offer novel knowledge on metastatic mechanisms as well as potential markers and therapeutic targets for CRC diagnosis and therapy. The high-precision single-cell RNA-seq dataset of matched adjacent normal tissues, primary tumors, and metastases from CRCs may serve as a rich resource for further studies.

Advances in Single-Cell Toxicogenomics in Environmental Toxicology

The toxicity evaluation system of environmental pollutants has undergone numerous changes due to the application of new technologies. Single-cell toxicogenomics is rapidly changing our view on environmental toxicology by increasing the resolution of our analysis to the level of a single cell. Applications of this technology in environmental toxicology have begun to emerge and are rapidly expanding the portfolio of existing technologies and applications. Here, we first summarized different methods involved in single-cell isolation and amplification in single-cell sequencing process, compared the advantages and disadvantages of different methods, and analyzed their development trends. Then, we reviewed the main advances of single-cell toxicogenomics in environmental toxicology, emphatically analyzed the application prospects of this technology in identifying the target cells of pollutants in early embryos, clarifying the heterogeneous response of cell subtypes to pollutants, and finding pathogenic bacteria in unknown microbes, and highlighted the unique characteristics of this approach with high resolution, high throughput, and high specificity by examples. We also offered a prediction of the further application of this technology and the revolution it brings in environmental toxicology. Overall, these advances will provide practical solutions for controlling or mitigating exogenous toxicological effects that threaten human and ecosystem health, contribute to improving our understanding of the physiological processes affected by pollutants, and lead to the emergence of new methods of pollution control.

Scelestial: Fast and accurate single-cell lineage tree inference based on a Steiner tree approximation algorithm

Single-cell genome sequencing provides a highly granular view of biological systems but is affected by high error rates, allelic amplification bias, and uneven genome coverage. This creates a need for data-specific computational methods, for purposes such as for cell lineage tree inference. The objective of cell lineage tree reconstruction is to infer the evolutionary process that generated a set of observed cell genomes. Lineage trees may enable a better understanding of tumor formation and growth, as well as of organ development for healthy body cells. We describe a method, Scelestial, for lineage tree reconstruction from single-cell data, which is based on an approximation algorithm for the Steiner tree problem and is a generalization of the neighbor-joining method. We adapt the algorithm to efficiently select a limited subset of potential sequences as internal nodes, in the presence of missing values, and to minimize cost by lineage tree-based missing value imputation. In a comparison against seven state-of-the-art single-cell lineage tree reconstruction algorithms—BitPhylogeny, OncoNEM, SCITE, SiFit, SASC, SCIPhI, and SiCloneFit—on simulated and real single-cell tumor samples, Scelestial performed best at reconstructing trees in terms of accuracy and run time. Scelestial has been implemented in C++. It is also available as an R package named RScelestial.

Application and prospect of single-cell sequencing in cancer metastasis

Cancer metastasis is a complicated process driven by internal genetic variations and developed through interactions with the external environment. This process usually causes therapeutic resistance and results in a low survival rate. In recent years, single-cell sequencing has become a popular method for revealing the tumor evolutionary genetic lineage, intra-tumoral heterogeneity and tumor microenvironment of the metastasis process. So as to find more therapeutic targets for clinical application, the spatial transcriptomics method has become a new rising field of cancer studies, which promotes the combination between clinical medicine and molecular biology. In future prospects, more accurate and personalized treatment models will come into reality.

Accounting for errors in data improves divergence time estimates in single-cell cancer evolution

Single-cell sequencing provides a new way to explore the evolutionary history of cells. Compared to traditional bulk sequencing, where a population of heterogeneous cells is pooled to form a single observation, single-cell sequencing isolates and amplifies genetic material from individual cells, thereby preserving the information about the origin of the sequences. However, single-cell data is more error-prone than bulk sequencing data due to the limited genomic material available per cell. Here, we present error and mutation models for evolutionary inference of single-cell data within a mature and extensible Bayesian framework, BEAST2. Our framework enables integration with biologically informative models such as relaxed molecular clocks and population dynamic models. Our simulations show that modeling errors increase the accuracy of relative divergence times and substitution parameters. We reconstruct the phylogenetic history of a colorectal cancer patient and a healthy patient from single-cell DNA sequencing data. We find that the estimated times of terminal splitting events are shifted forward in time compared to models which ignore errors. We observed that not accounting for errors can overestimate the phylogenetic diversity in single-cell DNA sequencing data. We estimate that 30-50% of the apparent diversity can be attributed to error. Our work enables a full Bayesian approach capable of accounting for errors in the data within the integrative Bayesian software framework BEAST2.

Applications and Achievements of Single-Cell Sequencing in Gastrointestinal Cancer

Gastrointestinal cancer represents a public health concern that seriously endangers human health. The emerging single-cell sequencing (SCS) technologies are different from the large-scale sequencing technologies which provide inaccurate data. SCS is a powerful tool for deciphering the single-cell resolutions of cellular and molecular landscapes, revealing the features of single-cell genomes, transcriptomes, and epigenomes. Recently, SCS has been applied in the field of gastrointestinal cancer research for clarifying the origin and heterogeneity of gastrointestinal cancer, acquiring micro-environmental information, and improving diagnostic and treatment methods. This review outlines the applications of SCS in gastrointestinal cancer research and summarizes the most recent advances in the field.

Identification, discrimination and heterogeneity of fibroblasts

Fibroblasts, the principal cell type of connective tissue, secrete extracellular matrix components during tissue development, homeostasis, repair and disease. Despite this crucial role, the identification and distinction of fibroblasts from other cell types are challenging and laden with caveats. Rapid progress in single-cell transcriptomics now yields detailed molecular portraits of fibroblasts and other cell types in our bodies, which complement and enrich classical histological and immunological descriptions, improve cell class definitions and guide further studies on the functional heterogeneity of cell subtypes and states, origins and fates in physiological and pathological processes. In this review, we summarize and discuss recent advances in the understanding of fibroblast identification and heterogeneity and how they discriminate from other cell types.

In this review, the authors look at how recent progress in single-cell transcriptomics complement and enrich the classical, largely morphological, portraits of fibroblasts. The detailed molecular information now available provides new insights into fibroblast identity, heterogeneity and function.

Somatic variant calling from single-cell DNA sequencing data

Single-cell sequencing has gained popularity in recent years. Despite its numerous applications, single-cell DNA sequencing data is highly error-prone due to technical biases arising from uneven sequencing coverage, allelic dropout, and amplification error. With these artifacts, the identification of somatic genomic variants becomes a challenging task, and over the years, several methods have been developed explicitly for this type of data. Single-cell variant callers implement distinct strategies, make different use of the data, and typically result in many discordant calls when applied to real data. Here, we review current approaches for single-cell variant calling, emphasizing single nucleotide variants. We highlight their potential benefits and shortcomings to help users choose a suitable tool for their data at hand.

Alignment of Fellowship Training with Practice Patterns for Molecular Pathologists: A Report of the Association for Molecular Pathology Training and Education Committee

In the two decades since Accreditation Council for Graduate Medical Education-accredited Molecular Genetic Pathology fellowships began, the field of clinical molecular pathology has evolved considerably. The American Board of Pathology gathered data from board-certified molecular genetic pathologists assessing the alignment of skills and knowledge gained during fellowship with current needs on the job. The Association of Molecular Pathology conducted a parallel survey of program directors, and included questions on how various topics were taught during fellowship, as well as ranking their importance. Both surveys showed that most training aligned well with the practice needs of former trainees. Genomic profiling of tumors by next-generation sequencing, bioinformatics, laboratory management, and regulatory issues were topics thought to require increased emphasis in training. Topics related to clinical genetics and microbiology were deemed less important by those in practice, perhaps reflecting the increasing subspecialization of molecular pathologists. Program directors still viewed these topics as important to provide foundational knowledge. Parentage, identity, and human leukocyte antigen testing were less important to both survey audiences. These data may be helpful in guiding future adjustments to the Molecular Genetic Pathology curriculum and Accreditation Council for Graduate Medical Education program requirements.

Pediatric Sarcomas: The Next Generation of Molecular Studies

Pediatric sarcomas constitute one of the largest groups of childhood cancers, following hematopoietic, neural, and renal lesions. Partly because of their diversity, they continue to offer challenges in diagnosis and treatment. In spite of the diagnostic, nosologic, and therapeutic gains made with genetic technology, newer means for investigation are needed. This article reviews emerging technology being used to study human neoplasia and how these methods might be applicable to pediatric sarcomas. Methods reviewed include single cell RNA sequencing (scRNAseq), spatial multi-omics, high-throughput functional genomics, and clustered regularly interspersed short palindromic sequence-Cas9 (CRISPR-Cas9) technology. In spite of these advances, the field continues to be challenged by a dearth of properly annotated materials, particularly from recurrences and metastases and pre- and post-treatment samples.

Breaking the Immune Complexity of the Tumor Microenvironment Using Single-Cell Technologies

Tumors are not a simple aggregate of transformed cells but rather a complicated ecosystem containing various components, including infiltrating immune cells, tumor-related stromal cells, endothelial cells, soluble factors, and extracellular matrix proteins. Profiling the immune contexture of this intricate framework is now mandatory to develop more effective cancer therapies and precise immunotherapeutic approaches by identifying exact targets or predictive biomarkers, respectively. Conventional technologies are limited in reaching this goal because they lack high resolution. Recent developments in single-cell technologies, such as single-cell RNA transcriptomics, mass cytometry, and multiparameter immunofluorescence, have revolutionized the cancer immunology field, capturing the heterogeneity of tumor-infiltrating immune cells and the dynamic complexity of tenets that regulate cell networks in the tumor microenvironment. In this review, we describe some of the current single-cell technologies and computational techniques applied for immune-profiling the cancer landscape and discuss future directions of how integrating multi-omics data can guide a new "precision oncology" advancement.

Single-cell sequencing: a promising approach for uncovering the mechanisms of tumor metastasis

Single-cell sequencing (SCS) is an emerging high-throughput technology that can be used to study the genomics, transcriptomics, and epigenetics at a single cell level. SCS is widely used in the diagnosis and treatment of various diseases, including cancer. Over the years, SCS has gradually become an effective clinical tool for the exploration of tumor metastasis mechanisms and the development of treatment strategies. Currently, SCS can be used not only to analyze metastasis-related malignant biological characteristics, such as tumor heterogeneity, drug resistance, and microenvironment, but also to construct metastasis-related cell maps for predicting and monitoring the dynamics of metastasis. SCS is also used to identify therapeutic targets related to metastasis as it provides insights into the distribution of tumor cell subsets and gene expression differences between primary and metastatic tumors. Additionally, SCS techniques in combination with artificial intelligence (AI) are used in liquid biopsy to identify circulating tumor cells (CTCs), thereby providing a novel strategy for treating tumor metastasis. In this review, we summarize the potential applications of SCS in the field of tumor metastasis and discuss the prospects and limitations of SCS to provide a theoretical basis for finding therapeutic targets and mechanisms of metastasis.

The Genetic Evolution of Metastasis

Cancer is an evolutionary process that is characterized by the emergence of multiple genetically distinct populations or clones within the primary tumor. Intratumor heterogeneity provides a substrate for the selection of adaptive clones, such as those that lead to metastasis. Comparative molecular studies of primary tumors and metastases have identified distinct genomic features associated with the development of metastases. In this review, we discuss how these insights could inform clinical decision-making and uncover rational antimetastasis treatment strategies.

Single-Cell Sequencing and Its Applications in Liver Cancer

As one of the most lethal cancers, primary liver cancer (PLC) has high tumor heterogeneity, including the heterogeneity between cancer cells. Traditional methods which have been used to identify tumor heterogeneity for a long time are based on large mixed cell samples, and the research results usually show average level of the cell population, ignoring the heterogeneity between cancer cells. In recent years, single-cell sequencing has been increasingly applied to the studies of PLCs. It can detect the heterogeneity between cancer cells, distinguish each cell subgroup in the tumor microenvironment (TME), and also reveal the clonal characteristics of cancer cells, contributing to understand the evolution of tumor. Here, we introduce the process of single-cell sequencing, review the applications of single-cell sequencing in the heterogeneity of cancer cells, TMEs, oncogenesis, and metastatic mechanisms of liver cancer, and discuss some of the current challenges in the field.

BiTSC 2: Bayesian inference of tumor clonal tree by joint analysis of single-cell SNV and CNA data

The rapid development of single-cell DNA sequencing (scDNA-seq) technology has greatly enhanced the resolution of tumor cell profiling, providing an unprecedented perspective in characterizing intra-tumoral heterogeneity and understanding tumor progression and metastasis. However, prominent algorithms for constructing tumor phylogeny based on scDNA-seq data usually only take single nucleotide variations (SNVs) as markers, failing to consider the effect caused by copy number alterations (CNAs). Here, we propose BiTSC$^2$, Bayesian inference of Tumor clonal Tree by joint analysis of Single-Cell SNV and CNA data. BiTSC$^2$ takes raw reads from scDNA-seq as input, accounts for the overlapping of CNA and SNV, models allelic dropout rate, sequencing errors and missing rate, as well as assigns single cells into subclones. By applying Markov Chain Monte Carlo sampling, BiTSC$^2$ can simultaneously estimate the subclonal scCNA and scSNV genotype matrices, subclonal assignments and tumor subclonal evolutionary tree. In comparison with existing methods on synthetic and real tumor data, BiTSC$^2$ shows high accuracy in genotype recovery, subclonal assignment and tree reconstruction. BiTSC$^2$ also performs robustly in dealing with scDNA-seq data with low sequencing depth and variant missing rate. BiTSC$^2$ software is available at https://github.com/ucasdp/BiTSC2.

Parsimonious Clone Tree Integration in cancer

BACKGROUND

Every tumor is composed of heterogeneous clones, each corresponding to a distinct subpopulation of cells that accumulated different types of somatic mutations, ranging from single-nucleotide variants (SNVs) to copy-number aberrations (CNAs). As the analysis of this intra-tumor heterogeneity has important clinical applications, several computational methods have been introduced to identify clones from DNA sequencing data. However, due to technological and methodological limitations, current analyses are restricted to identifying tumor clones only based on either SNVs or CNAs, preventing a comprehensive characterization of a tumor's clonal composition.

RESULTS

To overcome these challenges, we formulate the identification of clones in terms of both SNVs and CNAs as a integration problem while accounting for uncertainty in the input SNV and CNA proportions. We thus characterize the computational complexity of this problem and we introduce PACTION (PArsimonious Clone Tree integratION), an algorithm that solves the problem using a mixed integer linear programming formulation. On simulated data, we show that tumor clones can be identified reliably, especially when further taking into account the ancestral relationships that can be inferred from the input SNVs and CNAs. On 49 tumor samples from 10 prostate cancer patients, our integration approach provides a higher resolution view of tumor evolution than previous studies.

CONCLUSION

PACTION is an accurate and fast method that reconstructs clonal architecture of cancer tumors by integrating SNV and CNA clones inferred using existing methods.

SCClone: Accurate Clustering of Tumor Single-Cell DNA Sequencing Data

Single-cell DNA sequencing (scDNA-seq) enables high-resolution profiling of genetic diversity among single cells and is especially useful for deciphering the intra-tumor heterogeneity and evolutionary history of tumor. Specific technical issues such as allele dropout, false-positive errors, and doublets make scDNA-seq data incomplete and error-prone, giving rise to a severe challenge of accurately inferring clonal architecture of tumor. To effectively address these issues, we introduce a new computational method called SCClone for reasoning subclones from single nucleotide variation (SNV) data of single cells. Specifically, SCClone leverages a probability mixture model for binary data to cluster single cells into distinct subclones. To accurately decipher underlying clonal composition, a novel model selection scheme based on inter-cluster variance is employed to find the optimal number of subclones. Extensive evaluations on various simulated datasets suggest SCClone has strong robustness against different technical noises in scDNA-seq data and achieves better performance than the state-of-the-art methods in reasoning clonal composition. Further evaluations of SCClone on three real scDNA-seq datasets show that it can effectively find the underlying subclones from severely disturbed data. The SCClone software is freely available at https://github.com/qasimyu/scclone.

CellPhy: accurate and fast probabilistic inference of single-cell phylogenies from scDNA-seq data

We introduce CellPhy, a maximum likelihood framework for inferring phylogenetic trees from somatic single-cell single-nucleotide variants. CellPhy leverages a finite-site Markov genotype model with 16 diploid states and considers amplification error and allelic dropout. We implement CellPhy into RAxML-NG, a widely used phylogenetic inference package that provides statistical confidence measurements and scales well on large datasets with hundreds or thousands of cells. Comprehensive simulations suggest that CellPhy is more robust to single-cell genomics errors and outperforms state-of-the-art methods under realistic scenarios, both in accuracy and speed. CellPhy is freely available at https://github.com/amkozlov/cellphy .

Deciphering the Immune-Tumor Interplay During Early-Stage Lung Cancer Development via Single-Cell Technology

Lung cancer is the leading cause of cancer-related death worldwide. Cancer immunotherapy has shown great success in treating advanced-stage lung cancer but has yet been used to treat early-stage lung cancer, mostly due to lack of understanding of the tumor immune microenvironment in early-stage lung cancer. The immune system could both constrain and promote tumorigenesis in a process termed immune editing that can be divided into three phases, namely, elimination, equilibrium, and escape. Current understanding of the immune response toward tumor is mainly on the "escape" phase when the tumor is clinically detectable. The detailed mechanism by which tumor progenitor lesions was modulated by the immune system during early stage of lung cancer development remains elusive. The advent of single-cell sequencing technology enables tumor immunologists to address those fundamental questions. In this perspective, we will summarize our current understanding and big gaps about the immune response during early lung tumorigenesis. We will then present the state of the art of single-cell technology and then envision how single-cell technology could be used to address those questions. Advances in the understanding of the immune response and its dynamics during malignant transformation of pre-malignant lesion will shed light on how malignant cells interact with the immune system and evolve under immune selection. Such knowledge could then contribute to the development of precision and early intervention strategies toward lung malignancy.