Mitochondrial network determines intracellular ROS dynamics and sensitivity to oxidative stress through switching inter-mitochondrial messengers

Oxidative stresses caused by reactive oxygen species (ROS) can induce rapid depolarization of inner mitochondrial membrane potential and subsequent impairment of oxidative phosphorylation. Damaged mitochondria produce more ROS, especially the superoxide anion (O(2)(-)) and hydrogen peroxide (H(2)O(2)), which potentiate mitochondria-driven ROS propagation, so-called ROS-induced ROS release (RIRR), via activation of an inter-mitochondria signaling network. Therefore, loss of function in only a fraction of mitochondria might eventually affect cell viability through this positive feedback loop. Since ROS are very short-lived molecules in the biological milieu, mitochondrial network dynamics, such as density, number, and spatial distribution, can affect mitochondria-driven ROS propagation. To address this issue, we developed a mathematical model using an agent-based modeling approach, and tested the effect of mitochondrial network dynamics on RIRR for mitochondria under various conditions. Simulation results show that the intracellular ROS signaling pattern, such as ROS propagation speed and oxidative stress vulnerability, are critically affected by mitochondrial network dynamics. Mitochondrial network dynamics of mitochondrial distribution, density, activity, and size can mediate inter-mitochondrial signaling under certain conditions and determine the identity of the ROS signaling pattern. We further elucidated the potential mechanism of these actions, i.e., conversion of major messenger molecules involved in ROS signaling. If the average distance between neighboring mitochondria is large or mitochondrial distribution becomes randomized, messenger molecule of the ROS signaling network can be switched from O(2)(-) to H(2)O(2). In this case, mitochondria-driven ROS propagation is efficiently blocked by introduction of excess cytosolic glutathione peroxidase 1, while introduction of cytosolic superoxide dismutase has no effect. Together, these results suggest that mitochondrial network dynamics is a major determinant for cellular responses to RIRR through changing the key messenger molecules.     

The low conductance mitochondrial permeability transition pore confers excitability and CICR wave propagation in a computational model

Mitochondria have long been known to sequester cytosolic Ca²⁺ and even to shape intracellular patterns of endoplasmic reticulum-based Ca²⁺ signaling. Evidence suggests that the mitochondrial network is an excitable medium which can demonstrate independent Ca²⁺ induced Ca²⁺ release via the mitochondrial permeability transition. The role of this excitability remains unclear, but mitochondrial Ca²⁺ handling appears to be a crucial element in diverse diseases as diabetes, neurodegeneration and cardiac dysfunction that also have bioenergetic components. In this paper, we extend the modular Magnus-Keizer computational model for respiration-driven Ca²⁺ handling to include a permeability transition based on a channel-like pore mechanism. We demonstrate both excitability and Ca²⁺ wave propagation accompanied by depolarizations qualitatively similar to those reported in cell and isolated mitochondria preparations. These waves depend on the energy state of the mitochondria, as well as other elements of mitochondrial physiology. Our results support the concept that mitochondria can transmit state dependent signals about their function across the mitochondrial network. Our model provides the tools for predictions about the internal physiology that leads to this qualitatively different Ca²⁺ excitability seen in mitochondria.     

Mitochondrial Ca2+ dynamics reveals limited intramitochondrial Ca2+ diffusion

To reveal heterogeneity of mitochondrial function on the single-mitochondrion level we have studied the spatiotemporal dynamics of the mitochondrial Ca2+ signaling and the mitochondrial membrane potential using wide-field fluorescence imaging and digital image processing techniques. Here we demonstrate first-time discrete sites--intramitochondrial hotspots--of Ca2+ uptake after Ca2+ release from intracellular stores, and spreading of Ca2+ rise within the mitochondria. The phenomenon was characterized by comparison of observations in intact cells stimulated by ATP and in plasma membrane permeabilized or in ionophore-treated cells exposed to elevated buffer [Ca2+]. The findings indicate that Ca2+ diffuses laterally within the mitochondria, and that the diffusion is limited for shorter segments of the mitochondrial network. These observations were supported by mathematical simulation of buffered diffusion. The mitochondrial membrane potential was investigated using the potentiometric dye TMRM. Irradiation-induced fluctuations (flickering) of TMRM fluorescence showed synchronicity over large regions of the mitochondrial network, indicating that certain parts of this network form electrical syncytia. The spatial extension of these syncytia was decreased by 2-aminoethoxydiphenyl borate (2-APB) or by propranolol (blockers of nonclassical mitochondrial permeabilities). Our data suggest that mitochondria form syncytia of electrical conductance whereas the passage of Ca2+ is restricted to the individual organelle.     

