Resolving lineage relationships between cells is necessary to understand the fundamental mechanisms underlying normal development and the progression of disease. In recent years, new methods have emerged to enable cell lineage tracking with increasing resolution, leading to substantial biological insights (Woodworth et al., 2017). Specifically, genome editing of reporter constructs via CRISPR-Cas9 allowed synthetic reconstruction of cell lineage relationships in model organisms, and has been coupled with transcriptome profiling to inform cell fates (Spanjaard et al., 2018). These prospective ‘mutate-and-record’ methods provide powerful tools to resolve the developmental origin of cells in genetically engineered cells and organisms, but cannot be utilized in living humans, archival clinical samples, or any wild type organism (Woodworth et al., 2017). Given these limitations, retrospective lineage tracing using endogenous genetic markers is an alternative solution. Recent advances in sequencing enable naturally occurring somatic mutations to be used as lineage markers, which usually required single-cell genome sequencing to capture the sparse genetic information (Lee-Six et al., 2018; Wang et al., 2014). Regions with high mutation rates, such as microsatellite repeats, retrotransposons, and copy-number variants, has been used to resolve the lineage relationship for normal or cancerous tissue samples (Evrony et al., 2015; Biezuner et al., 2016). These methods reduce the cost of whole genome sequencing, but still lack information on cell phenotypes.

Simultaneous measurement of the lineage relationship and cell fates is ultimately required to address many biomedical questions. Here we describe EMBLEM (Epigenome and Mitochondrial Barcode of Lineage from Endogenous Mutations), a strategy to track cell lineage using endogenous mitochondrial DNA variants in ATAC-seq data. The end result of EMBLEM is single-cell lineage information and rich global epigenomic profile from the same individual cells (Figure 1A and Figure 1—figure supplement 1).

Figure 1 with 5 supplements see all

Download asset Open asset

We illustrate the utility of EMBLEM in human blood progenitor cells to clarify the process of pre-leukemic clonal evolution and the emerging biology of clonal hematopoiesis.

Assay of Transposase-Accessible Chromatin by sequencing (ATAC-seq) is a sensitive method used to study chromatin accessibility profiles in diverse cell types and organisms (Buenrostro et al., 2013). During DNA transposition and amplification in cells, mitochondrial DNA is also amplified at the same time (Figure 1A). Mitochondrial DNA (mtDNA) is a ~ 16 kb circular genome with ~10 fold higher mutation rate compared to the nuclear genome. Hence, mtDNA incrementally accumulates unique, irreversible genetic mutations that are passed on to daughter cells even in healthy humans and may be used for lineage tracing (Morris et al., 2017), (Fellous et al., 2009). The majority of somatic mtDNA mutations are noncoding and thought to be passenger (Ju et al., 2014). Importantly, the number of mitochondria (and therefore mtDNA) range from several hundreds to >10,000 per cell in different cell types, facilitating robust mtDNA analysis even from a single cell.

We first observed that ATAC-seq effectively enriches for mtDNA. While mtDNA is present in many kinds of DNA sequence libraries, it is substantially enriched in ATAC-seq libraries due to the fact that mtDNA is not chromatinized and is therefore highly accessible (Supplementary file 1). ATAC-seq enables a 17-fold or greater enrichment of mtDNA compared to exome sequencing or whole genome sequencing in GM12878 human B cells (Figure 1B), leading to an average ~18,000X coverage of mtDNA (Figure 1—figure supplement 2A). With this coverage, we detected 27 mitochondrial variants from GM12878 cells (Figure 1C). 13 of these variants have a variant allele frequency (VAF) greater than 90%, which are known as homoplasmic variants (Figure 1—figure supplement 2B). We also detected 14 low frequency mitochondrial DNA variants, with VAFs ranging from 0.1% to 24% (Figure 1C and Figure 1—figure supplement 2C). Similar results for mtDNA enrichment were observed in human K562 cells (Figure 1—figure supplement 3, Supplementary file 1).

The VAF from bulk ATAC-seq data represents the average of the allele frequencies of the cell population. A 25% VAF may arise from 25% of cells in the population with a homoplasmic variant, or alternatively arise from 100% of cells all having a quarter of their mitochondria with the variant allele (Figure 1D). To distinguish between these two models, we analyzed single-cell ATAC-seq data from GM12878. For four mtDNA variants (VAF between 0.5% ~ 24% at population level), we find that a mixture of both models is in action for different variants (Figure 1E). For instance, mtDNA mutation 3082 is widely spread among single cells, but at low frequency per cell. Because it is extremely unlikely (see Materials and method) that the identical mutation arose independently in every single cell, cells sharing the same mitochondrial mutations are inferred to have descended from the same ancestral cell. These results suggest that even low frequency heteroplasmic mtDNA mutations can be exploited for lineage tracing.

