  A quantum model on the chemically and physically induced pluripotency in stem
cells is proposed. Based on the conformational Hamiltonian and the idea of slow
variables (molecular torsions) slaving fast ones the conversion from the
differentiate state to pluripotent state is defined as the quantum transition
between conformational states. The transitional rate is calculated and an
analytical form for the rate formulas is deduced. Then the dependence of the
rate on the number of torsion angles of the gene and the magnitude of the rate
can be estimated by comparison with protein folding. The reaction equations of
the conformational change of the pluripotency genes in chemical reprogramming
are given. The characteristic time of the chemical reprogramming is calculated
and the result is consistent with experiments. The dependence of the transition
rate on physical factors such as temperature, PH value and the volume and shape
of the coherent domain is analyzed from the rate equation. It is suggested that
by decreasing the coherence degree of some pluripotency genes a more effective
approach to the physically induced pluripotency can be made.


  The multiple sequence alignment (MSA) of a protein family provides a wealth
of information in terms of the conservation pattern of amino acid residues not
only at each alignment site but also between distant sites. In order to
statistically model the MSA incorporating both short-range and long-range
correlations as well as insertions, I have derived a lattice gas model of the
MSA based on the principle of maximum entropy. The partition function, obtained
by the transfer matrix method with a mean-field approximation, accounts for all
possible alignments with all possible sequences. The model parameters for
short-range and long-range interactions were determined by a self-consistent
condition and by a Gaussian approximation, respectively. Using this model with
and without long-range interactions, I analyzed the globin and V-set domains by
increasing the "temperature" and by "mutating" a site. The correlations between
residue conservation and various measures of the system's stability indicate
that the long-range interactions make the conservation pattern more specific to
the structure, and increasingly stabilize better conserved residues.


  Classical age-structured mass-action models such as the McKendrick-von
Foerster equation have been extensively studied but they are structurally
unable to describe stochastic fluctuations or population-size-dependent birth
and death rates. Stochastic theories that treat semi-Markov age-dependent
processes using e.g., the Bellman-Harris equation, do not resolve a
population's age-structure and are unable to quantify population-size
dependencies. Conversely, current theories that include size-dependent
population dynamics (e.g., mathematical models that include carrying capacity
such as the Logistic equation) cannot be easily extended to take into account
age-dependent birth and death rates. In this paper, we present a systematic
derivation of a new fully stochastic kinetic theory for interacting
age-structured populations. By defining multiparticle probability density
functions, we derive a hierarchy of kinetic equations for the stochastic
evolution of an ageing population undergoing birth and death. We show that the
fully stochastic age-dependent birth-death process precludes factorization of
the corresponding probability densities, which then must be solved by using a
BBGKY-like hierarchy. However, explicit solutions are derived in two simple
limits and compared with their corresponding mean-field results. Our results
generalize both deterministic models and existing master equation approaches by
providing an intuitive and efficient way to simultaneously model age- and
population-dependent stochastic dynamics applicable to the study of demography,
stem cell dynamics, and disease evolution.


  Biology achieves novel functions like error correction, ultra-sensitivity and
accurate concentration measurement at the expense of free energy through
Maxwell Demon-like mechanisms. The design principles and free energy trade-offs
have been studied for a variety of such mechanisms. In this review, we
emphasize a perspective based on dynamical phases that can explain
commonalities shared by these mechanisms. Dynamical phases are defined by
typical trajectories executed by non-equilibrium systems in the space of
internal states. We find that coexistence of dynamical phases can have dramatic
consequences for function vs free energy cost trade-offs. Dynamical phases can
also provide an intuitive picture of the design principles behind such
biological Maxwell Demons.


  The study of correlated mutations in alignments of homologous proteins proved
to be succesful not only in the prediction of their native conformation, but
also in the developement of a two-body effective potential between pairs of
amino acids. In the present work we extend the effective potential, introducing
a many--body term based on the same theoretical framework, making use of a
principle of maximum entropy. The extended potential performs better than the
two--body one in predicting the energetic effect of 308 mutations in 14
proteins (including membrane proteins). The average value of the parameters of
the many-body term correlates with the degree of hydrophobicity of the
corresponding residues, suggesting that this term partly reflects the effect of
the solvent.


  Some genes can change their relative locations in a genome. Thus for
different individuals of the same species, the orders of genes might be
different. Such jumping genes are called transposons. A practical problem is to
determine transposons in given gene sequences. Through an intuitive rule, we
transform the biological problem of determining transposons into a rigorous
mathematical problem of determining the longest common subsequence. Depending
on whether the gene sequence is linear (each sequence has a fixed head and
tail) or circular (we can choose any gene as the head, and the previous one is
the tail), and whether genes have multiple copies, we classify the problem of
determining transposons into four scenarios: (1) linear sequences without
duplicated genes; (2) circular sequences without duplicated genes; (3) linear
sequences with duplicated genes; (4) circular sequences with duplicated genes.
With the help of graph theory, we design fast algorithms for different
scenarios. We also derive some results that might be of theoretical interests
in combinatorics.


  We study the behaviour of double-stranded RNA under twist and tension using
oxRNA, a recently developed coarse-grained model of RNA. Introducing explicit
salt-dependence into the model allows us to directly compare our results to
data from recent single-molecule experiments. The model reproduces extension
curves as a function of twist and stretching force, including the buckling
transition and the behaviour of plectoneme structures. For negative
supercoiling, we predict denaturation bubble formation in plectoneme end-loops,
suggesting preferential plectoneme localisation in weak base sequences. OxRNA
exhibits a positive twist-stretch coupling constant, in agreement with recent
experimental observations.


  The spread of invasive species can have far reaching environmental and
ecological consequences. Understanding invasion spread patterns and the
underlying process driving invasions are key to predicting and managing
invasions. We combine a set of statistical methods in a novel way to
characterize local spread properties and demonstrate their application using
simulated and historical data on invasive insects. Our method uses a Gaussian
process fit to the surface of waiting times to invasion in order to
characterize the vector field of spread. Using this method we estimate with
statistical uncertainties the speed and direction of spread at each location.
Simulations from a stratified diffusion model verify the accuracy of our
method. We show how we may link local rates of spread to environmental
covariates for two case studies: the spread of the gypsy moth (Lymantria
dispar), and hemlock wolly adelgid (Adelges tsugae) in North America. We
provide an R-package that automates the calculations for any spatially
referenced waiting time data.


  Constraints placed upon the phenotypes of organisms result from their
interactions with the environment. Over evolutionary timescales, these
constraints feed back onto smaller molecular subnetworks comprising the
organism. The evolution of biological networks is studied by considering a
network of a few nodes embedded in a larger context. Taking into account this
fact that any network under study is actually embedded in a larger context, we
define network architecture, not on the basis of physical interactions alone,
but rather as a specification of the manner in which constraints are placed
upon the states of its nodes. We show that such network architectures
possessing cycles in their topology, in contrast to those that do not, may be
subjected to unsatisfiable constraints. This may be a significant factor
leading to selection biased against those network architectures where such
inconsistent constraints are more likely to arise. We proceed to quantify the
likelihood of inconsistency arising as a function of network architecture
finding that, in the absence of sampling bias over the space of possible
constraints and for a given network size, networks with a larger number of
cycles are more likely to have unsatisfiable constraints placed upon them. Our
results identify a constraint that, at least in isolation, would contribute to
a bias in the evolutionary process toward more hierarchical-modular versus
completely connected network architectures. Together, these results highlight
the context-dependence of the functionality of biological networks.


  Background: Highly effective direct-acting antiviral (DAA) regimens (90%
efficacy) are becoming available for hepatitis C virus (HCV) treatment. This
therapeutic revolution leads us to consider possibility of eradicating the
virus. However, for this, an effective cascade of care is required. Methods: In
the context of the incoming DAAs, we used a dynamic individual-based model
including a model of the people who inject drugs (PWID) social network to
simulate the impact of improved testing, linkage to care, and adherence to
treatment, and of modified treatment recommendation on the transmission and on
the morbidity of HCV in PWID in France. Results: Under the current incidence
and cascade of care, with treatment initiated at fibrosis stage $\ge$F2, the
HCV prevalence decreased from 42.8% to 24.9% [95% confidence interval
24.8%--24.9%] after 10 years. Changing treatment initiation criteria to treat
from F0 was the only intervention leading to a substantial additional decrease
in the prevalence, which fell to 11.6% [11.6%--11.7%] at 10 years. Combining
this change with improved testing, linkage to care, and adherence to treatment
decreased HCV prevalence to 7% [7%--7.1%] at 10 years and avoided 15.3%
[14.0%-16.6%] and 29.0% [27.9%--30.1%] of cirrhosis complications over 10 and
40 years respectively. Conclusion: A high decrease in viral transmission occurs
only when treatment is initiated before liver disease progresses to severe
stages, suggesting that systematic treatment in PWID, where incidence remains
high, would be beneficial. However, eradication will be difficult to achieve.


  Identifying the physical basis of heterosis (or hybrid vigor) has remained
