Source author record

Johannes Berg

Johannes Berg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.stat-mech cond-mat.dis-nn Molecular Networks Quantitative Methods Populations and Evolution

Catalog footprint

What is connected

11works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Stochastic clonal dynamics and genetic turnover in exponentially growing populations

We consider an exponentially growing population of cells undergoing mutations and ask about the effect of reproductive fluctuations (genetic drift) on its long-term evolution. We combine first step analysis with the stochastic dynamics of a birth-death process to analytically calculate the probability that the parent of a given genotype will go extinct. We compare the results with numerical simulations and show how this turnover of genetic clones can be used to infer the rates underlying the population dynamics. Our work is motivated by growing populations of tumour cells, the epidemic spread of viruses, and bacterial growth.

preprint2016arXiv

Network inference in the non-equilibrium steady state

Non-equilibrium systems lack an explicit characterisation of their steady state like the Boltzmann distribution for equilibrium systems. This has drastic consequences for the inference of parameters of a model when its dynamics lacks detailed balance. Such non-equilibrium systems occur naturally in applications like neural networks or gene regulatory networks. Here, we focus on the paradigmatic asymmetric Ising model and show that we can learn its parameters from independent samples of the non-equilibrium steady state. We present both an exact inference algorithm and a computationally more efficient, approximate algorithm for weak interactions based on a systematic expansion around mean-field theory. Obtaining expressions for magnetisations, two- and three-point spin correlations, we establish that these observables are sufficient to infer the model parameters. Further, we discuss the symmetries characterising the different orders of the expansion around the mean field and show how different types of dynamics can be distinguished on the basis of samples from the non-equilibrium steady state.

preprint2015arXiv

Pervasive adaptation of gene expression in Drosophila

Gene expression levels are important molecular quantitative traits that link genotypes to molecular functions and fitness. In Drosophila, population-genetic studies in recent years have revealed substantial adaptive evolution at the genomic level. However, the evolutionary modes of gene expression have remained controversial. Here we present evidence that adaptation dominates the evolution of gene expression levels in flies. We show that 63% of the observed expression divergence across seven Drosophila species are adaptive changes driven by directional selection. Our results are derived from the variation of expression within species and the time-resolved divergence across a family of related species, using a new inference method for selection. We identify functional classes of adaptively regulated genes, as well as sex-specific adaptation occurring predominantly in males. Our analysis opens a new avenue to map system-wide selection on molecular quantitative traits independently of their genetic basis.

preprint2014arXiv

What makes the lac-pathway switch: identifying the fluctuations that trigger phenotype switching in gene regulatory systems

Multistable gene regulatory systems sustain different levels of gene expression under identical external conditions. Such multistability is used to encode phenotypic states in processes including nutrient uptake and persistence in bacteria, fate selection in viral infection, cell cycle control, and development. Stochastic switching between different phenotypes can occur as the result of random fluctuations in molecular copy numbers of mRNA and proteins arising in transcription, translation, transport, and binding. However, which component of a pathway triggers such a transition is generally not known. By linking single-cell experiments on the lactose-uptake pathway in E. coli to molecular simulations, we devise a general method to pinpoint the particular fluctuation driving phenotype switching and apply this method to the transition between the uninduced and induced states of the lac genes. We find that the transition to the induced state is not caused only by the single event of lac-repressor unbinding, but depends crucially on the time period over which the repressor remains unbound from the lac-operon. We confirm this notion in strains with a high expression level of the repressor (leading to shorter periods over which the lac-operon remains unbound), which show a reduced switching rate. Our techniques apply to multi-stable gene regulatory systems in general and allow to identify the molecular mechanisms behind stochastic transitions in gene regulatory circuits.

preprint2013arXiv

Can we always sweep the details of RNA-processing under the carpet?

