MQuad enables clonal substructure discovery using single cell mitochondrial variants

MitoSort: Robust Demultiplexing of Pooled Single-cell Genomic Data Using Endogenous Mitochondrial Variants

Multiplexing across donors has emerged as a popular strategy to increase throughput, reduce costs, overcome technical batch effects, and improve doublet detection in single-cell genomic studies. To eliminate additional experimental steps, endogenous nuclear genome variants are used for demultiplexing pooled single-cell RNA sequencing (scRNA-seq) data by several computational tools. However, these tools have limitations when applied to single-cell sequencing methods that do not cover nuclear genomic regions well, such as single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq). Here, we demonstrate that mitochondrial germline variants are an alternative, robust, and computationally efficient endogenous barcode for sample demultiplexing. We propose MitoSort, a tool that uses mitochondrial germline variants to assign cells to their donor origins and identify cross-genotype doublets in single-cell genomic datasets. We evaluate its performance by using in silico pooled mitochondrial scATAC-seq (mtscATAC-seq) libraries and experimentally multiplexed data with cell hashtags. MitoSort achieves high accuracy and efficiency in genotype clustering and doublet detection for mtscATAC-seq data, addressing the limitations of current computational techniques tailored for scRNA-seq data. Moreover, MitoSort exhibits versatility, and can be applied to various single-cell sequencing approaches beyond mtscATAC-seq provided that the mitochondrial variants are reliably detected. Furthermore, we demonstrate the application of MitoSort in a case study where B cells from eight donors were pooled and assayed by single-cell multi-omics sequencing. Altogether, our results demonstrate the accuracy and efficiency of MitoSort, which enables reliable sample demultiplexing in various single-cell genomic applications. MitoSort is available at https://github.com/tangzhj/MitoSort.

De novo detection of somatic variants in high-quality long-read single-cell RNA sequencing data

In cancer, genetic and transcriptomic variations generate clonal heterogeneity, leading to treatment resistance. Long-read single-cell RNA sequencing (LR scRNA-seq) has the potential to detect genetic and transcriptomic variations simultaneously. Here, we present LongSom, a computational workflow leveraging high-quality LR scRNA-seq data to call de novo somatic single-nucleotide variants (SNVs), including in mitochondria (mtSNVs), copy-number alterations (CNAs), and gene fusions, to reconstruct the tumor clonal heterogeneity. Before somatic variants calling, LongSom re-annotates marker gene based cell types using cell mutational profiles. LongSom distinguishes somatic SNVs from noise and germline polymorphisms by applying an extensive set of hard filters and statistical tests. Applying LongSom to human ovarian cancer samples, we detected clinically relevant somatic SNVs that were validated against matched DNA samples. Leveraging somatic SNVs and fusions, LongSom found subclones with different predicted treatment outcomes. In summary, LongSom enables de novo variant detection without the need for normal samples, facilitating the study of cancer evolution, clonal heterogeneity, and treatment resistance.

Identifying cancer cells from calling single-nucleotide variants in scRNA-seq data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) data are widely used to study cancer cell states and their heterogeneity. However, the tumour microenvironment is usually a mixture of healthy and cancerous cells and it can be difficult to fully separate these two populations based on transcriptomics alone. If available, somatic single-nucleotide variants (SNVs) observed in the scRNA-seq data could be used to identify the cancer population and match that information with the single cells' expression profile. However, calling somatic SNVs in scRNA-seq data is a challenging task, as most variants seen in the short-read data are not somatic, but can instead be germline variants, RNA edits or transcription, sequencing, or processing errors. In addition, only variants present in actively transcribed regions for each individual cell will be seen in the data.

RESULTS

To address these challenges, we develop CCLONE (Cancer Cell Labelling On Noisy Expression), an interpretable tool adapted to handle the uncertainty and sparsity of SNVs called from scRNA-seq data. CCLONE jointly identifies cancer clonal populations, and their associated variants. We apply CCLONE on two acute myeloid leukaemia datasets and one lung adenocarcinoma dataset and show that CCLONE captures both genetic clones and somatic events for multiple patients. These results show how CCLONE can be used to gather insight into the course of the disease and the origin of cancer cells in scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Source code is available at github.com/HaghverdiLab/CCLONE.

Inferring clonal somatic mutations directed by X chromosome inactivation status in single cells

Analysis of clonal dynamics in human tissues is enabled by somatic genetic variation. Here, we show that analysis of mitochondrial mutations in single cells is dramatically improved in females when using X chromosome inactivation to select informative clonal mutations. Applying this strategy to human peripheral mononuclear blood cells reveals clonal structures within T cells that otherwise are blurred by non-informative mutations, including the separation of gamma-delta T cells, suggesting this approach can be used to decipher clonal dynamics of cells in human tissues.