To prove the principle that somatic mitochondrial mutations can track cells from the same ancestor and to quantify the power of lineage mapping, we next applied EMBLEM to primary blood cells from patients with acute myeloid leukemia (AML). Human AML is organized as a hierarchy: a hematopoietic stem cell first acquires an initiating mutation in one of a number of chromatin modifier genes, previously termed as ‘pre-leukemic’ hematopoietic stem cell (pHSC) (Jan et al., 2012; Corces-Zimmerman et al., 2014). pHSCs are functionally normal and are not able to transplant AML, but upon accumulation of additional mutations, they give rise to leukemic stem cells (LSCs) that are able to self-renew and recapitulate AML disease upon transplantation (Jan et al., 2012; Thomas and Majeti, 2017). Finally, LSCs give rise to the bulk leukemic blast cells in AML (Thomas and Majeti, 2017). Targeted exome sequencing in these samples have identified somatic mutations in tumor suppressor genes and oncogenes that link the lineage relationship of pHSCs, LSCs and blasts, providing the ground truth for our analyses (Corces-Zimmerman et al., 2014).

We applied EMBLEM to the ATAC-seq profiles of FACS-purified LSCs and leukemic blasts first (Supplementary file 1). Using high-confidence mtDNA mutations, detected both from bulk ATAC-seq and single cell ATAC-seq, we found the LSC and blast populations not only shared the same heteroplasmic variants, but also showed similar distribution and allele frequency at the cellular level (Figure 1—figure supplement 4). These results indicate the two populations are identical at the genetic level, but divergent at the epigenomic level, consistent with previous studies (Corces et al., 2016; Schep et al., 2017). In patient SU353, we identified four diagnostic mtDNA mutations in the same cell (Figure 1F), which indicates these four mitochondrial variants already co-existed in the ancestral cell (see Materials and method). With the assumption that all these LSCs and blasts are clonal, we further quantified the detection rate of each mtDNA variant as a function of allele frequency and sequencing depth (Figure 1G). We found that when a single variant allele has a frequency greater than 20%, the detection rate can be up to 90% with >20X coverage (e.g. site 6776). In contrast, when the variant allele has a frequency lower than 1%, the detection rate drops to 20% when the coverage is below 100X (e.g. site 6705). While high drop-out rate is a common challenge for single-cell technologies (Kharchenko et al., 2014), computational imputation of the missing information from single cell data can address this problem (Zhang and Zhang, 2018). When multiple mtDNA variants are co-detected in multiple single cells, we can infer their origin and linkage in the ancestral cell (see Materials and method). Thus, cells containing any one of these variants will still inform their origin from the same lineage. With any combination of the four variants, 90% (sensitivity) of the cells can be unambiguously assigned to the correct lineage with just 20x mtDNA coverage (Figure 1H). Furthermore, two mtDNA mutations identified in other cells (e.g. pHSC specific site 2967,6268) were never detected (false positive = 0) in LSCs and blasts (Figure 1—figure supplement 4), showing a high specificity of the method. Similar performance of single cell lineage tracing for another patient (SU070) is shown in Figure 1—figure supplement 5. These results demonstrate that somatic DNA mutations in the mitochondrial genome are a powerful endogenous marker to identify clonal cell populations.

To expand on these findings to additional different cell lineages, we applied EMBLEM to bulk ATAC-seq data from sorted blood cells from healthy human donors and patients with AML(Supplementary file 1) (Corces et al., 2016). We identified heteroplasmic mtDNA mutations in multiple cell populations of primary blood cells from healthy donors and all AML patients (Figure 2—figure supplement 1A, Supplementary file 2). The heteroplasmic mtDNA mutations showed a similar mutant spectrum as observed by previous studies using cancer genomic data (Figure 2—figure supplement 1B and C) (Ju et al., 2014).

