Mitochondrial variant enrichment from high-throughput single-cell RNA sequencing resolves clonal populations

Abstract

Reconstructing lineage relationships in complex tissues can reveal mechanisms underlying development and disease. Recent methods combine single-cell transcriptomics with mitochondrial DNA variant detection to establish lineage relationships in primary human cells, but are not scalable to interrogate complex tissues. To overcome this limitation, here we develop a technology for highconfidence detection of mitochondrial mutations from high-throughput single-cell RNA-sequencing. We use the new method to identify skewed immune cell expansions in primary human clonal hematopoiesis. Main textSingle-cell RNA-sequencing (scRNA-seq) enables the unbiased assessment of cell states in health and disease 1,2 . Combined acquisition of cell state and genetic information can provide additional insight, such as targeted enrichment of cancer driver mutations from single-cell transcriptomes 3,4 . Separately, combining scRNA-seq with genetic cell barcodes is a powerful method to reveal clonal relationships and evolutionary dynamics of cells within organisms 5,6 . However, this has largely been limited to experimental model systems that can be genetically manipulated to insert cell barcodes. To infer clonal dynamics in primary human cells, recent methods have detected and utilized mitochondrial DNA (mtDNA) mutations as naturally occurring genetic cell barcodes [7][8][9] . The combination of scRNA-seq with mtDNA mutation detection can inform clonal relationships with high confidence, but is currently restricted to expensive, low-throughput, full-length transcript sequencing technologies like SmartSeq2 7,10 . To enable the reconstruction of clonal relationships in complex human tissues, we developed a method that captures genetic variants from high-throughput scRNA-seq platforms: MAESTER, or Mitochondrial Alteration Enrichment from Singlecell Transcriptomes to Establish Relatedness (Figure 1A). MAESTER is compatible with the most common high-throughput scRNA-seq platforms, including 10x Genomics 3' protocols, Seq-Well S 3 , and Drop-seq (Supplemental Figures 1-3) 11,12 . An intermediate step in each of these platforms yields full length cDNA transcripts, from which we enrich all 15 mitochondrial transcripts using pools of primers, while maintaining cell-identifying barcodes (Figure 1B, Supplemental Figure 4). Standard next-generation sequencing with 250 bp reads is then used to obtain the sequence of the amplified mitochondrial transcripts (Figure 1A). We developed a computational toolkit to call mtDNA variants from MAESTER data, the Mitochondrial Alteration Enrichment and Genome Analysis Toolkit (maegatk, Supplemental Figure 5, Methods). Building on previous tools that we developed 8 for mtDNA variant detection from single-cell ATAC or SmartSeq2, maegatk specifically handles technical biases implicit in high-throughput transcriptomic libraries. Critically, maegatk leverages unique molecular identifiers (UMIs) to collapse multiple sequencing reads of the same starting transcript, creating a consensus call for every

Clonally resolved single-cell multi-omics identifies routes of cellular differentiation in acute myeloid leukemia

“…To cover mitochondrial genomes at full-length, we used PCR primers located near the 5’ end of mitochondrial genes, and subsequently performed an optimized tagmentation protocol (Figure S1a,b). Full length coverage was substantially improved compared to default 10x, and similar to a recent report that optimized mitochondrial coverage using a purely PCR-based strategy 14 (Figure 1b). Overall, 55-85% of these libraries mapped to the mitochondrial genome, allowing for a cost-effective deep sequencing of mitochondrial genomes.…”

“…To answer these questions, we have developed a new strategy to add clonal resolution to high-throughput (droplet-based) single-cell RNA-seq data that robustly works across many of the heterogeneous AML genotypes. Existing approaches use single nucleotide variant (SNVs) or mitochondrial SNVs (mtSNVs) as qualitative markers to identify healthy and malignant cells from scRNAseq data [8][9][10]13,14 . However, these measurements are noisy and biased by cellular identity impeding quantitative analysis.…”

enclone: precision clonotyping and analysis of immune receptors

“…We puzzled over myriad examples where we were unsure of the true answer and became confident that (lacking ground truth) we could not solve the problem by simulation. A way out could be to generate single cell data capturing both VDJ and mitochondrial mRNA [Miller 2022]. One could also engineer cell-intrinsic barcodes [Bowling 2020: Leeper 2021], but not in humans.…”

Functional antibodies exhibit light chain coherence

“…However, in the absence of a sufficiently large dataset with multiple antigen specificities using cells from multiple humans or donors, it is currently impossible to assess the validity of any functional grouping scheme. Innovative methods such as mitochondrial lineage tracing could perhaps be used to validate computed clonotypes 26 .…”