Robust analysis of allele-specific copy number alterations from scRNA-seq data with XClone

Somatic copy number alterations (CNAs) are major mutations that contribute to the development and progression of various cancers. Despite a few computational methods proposed to detect CNAs from single-cell transcriptomic data, the technical sparsity of such data makes it challenging to identify allele-specific CNAs, particularly in complex clonal structures. In this study, we present a statistical method, XClone, that strengthens the signals of read depth and allelic imbalance by effective smoothing on cell neighborhood and gene coordinate graphs to detect haplotype-aware CNAs from scRNA-seq data. By applying XClone to multiple datasets with challenging compositions, we demonstrated its ability to robustly detect different types of allele-specific CNAs and potentially indicate whole genome duplication, therefore enabling the discovery of corresponding subclones and the dissection of their phenotypic impacts.

The bioenergetic landscape of cancer

BACKGROUND

Bioenergetic remodeling of core energy metabolism is essential to the initiation, survival, and progression of cancer cells through exergonic supply of adenosine triphosphate (ATP) and metabolic intermediates, as well as control of redox homeostasis. Mitochondria are evolutionarily conserved organelles that mediate cell survival by conferring energetic plasticity and adaptive potential. Mitochondrial ATP synthesis is coupled to the oxidation of a variety of substrates generated through diverse metabolic pathways. As such, inhibition of the mitochondrial bioenergetic system by restricting metabolite availability, direct inhibition of the respiratory Complexes, altering organelle structure, or coupling efficiency may restrict carcinogenic potential and cancer progression.

SCOPE OF REVIEW

Here, we review the role of bioenergetics as the principal conductor of energetic functions and carcinogenesis while highlighting the therapeutic potential of targeting mitochondrial functions.

MAJOR CONCLUSIONS

Mitochondrial bioenergetics significantly contribute to cancer initiation and survival. As a result, therapies designed to limit oxidative efficiency may reduce tumor burden and enhance the efficacy of currently available antineoplastic agents.

Joint inference of cell lineage and mitochondrial evolution from single-cell sequencing data

MOTIVATION

Eukaryotic cells contain organelles called mitochondria that have their own genome. Most cells contain thousands of mitochondria which replicate, even in nondividing cells, by means of a relatively error-prone process resulting in somatic mutations in their genome. Because of the higher mutation rate compared to the nuclear genome, mitochondrial mutations have been used to track cellular lineage, particularly using single-cell sequencing that measures mitochondrial mutations in individual cells. However, existing methods to infer the cell lineage tree from mitochondrial mutations do not model "heteroplasmy," which is the presence of multiple mitochondrial clones with distinct sets of mutations in an individual cell. Single-cell sequencing data thus provide a mixture of the mitochondrial clones in individual cells, with the ancestral relationships between these clones described by a mitochondrial clone tree. While deconvolution of somatic mutations from a mixture of evolutionarily related genomes has been extensively studied in the context of bulk sequencing of cancer tumor samples, the problem of mitochondrial deconvolution has the additional constraint that the mitochondrial clone tree must be concordant with the cell lineage tree.

RESULTS

We formalize the problem of inferring a concordant pair of a mitochondrial clone tree and a cell lineage tree from single-cell sequencing data as the Nested Perfect Phylogeny Mixture (NPPM) problem. We derive a combinatorial characterization of the solutions to the NPPM problem, and formulate an algorithm, MERLIN, to solve this problem exactly using a mixed integer linear program. We show on simulated data that MERLIN outperforms existing methods that do not model mitochondrial heteroplasmy nor the concordance between the mitochondrial clone tree and the cell lineage tree. We use MERLIN to analyze single-cell whole-genome sequencing data of 5220 cells of a gastric cancer cell line and show that MERLIN infers a more biologically plausible cell lineage tree and mitochondrial clone tree compared to existing methods.

AVAILABILITY AND IMPLEMENTATION

https://github.com/raphael-group/MERLIN.

Single-cell mtDNA dynamics in tumors is driven by coregulation of nuclear and mitochondrial genomes

The extent of cell-to-cell variation in tumor mitochondrial DNA (mtDNA) copy number and genotype, and the phenotypic and evolutionary consequences of such variation, are poorly characterized. Here we use amplification-free single-cell whole-genome sequencing (Direct Library Prep (DLP+)) to simultaneously assay mtDNA copy number and nuclear DNA (nuDNA) in 72,275 single cells derived from immortalized cell lines, patient-derived xenografts and primary human tumors. Cells typically contained thousands of mtDNA copies, but variation in mtDNA copy number was extensive and strongly associated with cell size. Pervasive whole-genome doubling events in nuDNA associated with stoichiometrically balanced adaptations in mtDNA copy number, implying that mtDNA-to-nuDNA ratio, rather than mtDNA copy number itself, mediated downstream phenotypes. Finally, multimodal analysis of DLP+ and single-cell RNA sequencing identified both somatic loss-of-function and germline noncoding variants in mtDNA linked to heteroplasmy-dependent changes in mtDNA copy number and mitochondrial transcription, revealing phenotypic adaptations to disrupted nuclear/mitochondrial balance.