elusive despite over a hundred years of research on the subject. The three main
theories of heterosis are dominance theory, overdominance theory, and epistasis
theory. Kacser and Burns (1981) identified the molecular basis of dominance,
which has greatly enhanced our understanding of its importance to heterosis.
This paper aims to explain how overdominance, and some features of epistasis,
can similarly emerge from the molecular dynamics of proteins. Possessing
multiple alleles at a gene locus results in the synthesis of different
allozymes at reduced concentrations. This in turn reduces the rate at which
each allozyme forms soluble oligomers, which are toxic and must be degraded,
because allozymes co-aggregate at low efficiencies. The model developed in this
paper will be used to explain how heterozygosity can impact the metabolic
efficiency of an organism. It can also explain why the viabilities of some
inbred lines seem to decline rapidly at high inbreeding coefficients (F > 0.5),
which may provide a physical basis for truncation selection for heterozygosity.
Finally, the model has implications for the ploidy level of organisms. It can
explain why polyploids are frequently found in environments where severe
physical stresses promote the formation of soluble oligomers. The model can
also explain why complex organisms, which need to synthesize aggregation-prone
proteins that contain intrinsically unstructured regions (IURs) and multiple
domains because they facilitate complex protein interaction networks (PINs),
tend to be diploid while haploidy tends to be restricted to relatively simple
organisms.


  On Monday January 19, 2015 a story broke that the National Football League
(NFL) had started an investigation into whether the New England Patriots
deliberately deflated the footballs they used during their championship win
over the Indianapolis Colts. Like an infectious disease, discussion regarding
Deflategate grew rapidly on social media sites in the hours and days after the
release of the story. However, after the Super Bowl was over, the scandal
slowly began to dissipate and lost much of the attention it had originally had,
as interest in the NFL wained at the completion of its season. We construct a
simple epidemic model for the infectiousness of the Deflategate news story. We
then use data from the social media site Twitter to estimate the parameters of
this model using standard techniques from the study of inverse problems. We
find that the infectiousness (as measured by the basic reproduction number) of
Deflategate rivals that of any infectious disease that we are aware of, and is
actually more infectious than recent news stories of greater importance - both
in terms of the basic reproduction number and in terms of the average amount of
time the average tweeter continued to tweet about the news story.


  We propose a novel molecular computing scheme for statistical inference. We
focus on the much-studied statistical inference problem of computing maximum
likelihood estimators for log-linear models. Our scheme takes log-linear models
to reaction systems, and the observed data to initial conditions, so that the
corresponding equilibrium of each reaction system encodes the corresponding
maximum likelihood estimator. The main idea is to exploit the coincidence
between thermodynamic entropy and statistical entropy. We map a Maximum Entropy
characterization of the maximum likelihood estimator onto a Maximum Entropy
characterization of the equilibrium concentrations for the reaction system.
This allows for an efficient encoding of the problem, and reveals that reaction
networks are superbly suited to statistical inference tasks. Such a scheme may
also provide a template to understanding how in vivo biochemical signaling
pathways integrate extensive information about their environment and history.


  Testing symptomatic individuals for a disease can deliver treatment
resources, if tests' results turn positive, which speeds up their treatment and
might also decrease individuals' contacts to other ones. An imperfect test,
however, might incorrectly consider susceptible individuals to be infected
(false positives). In this case, testing reduces the epidemic in the expense of
potentially misclassifying individuals. We present a mathematical model that
describes the dynamics of an infectious disease and its testing. Susceptible
individuals turn to "susceptible but deemed infected" at rate $\theta$.
Infected individuals go to a state "infected and tested positive" at rate
$\alpha$. Both of these rates are functions of test's sensitivity and
specificity. Analysis of the model permits us to derive an expression for $R_0$
and to find the conditions for reaching $R_0<1$, i.e., when the disease--free
equilibrium is stable. We find that under certain conditions it is possible to
get $R_0<1$, when originally, i.e., without testing, we would have $R_0>1$. We
also present numerical results to cover interesting scenarios such as using
different tests and to compare these results.


  Catastrophic transitions, where a system shifts abruptly between alternate
steady states, are a generic feature of many nonlinear systems. Recently these
regime shift were suggested as the mechanism underlies many ecological
catastrophes, such as desertification and coral reef collapses, which are
considered as a prominent threat to sustainability and to the well-being of
millions. Still, the methods proposed so far for the prediction of an imminent
transition are quite ineffective, and some empirical and theoretical studies
suggest that actual transitions may occur smoothly, without an abrupt shift.
Here we present a new diagnostic tool, based on monitoring the dynamics of
clusters through time. Our technique discriminates between systems with local
positive feedback, where the transition is abrupt, and systems with negative
density dependence, where the transition is smooth. Analyzing the spatial
dynamics of these two generic scenarios, we show that changes in the critical
cluster size provide a reliable early warning indicator for both transitions.
Our method may allow for the prediction, and thus hopefully the prevention of
such transitions, avoiding their destructive outcomes.


  For much of the last three decades Monte Carlo-simulation methods have been
the standard approach for accurately calculating the cyclization probability,
$J$, or J factor, for DNA models having sequence-dependent bends or
inhomogeneous bending flexibility. Within the last ten years, however,
approaches based on harmonic analysis of semi-flexible polymer models have been
introduced, which offer much greater computational efficiency than Monte Carlo
techniques. These methods consider the ensemble of molecular conformations in
terms of harmonic fluctuations about a well-defined elastic-energy minimum.
However, the harmonic approximation is only applicable for small systems,
because the accessible conformation space of larger systems is increasingly
dominated by anharmonic contributions. In the case of computed values of the J
factor, deviations of the harmonic approximation from the exact value of $J$ as
a function of DNA length have not been characterized. Using a recent,
numerically exact method that accounts for both anharmonic and harmonic
contributions to $J$ for wormlike chains of arbitrary size, we report here the
apparent error that results from neglecting anharmonic behavior. For wormlike
chains having contour lengths less than four times the persistence length the
error in $J$ arising from the harmonic approximation is generally small,
amounting to free energies less than the thermal energy, $k_B T$. For larger
systems, however, the deviations between harmonic and exact $J$ values increase
approximately linearly with size.


  Cooperation is a persistent behavioral pattern of entities pooling and
sharing resources. Its ubiquity in nature poses a conundrum. Whenever two
entities cooperate, one must willingly relinquish something of value to the
other. Why is this apparent altruism favored in evolution? Classical solutions
assume a net fitness gain in a cooperative transaction which, through
reciprocity or relatedness, finds its way back from recipient to donor. We seek
the source of this fitness gain. Our analysis rests on the insight that
evolutionary processes are typically multiplicative and noisy. Fluctuations
have a net negative effect on the long-time growth rate of resources but no
effect on the growth rate of their expectation value. This is an example of
non-ergodicity. By reducing the amplitude of fluctuations, pooling and sharing
increases the long-time growth rate for cooperating entities, meaning that
cooperators outgrow similar non-cooperators. We identify this increase in
growth rate as the net fitness gain, consistent with the concept of geometric
mean fitness in the biological literature. This constitutes a fundamental
mechanism for the evolution of cooperation. Its minimal assumptions make it a
candidate explanation of cooperation in settings too simple for other fitness
gains, such as emergent function and specialization, to be probable. One such
example is the transition from single cells to early multicellular life.


  Standard evolutionary theories of aging and mortality, implicitly based on
assumptions of spatial averaging, hold that natural selection cannot favor
shorter lifespan without direct compensating benefit to individual reproductive
success. Here we show that both theory and phenomenology are consistent with
programmed death. Spatial evolutionary models show that self-limited lifespan
robustly results in long-term benefit to a lineage; longer-lived variants may
have a reproductive advantage for many generations, but shorter lifespan
ultimately confers long-term reproductive advantage through environmental
feedback acting on much longer time scales. Numerous model variations produce
the same qualitative result, demonstrating insensitivity to detailed
assumptions; the key conditions under which self-limited lifespan is favored
are spatial extent and locally exhaustible resources. Numerous empirical
observations can parsimoniously be explained in terms of long-term selective
advantage for intrinsic mortality. Classically anomalous empirical data on
natural lifespans and intrinsic mortality, including observations of longer
lifespan associated with increased predation, and evidence of programmed death
in both unicellular and multicellular organisms, are consistent with specific
model predictions. The generic nature of the spatial model conditions under
which intrinsic mortality is favored suggests a firm theoretical basis for the
idea that evolution can quite generally select for shorter lifespan directly.


  The contributions of active deformation dynamics in cell nuclei to the
intra-nuclear positioning of hetero- and euchromatin are investigated. We
analyzed the behaviors of model chains containing two types of regions, one
with high and the other with low mobility, confined in a pulsating container.
Here, the regions with high and low mobility represent eu- and heterochromatic
regions, respectively, and the pulsating container simulates a nucleus
exhibiting dynamic deformations. The Brownian dynamics simulations of this
model show that the positioning of low mobility regions transition from sites
near the periphery to the central region of the container if the affinity
between low mobility regions and the container periphery disappears. Here, the
former and latter positioning are similar to the "conventional" and "inverted"
chromatin positioning observed in nuclei of normal differentiated cells and
cells lacking Lamin-related proteins like mouse rod photoreceptor cell.


  The opportunity to tell a white lie (i.e., a lie that benefits another
person) generates a moral conflict between two opposite moral dictates, one
pushing towards telling always the truth and the other pushing towards helping
others. Here we study how people resolve this moral conflict. What does telling
a white lie signal about a person's pro-social tendencies? To answer this
question, we conducted a two-stage 2x2 experiment. In the first stage, we used
a Deception Game to measure aversion to telling a Pareto white lie (i.e., a lie
that helps both the liar and the listener), and aversion to telling an
altruistic white lie (i.e., a lie that helps the listener at the expense of the
liar). In the second stage we measured altruistic tendencies using a Dictator
Game and cooperative tendencies using a Prisoner's dilemma. We found three
major results: (i) both altruism and cooperation are positively correlated with
aversion to telling a Pareto white lie; (ii) both altruism and cooperation are
negatively correlated with aversion to telling an altruistic white lie; (iii)
men are more likely than women to tell an altruistic white lie, but not to tell
a Pareto white lie. Our results shed light on the moral conflit between
pro-sociality and truth-telling. In particular, the first finding suggests that
a significant proportion of people have non-distributional notions of what the
right thing to do is: irrespective of their economic consequences, they tell
the truth, they cooperate, they share their money.


  The Kolmogorov-Arnold stochasticity parameter technique is applied for the