Furthermore, EMBLEM, not only confirmed the previous lineage hierarchy of AML, but also extended the previous model of pHSC heterogeneity (Figure 2A). In the AML cases with LSCs sequenced by ATAC-seq, the LSCs and their corresponding leukemic blasts have nearly identical heteroplasmic mtDNA mutations (Figure 2B–C and Figure 2—figure supplement 1D), suggesting a direct lineage relationship and short generation history between LSCs and blasts. We then examined whether any of the mtDNA variants present in LSCs can be seen in the pHSCs, where the first leukemia-associated protein-coding mutations have already occurred in functional normal hematopoietic stem cells (Corces-Zimmerman et al., 2014; Corces et al., 2016). We detected blast-associated mtDNA mutations in pHSCs in all 11 cases. Interestingly, we also detected additional heteroplasmic mtDNA mutations present specifically in pHSCs (Figure 2C). In the 11 cases we investigated, seven cases have pHSC-unique heteroplasmic mtDNA mutations (Figure 2—figure supplement 1D and E), a previously unrecognized level of pHSC heterogeneity. pHSCs are capable of long-term self-renewal and possess a clonal growth advantage, allowing them to clonally outcompete normal HSCs. Indeed, the clonal frequency of pHSCs is a poor prognostic factor for overall survival in AML (Corces et al., 2016). Our discovery of pHSCs with distinct heteroplasmic mtDNA mutations suggests the existence of multiple distinct sub-clones of pHSCs in AML patients.

Figure 2 with 4 supplements see all

Download asset Open asset

To validate the heterogeneity of pHSCs inferred from EMBLEM of bulk cell populations, we performed single-cell ATAC-seq of HSCs from AML patient SU353, which exhibited both a high burden of pre-leukemic somatic coding gene mutations and high frequency of pHSC-specific heteroplasmic mtDNA mutations (Corces et al., 2016). We identified the heteroplasmic mtDNA variants from each single cell, which separated the HSCs into three lineages: Two clonal subpopulations termed ‘clone 1’ (18 cells) and ‘clone 2’ (104 cells), and a third population with no mtDNA variants despite sufficient mtDNA coverage (pHSC with WT mtDNA, 31 cells) (Figure 2D). Notably, clone 2 possessed pHSC-specific mtDNA mutations, while clone 1 possessed mtDNA mutations shared with LSCs, indicating clone1 is the lineage precursor of AML. These results confirm that multiple pHSC clones arise in AML patients, and one subclone eventually evolved to become the LSC (Figure 2—figure supplement 2A).

Finally, we related the clonotype of pHSCs to their single-cell chromatin accessibility profiles. We interrogated the patterns of active DNA elements and enriched transcription factor motifs in sequential stages of AML development from the same patient, and contrasted with HSCs from normal donors using ChromVAR (Schep et al., 2017) (Figure 2E and Figure 2—figure supplement 2B). The chromatin accessibility profiles of pHSCs are more similar to HSCs than to LSCs or leukemic blasts. The greatest deviation between HSC and other cell types occurred at DNA binding motifs of the transcription factor Jun/Fos, a known key regulator of HSC biology (Santaguida et al., 2009) (Figure 2F). Furthermore, the three lineages of pHSCs revealed by mtDNA mutations also showed distinctive chromatin profiles (Figure 2G). Clone 1 pHSC, which gives rise to the LSC and AML leukemia, is already more similar to LSCs and blasts in its chromatin accessibility. In contrast, clone 2 that comprises the larger fraction of pHSCs exhibited variable chromatin profiles at the single-cell level that spanned the range of normal HSCs, pHSC with WT mtDNA(WT cells), is also diverged from normal HSCs. Thus, both lineage tracing and single cell epigenomic states indicate clone 1 as the original stem cell of the AML in patient SU353. Supervised comparison of the chromatin accessibility profiles among these clonal sub-populations further identified distinct and significantly enriched transcription factor motifs (Figure 2H and Figure 2—figure supplement 2C–E). These results indicate the heterogeneity of HSCs from AML patients both on a genetic and epigenomic level.

We present a computational strategy to combine cell lineages tracing by endogenous mtDNA mutations and chromatin accessibility profiling in the same cell using single-cell ATAC-seq data. This approach is applicable to any eukaryote, does not require genetic engineering or genome editing, and is cost effective as the lineage information comes ‘for free’ on top of epigenomic insights. The relative merits of mtDNA vs other genetic markers for lineage tracing are outlined in Supplementary file 3. An important advantage of EMBLEM is that we enable clonotype tracing in existing ATAC-seq data sets and hierarchical lineage construction from ATAC-seq that thousands of labs have already generated. All future ATAC-seq data acquired for other inquiries will also have the benefit of lineage information. EMBLEM may also be extended to other single cell technologies, in which mtDNA is sequenced. We show that EMBLEM is successful even with low frequency heteroplasmic mutations, detection of rare clones in a population, and authentic clinical samples. With advances in the throughput and depth of single-cell genomic technologies, we believe EMBLEM may be a powerful tool to bring insight for many biomedical questions, including development, regeneration, immunity, and cancer with integration of genotype and phenotype information from the same cell. During revision of this work, Ludwig et al. (2019) reported the feasibility of using mtDNA and single cell genomics for lineage tracing, which independently validates the potentially broad utility of this approach.

