Pervasive within-Mitochondrion Single-Nucleotide Variant Heteroplasmy as Revealed by Single-Mitochondrion Sequencing

Mitochondrial heterogeneity and crosstalk in aging: Time for a paradigm shift?

The hallmarks of aging have been influential in guiding the biology of aging research, with more recent and growing recognition of the interdependence of these hallmarks on age-related health outcomes. However, a current challenge is personalizing aging trajectories to promote healthy aging, given the diversity of genotypes and lived experience. We suggest that incorporating heterogeneity-including intrinsic (e.g., genetic and structural) and extrinsic (e.g., environmental and exposome) factors and their interdependence of hallmarks-may move the dial. This editorial perspective will focus on one hallmark, namely mitochondrial dysfunction, to exemplify how consideration of heterogeneity and interdependence or crosstalk may reveal new perspectives and opportunities for personalizing aging research. To this end, we highlight heterogeneity within mitochondria as a model.

Mitochondrial genetic variations in leukemia: a comprehensive overview

Leukemias are a group of heterogeneous hematological malignancies driven by diverse genetic variations, and the advent of genomic sequencing technologies facilitates the investigation of genetic abnormalities in leukemia. However, these sequencing-based studies mainly focus on nuclear DNAs. Increasing evidence indicates that mitochondrial dysfunction is an important mechanism of leukemia pathogenesis, which is closely related to the mitochondrial genome variations. Here, we provide an overview of current research progress concerning mitochondrial genetic variations in leukemia, encompassing gene mutations and copy number variations. We also summarize currently accessible mitochondrial DNA (mtDNA) sequencing methods. Notably, somatic mtDNA mutations may serve as natural genetic barcodes for lineage tracing and longitudinal assessment of clonal dynamics. Collectively, these findings enhance our understanding of leukemia pathogenesis and foster the identification of novel therapeutic targets and interventions.

Unveiling aging dynamics in the hematopoietic system insights from single-cell technologies

As the demographic structure shifts towards an aging society, strategies aimed at slowing down or reversing the aging process become increasingly essential. Aging is a major predisposing factor for many chronic diseases in humans. The hematopoietic system, comprising blood cells and their associated bone marrow microenvironment, intricately participates in hematopoiesis, coagulation, immune regulation and other physiological phenomena. The aging process triggers various alterations within the hematopoietic system, serving as a spectrum of risk factors for hematopoietic disorders, including clonal hematopoiesis, immune senescence, myeloproliferative neoplasms and leukemia. The emerging single-cell technologies provide novel insights into age-related changes in the hematopoietic system. In this review, we summarize recent studies dissecting hematopoietic system aging using single-cell technologies. We discuss cellular changes occurring during aging in the hematopoietic system at the levels of the genomics, transcriptomics, epigenomics, proteomics, metabolomics and spatial multi-omics. Finally, we contemplate the future prospects of single-cell technologies, emphasizing the impact they may bring to the field of hematopoietic system aging research.

Single-mitochondrion sequencing uncovers distinct mutational patterns and heteroplasmy landscape in mouse astrocytes and neurons

BACKGROUND

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting "normal" mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes, we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes. In this study, we amplified mt-genomes from 1645 single mitochondria isolated from mouse single astrocytes and neurons to (1) determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, (2) assess differences in mtDNA SNVs between neurons and astrocytes, and (3) study co-segregation of variants in the mouse mtDNA.

RESULTS

(1) The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. (2) Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G > A and 9419:C > T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. (3) Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

CONCLUSIONS

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Single-Mitochondrion Sequencing Uncovers Distinct Mutational Patterns and Heteroplasmy Landscape in Mouse Astrocytes and Neurons

Background

Mitochondrial (mt) heteroplasmy can cause adverse biological consequences when deleterious mtDNA mutations accumulate disrupting 'normal' mt-driven processes and cellular functions. To investigate the heteroplasmy of such mtDNA changes we developed a moderate throughput mt isolation procedure to quantify the mt single-nucleotide variant (SNV) landscape in individual mouse neurons and astrocytes In this study we amplified mt-genomes from 1,645 single mitochondria (mts) isolated from mouse single astrocytes and neurons to 1. determine the distribution and proportion of mt-SNVs as well as mutation pattern in specific target regions across the mt-genome, 2. assess differences in mtDNA SNVs between neurons and astrocytes, and 3. Study cosegregation of variants in the mouse mtDNA.

Results

1. The data show that specific sites of the mt-genome are permissive to SNV presentation while others appear to be under stringent purifying selection. Nested hierarchical analysis at the levels of mitochondrion, cell, and mouse reveals distinct patterns of inter- and intra-cellular variation for mt-SNVs at different sites. 2. Further, differences in the SNV incidence were observed between mouse neurons and astrocytes for two mt-SNV 9027:G>A and 9419:C>T showing variation in the mutational propensity between these cell types. Purifying selection was observed in neurons as shown by the Ka/Ks statistic, suggesting that neurons are under stronger evolutionary constraint as compared to astrocytes. 3. Intriguingly, these data show strong linkage between the SNV sites at nucleotide positions 9027 and 9461.

Conclusion

This study suggests that segregation as well as clonal expansion of mt-SNVs is specific to individual genomic loci, which is important foundational data in understanding of heteroplasmy and disease thresholds for mutation of pathogenic variants.

Nanobiopsy investigation of the subcellular mtDNA heteroplasmy in human tissues

Mitochondrial function is critical to continued cellular vitality and is an important contributor to a growing number of human diseases. Mitochondrial dysfunction is typically heterogeneous, mediated through the clonal expansion of mitochondrial DNA (mtDNA) variants in a subset of cells in a given tissue. To date, our understanding of the dynamics of clonal expansion of mtDNA variants has been technically limited to the single cell-level. Here, we report the use of nanobiopsy for subcellular sampling from human tissues, combined with next-generation sequencing to assess subcellular mtDNA mutation load in human tissue from mitochondrial disease patients. The ability to map mitochondrial mutation loads within individual cells of diseased tissue samples will further our understanding of mitochondrial genetic diseases.

