Go to The Journal of Clinical Investigation
  • About
  • Editors
  • Consulting Editors
  • For authors
  • Journal stats
  • Publication ethics
  • Publication alerts by email
  • Transfers
  • Advertising
  • Job board
  • Contact
  • Physician-Scientist Development
  • Current issue
  • Past issues
  • By specialty
    • COVID-19
    • Cardiology
    • Immunology
    • Metabolism
    • Nephrology
    • Oncology
    • Pulmonology
    • All ...
  • Videos
  • Collections
    • In-Press Preview
    • Resource and Technical Advances
    • Clinical Research and Public Health
    • Research Letters
    • Editorials
    • Perspectives
    • Physician-Scientist Development
    • Reviews
    • Top read articles

  • Current issue
  • Past issues
  • Specialties
  • In-Press Preview
  • Resource and Technical Advances
  • Clinical Research and Public Health
  • Research Letters
  • Editorials
  • Perspectives
  • Physician-Scientist Development
  • Reviews
  • Top read articles
  • About
  • Editors
  • Consulting Editors
  • For authors
  • Journal stats
  • Publication ethics
  • Publication alerts by email
  • Transfers
  • Advertising
  • Job board
  • Contact
Multimodal single-cell analysis of nonrandom heteroplasmy distribution in human retinal mitochondrial disease
Nathaniel K. Mullin, Andrew P. Voigt, Miles J. Flamme-Wiese, Xiuying Liu, Megan J. Riker, Katayoun Varzavand, Edwin M. Stone, Budd A. Tucker, Robert F. Mullins
Nathaniel K. Mullin, Andrew P. Voigt, Miles J. Flamme-Wiese, Xiuying Liu, Megan J. Riker, Katayoun Varzavand, Edwin M. Stone, Budd A. Tucker, Robert F. Mullins
View: Text | PDF
Research Article Genetics Ophthalmology

Multimodal single-cell analysis of nonrandom heteroplasmy distribution in human retinal mitochondrial disease

  • Text
  • PDF
Abstract

Variants within the high copy number mitochondrial genome (mtDNA) can disrupt organelle function and lead to severe multisystem disease. The wide range of manifestations observed in patients with mitochondrial disease 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 nonrandom distribution of a pathogenic mtDNA variant across a complex tissue using single-cell RNA-Seq, mitochondrial single-cell ATAC sequencing, and multimodal single-cell sequencing. We profiled 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 from healthy control donors. Utilizing the retina as a model for complex multilineage 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 WT 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 nonrandom nature of mitochondrial variant partitioning in human mitochondrial disease and underscore its implications for mitochondrial disease pathogenesis and treatment.

Authors

Nathaniel K. Mullin, Andrew P. Voigt, Miles J. Flamme-Wiese, Xiuying Liu, Megan J. Riker, Katayoun Varzavand, Edwin M. Stone, Budd A. Tucker, Robert F. Mullins

×

Figure 5

Differential accessibility of transcription factor binding motifs in MELAS choroid.

Options: View larger image (or click on image) Download as PowerPoint
Differential accessibility of transcription factor binding motifs in MEL...
(A) Top transcription factor binding motifs enriched in MELAS or control T cells. The YY1 and ZFP42 motifs are deenriched in MELAS T cells. (B) Motif logo of the JASPAR YY1 motif (MA0095.2). (C) Motif foot printing plot shows Tn5 transposase insertion enrichment in the locus flanking the YY1 binding motif. Traces for 4 cell types with lower YY1 motif enrichment in MELAS are shown. Foot printing traces represent accessibility in MELAS and control cells. (D) YY1 binding motif enrichment score shown for all choroidal cell types. Four cell types with low m.3243A>G heteroplasmy have lower YY1 binding motif enrichment in MELAS versus control. (E) Enriched motifs in choroidal endothelial cells. (F) RFX3 motif logo, one of several RFX factor motifs enriched in control endothelial cells versus controls. (G) RFX3 motif footprint plot showing traces for T cells and endothelial cells. (H) ChromVAR activity for the RFX3 motif is lower in MELAS endothelial and T cells in MELAS sample versus control. (I) Enriched motifs in the pericyte and SMC population of MELAS versus control. (J) The NR3C2 binding motif. (K) NR3C2 footprint plot indicating transcription factor-chromatin interaction in cell types of interest. (L) ChromVAR activity of NR3C2 is elevated in the low-heteroplasmy pericyte/SMC, endothelial, and Schwann cell populations in MELAS versus control.

Copyright © 2026 American Society for Clinical Investigation
ISSN 2379-3708

Sign up for email alerts