Stochastic survival of the densest and mitochondrial DNA clonal expansion in aging.

The expansion of mitochondrial DNA molecules with deletions has been associated with aging, particularly in skeletal muscle fibers; its mechanism has remained unclear for three decades. Previous accounts have assigned a replicative advantage (RA) to mitochondrial DNA containing deletion mutations, but there is also evidence that cells can selectively remove defective mitochondrial DNA. Here we present a spatial model that, without an RA, but instead through a combination of enhanced density for mutants and noise, produces a wave of expanding mutations with speeds consistent with experimental data. A standard model based on RA yields waves that are too fast. We provide a formula that predicts that wave speed drops with copy number, consonant with experimental data. Crucially, our model yields traveling waves of mutants even if mutants are preferentially eliminated. Additionally, we predict that mutant loads observed in single-cell experiments can be produced by de novo mutation rates that are drastically lower than previously thought for neutral models. Given this exemplar of how spatial structure (multiple linked mtDNA populations), noise, and density affect muscle cell aging, we introduce the mechanism of stochastic survival of the densest (SSD), an alternative to RA, that may underpin other evolutionary phenomena.

Influencing dynamics on social networks without knowledge of network microstructure.

Social network-based information campaigns can be used for promoting beneficial health behaviours and mitigating polarization (e.g. regarding climate change or vaccines). Network-based intervention strategies typically rely on full knowledge of network structure. It is largely not possible or desirable to obtain population-level social network data due to availability and privacy issues. It is easier to obtain information about individuals' attributes (e.g. age, income), which are jointly informative of an individual's opinions and their social network position. We investigate strategies for influencing the system state in a statistical mechanics based model of opinion formation. Using synthetic and data-based examples we illustrate the advantages of implementing coarse-grained influence strategies on Ising models with modular structure in the presence of external fields. Our work provides a scalable methodology for influencing Ising systems on large graphs and the first exploration of the Ising influence problem in the presence of ambient (social) fields. By exploiting the observation that strong ambient fields can simplify control of networked dynamics, our findings open the possibility of efficiently computing and implementing public information campaigns using insights from social network theory without costly or invasive levels of data collection.

Discovering Cellular Mitochondrial Heteroplasmy Heterogeneity with Single Cell RNA and ATAC Sequencing.

Next-generation sequencing technologies have revolutionised the study of biological systems by enabling the examination of a broad range of tissues. Its application to single-cell genomics has generated a dynamic and evolving field with a vast amount of research highlighting heterogeneity in transcriptional, genetic and epigenomic state between cells. However, compared to these aspects of cellular heterogeneity, relatively little has been gleaned from single-cell datasets regarding cellular mitochondrial heterogeneity. Single-cell sequencing techniques can provide coverage of the mitochondrial genome which allows researchers to probe heteroplasmies at the level of the single cell, and observe interactions with cellular function. In this review, we give an overview of two popular single-cell modalities-single-cell RNA sequencing and single-cell ATAC sequencing-whose throughput and widespread usage offers researchers the chance to probe heteroplasmy combined with cell state in detailed resolution across thousands of cells. After summarising these technologies in the context of mitochondrial research, we give an overview of recent methods which have used these approaches for discovering mitochondrial heterogeneity. We conclude by highlighting current limitations of these approaches and open problems for future consideration.

Inference and influence of network structure using snapshot social behavior without network data.

Population behavior, like voting and vaccination, depends on the structure of social networks. This structure can differ depending on behavior type and is typically hidden. However, we do often have behavioral data, albeit only snapshots taken at one time point. We present a method jointly inferring a model for both network structure and human behavior using only snapshot population-level behavioral data. This exploits the simplicity of a few parameter model, geometric sociodemographic network model, and a spin-based model of behavior. We illustrate, for the European Union referendum and two London mayoral elections, how the model offers both prediction and the interpretation of the homophilic inclinations of the population. Beyond extracting behavior-specific network structure from behavioral datasets, our approach yields a framework linking inequalities and social preferences to behavioral outcomes. We illustrate potential network-sensitive policies: How changes to income inequality, social temperature, and homophilic preferences might have reduced polarization in a recent election.

Bladder pressure encoding by sacral dorsal root ganglion fibres: implications for decoding.

