Abstract
Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multi-system disease. The wide range of manifestations observed in mitochondrial disease patients results from varying fractions of abnormal mtDNA molecules in different cells and tissues, a phenomenon termed heteroplasmy. However, the landscape of heteroplasmy across cell types within tissues and its influence on phenotype expression in affected patients remains largely unexplored. Here, we identify non- random distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA sequencing, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profile the transcriptome, chromatin accessibility state, and heteroplasmy in cells from the eyes of a patient with mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes (MELAS) and healthy control donors. Utilizing the retina as a model for complex multi-lineage tissues, we found that the proportion of the pathogenic m.3243A>G allele was neither evenly nor randomly distributed across diverse cell types. All neuroectoderm- derived neural cells exhibited a high percentage of the mutant variant. However, a subset of mesoderm- derived lineage, namely the vasculature of the choroid, was near homoplasmic for the wildtype allele. Gene expression and chromatin accessibility profiles of cell types with high and low proportions of m.3243A>G implicate mTOR signaling in the cellular response to heteroplasmy. We further found by multimodal single-cell sequencing of retinal pigment epithelial cells that a high proportion of the pathogenic mtDNA variant was associated with transcriptionally and morphologically abnormal cells. Together, these findings show the non- random nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.
Publisher
Cold Spring Harbor Laboratory