High-Throughput Single-Cell, Single-Mitochondrial DNA Assay Using Hydrogel Droplet Microfluidics

There is growing interest in understanding the biological implications of single cell heterogeneity and heteroplasmy of mitochondrial DNA (mtDNA), but current methodologies for single-cell mtDNA analysis limit the scale of analysis to small cell populations. Although droplet microfluidics have increased the throughput of single-cell genomic, RNA, and protein analysis, their application to sub-cellular organelle analysis has remained a largely unsolved challenge. Here, we introduce an agarose-based droplet microfluidic approach for single-cell, single-mtDNA analysis, which allows simultaneous processing of hundreds of individual mtDNA molecules within >10,000 individual cells. Our microfluidic chip encapsulates individual cells in agarose beads, designed to have a sufficiently dense hydrogel network to retain mtDNA after lysis and provide a robust scaffold for subsequent multi-step processing and analysis. To mitigate the impact of the high viscosity of agarose required for mtDNA retention on the throughput of microfluidics, we developed a parallelized device, successfully achieving ~95 % mtDNA retention from single cells within our microbeads at >700,000 drops/minute. To demonstrate utility, we analyzed specific regions of the single-mtDNA using a multiplexed rolling circle amplification (RCA) assay. We demonstrated compatibility with both microscopy, for digital counting of individual RCA products, and flow cytometry for higher throughput analysis.

Bioinformatics analyses and experimental validation of the role of phagocytosis in low-grade glioma

BACKGROUND

Phagocytosis is of vital importance in tumor immune response. The alteration of phagocytosis in low-grade glioma (LGG) has not been investigated.

METHODS

The mRNA, copy number variation, single nucleotide variation, and methylation levels of phagocytosis-related genes were summarized in pan-cancer. Non-negative matrix factorization clustering was utilized to identify two LGG subtypes. LASSO regression analysis was performed to construct a phagocytosis-related prognostic signature (PRPS). Immune characteristics, immunotherapy response, and targeted-drug sensitivity were further explored. The phagocytosis activity in glioma was evaluated using scRNA-seq data. Multiplex immunohistochemical (m-IHC) technology was performed to identify the tumor-infiltrating immune cells in LGG.

RESULTS

The phagocytosis-related genes altered obviously in pan-cancer compared with corresponding normal tissues. Two LGG subtypes were obtained and the subtype with poor prognosis was combined with lower tumor purity, more active immune-related pathways, increasing infiltration of CD4+ T cells, CD8+ T cells, and natural killer (NK) cells, decreasing infiltration of macrophages, mast cells, and neutrophils, distinct pathway activity and cell death status, greater response to immunotherapy, and higher sensitivity to cyclophosphamide, erlotinib, gefitinib, lapatinib, and sorafenib. In addition, a PRPS involving 10 genes (i.e., SLC11A1, CAMK1D, PLA2G5, STAP1, ALOX15, PLCG2, SFTPD, AZU1, RAB27A, and LAMTOR2) was constructed to estimate the risk level of each LGG sample and high risk LGG patients had poor prognosis, upregulated infiltration of neutrophil, macrophage, Treg, and myeloid dendritic cell, down regulated infiltration of monocyte and NK cell, and increasing expression of large number of immune checkpoint genes. The phagocytosis activity is notably active in monocyte/macrophage. The m-IHC results confirmed increased infiltration of macrophages and neutrophils in LGG samples with high SLC11A1 expression.

CONCLUSION

The molecular characteristics of phagocytosis were revealed and the PRPS laid the foundation for personalized therapy in LGG.

High-throughput single-cell, single-mitochondrial DNA assay using hydrogel droplet microfluidics

There is growing interest in understanding the biological implications of single cell heterogeneity and intracellular heteroplasmy of mtDNA, but current methodologies for single-cell mtDNA analysis limit the scale of analysis to small cell populations. Although droplet microfluidics have increased the throughput of single-cell genomic, RNA, and protein analysis, their application to sub-cellular organelle analysis has remained a largely unsolved challenge. Here, we introduce an agarose-based droplet microfluidic approach for single-cell, single-mtDNA analysis, which allows simultaneous processing of hundreds of individual mtDNA molecules within >10,000 individual cells. Our microfluidic chip encapsulates individual cells in agarose beads, designed to have a sufficiently dense hydrogel network to retain mtDNA after lysis and provide a robust scaffold for subsequent multi-step processing and analysis. To mitigate the impact of the high viscosity of agarose required for mtDNA retention on the throughput of microfluidics, we developed a parallelized device, successfully achieving ~95% mtDNA retention from single cells within our microbeads at >700,000 drops/minute. To demonstrate utility, we analyzed specific regions of the single mtDNA using a multiplexed rolling circle amplification (RCA) assay. We demonstrated compatibility with both microscopy, for digital counting of individual RCA products, and flow cytometry for higher throughput analysis.

Cellular barcoding: From developmental tracing to anti-tumor drug discovery

Clonal evolution has gained immense attention in explaining cancer cell status, history, and fate during cancer progression. Current single-cell or spatial transcriptome technologies have broadened our understanding of various mechanisms underlying cancer initiation, relapse, and drug resistance. However, technical challenges still hinder a better understanding of the dynamics of distinctive phenotypic states and abnormal trajectories from normal physiological transition to malignant stages. Cellular barcoding enabled lineage tracing on parallelly massive cells at single-cell resolution through different mechanisms lately, enabling new insights into exploring developmental trajectories, cancer progression, and targeted therapies. This review summarizes the latest noteworthy and robust strategies for different types of cellular barcodes. To introduce the major characteristics, advantages and limitations of these different strategies, this review will further guide in choosing or improving cellular barcoding technologies and their applications in cancer research.