Objective.We aim at characterising the encoding of bladder pressure (intravesical pressure) by a population of sensory fibres. This research is motivated by the possibility to restore bladder function in elderly patients or after spinal cord injury using implanted devices, so called bioelectronic medicines. For these devices, nerve-based estimation of intravesical pressure can enable a personalized and on-demand stimulation paradigm, which has promise of being more effective and efficient. In this context, a better understanding of the encoding strategies employed by the body might in the future be exploited by informed decoding algorithms that enable a precise and robust bladder-pressure estimation.Approach.To this end, we apply information theory to microelectrode-array recordings from the cat sacral dorsal root ganglion while filling the bladder, conduct surrogate data studies to augment the data we have, and finally decode pressure in a simple informed approach.Main results.We find an encoding scheme by different main bladder neuron types that we divide into three response types (slow tonic, phasic, and derivative fibres). We show that an encoding by different bladder neuron types, each represented by multiple cells, offers reliability through within-type redundancy and high information rates through semi-independence of different types. Our subsequent decoding study shows a more robust decoding from mean responses of homogeneous cell pools.Significance.We have here, for the first time, established a link between an information theoretic analysis of the encoding of intravesical pressure by a population of sensory neurons to an informed decoding paradigm. We show that even a simple adapted decoder can exploit the redundancy in the population to be more robust against cell loss. This work thus paves the way towards principled encoding studies in the periphery and towards a new generation of informed peripheral nerve decoders for bioelectronic medicines.

Inference of a universal social scale and segregation measures using social connectivity kernels.

How people connect with one another is a fundamental question in the social sciences, and the resulting social networks can have a profound impact on our daily lives. Blau offered a powerful explanation: people connect with one another based on their positions in a social space. Yet a principled measure of social distance, allowing comparison within and between societies, remains elusive. We use the connectivity kernel of conditionally independent edge models to develop a family of segregation statistics with desirable properties: they offer an intuitive and universal characteristic scale on social space (facilitating comparison across datasets and societies), are applicable to multivariate and mixed node attributes, and capture segregation at the level of individuals, pairs of individuals and society as a whole. We show that the segregation statistics can induce a metric on Blau space (a space spanned by the attributes of the members of society) and provide maps of two societies. Under a Bayesian paradigm, we infer the parameters of the connectivity kernel from 11 ego-network datasets collected in four surveys in the UK and USA. The importance of different dimensions of Blau space is similar across time and location, suggesting a macroscopically stable social fabric. Physical separation and age differences have the most significant impact on segregation within friendship networks with implications for intergenerational mixing and isolation in later stages of life.

Cell identity and nucleo-mitochondrial genetic context modulate OXPHOS performance and determine somatic heteroplasmy dynamics.

Heteroplasmy, multiple variants of mitochondrial DNA (mtDNA) in the same cytoplasm, may be naturally generated by mutations but is counteracted by a genetic mtDNA bottleneck during oocyte development. Engineered heteroplasmic mice with nonpathological mtDNA variants reveal a nonrandom tissue-specific mtDNA segregation pattern, with few tissues that do not show segregation. The driving force for this dynamic complex pattern has remained unexplained for decades, challenging our understanding of this fundamental biological problem and hindering clinical planning for inherited diseases. Here, we demonstrate that the nonrandom mtDNA segregation is an intracellular process based on organelle selection. This cell type-specific decision arises jointly from the impact of mtDNA haplotypes on the oxidative phosphorylation (OXPHOS) system and the cell metabolic requirements and is strongly sensitive to the nuclear context and to environmental cues.

Characterizing soundscapes across diverse ecosystems using a universal acoustic feature set.

Natural habitats are being impacted by human pressures at an alarming rate. Monitoring these ecosystem-level changes often requires labor-intensive surveys that are unable to detect rapid or unanticipated environmental changes. Here we have developed a generalizable, data-driven solution to this challenge using eco-acoustic data. We exploited a convolutional neural network to embed soundscapes from a variety of ecosystems into a common acoustic space. In both supervised and unsupervised modes, this allowed us to accurately quantify variation in habitat quality across space and in biodiversity through time. On the scale of seconds, we learned a typical soundscape model that allowed automatic identification of anomalous sounds in playback experiments, providing a potential route for real-time automated detection of irregular environmental behavior including illegal logging and hunting. Our highly generalizable approach, and the common set of features, will enable scientists to unlock previously hidden insights from acoustic data and offers promise as a backbone technology for global collaborative autonomous ecosystem monitoring efforts.

A self-organizing, living library of time-series data.

Time-series data are measured across the sciences, from astronomy to biomedicine, but meaningful cross-disciplinary interactions are limited by the challenge of identifying fruitful connections. Here we introduce the web platform, CompEngine, a self-organizing, living library of time-series data, that lowers the barrier to forming meaningful interdisciplinary connections between time series. Using a canonical feature-based representation, CompEngine places all time series in a common feature space, regardless of their origin, allowing users to upload their data and immediately explore diverse data with similar properties, and be alerted when similar data is uploaded in future. In contrast to conventional databases which are organized by assigned metadata, CompEngine incentivizes data sharing by automatically connecting experimental and theoretical scientists across disciplines based on the empirical structure of the data they measure. CompEngine's growing library of interdisciplinary time-series data also enables the comprehensive characterization of time-series analysis algorithms across diverse types of empirical data.

