- Department of Molecular Energetics of Microorganisms, Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia
Introduction
Individual cells within a clonal population can display pronounced phenotypic differences. In multicellular organisms, such differentiation is controlled by numerous microenvironmental, positional, and genetic cues. However, even unicellular Saccharomyces cerevisiae cells show significant levels of cell-to-cell variability. This variability arises from de novo mutations, gene expression noise and asymmetric division. Here, I suggest that another contributing factor could be genetic drift in mitochondrial DNA (mtDNA) heteroplasmy levels. Moreover, because of the activation of general stress response pathways in cells with insufficient levels of wild type mtDNA, heteroplasmic cells can withstand certain stresses better than wild type cells. As a result, the heterogeneity within the heteroplasmic cell population could be a bet-hedging mechanism, increasing the chances of cell line survival.
Main text
Eukaryotic cells reproduce asexually most of the time. In the species with obligate sexual reproduction, there are thousands of mitotic divisions per meiotic division. For instance, human body comprises 3*1013 clonal somatic cells (Sender et al., 2016) that correspond to at average 13.4 divisions of zygote descendants. The actual number of cell generations is even higher due to cell turnover. Moreover, certain species have facultative sexual reproduction, so baker’s yeast cells go through a sexual reproduction phase only once in every ∼10,000 generations (Nieuwenhuis and James, 2016). Therefore, a majority of the living eukaryotic cells have multiple clonal copies of their own.
Clonal cells, despite their genetic identity, can pronouncedly differ from each other by the phenotype. This is true, even in the case of unicellular organisms, such as baker’s yeast S. cerevisiae. This is evidenced by differing protein concentrations and stress response readiness in genetically identical cells (Levy et al., 2012). The driving forces behind this heterogeneity include gene expression noise, cyclic processes, and asymmetry in cell divisions (Knorre et al., 2018). Recent research has shown that molecular stochasticity of 5′-UTR scanning by ribosome contributes to yeast population heterogeneity in cell response to nutrient limitation (Meng et al., 2023). In certain cases, yeast cells suspensions exhibit phenotypic bistability which is characterised by the coexistence of cells with two distinguishable phenotypes in a clonal cell line. For example, introns mediate bimodal distribution of the yeast ribosomal protein Rps22B in yeast cells (Lukačišin et al., 2022).
Cellular regulatory systems serve to minimise phenotypic variation. Therefore, gene promoter regions are subject to purifying selection pressure which mitigates gene expression noise (Metzger et al., 2015). Indeed, an elevation in phenotypic variance means that a growing percentage of cells within a subpopulation deviates from an optimal phenotype. However, in some cases there is no single optimal phenotypic value, particularly when the optimal protein concentration within a cell depends on fluctuating environmental conditions. In particular, the growth rate of yeast negatively correlates with its cellular resistance to stress (Ho and Gasch, 2015). At the same time, a genetically homogeneous suspension culture of yeast is usually heterogeneous in the rate of division of individual cells (Arabaciyan et al., 2021) where slow-dividing cells upregulate general stress-response transcription factors Msn2p/Msn4p and exhibit a higher trehalose synthase Tsl1 concentration (Li et al., 2018). The high trehalose synthase Tsl1p concentration is beneficial during heat stress but associated with decreased growth rate under normal conditions (Levy et al., 2012). We have also recently demonstrated that post-diauxic yeast cells form two subpopulations characterised by high and low concentrations of mitochondrial ATP-ase inhibitor proteins (Inh1p and Stf1p). Cells with an elevated inhibitor level recovered more rapidly from the stationary phase than their low-level counterparts, at the same time a high concentration of the inhibitor proteins was deleterious under the conditions of mitochondrial dysfunction (Galkina et al., 2022). These examples show that cell culture phenotypic heterogeneity can be detrimental or adaptive depending on which specific phenotypic properties are considered.
Eukaryotic cells harbour mitochondria, the cellular organelle with its own genome; refer to (Janouškovec et al., 2017) for the rare exceptions. The number of mtDNA molecules in cells varies depending on conditions, cell type and the previous cell’s history (Medeiros et al., 2018; Zhang et al., 2021; Galeota-Sprung et al., 2022). Meanwhile, heteroplasmy, a condition characterised by the presence of multiple mtDNA variants within a single cell, may occur as a result of mutations or cell fusion (Figure). The heteroplasmy is a common phenomenon. Many human genetic mitochondrial disorders are linked to the coexistence of a wild type mtDNA molecule and a variant bearing a deleterious mutation (Wallace and Chalkia, 2013). Furthermore, heteroplasmy in germline can be transmitted for several generations (Burgstaller et al., 2018).
Several forces shape mtDNA allele frequency within a cell line (Rand, 2001). 1) The random sampling of mtDNA during cell division makes the relative number of two mitochondrial genotypes drift over generations (Figure). The drift continues until one of the variants becomes fixed. Given that mtDNA can replicate under cell cycle arrest (Newlon and Fangman, 1975), a proportion of different mtDNAs can also change over time in post-mitotic heteroplasmic cells (Taylor and Turnbull, 2005). 2) At the same time, dividing cells with a high proportion of the pathogenic mtDNA molecules become unable to proliferate and, therefore, they are outcompeted by other cells with functional mitochondria. As accumulation of mutant mtDNA in a cell has a threshold effect on a cell phenotype (Rossignol et al., 2003), pathogenic mtDNA variants can persist in a cell line at a low proportion for an extended period. The pathogenicity threshold usually lies in the range 50%–90% of mutant mtDNA allele frequency; the number depends on cell type and the nature of mutation (Rossignol et al., 2003). This suggests that the wild type mtDNA is usually present in a moderate excess in normal cells. 3) Due to their superior intracellular fitness, certain deleterious mtDNAs variants can displace wild type variants within the cell (Gitschlag et al., 2020).
These three factors contribute to the ability of harmful mtDNA variants to persist in clonal populations along with wild type mtDNAs. Such populations should sporadically produce cells with a high proportion of mutant mtDNA and dysfunctional mitochondria. Indeed, multicellular organisms develop from a single cell, which experiences a significant genetic bottleneck in mtDNA quantity. This bottleneck decreases the diversity of mtDNA variants (Krakauer and Mira, 1999; Zaidi et al., 2019). Despite this, the proportion of respiration-deficient cells in some tissues increases with age (Taylor and Turnbull, 2005; Payne et al., 2011; Hipps et al., 2022).
Eukaryotic cells are equipped with the mechanisms capable of compensating deleterious effects of mitochondrial dysfunction. In nematodes, for instance, the depletion of mitochondrial DNA triggers the expression of mitochondrial chaperones (Nargund et al., 2012). In addition, mitochondrial dysfunction induces the general stress response mechanisms, which encompass both mitochondrial and cytosolic heat-shock response proteins (Poveda-Huertes et al., 2020). Furthermore, yeast cells with dysfunctional mitochondria demonstrate a surprising resilience to certain stress factors, outperforming cells with functional oxidative phosphorylation. Such yeast cells display enhanced resistance to a range of stressors, including the protein synthesis inhibitor cycloheximide (Zhang and Moye-Rowley, 2001), cytoplasm acidification (Ludovico et al., 2002), mating pheromone-induced death (Severin and Hyman, 2002), lipophilic cations (Antonenko et al., 2011), and the prooxidant paraquat (Stenberg et al., 2022). The resistance to these stressors is facilitated, in part, by the activation of pleiotropic drug resistance (PDR) transporter genes in yeast cells without functional mtDNA. These genes encode ABC-transporters with broad substrate specificity capable of effluxing xenobiotics from the cytoplasm into the environment (Zhang and Moye-Rowley, 2001).
Taken together, it turns out that, on the one hand, in clonal populations mitochondrial heteroplasmy causes the sporadic occurrence of the cells with impaired mitochondrial function; and mitochondrial dysfunction activates protective mechanisms. On the other hand, the activation of such a mechanism can be useful for cells in a wide range of stressful conditions. Importantly, in a heteroplasmic cell with a high load of mutant mtDNA but retaining wild-type mtDNA the proportion of full-length wild-type mtDNA can be increased over several generations in the absence of stress. Therefore, mitochondrial heteroplasmy can be a natural way to generate phenotypic bistability in clonal populations. I suggest that bistability mediated by mitochondrial heteroplasmy can be a mechanism betting the risk in the fluctuating environment (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1276629/fcell-11-1276629-HTML/image_m/fcell-11-1276629-g001.jpg)
FIGURE 1. Mitochondrial DNA heteroplasmy provides bistability in mitochondrial dysfunction phenotypes. Random mutations and genetic drift generate a fraction of cells with dysfunctional mitochondria (A). These cells are outcompeted by the cells with low burden of mutant mtDNA under normal conditions (B), but better survive under certain stresses (C), and wild type mtDNA can repopulate cells again when the stress is relieved. (D) Dotted arrows indicate the direction of change in mtDNA allele frequencies.
Author contributions
DK: Conceptualization, Visualization, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Russian Science Foundation (project 22-14-00108).
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Antonenko, Y. N., Avetisyan, A. V., Cherepanov, D. A., Knorre, D. A., Korshunova, G. A., Markova, O. V., et al. (2011). Derivatives of rhodamine 19 as mild mitochondria-targeted cationic uncouplers. J. Biol. Chem. 286, 17831–17840. doi:10.1074/jbc.M110.212837
Arabaciyan, S., Saint-Antoine, M., Maugis-Rabusseau, C., François, J.-M., Singh, A., Parrou, J.-L., et al. (2021). Insights on the control of yeast single-cell growth variability by members of the trehalose phosphate synthase (TPS) complex. Front. Cell Dev. Biol. 9, 607628. doi:10.3389/fcell.2021.607628
Burgstaller, J. P., Kolbe, T., Havlicek, V., Hembach, S., Poulton, J., Piálek, J., et al. (2018). Large-scale genetic analysis reveals mammalian mtDNA heteroplasmy dynamics and variance increase through lifetimes and generations. Nat. Commun. 9, 2488. doi:10.1038/s41467-018-04797-2
Galeota-Sprung, B., Fernandez, A., and Sniegowski, P. (2022). Changes to the mtDNA copy number during yeast culture growth. R. Soc. Open Sci. 9, 211842. doi:10.1098/rsos.211842
Galkina, K., Zubareva, V., Kashko, N., Lapashina, A., Markova, O., Feniouk, B., et al. (2022). Heterogeneity of starved yeast cells in IF1 levels suggests the role of this protein in vivo. Front. Microbiol. 13, 816622. doi:10.3389/fmicb.2022.816622
Gitschlag, B. L., Tate, A. T., and Patel, M. R. (2020). Nutrient status shapes selfish mitochondrial genome dynamics across different levels of selection. Elife 9, e56686. doi:10.7554/eLife.56686
Hipps, D., Dobson, P. F., Warren, C., McDonald, D., Fuller, A., Filby, A., et al. (2022). Detecting respiratory chain defects in osteoblasts from osteoarthritic patients using imaging mass cytometry. Bone 158, 116371. doi:10.1016/j.bone.2022.116371
Ho, Y.-H., and Gasch, A. P. (2015). Exploiting the yeast stress-activated signaling network to inform on stress biology and disease signaling. Curr. Genet. 61, 503–511. doi:10.1007/s00294-015-0491-0
Janouškovec, J., Tikhonenkov, D. V., Burki, F., Howe, A. T., Rohwer, F. L., Mylnikov, A. P., et al. (2017). A new lineage of eukaryotes illuminates early mitochondrial genome reduction. Curr. Biol. 27, 3717–3724. doi:10.1016/j.cub.2017.10.051
Knorre, D. A., Azbarova, A. V., Galkina, K. V., Feniouk, B. A., and Severin, F. F. (2018). Replicative aging as a source of cell heterogeneity in budding yeast. Mech. Ageing Dev. 176, 24–31. doi:10.1016/j.mad.2018.09.001
Krakauer, D. C., and Mira, A. (1999). Mitochondria and germ-cell death. Nature 400, 125–126. doi:10.1038/22026
Levy, S. F., Ziv, N., and Siegal, M. L. (2012). Bet hedging in yeast by heterogeneous, age-correlated expression of a stress protectant. PLoS Biol. 10, e1001325. doi:10.1371/journal.pbio.1001325
Li, S., Giardina, D. M., and Siegal, M. L. (2018). Control of nongenetic heterogeneity in growth rate and stress tolerance of Saccharomyces cerevisiae by cyclic AMP-regulated transcription factors. PLoS Genet. 14, e1007744. doi:10.1371/journal.pgen.1007744
Ludovico, P., Rodrigues, F., Almeida, A., Silva, M. T., Barrientos, A., and Côrte-Real, M. (2002). Cytochrome c release and mitochondria involvement in programmed cell death induced by acetic acid in Saccharomyces cerevisiae. Mol. Biol. Cell 13, 2598–2606. doi:10.1091/mbc.e01-12-0161
Lukačišin, M., Espinosa-Cantú, A., and Bollenbach, T. (2022). Intron-mediated induction of phenotypic heterogeneity. Nature 605, 113–118. doi:10.1038/s41586-022-04633-0
Medeiros, T. C., Thomas, R. L., Ghillebert, R., and Graef, M. (2018). Autophagy balances mtDNA synthesis and degradation by DNA polymerase POLG during starvation. J. Cell Biol. 217, 1601–1611. doi:10.1083/jcb.201801168
Meng, X., Reed, A., Lai, S., Szavits-Nossan, J., and McCarthy, J. E. G. (2023). Stochastic scanning events on the GCN4 mRNA 5’ untranslated region generate cell-to-cell heterogeneity in the yeast nutritional stress response. Nucleic Acids Res. 51, 6609–6621. doi:10.1093/nar/gkad433
Metzger, B. P. H., Yuan, D. C., Gruber, J. D., Duveau, F., and Wittkopp, P. J. (2015). Selection on noise constrains variation in a eukaryotic promoter. Nature 521, 344–347. doi:10.1038/nature14244
Nargund, A. M., Pellegrino, M. W., Fiorese, C. J., Baker, B. M., and Haynes, C. M. (2012). Mitochondrial import efficiency of ATFS-1 regulates mitochondrial UPR activation. Science 337, 587–590. doi:10.1126/science.1223560
Newlon, C. S., and Fangman, W. L. (1975). Mitochondrial DNA synthesis in cell cycle mutants of Saccharomyces cerevisiae. Cell 5, 423–428. doi:10.1016/0092-8674(75)90061-6
Nieuwenhuis, B. P. S., and James, T. Y. (2016). The frequency of sex in fungi. Philos. Trans. R. Soc. Lond. B Biol. Sci. 371, 20150540. doi:10.1098/rstb.2015.0540
Payne, B. A. I., Wilson, I. J., Hateley, C. A., Horvath, R., Santibanez-Koref, M., Samuels, D. C., et al. (2011). Mitochondrial aging is accelerated by anti-retroviral therapy through the clonal expansion of mtDNA mutations. Nat. Genet. 43, 806–810. doi:10.1038/ng.863
Poveda-Huertes, D., Matic, S., Marada, A., Habernig, L., Licheva, M., Myketin, L., et al. (2020). An early mtUPR: redistribution of the nuclear transcription factor Rox1 to mitochondria protects against intramitochondrial proteotoxic aggregates. Mol. Cell 77, 180–188. doi:10.1016/j.molcel.2019.09.026
Rand, D. M. (2001). The units of selection on mitochondrial DNA. Annu. Rev. Ecol. Syst. 32, 415–448. doi:10.1146/annurev.ecolsys.32.081501.114109
Rossignol, R., Faustin, B., Rocher, C., Malgat, M., Mazat, J.-P., and Letellier, T. (2003). Mitochondrial threshold effects. Biochem. J. 370, 751–762. doi:10.1042/BJ20021594
Sender, R., Fuchs, S., and Milo, R. (2016). Revised estimates for the number of human and bacteria cells in the body. PLoS Biol. 14, e1002533. doi:10.1371/journal.pbio.1002533
Severin, F. F., and Hyman, A. A. (2002). Pheromone induces programmed cell death in S. cerevisiae. Curr. Biol. 12, R233–R235. doi:10.1016/s0960-9822(02)00776-5
Stenberg, S., Li, J., Gjuvsland, A. B., Persson, K., Demitz-Helin, E., González Peña, C., et al. (2022). Genetically controlled mtDNA deletions prevent ROS damage by arresting oxidative phosphorylation. Elife 11, e76095. doi:10.7554/eLife.76095
Taylor, R. W., and Turnbull, D. M. (2005). Mitochondrial DNA mutations in human disease. Nat. Rev. Genet. 6, 389–402. doi:10.1038/nrg1606
Wallace, D. C., and Chalkia, D. (2013). Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease. Cold Spring Harb. Perspect. Biol. 5, a021220. doi:10.1101/cshperspect.a021220
Zhang, X., and Moye-Rowley, W. S. (2001). Saccharomyces cerevisiae multidrug resistance gene expression inversely correlates with the status of the F(0) component of the mitochondrial ATPase. J. Biol. Chem. 276, 47844–47852. doi:10.1074/jbc.M106285200
Zaidi, A. A., Wilton, P. R., Su, M. S.-W., Paul, I. M., Arbeithuber, B., Anthony, K., et al. (2019). Bottleneck and selection in the germline and maternal age influence transmission of mitochondrial DNA in human pedigrees. Proc. Natl. Acad. Sci. U. S. A. 116, 25172–25178. doi:10.1073/pnas.1906331116
Keywords: yeast, heterogeneity, heteroplasmy, pathogenicity threshold, mtDNA, retrograde signaling, mitochondrial dysfuncion
Citation: Knorre DA (2023) Mitochondrial heteroplasmy as a cause of cell-to-cell phenotypic heterogeneity in clonal populations. Front. Cell Dev. Biol. 11:1276629. doi: 10.3389/fcell.2023.1276629
Received: 29 August 2023; Accepted: 26 September 2023;
Published: 09 October 2023.
Edited by:
Dominic C. Voon, Kanazawa University, JapanReviewed by:
Conor Lawless, Newcastle University, United KingdomCopyright © 2023 Knorre. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dmitry A. Knorre, knorre@belozersky.msu.ru