Although powerful and broadly applicable, mtDNA lineage tracing also has its limitations. One limitation of this method is absence of mtDNA mutations in cells and tissues of embryos and young animals, which precluded us from applying EMBLEM to published scATAC-seq data of early animal development. Moreover, the possibility of selective mitochondrial inheritance or intercellular mitochondria transfer may affect the accuracy of inferred lineages (Mishra and Chan, 2014; Moschoi et al., 2016; Marlein et al., 2017; Hayakawa et al., 2016). On the other hand, asymmetric transmission of mitochondria would not necessarily affect cellular lineage tracing, as long as the variant alleles are randomly segregated. Using scATAC-seq data from a mixing experiment with human and mouse cells (Satpathy et al., 2018), we found species-specific mtDNA always paired with species-specific nuclear genomic DNA (Figure 2—figure supplement 4). These results suggest that mitochondrial horizontal transfer is not a confounder of our study and does not universally occur between cells. The aforementioned two scenarios reflect the potential uncoupling of nuclear and mitochondrial genomes, which would be of interest to investigate by EMBLEM in combination with other nucDNA tracing methods.

mtDNA lineage tracing produced new insights concerning the pHSC, the human hematopoietic stem cell that suffers the first oncogenic mutation in AML evolution. Our results add to the evidence that the pHSC population is heterogeneous, with evidence of multiple mtDNA clones. Unexpectedly, the pHSC lineage that gives rise to the subsequent acute myeloid leukemia is not the lineage with the best competitive potential among pHSCs, as the leukemogenic lineage is often in the minority. pHSC burden is a strong poor prognostic predictor of AML survival (Corces et al., 2016). It is widely believed that the association between high pHSC burden and poor AML patient prognosis reflected the enhanced self-renewal and competitive ability of the mutant pHSC. Our analysis suggests that high pHSC burden may reflect the diversity of pHSCs or the underlying mutational processes. These alternative interpretations of the link between pHSC burden and poor clinical prognosis should be addressed in future studies.

Sequencing data have been deposited in GEO under accession codes GSE122576 and GSE122577.

1. 1. Xu J

   (2018) NCBI Gene Expression Omnibus

   ID GSE122576. Single cell lineage tracing by endogenous mitochondrial DNA mutations in ATAC-seq data.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122576
2. 1. Xu J
   2. Chang HY

   (2018) NCBI Gene Expression Omnibus

   ID GSE122577. Single cell lineage tracing by endogenous mitochondrial DNA mutations in ATAC-seq data.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE122577

1. 1. Buenrostro JD
   2. Giresi PG
   3. Zaba LC
   4. Chang HY
   5. Greenleaf WJ

   (2013) NCBI Gene Expression Omnibus

   ID GSE47753. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE47753
2. 1. Buenrostro JD

   (2015) NCBI Gene Expression Omnibus

   ID GSE65360. Single-cell chromatin accessibility data using scATAC-seq.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE65360
3. 1. Buenrostro JD

   (2016) NCBI Gene Expression Omnibus

   ID GSE74912. ATAC-seq data.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE74912
4. 1. Buenrostro JD

   (2016) NCBI Gene Expression Omnibus

   ID GSE74310. Single-cell chromatin accessibility data using scATAC-seq.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE74310
5. 1. Buenrostro JD

   (2018) NCBI Gene Expression Omnibus

   ID GSE96772. Single-cell epigenomics maps the continuous regulatory landscape of human hematopoietic differentiation.

   https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE96772

1. 1. Biezuner T
   2. Spiro A
   3. Raz O
   4. Amir S
   5. Milo L
   6. Adar R
   7. Chapal-Ilani N
   8. Berman V
   9. Fried Y
   10. Ainbinder E
   11. Cohen G
   12. Barr HM
   13. Halaban R
   14. Shapiro E

   (2016) A generic, cost-effective, and scalable cell lineage analysis platform

   Genome Research 26:1588–1599.

   https://doi.org/10.1101/gr.202903.115

   • PubMed
   • Google Scholar
2. 1. Buenrostro JD
   2. Giresi PG
   3. Zaba LC
   4. Chang HY
   5. Greenleaf WJ

   (2013) Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position

   Nature Methods 10:1213–1218.

   https://doi.org/10.1038/nmeth.2688

   • PubMed
   • Google Scholar
3. 1. Buenrostro JD
   2. Wu B
   3. Litzenburger UM
   4. Ruff D
   5. Gonzales ML
   6. Snyder MP
   7. Chang HY
   8. Greenleaf WJ

   (2015) Single-cell chromatin accessibility reveals principles of regulatory variation

   Nature 523:486–490.

   https://doi.org/10.1038/nature14590

   • PubMed
   • Google Scholar
