ARTICLE: Cellular energy budgets and antimicrobial resistance

Dynamic Boolean modelling reveals the influence of energy supply on bacterial efflux pump expression
Ryan Kerr, Sara Jabbari, Jessica MA Blair, Iain G Johnston
Journal of the Royal Society Interface 19 20210771 (2021)

Antimicrobial resistance or AMR is a major global health issue, with disease-causing organisms like bacteria acquiring resistance to the drugs we use to kill them. One way that bacteria acquire this resistance is through so-called efflux pumps -- cellular machinery that removes chemicals like drugs from inside the bacterium. Bacteria produce these pumps when faced with drug treatments, but not all cells produce the same amount or at the same time. Understanding this variability could help the theory behind future treatments.

After finding that the available "energy budget" influences the behaviour of cellular programs, we asked whether energy variability could be a cause of these differences. Using lots of diverse experimental observations, we built a theoretical model of the signals that tell a bacterium to produce efflux pumps in response to sensing a drug, with a new and simple way of accounting for how energy affects these signals. We then simulated this model in a computer to see how model cells with different amounts of available energy (as we see in real bacterial populations) behaved.

We found that differences in cellular energy budgets can have a profound effect on when, and how much, efflux machinery is produced. This variability builds on the natural randomness of the system, leading to several interesting results: energy changes the dynamics of how signalling programs work in the cell, alters timescales, and affects the "priming" of a population of cells to anticipate future stress. The approach we developed is quite general and can be used to explore energy influence on any other cell programs and signals too.

Including ATP, an important cellular energy currency, in models of how bacteria express efflux machinery helps us understand how cell-to-cell differences in energy budget may influence AMR.

ARTICLE: What makes mitochondria selfish, and when do selfish ones win?

MtDNA sequence features associated with ‘selfish genomes’ predict tissue-specific segregation and reversion, Nucleic Acids Research 48 8290 (202)

Mitochondria, the power stations of the cell, are in some senses like people in a company. The company needs several people to contribute if it is to survive. Some people may work hard and contribute lots to the company. Others may selfishly slack off and rely on others doing the work.

The cell needs mitochondria to produce ATP, the chemical that powers many important processes. But there is evidence that some mitochondria are more selfish, and some less so, than others. Unselfish mitochondria produce machinery which helps produce ATP. Selfish mitochondria prefer to replicate, copying their DNA and contributing less to the cell. Interestingly, at the molecular level, there is something like a "switch": a mitochondrion either takes steps to produce useful machinery, or takes steps that will help it replicate.

We were interested in why different mitochondria choose different positions of this switch, and what is a good "strategy" for mitochondria under different conditions. We built a simple model of this behaviour to understand it. Unsurprisingly, we found that selfish mitochondria -- favouring replication, and contributing less to the cell -- profilerate over unselfish ones in cells where there's little pressure to co-operate. As they replicate more, selfish mitochondria eventually come to dominate such cells. Where there is cellular pressure, however, unselfish mitochondria may win out. This is because cells full of selfish mitochondria won't perform adequately, and the whole cell and all its mitochondria will die -- leaving those cells with more unselfish mitochondria remaining.

What do we mean by "cellular pressure"? Well, if some type of cells never die, clearly the latter event can't happen, and we'd expect selfish mitochondria to win. If cells die regularly, perhaps there's more capacity to select those filled with unselfish mitochondria. We looked at different tissues where cells die with different rates, in mice where cells had two different types of mitochondrial DNA (mtDNA). We found a consistent pattern where one type of mtDNA proliferated in slow-dying cells and the other proliferated in fast-dying cells. Why different mtDNA types "win" in different tissues is a big question (which we've looked at before!), and it looks like this might help explain some of these differences.

(A) Sequence features may make different mtDNA types more "selfish" (favouring replication) or "unselfish" (favouring the production of useful machinery). (B) Our theory shows how, depending on cellular pressures, one or the other strategy can be favoured, leading to selection for one or the other mtDNA type.

We also asked what it is about a particular mtDNA sequence that might make it more or less selfish. Based on how mtDNA produces useful machinery, and how it replicates, we hypothesised that some features in the so-called "control region" of mtDNA may influence selfishness. Using sequence information, we found that these features tied quite neatly in with the observations in these mouse models, and also in (more limited) observations from human cells. While certainly not resolved, this picture suggests a link between sequence features of mtDNA, cellular selfishness, and proliferation differences across different tissues. You can read more (for free) in Nucleic Acids Research here.

ARTICLE: Powering cellular decision-making

Intracellular energy Variability Modulates cellular Decision-Making capacity
Ryan Kerr, Sara Jabbari, Iain G Johnston
Scientific Reports 9 1 (2019)

The ability to process information and make decisions is fundamental to life. Intelligent organisms use their brains to do this, but individual cells are also constantly making decisions, changing their behaviour in response to microscopic stimuli. Examples of this cellular decision-making abound in biology: stem cells decide which type of cell to become; some bacteria decide to become robust "persister" cells that can survive drug treatments; cells in plant seeds decide when to germinate.

Often, the "decisions" that cells make involve which genes to express. Genes contain information on how to build cellular machinery, and "expressing" a gene in a sense means turning it on so that its machinery gets built in the cell. We often see that two genes, say A and B, build proteins that switch each other's genes off. So if we have lots of A, it's very hard to produce B, and vice versa. These genes can determine the type of cell we have -- for example, cells with lots of A might be white blood cells, and cells with lots of B might be red blood cells. A blood stem cell could then become a white or a red cell depending on how the interaction between A and B plays out.