Single-cell multi-omics of mitochondrial DNA disorders reveals dynamics of purifying selection across human immune cells

Pathogenic mutations in mitochondrial DNA (mtDNA) compromise cellular metabolism, contributing to cellular heterogeneity and disease. Diverse mutations are associated with diverse clinical phenotypes, suggesting distinct organ- and cell-type-specific metabolic vulnerabilities. Here we establish a multi-omics approach to quantify deletions in mtDNA alongside cell state features in single cells derived from six patients across the phenotypic spectrum of single large-scale mtDNA deletions (SLSMDs). By profiling 206,663 cells, we reveal the dynamics of pathogenic mtDNA deletion heteroplasmy consistent with purifying selection and distinct metabolic vulnerabilities across T-cell states in vivo and validate these observations in vitro. By extending analyses to hematopoietic and erythroid progenitors, we reveal mtDNA dynamics and cell-type-specific gene regulatory adaptations, demonstrating the context-dependence of perturbing mitochondrial genomic integrity. Collectively, we report pathogenic mtDNA heteroplasmy dynamics of individual blood and immune cells across lineages, demonstrating the power of single-cell multi-omics for revealing fundamental properties of mitochondrial genetics.

Quantitative haplotype-resolved analysis of mitochondrial DNA heteroplasmy in Human single oocytes, blastoids, and pluripotent stem cells

Maternal mitochondria are the sole source of mtDNA for every cell of the offspring. Heteroplasmic mtDNA mutations inherited from the oocyte are a common cause of metabolic diseases and associated with late-onset diseases. However, the origin and dynamics of mtDNA heteroplasmy remain unclear. We used our individual Mitochondrial Genome sequencing (iMiGseq) technology to study mtDNA heterogeneity, quantitate single nucleotide variants (SNVs) and large structural variants (SVs), track heteroplasmy dynamics, and analyze genetic linkage between variants at the individual mtDNA molecule level in single oocytes and human blastoids. Our study presented the first single-mtDNA analysis of the comprehensive heteroplasmy landscape in single human oocytes. Unappreciated levels of rare heteroplasmic variants well below the detection limit of conventional methods were identified in healthy human oocytes, of which many are reported to be deleterious and associated with mitochondrial disease and cancer. Quantitative genetic linkage analysis revealed dramatic shifts of variant frequency and clonal expansions of large SVs during oogenesis in single-donor oocytes. iMiGseq of a single human blastoid suggested stable heteroplasmy levels during early lineage differentiation of naïve pluripotent stem cells. Therefore, our data provided new insights of mtDNA genetics and laid a foundation for understanding mtDNA heteroplasmy at early stages of life.

High-frequency and functional mitochondrial DNA mutations at the single-cell level

Decline in mitochondrial function underlies aging and age-related diseases, but the role of mitochondrial DNA (mtDNA) mutations in these processes remains elusive. To investigate patterns of mtDNA mutations, it is particularly important to quantify mtDNA mutations and their associated pathogenic effects at the single-cell level. However, existing single-cell mtDNA sequencing approaches remain inefficient due to high cost and low mtDNA on-target rates. In this study, we developed a cost-effective mtDNA targeted-sequencing protocol called single-cell sequencing by targeted amplification of multiplex probes (scSTAMP) and experimentally validated its reliability. We then applied our method to assess single-cell mtDNA mutations in 768 B lymphocytes and 768 monocytes from a 76-y-old female. Across 632 B lymphocyte and 617 monocytes with medium mtDNA coverage over >100×, our results indicated that over 50% of cells carried at least one mtDNA mutation with variant allele frequencies (VAFs) over 20%, and that cells carried an average of 0.658 and 0.712 such mutation for B lymphocytes and monocytes, respectively. Surprisingly, more than 20% of the observed mutations had VAFs of over 90% in either cell population. In addition, over 60% of the mutations were in protein-coding genes, of which over 70% were nonsynonymous, and more than 50% of the nonsynonymous mutations were predicted to be highly pathogenic. Interestingly, about 80% of the observed mutations were singletons in the respective cell populations. Our results revealed mtDNA mutations with functional significance might be prevalent at advanced age, calling further investigation on age-related mtDNA mutation dynamics at the single-cell level.

Organelle bottlenecks facilitate evolvability by traversing heteroplasmic fitness valleys

Bioenergetic organelles-mitochondria and plastids-retain their own genomes (mtDNA and ptDNA), and these organelle DNA (oDNA) molecules are vital for eukaryotic life. Like all genomes, oDNA must be able to evolve to suit new environmental challenges. However, mixed oDNA populations in cells can challenge cellular bioenergetics, providing a penalty to the appearance and adaptation of new mutations. Here we show that organelle "bottlenecks," mechanisms increasing cell-to-cell oDNA variability during development, can overcome this mixture penalty and facilitate the adaptation of beneficial mutations. We show that oDNA heteroplasmy and bottlenecks naturally emerge in evolutionary simulations subjected to fluctuating environments, demonstrating that this evolvability is itself evolvable. Usually thought of as a mechanism to clear damaging mutations, organelle bottlenecks therefore also resolve the tension between intracellular selection for pure cellular oDNA populations and the "bet-hedging" need for evolvability and adaptation to new environments. This general theory suggests a reason for the maintenance of organelle heteroplasmy in cells, and may explain some of the observed diversity in organelle maintenance and inheritance across taxa.

Future of Digital Assays to Resolve Clinical Heterogeneity of Single Extracellular Vesicles

Extracellular vesicles (EVs) are complex lipid membrane vehicles with variable expressions of molecular cargo, composed of diverse subpopulations that participate in the intercellular signaling of biological responses in disease. EV-based liquid biopsies demonstrate invaluable clinical potential for overhauling current practices of disease management. Yet, EV heterogeneity is a major needle-in-a-haystack challenge to translate their use into clinical practice. In this review, existing digital assays will be discussed to analyze EVs at a single vesicle resolution, and future opportunities to optimize the throughput, multiplexing, and sensitivity of current digital EV assays will be highlighted. Furthermore, this review will outline the challenges and opportunities that impact the clinical translation of single EV technologies for disease diagnostics and treatment monitoring.