first time to the study of cancer genome sequencing, to reveal mutations. Using
data generated by next generation sequencing technologies, we have analyzed the
exome sequences of brain tumor patients with matched tumor and normal blood. We
show that mutations contained in sequencing data can be revealed using this
technique thus providing a new methodology for determining subsequences of
given length containing mutations i.e. its value differs from those of
subsequences without mutations. A potential application for this technique
involves simplifying the procedure of finding segments with mutations, speeding
up genomic research, and accelerating its implementation in clinical
diagnostic. Moreover, the prediction of a mutation associated to a family of
frequent mutations in numerous types of cancers based purely on the value of
the Kolmogorov function, indicates that this applied marker may recognize
genomic sequences that are in extremely low abundance and can be used in
revealing new types of mutations.


  A model of two consumer-resource systems linked by interspecific interference
competition of consumers is considered. The basic assumption of the model is
that the dynamics of the resource is much slower than that of the consumer. In
the absence of interaction each consumer-resource pair has a unique stable
steady state which is completely nonoscillatory. When weakly coupled, the
consumer-resource pairs are shown to exhibit sustained low-frequency
synchronous antiphase relaxation oscillations.


  The asymmetrical inheritance of plasmid DNA, as well as other cellular
components, has been shown to be involved in replicative aging. In
Saccharomyces cerevisiae, there is an ongoing debate regarding the mechanisms
underlying this important asymmetry. Currently proposed models suggest it is
established via diffusion, but differ on whether a diffusion barrier is
necessary or not. However, no study so far incorporated key aspects to
segregation, such as dynamic morphology changes throughout anaphase or plasmids
size. Here, we determine the distinct effects and contributions of individual
cellular variability, plasmid volume and moving boundaries in the asymmetric
segregation of plasmids. We do this by measuring cellular nuclear geometries
and plasmid diffusion rates with confocal microscopy, subsequently
incorporating this data into a growing domain stochastic spatial simulator. Our
modelling and simulations confirms that plasmid asymmetrical inheritance does
not require an active barrier to diffusion, and provides a full analysis on
plasmid size effects.


  We make use of ideas from the theory of complex networks to implement a
machine learning classification of human DNA methylation data, that carry
signatures of cancer development. The data were obtained from patients with
various kinds of cancers and represented as parenclictic networks, wherein
nodes correspond to genes, and edges are weighted according to pairwise
variation from control group subjects. We demonstrate that for the $10$ types
of cancer under study, it is possible to obtain a high performance of binary
classification between cancer-positive and negative samples based on network
measures. Remarkably, an accuracy as high as $93-99\%$ is achieved with only
$12$ network topology indices, in a dramatic reduction of complexity from the
original $15295$ gene methylation levels. Moreover, it was found that the
parenclictic networks are scale-free in cancer-negative subjects, and deviate
from the power-law node degree distribution in cancer. The node centrality
ranking and arising modular structure could provide insights into the systems
biology of cancer.


  Understanding the mechanisms by which cells coordinate their size with their
ability to divide has long attracted the interest of biologists. The Target of
Rapamycin (TOR) pathway is becoming increasingly recognized as a master
regulator of cell size, however less is known how TOR activity might be coupled
with the cell cycle. Here, we establish that mTOR complex 1 (mTORC1) promotes
cytokinesis through activation of a Rho GTPase-Rho Kinase (ROCK) signaling
cascade. Hyperactivation of mTORC1 signaling by depletion of any of its
negative regulators: TSC1, TSC2, PTEN, or DEPTOR, induces polyploidy in a
rapamycin-sensitive manner. mTORC1 hyperactivation-mediated polyploidization
occurs by a prolonged, but ultimately failed attempt at abcission followed by
re-fusion. Similar to the effects of ROCK2 overexpression, these mTORC1-driven
aberrant cytokinesis events are accompanied by increased Rho-GTP loading,
extensive plasma membrane blebbing, and increased actin-myosin contractility,
all of which can be rescued by either mTORC1 or ROCK inhibition. These results
provide evidence for the existence of a novel mTORC1-Rho-ROCK pathway during
cytokinesis and suggest that mTORC1 might play a critical role in setting the
size at which a mammalian cell divides.


  Phylogenetic networks are rooted, labelled directed acyclic graphs which are
commonly used to represent reticulate evolution. There is a close relationship
between phylogenetic networks and multi-labelled trees (MUL-trees). Indeed, any
phylogenetic network $N$ can be 'unfolded' to obtain a MUL-tree $U(N)$ and,
conversely, a MUL-tree $T$ can in certain circumstances be 'folded' to obtain a
phylogenetic network $F(T)$ that exhibits $T$. In this paper, we study
properties of the operations $U$ and $F$ in more detail. In particular, we
introduce the class of stable networks, phylogenetic networks $N$ for which
$F(U(N))$ is isomorphic to $N$, characterise such networks, and show that that
they are related to the well-known class of tree-sibling networks. We also
explore how the concept of displaying a tree in a network $N$ can be related to
displaying the tree in the MUL-tree $U(N)$. To do this, we develop a
phylogenetic analogue of graph fibrations. This allows us to view $U(N)$ as the
analogue of the universal cover of a digraph, and to establish a close
connection between displaying trees in $U(N)$ and reconciling phylogenetic
trees with networks.


  Cooperation is ubiquitous in biological and social systems. Previous studies
revealed that a preference toward similar appearance promotes cooperation, a
phenomenon called tag-mediated cooperation or communitarian cooperation. This
effect is enhanced when a spatial structure is incorporated, because space
allows agents sharing an identical tag to regroup to form locally cooperative
clusters. In spatially distributed settings, one can also consider migration of
organisms, which has a potential to further promote evolution of cooperation by
facilitating spatial clustering. However, it has not yet been considered in
spatial tag-mediated cooperation models. Here we show, using computer
simulations of a spatial model of evolutionary games with organismal migration,
that tag-based segregation and homophilic cooperation arise for a wide range of
parameters. In the meantime, our results also show another evolutionarily
stable outcome, where a high level of heterophilic cooperation is maintained in
spatially well-mixed patterns. We found that these two different forms of
tag-mediated cooperation appear alternately as the parameter for temptation to
defect is increased.


  We explain how hierarchical organization of biological systems emerges
naturally during evolution, through a transition in the units of individuality.
We will show how these transitions are the result of competing selective forces
operating at different levels of organization, each level having different
units of individuality. Such a transition represents a singular point in the
evolutionary process, which we will show corresponds to a phase transition in
the way information is encoded, with the formation of self-reinforcing
information bottlenecks. We present an abstract model for characterizing these
transitions that is quite general, applicable to many different versions of
such transitions. As a concrete example, we consider the transition to
multicellularity. Specifically, we study a stochastic model where isolated
communities of interacting individuals (e.g. cells) undergo a transition to
higher-order individuality (e.g. multicellularity). This transition is
indicated by the marked decrease in the number of cells utilized to generate
new communities from pre-existing ones. In this sense, the community begins to
reproduce as a whole via a decreasing number of cells. We show that the fitness
barrier to this transition is strongly reduced by horizontal gene transfer.
These features capture two of the most prominent aspects of the transition to
multicellularity: the evolution of a developmental process and reproduction
through a unicellular bottleneck.


  The total-evidence approach to divergence-time dating uses molecular and
morphological data from extant and fossil species to infer phylogenetic
relationships, species divergence times, and macroevolutionary parameters in a
single coherent framework. Current model-based implementations of this approach
lack an appropriate model for the tree describing the diversification and
fossilization process and can produce estimates that lead to erroneous
conclusions. We address this shortcoming by providing a total-evidence method
implemented in a Bayesian framework. This approach uses a mechanistic tree
prior to describe the underlying diversification process that generated the
tree of extant and fossil taxa. Previous attempts to apply the total-evidence
approach have used tree priors that do not account for the possibility that
fossil samples may be direct ancestors of other samples. The fossilized
birth-death (FBD) process explicitly models the diversification, fossilization,
and sampling processes and naturally allows for sampled ancestors. This model
was recently applied to estimate divergence times based on molecular data and
fossil occurrence dates. We incorporate the FBD model and a model of
morphological trait evolution into a Bayesian total-evidence approach to dating
species phylogenies. We apply this method to extant and fossil penguins and
show that the modern penguins radiated much more recently than has been
previously estimated, with the basal divergence in the crown clade occurring at
~12.7 Ma and most splits leading to extant species occurring in the last 2
million years. Our results demonstrate that including stem-fossil diversity can
greatly improve the estimates of the divergence times of crown taxa. The method
is available in BEAST2 (v. 2.4) www.beast2.org with packages SA (v. at least
1.1.4) and morph-models (v. at least 1.0.4).


  We study the dynamics of the nonlinear $q$-voter model with inflexible
zealots in a finite well-mixed population. In this system, each individual
supports one of two parties and is either a susceptible voter or an inflexible
zealot. At each time step, a susceptible adopts the opinion of a neighbor if
this belongs to a group of $q\geq 2$ neighbors all in the same state, whereas
inflexible zealots never change their opinion. In the presence of zealots of
both parties the model is characterized by a fluctuating stationary state and,
below a zealotry density threshold, the distribution of opinions is bimodal.
After a characteristic time, most susceptibles become supporters of the party
having more zealots and the opinion distribution is asymmetric. When the number
of zealots of both parties is the same, the opinion distribution is symmetric
and, in the long run, susceptibles endlessly swing from the state where they
all support one party to the opposite state. Above the zealotry density
threshold, when there is an unequal number of zealots of each type, the
probability distribution is single-peaked and non-Gaussian. These properties
are investigated analytically and with stochastic simulations. We also study
the mean time to reach a consensus when zealots support only one party.


  Based on the concept of DNA strand displacement and DNA strand algebra we