All this is reasonably common knowledge (though rather simplified!). But we got interested in how energy plays a role in these decisions. Gene expression requires energy, which in the cell is provided by a molecule called ATP. Different cells have different amounts of ATP, so the processes involved in the interaction of our genes A and B can take place at different rates. Following some ideas we laid out here, we asked, using maths, how this energy dependence might affect the decisions that cells make.

We found, in a new paper free to read in Scientific Reports, that ATP levels strongly influence the decision-making capacity of a cell. Consider the simple A-B case above. Four states are possible: no A or B (state 0), more A than B (state A), more B than A (state B), and high A and B (state AB). We found that, at low ATP, only state 0 is possible (the cell can't make any decisions). As ATP increases, states A and B become possible, and for high ATP the state AB also appears. So, the number of states a cell can choose between (for example, white, red, or stem blood cell) depends strongly on how much energy that cell has available to power these genetic interactions.

We also found that more energy stabilised the decisions that could be made (cells are noisy, so decisions can be randomly "overturned" by gene expression fluctuations), and mapped out the "landscape" of decisions that can be made as the biochemical features of the genes involved change. We're now going to the lab to explore these mathematical predictions in real cells -- particularly in bacterial persister cells -- and developing the theory further for more complicated decision-making circuits.

ARTICLE: The inheritance of mtDNA

Regulation of mother-to-offspring transmission of mtDNA heteroplasmy
Ana Latorre-Pellicer, Ana Victoria Lechuga-Vieco, Iain G Johnston, Riikka H Hämäläinen, Juan Pellico, Raquel Justo-Méndez, Jose María Fernández-Toro, Cristina Clavería, Adela Guaras, Rocío Sierra, Jordi Llop, Miguel Torres, Luis Miguel Criado, Anu Suomalainen, Nick S Jones, Jesús Ruíz-Cabello, José Antonio Enríquez
Cell Metabolism 30 1120 (2019)

Mitochondrial DNA (mtDNA) is inherited from mothers to children. If two or more types of mtDNA exist in a cell, the cell is called "heteroplasmic". Mothers may carry a heteroplasmic mixture of different mtDNA types in each of their oocytes (egg cells), so a mixture of different types may be passed on to children. Different oocytes may have different mixtures -- for example, one cel may have 50% type A and 50% type B, and another may have 70% A and 30% B. 

The inheritance of heteroplasmy depends both on a complicated "bottleneck" (see here and here) and whether either type has some advantage over the other -- a question that is hotly debated. The mechanisms that shape the inheritance of mtDNA populations remain poorly understood, so it's hard to predict which offspring will inherit which mixture. It's often the case that a disease is caused by a particular mixture -- for example, over 60% of type B -- so this complex inheritance makes it hard to plan fertility treatments too.

In a new paper in Cell Metabolism, we looked at the inheritance and consequences of heteroplasmy in mice. Strikingly, we found that the presence of any heteroplasmy has generally negative consequences for the cell. This is perhaps surprising, given the above view that we normally need a certain amount of a dangerous mtDNA type to cause disease. But it does match a prediction that we recently made by mathematically considering how a cell must invest energy in controlling mixed mtDNA populations. Correspondingly, we found that regardless of how much type A and type B there is, having a heteroplasmic mixture challenges metabolism in the embryo, and affects how readily induced pluripotent stem cells (iPSCs) can be produced from cells.

Given that heteroplasmy is a challenge, it seems that cells have evolved mechanisms to sense and address the inheritance of heteroplasmy. In addition to the bottleneck, we found (as in our previous work) that cell-to-cell variance of heteroplasmy increased in oocytes with age -- which will have the eventual effect of reducing heteroplasmy. We also found that particular mtDNA types had a selective advantage through inheritance, and identified a set of genes that shape this advantage. Variability in the expression of these genes, and variability in metabolic factors, led to differences in the strength of selection.

Some key findings from this paper, and links to our previous work.

This work was exciting because it provided some insights into the mechanisms that shape mtDNA populations between generations -- but also because it validated several predictions that our theoretical work had proposed in the past:

• Increasing heteroplasmy variance with age (predicted here, observed here)
• MtDNA selection occurs at different developmental stages (as we found here and here)
• Mixed mtDNA populations challenge the cell (predicted here)
• Genes related to mitochondrial dynamics shape mtDNA genetic makeup (predicted here)

We're continuing this exciting collaboration and looking in more depth at the behaviour of mtDNA over time.

ARTICLES: Evolving cellular populations of mtDNA

Evolving mtDNA populations within cells
Iain G Johnston, Joerg P Burgstaller
Biochemical Society Transactions 47 1367 (2019)

and

Varied mechanisms and models for the varying mitochondrial bottleneck
Iain G Johnston
Frontiers in Cell and Developmental Biology 7 294 (2019)

We've recently written two review papers looking at the dynamics of mitochondrial DNA (mtDNA) in cells. As we've written about before, cells contain populations of hundreds or thousands of mtDNA molecules. These molecules replicate and degrade, so that over time, cellular populations of mtDNA change and evolve. The amount of disease-causing mutations, and the number and structure of mtDNA molecules, may all change as organisms develop and age, with different consequences.

The first article, in Biochemical Society Transactions, takes a broad look at how cellular mtDNA populations change over time, considering a range of organisms from humans and other animals to plants and fungi. We look at the different processes that change mtDNA populations, which include replication and degradation but may also include recombination (particularly in plants) and cell-to-cell exchange of mitochondria. The review particularly highlights the importance of understanding cell-to-cell variability in mtDNA populations -- as it only takes a few cells with lots of mutant mtDNA to cause disease, it's important to understand the statistics of mtDNA populations across cells. We review experiments and theory aiming to do so, including our recent work showing that cell-to-cell variability of mtDNA mutant load increases over time in a wide variety of circumstances.

The second article, in Frontiers in Cell and Developmental Biology, focuses on the so-called "mtDNA bottleneck", a process that shapes mtDNA populations in early mammalian development, and helps prevent the inheritance of mutant mtDNA. Specifically, mtDNA undergoes a "genetic bottleneck" between generations, meaning that mothers' egg cells, and offspring, often have dramatically different mtDNA populations. The review emphasises that this "genetic bottleneck" is an effective quantity, not a directly measurable observation, that arises from several physical processes, including but not limited to a "physical bottleneck" or mtDNA depletion during development. Different ways of modelling, analysing, and explaining the "genetic bottleneck" are reviewed, from human populations to mouse egg cells and Adélie penguins. We invest some time in trying to explain the different assumptions, symbols, and methods that researchers have used to quantify the bottleneck over the years. Again, the importance of understanding cell-to-cell variance in mtDNA populations is a core theme.

A. Different processes shaping mixed mtDNA populations inside cells. B. The "genetic bottleneck", increasing mtDNA variance between egg cells and offspring.  

Like all reviews, these articles don't have new results, but attempt to summarise existing knowledge and thinking on these topics. We hope that both papers provide some interesting insights, references, and (in the case of the bottleneck paper) visualisations that may help understand these often confusing topics.

ARTICLE: The cell's power station policies

Energetic costs of cellular and therapeutic control of stochastic mitochondrial DNA populations

Hanne Hoitzing, Payam A Gammage, Lindsey van Haute, Michal Minczuk, Iain G Johnston, Nick S Jones

PLoS Computational Biology 15 e1007023 (2019)

(Hanne's also written a post about this paper, you can read it here)

Our cells are filled with populations of mitochondrial DNA (mtDNA) molecules, which encode vital cellular machinery that supports our energy requirements. The cell invests energy in maintaining its mtDNA population, like us using electricity-powered tools to help maintain our power stations. Our cellular power stations can vary in quality (for example, mutations can damage mtDNA), and are subject to random influences. How should the cell best invest energy in controlling and maintaining its power stations? And can we use this answer to design better therapies to address damaged mtDNA?

In a new paper here in PLoS Computational Biology, we attempt to answer this question using mathematical modelling, linking with genetic experiments done by our excellent collaborators at Cambridge (Payam Gammage, Lindsey Van Haute and Michal Minczuk). We first expand a mathematical model for how diverse mtDNA populations within cells change over time – building new power stations and decommissioning old ones, under the “governance” of the cell. We then produce an “energy budget” for the cellular “society” – describing the costs of building, decommissioning, and maintaining different power stations, and the corresponding profits of energy generation.

We find some surprising results. First, it can get harder to maintain a good energy budget in a tissue (a collection of individual cellular “societies”) over time, even if demands stay the same and average mtDNA quality doesn’t change. This is because the cell-to-cell variability in mtDNA quality does increase, carrying with it an added energetic challenge. This increased challenge could be a contributing factor to the collection of problems involved in ageing.

An overview of our approach. A mathematical model for the processes and "budget" involved in controlling mtDNA populations makes a general set of biological predictions and explains gene-therapy observations

Next, we found that cells with only low-quality mtDNA can perform worse than cells with a mix of low- and high-quality mtDNA. This is because low-quality mtDNA may consume less cellular resource, although global efficiency is decreased. Linked to this, removal of low-quality mtDNA (decommissioning bad power stations) alone is not always the best strategy to improve performance. Instead, jointly elevating low- and high-quality mtDNA levels, avoiding this detrimental mixed regime, is the best strategy for some situations. These insights may help explain some of the negative effects recently observed in cells with mixed mtDNA populations.

Our theory suggests that mixed mtDNA populations may do worse than pure ones, even if the pure population is a low-functionality mutant. Image from Hanne's post here 

We identified how best to control cellular mtDNA populations across the full range of possible populations, and used this insight to link with exciting gene therapies where low-quality mtDNA is preferentially removed through an experimental intervention (using so-called “endonucleases” to cut particular mtDNA molecules). We found that strong, single treatments will be outperformed by weaker, longer-term treatments, and identified how the mtDNA variability we know is present can practically effect the outcome of these therapies. We hope that the principles found in this work both add to our basic understanding of ageing and mixed (“heteroplasmic”) mitochondrial populations, and may inform more efficient therapeutic approaches in the future. Iain, Hanne, Nick

ARTICLE: Tension and Resolution

Tension and resolution: dynamic, evolving populations of organelle genomes within plant cells
IG Johnston
Molecular Plant 12 764 (2019)

Mitochondria and chloroplasts are compartments in cells that power complex life. Both started out billions of years ago as independent organisms with complete genomes, that were acquired by ancestral cells. Since these endosymbioses, the genomes of mitochondria (mt) and chloroplasts (cp) have become stripped down. Modern mt and cp have lost lots of genes either completely or the “host” cell nucleus. Mt and cp now exist in dynamic populations within the cells of modern organisms. In plants and algae, the two co-exist, sharing responsibility for the energy balance of the organism – and hence ultimately powering and feeding life, including the human population.

Plant mt and cp populations are weird. Different plants and algae have very different mt and cp genomes – some huge (many megabases, several chromosomes in the case of some mt) and some tiny. Unlike the more familiar animal (and human) case, plant mt genomes readily recombine, mixing up their structures and genetic content within the cell. Both mt and cp move around plant cells rapidly – we’re not sure why, particular for mt. Again, unlike animal mt, neither plant mt not cp are particularly prone to meet up and fuse into big networks – they usually stay as individual compartments, except for short interactions. We do know that if we perturb the physical or genetic dynamics of organelles, the plant suffers – which we can sometimes exploit in breeding efficient crops.

 Populations of mitochondria (A green, B) and chloroplasts (A blue, C) moving in the plant cell

In a recent review article here in Molecular Plant, we reviewed current knowledge about these dynamics and speculated about what principles these populations of mt and cp may be responding to. We first asked why mt and cp may retain different sets of genes in different species – a question we’ve touched upon before here (blog). Retaining more genes in organelles may have the “pro” of making individual organelles more independent, and better at responding to demands (see John Allen’s CoRR hypothesis, e.g. here). But there’s the “con” that organelles are dangerous places, and genes retained there may be more subject to damage than in the safe haven of the nucleus. So individual plants may choose to retain mt and cp genes for dynamism, or shift them to the nucleus for robustness. Neither extreme is perfect – there are always pros and cons – leading to a tension to which different plants have selected different resolutions.

Pursuing this line, we next speculated that because plants are immobile (and hence unable to move away from challenging conditions), they may favour the “dynamism” side over the “robustness” side. This would explain why they often retain more organelle genes than motile organisms, but would also predict that they face a double challenge: (i) more organelle genes and (ii) exposure to more challenging environments, both of which may lead to genetic damage. This could be a reason why plant organelles undergo recombination – as a way of ameliorating genetic damage. But again, there are pros and cons: the “pro” of fixing genetic damage is balanced by the “con” of recombination mixing and confusing genetic structure. Perhaps this is why the physical behaviour of plant organelles is different to that in animals – keeping mt and cp separate may limit the amount of recombination that can take place, allowing the plant to control this second pro-con tradeoff.

 

(left) the proposed tension between robustness (i) and dynamism (ii). Perhaps plants are more (ii)-like because they need to respond to fluctuating conditions... because of their immobility (right) with hypothesised knock-on consequences.

All of these ideas are presented as hypotheses, and we proposed some ways that a combination of new experiment and theory can help make progress understanding these complex, vital systems in future. Watch this space! Iain

ARTICLE: Coupling mitochondrial physics and genetics

Mitochondrial Network State Scales mtDNA Genetic Dynamics
Juvid Aryaman, Charlotte Bowles, Nick S. Jones and Iain G. Johnston
Genetics Early online July 10, 2019; https://doi.org/10.1534/genetics.119.302423

Mitochondrial DNA (mtDNA) populations within our cells encode vital energetic machinery. MtDNA is housed within mitochondria, cellular compartments lined by two membranes, that lead a very dynamic life. Individual mitochondria can fuse when they meet, and fused mitochondria can fragment to become individual smaller mitochondria, all the while moving throughout the cell. The reasons for this dynamic activity remain unclear (we’ve compared hypotheses about them before here and here, with blog articles here). But what influence do these physical mitochondrial dynamics have on the genetic composition of mtDNA populations?

MtDNA populations can, naturally or as a result of gene therapies, consist of a mixture of different mtDNA types. Typically, different cells will have different proportions of, say, type A and type B. For example, one cell may be 20% type A, another cell may be 40% type A, and a third may be 70% type A. This variability matters because when a certain threshold (often around 60%) is crossed for some mtDNA types, we get devastating diseases.

We previously showed mathematically (blog) and experimentally (blog) that this cell-to-cell variability in mtDNA proportions (often called “heteroplasmy variance” and sometimes referred to via the “mtDNA bottleneck”) is expected to increase linearly over time. However, this analysis pictured mtDNAs as individual molecules, outside of their mitochondrial compartments. When mitochondria fuse to form larger compartments, their mtDNA is more protected: smaller mitochondria (and their internal mtDNA) are subject to greater degradation. More degradation means more replication, and more opportunities for the fraction of a particular type of mtDNA to change per unit time. In a new paper here in Genetics, we show (using a mathematical tour de force by Juvid) that this protection can dramatically influence cell-to-cell mtDNA variability. Specifically, the rate of heteroplasmy variance increase is scaled by the proportion of mitochondria that exist in a fragmented state. (It turns out that it's the proportion of itochondria that are fragmented that's important -- not whether the rate of fission-fusion is fast or slow).

This has knock-on effects for how the cell can best get rid of low-quality mutant mtDNA. In particular, if mitochondria are allowed to fuse based on their quality (“selective fusion”), we show that intermediate rates of fusion are best for removing mutants. Too much fusion, and all mtDNA is protected; too little, and good mtDNA cannot be sorted from bad mtDNA using the mitochondrial network. This mechanism could help explain why we see different levels of mitochondrial fusion in different conditions. More broadly, this link between mitochondrial physics and genetics (which we’ve also speculated about here (blog) and here) suggests one way that selective pressures and tradeoffs could influence mitochondrial dynamics, giving rise to the wide variety of behaviours that remain unexplained. Juvid, Nick, and Iain

ARTICLE: How mitochondria can vary, and consequences for human health

(cross-posted from Imperial Mitochondriacs)

Mitochondria are components of the cell which are involved in generating “energy currency” molecules called ATP across much of complex life. Since many mitochondria exist within single cells (often hundreds or thousands), it is possible for the characteristics of individual mitochondria to vary within cells, and within tissues. This variation of mitochondrial characteristics can affect biological function and human health.

Since mitochondria possess their own, small, circular, DNA molecules (mtDNA), we can split mitochondrial characteristics into two categories: genetic and non-genetic. In our review, we discuss a number of aspects in which mitochondria vary, from both genetic and non-genetic perspectives. 

In terms of mitochondrial genetics, the amount of mtDNA per cell is variable. When a cell divides, its daughters receive a share of its parents mtDNA, but the split isn’t precisely 50/50, so cell division can cause variability in the number of mtDNAs per cell. As mtDNAs are replicated and degraded over time, errors in the copying process may give rise to mtDNA mutations, which may spread throughout a cell. Factors such as: the total amount, the rate of degradation/replication, the mean fraction of mutants, and the extent of fragmentation in the mitochondrial network, can all influence how variable the fraction of mutated mtDNAs becomes through time (see here for a preview of some upcoming work on this topic). The total amount, and mutated fraction of mtDNAs, are implicated in diseases such as neurodegeneration, as well as the ageing process.

Apart from genetic variations, there are many non-genetic features of mitochondria which also vary within and between cells. Changes in mtDNA sequence can change the amino-acid sequence of the proteins encoded by mtDNA, causing structural changes in the molecular machines which generate ATP. The shape of the membranes of mitochondria are also highly variable, and respond to mitochondrial activity through quantities such as pH, where mitochondrial activity itself may depend on mtDNA sequence. The previous two examples (mitochondrial protein and membrane structure) demonstrate how the genetic state of mitochondria may influence their non-genetic characteristics. Mitochondrial non-genetic characteristics may also influence the genetic state: for instance, mitochondrial membrane potential can influence the probability of a mitochondria being degraded, along with its mtDNA.

The inter-dependence of genetic and non-genetic characteristics demonstrate the complex feedback loops linking these two aspects of mitochondrial physiology. We suggest here that, since changes in mitochondrial genetics occur more slowly than most physical aspects of mitochondrial physiology, understanding mitochondrial genetics may be especially important in explaining phenomena such as ageing, which appears to be closely related to mitochondrial heterogeneity. You can freely access our work, which has recently been published in Frontiers in Genetics, as “Mitochondrial Heterogeneity” https://www.frontiersin.org/articles/10.3389/fgene.2018.00718/full Juvid, Iain and Nick.

 

ARTICLE: Time marches on -- mitochondria, ageing, and disease

Burgstaller, J.P., Kolbe, T., Havlicek, V., Hembach, S., Poulton, J., Piálek, J., Steinborn, R., Rülicke, T., Brem, G., Jones, N.S. and Johnston, I.G. Large-scale genetic analysis reveals mammalian mtDNA heteroplasmy dynamics and variance increase through lifetimes and generations. Nature communications, 9 2488 (2018)

DNA in mitochondria, the powerhouses of the cell, is passed down from mother to child. But there are many mitochondria in each cell, and these mitochondria may have different genetic features. If a mother carries a mixture of mitochondrial DNA (mtDNA) types, this can make it hard to say which features their children will inherit. For mothers carrying a disease-causing mtDNA mutation, this makes family planning and clinical therapies challenging.

In particular, the role of a mother's age has long been a mystery. Is the probability of a child inheriting a particular mtDNA feature higher when mothers are younger or older? An answer to this question could help plan clinical strategies to improve fertility and prevent the inheritance of deadly mitochondrial disease.

To address this, we worked with our excellent collaborators with a combination of maths, statistics, and experiment. Our collaborators used cutting-edge technology to reveal the mixtures of mtDNA in the egg cells of mother mice at a wide range of ages, and in the litters of offspring the mothers produced. This experimental work was the largest-scale study of mammalian mtDNA that we're aware of, involving thousands of observations throughout lifetimes and between generations. In concert, we developed a mathematical model describing the changes to, and inheritance of, mtDNA from mother to offspring. We combined the model and data to learn how different biological processes affect mtDNA through and between generations.

Cells contain populations of mitochondria, and these populations change over time. In European mice, we observed how variability in these populations evolves as mammals age and reproduce. We found that older mother have more varied mitochondria and pass this variance on to their offspring -- of central importance in the inheritance of genetic disease. 

We found that the variability of mtDNA dramatically increased as mothers aged. This means that the probability of inheriting more extreme -- both lower and higher -- levels of a genetic feature increases for older mothers. We also found that different mtDNA mixtures were inherited in different ways - with some mtDNA types favoured for inheritance and some disfavoured. We used our findings to create a way to predict how the risk that offspring would inherit disease-causing mtDNA features changes over time. Moving forward, we're aiming to harness these powerful ways of using large datasets to describe and predict the dynamics of mtDNA inheritance in humans, and to learn what it is about these mtDNA types that predicts their evolution across generations. You can read the article for free in Nature Communications here.

ARTICLE: How cells adapt to progressive increase in mitochondrial mutation

Aryaman, J., Johnston, I.G. and Jones, N.S. Mitochondrial DNA density homeostasis accounts for a threshold effect in a cybrid model of a human mitochondrial disease. Biochemical Journal, 474 4019 (2017).

Mitochondria produce the cell's major energy currency: ATP. If mitochondria become dysfunctional, this can be associated with a variety of devastating diseases, from Parkinson's disease to cancer. Technological advances have allowed us to generate huge volumes of data about these diseases. However, it can be a challenge to turn these large, complicated, datasets into basic understanding of how these diseases work, so that we can come up with rational treatments.

We were interested in a dataset (see here) which measured what happened to cells as their mitochondria became progressively more dysfunctional. A typical cell has roughly 1000 copies of mitochondrial DNA (mtDNA), which contains information on how to build some of the most important parts of the machinery responsible for making ATP in your cells. When mitochondrial DNA becomes mutated, these instructions accumulate errors, preventing the cell's energy machinery from working properly. Since your cells each contain about 1000 copies of mitochondrial DNA, it is interesting to think about what happens to a cell as the fraction of mutated mitochondrial DNA (called 'heteroplasmy') gradually increases.  We used maths to try and explain how a cell attempts to cope with increasing levels of heteroplasmy, resulting in a wealth of hypotheses which we hope to explore experimentally in the future.

The central idea arising from our analysis of this large dataset is that cells seem to attempt to maintain the number of normal mtDNAs per cell volume as heteroplasmy initially increases from 0% mutant. We suggest they do this by shrinking their size. By getting smaller, cells are able to reduce their energy demands as the fraction of mutant mtDNA increases, allowing them to balance their energy budget and maintain energy supply = demand. However, cells can only get so small and eventually the cell must change its strategy. At a critical fraction of mutated mtDNA (h* in the cartoon above), we suggest that cells switch on an alternative energy production mode called glycolysis. This causes energy supply to increase, and as a result, cells grow larger in size again. These ideas, as well as experimental proposals to test them, are freely available in the Biochemical Journal "Mitochondrial DNA Density Homeostasis Accounts for a Threshold Effect in a Cybrid Model of a Human Mitochondrial Disease". Juvid, Iain and Nick

ARTICLE: Supply, demand, energy, and death

Mitochondrial heterogeneity, metabolic scaling and cell death
J Aryaman, H Hoitzing, JP Burgstaller, IG Johnston, NS Jones
BioEssays e201700001; doi:10.1002/bies.201700001 (2017)

•  The links between mitochondrial functionality and various aspects of cell physiology remain unclear; we combine recent experimental insights with mathematical modelling to produce quantitative hypotheses linking metabolism, cell proliferation, and mitochondria.

Cells need energy to produce functional machinery, deal with challenges, and continue to grow and divide -- these activities and others are collectively referred to as "cell physiology". Mitochondria are the dominant energy sources in most of our cells, so we'd expect a strong link between how well mitochondria perform and cell physiology. Indeed, when mitochondrial energy production is compromised, deadly diseases can result -- as we've written about before.

The details of this link -- how cells with different mitochondrial populations may differ physiologically -- is not well understood. A recent article shed new light on this link by looking at a measure of mitochondrial functionality in cells of different sizes. They found what we'll call the "mitopeak" -- mitochondrial functionality peaks at intermediate cell sizes, with larger and smaller cells having less functional mitochondria. The subsequent interpretation was that there is an “optimal”, intermediate, size for cells. Above this size, it was suggested that a proposed universal relationship between the energy demands of organisms (from microorganisms to elephants) and their size predicts the reduction in the function of mitochondria. Smaller cells, which result from a large cell having divided, were suggested to have inherited their parent's low mitochondrial functionality. Cells were predicted to “reset” their mitochondrial activity as they initially grow and reach an “optimal” size.

We were interested in the mitopeak, and wondered if scientifically simpler hypotheses could account for it. Using mathematical modelling, our idea was to use the observation that as a cell becomes larger in volume, the size of its mitochondrial population (and hence power supply) increases in concert. We considered that a cell has power demands which also track its volume, as well as demands which are proportional to surface area and power demands which do not depend on cell size at all (such as the energetic cost of replicating the genome at cell division, since the size of a cell's genome does not depend on how big the cell is). Assuming that power supply = demand in a cell, then bigger cells may more easily satisfy e.g. the constant power demands. This is because the number of mitochondria increases with cell volume yet the constant demands remain the same regardless of cell size. In other words, if a cell has more mitochondria as it gets larger, then each mitochondrion has to work less hard to satisfy power demand.

To explain why the smallest cells also have mitochondria which do not appear to work hard, we suggested that some smaller cells could be in the process of dying. If smaller cells are more likely to die, and if dying cells have low mitochondrial functionality (both of these ideas are biologically supported), then, by combining this with the power supply/demand picture above, the observed mitopeak naturally emerges from our mathematical model.

As an alternative model, we also suggested that the mitopeak could come entirely from a nonlinear relationship between cell size and cell death, with mitochondrial functionality as a passive indicator of how healthy a cell is. This indicates the existence of multiple hypotheses which could explain this new dataset.

A recent study has provided new data for the relationship between cell physiology and mitochondrial functionality. We have used mathematical modelling to suggest that a mixture of cellular power demand scaling, as well as cell death, could intuitively account for these new data. However, a nonlinear relationship between cell death and cell size could also account for these data, as well as a nonlinear relationship between mitochondrial functionality and cell size, as proposed by the original authors of the dataset. By integrating such a relationship between cell size and mitochondrial functionality into one of our existing models, we found that this “mitopeak” helps explain a wider set of cell physiological data. Using our model to highlight these competing hypotheses, we suggest future experiments to gather further support for these potential explanations.

Interestingly, we also found that the mitopeak could be an alternative to one aspect of a model we used some time ago to explain a different dataset, looking at the physiological influence of mitochondrial variability. Then, we modelled the activity of mitochondria as a quantity that is inherited identically by each daughter cell from its parent, plus some noise -- noting that this was a guess at the true behaviour because we didn't have the data to make a firm statement. We needed this relationship because observed functionality varied comparatively little between sister cells but substantially across a population. The mitopeak induces this variability without needing random inheritance of functionality, and may thus be the refined picture we've been looking for. These ideas, and suggestions for future strategies to explore the link between mitochondria and cell physiology in more detail, are in our new BioEssays article here. Juvid, Nick, and Iain.

ARTICLE: Who keeps the plans for our power stations?

Evolutionary Inference across Eukaryotes Identifies Specific Pressures Favoring Mitochondrial Gene Retention

IG Johnston, BP Williams

Cell Systems 2 (2), 101-111 (2016)

• Why some genes are retained in mitochondria, where they are prone to disease-causing mutation, is a much-debated evolutionary question: we use new and generalisable maths and statistics to harness a large volume of sequence data and find the features of genes that predict the patterns of mitochondrial evolution that we observe.

Billions of years ago, a single-celled organism that would become our ancestor engulfed another smaller single-celled organism. The engulfed cell was probably intended to be lunch, but for reasons that remain mysterious (though recently explored here), it remained intact within our ancestor. It produced valuable chemicals that our ancestor could make use of, and was protected within the larger cell. This started a mutually beneficial relationship that evolved over billions of years to give rise to our situation today -- we are the descendants of the big cell, and our mitochondria are the descendants of the small, engulfed cell. 

As they were once independent organisms, mitochondria possess their own genomes (mitochondrial DNA, or mtDNA, which we've written about before). However, unlike the genomes of independent single-celled organisms like bacteria, mtDNA has only a handful of genes: why? Over evolutionary time, the majority of genes have either vanished from mtDNA or been transferred to the nucleus of the host cell. The reasons for transferring these genes to the nucleus are quite well understood; the nucleus is a safer environment for genes, less prone to mutation, and has several other evolutionary advantages.  But, given that transfer to the nucleus is possible, and genes in mtDNA are susceptible to mutation and damage (often giving rise to devastating diseases, which we study and try to prevent), why have mitochondria retained any genes at all? 

This question has been asked for decades, but until recently we lacked the data and the mathematical language to answer it quantitatively. Scientists energetically debate several different hypotheses: our approach attempts to let the data speak for itself without any preconceived ideas about which hypotheses are most likely. To this end, we built a mathematical model encompassing the evolutionary history of organisms with mitochondria, and a powerful statistical framework to amalgamate all the data that has been collected in recent years -- thousands of mitochondrial genomes from organisms from plants to protists (and humans) -- and harness it to compare the many disputed hypotheses addressing this question. 

Our mathematical approach allows us to "rewind the tape of evolution" and explore how mitochondrial genes have evolved. We're looking at Complex I -- an important protein complex involved in respiration -- over time, and watching the number of its subunits encoded in mitochondrial DNA (coloured black) decrease over evolutionary time, according to rules which we identify. The skyscrapers in the background are part of a graph describing how more mtDNA genes are lost over evolutionary history.

In a new paper in Cell Systems here (free here) we found several features that are most related to whether a gene is retained in mtDNA. Before discussing what they were, note that this picture -- several different features each with some influence -- explains and justifies the existing scientific debate. If hypothesis X and hypothesis Y both represent parts of the underlying "truth", then scientists advocating X alone and scientists advocating Y alone are neither completely wrong nor necessarily at odds -- everyone's partly right and the truth lies in the combination of the two arguments. 

The features that predict mtDNA gene retention are how central a gene's product is in its protein complex, the hydrophobicity of the protein the gene encodes, and the proportion of G's and C's in the gene's sequence. This suggests that genes are retained in mtDNA:  

(a) To allow local control of mitochondrial machinery (individual mitochondria can be controlled in response to cellular demands, rather than having to apply changes to the entire cellular population of mitochondria at once).  

(b) To prevent hydrophobic proteins ending up in the wrong place in the cell (if encoded by the far-away nucleus, these proteins may not be able to reach or enter the mitochondrion). 

(c and most speculatively) Because they are capable of withstanding the damaging environment of the mitochondria (GC-rich DNA and RNA is chemically more robust than GC-poor molecules). 

We found that the combination of the features we identified also predicted the success of experiments where scientists have attempted to mimic evolution and artificially transfer genes from the mitochondrion to the nucleus. Our results, as well as addressing a central mystery of evolutionary biology, thus also have the potential to inform synthetic biology approaches to tailor the genetics and bioenergetics of organisms. One final but important point is that the mathematical and statistical machinery we built for this project is highly generalisable and an efficient way of harnessing large sets of data about evolutionary and progressive processes -- we hope to use it to explore lots of other questions, including figuring out the pathways of disease progression and suggesting personalised medicine strategies in the clinic. Iain and Ben

ARTICLE: Go green -- recycle mitochondria

A novel quantitative assay of mitophagy: Combining high content fluorescence microscopy and mitochondrial DNA load to quantify mitophagy and identify novel pharmacological tools against pathogenic heteroplasmic mtDNA

• Mitophagy degrades mitochondria, and likely plays important roles in the cell's responses to mitochondrial disease, but is hard to measure and thus poorly understood: we propose new ways of measuring mitophagy and use them to explore drugs that may help change damaged mitochondrial populations

Mitochondria, as we've written about before, are important entities in our cells that produce energy and take part in many other vital processes. Mitochondrial DNA (mtDNA), inherited from our mothers, contains instructions on how to build important mitochondrial machinery. MtDNA is sometimes mutated, leading to problems with our mitochondria. How do our cells cope?

Mitophagy (from mito-(chondria) and -phagy (eating)) is a process by which cells degrade and recycle mitochondria, allowing dysfunctional mitochondria to be removed and replaced. Mitophagy is one of a number of cellular mechanisms that maintain a healthy population of mitochondria, and appears to play a central role in determining the inheritance and evolution of mtDNA over our lifetimes. However, our understanding of mitophagy is limited because it is hard to observe.

In a recent and epically-titled paper in Pharmacological Research here, we explore two different approaches for measuring mitophagy in cells. The first is physical. We used chemicals to make mitochondria glow red, and autophagosomes (the cellular machines responsible for the degradation of mitochondria) glow green. We then used a microscope to examine large numbers of cells and recorded how often red (mitochondria) and green (autophagosomes) were seen together, which we took to imply that mitophagy may be occurring. We confirmed that various drugs and chemicals known to affect mitophagy had the expected effects on this estimate of mitophagy, and that perturbing ATG7 (an essential part of the autophagic machinery) sustantially reduced our observed mitophagy levels.

We also subjected cells to stress by growing them with a less plentiful supply of energy. We found that this energy stress increased the amount of mitophagy (perhaps as cells struggle to make the very best of their mitochondrial populations). We also found that mitophagy broadly decreased in cells from older people, and was increased in cells from people carrying an mtDNA disease (negatively affecting mitochondrial functionality).

The second approach is genetic. In cells from patients with mtDNA disease, some mtDNA is normal and some is mutated -- we used genetic tools to measure the proportion of mutant mtDNA in cells. We observed that when we stressed patients' cells, levels of mutant mtDNA decreased while our physically observed measure of mitophagy increased, supporting a picture in which mitophagy removes dysfunctional mitochondria when energy output is of central importance. We also found evidence for undirected mitophagy, where mtDNA copy number is depleted with no preference for mutant or wildtype.

Observing the colocalisation of autophagosomes (green) and mitochondria (red), as well as the proportion of mutant mtDNA (white stars), allows a bilateral characterisation of mitophagy. The patterns of changes in these observations tell us about how drug treatments and different environments change mitochondrial populations.

The physical and genetic approaches give us two largely independent means to estimate mitophagy, placing our understanding of this vital process on a solid analytical foundation. We used these tools to assess the effects of various drugs on mitophagy, allowing us to characterise the effects of drugs like metformin (inhibiting mitophagy) and phenanthroline (inducing undirected mitophagy) in unprecedented detail and facilitating more precise statements about their utility in clinical contexts. Iain

ARTICLE: Warburg Ensemble

Monitoring Intracellular Oxygen Concentration: Implications for Hypoxia Studies and Real-time Oxygen Monitoring

• Cancer cells vary in how they produce their energy: we make progress understanding this variability, which may eventually help scientists design better therapies.

Cells can produce energy through several processes. We'll consider two – process "O" (for "oxidative phosphorylation"), and process "G" (for "glycolysis"). "O" uses oxygen, and harnesses the cell's mitochondria to produce energy. "G" does not use oxygen and produces energy without directly using mitochondria.

Healthy cells use both “O” and “G”, but cancer cells are often observed to rely on "G" much more. The shift away from "O+G" towards just "G" in cancer is often called the "Warburg effect", after Otto Warburg, who wrote about the shift in the 1950s. It remains unclear, however, whether the Warburg effect applies to all cancer cells under all conditions, or if different cells and different environments experience different shifts. This is important because understanding how cancer cells get their energy -- and, more generally, what changes occur in cancer cells compared to healthy cells -- may allow us to design therapies that challenge cancer cells while leaving healthy cells undamaged.

We used some fancy modern technology (focussed around the MitoXpress-Intra probe) to measure the difference between oxygen levels within a cell and oxygen levels in the cell's environment. We developed a mathematical way of producing "calibration curves", directly linking the observed MitoXpress behaviour to oxygen concentrations. If cells are using "G" alone, these levels are similar, as no oxygen is being consumed by the cells. If cells are also using "O", oxygen levels within cells should be rather lower than in their environment.

We found that two different cancer cell lines (with the rather jargon-y names "RD" and "U87MG") behaved surprisingly differently. When grown on glucose, U87MG looks quite "G", with oxygen levels within cells similar to those in the environment (e.g. 17.1% in cells, 18% outside). RD looks much more "O+G", with substantial differences between in-cell and outside-cell oxygen levels (e.g 13.2% in cells, 18% outside). Importantly, these findings were reproduced across a range of environmental oxygen levels (18% to 5%), modelling the range of conditions that cancer cells experience in tumours in the body. The two cancer cell lines thus seem to produce their energy in rather different ways, underlining that the Warburg effect is not an invariant across all cancers, and that treatments may be improved by taking this into account. We also showed that treating a different cancer cell line ("786-0") with phenformin, a drug inhibiting mitochondria, shifts cells away from "O+G" to "G", and that this shift can be monitored in real time with MitoXpress.

Different cancer cell lines (U87MG and RD) produce energy through different pathways, engaging more “G” (glycolysis) or “O” (oxidative phosphorylation). “O” uses oxygen (O2), lowering oxygen levels in cells compared to their environment. The different balance of “G” and “O” in different cases is important for understanding the heterogeneity of cancer.

Our paper appears in a book with the catchy title "Oxygen Transport to Tissue XXXVII", associated with the journal Advances in Experimental Medicine and Biology. You can get a sneak peek here and we'll update with a link when possible. Iain