Analysis of Mitochondria by Single-Organelle Resolution

Cellular organelles are highly specialized compartments with distinct functions. With the increasing resolution of detection methods, it is becoming clearer that same organelles may have different functions or properties not only within different cell populations of a tissue but also within the same cell. Dysfunction or altered function affects the organelle itself and may also lead to malignancies or undesirable cell death. To understand cellular function or dysfunction, it is therefore necessary to analyze cellular components at the single-organelle level. Here, we review the recent advances in analyzing cellular function at single-organelle resolution using high-parameter flow cytometry or multicolor confocal microscopy. We focus on the analysis of mitochondria, as they are organelles at the crossroads of various cellular signaling pathways and functions. However, most of the applied methods/technologies are transferable to any other organelle, such as the endoplasmic reticulum, lysosomes, or peroxisomes.

A subcellular cookie cutter for spatial genomics in human tissue

Intracellular heterogeneity contributes significantly to cellular physiology and, in a number of debilitating diseases, cellular pathophysiology. This is greatly influenced by distinct organelle populations and to understand the aetiology of disease, it is important to have tools able to isolate and differentially analyse organelles from precise location within tissues. Here, we report the development of a subcellular biopsy technology that facilitates the isolation of organelles, such as mitochondria, from human tissue. We compared the subcellular biopsy technology to laser capture microdissection (LCM) that is the state-of-the-art technique for the isolation of cells from their surrounding tissues. We demonstrate an operational limit of  >20 µm for LCM and then, for the first time in human tissue, show that subcellular biopsy can be used to isolate mitochondria beyond this limit.

Mitochondrial DNA Mutagenesis: Feature of and Biomarker for Environmental Exposures and Aging

PURPOSE OF REVIEW

Mitochondrial dysfunction is a hallmark of aging. Mitochondrial genome (mtDNA) instability contributes to mitochondrial dysfunction, and mtDNA mutagenesis may contribute to aging. However, the origin of mtDNA mutations remains somewhat controversial. The goals of this review are to introduce and review recent literature on mtDNA mutagenesis and aging, address recent animal and epidemiological evidence for the effects of chemicals on mtDNA damage and mutagenesis, propose hypotheses regarding the contribution of environmental toxicant exposure to mtDNA mutagenesis in the context of aging, and suggest future directions and approaches for environmental health researchers.

RECENT FINDINGS

Stressors such as pollutants, pharmaceuticals, and ultraviolet radiation can damage the mitochondrial genome or disrupt mtDNA replication, repair, and organelle homeostatic processes, potentially influencing the rate of accumulation of mtDNA mutations. Accelerated mtDNA mutagenesis could contribute to aging, diseases of aging, and sensitize individuals with pathogenic mtDNA variants to stressors. We propose three potential mechanisms of toxicant-induced effects on mtDNA mutagenesis over lifespan: (1) increased de novo mtDNA mutations, (2) altered frequencies of mtDNA mutations, or (3) both. There are remarkably few studies that have investigated the impact of environmental chemical exposures on mtDNA instability and mutagenesis, and even fewer in the context of aging. More studies are warranted because people are exposed to tens of thousands of chemicals, and are living longer. Finally, we suggest that toxicant-induced mtDNA damage and mutational signatures may be a sensitive biomarker for some exposures.

A coarse-grained NADH redox model enables inference of subcellular metabolic fluxes from fluorescence lifetime imaging

Mitochondrial metabolism is of central importance to diverse aspects of cell and developmental biology. Defects in mitochondria are associated with many diseases, including cancer, neuropathology, and infertility. Our understanding of mitochondrial metabolism in situ and dysfunction in diseases are limited by the lack of techniques to measure mitochondrial metabolic fluxes with sufficient spatiotemporal resolution. Herein, we developed a new method to infer mitochondrial metabolic fluxes in living cells with subcellular resolution from fluorescence lifetime imaging of NADH. This result is based on the use of a generic coarse-grained NADH redox model. We tested the model in mouse oocytes and human tissue culture cells subject to a wide variety of perturbations by comparing predicted fluxes through the electron transport chain (ETC) to direct measurements of oxygen consumption rate. Interpreting the fluorescence lifetime imaging microscopy measurements of NADH using this model, we discovered a homeostasis of ETC flux in mouse oocytes: perturbations of nutrient supply and energy demand of the cell do not change ETC flux despite significantly impacting NADH metabolic state. Furthermore, we observed a subcellular spatial gradient of ETC flux in mouse oocytes and found that this gradient is primarily a result of a spatially heterogeneous mitochondrial proton leak. We concluded from these observations that ETC flux in mouse oocytes is not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration.

New insights into mitophagy and stem cells

Mitophagy is a specific autophagic phenomenon in which damaged or redundant mitochondria are selectively cleared by autophagic lysosomes. A decrease in mitophagy can accelerate the aging process. Mitophagy is related to health and longevity and is the key to protecting stem cells from metabolic stress damage. Mitophagy decreases the metabolic level of stem cells by clearing active mitochondria, so mitophagy is becoming increasingly necessary to maintain the regenerative capacity of old stem cells. Stem cell senescence is the core problem of tissue aging, and tissue aging occurs not only in stem cells but also in transport amplifying cell chambers and the stem cell environment. The loss of the autophagic ability of stem cells can cause the accumulation of mitochondria and the activation of the metabolic state as well as damage the self-renewal ability and regeneration potential of stem cells. However, the claim remains controversial. Mitophagy is an important survival strategy against nutrient deficiency and starvation, and mitochondrial function and integrity may affect the viability, proliferation and differentiation potential, and longevity of normal stem cells. Mitophagy can affect the health and longevity of the human body, so the number of studies in this field has increased, but the mechanism by which mitophagy participates in stem cell development is still not fully understood. This review describes the potential significance of mitophagy in stem cell developmental processes, such as self-renewal, differentiation and aging. Through this work, we discovered the role and mechanism of mitophagy in different types of stem cells, identified novel targets for killing cancer stem cells and curing cancer, and provided new insights for future research in this field.