have developed a method for logical inference which is not based on silicon
based computing. Essentially, it is a paradigm shift from silicon to carbon. In
this paper we have considered the inference mechanism, viz. modus ponens, to
draw conclusion from any observed fact. Thus, the present approach to logical
inference based on DNA strand algebra is basically an attempt to develop expert
system design in the domain of DNA computing. We have illustrated our
methodology with respect to worked out example. Our methodology is very
flexible for implementation of different expert system applications.


  We present a model for the growth of the transportation network inside nests
of the social insect subfamily Termitinae (Isoptera, termitidae). These nests
consist of large chambers (nodes) connected by tunnels (edges). The model based
on the empirical analysis of the real nest networks combined with pruning (edge
removal, either random or weighted by betweenness centrality) and a memory
effect (preferential growth from the latest added chambers) successfully
predicts emergent nest properties (degree distribution, size of the largest
connected component, average path lengths, backbone link ratios, and local
graph redundancy). The two pruning alternatives can be associated with
different genuses in the subfamily. A sensitivity analysis on the pruning and
memory parameters indicates that Termitinae networks favor fast internal
transportation over efficient defense strategies against ant predators. Our
results provide an example of how complex network organization and efficient
network properties can be generated from simple building rules based on local
interactions and contribute to our understanding of the mechanisms that come
into play for the formation of termite networks and of biological
transportation networks in general.


  Evolution of gene regulation is crucial for our understanding of the
phenotypic differences between species, populations and individuals.
Sequence-specific binding of transcription factors to the regulatory regions on
the DNA is a key regulatory mechanism that determines gene expression and hence
heritable phenotypic variation. We use a biophysical model for directional
selection on gene expression to estimate the rates of gain and loss of
transcription factor binding sites (TFBS) in finite populations under both
point and insertion/deletion mutations. Our results show that these rates are
typically slow for a single TFBS in an isolated DNA region, unless the
selection is extremely strong. These rates decrease drastically with increasing
TFBS length or increasingly specific protein-DNA interactions, making the
evolution of sites longer than ~10 bp unlikely on typical eukaryotic speciation
timescales. Similarly, evolution converges to the stationary distribution of
binding sequences very slowly, making the equilibrium assumption questionable.
The availability of longer regulatory sequences in which multiple binding sites
can evolve simultaneously, the presence of "pre-sites" or partially decayed old
sites in the initial sequence, and biophysical cooperativity between
transcription factors, can all facilitate gain of TFBS and reconcile
theoretical calculations with timescales inferred from comparative genetics.


  The recent super-exponential growth in the amount of sequencing data
generated worldwide has put techniques for compressed storage into the focus.
Most available solutions, however, are strictly tied to specific bioinformatics
formats, sometimes inheriting from them suboptimal design choices; this hinders
flexible and effective data sharing. Here we present CARGO (Compressed
ARchiving for GenOmics), a high-level framework to automatically generate
software systems optimized for the compressed storage of arbitrary types of
large genomic data collections. Straightforward applications of our approach to
FASTQ and SAM archives require a few lines of code, produce solutions that
match and sometimes outperform specialized format-tailored compressors, and
scale well to multi-TB datasets.


  In this paper we propose network methodology to infer prognostic cancer