Gene-expression memory-based prediction of cell lineages from scRNA-seq datasets

Assigning single cell transcriptomes to cellular lineage trees by lineage tracing has transformed our understanding of differentiation during development, regeneration, and disease. However, lineage tracing is technically demanding, often restricted in time-resolution, and most scRNA-seq datasets are devoid of lineage information. Here we introduce Gene Expression Memory-based Lineage Inference (GEMLI), a computational tool allowing to robustly identify small to medium-sized cell lineages solely from scRNA-seq datasets. GEMLI allows to study heritable gene expression, to discriminate symmetric and asymmetric cell fate decisions and to reconstruct individual multicellular structures from pooled scRNA-seq datasets. In human breast cancer biopsies, GEMLI reveals previously unknown gene expression changes at the onset of cancer invasiveness. The universal applicability of GEMLI allows studying the role of small cell lineages in a wide range of physiological and pathological contexts, notably in vivo. GEMLI is available as an R package on GitHub ( https://github.com/UPSUTER/GEMLI ).

Deciphering cell states and genealogies of human haematopoiesis

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Single-cell lineage tracing with endogenous markers

Resolving lineage relationships between cells in an organism provides key insights into the fate of individual cells and drives a fundamental understanding of the process of development and disease. A recent rapid increase in experimental and computational advances for detecting naturally occurring somatic nuclear and mitochondrial mutation at single-cell resolution has expanded lineage tracing from model organisms to humans. This review discusses the advantages and challenges of experimental and computational techniques for cell lineage tracing using somatic mutation as endogenous DNA barcodes to decipher the relationships between cells during development and tumour evolution. We outlook the advantages of spatial clonal evolution analysis and single-cell lineage tracing using endogenous genetic markers.

MitoTracer facilitates the identification of informative mitochondrial mutations for precise lineage reconstruction

Mitochondrial (MT) mutations serve as natural genetic markers for inferring clonal relationships using single cell sequencing data. However, the fundamental challenge of MT mutation-based lineage tracing is automated identification of informative MT mutations. Here, we introduced an open-source computational algorithm called "MitoTracer", which accurately identified clonally informative MT mutations and inferred evolutionary lineage from scRNA-seq or scATAC-seq samples. We benchmarked MitoTracer using the ground-truth experimental lineage sequencing data and demonstrated its superior performance over the existing methods measured by high sensitivity and specificity. MitoTracer is compatible with multiple single cell sequencing platforms. Its application to a cancer evolution dataset revealed the genes related to primary BRAF-inhibitor resistance from scRNA-seq data of BRAF-mutated cancer cells. Overall, our work provided a valuable tool for capturing real informative MT mutations and tracing the lineages among cells.

Teaser

MitoTracer enables automatically and accurately discover informative mitochondrial mutations for lineage tracing.

Systematic investigation of mitochondrial transfer between cancer cells and T cells at single-cell resolution

Mitochondria (MT) participate in most metabolic activities of mammalian cells. A near-unidirectional mitochondrial transfer from T cells to cancer cells was recently observed to "metabolically empower" cancer cells while "depleting immune cells," providing new insights into tumor-T cell interaction and immune evasion. Here, we leverage single-cell RNA-seq technology and introduce MERCI, a statistical deconvolution method for tracing and quantifying mitochondrial trafficking between cancer and T cells. Through rigorous benchmarking and validation, MERCI accurately predicts the recipient cells and their relative mitochondrial compositions. Application of MERCI to human cancer samples identifies a reproducible MT transfer phenotype, with its signature genes involved in cytoskeleton remodeling, energy production, and TNF-α signaling pathways. Moreover, MT transfer is associated with increased cell cycle activity and poor clinical outcome across different cancer types. In summary, MERCI enables systematic investigation of an understudied aspect of tumor-T cell interactions that may lead to the development of therapeutic opportunities.

MitoTrace: A Computational Framework for Analyzing Mitochondrial Variation in Single-Cell RNA Sequencing Data

Genetic variation in the mitochondrial genome is linked to important biological functions and various human diseases. Recent progress in single-cell genomics has established single-cell RNA sequencing (scRNAseq) as a popular and powerful technique to profile transcriptomics at the cellular level. While most studies focus on deciphering gene expression, polymorphisms including mitochondrial variants can also be readily inferred from scRNAseq. However, limited attention has been paid to investigate the single-cell landscape of mitochondrial variants, despite the rapid accumulation of scRNAseq data in the community. In addition, a diploid context is assumed for most variant calling tools, which is not appropriate for mitochondrial heteroplasmies. Here, we introduce MitoTrace, an R package for the analysis of mitochondrial genetic variation in bulk and scRNAseq data. We applied MitoTrace to several publicly accessible data sets and demonstrated its ability to robustly recover genetic variants from scRNAseq data. We also validated the applicability of MitoTrace to scRNAseq data from diverse platforms. Overall, MitoTrace is a powerful and user-friendly tool to investigate mitochondrial variants from scRNAseq data.

Haplotype-aware analysis of somatic copy number variations from single-cell transcriptomes

Genome instability and aberrant alterations of transcriptional programs both play important roles in cancer. Single-cell RNA sequencing (scRNA-seq) has the potential to investigate both genetic and nongenetic sources of tumor heterogeneity in a single assay. Here we present a computational method, Numbat, that integrates haplotype information obtained from population-based phasing with allele and expression signals to enhance detection of copy number variations from scRNA-seq. Numbat exploits the evolutionary relationships between subclones to iteratively infer single-cell copy number profiles and tumor clonal phylogeny. Analysis of 22 tumor samples, including multiple myeloma, gastric, breast and thyroid cancers, shows that Numbat can reconstruct the tumor copy number profile and precisely identify malignant cells in the tumor microenvironment. We identify genetic subpopulations with transcriptional signatures relevant to tumor progression and therapy resistance. Numbat requires neither sample-matched DNA data nor a priori genotyping, and is applicable to a wide range of experimental settings and cancer types.

Mitochondrial DNA Mutations and Ageing

Mitochondria are subcellular organelles present in most eukaryotic cells which play a significant role in numerous aspects of cell biology. These include carbohydrate and fatty acid metabolism to generate cellular energy through oxidative phosphorylation, apoptosis, cell signalling, haem biosynthesis and reactive oxygen species production. Mitochondrial dysfunction is a feature of many human ageing tissues, and since the discovery that mitochondrial DNA mutations were a major underlying cause of changes in oxidative phosphorylation capacity, it has been proposed that they have a role in human ageing. However, there is still much debate on whether mitochondrial DNA mutations play a causal role in ageing or are simply a consequence of the ageing process. This chapter describes the structure of mammalian mitochondria, and the unique features of mitochondrial genetics, and reviews the current evidence surrounding the role of mitochondrial DNA mutations in the ageing process. It then focusses on more recent discoveries regarding the role of mitochondrial dysfunction in stem cell ageing and age-related inflammation.

Mitochondrial DNA is a major source of driver mutations in cancer

Mitochondrial DNA (mtDNA) mutations are among the most common genetic events in all tumors and directly impact metabolic homeostasis. Despite the central role mitochondria play in energy metabolism and cellular physiology, the role of mutations in the mitochondrial genomes of tumors has been contentious. Until recently, genomic and functional studies of mtDNA variants were impeded by a lack of adequate tumor mtDNA sequencing data and available methods for mitochondrial genome engineering. These barriers and a conceptual fog surrounding the functional impact of mtDNA mutations in tumors have begun to lift, revealing a path to understanding the role of this essential metabolic genome in cancer initiation and progression. Here we discuss the history, recent developments, and challenges that remain for mitochondrial oncogenetics as the impact of a major new class of cancer-associated mutations is unveiled.

Deciphering Innate Immune Cell-Tumor Microenvironment Crosstalk at a Single-Cell Level

The tumor microenvironment encompasses various innate immune cells which regulate tumor progression. Exploiting innate immune cells is a new frontier of cancer immunotherapy. However, the classical surface markers for cell-type classification cannot always well-conclude the phenotype, which will further hinge our understanding. The innate immune cells include dendritic cells, monocytes/macrophages, natural killer cells, and innate lymphoid cells. They play important roles in tumor growth and survival, in some cases promoting cancer, in other cases negating cancer. The precise characterization of innate immune cells at the single-cell level will boost the potential of cancer immunotherapy. With the development of single-cell RNA sequencing technology, the transcriptome of each cell in the tumor microenvironment can be dissected at a single-cell level, which paves a way for a better understanding of the cell type and its functions. Here, we summarize the subtypes and functions of innate immune cells in the tumor microenvironment based on recent literature on single-cell technology. We provide updates on recent achievements and prospects for how to exploit novel functions of tumor-associated innate immune cells and target them for cancer immunotherapy.