Organellar vs cellular control of mitochondrial dynamics

Mitochondrial dynamics, the fusion and fission of individual mitochondrial units, is critical to the exchange of the metabolic, genetic and proteomic contents of individual mitochondria. In this regard, fusion and fission events have been shown to modulate mitochondrial bioenergetics, as well as several cellular processes including fuel sensing, ATP production, autophagy, apoptosis, and the cell cycle. Regulation of the dynamic events of fusion and fission occur at two redundant and interactive levels. Locally, the microenvironment of the individual mitochondrion can alter its ability to fuse, divide or move through the cell. Globally, nuclear-encoded processes and cellular ionic and second messenger systems can alter or activate mitochondrial proteins, regulate mitochondrial dynamics and concomitantly change the condition of the mitochondrial population. In this review we investigate the different global and local signals that control mitochondrial biology. This discussion is carried out to clarify the different signals that impact the status of the mitochondrial population.     

From mitochondrial dynamics to arrhythmias

The reactive oxygen species (ROS)-dependent mitochondrial oscillator described in cardiac cells exhibits at least two modes of function under physiological conditions or in response to metabolic and oxidative stress. Both modes depend upon network behavior of mitochondria. Under physiological conditions cardiac mitochondria behave as a network of coupled oscillators with a broad range of frequencies. ROS weakly couples mitochondria under normal conditions but becomes a strong coupling messenger when, under oxidative stress, the mitochondrial network attains criticality. Mitochondrial criticality is achieved when a threshold of ROS is overcome and a certain density of mitochondria forms a cluster that spans the whole cell. Under these conditions, the slightest perturbation triggers a cell-wide collapse of the mitochondrial membrane potential, Δψ_m, visualized as a depolarization wave throughout the cell which is followed by whole cell synchronized oscillations in Δψ_m, NADH, ROS, and GSH. This dynamic behavior scales from the mitochondrion to the cell by driving cellular excitability and the whole heart into catastrophic arrhythmias. A network collapse of Δψ_m under criticality leads to: (i) energetic failure, (ii) temporal and regional alterations in action potential (AP), (iii) development of zones of impaired conduction in the myocardium, and, ultimately, (iv) a fatal ventricular arrhythmia.     

Detection and manipulation of mitochondrial reactive oxygen species in mammalian cells

Reactive oxygen species (ROS) are formed upon incomplete reduction of molecular oxygen (O₂) as an inevitable consequence of mitochondrial metabolism. Because ROS can damage biomolecules, cells contain elaborate antioxidant defense systems to prevent oxidative stress. In addition to their damaging effect, ROS can also operate as intracellular signaling molecules. Given the fact that mitochondrial ROS appear to be only generated at specific sites and that particular ROS species display a unique chemistry and have specific molecular targets, mitochondria-derived ROS might constitute local regulatory signals. The latter would allow individual mitochondria to auto-regulate their metabolism, shape and motility, enabling them to respond autonomously to the metabolic requirements of the cell. In this review we first summarize how mitochondrial ROS can be generated and removed in the living cell. Then we discuss experimental strategies for (local) detection of ROS by combining chemical or proteinaceous reporter molecules with quantitative live cell microscopy. Finally, approaches involving targeted pro- and antioxidants are presented, which allow the local manipulation of ROS levels.     

Ca(2+) homeostasis during mitochondrial fragmentation and perinuclear clustering induced by hFis1