biomarkers based on the epigenetic pattern DNA methylation. Epigenetic
processes such as DNA methylation reflect environmental risk factors, and are
increasingly recognised for their fundamental role in diseases such as cancer.
DNA methylation is a gene-regulatory pattern, and hence provides a means by
which to assess genomic regulatory interactions. Network models are a natural
way to represent and analyse groups of such interactions. The utility of
network models also increases as the quantity of data and number of variables
increase, making them increasingly relevant to large-scale genomic studies. We
propose methodology to infer prognostic genomic networks from a DNA
methylation-based measure of genomic interaction and association. We then show
how to identify prognostic biomarkers from such networks, which we term
`network community oncomarkers'. We illustrate the power of our proposed
methodology in the context of a large publicly available breast cancer dataset.


  The prevalence of neutral mutations implies that biological systems typically
have many more genotypes than phenotypes. But can the way that genotypes are
distributed over phenotypes determine evolutionary outcomes? Answering such
questions is difficult because the number of genotypes can be
hyper-astronomically large. By solving the genotype-phenotype (GP) map for RNA
secondary structure for systems up to length $L=126$ nucleotides (where the set
of all possible RNA strands would weigh more than the mass of the visible
universe) we show that the GP map strongly constrains the evolution of
non-coding RNA (ncRNA). Simple random sampling over genotypes predicts the
distribution of properties such as the mutational robustness or the number of
stems per secondary structure found in naturally occurring ncRNA with
surprising accuracy. Since we ignore natural selection, this strikingly close
correspondence with the mapping suggests that structures allowing for
functionality are easily discovered, despite the enormous size of the genetic
spaces. The mapping is extremely biased: the majority of genotypes map to an
exponentially small portion of the morphospace of all biophysically possible
structures. Such strong constraints provide a non-adaptive explanation for the
convergent evolution of structures such as the hammerhead ribozyme. These
results presents a particularly clear example of bias in the arrival of
variation strongly shaping evolutionary outcomes and may be relevant to Mayr's
distinction between proximate and ultimate causes in evolutionary biology.


  Motivation: RNA-seq has made feasible the analysis of a whole set of
expressed mRNAs. Mapping-based assembly of RNA-seq reads sometimes is
infeasible due to lack of high-quality references. However, de novo assembly is
very challenging due to uneven expression levels among transcripts and also the
read coverage variation within a single transcript. Existing methods either
apply de Bruijn graphs of single-sized k-mers to assemble the full set of
transcripts, or conduct multiple runs of assembly, but still apply graphs of
single-sized k-mers at each run. However, a single k-mer size is not suitable
for all the regions of the transcripts with varied coverage. Contribution: This
paper presents a de novo assembler Bermuda with new insights for handling
uneven coverage. Opposed to existing methods that use a single k-mer size for
all the transcripts in each run of assembly, Bermuda self-adaptively uses a few
k-mer sizes to assemble different regions of a single transcript according to
their local coverage. As such, Bermuda can deal with uneven expression levels
and coverage not only among transcripts, but also within a single transcript.
Extensive tests show that Bermuda outperforms popular de novo assemblers in
reconstructing unevenly-expressed transcripts with longer length, better
contiguity and lower redundancy. Further, Bermuda is computationally efficient
with moderate memory consumption.


  We present an analytical treatment of a genetic switch model consisting of
two mutually inhibiting genes operating without cooperative binding of the
corresponding transcription factors. Previous studies have numerically shown
that these systems can exhibit bimodal dynamics without possessing two stable
fixed points at the deterministic level. We analytically show that bimodality
is induced by the noise and find the critical repression strength that controls
a transition between the bimodal and non-bimodal regimes. We also identify
characteristic polynomial scaling laws of the mean switching time between
bimodal states. These results, independent of the model under study, reveal
essential differences between these systems and systems with cooperative
binding, where there is no critical threshold for bimodality and the mean
switching time scales exponentially with the system size.


  A major issue in the modelling of weed resistance to herbicide lies in
effectively handling the wild dynamics, that is, the allele frequency prior to
the herbicide application, and in particular when starting to use the
herbicide. The wild allele frequency is a key variable as the resistance
evolution is highly sensitive to it, moreover it is extremely difficult to
measure in the agricultural field. In this paper we propose a model for weed
control that handles the allele frequency in a direct manner, grounded on the
very dynamics of the weed life cycle, with no need of a priori distributions
nor the use of the Hardy-Weinberg equilibrium. The proposed model is individual
based, stochastic, and considers some phenomena like the relative fitnesses and
mutation that are prominent in the resistance dynamics without herbicide. A
case study is presented for the herbicide nicosulfuron in a field with the weed
Bidens pilosa. Another two models having a standard deterministic dynamics are
compared with ours in terms of the initial allele frequency, its time
evolution, and the resistance visualization in the field, indicating that the
proposed model is effective to provide more realistic simulations for the weed
resistance.


  Oscillation is an important cellular process that regulates timing of
different vital life cycles. However, in the noisy cellular environment,
oscillations can be highly inaccurate due to phase fluctuations. It remains
poorly understood how biochemical circuits suppress phase fluctuations and what
is the incurred thermodynamic cost. Here, we study four different types of
biochemical oscillations representing three basic oscillation motifs shared by
all known oscillatory systems. We find that the phase diffusion constant
follows the same inverse dependence on the free energy dissipation per period
for all systems studied. This relationship between the phase diffusion and
energy dissipation is shown analytically in a model of noisy oscillation.
Microscopically, we find that the oscillation is driven by multiple
irreversible cycles that hydrolyze the fuel molecules such as ATP; the number
of phase coherent periods is proportional to the free energy consumed per
period. Experimental evidence in support of this universal relationship and
testable predictions are also presented.


  We wish to formally test for changes in the taxonomic diversity of a
community, especially in the presence of high latent diversity. Drawing on the
meta-analysis literature, we construct a model for diversity that accounts for
covariate effects as well as sampling variability. This permits inference for
changes in richness with covariates and also a test for homogeneity. We argue
that we can use the principles of shrinkage estimation to improve richness
estimation in this nonstandard context, which is especially important given the
high variance of richness estimators and the increasing abundance of community
composition data. We demonstrate the methodology under simulation, in a gut
microbiome study (testing for a decrease in richness with antibiotics), and in
a soil microbiome study (testing for homogeneity of replicates). We believe
that this is the first formal procedure for analyzing changes in species
richness.


  RNA-Seq and gene expression microarrays provide comprehensive profiles of
gene activity, but lack of reproducibility has hindered their application. A
key challenge in the data analysis is the normalization of gene expression
levels, which is currently performed following the implicit assumption that
most genes are not differentially expressed. Here, we present a mathematical
approach to normalization that makes no assumption of this sort. We have found
that variation in gene expression is much larger than currently believed, and
that it can be measured with available assays. Our results also explain, at
least partially, the reproducibility problems encountered in transcriptomics
studies. We expect that this improvement in detection will help efforts to
realize the full potential of gene expression profiling, especially in analyses
of cellular processes involving complex modulations of gene expression.


  Nuclear magnetic resonance (NMR) spectroscopy is a powerful method for the
investigation of three-dimensional structures of biological molecules such as
proteins. Determining a protein structure is essential for understanding its
function and alterations in function which lead to disease. One of the major
challenges of the post-genomic era is to obtain structural and functional
information on the many unknown proteins encoded by thousands of newly
identified genes. The goal of this research is to design an algorithm capable
of automating the analysis of backbone protein NMR data by implementing AI
strategies such as greedy and A* search.


  Comparative analyses of protein-protein interaction networks play important
roles in the understanding of biological processes. However, the growing
enormity of available data on the networks becomes a computational challenge
for the conventional alignment algorithms. Quantum algorithms generally provide
greater efficiency over their classical counterparts in solving various
problems.
  One of such algorithms is the quantum phase estimation algorithm which
generates the principal eigenvector of a stochastic matrix with probability
one.
  Using the quantum phase estimation algorithm, we introduce a quantum
computing approach for the alignment of protein-protein interaction networks by
following the classical algorithm IsoRank which uses the principal eigenvector
of the stochastic matrix representing the Kronecker product of the normalized
adjacency matrices of networks for the pairwise alignment. We also present a
greedy quantum measurement scheme to efficiently procure the alignment from the
output state of the phase estimation algorithm where the eigenvector is encoded
as the amplitudes of this state. The complexity of the quantum approach
outperforms the classical running time.


  Numerical methods are proposed for an advanced Poisson-Nernst-Planck-Fermi
(PNPF) model for studying ion transport through biological ion channels. PNPF
contains many more correlations than most models and simulations of channels,
because it includes water and calculates dielectric properties consistently as
outputs. This model accounts for the steric effect of ions and water molecules
with different sizes and interstitial voids, the correlation effect of crowded
ions with different valences, and the screening effect of polarized water
molecules in an inhomogeneous aqueous electrolyte. The steric energy is shown
to be comparable to the electrical energy under physiological conditions,
demonstrating the crucial role of the excluded volume of particles and the
voids in the natural function of channel proteins. Water is shown to play a
critical role in both correlation and steric effects in the model. We extend
the classical Scharfetter-Gummel (SG) method for semiconductor devices to
include the steric potential for ion channels, which is a fundamental physical
property not present in semiconductors. Together with a simplified matched
interface and boundary (SMIB) method for treating molecular surfaces and
singular charges of channel proteins, the extended SG method is shown to
exhibit important features in flow simulations such as optimal convergence,
efficient nonlinear iterations, and physical conservation. Two different
methods --- called the SMIB and multiscale methods --- are proposed for two
different types of channels, namely, the gramicidin A channel and an L-type
calcium channel, depending on whether water is allowed to pass through the
channel. The PNPF currents are in accord with the experimental I-V (V for
applied voltages) data of the gramicidin A channel and I-C (C for bath
concentrations) data of the calcium channel.


  An epidemic model with births and deaths is considered on a two dimensional
LxL lattice. Each individual can have global infective contacts according to
the standard SIR model rules or local infective contacts with its nearest
neighbors. We propose a deterministic approach to this model and verified that
there is a good agreement with the stochastic simulations for different
situations of the disease transmission and parameters corresponding to
pertussis and rubella in the prevaccine era.


  In complex environments, there are costs to both ignorance and perception. An
organism needs to track fitness-relevant information about its world, but the
more information it tracks, the more resources it must devote to memory and
processing. Rate-distortion theory shows that, when errors are allowed,
remarkably efficient internal representations can be found by
biologically-plausible hill-climbing mechanisms. We identify two regimes: a
high-fidelity regime where perceptual costs scale logarithmically with
environmental complexity, and a low-fidelity regime where perceptual costs are,
remarkably, independent of the environment. When environmental complexity is
rising, Darwinian evolution should drive organisms to the threshold between the
high- and low-fidelity regimes. Organisms that code efficiently will find
themselves able to make, just barely, the most subtle distinctions in their
environment.


  The article presents the theoretical foundations of the algorithm for
calculating the number of different genomes in the medium under study and of
two algorithms for determining the presence of a particular (known) genome in
this medium. The approach is based on the analysis of the compositional spectra
of subsequently sequenced samples of the medium. The theoretical estimations
required for the implementation of the algorithms are obtained.


  The Bak--Sneppen model is an abstract representation of a biological system
that evolves according to the Darwinian principles of random mutation and
selection. The species in the system are characterized by a numerical fitness
value between zero and one. We show that in the case of five species the
steady-state fitness distribution can be obtained as a solution to a linear
differential equation of order five with hypergeometric coefficients. Similar
representations for the asymptotic fitness distribution in larger systems may
help pave the way towards a resolution of the question of whether or not, in
the limit of infinitely many species, the fitness is asymptotically uniformly
distributed on the interval [f,1] with f approximately equal to 2/3.


  Characterizing macromolecular kinetics from molecular dynamics (MD)
simulations requires a distance metric that can distinguish
slowly-interconverting states. Here we build upon diffusion map theory and
define a kinetic distance for irreducible Markov processes that quantifies how
slowly molecular conformations interconvert. The kinetic distance can be
computed given a model that approximates the eigenvalues and eigenvectors
(reaction coordinates) of the MD Markov operator. Here we employ the
time-lagged independent component analysis (TICA). The TICA components can be
scaled to provide a kinetic map in which the Euclidean distance corresponds to
the kinetic distance. As a result, the question of how many TICA dimensions
should be kept in a dimensionality reduction approach becomes obsolete, and one
parameter less needs to be specified in the kinetic model construction. We
demonstrate the approach using TICA and Markov state model (MSM) analyses for
illustrative models, protein conformation dynamics in bovine pancreatic trypsin
inhibitor and protein-inhibitor association in trypsin and benzamidine.


  Ecological studies have now gone beyond measures of species turnover towards
measures of phylogenetic and functional dissimilarity with a main objective:
disentangling the processes that drive species distributions from local to
broad scales. A fundamental difference between phylogenetic and functional
analyses is that phylogeny is intrinsically dependent on a tree-like structure.
When the branches of a phylogenetic tree have lengths, then each evolutionary
unit on these branches can be considered as a basic entity on which
dissimilarities among sites should be measured. Several of the recent measures
of phylogenetic dissimilarities among sites thus are traditional dissimilarity
indices where species are replaced by evolutionary units. The resulting indices
were named PD-dissimilarity indices. Here I review and compare indices and
ordination approaches that, although first developed to analyse the differences
in the species compositions of sites, can be adapted to describe
PD-dissimilarities among sites, thus revealing how lineages are distributed
along environmental gradients, or among habitats or regions. As an
illustration, I show that the amount of PD-dissimilarities among the main
habitats of a disturbance gradient in Selva Lacandona of Chiapas, Mexico is
strongly dependent on whether species are weighted by their abundance or not,
and on the index used to measure PD-dissimilarity. Overall, the family of
PD-dissimilarity indices has a critical potential for future analyses of
phylogenetic diversity as it benefits from decades of research on the measure
of species dissimilarity. I provide clues to help to choose among many
potential indices, identifying which indices satisfy minimal basis properties,
and analysing their sensitivity to abundance, size, diversity, and joint
absences.


  Understanding how the structure of community interactions is modified by
coevolution is vital for understanding system responses to change at all
scales. However, in absence of a group selection process, collective community
behaviours cannot be organised or adapted in a Darwinian sense. An open
question thus persists: are there alternative organising principles that enable
us to understand how coevolution of component species creates complex
collective behaviours exhibited at the community level? We address this issue
using principles from connectionist learning, a discipline with well-developed
theories of emergent behaviours in simple networks. We identify conditions
where selection on ecological interactions is equivalent to 'unsupervised
learning' (a simple type of connectionist learning) and observe that this
enables communities to self organize without community-level selection. Despite
not being a Darwinian unit, ecological communities can behave like
connectionist learning systems, creating internal organisation that habituates
to past environmental conditions and actively recalling those conditions.


  Biological processes involve a variety of spatial and temporal scales. A
holistic understanding of many biological processes therefore requires
multi-scale models which capture the relevant properties on all these scales.
In this manuscript we review mathematical modelling approaches used to describe
the individual spatial scales and how they are integrated into holistic models.
We discuss the relation between spatial and temporal scales and the implication
of that on multi-scale modelling. Based upon this overview over
state-of-the-art modelling approaches, we formulate key challenges in
mathematical and computational modelling of biological multi-scale and
multi-physics processes. In particular, we considered the availability of
analysis tools for multi-scale models and model-based multi-scale data
integration. We provide a compact review of methods for model-based data
integration and model-based hypothesis testing. Furthermore, novel approaches
and recent trends are discussed, including computation time reduction using
reduced order and surrogate models, which contribute to the solution of
inference problems. We conclude the manuscript by providing a few ideas for the
development of tailored multi-scale inference methods.


  Given two phylogenetic trees on the same set of taxa X, the maximum parsimony
distance d_MP is defined as the maximum, ranging over all characters c on X, of
the absolute difference in parsimony score induced by c on the two trees. In
this note we prove that for binary trees there exists a character achieving
this maximum that is convex on one of the trees (i.e. the parsimony score
induced on that tree is equal to the number of states in the character minus 1)
and such that the number of states in the character is at most 7d_MP - 5. This
is the first non-trivial bound on the number of states required by optimal
characters, convex or otherwise. The result potentially has algorithmic
significance because, unlike general characters, convex characters with a
bounded number of states can be enumerated in polynomial time.


  In lymphoma, mutations in genes of histone modifying proteins are frequently
observed. Notably, somatic mutations in the activatory histone modification
writing protein MLL2 and the repressive modification writer EZH2 are the most
frequent. Gain of function mutations are typically detected in EZH2 whilst MLL2
mutations are usually observed as conferring a homozygous loss of function. The
gain-of-function mutations in EZH2 provide an obvious target for the
development of inhibitors with therapeutic potential. To counter the loss of
functional MLL2 protein, we computationally predicted compounds that are able
to modulate the reader of the corresponding modifications, BPTF, to recognize
other forms of the histone H3 lysine 4, instead of the tri-methylated form
normally produced by MLL2. By forming a synthetic triple-complex of a compound,
the histone H3 tail and BPTF we potentially circumvent the requirement for
functional MLL2 methyl-transferase through the modulation of BPTF activity.
Here we show a proof-of-principle that special compounds, named variators, can
reprogram selectivity of protein binding and thus create artificial regulatory
pathways which can have a potential therapeutic role. A therapeutic role of
BPTF variators may extend to other diseases that involve loss of MLL2 function,
such as Kabuki syndrome or the aberrant functioning of H3K4 modification as
observed in Huntington disease and in memory formation.


  Under the multispecies coalescent model of molecular evolution, gene trees
have independent evolutionary histories within a shared species tree. In
comparison, supermatrix concatenation methods assume that gene trees share a
single common genealogical history, thereby equating gene coalescence with
species divergence. The multispecies coalescent is supported by previous
studies which found that its predicted distributions fit empirical data, and
that concatenation is not a consistent estimator of the species tree. *BEAST, a
fully Bayesian implementation of the multispecies coalescent, is popular but
computationally intensive, so the increasing size of phylogenetic data sets is
both a computational challenge and an opportunity for better systematics. Using
simulation studies, we characterize the scaling behaviour of *BEAST, and enable
quantitative prediction of the impact increasing the number of loci has on both
computational performance and statistical accuracy. Follow up simulations over
a wide range of parameters show that the statistical performance of *BEAST
relative to concatenation improves both as branch length is reduced and as the
number of loci is increased. Finally, using simulations based on estimated
parameters from two phylogenomic data sets, we compare the performance of a
range of species tree and concatenation methods to show that using *BEAST with
tens of loci can be preferable to using concatenation with thousands of loci.
Our results provide insight into the practicalities of Bayesian species tree
estimation, the number of loci required to obtain a given level of accuracy and
the situations in which supermatrix or summary methods will be outperformed by
the fully Bayesian multispecies coalescent.


  Frequency dependent selection and demographic fluctuations play important
roles in evolutionary and ecological processes. Under frequency dependent
selection, the average fitness of the population may increase or decrease based
on interactions between individuals within the population. This should be
reflected in fluctuations of the population size even in constant environments.
Here, we propose a stochastic model, which naturally combines these two
evolutionary ingredients by assuming frequency dependent competition between
different types in an individual-based model. In contrast to previous game
theoretic models, the carrying capacity of the population and thus the
population size is determined by pairwise competition of individuals mediated
by evolutionary games and demographic stochasticity. In the limit of infinite
population size, the averaged stochastic dynamics is captured by the
deterministic competitive Lotka-Volterra equations. In small populations,
demographic stochasticity may instead lead to the extinction of the entire
population. As the population size is driven by the fitness in evolutionary
games, a population of cooperators is less prone to go extinct than a
population of defectors, whereas in the usual systems of fixed size, the
population would thrive regardless of its average payoff.


  Computational prediction of protein structures is a difficult task, which
involves fast and accurate evaluation of candidate model structures. We propose
to enhance single model quality assessment with a functionality evaluation
phase for proteins whose quantitative functional characteristics are known. In
particular, this idea can be applied to evaluation of structural models of ion
channels, whose main function - conducting ions - can be quantitatively
measured with the patch-clamp technique providing the current-voltage
characteristics. The study was performed on a set of KcsA channel models
obtained from complete and incomplete contact maps. A fast continuous
electrodiffusion model was used for calculating the current-voltage
characteristics of structural models. We found that the computed charge
selectivity and total current were sensitive to structural and electrostatic
quality of models. In practical terms, we show that evaluating predicted
conductance values is an appropriate method to eliminate modes with an occluded
pore or with multiple erroneously created pores. Moreover, filtering models on
the basis of their predicted charge selectivity results in a substantial
enrichment of the candidate set in highly accurate models. In addition to being
a proof of the concept, our function-oriented single model quality assessment
tool can be directly applied for evaluation of structural models of
strongly-selective protein channels. Finally, our work raises an important
question whether a computational validation of functionality should not be
included in the evaluation process of structural models, whenever possible.


  Previously reported [http://arxiv.org/abs/1506.06433] reprogramming of
substrate specificity of H3K4Me3 epigenetic marks reading PHD domain of BPTF
protein illustrates therapeutic potential of a new class of non-inhibitor small
organic compounds - variators. Here we address the question about
reproducibility of rational design of variators by reprogramming of the second
epigenetic marks reading domain of BPTF protein - bromodomain. Bromodomain of
BPTF binds to epigenetic marks in form of acetylated lysine of histone H4
(H4K12Ac, H4K16Ac and H4K20Ac), which physicochemical properties and binding
mode differs considerably from those of methylated H3K4 marks. Thus, detailed
description of computational approach for reprogramming of bromodomain
substrate specificity illustrates both general and target specific attributes
of computer aided variators design.


  We present a novel experimental setup in which magnetic and optical tweezers
are combined for torque and force transduction onto single filamentous
molecules in a transverse configuration to allow simultaneous mechanical
measurement and manipulation. Previously we have developed a super-resolution
imaging module which in conjunction with advanced imaging techniques such as
Blinking assisted Localisation Microscopy (BaLM) achieves localisation
precision of single fluorescent dye molecules bound to DNA of ~30 nm along the
contour of the molecule; our work here describes developments in producing a
system which combines tweezing and super-resolution fluorescence imaging. The
instrument also features an acousto-optic deflector that temporally divides the
laser beam to form multiple traps for high throughput statistics collection.
Our motivation for developing the new tool is to enable direct observation of
detailed molecular topological transformation and protein binding event
localisation in a stretching/twisting mechanical assay that previously could
hitherto only be deduced indirectly from the end-to-end length variation of
DNA. Our approach is simple and robust enough for reproduction in the lab
without the requirement of precise hardware engineering, yet is capable of
unveiling the elastic and dynamic properties of filamentous molecules that have
been hidden using traditional tools.


  Gene regulation relies on the specificity of transcription factor (TF) - DNA
interactions. In equilibrium, limited specificity may lead to crosstalk: a
regulatory state in which a gene is either incorrectly activated due to
noncognate TF-DNA interactions or remains erroneously inactive. We present a
tractable biophysical model of global crosstalk, where many genes are
simultaneously regulated by many TFs. We show that in the simplest regulatory
scenario, a lower bound on crosstalk severity can be analytically derived
solely from the number of (co)regulated genes and a suitable parameter that
describes binding site similarity. Estimates show that crosstalk could present
a significant challenge for organisms with low-specificity TFs, such as
metazoans, unless they use appropriate regulation schemes. Strong cooperativity
substantially decreases crosstalk, while joint regulation by activators and
repressors, surprisingly, does not; moreover, certain microscopic details about
promoter architecture emerge as globally important determinants of crosstalk
strength. Our results suggest that crosstalk imposes a new type of global
constraint on the functioning and evolution of regulatory networks, which is
qualitatively distinct from the known constraints acting at the level of
individual gene regulatory elements.


  Stability is a desirable property of complex ecosystems. If a community of
interacting species is at a stable equilibrium point then it is able to
withstand small perturbations to component species' abundances without
suffering adverse effects. In ecology, the Jacobian matrix evaluated at an
equilibrium point is known as the community matrix, which describes the
population dynamics of interacting species. A system's asymptotic short- and
long-term behaviour can be determined from eigenvalues derived from the
community matrix. Here we use results from the theory of pseudospectra to
describe intermediate, transient dynamics. We first recover the established
result that the transition from stable to unstable dynamics includes a region
of `transient instability', where the effect of a small perturbation to
species' abundances---to the population vector---is amplified before ultimately
decaying. Then we show that the shift from stability to transient instability
can be affected by uncertainty in, or small changes to, entries in the
community matrix, and determine lower and upper bounds to the maximum amplitude
of perturbations to the population vector. Of five different types of community
matrix, we find that amplification is least severe when predator-prey
interactions dominate. This analysis is relevant to other systems whose
dynamics can be expressed in terms of the Jacobian matrix. Our results will
lead to improved understanding of how multiple perturbations to a complex
system may irrecoverably break stability.


  The origin of life can be understood mathematically to be the origin of