4. 1. Buenrostro JD
   2. Corces MR
   3. Lareau CA
   4. Wu B
   5. Schep AN
   6. Aryee MJ
   7. Majeti R
   8. Chang HY
   9. Greenleaf WJ

   (2018) Integrated Single-Cell analysis maps the continuous regulatory landscape of human hematopoietic differentiation

   Cell 173:1535–1548.

   https://doi.org/10.1016/j.cell.2018.03.074

   • PubMed
   • Google Scholar
5. 1. Chen L
   2. Liu P
   3. Evans TC
   4. Ettwiller LM

   (2017) DNA damage is a pervasive cause of sequencing errors, directly confounding variant identification

   Science 355:752–756.

   https://doi.org/10.1126/science.aai8690

   • PubMed
   • Google Scholar
6. 1. Coller HA
   2. Khrapko K
   3. Bodyak ND
   4. Nekhaeva E
   5. Herrero-Jimenez P
   6. Thilly WG

   (2001) High frequency of homoplasmic mitochondrial DNA mutations in human tumors can be explained without selection

   Nature Genetics 28:147–150.

   https://doi.org/10.1038/88859

   • PubMed
   • Google Scholar
7. 1. Conway JR
   2. Lex A
   3. Gehlenborg N

   (2017) UpSetR: an R package for the visualization of intersecting sets and their properties

   Bioinformatics 33:2938–2940.

   https://doi.org/10.1093/bioinformatics/btx364

   • PubMed
   • Google Scholar
8. 1. Corces MR
   2. Buenrostro JD
   3. Wu B
   4. Greenside PG
   5. Chan SM
   6. Koenig JL
   7. Snyder MP
   8. Pritchard JK
   9. Kundaje A
   10. Greenleaf WJ
   11. Majeti R
   12. Chang HY

   (2016) Lineage-specific and single-cell chromatin accessibility charts human hematopoiesis and leukemia evolution

   Nature Genetics 48:1193–1203.

   https://doi.org/10.1038/ng.3646

   • PubMed
   • Google Scholar
9. 1. Corces-Zimmerman MR
   2. Hong WJ
   3. Weissman IL
   4. Medeiros BC
   5. Majeti R

   (2014) Preleukemic mutations in human acute myeloid leukemia affect epigenetic regulators and persist in remission

   PNAS 111:2548–2553.

   https://doi.org/10.1073/pnas.1324297111

   • PubMed
   • Google Scholar
10. 1. Evrony GD
    2. Lee E
    3. Mehta BK
    4. Benjamini Y
    5. Johnson RM
    6. Cai X
    7. Yang L
    8. Haseley P
    9. Lehmann HS
    10. Park PJ
    11. Walsh CA

    (2015) Cell lineage analysis in human brain using endogenous retroelements

    Neuron 85:49–59.

    https://doi.org/10.1016/j.neuron.2014.12.028

    • PubMed
    • Google Scholar
11. 1. Fellous TG
    2. McDonald SA
    3. Burkert J
    4. Humphries A
    5. Islam S
    6. De-Alwis NM
    7. Gutierrez-Gonzalez L
    8. Tadrous PJ
    9. Elia G
    10. Kocher HM
    11. Bhattacharya S
    12. Mears L
    13. El-Bahrawy M
    14. Turnbull DM
    15. Taylor RW
    16. Greaves LC
    17. Chinnery PF
    18. Day CP
    19. Wright NA
    20. Alison MR

    (2009) A methodological approach to tracing cell lineage in human epithelial tissues

    Stem Cells 27:1410–1420.

    https://doi.org/10.1002/stem.67

    • PubMed
    • Google Scholar
12. 1. Hayakawa K
    2. Esposito E
    3. Wang X
    4. Terasaki Y
    5. Liu Y
    6. Xing C
    7. Ji X
    8. Lo EH

    (2016) Transfer of mitochondria from astrocytes to neurons after stroke

    Nature 535:551–555.

    https://doi.org/10.1038/nature18928

    • PubMed
    • Google Scholar
13. 1. Jan M
    2. Snyder TM
    3. Corces-Zimmerman MR
    4. Vyas P
    5. Weissman IL
    6. Quake SR
    7. Majeti R

    (2012) Clonal evolution of preleukemic hematopoietic stem cells precedes human acute myeloid leukemia

    Science Translational Medicine 4:ra118.

    https://doi.org/10.1126/scitranslmed.3004315

    • Google Scholar