Mitochondria modulate Ca(2+) signals by taking up, buffering, and releasing Ca(2+) at key locations near Ca(2+) release or influx channels. The role of such local interactions between channels and organelles is difficult to establish in living cells because mitochondria form an interconnected network constantly remodeled by coordinated fusion and fission reactions. To study the effect of a controlled disruption of the mitochondrial network on Ca(2+) homeostasis, we took advantage of hFis1, a protein that promotes mitochondrial fission by recruiting the dynamin-related protein, Drp1. hFis1 expression in HeLa cells induced a rapid and complete fragmentation of mitochondria, which redistributed away from the plasma membrane and clustered around the nucleus. Despite the dramatic morphological alteration, hFis1-fragmented mitochondria maintained a normal transmembrane potential and pH and took up normally the Ca(2+) released from intracellular stores upon agonist stimulation, as measured with a targeted ratiometric pericam probe. In contrast, hFis1-fragmented mitochondria took up more slowly the Ca(2+) entering across plasma membrane channels, because the Ca(2+) ions reaching mitochondria propagated faster and in a more coordinated manner in interconnected than in fragmented mitochondria. In parallel cytosolic fura-2 measurements, the capacitative Ca(2+) entry (CCE) elicited by store depletion was only marginally reduced by hFis1 expression. Regardless of mitochondria shape and location, disruption of mitochondrial potential with uncouplers or oligomycin/rotenone reduced CCE by approximately 35%. These observations indicate that close contact to Ca(2+) influx channels is not required for CCE modulation and that the formation of a mitochondrial network facilitates Ca(2+) propagation within interconnected mitochondria.     

Mitochondrial quality control in cardiac cells: Mechanisms and role in cardiac cell injury and disease

Mitochondria play an important role in maintaining cardiac homeostasis by supplying the major energy required for cardiac excitation–contraction coupling as well as controlling the key intracellular survival and death pathways. Healthy mitochondria generate ATP molecules through an aerobic process known as oxidative phosphorylation (OXPHOS). Mitochondrial injury during myocardial infarction (MI) impairs OXPHOS and results in the excessive production of reactive oxygen species (ROS), bioenergetic insufficiency, and contributes to the development of cardiovascular diseases. Therefore, mitochondrial biogenesis along with proper mitochondrial quality control machinery, which removes unhealthy mitochondria is pivotal for mitochondrial homeostasis and cardiac health. Upon damage to the mitochondrial network, mitochondrial quality control components are recruited to segregate the unhealthy mitochondria and target aberrant mitochondrial proteins for degradation and elimination. Impairment of mitochondrial quality control and accumulation of abnormal mitochondria have been reported in the pathogenesis of various cardiac disorders and heart failure. Here, we provide an overview of the recent studies describing various mechanistic pathways underlying mitochondrial homeostasis with the main focus on cardiac cells. In addition, this review demonstrates the potential effects of mitochondrial quality control dysregulation in the development of cardiovascular disease.     

Differential mitochondrial calcium responses in different cell types detected with a mitochondrial calcium fluorescent indicator, mito-GCaMP2

Mitochondrial calcium plays a crucial role in mitochondrial metabolism, cell calcium handling, and cell death. However, some mechanisms concerning mitochondrial calcium regulation are still unknown, especially how mitochondrial calcium couples with cytosolic calcium. In this work, we constructed a novel mitochondrial calcium fluorescent indicator (mito-GCaMP2) by genetic manipulation. Mito-GCaMP2 was imported into mitochondria with high efficiency and the fluorescent signals co-localized with that of tetramethyl rhodamine methyl ester, a mitochondrial membrane potential indicator. The mitochondrial inhibitors specifically decreased the signals of mito-GCaMP2. The apparent K(d) of mito-GCaMP2 was 195.0 nmol/L at pH 8.0 in adult rat cardiomyocytes. Furthermore, we observed that mito-GCaMP2 preferred the alkaline pH surrounding of mitochondria. In HeLa cells, we found that mitochondrial calcium ([Ca(2+)](mito)) responded to the changes of cytosolic calcium ([Ca(2+)](cyto)) induced by histamine or thapasigargin. Moreover, external Ca(2+) (100 μmol/L) directly induced an increase of [Ca(2+)](mito) in permeabilized HeLa cells. However, in rat cardiomyocytes [Ca(2+)](mito) did not respond to cytosolic calcium transients stimulated by electric pacing or caffeine. In permeabilized cardiomyocytes, 600 nmol/L free Ca(2+) repeatedly increased the fluorescent signals of mito-GCaMP2, which excluded the possibility that mito-GCaMP2 lost its function in cardiomyocytes mitochondria. These results showed that the response of mitochondrial calcium is diverse in different cell lineages and suggested that mitochondria in cardiomyocytes may have a special defense mechanism to control calcium flux.     

A phosphodiesterase 2A isoform localized to mitochondria regulates respiration