information that can replicate. The likelihood that entropy spontaneously
becomes information can be calculated from first principles, and depends
exponentially on the amount of information that is necessary for replication.
We do not know what the minimum amount of information for self-replication is
because it must depend on the local chemistry, but we can study how this
likelihood behaves in different known chemistries, and we can study ways in
which this likelihood can be enhanced. Here we present evidence from numerical
simulations (using the digital life chemistry "Avida") that using a biased
probability distribution for the creation of monomers (the "biased typewriter")
can exponentially increase the likelihood of spontaneous emergence of
information from entropy. We show that this likelihood may depend on the length
of the sequence that the information is embedded in, but in a non-trivial
manner: there may be an optimum sequence length that maximizes the likelihood.
We conclude that the likelihood of spontaneous emergence of self-replication is
much more malleable than previously thought, and that the biased probability
distributions of monomers that are the norm in biochemistry may significantly
enhance these likelihoods


  The Wright-Fisher family of diffusion processes is a widely used class of
evolutionary models. However, simulation is difficult because there is no known
closed-form formula for its transition function. In this article we demonstrate
that it is in fact possible to simulate exactly from a broad class of
Wright-Fisher diffusion processes and their bridges. For those diffusions
corresponding to reversible, neutral evolution, our key idea is to exploit an
eigenfunction expansion of the transition function; this approach even applies
to its infinite-dimensional analogue, the Fleming-Viot process. We then develop
an exact rejection algorithm for processes with more general drift functions,
including those modelling natural selection, using ideas from retrospective
simulation. Our approach also yields methods for exact simulation of the moment
dual of the Wright-Fisher diffusion, the ancestral process of an infinite-leaf
Kingman coalescent tree. We believe our new perspective on diffusion simulation
holds promise for other models admitting a transition eigenfunction expansion.


  Biopolymers are characterized by heterogeneous interactions, and usually