14. 1. Ju YS
    2. Alexandrov LB
    3. Gerstung M
    4. Martincorena I
    5. Nik-Zainal S
    6. Ramakrishna M
    7. Davies HR
    8. Papaemmanuil E
    9. Gundem G
    10. Shlien A
    11. Bolli N
    12. Behjati S
    13. Tarpey PS
    14. Nangalia J
    15. Massie CE
    16. Butler AP
    17. Teague JW
    18. Vassiliou GS
    19. Green AR
    20. Du MQ
    21. Unnikrishnan A
    22. Pimanda JE
    23. Teh BT
    24. Munshi N
    25. Greaves M
    26. Vyas P
    27. El-Naggar AK
    28. Santarius T
    29. Collins VP
    30. Grundy R
    31. Taylor JA
    32. Hayes DN
    33. Malkin D
    34. Foster CS
    35. Warren AY
    36. Whitaker HC
    37. Brewer D
    38. Eeles R
    39. Cooper C
    40. Neal D
    41. Visakorpi T
    42. Isaacs WB
    43. Bova GS
    44. Flanagan AM
    45. Futreal PA
    46. Lynch AG
    47. Chinnery PF
    48. McDermott U
    49. Stratton MR
    50. Campbell PJ
    51. ICGC Breast Cancer Group
    52. ICGC Chronic Myeloid Disorders Group
    53. ICGC Prostate Cancer Group

    (2014) Origins and functional consequences of somatic mitochondrial DNA mutations in human cancer

    eLife 3:e02935.

    https://doi.org/10.7554/eLife.02935

    • Google Scholar
15. 1. Kharchenko PV
    2. Silberstein L
    3. Scadden DT

    (2014) Bayesian approach to single-cell differential expression analysis

    Nature Methods 11:740–742.

    https://doi.org/10.1038/nmeth.2967

    • PubMed
    • Google Scholar
16. 1. Koboldt DC
    2. Zhang Q
    3. Larson DE
    4. Shen D
    5. McLellan MD
    6. Lin L
    7. Miller CA
    8. Mardis ER
    9. Ding L
    10. Wilson RK

    (2012) VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing

    Genome Research 22:568–576.

    https://doi.org/10.1101/gr.129684.111

    • PubMed
    • Google Scholar
17. 1. Lee-Six H
    2. Øbro NF
    3. Shepherd MS
    4. Grossmann S
    5. Dawson K
    6. Belmonte M
    7. Osborne RJ
    8. Huntly BJP
    9. Martincorena I
    10. Anderson E
    11. O'Neill L
    12. Stratton MR
    13. Laurenti E
    14. Green AR
    15. Kent DG
    16. Campbell PJ

    (2018) Population dynamics of normal human blood inferred from somatic mutations

    Nature 561:473–478.

    https://doi.org/10.1038/s41586-018-0497-0

    • PubMed
    • Google Scholar
18. 1. Li H
    2. Handsaker B
    3. Wysoker A
    4. Fennell T
    5. Ruan J
    6. Homer N
    7. Marth G
    8. Abecasis G
    9. Durbin R
    10. 1000 Genome Project Data Processing Subgroup

    (2009b) The sequence alignment/Map format and SAMtools

    Bioinformatics 25:2078–2079.

    https://doi.org/10.1093/bioinformatics/btp352

    • PubMed
    • Google Scholar
19. 1. Li H
    2. Durbin R

    (2009a) Fast and accurate short read alignment with Burrows-Wheeler transform

    Bioinformatics 25:1754–1760.

    https://doi.org/10.1093/bioinformatics/btp324

    • PubMed
    • Google Scholar
20. 1. Ludwig LS
    2. Lareau CA
    3. Ulirsch JC
    4. Christian E
    5. Muus C
    6. Li LH
    7. Pelka K
    8. Ge W
    9. Oren Y
    10. Brack A
    11. Law T
    12. Rodman C
    13. Chen JH
    14. Boland GM
    15. Hacohen N
    16. Rozenblatt-Rosen O
    17. Aryee MJ
    18. Buenrostro JD
    19. Regev A
    20. Sankaran VG

    (2019) Lineage tracing in humans enabled by mitochondrial mutations and Single-Cell genomics

    Cell 176:1325–1339.

    https://doi.org/10.1016/j.cell.2019.01.022

    • PubMed
    • Google Scholar
21. 1. Marlein CR
    2. Zaitseva L
    3. Piddock RE
    4. Robinson SD
    5. Edwards DR
    6. Shafat MS
    7. Zhou Z
    8. Lawes M
    9. Bowles KM
    10. Rushworth SA

    (2017) NADPH oxidase-2 derived superoxide drives mitochondrial transfer from bone marrow stromal cells to leukemic blasts

    Blood 130:blood-2017-03-772939.

    https://doi.org/10.1182/blood-2017-03-772939

    • PubMed
    • Google Scholar