Mitochondria are central organelles in cellular energy metabolism, apoptosis, and aging processes. A signaling network regulating these functions was recently shown to include soluble adenylyl cyclase as a local source of the second messenger cAMP in the mitochondrial matrix. However, a mitochondrial cAMP-degrading phosphodiesterase (PDE) necessary for switching off this cAMP signal has not yet been identified. Here, we describe the identification and characterization of a PDE2A isoform in mitochondria from rodent liver and brain. We find that mitochondrial PDE2A is located in the matrix and that the unique N terminus of PDE2A isoform 2 specifically leads to mitochondrial localization of this isoform. Functional assays show that mitochondrial PDE2A forms a local signaling system with soluble adenylyl cyclase in the matrix, which regulates the activity of the respiratory chain. Our findings complete a cAMP signaling cascade in mitochondria and have implications for understanding the regulation of mitochondrial processes and for their pharmacological modulation.     

Inhibiting the mitochondrial fission machinery does not prevent Bax/Bak-dependent apoptosis

Apoptosis, induced by a number of death stimuli, is associated with a fragmentation of the mitochondrial network. These morphological changes in mitochondria have been shown to require proteins, such as Drp1 or hFis1, which are involved in regulating the fission of mitochondria. However, the precise role of mitochondrial fission during apoptosis remains elusive. Here we report that inhibiting the fission machinery in Bax/Bak-mediated apoptosis, by down-regulating of Drp1 or hFis1, prevents the fragmentation of the mitochondrial network and partially inhibits the release of cytochrome c from the mitochondria but fails to block the efflux of Smac/DIABLO. In addition, preventing mitochondrial fragmentation does not inhibit cell death induced by Bax/Bak-dependent death stimuli, in contrast to the effects of Bcl-xL or caspase inhibition. Therefore, the fission of mitochondria is a dispensable event in Bax/Bak-dependent apoptosis.     

The anticodon and the D-domain sequences are essential determinants for plant cytosolic tRNA(Val) import into mitochondria

In higher plants, one-third to one-half of the mitochondrial tRNAs are encoded in the nucleus and are imported into mitochondria. This process appears to be highly specific for some tRNAs, but the factors that interact with tRNAs before and/or during import, as well as the signals present on the tRNAs, still need to be identified. The rare experiments performed so far suggest that, besides the probable implication of aminoacyl-tRNA synthetases, at least one additional import factor and/or structural features shared by imported tRNAs must be involved in plant mitochondrial tRNA import. To look for determinants that direct tRNA import into higher plant mitochondria, we have transformed BY2 tobacco cells with Arabidopsis thaliana cytosolic tRNA(Val)(AAC) carrying various mutations. The nucleotide replacements introduced in this naturally imported tRNA correspond to the anticodon and/or D-domain of the non-imported cytosolic tRNA(Met-e). Unlike the wild-type tRNA(Val)(AAC), a mutant tRNA(Val) carrying a methionine CAU anticodon that switches the aminoacylation of this tRNA from valine to methionine is not present in the mitochondrial fraction. Furthermore, mutant tRNAs(Val) carrying the D-domain of the tRNA(Met-e), although still efficiently recognized by the valyl-tRNA synthetase, are not imported any more into mitochondria. These data demonstrate that in plants, besides identity elements required for the recognition by the cognate aminoacyl-tRNA synthetase, tRNA molecules contain other determinants that are essential for mitochondrial import selectivity. Indeed, this suggests that the tRNA import mechanism occurring in plant mitochondria may be different from what has been described so far in yeast or in protozoa.     

A Reaction-Diffusion Model of ROS-Induced ROS Release in a Mitochondrial Network