perform their biological tasks forming contacts within domains of limited size.
Combining polymer theory with a replica approach, we study the scaling
properties of the probability of contact formation in random heteropolymers as
a function of their linear distance. It is found that close or above the
theta--point, it is possible to define a contact probability which is typical
(i.e. "self-averaging") for different realizations of the heterogeneous
interactions, and which displays an exponential cut--off, dependent on
temperature and on the interaction range. In many cases this cut--off is
comparable with the typical sizes of domains in biopolymers. While it is well
known that disorder causes interesting effects at low temperature, the behavior
elucidated in the present study is an example of a non--trivial effect at high
temperature.


  The Price equation shows the unity between the fundamental expressions of
change in biology, in information and entropy descriptions of populations, and
in aspects of thermodynamics. The Price equation partitions the change in the
average value of a metric between two populations. A population may be composed
of organisms or particles or any members of a set to which we can assign
probabilities. A metric may be biological fitness or physical energy or the
output of an arbitrarily complicated function that assigns quantitative values
to members of the population. The first part of the Price equation describes
how directly applied forces change the probabilities assigned to members of the
population when holding constant the metrical values of the members---a fixed
metrical frame of reference. The second part describes how the metrical values
change, altering the metrical frame of reference. In canonical examples, the
direct forces balance the changing metrical frame of reference, leaving the
average or total metrical values unchanged. In biology, relative reproductive
success (fitness) remains invariant as a simple consequence of the conservation
of total probability. In physics, systems often conserve total energy.
Nonconservative metrics can be described by starting with conserved metrics,
and then studying how coordinate transformations between conserved and
nonconserved metrics alter the geometry of the dynamics and the aggregate
values of populations. From this abstract perspective, key results from
different subjects appear more simply as universal geometric principles for the
dynamics of populations subject to the constraints of particular conserved
quantities


  I develop a framework for interpreting the forces that act on any population
described by frequencies. The conservation of total frequency, or total
probability, shapes the characteristics of force. I begin with Fisher's
fundamental theorem of natural selection. That theorem partitions the total
evolutionary change of a population into two components. The first component is
the partial change caused by the direct force of natural selection, holding
constant all aspects of the environment. The second component is the partial
change caused by the changing environment. I demonstrate that Fisher's
partition of total change into the direct force of selection and the forces
from the changing environmental frame of reference is identical to d'Alembert's
principle of mechanics, which separates the work done by the direct forces from
the work done by the inertial forces associated with the changing frame of
reference. In d'Alembert's principle, there exist inertial forces from a change
in the frame of reference that exactly balance the direct forces. I show that
the conservation of total probability strongly shapes the form of the balance
between the direct and inertial forces. I then use the strong results for
conserved probability to obtain general results for the change in any system
quantity, such as biological fitness or energy. Those general results derive
from simple coordinate changes between frequencies and system quantities.
Ultimately, d'Alembert's separation of direct and inertial forces provides deep
conceptual insight into the interpretation of forces and the unification of
disparate fields of study.


  Cancer research reflects an implicit conflict. On the one hand, there is an
overwhelming desire to control the disease. We all wish that. On the other
hand, we would like to understand why cancer follows so many clearly defined
yet puzzling patterns. Why is there such regularity in the rates of
progression? Why do different tissues vary so much? There should, of course, be
no conflict between control and understanding. But the history of cancer
research seems to say that those different goals remain oddly estranged. Peto's
1977 article locates the seeds of this conflict most clearly. He describes what
is still the most powerful theoretical perspective for analyzing the causes of
cancer. He presents many key unsolved puzzles within that context. He also says
why most cancer researchers are not interested in these fundamental issues. The
subsequent decades of research grew around this rift, blindly, in the way that
research disciplines often grow. Let us revisit Peto, almost 40 years ago. We
can learn much about the current nature of cancer research.


  Evolutionary game theory is a powerful framework for studying evolution in
populations of interacting individuals. A common assumption in evolutionary
game theory is that interactions are symmetric, which means that the players
are distinguished by only their strategies. In nature, however, the microscopic
interactions between players are nearly always asymmetric due to environmental
effects, differing baseline characteristics, and other possible sources of
heterogeneity. To model these phenomena, we introduce into evolutionary game
theory two broad classes of asymmetric interactions: ecological and genotypic.
Ecological asymmetry results from variation in the environments of the players,
while genotypic asymmetry is a consequence of the players having differing
baseline genotypes. We develop a theory of these forms of asymmetry for games
in structured populations and use the classical social dilemmas, the Prisoner's
Dilemma and the Snowdrift Game, for illustrations. Interestingly, asymmetric
games reveal essential differences between models of genetic evolution based on
reproduction and models of cultural evolution based on imitation that are not
apparent in symmetric games.


  Spatial patterning can be crucially important for understanding the behavior
of interacting populations. Here we investigate a simple model of parasite and
host populations in which parasites are random walkers that must come into
contact with a host in order to reproduce. We focus on the spatial arrangement
of parasites around a single host, and we derive using analytics and numerical
simulations the necessary conditions placed on the parasite fecundity and
lifetime for the populations long-term survival. We also show that the parasite
population can be pushed to extinction by a large drift velocity, but,
counterintuitively, a small drift velocity generally increases the parasite
population.


  Recently, Hajirasouliha and Raphael (WABI 2014) proposed a model for
deconvoluting mixed tumor samples measured from a collection of high-throughput
sequencing reads. This is related to understanding tumor evolution and critical
cancer mutations. In short, their formulation asks to split each row of a
binary matrix so that the resulting matrix corresponds to a perfect phylogeny
and has the minimum number of rows among all matrices with this property. In
this paper we disprove several claims about this problem, including an
NP-hardness proof of it. However, we show that the problem is indeed NP-hard,
by providing a different proof. We also prove NP-completeness of a variant of
this problem proposed in the same paper. On the positive side, we propose an
efficient (though not necessarily optimal) heuristic algorithm based on
coloring co-comparability graphs, and a polynomial time algorithm for solving
the problem optimally on matrix instances in which no column is contained in
both columns of a pair of conflicting columns. Implementations of these
algorithms are freely available at
https://github.com/alexandrutomescu/MixedPerfectPhylogeny


  Phylogenetic comparative analysis is an approach to inferring evolutionary
process from a combination of phylogenetic and phenotypic data. The last few
years have seen increasingly sophisticated models employed in the evaluation of
more and more detailed evolutionary hypotheses, including adaptive hypotheses
with multiple selective optima and hypotheses with rate variation within and
across lineages. The statistical performance of these sophisticated models has
received relatively little systematic attention, however. We conducted an
extensive simulation study to quantify the statistical properties of a class of
models toward the simpler end of the spectrum that model phenotypic evolution
using Ornstein-Uhlenbeck processes. We focused on identifying where, how, and
why these methods break down so that users can apply them with greater
understanding of their strengths and weaknesses. Our analysis identifies three
key determinants of performance: a discriminability ratio, a signal-to-noise
ratio, and the number of taxa sampled. Interestingly, we find that
model-selection power can be high even in regions that were previously thought
to be difficult, such as when tree size is small. On the other hand, we find
that model parameters are in many circumstances difficult to estimate
accurately, indicating a relative paucity of information in the data relative
to these parameters. Nevertheless, we note that accurate model selection is
often possible when parameters are only weakly identified. Our results have
implications for more sophisticated methods inasmuch as the latter are
generalizations of the case we study.


  The ionic atmospheres around nucleic acids play important roles in biological