22. 1. McKenna A
    2. Hanna M
    3. Banks E
    4. Sivachenko A
    5. Cibulskis K
    6. Kernytsky A
    7. Garimella K
    8. Altshuler D
    9. Gabriel S
    10. Daly M
    11. DePristo MA

    (2010) The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data

    Genome Research 20:1297–1303.

    https://doi.org/10.1101/gr.107524.110

    • PubMed
    • Google Scholar
23. 1. Miller FJ
    2. Rosenfeldt FL
    3. Zhang C
    4. Linnane AW
    5. Nagley P

    (2003) Precise determination of mitochondrial DNA copy number in human skeletal and cardiac muscle by a PCR-based assay: lack of change of copy number with age

    Nucleic Acids Research 31:e61.

    https://doi.org/10.1093/nar/gng060

    • PubMed
    • Google Scholar
24. 1. Mishra P
    2. Chan DC

    (2014) Mitochondrial dynamics and inheritance during cell division, development and disease

    Nature Reviews Molecular Cell Biology 15:634–646.

    https://doi.org/10.1038/nrm3877

    • PubMed
    • Google Scholar
25. 1. Morris J
    2. Na YJ
    3. Zhu H
    4. Lee JH
    5. Giang H
    6. Ulyanova AV
    7. Baltuch GH
    8. Brem S
    9. Chen HI
    10. Kung DK
    11. Lucas TH
    12. O'Rourke DM
    13. Wolf JA
    14. Grady MS
    15. Sul JY
    16. Kim J
    17. Eberwine J

    (2017) Pervasive within-Mitochondrion Single-Nucleotide variant heteroplasmy as revealed by Single-Mitochondrion sequencing

    Cell Reports 21:2706–2713.

    https://doi.org/10.1016/j.celrep.2017.11.031

    • PubMed
    • Google Scholar
26. 1. Moschoi R
    2. Imbert V
    3. Nebout M
    4. Chiche J
    5. Mary D
    6. Prebet T
    7. Saland E
    8. Castellano R
    9. Pouyet L
    10. Collette Y
    11. Vey N
    12. Chabannon C
    13. Recher C
    14. Sarry JE
    15. Alcor D
    16. Peyron JF
    17. Griessinger E

    (2016) Protective mitochondrial transfer from bone marrow stromal cells to acute myeloid leukemic cells during chemotherapy

    Blood 128:253–264.

    https://doi.org/10.1182/blood-2015-07-655860

    • PubMed
    • Google Scholar
27. 1. Santaguida M
    2. Schepers K
    3. King B
    4. Sabnis AJ
    5. Forsberg EC
    6. Attema JL
    7. Braun BS
    8. Passegué E

    (2009) JunB protects against myeloid malignancies by limiting hematopoietic stem cell proliferation and differentiation without affecting self-renewal

    Cancer Cell 15:341–352.

    https://doi.org/10.1016/j.ccr.2009.02.016

    • PubMed
    • Google Scholar
28. 1. Satpathy AT
    2. Saligrama N
    3. Buenrostro JD
    4. Wei Y
    5. Wu B
    6. Rubin AJ
    7. Granja JM
    8. Lareau CA
    9. Li R
    10. Qi Y
    11. Parker KR
    12. Mumbach MR
    13. Serratelli WS
    14. Gennert DG
    15. Schep AN
    16. Corces MR
    17. Khodadoust MS
    18. Kim YH
    19. Khavari PA
    20. Greenleaf WJ
    21. Davis MM
    22. Chang HY

    (2018) Transcript-indexed ATAC-seq for precision immune profiling

    Nature Medicine 24:580–590.

    https://doi.org/10.1038/s41591-018-0008-8

    • PubMed
    • Google Scholar
29. 1. Schep AN
    2. Wu B
    3. Buenrostro JD
    4. Greenleaf WJ

    (2017) chromVAR: inferring transcription-factor-associated accessibility from single-cell epigenomic data

    Nature Methods 14:975–978.

    https://doi.org/10.1038/nmeth.4401

    • PubMed
    • Google Scholar
30. 1. Spanjaard B
    2. Hu B
    3. Mitic N
    4. Olivares-Chauvet P
    5. Janjuha S
    6. Ninov N
    7. Junker JP

    (2018) Simultaneous lineage tracing and cell-type identification using CRISPR-Cas9-induced genetic scars

    Nature Biotechnology 36:469–473.

    https://doi.org/10.1038/nbt.4124

    • PubMed
    • Google Scholar
31. 1. Thomas D
    2. Majeti R

    (2017) Biology and relevance of human acute myeloid leukemia stem cells

    Blood 129:1577–1585.

    https://doi.org/10.1182/blood-2016-10-696054

    • Google Scholar