Physical bioenergetics: Energy fluxes, budgets, and constraints in cells

Cells are the basic units of all living matter which harness the flow of energy to drive the processes of life. While the biochemical networks involved in energy transduction are well-characterized, the energetic costs and constraints for specific cellular processes remain largely unknown. In particular, what are the energy budgets of cells? What are the constraints and limits energy flows impose on cellular processes? Do cells operate near these limits, and if so how do energetic constraints impact cellular functions? Physics has provided many tools to study nonequilibrium systems and to define physical limits, but applying these tools to cell biology remains a challenge. Physical bioenergetics, which resides at the interface of nonequilibrium physics, energy metabolism, and cell biology, seeks to understand how much energy cells are using, how they partition this energy between different cellular processes, and the associated energetic constraints. Here we review recent advances and discuss open questions and challenges in physical bioenergetics.

Massively parallel single-cell mitochondrial DNA genotyping and chromatin profiling

Natural mitochondrial DNA (mtDNA) mutations enable the inference of clonal relationships among cells. mtDNA can be profiled along with measures of cell state, but has not yet been combined with the massively parallel approaches needed to tackle the complexity of human tissue. Here, we introduce a high-throughput, droplet-based mitochondrial single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), a method that combines high-confidence mtDNA mutation calling in thousands of single cells with their concomitant high-quality accessible chromatin profile. This enables the inference of mtDNA heteroplasmy, clonal relationships, cell state and accessible chromatin variation in individual cells. We reveal single-cell variation in heteroplasmy of a pathologic mtDNA variant, which we associate with intra-individual chromatin variability and clonal evolution. We clonally trace thousands of cells from cancers, linking epigenomic variability to subclonal evolution, and infer cellular dynamics of differentiating hematopoietic cells in vitro and in vivo. Taken together, our approach enables the study of cellular population dynamics and clonal properties in vivo.

Mitochondrial isolation: when size matters

Mitochondrial vitality is critical to cellular function, with mitochondrial dysfunction linked to a growing number of human diseases. Tissue and cellular heterogeneity, in terms of genetics, dynamics and function means that increasingly mitochondrial research is conducted at the single cell level. Whilst there are several technologies that are currently available for single-cell analysis, each with their advantages, they cannot be easily adapted to study mitochondria with subcellular resolution. Here we review the current techniques and strategies for mitochondrial isolation, critically discussing each technology's limitations for future mitochondrial research. Finally, we highlight and discuss the recent breakthroughs in sub-cellular isolation techniques, with a particular focus on nanotechnologies that enable the isolation of mitochondria from subcellular compartments. This allows isolation of mitochondria with unprecedented spatial precision with minimal disruption to mitochondria and their immediate cellular environment.

Micro- and Nano-Devices for Studying Subcellular Biology

Cells are complex machines whose behaviors arise from their internal collection of dynamically interacting organelles, supramolecular complexes, and cytoplasmic chemicals. The current understanding of the nature by which subcellular biology produces cell-level behaviors is limited by the technological hurdle of measuring the large number (>103 ) of small-sized (<1 μm) heterogeneous organelles and subcellular structures found within each cell. In this review, the emergence of a suite of micro- and nano-technologies for studying intracellular biology on the scale of organelles is described. Devices that use microfluidic and microelectronic components for 1) extracting and isolating subcellular structures from cells and lysate; 2) analyzing the physiology of individual organelles; and 3) recreating subcellular assembly and functions in vitro, are described. The authors envision that the continued development of single organelle technologies and analyses will serve as a foundation for organelle systems biology and will allow new insight into fundamental and clinically relevant biological questions.

Mitochondrial isolation: when size matters

Mitochondrial vitality is critical to cellular function, with mitochondrial dysfunction linked to a growing number of human diseases. Tissue and cellular heterogeneity, in terms of genetics, dynamics and function means that increasingly mitochondrial research is conducted at the single cell level. Whilst there are several technologies that are currently available for single-cell analysis, each with their advantages, they cannot be easily adapted to study mitochondria with subcellular resolution. Here we review the current techniques and strategies for mitochondrial isolation, critically discussing each technology's limitations for future mitochondrial research. Finally, we highlight and discuss the recent breakthroughs in sub-cellular isolation techniques, with a particular focus on nanotechnologies that enable the isolation of mitochondria from subcellular compartments. This allows isolation of mitochondria with unprecedented spatial precision with minimal disruption to mitochondria and their immediate cellular environment.

STAMP: a multiplex sequencing method for simultaneous evaluation of mitochondrial DNA heteroplasmies and content

Human mitochondrial genome (mtDNA) variations, such as mtDNA heteroplasmies (the co-existence of mutated and wild-type mtDNA), have received increasing attention in recent years for their clinical relevance to numerous diseases. But large-scale population studies of mtDNA heteroplasmies have been lagging due to the lack of a labor- and cost-effective method. Here, we present a novel human mtDNA sequencing method called STAMP (sequencing by targeted amplification of multiplex probes) for measuring mtDNA heteroplasmies and content in a streamlined workflow. We show that STAMP has high-mapping rates to mtDNA, deep coverage of unique reads and high tolerance to sequencing and polymerase chain reaction errors when applied to human samples. STAMP also has high sensitivity and low false positive rates in identifying artificial mtDNA variants at fractions as low as 0.5% in genomic DNA samples. We further extend STAMP, by including nuclear DNA-targeting probes, to enable assessment of relative mtDNA content in the same assay. The high cost-effectiveness of STAMP, along with the flexibility of using it for measuring various aspects of mtDNA variations, will accelerate the research of mtDNA heteroplasmies and content in large population cohorts, and in the context of human diseases and aging.