A mechanistic explanation of the transition to simple multicellularity in fungi.

Development of multicellularity was one of the major transitions in evolution and occurred independently multiple times in algae, plants, animals, and fungi. However recent comparative genome analyses suggest that fungi followed a different route to other eukaryotic lineages. To understand the driving forces behind the transition from unicellular fungi to hyphal forms of growth, we develop a comparative model of osmotrophic resource acquisition. This predicts that whenever the local resource is immobile, hard-to-digest, and nutrient poor, hyphal osmotrophs outcompete motile or autolytic unicellular osmotrophs. This hyphal advantage arises because transporting nutrients via a contiguous cytoplasm enables continued exploitation of remaining resources after local depletion of essential nutrients, and more efficient use of costly exoenzymes. The model provides a mechanistic explanation for the origins of multicellular hyphal organisms, and explains why fungi, rather than unicellular bacteria, evolved to dominate decay of recalcitrant, nutrient poor substrates such as leaf litter or wood.

Community detection in networks without observing edges.

We develop a Bayesian hierarchical model to identify communities of time series. Fitting the model provides an end-to-end community detection algorithm that does not extract information as a sequence of point estimates but propagates uncertainties from the raw data to the community labels. Our approach naturally supports multiscale community detection and the selection of an optimal scale using model comparison. We study the properties of the algorithm using synthetic data and apply it to daily returns of constituents of the S&P100 index and climate data from U.S. cities.

Visualizing, quantifying, and manipulating mitochondrial DNA in vivo.

Mitochondrial DNA (mtDNA) encodes proteins and RNAs that support the functions of mitochondria and thereby numerous physiological processes. Mutations of mtDNA can cause mitochondrial diseases and are implicated in aging. The mtDNA within cells is organized into nucleoids within the mitochondrial matrix, but how mtDNA nucleoids are formed and regulated within cells remains incompletely resolved. Visualization of mtDNA within cells is a powerful means by which mechanistic insight can be gained. Manipulation of the amount and sequence of mtDNA within cells is important experimentally and for developing therapeutic interventions to treat mitochondrial disease. This review details recent developments and opportunities for improvements in the experimental tools and techniques that can be used to visualize, quantify, and manipulate the properties of mtDNA within cells.

The homeostatic dynamics of feeding behaviour identify novel mechanisms of anorectic agents.

Better understanding of feeding behaviour will be vital in reducing obesity and metabolic syndrome, but we lack a standard model that captures the complexity of feeding behaviour. We construct an accurate stochastic model of rodent feeding at the bout level in order to perform quantitative behavioural analysis. Analysing the different effects on feeding behaviour of peptide YY3-36 (PYY3-36), lithium chloride, glucagon-like peptide 1 (GLP-1), and leptin shows the precise behavioural changes caused by each anorectic agent. Our analysis demonstrates that the changes in feeding behaviour evoked by the anorectic agents investigated do not mimic the behaviour of well-fed animals and that the intermeal interval is influenced by fullness. We show how robust homeostatic control of feeding thwarts attempts to reduce food intake and how this might be overcome. In silico experiments suggest that introducing a minimum intermeal interval or modulating upper gut emptying can be as effective as anorectic drug administration.

Regulation of Mother-to-Offspring Transmission of mtDNA Heteroplasmy.

mtDNA is present in multiple copies in each cell derived from the expansions of those in the oocyte. Heteroplasmy, more than one mtDNA variant, may be generated by mutagenesis, paternal mtDNA leakage, and novel medical technologies aiming to prevent inheritance of mtDNA-linked diseases. Heteroplasmy phenotypic impact remains poorly understood. Mouse studies led to contradictory models of random drift or haplotype selection for mother-to-offspring transmission of mtDNA heteroplasmy. Here, we show that mtDNA heteroplasmy affects embryo metabolism, cell fitness, and induced pluripotent stem cell (iPSC) generation. Thus, genetic and pharmacological interventions affecting oxidative phosphorylation (OXPHOS) modify competition among mtDNA haplotypes during oocyte development and/or at early embryonic stages. We show that heteroplasmy behavior can fall on a spectrum from random drift to strong selection, depending on mito-nuclear interactions and metabolic factors. Understanding heteroplasmy dynamics and its mechanisms provide novel knowledge of a fundamental biological process and enhance our ability to mitigate risks in clinical applications affecting mtDNA transmission.