32. 1. Wang Y
    2. Waters J
    3. Leung ML
    4. Unruh A
    5. Roh W
    6. Shi X
    7. Chen K
    8. Scheet P
    9. Vattathil S
    10. Liang H
    11. Multani A
    12. Zhang H
    13. Zhao R
    14. Michor F
    15. Meric-Bernstam F
    16. Navin NE

    (2014) Clonal evolution in breast cancer revealed by single nucleus genome sequencing

    Nature 512:155–160.

    https://doi.org/10.1038/nature13600

    • PubMed
    • Google Scholar
33. 1. Woodworth MB
    2. Girskis KM
    3. Walsh CA

    (2017) Building a lineage from single cells: genetic techniques for cell lineage tracking

    Nature Reviews Genetics 18:230–244.

    https://doi.org/10.1038/nrg.2016.159

    • PubMed
    • Google Scholar
34. Software

    1. Xu J

    (2019) ATAC_mito_sc

    Github.

    https://github.com/ChangLab/ATAC_mito_sc
35. 1. Zhang L
    2. Zhang S

    (2018) Comparison of computational methods for imputing single-cell RNA-sequencing data

    IEEE/ACM Transactions on Computational Biology and Bioinformatics 1:.

    https://doi.org/10.1109/TCBB.2018.2848633

    • Google Scholar

Author details

1. Jin Xu

   1. Center for Personal Dynamic Regulomes, Stanford, United States
   2. Department of Dermatology, Stanford University School of Medicine, Stanford, United States
   3. Department of Genetics, Stanford University School of Medicine, Stanford, United States

   Contribution

   Conceptualization, Data curation, Software, Formal analysis, Investigation, Methodology, Writing—original draft, Writing—review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-0944-9835

2. Kevin Nuno

   1. Institute for Stem Cell Biology and Regenerative Medicine, Stanford University School of Medicine, Stanford, United States
   2. Division of Hematology, Department of Medicine, Stanford University School of Medicine, Stanford, United States

   Contribution

   Resources, Data curation, Investigation, Methodology, Writing—review and editing

   Competing interests

   No competing interests declared

3. Ulrike M Litzenburger

   1. Center for Personal Dynamic Regulomes, Stanford, United States
   2. Department of Dermatology, Stanford University School of Medicine, Stanford, United States
   3. Department of Genetics, Stanford University School of Medicine, Stanford, United States

   Contribution

   Resources, Data curation, Investigation, Methodology

   Competing interests

   No competing interests declared

4. Yanyan Qi

   1. Center for Personal Dynamic Regulomes, Stanford, United States
   2. Department of Dermatology, Stanford University School of Medicine, Stanford, United States
   3. Department of Genetics, Stanford University School of Medicine, Stanford, United States

   Contribution

   Resources, Data curation, Investigation

   Competing interests

   No competing interests declared

5. M Ryan Corces

   1. Center for Personal Dynamic Regulomes, Stanford, United States
   2. Department of Dermatology, Stanford University School of Medicine, Stanford, United States
   3. Department of Genetics, Stanford University School of Medicine, Stanford, United States

   Contribution

   Resources, Data curation

   Competing interests

   No competing interests declared

6. Ravindra Majeti

   1. Institute for Stem Cell Biology and Regenerative Medicine, Stanford University School of Medicine, Stanford, United States
   2. Division of Hematology, Department of Medicine, Stanford University School of Medicine, Stanford, United States

   Contribution

   Resources, Data curation, Supervision, Funding acquisition, Methodology, Writing—review and editing

   For correspondence

   rmajeti@stanford.edu

   Competing interests

   is a co-founder, consultant, equity holder, and serves on the Board of Directors of Forty Seven Inc. Ravindra Majeti is Reviewing editor, eLife

7. Howard Y Chang

   1. Center for Personal Dynamic Regulomes, Stanford, United States
   2. Department of Dermatology, Stanford University School of Medicine, Stanford, United States
   3. Department of Genetics, Stanford University School of Medicine, Stanford, United States
   4. Howard Hughes Medical Institute, Stanford University, Stanford, United States

   Contribution

   Conceptualization, Supervision, Funding acquisition, Writing—original draft, Writing—review and editing

   For correspondence

   howchang@stanford.edu

   Competing interests

   is a co-founder of Accent Therapeutics and an advisor for 10x Genomics and Spring Discovery. Stanford University has filed a patent on ATAC-seq(US20160060691A1), on which HYC is named as an inventor.

   "This ORCID iD identifies the author of this article:" 0000-0002-9459-4393

• 12,503

  views

• 1,548

  downloads

• 104

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.