The mitochondria regulation of stem cell aging

Mitochondrial dysfunction and stem cell exhaustion are among the nine separate hallmarks of aging. Emerging evidence however suggests that mitochondrial activity can have a profound influence on the self-renewal and function of stem cells, thus mechanistically linking mitochondrial function and stem cell decline. In this review, we discuss how accumulation of mtDNA mutations or alterations in mitochondrial dynamics, turnover, and signaling can modulate age-dependent stem cell function. Finally, we also describe how mitochondrial substrate utilization influences stem and progenitor activity. Together, this growing body of evidence suggests that modulation of mitochondrial activity might provide a strategy to slow or reverse age-dependent stem cell decline, and potentially, slow or reverse human aging.

Marked mitochondrial genetic variation in individuals and populations of the carcinogenic liver fluke Clonorchis sinensis

Clonorchiasis is a neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, and is often associated with a malignant form of bile duct cancer (cholangiocarcinoma). Although some aspects of the epidemiology of clonorchiasis are understood, little is known about the genetics of C. sinensis populations. Here, we conducted a comprehensive genetic exploration of C. sinensis from endemic geographic regions using complete mitochondrial protein gene sets. Genomic DNA samples from C. sinensis individuals (n = 183) collected from cats and dogs in China (provinces of Guangdong, Guangxi, Hunan, Heilongjiang and Jilin) as well as from rats infected with metacercariae from cyprinid fish from the Russian Far East (Primorsky Krai region) were deep sequenced using the BGISEQ-500 platform. Informatic analyses of mitochondrial protein gene data sets revealed marked genetic variation within C. sinensis; significant variation was identified within and among individual worms from distinct geographical locations. No clear affiliation with a particular location or host species was evident, suggesting a high rate of dispersal of the parasite across endemic regions. The present work provides a foundation for future biological, epidemiological and ecological studies using mitochondrial protein gene data sets, which could aid in elucidating associations between particular C. sinensis genotypes/haplotypes and the pathogenesis or severity of clonorchiasis and its complications (including cholangiocarcinoma) in humans.

Clonorchiasis is an important neglected tropical disease caused by the Chinese liver fluke, Clonorchis sinensis, which can induce malignant bile duct cancer (cholangiocarcinoma). Little precise information is available on the biology, epidemiology and population genetics of C. sinensis. For this reason, we explored here the genetic composition of C. sinensis populations in distinct endemic areas in China and Russia. Using a deep sequencing-informatic approach, we revealed marked mitochondrial genetic variation within and between individuals and populations of C. sinensis, with no particular affiliation with geographic or host origin. These molecular findings and the methodology established should underpin future genetic studies of C. sinensis causing human clonorchiasis and associated complications (cancer) as well as transmission patterns in endemic regions.

Identification of unique and shared mitochondrial DNA mutations in neurodegeneration and cancer by single-cell mitochondrial DNA structural variation sequencing (MitoSV-seq)

Background

Point mutations and structural variations (SVs) in mitochondrial DNA (mtDNA) contribute to many neurodegenerative diseases. Technical limitations and heteroplasmy, however, have impeded their identification, preventing these changes from being examined in neurons in healthy and disease states.

Methods

We have developed a high-resolution technique—Mitochondrial DNA Structural Variation Sequencing (MitoSV-seq)—that identifies all types of mtDNA SVs and single-nucleotide variations (SNVs) in single neurons and novel variations that have been undetectable with conventional techniques.

Findings

Using MitoSV-seq, we discovered SVs/SNVs in dopaminergic neurons in the Ifnar1^−/− murine model of Parkinson disease. Further, MitoSV-seq was found to have broad applicability, delivering high-quality, full-length mtDNA sequences in a species-independent manner from human PBMCs, haematological cancers, and tumour cell lines, regardless of heteroplasmy. We characterised several common SVs in haematological cancers (AML and MDS) that were linked to the same mtDNA region, MT-ND5, using only 10 cells, indicating the power of MitoSV-seq in determining single-cancer-cell ontologies. Notably, the MT-ND5 hotspot, shared between all examined cancers and Ifnar1^−/− dopaminergic neurons, suggests that its mutations have clinical value as disease biomarkers.

Interpretation

MitoSV-seq identifies disease-relevant mtDNA mutations in single cells with high resolution, rendering it a potential drug screening platform in neurodegenerative diseases and cancers.

Funding

The Lundbeck Foundation, Danish Council for Independent Research-Medicine, and European Union Horizon 2020 Research and Innovation Programme.

High throughput single cell analysis of mitochondrial heteroplasmy in mitochondrial diseases

Mitochondrial heteroplasmy, which fundamentally means intracellular heterogeneity of mitochondrial DNA (mtDNA), has been measured in a group of cells, regardless of intercellular heterogeneity. Ordinal methods for mitochondrial heteroplasmy cannot discriminate between an intercellular homogenic population composed of cells with similar intracellular heterogeneity for mtDNA and an intercellular heterogenic population composed of cells with different rates of mutated mtDNA. A high-throughput method to determine mitochondrial heteroplasmy in a single cell was developed by using droplet digital PCR with TaqMan polymerase in this study. This technique revealed that there are three different cell populations of cultured fibroblasts derived from patients with mitochondrial disease carrying a mutation in the mtDNA; cells with homoplasmy of either mutated or healthy mtDNA; and cells mixed with mutated and healthy mtDNA. The presence of intercellular heterogeneity, even in uniformed cultured fibroblasts, suggests that heterogeneity should exist among different kinds of cells. The diagnosis of intercellular heterogeneity with respect to mitochondrial heteroplasmy by this methodology could provide novel insight into developing a treatment strategy for mitochondrial diseases.

Age-related accumulation of de novo mitochondrial mutations in mammalian oocytes and somatic tissues