Precision identification of high-risk phenotypes and progression pathways in severe malaria without requiring longitudinal data.

More than 400,000 deaths from severe malaria (SM) are reported every year, mainly in African children. The diversity of clinical presentations associated with SM indicates important differences in disease pathogenesis that require specific treatment, and this clinical heterogeneity of SM remains poorly understood. Here, we apply tools from machine learning and model-based inference to harness large-scale data and dissect the heterogeneity in patterns of clinical features associated with SM in 2904 Gambian children admitted to hospital with malaria. This quantitative analysis reveals features predicting the severity of individual patient outcomes, and the dynamic pathways of SM progression, notably inferred without requiring longitudinal observations. Bayesian inference of these pathways allows us assign quantitative mortality risks to individual patients. By independently surveying expert practitioners, we show that this data-driven approach agrees with and expands the current state of knowledge on malaria progression, while simultaneously providing a data-supported framework for predicting clinical risk.

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

The dynamics of the cellular proportion of mutant mtDNA molecules is crucial for mitochondrial diseases. Cellular populations of mitochondria are under homeostatic control, but the details of the control mechanisms involved remain elusive. Here, we use stochastic modelling to derive general results for the impact of cellular control on mtDNA populations, the cost to the cell of different mtDNA states, and the optimisation of therapeutic control of mtDNA populations. This formalism yields a wealth of biological results, including that an increasing mtDNA variance can increase the energetic cost of maintaining a tissue, that intermediate levels of heteroplasmy can be more detrimental than homoplasmy even for a dysfunctional mutant, that heteroplasmy distribution (not mean alone) is crucial for the success of gene therapies, and that long-term rather than short intense gene therapies are more likely to beneficially impact mtDNA populations.

Mitochondrial Network State Scales mtDNA Genetic Dynamics

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

Correction to: PyPNS: Multiscale Simulation of a Peripheral Nerve in Python.

The original version of this article unfortunately contained a mistake. The following text: "This project has received funding from European Research Council (ERC) Synergy Grant no. 319818." is missing in the Acknowledgments.

Frequency and signature of somatic variants in 1461 human brain exomes.

PURPOSE: To systematically study somatic variants arising during development in the human brain across a spectrum of neurodegenerative disorders. METHODS: In this study we developed a pipeline to identify somatic variants from exome sequencing data in 1461 diseased and control human brains. Eighty-eight percent of the DNA samples were extracted from the cerebellum. Identified somatic variants were validated by targeted amplicon sequencing and/or PyroMark® Q24. RESULTS: We observed somatic coding variants present in >10% of sampled cells in at least 1% of brains. The mutational signature of the detected variants showed a predominance of C>T variants most consistent with arising from DNA mismatch repair, occurred frequently in genes that are highly expressed within the central nervous system, and with a minimum somatic mutation rate of 4.25 × 10-10 per base pair per individual. CONCLUSION: These findings provide proof-of-principle that deleterious somatic variants can affect sizeable brain regions in at least 1% of the population, and thus have the potential to contribute to the pathogenesis of common neurodegenerative diseases.

PyPNS: Multiscale Simulation of a Peripheral Nerve in Python.

Bioelectronic Medicines that modulate the activity patterns on peripheral nerves have promise as a new way of treating diverse medical conditions from epilepsy to rheumatism. Progress in the field builds upon time consuming and expensive experiments in living organisms. To reduce experimentation load and allow for a faster, more detailed analysis of peripheral nerve stimulation and recording, computational models incorporating experimental insights will be of great help. We present a peripheral nerve simulator that combines biophysical axon models and numerically solved and idealised extracellular space models in one environment. We modelled the extracellular space as a three-dimensional resistive continuum governed by the electro-quasistatic approximation of the Maxwell equations. Potential distributions were precomputed in finite element models for different media (homogeneous, nerve in saline, nerve in cuff) and imported into our simulator. Axons, on the other hand, were modelled more abstractly as one-dimensional chains of compartments. Unmyelinated fibres were based on the Hodgkin-Huxley model; for myelinated fibres, we adapted the model proposed by McIntyre et al. in 2002 to smaller diameters. To obtain realistic axon shapes, an iterative algorithm positioned fibres along the nerve with a variable tortuosity fit to imaged trajectories. We validated our model with data from the stimulated rat vagus nerve. Simulation results predicted that tortuosity alters recorded signal shapes and increases stimulation thresholds. The model we developed can easily be adapted to different nerves, and may be of use for Bioelectronic Medicine research in the future.