RNA molecules follow a succession of enzyme-mediated processing steps from transcription until maturation. The participating enzymes, for example the spliceosome for mRNAs and Drosha and Dicer for microRNAs, are also produced in the cell and their copy-numbers fluctuate over time. Enzyme copy-number changes affect the processing rate of the substrate molecules; high enzyme numbers increase the processing probability, low enzyme numbers decrease it. We study different RNA processing cascades where enzyme copy-numbers are either fixed or fluctuate. We find that for fixed enzyme-copy numbers the substrates at steady-state are Poisson-distributed, and the whole RNA cascade dynamics can be understood as a single birth-death process of the mature RNA product. In this case, solely fluctuations in the timing of RNA processing lead to variation in the number of RNA molecules. However, we show analytically and numerically that when enzyme copy-numbers fluctuate, the strength of RNA fluctuations increases linearly with the RNA transcription rate. This linear effect becomes stronger as the speed of enzyme dynamics decreases relative to the speed of RNA dynamics. Interestingly, we find that under certain conditions, the RNA cascade can reduce the strength of fluctuations in the expression level of the mature RNA product. Finally, by investigating the effects of processing polymorphisms we show that it is possible for the effects of transcriptional polymorphisms to be enhanced, reduced, or even reversed. Our results provide a framework to understand the dynamics of RNA processing.

preprint2012arXiv

Bethe-Peierls approximation and the inverse Ising model

We apply the Bethe-Peierls approximation to the problem of the inverse Ising model and show how the linear response relation leads to a simple method to reconstruct couplings and fields of the Ising model. This reconstruction is exact on tree graphs, yet its computational expense is comparable to other mean-field methods. We compare the performance of this method to the independent-pair, naive mean- field, Thouless-Anderson-Palmer approximations, the Sessak-Monasson expansion, and susceptibility propagation in the Cayley tree, SK-model and random graph with fixed connectivity. At low temperatures, Bethe reconstruction outperforms all these methods, while at high temperatures it is comparable to the best method available so far (Sessak-Monasson). The relationship between Bethe reconstruction and other mean- field methods is discussed.

preprint2012arXiv

Mean-field theory for the inverse Ising problem at low temperatures

The large amounts of data from molecular biology and neuroscience have lead to a renewed interest in the inverse Ising problem: how to reconstruct parameters of the Ising model (couplings between spins and external fields) from a number of spin configurations sampled from the Boltzmann measure. To invert the relationship between model parameters and observables (magnetisations and correlations) mean-field approximations are often used, allowing to determine model parameters from data. However, all known mean-field methods fail at low temperatures with the emergence of multiple thermodynamic states. Here we show how clustering spin configurations can approximate these thermodynamic states, and how mean-field methods applied to thermodynamic states allow an efficient reconstruction of Ising models also at low temperatures.

preprint2010arXiv

Significance analysis and statistical mechanics: an application to clustering

This paper addresses the statistical significance of structures in random data: Given a set of vectors and a measure of mutual similarity, how likely does a subset of these vectors form a cluster with enhanced similarity among its elements? The computation of this cluster p-value for randomly distributed vectors is mapped onto a well-defined problem of statistical mechanics. We solve this problem analytically, establishing a connection between the physics of quenched disorder and multiple testing statistics in clustering and related problems. In an application to gene expression data, we find a remarkable link between the statistical significance of a cluster and the functional relationships between its genes.

preprint2009arXiv

Adaptive gene regulatory networks

Regulatory interactions between genes show a large amount of cross-species variability, even when the underlying functions are conserved: There are many ways to achieve the same function. Here we investigate the ability of regulatory networks to reproduce given expression levels within a simple model of gene regulation. We find an exponentially large space of regulatory networks compatible with a given set of expression levels, giving rise to an extensive entropy of networks. Typical realisations of regulatory networks are found to share a bias towards symmetric interactions, in line with empirical evidence.

preprint2001arXiv

On random graphs and the statistical mechanics of granular matter

The dynamics of spins on a random graph with ferromagnetic three-spin interactions is used to model the compaction of granular matter under a series of taps. Taps are modelled as the random flipping of a small fraction of the spins followed by a quench at zero temperature. We find that the density approached during a logarithmically slow compaction - the random-close-packing density - corresponds to a dynamical phase transition. We discuss the the role of cascades of successive spin-flips in this model and link them with density-noise power fluctuations observed in recent experiments.