Mutations create genetic variation for other evolutionary forces to operate on and cause numerous genetic diseases. Nevertheless, how de novo mutations arise remains poorly understood. Progress in the area is hindered by the fact that error rates of conventional sequencing technologies (1 in 100 or 1,000 base pairs) are several orders of magnitude higher than de novo mutation rates (1 in 10,000,000 or 100,000,000 base pairs per generation). Moreover, previous analyses of germline de novo mutations examined pedigrees (and not germ cells) and thus were likely affected by selection. Here, we applied highly accurate duplex sequencing to detect low-frequency, de novo mutations in mitochondrial DNA (mtDNA) directly from oocytes and from somatic tissues (brain and muscle) of 36 mice from two independent pedigrees. We found mtDNA mutation frequencies 2- to 3-fold higher in 10-month-old than in 1-month-old mice, demonstrating mutation accumulation during the period of only 9 mo. Mutation frequencies and patterns differed between germline and somatic tissues and among mtDNA regions, suggestive of distinct mutagenesis mechanisms. Additionally, we discovered a more pronounced genetic drift of mitochondrial genetic variants in the germline of older versus younger mice, arguing for mtDNA turnover during oocyte meiotic arrest. Our study deciphered for the first time the intricacies of germline de novo mutagenesis using duplex sequencing directly in oocytes, which provided unprecedented resolution and minimized selection effects present in pedigree studies. Moreover, our work provides important information about the origins and accumulation of mutations with aging/maturation and has implications for delayed reproduction in modern human societies. Furthermore, the duplex sequencing method we optimized for single cells opens avenues for investigating low-frequency mutations in other studies.

Disentangling the intertwined roles of mutation, selection and drift in the mitochondrial genome

Understanding and quantifying the rates of change in the mitochondrial genome is a major component of many areas of biological inquiry, from phylogenetics to human health. A critical parameter in understanding rates of change is estimating the mitochondrial mutation rate (mtDNA MR). Although the first direct estimates of mtDNA MRs were reported almost 20 years ago, the number of estimates has not grown markedly since that time. This is largely owing to the challenges associated with time- and labour-intensive mutation accumulation (MA) experiments. But even MA experiments do not solve a major problem with estimating mtDNA MRs-the challenge of disentangling the role of mutation from other evolutionary forces acting within the cell. Now that it is widely understood that any newly generated mutant allele in the mitochondria will initially be at very low frequency (1/N, where N is the number of mtDNA molecules in the cell), the importance of understanding the effective population size (Ne) of the mtDNA and the size of genetic bottlenecks during gametogenesis and development has come into the spotlight. In addition to these factors regulating the role of genetic drift, advances in our understanding of mitochondrial replication and turnover allow us to more easily envision how natural selection within the cell might favour or purge mutations in multi-copy organellar genomes. Here, we review the unique features of the mitochondrial genome that pose a challenge for accurate MR estimation and discuss ways to overcome those challenges. Estimates of mtDNA MRs remain one of the most widely used parameters in biology, thus accurate quantification and a deeper understanding of how and why they may vary within and between individuals, populations and species is an important goal. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Mitochondrial Network State Scales mtDNA Genetic Dynamics

Mitochondrial DNA (mtDNA) mutations cause severe congenital diseases but may also be associated with healthy aging. mtDNA is stochastically replicated and degraded, and exists within organelles which undergo dynamic fusion and fission. The role of the resulting mitochondrial networks in the time evolution of the cellular proportion of mutated mtDNA molecules (heteroplasmy), and cell-to-cell variability in heteroplasmy (heteroplasmy variance), remains incompletely understood. Heteroplasmy variance is particularly important since it modulates the number of pathological cells in a tissue. Here, we provide the first wide-reaching theoretical framework which bridges mitochondrial network and genetic states. We show that, under a range of conditions, the (genetic) rate of increase in heteroplasmy variance and de novo mutation are proportionally modulated by the (physical) fraction of unfused mitochondria, independently of the absolute fission-fusion rate. In the context of selective fusion, we show that intermediate fusion:fission ratios are optimal for the clearance of mtDNA mutants. Our findings imply that modulating network state, mitophagy rate, and copy number to slow down heteroplasmy dynamics when mean heteroplasmy is low could have therapeutic advantages for mitochondrial disease and healthy aging.

A nanoscale, multi-parametric flow cytometry-based platform to study mitochondrial heterogeneity and mitochondrial DNA dynamics

Mitochondria are well-characterized regarding their function in both energy production and regulation of cell death; however, the heterogeneity that exists within mitochondrial populations is poorly understood. Typically analyzed as pooled samples comprised of millions of individual mitochondria, there is little information regarding potentially different functionality across subpopulations of mitochondria. Herein we present a new methodology to analyze mitochondria as individual components of a complex and heterogeneous network, using a nanoscale and multi-parametric flow cytometry-based platform. We validate the platform using multiple downstream assays, including electron microscopy, ATP generation, quantitative mass-spectrometry proteomic profiling, and mtDNA analysis at the level of single organelles. These strategies allow robust analysis and isolation of mitochondrial subpopulations to more broadly elucidate the underlying complexities of mitochondria as these organelles function collectively within a cell.

Single-cell lineage tracing by endogenous mutations enriched in transposase accessible mitochondrial DNA

Simultaneous measurement of cell lineage and cell fates is a longstanding goal in biomedicine. Here we describe EMBLEM, a strategy to track cell lineage using endogenous mitochondrial DNA variants in ATAC-seq data. We show that somatic mutations in mitochondrial DNA can reconstruct cell lineage relationships at single cell resolution with high sensitivity and specificity. Using EMBLEM, we define the genetic and epigenomic clonal evolution of hematopoietic stem cells and their progenies in patients with acute myeloid leukemia. EMBLEM extends lineage tracing to any eukaryotic organism without genetic engineering.

Mitochondrial Heterogeneity