Loss of mitochondrial function is a fundamental determinant of cell injury and death. In heart cells under metabolic stress, we have previously described how the abrupt collapse or oscillation of the mitochondrial energy state is synchronized across the mitochondrial network by local interactions dependent upon reactive oxygen species (ROS). Here, we develop a mathematical model of ROS-induced ROS release (RIRR) based on reaction-diffusion (RD-RIRR) in one- and two-dimensional mitochondrial networks. The nodes of the RD-RIRR network are comprised of models of individual mitochondria that include a mechanism of ROS-dependent oscillation based on the interplay between ROS production, transport, and scavenging; and incorporating the tricarboxylic acid (TCA) cycle, oxidative phosphorylation, and Ca²⁺ handling. Local mitochondrial interaction is mediated by superoxide (O₂ ^.−) diffusion and the O₂ ^.−-dependent activation of an inner membrane anion channel (IMAC). In a 2D network composed of 500 mitochondria, model simulations reveal ΔΨ_m depolarization waves similar to those observed when isolated guinea pig cardiomyocytes are subjected to a localized laser-flash or antioxidant depletion. The sensitivity of the propagation rate of the depolarization wave to O₂^.− diffusion, production, and scavenging in the reaction-diffusion model is similar to that observed experimentally. In addition, we present novel experimental evidence, obtained in permeabilized cardiomyocytes, confirming that ΔΨ_m depolarization is mediated specifically by O₂ ^.−. The present work demonstrates that the observed emergent macroscopic properties of the mitochondrial network can be reproduced in a reaction-diffusion model of RIRR. Moreover, the findings have uncovered a novel aspect of the synchronization mechanism, which is that clusters of mitochondria that are oscillating can entrain mitochondria that would otherwise display stable dynamics. The work identifies the fundamental mechanisms leading from the failure of individual organelles to the whole cell, thus it has important implications for understanding cell death during the progression of heart disease.     

Analysis of mtDNA copy number and composition of single mitochondrial particles using flow cytometry and PCR

Mitochondrial DNA (mtDNA) is a multicopy, maternally inherited, genome. Individuals frequently carry a mixture of genetically distinct mtDNA molecules whose proportions may vary between sexual generations or among tissues from the same individual. Analyses of the genetic composition of mitochondria have previously relied on electron microscopy and have not permitted the genotype of single mitochondria to be determined. We have developed flow cytometry techniques to isolate single mitochondrial particles and PCR-based assays to determine the mtDNA copy number and composition of individual particles. In a first application of this method, we studied mitochondrial particles from fibroblast cells heteroplasmic for the tRNA lys(8344) point mutation, associated with myoclonus epilepsy and ragged red fiber (MERRF). Individual mitochondrial particles contained between 0 and 11 mtDNA molecules with a mean of 2.0 (95% CI 1.6-2.4). The majority (75%) of the mitochondrial particles from which a PCR product was obtained contained only one type of mtDNA, consistent with the low mean mtDNA copy number. The method developed may be applied to studies of the copy number and distribution of mtDNA genomes in different cell types.     

Inherited complex I deficiency is associated with faster protein diffusion in the matrix of moving mitochondria

Mitochondria continuously change shape, position, and matrix configuration for optimal metabolite exchange. It is well established that changes in mitochondrial metabolism influence mitochondrial shape and matrix configuration. We demonstrated previously that inhibition of mitochondrial complex I (CI or NADH:ubiquinone oxidoreductase) by rotenone accelerated matrix protein diffusion and decreased the fraction and velocity of moving mitochondria. In the present study, we investigated the relationship between inherited CI deficiency, mitochondrial shape, mobility, and matrix protein diffusion. To this end, we analyzed fibroblasts of two children that represented opposite extremes in a cohort of 16 patients, with respect to their residual CI activity and mitochondrial shape. Fluorescence correlation spectroscopy (FCS) revealed no relationship between residual CI activity, mitochondrial shape, the fraction of moving mitochondria, their velocity, and the rate of matrix-targeted enhanced yellow fluorescent protein (mitoEYFP) diffusion. However, mitochondrial velocity and matrix protein diffusion in moving mitochondria were two to three times higher in patient cells than in control cells. Nocodazole inhibited mitochondrial movement without altering matrix EYFP diffusion, suggesting that both activities are mutually independent. Unexpectedly, electron microscopy analysis revealed no differences in mitochondrial ultrastructure between control and patient cells. It is discussed that the matrix of a moving mitochondrion in the CI-deficient state becomes less dense, allowing faster metabolite diffusion, and that fibroblasts of CI-deficient patients become more glycolytic, allowing a higher mitochondrial velocity.     

Agonist-evoked mitochondrial Ca2+ signals in mouse pancreatic acinar cells