function. Large-scale explicit solvent simulations coupled to experimental
assays such as anomalous small-angle X-ray scattering (ASAXS) can provide
important insights into the structure and energetics of such atmospheres but
are time- and resource-intensive. In this paper, we use classical density
functional theory (cDFT) to explore the balance between ion-DNA, ion-water, and
ion-ion interactions in ionic atmospheres of RbCl, SrCl$_2$, and CoHexCl$_3$
(cobalt hexammine chloride) around a B-form DNA molecule. The accuracy of the
cDFT calculations was assessed by comparison between simulated and experimental
ASAXS curves, demonstrating that an accurate model should take into account
ion-ion correlation and ion hydration forces, DNA topology, and the discrete
distribution of charges on DNA strands. As expected, these calculations
revealed significant differences between monovalent, divalent, and trivalent
cation distributions around DNA. About half of the DNA-bound Rb$^+$ ions
penetrate into the minor groove of the DNA and half adsorb on the DNA strands.
The fraction of cations in the minor groove decreases for the larger Sr$^{2+}$
ions and becomes zero for CoHex$^{3+}$ ions, which all adsorb on the DNA
strands. The distribution of CoHex$^{3+}$ ions is mainly determined by Coulomb
and steric interactions, while ion-correlation forces play a central role in
the monovalent Rb$^+$ distribution and a combination of ion-correlation and
hydration forces affect the Sr$^{2+}$ distribution around DNA.


  Summary: Both theory and empirical evidence indicate that phylogenies (trees)
of different genes (loci) do not display precisely matched topologies. This
phylogenetic incongruence is attributed to the reticulated evolutionary history
of most species due to meiotic sexual recombination in eukaryotes, or
horizontal transfers of genetic materials in prokaryotes. Nonetheless, most
genes do display topologically related phylogenies; this implies they form
cohesive subsets (clusters). In this work, we compare popular clustering
methods, and show how the performance of the normalized cut framework is
efficient and statistically accurate when obtaining clusters on the set of gene
trees based on the geodesic distance between them over the
Billera-Holmes-Vogtmann (BHV) tree space. We proceed to present a computational
study on the performance of different clustering methods with and without
preprocessing under different distance metrics and using a series of dimension
reduction techniques. Results: First, we show using simulated data that indeed
the Ncut framework accurately clusters the set of gene trees given a species
tree under the coalescent process. We then depict the success of our framework
by comparing its performance to other clustering techniques, including k-means
and hierarchical clustering. The main computational results can be summarized
to the stellar performance of the Ncut framework even without dimension
reduction, the similar performance portrayed by k-means and Ncut under most
dimension reduction schemes, the utter failure of hierarchical clustering to
accurately capture clusters, as well as the significantly better performance of
the NJp method, as compared to MLE.


  Cells moving on a two dimensional substrate generate motion by polymerizing
actin filament networks inside a flat membrane protrusion. New filaments are
generated by branching off existing ones, giving rise to branched network
structures. We investigate the force-extension relation of branched filaments,
grafted on an elastic structure at one end and pushing with the free ends
against the leading edge cell membrane. Single filaments are modeled as
worm-like chains, whose thermal bending fluctuations are restricted by the
leading edge cell membrane, resulting in an effective force. Branching can
increase the stiffness considerably; however the effect depends on branch point
position and filament orientation, being most pronounced for intermediate tilt
angles and intermediate branch point positions. We describe filament networks
without cross-linkers to focus on the effect of branching. We use randomly
positioned branch points, as generated in the process of treadmilling, and
orientation distributions as measured in lamellipodia. These networks reproduce
both the weak and strong force response of lamellipodia as measured in
force-velocity experiments. We compare properties of branched and unbranched
networks. The ratio of the network average of the force per branched filament
to the average force per unbranched filament depends on the orientation
distribution of the filaments. The ratio exhibits compression dependence and
may go up to about 4.5 in networks with a narrow orientation distribution. With
orientation distributions measured in lamellipodia, it is about two and
essentially independent from network compression, graft elasticity and filament
persistence length.


  Motivation: Optimizing seed selection is an important problem in read
mapping. The number of non-overlapping seeds a mapper selects determines the
sensitivity of the mapper while the total frequency of all selected seeds
determines the speed of the mapper. Modern seed-and-extend mappers usually
select seeds with either an equal and fixed-length scheme or with an inflexible
placement scheme, both of which limit the potential of the mapper to select
less frequent seeds to speed up the mapping process. Therefore, it is crucial
to develop a new algorithm that can adjust both the individual seed length and
the seed placement, as well as derive less frequent seeds.
  Results: We present the Optimal Seed Solver (OSS), a dynamic programming
algorithm that discovers the least frequently-occurring set of x seeds in an
L-bp read in $O(x \times L)$ operations on average and in $O(x \times L^{2})$
operations in the worst case. We compared OSS against four state-of-the-art
seed selection schemes and observed that OSS provides a 3-fold reduction of
average seed frequency over the best previous seed selection optimizations.


  Summary: BGT is a compact format, a fast command line tool and a simple web
application for efficient and convenient query of whole-genome genotypes and
frequencies across tens to hundreds of thousands of samples. On real data, it
encodes the haplotypes of 32,488 samples across 39.2 million SNPs into a 7.4GB
database and decodes a couple of hundred million genotypes per CPU second. The
high performance enables real-time responses to complex queries.
  Availability and implementation: https://github.com/lh3/bgt
  Contact: hengli@broadinstitute.org


  This article shows how to specify and construct a discrete, stochastic,
continuous-time model specifically for ecological systems. The model is more
broad than typical chemical kinetics models in two ways. First, using
time-dependent hazard rates simplifies the process of making models more
faithful. Second, the state of the system includes individual traits and use of
environmental resources. The models defined here focus on taking survival
analysis of observations in the field and using the measured hazard rates to
generate simulations which match exactly what was measured.


  The linking genotype to phenotype is the fundamental aim of modern genetics.
We focus on study of links between gene expression data and phenotype data
through integrative analysis. We propose three approaches.
  1) The inherent complexity of phenotypes makes high-throughput phenotype
profiling a very difficult and laborious process. We propose a method of
automated multi-dimensional profiling which uses gene expression similarity.
Large-scale analysis show that our method can provide robust profiling that
reveals different phenotypic aspects of samples. This profiling technique is
also capable of interpolation and extrapolation beyond the phenotype
information given in training data. It can be used in many applications,
including facilitating experimental design and detecting confounding factors.
  2) Phenotype association analysis problems are complicated by small sample
size and high dimensionality. Consequently, phenotype-associated gene subsets
obtained from training data are very sensitive to selection of training
samples, and the constructed sample phenotype classifiers tend to have poor
generalization properties. To eliminate these obstacles, we propose a novel
approach that generates sequences of increasingly discriminative gene cluster
combinations. Our experiments on both simulated and real datasets show robust
and accurate classification performance.
  3) Many complex phenotypes, such as cancer, are the product of not only gene
expression, but also gene interaction. We propose an integrative approach to
find gene network modules that activate under different phenotype conditions.
Using our method, we discovered cancer subtype-specific network modules, as
well as the ways in which these modules coordinate. In particular, we detected
a breast-cancer specific tumor suppressor network module with a hub gene,
PDGFRL, which may play an important role in this module.


  Sexually transmitted infections (STIs) continue to present a complex and
costly challenge to public health programs. The preferences and social dynamics
of a population can have a large impact on the course of an outbreak as well as
the effectiveness of interventions intended to influence individual behavior.
In addition, individuals may alter their sexual behavior in response to the
presence of STIs, creating a feedback loop between transmission and behavior.
We investigate the consequences of modeling the interaction between STI
transmission and prophylactic use with a model that links a
Susceptible-Infectious-Susceptible (SIS) system to evolutionary game dynamics
that determine the effective contact rate. The combined model framework allows
us to address protective behavior by both infected and susceptible individuals.
Feedback between behavioral adaptation and prevalence creates a wide range of
dynamic behaviors in the combined model, including damped and sustained
oscillations as well as bistability, depending on the behavioral parameters and
disease growth rate. We found that disease extinction is possible for multiple
regions where R0 > 1, due to behavior adaptation driving the epidemic downward,
although conversely endemic prevalence for arbitrarily low R0 is also possible
if contact rates are sufficiently high. We also tested how model
misspecification might affect disease forecasting and estimation of the model
parameters and R0. We found that alternative models that neglect the behavioral
feedback or only consider behavior adaptation by susceptible individuals can
potentially yield misleading parameter estimates or omit significant features
of the disease trajectory.


  The dynamics of adaptation is difficult to predict because it is highly
stochastic even in large populations. The uncertainty emerges from number
fluctuations, called genetic drift, arising in the small number of particularly
fit individuals of the population. Random genetic drift in this evolutionary
vanguard also limits the speed of adaptation, which diverges in deterministic
models that ignore these chance effects. Several approaches have been developed
to analyze the crucial role of noise on the expected dynamics of adaptation,
including the mean fitness of the entire population, or the fate of newly
arising beneficial deleterious mutations. However, very little is known about
how genetic drift causes fluctuations to emerge on the population level,
including fitness distribution variations and speed variations. Yet, these
phenomena control the replicability of experimental evolution experiments and
are key to a truly predictive understanding of evolutionary processes. Here, we
develop an exact approach to these emergent fluctuations by a combination of
computational and analytical methods. We show, analytically, that the infinite
hierarchy of moment equations can be closed at any arbitrary order by a
suitable choice of a dynamical constraint. This constraint regulates (rather
than fixes) the population size, accounting for resource limitations. The
resulting linear equations, which can be accurately solved numerically, exhibit
fluctuation-induced terms that amplify short-distance correlations and suppress
long-distance ones. Importantly, by accounting for the dynamics of
sub-populations, we provide a systematic route to key population genetic
quantities, such as fixation probabilities and decay rates of the genetic
diversity.