Cell-to-cell heterogeneity drives a range of (patho)physiologically important phenomena, such as cell fate and chemotherapeutic resistance. The role of metabolism, and particularly of mitochondria, is increasingly being recognized as an important explanatory factor in cell-to-cell heterogeneity. Most eukaryotic cells possess a population of mitochondria, in the sense that mitochondrial DNA (mtDNA) is held in multiple copies per cell, where the sequence of each molecule can vary. Hence, intra-cellular mitochondrial heterogeneity is possible, which can induce inter-cellular mitochondrial heterogeneity, and may drive aspects of cellular noise. In this review, we discuss sources of mitochondrial heterogeneity (variations between mitochondria in the same cell, and mitochondrial variations between supposedly identical cells) from both genetic and non-genetic perspectives, and mitochondrial genotype-phenotype links. We discuss the apparent homeostasis of mtDNA copy number, the observation of pervasive intra-cellular mtDNA mutation (which is termed "microheteroplasmy"), and developments in the understanding of inter-cellular mtDNA mutation ("macroheteroplasmy"). We point to the relationship between mitochondrial supercomplexes, cristal structure, pH, and cardiolipin as a potential amplifier of the mitochondrial genotype-phenotype link. We also discuss mitochondrial membrane potential and networks as sources of mitochondrial heterogeneity, and their influence upon the mitochondrial genome. Finally, we revisit the idea of mitochondrial complementation as a means of dampening mitochondrial genotype-phenotype links in light of recent experimental developments. The diverse sources of mitochondrial heterogeneity, as well as their increasingly recognized role in contributing to cellular heterogeneity, highlights the need for future single-cell mitochondrial measurements in the context of cellular noise studies.

The Cellular Mitochondrial Genome Landscape in Disease

Mitochondrial genome (mitochondrial DNA, mtDNA) lesions that unbalance bioenergetic and oxidative outputs are an important cause of human disease. A major impediment in our understanding of the pathophysiology of mitochondrial disorders is the complexity with which mtDNA mutations are spatiotemporally distributed and managed within individual cells, tissues, and organs. Unlike the comparatively static nuclear genome, accumulating evidence highlights the variability, dynamism, and modifiability of the mtDNA nucleotide sequence between individual cells over time. In this review, we summarize and discuss the impact of mtDNA defects on disease within the context of a mosaic and shifting mutational landscape.

Multilevel heterogeneity of mitochondrial respiratory chain deficiency

Mitochondrial diseases are heterogeneous multisystem disorders that show a mosaic pattern of mitochondrial respiratory chain dysfunction. The mitochondrial DNA (mtDNA) mutation load is heterogeneous at multiple levels: across organs, between cells, and between subcellular compartments. Such heterogeneity poses a diagnostic challenge, but also provides a scientific opportunity to explore the biological mechanisms underlying the onset and progression of these disorders. A recent article in The Journal of Pathology described a novel histochemical technique - nitro blue tetrazolium exclusion assay (NBTx) - to quantify mitochondrial cytochrome c oxidase (COX, or complex IV) deficiency. This technique is rapid, cost-effective, and quantitative, and is more sensitive than previous histochemical methods. It can also be applied across model organisms and human tissues. The NBTx method should therefore be a useful diagnostic tool, and may catalyze research examining the cellular and subcellular mechanisms that drive the onset and progression of inherited and acquired mtDNA disorders. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Germline and somatic mtDNA mutations in mouse aging

The accumulation of acquired mitochondrial genome (mtDNA) mutations with aging in somatic cells has been implicated in mitochondrial dysfunction and linked to age-onset diseases in humans. Here, we asked if somatic mtDNA mutations are also associated with aging in the mouse. MtDNA integrity in multiple organs and tissues in young and old (2-34 months) wild type (wt) mice was investigated by whole genome sequencing. Remarkably, no acquired somatic mutations were detected in tested tissues. However, we identified several non-synonymous germline mtDNA variants whose heteroplasmy levels (ratio of normal to mutant mtDNA) increased significantly with aging suggesting clonal expansion of inherited mtDNA mutations. Polg mutator mice, a model for premature aging, exhibited both germline and somatic mtDNA mutations whose numbers and heteroplasmy levels increased significantly with age implicating involvement in premature aging. Our results suggest that, in contrast to humans, acquired somatic mtDNA mutations do not accompany the aging process in wt mice.

Study of spontaneous mutations in the transmission of poplar chloroplast genomes from mother to offspring

Background

Chloroplasts have their own genomes, independent from nuclear genomes, that play vital roles in growth, which is a major targeted trait for genetic improvement in Populus. Angiosperm chloroplast genomes are maternally inherited, but the chloroplast’ variation pattern of poplar at the single-base level during the transmission from mother to offspring remains unknown.

Results

Here, we constructed high-quality and almost complete chloroplast genomes for three poplar clones, ‘NL895’ and its parents, ‘I69’ and ‘I45’, from the short-read datasets using multi-pass sequencing (15–16 times per clone) and ultra-high coverage (at least 8500× per clone), with the four-step strategy of Simulation–Assembly–Merging–Correction. Each of the three resulting chloroplast assemblies contained contigs covering > 99% of Populus trichocarpa chloroplast DNA as a reference. A total of 401 variant loci were identified by a hybrid strategy of genome comparison-based and mapping-based single nucleotide polymorphism calling. The genotypes of 94 variant loci were different among the three poplar clones. However, only 1 of the 94 loci was a missense mutation, which was located in the exon region of rpoC1 encoding the β’ subunit of plastid-encoded RNA polymerase. The genotype of the loci in NL895 and its female parent (I69) was different from that of its male parent (I45).

Conclusions

This research provides resources for further chloroplast genomic studies of a F₁ full-sibling family derived from a cross between I69 and I45, and will improve the application of chloroplast genomic information in modern Populus breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4813-8) contains supplementary material, which is available to authorized users.