In the present study we have investigated cytosolic and mitochondrial Ca(2+) signals in isolated mouse pancreatic acinar cells double-loaded with the fluorescent probes fluo-3 and rhod-2. Stimulation of pancreatic acinar cells with 500 nm acetylcholine caused release of Ca(2+) from intracellular stores and produced cytosolic Ca(2+) signals in form of Ca(2+) waves propagating from the luminal to the basal cell pole. The increase in the cytosolic Ca(2+) concentration was followed by Ca(2+) uptake into mitochondria. Between onset of cytosolic and mitochondrial Ca(2+) signals there was a delay of 10.7 +/- 0.4 s. Ca(2+) uptake into mitochondria could be inhibited with Ruthenium Red and carbonyl cyanide m-chlorophenylhydrazone, whereas 2,5-di-tert-butylhydroquinone, which inhibits sarco(endo)plasmic reticulum Ca(2+) ATPases, did not prevent Ca(2+) accumulation in mitochondria. Carbonyl cyanide m-chlorophenylhydrazone-induced Ca(2+) release from mitochondria could only be observed after a preceding stimulation of the cell with a physiological agonist or by treatment with 2, 5-di-tert-butylhydroquinone, indicating that under resting conditions mitochondria do not contain releasable Ca(2+) ions. Analysis of the propagation rate of acetylcholine-induced Ca(2+) waves revealed that inhibition of mitochondrial Ca(2+) uptake did not accelerate spreading of cytosolic Ca(2+) signals. Our experiments indicate that in the early phase of secretagogue-induced Ca(2+) signals, mitochondria behave as passive Ca(2+)-buffering elements and do not actively suppress spreading of Ca(2+) signals in pancreatic acinar cells.     

Human mitochondria and mitochondrial genome function as a single dynamic cellular unit

rho 0 HeLa cells entirely lacking mitochondrial DNA (mtDNA) and mitochondrial transfection techniques were used to examine intermitochondrial interactions between mitochondria with and without mtDNA, and also between those with wild-type (wt) and mutant-type mtDNA in living human cells. First, unambiguous evidence was obtained that the DNA-binding dyes ethidium bromide (EtBr) and 4',6-diamidino-2- phenylindole (DAPI) exclusively stained mitochondria containing mtDNA in living human cells. Then, using EtBr or DAPI fluorescence as a probe, mtDNA was shown to spread rapidly to all rho 0 HeLa mitochondria when EtBr- or DAPI-stained HeLa mitochondria were introduced into rho 0 HeLa cells. Moreover, coexisting wt-mtDNA and mutant mtDNA with a large deletion (delta-mtDNA) were shown to mix homogeneously throughout mitochondria, not to remain segregated by use of electron microscopic analysis of cytochrome c oxidase activities of individual mitochondria as a probe to identify mitochondria with predominantly wt- or delta- mtDNA in single cells. This rapid diffusion of mtDNA and the resultant homogeneous distribution of the heteroplasmic wt- and delta-mtDNA molecules throughout mitochondria in a cell suggest that the mitochondria in living human cells have lost their individuality. Thus, the actual number of mitochondria per cell is not of crucial importance, and mitochondria in a cell should be considered as a virtually single dynamic unit.     

Mathematical Modeling of the Role of Mitochondrial Fusion and Fission in Mitochondrial DNA Maintenance

Accumulation of mitochondrial DNA (mtDNA) mutations has been implicated in a wide range of human pathologies, including neurodegenerative diseases, sarcopenia, and the aging process itself. In cells, mtDNA molecules are constantly turned over (i.e. replicated and degraded) and are also exchanged among mitochondria during the fusion and fission of these organelles. While the expansion of a mutant mtDNA population is believed to occur by random segregation of these molecules during turnover, the role of mitochondrial fusion-fission in this context is currently not well understood. In this study, an in silico modeling approach is taken to investigate the effects of mitochondrial fusion and fission dynamics on mutant mtDNA accumulation. Here we report model simulations suggesting that when mitochondrial fusion-fission rate is low, the slow mtDNA mixing can lead to an uneven distribution of mutant mtDNA among mitochondria in between two mitochondrial autophagic events leading to more stochasticity in the outcomes from a single random autophagic event. Consequently, slower mitochondrial fusion-fission results in higher variability in the mtDNA mutation burden among cells in a tissue over time, and mtDNA mutations have a higher propensity to clonally expand due to the increased stochasticity. When these mutations affect cellular energetics, nuclear retrograde signalling can upregulate mtDNA replication, which is expected to slow clonal expansion of these mutant mtDNA. However, our simulations suggest that the protective ability of retrograde signalling depends on the efficiency of fusion-fission process. Our results thus shed light on the interplay between mitochondrial fusion-fission and mtDNA turnover and may explain the mechanism underlying the experimentally observed increase in the accumulation of mtDNA mutations when either mitochondrial fusion or fission is inhibited.