  Vector-borne diseases with reservoir cycles are complex to understand because
new infections come from contacts of the vector with humans and different
reservoirs. In this scenario, the basic reproductive number $\mathcal{R}^h_0$
of the system where the reservoirs are not included could turn out to be less
than one, yet, an endemic equilibrium be observed. Indeed, when the reservoirs
are taken back into account, the basic reproductive number $\mathcal{R}_0^r$,
of only vectors and reservoirs, explains the endemic state. Furthermore,
reservoirs cycles with a small basic reproductive number could contribute to
reach an endemic state in the human cycle. Therefore, when controlling for the
spread of a disease, it could not be enough to focus on specific reservoir
cycles or only on the vector. In this work, we created a simple epidemiological
model with a network of reservoirs where $\mathcal{R}_0^r$ is a bifurcation
parameter of the system, explaining disease endemicity in the absence of a
strong reservoir cycle. This simple model may help to explain transmission
dynamics of diseases such as Chagas, Leishmaniasis and Dengue.


  Since the late `60s, various genome evolutionary models have been proposed to
predict the evolution of a DNA sequence as the generations pass. Most of these
models are based on nucleotides evolution, so they use a mutation matrix of
size 4x4. They encompass for instance the well-known models of Jukes and
Cantor, Kimura, and Tamura. By essence, all of these models relate the
evolution of DNA sequences to the computation of the successive powers of a
mutation matrix. To make this computation possible, particular forms for the
mutation matrix are assumed, which are not compatible with mutation rates that
have been recently obtained experimentally on gene ura3 of the Yeast
Saccharomyces cerevisiae. Using this experimental study, authors of this paper
have deduced a simple mutation matrice, compute the future evolution of the
rate purine/pyrimidine for ura3, investigate the particular behavior of
cytosines and thymines compared to purines, and simulate the evolution of each
nucleotide.


  We consider a discrete-time Markov chain, called fragmentation process, that
describes a specific way of successively removing objects from a linear
arrangement. The process arises in population genetics and describes the
ancestry of the genetic material of individuals in a population experiencing
recombination. We aim at the law of the process over time. To this end, we
investigate sets of realisations of this process that agree with respect to a
specific order of events and represent each such set by a rooted (binary) tree.
The probability of each tree is, in turn, obtained by M\"obius inversion on a
suitable poset of all rooted forests that can be obtained from the tree by edge
deletion; we call this poset the \textit{pruning poset}. Dependencies within
the fragments make it difficult to obtain explicit expressions for the
probabilities of the trees. We therefore construct an auxiliary process for
every given tree, which is i.i.d. over time, and which allows to give a
pathwise construction of realisations that match the tree.


  A key ecological parameter for planktonic copepods studies is their
interspecies encounter rate which is driven by their behaviour and is strongly
influenced by turbulence of the surrounding environment. A distinctive feature
of copepods motility is their ability to perform quick displacements, often
dubbed jumps, by means of powerful swimming strokes. Such a reaction has been
associated to an escape behaviour from flow disturbances due to predators or
other external dangers. In the present study, the encounter rate of copepods in
a developed turbulent flow with intensity comparable to the one found in
copepods' habitat is numerically investigated. This is done by means of a
Lagrangian copepod (LC) model that mimics the jump escape reaction behaviour
from localised high-shear rate fluctuations in the turbulent flows. Our
analysis shows that the encounter rate for copepods of typical perception
radius of ~ {\eta}, where {\eta} is the dissipative scale of turbulence, can be
increased by a factor up to ~ 100 compared to the one experienced by passively
transported fluid tracers. Furthermore, we address the effect of introducing in
the LC model a minimal waiting time between consecutive jumps. It is shown that
any encounter-rate enhancement is lost if such time goes beyond the dissipative
time-scale of turbulence, {\tau}_{\eta}. Because typically in the ocean {\eta}
~ 0.001m and {\tau}_{\eta} ~ 1s, this provides stringent constraints on the
turbulent-driven enhancement of encounter-rate due to a purely mechanical
induced escape reaction.


  Recently developed deep learning techniques have significantly improved the
accuracy of various speech and image recognition systems. In this paper we show
how to adapt some of these techniques to create a novel chained convolutional
architecture with next-step conditioning for improving performance on protein
sequence prediction problems. We explore its value by demonstrating its ability
to improve performance on eight-class secondary structure prediction. We first
establish a state-of-the-art baseline by adapting recent advances in
convolutional neural networks which were developed for vision tasks. This model
achieves 70.0% per amino acid accuracy on the CB513 benchmark dataset without
use of standard performance-boosting techniques such as ensembling or multitask
learning. We then improve upon this state-of-the-art result using a novel
chained prediction approach which frames the secondary structure prediction as
a next-step prediction problem. This sequential model achieves 70.3% Q8
accuracy on CB513 with a single model; an ensemble of these models produces
71.4% Q8 accuracy on the same test set, improving upon the previous overall
state of the art for the eight-class secondary structure problem. Our models
are implemented using TensorFlow, an open-source machine learning software
library available at TensorFlow.org; we aim to release the code for these
experiments as part of the TensorFlow repository.


  Tree containment problem is a fundamental problem in phylogenetic study, as
it is used to verify a network model. It asks whether a given network contain a
subtree that resembles a binary tree. The problem is NP-complete in general,
even in the class of binary network. Recently, it was proven to be solvable in
cubic time, and later in quadratic time for the class of general reticulation
visible networks. In this paper, we further improve the time complexity into
linear time.


  In this work we explore the dissimilarity between symmetric word pairs, by
comparing the inter-word distance distribution of a word to that of its
reversed complement. We propose a new measure of dissimilarity between such
distributions. Since symmetric pairs with different patterns could point to
evolutionary features, we search for the pairs with the most dissimilar
behaviour. We focus our study on the complete human genome and its
repeat-masked version.


  We consider a class of branching processes called Markovian binary trees, in
which the individuals lifetime and reproduction epochs are modeled using a
transient Markovian arrival process (TMAP). We estimate the parameters of the
TMAP based on population data containing information on age-specific fertility
and mortality rates. Depending on the degree of detail of the available data, a
weighted non-linear regression method or a maximum likelihood method is
applied. We discuss the optimal choice of the number of phases in the TMAP, and
we provide confidence intervals for the model outputs. The results are
illustrated using real data on human and bird populations.


  Recent 60Fe results have suggested that the estimated distances of supernovae
in the last few million years should be reduced from 100 pc to 50 pc. Two
events or series of events are suggested, one about 2.7 million years to 1.7
million years ago, and another may at 6.5 to 8.7 million years ago. We ask what
effects such supernovae are expected to have on the terrestrial atmosphere and
biota. Assuming that the Local Bubble was formed before the event being
considered, and that the supernova and the Earth were both inside a weak,
disordered magnetic field at that time, TeV-PeV cosmic rays at Earth will
increase by a factor of a few hundred. Tropospheric ionization will increase
proportionately, and the overall muon radiation load on terrestrial organisms
will increase by a factor of 150. All return to pre-burst levels within 10kyr.
In the case of an ordered magnetic field, effects depend strongly on the field
orientation. The upper bound in this case is with a largely coherent field
aligned along the line of sight to the supernova, in which case TeV-PeV cosmic
ray flux increases are 10^4; in the case of a transverse field they are below
current levels. We suggest a substantial increase in the extended effects of
supernovae on Earth and in the lethal distance estimate; more work is
needed.This paper is an explicit followup to Thomas et al. (2016). We also here
provide more detail on the computational procedures used in both works.


  We demonstrate the presence of chaos in stochastic simulations that are
widely used to study biodiversity in nature. The investigation deals with a set
of three distinct species that evolve according to the standard rules of
mobility, reproduction and predation, with predation following the cyclic rules
of the popular rock, paper and scissors game. The study uncovers the
possibility to distinguish between time evolutions that start from slightly
different initial states, guided by the Hamming distance which heuristically
unveils the chaotic behavior. The finding opens up a quantitative approach that
relates the correlation length to the average density of maxima of a typical
species, and an ensemble of stochastic simulations is implemented to support
the procedure. The main result of the work shows how a single and simple
experimental realization that counts the density of maxima associated with the
chaotic evolution of the species serves to infer its correlation length. We use
the result to investigate others distinct complex systems, one dealing with a
set of differential equations that can be used to model a diversity of natural
and artificial chaotic systems, and another one, focusing on the ocean water
level.


  Stochastic reaction-diffusion equations are a popular modelling approach for
studying interacting populations in a heterogeneous environment under the
influence of environmental fluctuations. Although the theoretical basis of
alternative models such as Fokker-Planck diffusion is not less convincing,
movement of populations is commonly modelled using the diffusion law due to
Fick. It is an interesting feature of Fokker-Planck diffusion that for
spatially varying diffusion coefficients the stationary solution is not a
homogeneous distribution; in contrast to Fickian diffusion. Instead,
concentration accumulates in regions of low diffusivity and tends to lower
levels for areas of high diffusivity. Thus, we may interpret the stationary
distribution of the Fokker-Planck diffusion as a reflection of different levels
of habitat quality. Moreover, the most common model for environmental
fluctuations, linear multiplicative noise, is based on the assumption that
individuals respond independently to stochastic environmental fluctuations. For
large population densities the assumption of independence is debatable and the
model further implies that noise intensities can increase to arbitrarily high
levels. Therefore, instead of the commonly used linear multiplicative noise
model, we implement environmental variability by an alternative nonlinear noise
term which never exceeds a certain maximum noise intensity. With Fokker-Planck
diffusion and the nonlinear noise model replacing the classical approaches we
investigate a simple invasive system based on the Lotka-Volterra competition
model. We observe that the heterogeneous stationary distribution generated by
Fokker-Planck diffusion generally facilitates the formation of segregated
habitats of resident and invader. However, this segregation can be broken by
nonlinear noise leading to coexistence of resident and invader across the whole
spatial domain.


  A trademark of eusocial insect species is reproductive division of labor, in
which workers forego their own reproduction while the queen produces almost all
offspring. The presence of the queen is key for maintaining social harmony, but
the specific role of the queen in the evolution of eusociality remains unclear.
A long-discussed scenario is that a queen either behaviorally or chemically
sterilizes her workers. However, the demographic and ecological conditions that
enable such manipulation are unknown. Accordingly, we propose a simple model of
evolutionary dynamics that is based on haplodiploid genetics. We consider a
mutation that acts in a queen, causing her to control the reproductive behavior
of her workers. Our mathematical analysis yields precise conditions for the
evolutionary emergence and stability of queen-induced worker sterility. These
conditions do not depend on the queen's mating frequency. Moreover, we find
that queen control is always established if it increases colony reproductive
efficiency and can evolve even if it decreases colony efficiency. We further
outline the conditions under which queen control is evolutionarily stable
against invasion by mutant, reproductive workers.


  This paper analyses an SIRS-type model for infectious diseases with account
for behavioural changes associated with the simultaneous spread of awareness in
the population. Two types of awareness are included into the model: private
awareness associated with direct contacts between unaware and aware
populations, and public information campaign. Stability analysis of different
steady states in the model provides information about potential spread of
disease in a population, and well as about how the disease dynamics is affected
by the two types of awareness. Numerical simulations are performed to
illustrate the behaviour of the system in different dynamical regimes.


  We review the status of protein-based molecular electronics. First we discuss
fundamental concepts of electron transfer and transport in and across proteins
and proposed mechanisms for these processes. We then describe the
immobilization of proteins to solid-state surfaces in both nanoscale and
macroscopic approaches, and highlight how different methodologies can alter
protein electronic properties. Because immobilizing proteins while retaining
biological activity is crucial to the successful development of bioelectronic
devices, we discuss this process at length. We briefly discuss computational
predictions and their link to experimental results. We then summarize how the
biological activity of immobilized proteins is beneficial for bioelectronics
devices, and how conductance measurements can shed light on protein properties.
Finally, we consider how the research to date could influence the development
of future bioelectronics devices.


  Signal transduction within biological cells is governed by networks of
interacting proteins. Communication between these proteins is mediated by
signaling molecules which bind to receptors and induce stochastic transitions
between different conformational states. Signaling is typically a cooperative
process which requires the occurrence of multiple binding events so that
reaction rates have a nonlinear dependence on the amount of signaling molecule.
It is this nonlinearity that endows biological signaling networks with robust
switch-like properties which are critical to their biological function. In this
study, we investigate how the properties of these signaling systems depend on
the network architecture. Our main result is that these nonlinear networks
exhibit bistability where the network activity can switch between states that
correspond to a low and high activity level. We show that this bistable regime
emerges at a critical coupling strength that is determined by the spectral
structure of the network. In particular, the set of nodes that correspond to
large components of the leading eigenvector of the adjacency matrix determines
the onset of bistability. Above this transition, the eigenvectors of the
adjacency matrix determine a hierarchy of clusters, defined by its spectral
properties, which are activated sequentially with increasing network activity.
We argue further that the onset of bistability occurs either continuously or
discontinuously depending upon whether the leading eigenvector is localized or
delocalized. Finally, we show that at low network coupling stochastic
transitions to the active branch are also driven by the set of nodes that
contribute more strongly to the leading eigenvector.


  We investigate theoretically if treatment alone can reduce the
schistosomiasis's prevalence in an infected population, in a long-lasting
sustainable way. We use a non-linear system of ordinary differential equations
(a SI system combined with a logistic population growth) which describes the
time evolution of the non-infected and infected populations, in terms of the
recovering, infection, and demographic rates.
  Our model leads to the conclusion that the only way to eliminate this endemic
disease is to implement public health policies aimed at both treatment and
environment.
  We apply our model to the endemic area of Virgem das Gra\c{c}as, in Brazil,
where the prevalence of Schistosiamisis in 2001 was greater than 50\%. The
epidemiological data are extracted from a longitudinal study carried on the
region between 2001 and 2010 and the demographic parameter from official
Brazilian population data. When these estimated parameters are entered, our
model gives a limit prevalence for Virgem das Gra\c{c}as of 11\%, which is
still significantly high even though treatment in accordance to the government
regulations is systematically performed. This estimative reinforces once again,
that in order to eliminate Schistomiasis, public health policies aimed at
treatment, sanitation, including snail control, and health education programs
are mandatory.


  We study the population size time series of a Neotropical small mammal with
the intent of detecting and modelling population regulation processes generated
by density-dependent factors and their possible delayed effects. The
application of analysis tools based on principles of statistical generality are
nowadays a common practice for describing these phenomena, but, in general,
they are more capable of generating clear diagnosis rather than granting
valuable modelling. For this reason, in our approach, we detect the principal
temporal structures on the bases of different correlation measures, and from
these results we build an ad-hoc minimalist autoregressive model that
incorporates the main drivers of the dynamics. Surprisingly our model is
capable of reproducing very well the time patterns of the empirical series and,
for the first time, clearly outlines the importance of the time of attaining
sexual maturity as a central temporal scale for the dynamics of this species.
In fact, an important advantage of this analysis scheme is that all the model
parameters are directly biologically interpretable and potentially measurable,
allowing a consistency check between model outputs and independent
measurements.


  We consider a class of birth-and-death processes describing a population made
of $d$ sub-populations of different types which interact with one another. The
state space is $\mathbb{Z}_+^d$ (unbounded). We assume that the population goes
almost surely to extinction, so that the unique stationary distribution is the
Dirac measure at the origin. These processes are parametrized by a scaling
parameter $K$ which can be thought as the order of magnitude of the total size
of the population at time $0$. For any fixed finite time span, it is well-known
that such processes, when renormalized by $K$, are close, in the limit
$K\to+\infty$, to the solutions of a certain differential equation in
$\mathbb{R}_+^d$ whose vector field is determined by the birth and death rates.
We consider the case where there is a unique attractive fixed point (off the
boundary of the positive orthant) for the vector field (while the origin is
repulsive). What is expected is that, for $K$ large, the process will stay in
the vicinity of the fixed point for a very long time before being absorbed at
the origin. To precisely describe this behavior, we prove the existence of a
quasi-stationary distribution (qsd). In fact, we establish a bound for the
total variation distance between the process conditioned to non-extinction
before time $t$ and the qsd. This bound is exponentially small in $t$, for
$t\gg \log K$. As a by-product, we obtain an estimate for the mean time to
extinction in the qsd. We also quantify how close is the law of the process
(not conditioned to non-extinction) either to the Dirac measure at the origin
or to the qsd, for times much larger than $\log K$ and much smaller than the
mean time to extinction, which is exponentially large as a function of $K$. Let
us stress that we are interested in what happens for finite $K$. We obtain
results much beyond what large deviation techniques could provide.


  Assuming a steady-state condition within a cell, metabolic fluxes satisfy an
under-determined linear system of stoichiometric equations. Characterizing the
space of fluxes that satisfy such equations along with given bounds (and
possibly additional relevant constraints) is considered of utmost importance
for the understanding of cellular metabolism. Extreme values for each
individual flux can be computed with Linear Programming (as Flux Balance
Analysis), and their marginal distributions can be approximately computed with
Monte-Carlo sampling. Here we present an approximate analytic method for the
latter task based on Expectation Propagation equations that does not involve
sampling and can achieve much better predictions than other existing analytic
methods. The method is iterative, and its computation time is dominated by one
matrix inversion per iteration. With respect to sampling, we show through
extensive simulation that it has some advantages including computation time,
and the ability to efficiently fix empirically estimated distributions of
fluxes.


  The stochastic dynamics of biochemical networks are usually modelled with the
chemical master equation (CME). The stationary distributions of CMEs are seldom
solvable analytically, and numerical methods typically produce estimates with
uncontrolled errors. Here, we introduce mathematical programming approaches
that yield approximations of these distributions with computable error bounds
which enable the verification of their accuracy. First, we use semidefinite
programming to compute increasingly tighter upper and lower bounds on the
moments of the stationary distributions for networks with rational
propensities. Second, we use these moment bounds to formulate linear programs
that yield convergent upper and lower bounds on the stationary distributions
themselves, their marginals and stationary averages. The bounds obtained also
provide a computational test for the uniqueness of the distribution. In the
unique case, the bounds form an approximation of the stationary distribution
with a computable bound on its error. In the non-unique case, our approach
yields converging approximations of the ergodic distributions. We illustrate
our methodology through several biochemical examples taken from the literature:
Schl\"ogl's model for a chemical bifurcation, a two-dimensional toggle switch,
a model for bursty gene expression, and a dimerisation model with multiple
stationary distributions.


  Phylogenetic networks are a generalization of phylogenetic trees that allow
for representation of reticulate evolution. Recently, a space of unrooted
phylogenetic networks was introduced, where such a network is a connected graph
in which every vertex has degree 1 or 3 and whose leaf-set is a fixed set $X$
of taxa. This space, denoted $\mathcal{N}(X)$, is defined in terms of two
operations on networks -- the nearest neighbor interchange and triangle
operations -- which can be used to transform any network with leaf set $X$ into
any other network with that leaf set. In particular, it gives rise to a metric
$d$ on $\mathcal N(X)$ which is given by the smallest number of operations
required to transform one network in $\mathcal N(X)$ into another in $\mathcal
N(X)$. The metric generalizes the well-known NNI-metric on phylogenetic trees
which has been intensively studied in the literature. In this paper, we derive
a bound for the metric $d$ as well as a related metric $d_{N\!N\!I}$ which
arises when restricting $d$ to the subset of $\mathcal{N}(X)$ consisting of all
networks with $2(|X|-1+i)$ vertices, $i \ge 1$. We also introduce two new
metrics on networks -- the SPR and TBR metrics -- which generalize the metrics
on phylogenetic trees with the same name and give bounds for these new metrics.
We expect our results to eventually have applications to the development and
understanding of network search algorithms.


  Cytosine methylation has been found to play a crucial role in various
biological processes, including a number of human diseases. The detection of
this small modification remains challenging. In this work, we computationally
explore the possibility of detecting methylated DNA strands through direct
electrical conductance measurements. Using density functional theory and the
Landauer-Buttiker method, we study the electronic properties and charge
transport through an eight base-pair methylated DNA strand and its native
counterpart. We first analyze the effect of cytosine methylation on the
tight-binding parameters of two DNA strands and then model the transmission of
the electrons and conductance through the strands both with and without
decoherence. We find that the main difference of the tight-binding parameters
between the native DNA and the methylated DNA lies in the on-site energies of
(methylated) cytosine bases. The intra- and inter- strand hopping integrals
between two nearest neighboring guanine base and (methylated) cytosine base
also change with the addition of the methyl groups. Our calculations show that
in the phase-coherent limit, the transmission of the methylated strand is close
to the native strand when the energy is nearby the highest occupied molecular
orbital level and larger than the native strand by 5 times in the bandgap. The
trend in transmission also holds in the presence of the decoherence with the
same rate. The lower conductance for the methylated strand in the experiment is
suggested to be caused by the more stable structure due to the introduction of
the methyl groups. We also study the role of the exchangecorrelation functional
and the effect of contact coupling by choosing coupling strengths ranging from
weak to strong coupling limit.


  We investigate the effects of social interactions in task al- location using
Evolutionary Game Theory (EGT). We propose a simple task-allocation game and
study how different learning mechanisms can give rise to specialised and non-
specialised colonies under different ecological conditions. By combining
agent-based simulations and adaptive dynamics we show that social learning can
result in colonies of generalists or specialists, depending on ecological
parameters. Agent-based simulations further show that learning dynamics play a
crucial role in task allocation. In particular, introspective individual
learning readily favours the emergence of specialists, while a process
resembling task recruitment favours the emergence of generalists.


  The learning rate is an information-theoretical quantity for bipartite Markov
chains describing two coupled subsystems. It is defined as the rate at which
transitions in the downstream subsystem tend to increase the mutual information
between the two subsystems, and is bounded by the dissipation arising from
these transitions. Its physical interpretation, however, is unclear, although
it has been used as a metric for the sensing performance of the downstream
subsystem. In this paper, we explore the behaviour of the learning rate for a
number of simple model systems, establishing when and how its behaviour is
distinct from the instantaneous mutual information between subsystems. In the
simplest case, the two are almost equivalent. In more complex steady-state
systems, the mutual information and the learning rate behave qualitatively
distinctly, with the learning rate clearly now reflecting the rate at which the
downstream system must update its information in response to changes in the
upstream system. It is not clear whether this quantity is the most natural
measure for sensor performance, and, indeed, we provide an example in which
optimising the learning rate over a region of parameter space of the downstream
system yields an apparently sub-optimal sensor.


  Boolean matrix factorisation aims to decompose a binary data matrix into an
approximate Boolean product of two low rank, binary matrices: one containing
meaningful patterns, the other quantifying how the observations can be
expressed as a combination of these patterns. We introduce the OrMachine, a
probabilistic generative model for Boolean matrix factorisation and derive a
Metropolised Gibbs sampler that facilitates efficient parallel posterior
inference. On real world and simulated data, our method outperforms all
currently existing approaches for Boolean matrix factorisation and completion.
This is the first method to provide full posterior inference for Boolean Matrix
factorisation which is relevant in applications, e.g. for controlling false
positive rates in collaborative filtering and, crucially, improves the
interpretability of the inferred patterns. The proposed algorithm scales to
large datasets as we demonstrate by analysing single cell gene expression data
in 1.3 million mouse brain cells across 11 thousand genes on commodity
hardware.


  We analyse the asymptotic behaviour of integro-differential equations
modelling $N$ populations in interaction, all structured by different traits.
Interactions are modelled by non-local terms involving linear combinations of
the total number of individuals in each population. These models have already
been shown to be suitable for the modelling of drug resistance in cancer, and
they generalise the usual Lotka-Volterra ordinary differential equations. Our
aim is to give conditions under which there is persistence of all species.
Through the analysis of a Lyapunov function, our first main result gives a
simple and general condition on the matrix of interactions, together with a
convergence rate. The second main result establishes another type of condition
in the specific case of mutualistic interactions. When either of these
conditions is met, we describe which traits are asymptotically selected.


  Zika virus (ZIKV) exhibits unique transmission dynamics in that it is
concurrently spread by a mosquito vector and through sexual contact. We show
that this sexual component of ZIKV transmission induces novel processes on
networks through the highly asymmetric durations of infectiousness between
males and females -- it is estimated that males are infectious for periods up
to ten times longer than females -- leading to an asymmetric percolation
process on the network of sexual contacts. We exactly solve the properties of
this asymmetric percolation on random sexual contact networks and show that
this process exhibits two epidemic transitions corresponding to a
core-periphery structure. This structure is not present in the underlying
contact networks, which are not distinguishable from random networks, and
emerges because of the asymmetric percolation. We provide an exact analytical
description of this double transition and discuss the implications of our
results in the context of ZIKV epidemics. Most importantly, our study suggests
a bias in our current ZIKV surveillance as the community most at risk is also
one of the least likely to get tested.


  Whole genome duplication (WGD) is one of the most important events in the
molecular evolution of organisms. In fish species, a WGD is considered to have
occurred in the ancestral lineage of teleosts. Recent comprehensive ortholog
comparisons among teleost genomes have provided useful data and insights into
the fate of redundant genes generated by WGD. Based on these data, a
mathematical model is proposed to explain the evolutionary scenario of genes
after WGD. The model is parameterized taking into account an equilibrium
between i) rapid loss of either of the duplicate genes and ii) moderate
functional differentiation of each of duplicate genes, both of which are
followed by slow gene loss under purifying selection. This model predicts that,
in the teleost lineage, a maximum of about 3000 gene pairs may have
differentiated functionally during 90 million years after WGD. Thus, the
present study provides a possibility that the whole impact of WGD can be
quantitatively assessed according to the model parameters, before details of
genomic structural changes or functional differentiation are investigated. If
the equilibrium model is valid not only for teleosts but also for other
lineages that have undergone WGDs, correlations between the assessment indices
and evolutionarily significant events, such as the diversification of species
or the occurrence of novel phenotypes, could be tested and compared among those
lineages.


  Cancer is a disease of cellular regulation, often initiated by genetic
mutation within cells, and leading to a heterogeneous cell population within
tissues. In the competition for nutrients and growth space within the tumors
the phenotype of each cell determines its success. Selection in this process is
imposed by both the microenvironment (neighboring cells, extracellular matrix,
and diffusing substances), and the whole of the organism through for example
the blood supply. In this view, the development of tumor cells is in close
interaction with their increasingly changing environment: the more cells can
change, the more their environment will change. Furthermore, instabilities are
also introduced on the organism level: blood supply can be blocked by increased
tissue pressure or the tortuosity of the tumor-neovascular vessels. This
coupling between cell, microenvironment, and organism results in behavior that
is hard to predict. Here we introduce a cell-based computational model to study
the effect of blood flow obstruction on the micro-evolution of cells within a
cancerous tissue. We demonstrate that stages of tumor development emerge
naturally, without the need for sequential mutation of specific genes.
Secondly, we show that instabilities in blood supply can impact the overall
development of tumors and lead to the extinction of the dominant aggressive
phenotype, showing a clear distinction between the fitness at the cell level
and survival of the population. This provides new insights into potential side
effects of recent tumor vasculature renormalization approaches.


  Most of our current understanding of mechanisms of photosynthesis comes from
spectroscopy. However, classical definition of radio-antenna can be extended to
optical regime to discuss the function of light-harvesting antennae. Further to
our previously proposed model of a loop antenna we provide several more
physical explanations on considering the non-reciprocal properties of the light
harvesters of bacteria. We explained the function of the non-heme iron at the
reaction center, and presented reasons for each module of the light harvester
being composed of one carotenoid, two short $\alpha$-helical polypeptides and
three bacteriochlorophylls; we explained also the toroidal shape of the light
harvester, the upper bound of the characteristic length of the light harvester,
the functional role played by the long-lasting spectrometric signal observed,
and the photon anti-bunching observed. Based on these analyses, two mechanisms
might be used by radiation-durable bacteria, {\it Deinococcus radiodurans}; and
the non-reciprocity of an archaeon, {\it Haloquadratum walsbyi}, are analyzed.
The physical lessons involved are useful for designing artificial light
harvesters, optical sensors, wireless power chargers, passive super-Planckian
heat radiators, photocatalytic hydrogen generators, and radiation protective
cloaks. In particular it can predict what kind of particles should be used to
separate sunlight into a photovoltaically and thermally useful range to enhance
the efficiency of solar cells.


  When analyzing the genome, researchers have discovered that proteins bind to
DNA based on certain patterns of the DNA sequence known as "motifs". However,
it is difficult to manually construct motifs due to their complexity. Recently,
externally learned memory models have proven to be effective methods for
reasoning over inputs and supporting sets. In this work, we present memory
matching networks (MMN) for classifying DNA sequences as protein binding sites.
Our model learns a memory bank of encoded motifs, which are dynamic memory
modules, and then matches a new test sequence to each of the motifs to classify
the sequence as a binding or nonbinding site.


  The most profound change in the relationship between humans and their
environment was the introduction of agriculture and pastoralism. [....] For an
understanding of the expansion process, it appears appropriate to apply a
diffusive model. Broadly, these numerical modeling approaches can be catego-
rized in correlative, continuous and discrete. Common to all approaches is the
comparison to collections of radiocarbon data that show the apparent wave of
advance of the transition to farming. However, these data sets differ in entry
density and data quality. Often they disregard local and regional specifics and
research gaps, or dating uncertainties. Thus, most of these data bases may only
be used on a very general, broad scale. One of the pitfalls of using
irregularly spaced or irregularly documented radiocarbon data becomes evident
from the map generated by Fort (this volume, Chapter 16): while the general
east-west and south-north trends become evident, some areas appear as having
undergone anomalously early transitions to farming. This may be due to faulty
entries into the data base or regional problems with radiocarbon dating, if not
unnoticed or undocumented laboratory mistakes.


  The advent of rapid and inexpensive DNA sequencing has led to an explosion of
data waiting to be transformed into knowledge about genome organization and
function. Gene prediction is customarily the starting point for genome
analysis. This paper presents a bioinformatics study of the oil palm genome,
including comparative genomics analysis, database and tools development, and
mining of biological data for genes of interest. We have annotated 26,059 oil
palm genes integrated from two independent gene-prediction pipelines, Fgenesh++
and Seqping. This integrated annotation constitutes a significant improvement
in comparison to the preliminary annotation published in 2013. We conducted a
comprehensive analysis of intronless, resistance and fatty acid biosynthesis
genes, and demonstrated that the high quality of the current genome annotation.
3,658 intronless genes were identified in the oil palm genome, an important
resource for evolutionary study. Further analysis of the oil palm genes
revealed 210 candidate resistance genes involved in pathogen defense. Fatty
acids have diverse applications ranging from food to industrial feedstocks, and
we identified 42 key genes involved in fatty acid biosynthesis in oil palm.
These results provide an important resource for studies of plant genomes and a
theoretical foundation for marker-assisted breeding of oil palm and related
crops.


  Partially-observed Boolean dynamical systems (POBDS) are a general class of
nonlinear models with application in estimation and control of Boolean
processes based on noisy and incomplete measurements. The optimal minimum mean
square error (MMSE) algorithms for POBDS state estimation, namely, the Boolean
Kalman filter (BKF) and Boolean Kalman smoother (BKS), are intractable in the
case of large systems, due to computational and memory requirements. To address
this, we propose approximate MMSE filtering and smoothing algorithms based on
the auxiliary particle filter (APF) method from sequential Monte-Carlo theory.
These algorithms are used jointly with maximum-likelihood (ML) methods for
simultaneous state and parameter estimation in POBDS models. In the presence of
continuous parameters, ML estimation is performed using the
expectation-maximization (EM) algorithm; we develop for this purpose a special
smoother which reduces the computational complexity of the EM algorithm. The
resulting particle-based adaptive filter is applied to a POBDS model of Boolean
gene regulatory networks observed through noisy RNA-Seq time series data, and
performance is assessed through a series of numerical experiments using the
well-known cell cycle gene regulatory model.


  Efficient bacterial chromosome segregation typically requires the coordinated
action of a three-components, fueled by adenosine triphosphate machinery called
the partition complex. We present a phenomenological model accounting for the
dynamic activity of this system. The model is obtained by coupling simple
linear reaction-diffusion equations with a proteophoresis, or "volumetric"
chemophoresis, force field. This minimal description, in the sense of Occam's
Razor principle, captures most known experimental observations: dynamic
oscillations of complex components, complex separation and subsequent
symmetrical positioning. The predictions of our model are in phenomenological
agreement with and provide substantial insights into recent experiments. From a
non-linear physics view point, this system explores the active separation of
matter at micrometric scales with a dynamical instability between static
positioning and travelling wave regimes.


  Allosteric regulation is found across all domains of life, yet we still lack
simple, predictive theories that directly link the experimentally tunable
parameters of a system to its input-output response. To that end, we present a
general theory of allosteric transcriptional regulation using the
Monod-Wyman-Changeux model. We rigorously test this model using the ubiquitous
simple repression motif in bacteria by first predicting the behavior of strains
that span a large range of repressor copy numbers and DNA binding strengths and
then constructing and measuring their response. Our model not only accurately
captures the induction profiles of these strains but also enables us to derive
analytic expressions for key properties such as the dynamic range and
$[EC_{50}]$. Finally, we derive an expression for the free energy of allosteric
repressors which enables us to collapse our experimental data onto a single
master curve that captures the diverse phenomenology of the induction profiles.


  This paper is concerned with the problem of stochastic control of gene
regulatory networks (GRNs) observed indirectly through noisy measurements and
with uncertainty in the intervention inputs. The partial observability of the
gene states and uncertainty in the intervention process are accounted for by
modeling GRNs using the partially-observed Boolean dynamical system (POBDS)
signal model with noisy gene expression measurements. Obtaining the optimal
infinite-horizon control strategy for this problem is not attainable in
general, and we apply reinforcement learning and Gaussian process techniques to
find a near-optimal solution. The POBDS is first transformed to a
directly-observed Markov Decision Process in a continuous belief space, and the
Gaussian process is used for modeling the cost function over the belief and
intervention spaces. Reinforcement learning then is used to learn the cost
function from the available gene expression data. In addition, we employ
sparsification, which enables the control of large partially-observed GRNs. The
performance of the resulting algorithm is studied through a comprehensive set
of numerical experiments using synthetic gene expression data generated from a
melanoma gene regulatory network.


  Hamiltonian Monte Carlo (HMC) is an efficient and effective means of sampling
posterior distributions on Euclidean space, which has been extended to
manifolds with boundary. However, some applications require an extension to
more general spaces. For example, phylogenetic (evolutionary) trees are defined
in terms of both a discrete graph and associated continuous parameters;
although one can represent these aspects using a single connected space, this
rather complex space is not suitable for existing HMC algorithms. In this
paper, we develop Probabilistic Path HMC (PPHMC) as a first step to sampling
distributions on spaces with intricate combinatorial structure. We define PPHMC
on orthant complexes, show that the resulting Markov chain is ergodic, and
provide a promising implementation for the case of phylogenetic trees in
open-source software. We also show that a surrogate function to ease the
transition across a boundary on which the log-posterior has discontinuous
derivatives can greatly improve efficiency.


  Ultrametric trees are trees whose leaves lie at the same distance from the
root. They are used to model the genealogy of a population of particles
co-existing at the same point in time. We show how the boundary of an
ultrametric tree, like any compact ultrametric space, can be represented in a
simple way via the so-called comb metric. We display a variety of examples of
random combs and explain how they can be used in applications. In particular,
we review some old and recent results regarding the genetic structure of the
population when throwing neutral mutations on the skeleton of the tree.


  Stochastic exponential growth is observed in a variety of contexts, including
molecular autocatalysis, nuclear fission, population growth, inflation of the
universe, viral social media posts, and financial markets. Yet literature on
modeling the phenomenology of these stochastic dynamics has predominantly
focused on one model, Geometric Brownian Motion (GBM), which can be described
as the solution of a Langevin equation with linear drift and linear
multiplicative noise. Using recent experimental results on stochastic
exponential growth of individual bacterial cell sizes, we motivate the need for
a more general class of phenomenological models of stochastic exponential
growth, which are consistent with the observation that the mean-rescaled
distributions are approximately stationary at long times. We show that this
behavior is not consistent with GBM, instead it is consistent with power law
multiplicative noise with positive fractional powers. Therefore, we consider
this general class of phenomenological models for stochastic exponential
growth, provide analytical solutions, and identify the important dimensionless
combination of model parameters which determines the shape of the mean-rescaled
distribution. We also provide a prescription for robustly inferring model
parameters from experimentally observed stochastic growth trajectories.


  The key finding in the DNA double helix model is the specific pairing or
binding between nucleotides A-T and C-G, and the pairing rules are the molecule
basis of genetic code. Unfortunately, no such rules have been discovered for
proteins. Here we show that similar rules and intrinsic sequence patterns
between intra-protein binding peptide fragments do exist, and they can be
extracted using a deep learning algorithm. Multi-millions of binding and
non-binding peptide fragments from currently available protein X-ray structures
are classified with an accuracy of up to 93%. This discovery has the potential
in helping solve protein folding and protein-protein interaction problems, two
open and fundamental problems in molecular biology.


  Optimal subset selection is an important task that has numerous algorithms
designed for it and has many application areas. STPGA contains a special
genetic algorithm supplemented with a tabu memory property (that keeps track of
previously tried solutions and their fitness for a number of iterations), and
with a regression of the fitness of the solutions on their coding that is used
to form the ideal estimated solution (look ahead property) to search for
solutions of generic optimal subset selection problems. I have initially
developed the programs for the specific problem of selecting training
populations for genomic prediction or association problems, therefore I give
discussion of the theory behind optimal design of experiments to explain the
default optimization criteria in STPGA, and illustrate the use of the programs
in this endeavor. Nevertheless, I have picked a few other areas of application:
supervised and unsupervised variable selection based on kernel alignment,
supervised variable selection with design criteria, influential observation
identification for regression, solving mixed integer quadratic optimization
problems, balancing gains and inbreeding in a breeding population. Some of
these illustrations pertain new statistical approaches.


  Using state-of-the-art techniques combining imaging methods and
high-throughput genomic mapping tools leaded to the significant progress in
detailing chromosome architecture of various organisms. However, a gap still
remains between the rapidly growing structural data on the chromosome folding
and the large-scale genome organization. Could a part of information on the
chromosome folding be obtained directly from underlying genomic DNA sequences
abundantly stored in the databanks? To answer this question, we developed an
original discrete double Fourier transform (DDFT). DDFT serves for the
detection of large-scale genome regularities associated with domains/units at
the different levels of hierarchical chromosome folding. The method is
versatile and can be applied to both genomic DNA sequences and corresponding
physico-chemical parameters such as base-pairing free energy. The latter
characteristic is closely related to the replication and transcription and can
also be used for the assessment of temperature or supercoiling effects on the
chromosome folding. We tested the method on the genome of Escherichia coli K-12
and found good correspondence with the annotated domains/units established
experimentally. As a brief illustration of further abilities of DDFT, the study
of large-scale genome organization for bacteriophage PHIX174 and bacterium
Caulobacter crescentus was also added. The combined experimental, modeling, and
bioinformatic DDFT analysis should yield more complete knowledge on the
chromosome architecture and genome organization.


  Allosteric effects are often underlying the activity of proteins and
elucidating generic design aspects and functional principles which are unique
to allosteric phenomena represents a major challenge. Here an approach which
consists in the in silico design of synthetic structures which, as the
principal element of allostery, encode dynamical long-range coupling among two
sites is presented. The structures are represented by elastic networks, similar
to coarse-grained models of real proteins. A strategy of evolutionary
optimization was implemented to iteratively improve allosteric coupling. In the
designed structures allosteric interactions were analyzed in terms of strain
propagation and simple pathways which emerged during evolution were identified
as signatures through which long-range communication was established. Moreover,
robustness of allosteric performance with respect to mutations was
demonstrated. As it turned out, the designed prototype structures reveal
dynamical properties resembling those found in real allosteric proteins. Hence,
they may serve as toy models of complex allosteric systems, such as proteins.


  The idea that a genetically fixed behavior evolved from the once differential
learning ability of individuals that performed the behavior is known as the
Baldwin effect. A highly influential paper [Hinton G.E., Nowlan S.J., 1987. How
learning can guide evolution. Complex Syst. 1, 495-502] claimed that this
effect can be observed in silico, but here we argue that what was actually
shown is that the learning ability is easily selected for. Then we demonstrate
the Baldwin effect to happen in the in silico scenario by estimating the
probability and waiting times for the learned behavior to become innate.
Depending on parameter values, we find that learning can increase the chance of
fixation of the learned behavior by several orders of magnitude compared with
the non-learning situation.


  Efficient allocation of limited medical resources is crucial for controlling
epidemic spreading on networks. Based on the susceptible-infected-susceptible
model, we solve an optimization problem as how best to allocate the limited
resources so as to minimize the prevalence, providing that the curing rate of
each node is positively correlated to its medical resource. By quenched
mean-field theory and heterogeneous mean-field (HMF) theory, we prove that
epidemic outbreak will be suppressed to the greatest extent if the curing rate
of each node is directly proportional to its degree, under which the effective
infection rate $\lambda$ has a maximal threshold
$\lambda_c^{opt}=1/\left\langle k \right\rangle$ where $\left\langle k
\right\rangle$ is average degree of the underlying network. For weak infection
region ($\lambda\gtrsim\lambda_c^{opt}$), we combine a perturbation theory with
Lagrange multiplier method (LMM) to derive the analytical expression of optimal
allocation of the curing rates and the corresponding minimized prevalence. For
general infection region ($\lambda>\lambda_c^{opt}$), the high-dimensional
optimization problem is converted into numerically solving low-dimensional
nonlinear equations by the HMF theory and LMM. Counterintuitively, in the
strong infection region the low-degree nodes should be allocated more medical
resources than the high-degree nodes to minimize the prevalence. Finally, we
use simulated annealing to validate the theoretical results.


  Neutral evolution assumes that there are no selective forces distinguishing
different variants in a population. Despite this striking assumption, many
recent studies have sought to assess whether neutrality can provide a good
description of different episodes of cultural change. One approach has been to
test whether neutral predictions are consistent with observed progeny
distributions, recording the number of variants that have produced a given
number of new instances within a specified time interval: a classic example is
the distribution of baby names. Using an overlapping generations model we show
that these distributions consist of two phases: a power law phase with a
constant exponent of -3/2, followed by an exponential cut-off for variants with
very large numbers of progeny. Maximum likelihood estimations of the model
parameters provide a direct way to establish whether observed empirical
patterns are consistent with neutral evolution. We apply our approach to a
complete data set of baby names from Australia. Crucially we show that analyses
based on only the most popular variants, as is often the case in studies of
cultural evolution, can provide misleading evidence for underlying transmission
hypotheses. While neutrality provides a plausible description of progeny
distributions of abundant variants, rare variants deviate from neutrality.
Further, we develop a simulation framework that allows for the detection of
alternative cultural transmission processes. We show that anti-novelty bias is
able to replicate the complete progeny distribution of the Australian data set.


  "Three is a crowd" is an old proverb that applies as much to social
interactions, as it does to frustrated configurations in statistical physics
models. Accordingly, social relations within a triangle deserve special
attention. With this motivation, we explore the impact of topological
frustration on the evolutionary dynamics of the snowdrift game on a triangular
lattice. This topology provides an irreconcilable frustration, which prevents
anti-coordination of competing strategies that would be needed for an optimal
outcome of the game. By using different strategy updating protocols, we observe
complex spatial patterns in dependence on payoff values that are reminiscent to
a honeycomb-like organization, which helps to minimize the negative consequence
of the topological frustration. We relate the emergence of these patterns to
the microscopic dynamics of the evolutionary process, both by means of
mean-field approximations and Monte Carlo simulations. For comparison, we also
consider the same evolutionary dynamics on the square lattice, where of course
the topological frustration is absent. However, with the deletion of diagonal
links of the triangular lattice, we can gradually bridge the gap to the square
lattice. Interestingly, in this case the level of cooperation in the system is
a direct indicator of the level of topological frustration, thus providing a
method to determine frustration levels in an arbitrary interaction network.


  Simplified mechanistic models in ecology have been criticized for the fact
that a good fit to data does not imply the mechanism is true: pattern does not
equal process. In parallel, the maximum entropy principle (MaxEnt) has been
applied in ecology to make predictions constrained by just a handful of state
variables, like total abundance or species richness. But an outstanding
question remains: what principle tells us which state variables to constrain?
Here we attempt to solve both problems simultaneously, by translating a given
set of mechanisms into the state variables to be used in MaxEnt, and then using
this MaxEnt theory as a null model against which to compare mechanistic
predictions. In particular, we identify the sufficient statistics needed to
parametrize a given mechanistic model from data and use them as MaxEnt
constraints. Our approach isolates exactly what mechanism is telling us over
and above the state variables alone.


  As South and Central American countries prepare for increased birth defects
from Zika virus outbreaks and plan for mitigation strategies to minimize
ongoing and future outbreaks, understanding important characteristics of Zika
outbreaks and how they vary across regions is a challenging and important
problem. We developed a mathematical model for the 2015 Zika virus outbreak
dynamics in Colombia, El Salvador, and Suriname. We fit the model to publicly
available data provided by the Pan American Health Organization, using
Approximate Bayesian Computation to estimate parameter distributions and
provide uncertainty quantification. An important model input is the at-risk
susceptible population, which can vary with a number of factors including
climate, elevation, population density, and socio-economic status. We informed
this initial condition using the highest historically reported dengue incidence
modified by the probable dengue reporting rates in the chosen countries. The
model indicated that a country-level analysis was not appropriate for Colombia.
We then estimated the basic reproduction number, or the expected number of new
human infections arising from a single infected human, to range between 4 and 6
for El Salvador and Suriname with a median of 4.3 and 5.3, respectively. We
estimated the reporting rate to be around 16% in El Salvador and 18% in
Suriname with estimated total outbreak sizes of 73,395 and 21,647 people,
respectively. The uncertainty in parameter estimates highlights a need for
research and data collection that will better constrain parameter ranges.


  We study a generic one-dimensional model for an intracellular cargo driven by
N motor proteins against an external applied force. The model includes
motor-cargo and motor-motor interactions. The cargo motion is described by an
over-damped Langevin equation, while motor dynamics is specified by hopping
rates which follow a local detailed balance condition with respect to change in
energy per hopping event. Based on this model, we show that the stall force,
the mean external force corresponding to zero mean cargo velocity, is
completely independent of the details of the interactions and is, therefore,
always equal to the sum of the stall forces of the individual motors. This
exact result is arrived on the basis of a simple assumption: the (macroscopic)
state of stall of the cargo is analogous to a state of thermodynamic
equilibrium, and is characterized by vanishing net probability current between
any two microstates, with the latter specified by motor positions relative to
the cargo. The corresponding probability distribution of the microstates under
stall is also determined. These predictions are in complete agreement with
numerical simulations, carried out using specific forms of interaction
potentials.


  Oscillations in a stochastic dynamical system, whose deterministic
counterpart has a stable steady state, are a widely reported phenomenon.
Traditional methods of finding parameter regimes for stochastically-driven
resonances are, however, cumbersome for any but the smallest networks. In this
letter we show by example of the Brusselator how to use real root counting
algorithms and graph theoretic tools to efficiently determine the number of
resonant modes and parameter ranges for stochastic oscillations. We argue that
stochastic resonance is a network property by showing that resonant modes only
depend on the squared Jacobian matrix $J^2$ , unlike deterministic oscillations
which are determined by $J$. By using graph theoretic tools, analysis of
stochastic behaviour for larger networks is simplified and chemical reaction
networks with multiple resonant modes can be identified easily.


  Adaptation in a fluctuating environment is a process of fueling environmental
information to gain fitness. Living systems have gradually developed strategies
for adaptation from random and passive diversification of the phenotype to more
proactive decision making, in which environmental information is sensed and
exploited more actively and effectively. Understanding the fundamental relation
between fitness and information is therefore crucial to clarify the limits and
universal properties of adaptation. In this work, we elucidate the underlying
stochastic and information-thermodynamic structure in this process, by deriving
causal fluctuation relations (FRs) of fitness and information. Combined with a
duality between phenotypic and environmental dynamics, the FRs reveal the limit
of fitness gain, the relation of time reversibility with the achievability of
the limit, and the possibility and condition for gaining excess fitness due to
environmental fluctuation. The loss of fitness due to causal constraints and
the limited capacity of real organisms is shown to be the difference between
time-forward and time-backward path probabilities of phenotypic and
environmental dynamics. Furthermore, the FRs generalize the concept of
evolutionary stable state (ESS) for fluctuating environment by giving the
probability that the optimal strategy on average can be invaded by a suboptimal
one owing to rare environmental fluctuation. These results clarify the
information thermodynamic structures in adaptation and evolution.


  We quantify the amount of regulation required to control growth in living
cells by a Maximum Entropy approach to the space of underlying metabolic states
described by genome-scale models. Results obtained for E. coli and human cells
are consistent with experiments and point to different regulatory strategies by
which growth can be fostered or repressed. Moreover we explicitly connect the
`inverse temperature' that controls MaxEnt distributions to the growth
dynamics, showing that the initial size of a colony may be crucial in
determining how an exponentially growing population organizes the phenotypic
space.


  Today, the CRISPR (clustered regularly interspaced short palindromic repeats)
region within bacterial and archaeal genomes is known to encode an adaptive
immune system. We rely on previous results on the evolution of the CRISPR
arrays, which led to the ordered independent loss model, introduced by Kupczok
and Bollback (2013). When focusing on the spacers (between the repeats), new
elements enter a CRISPR array at rate $\theta$ at the leader end of the array,
while all spacers present are lost at rate $\rho$ along the phylogeny relating
the sample. Within this model, we compute the distribution of distances of
spacers which are present in all arrays in samples of size $n=2$ and $n=3$. We
use these results to estimate the loss rate $\rho$ from spacer array data.


  The evolution of the adaptive immune system is characterized by changes in
the relative abundances of the B- and T-cell clones that make up its
repertoires. To fully capture this evolution, we need to describe the complex
dynamics of the response to pathogenic and self-antigenic stimulations, as well
as the statistics of novel lymphocyte receptors introduced throughout life.
Recent experiments, ranging from high-throughput immune repertoire sequencing
to quantification of the response to specific antigens, can help us
characterize the effective dynamics of the immune response. Here we describe
mathematical models informed by experiments that lead to a picture of clonal
competition in a highly stochastic context. We discuss how different types of
competition, noise and selection shape the observed clone-size distributions,
and contrast them with predictions of a neutral theory of clonal evolution.
These mathematical models show that memory and effector immune repertoire
evolution is far from neutral, and is driven by the history of the pathogenic
environment, while naive repertoire dynamics are consistent with neutral theory
and competition in a fixed antigenic environment. Lastly, we investigate the
effect of long-term clonal selection on repertoire aging.


  Geographic structure can affect patterns of genetic differentiation and
speciation rates. In this article, we investigate the dynamics of genetic
distances in a geographically structured metapopulation. We model the
metapopulation as a weighted directed graph, with d vertices corresponding to d
subpopulations that evolve according to an individual based model. The dynamics
of the genetic distances is then controlled by two types of transitions
-mutation and migration events. We show that, under a rare mutation - rare
migration regime, intra subpopulation diversity can be neglected and our model
can be approximated by a population based model. We show that under a large
population - large number of loci limit, the genetic distance between two
subpopulations converges to a deterministic quantity that can asymptotically be
expressed in terms of the hitting time between two random walks in the
metapopulation graph. Our result shows that the genetic distance between two
subpopulations does not only depend on the direct migration rates between them
but on the whole metapopulation structure.


  Despite its radical assumption of ecological equivalence between species,
neutral biodiversity theory can often provide good fits to species abundance
distributions observed in nature. Major criticisms of neutral theory have
focused on interspecific differences, which are in conflict with ecological
equivalence. However, neutrality in nature is also broken by differences
between conspecific individuals at different life stages, which in many
communities may vastly exceed interspecific differences between individuals at
similar stages. These within-species asymmetries have not been fully explored
in species-neutral models, and it is not known whether demographic stage
structure affects macroecological patterns in neutral theory. Here we present a
two-stage neutral model where fecundity and mortality change as an individual
transitions from one stage to the other. We explore several qualitatively
different scenarios, and compare numerically obtained species abundance
distributions to the predictions of unstructured neutral theory. We find that
abundance distributions are generally robust to this kind of stage structure,
but significant departures from unstructured predictions occur if adults have
sufficiently low fecundity and mortality. In addition, we show that the
cumulative number of births per species, which is distributed as a power law
with a 3/2 exponent, is invariant even when the abundance distribution departs
from unstructured model predictions. Our findings potentially explain power
law-like abundance distributions in organisms with strong demographic
structure, such as eusocial insects and humans, and partially rehabilitate
species abundance distributions from past criticisms as to their inability to
distinguish between biological mechanisms.


  High hydrostatic pressure is commonly encountered in many environments, but
the effects of high pressure on eukaryotic cells have been understudied. To
understand the effects of hydrostatic pressure in the model eukaryote,
Saccharomyces cerevisiae, we have performed quantitative experiments of cell
division, cell morphology, and cell death under a wide range of pressures. We
developed an automated image analysis method for quantification of the yeast
budding index - a measure of cell cycle state - as well as a continuum model of
budding to investigate the effect of pressure on cell division and cell
morphology. We find that the budding index, the average cell size, and the
eccentricity - a measure of how much the cell morphology varies from the being
elliptical - of the cells decrease with increasing pressure. Furthermore, high
hydrostatic pressure led to the small but finite probability of cell death via
both apoptosis and necrosis. Our experiments suggest that decrease of budding
index arises from cellular arrest or death at the cell cycle checkpoints during
different stages of cell division.


  We introduce an improved version of RECKONER, an error corrector for Illumina
whole genome sequencing data. By modifying its workflow we reduce the
computation time even 10 times. We also propose a new method of determination
of $k$-mer length, the key parameter of $k$-spectrum-based family of
correctors. The correction algorithms are examined on huge data sets, i.e.,
human and maize genomes for both Illumina HiSeq and MiSeq instruments.


  We present *K-means clustering algorithm and source code by expanding
statistical clustering methods applied in https://ssrn.com/abstract=2802753 to
quantitative finance. *K-means is statistically deterministic without
specifying initial centers, etc. We apply *K-means to extracting cancer
signatures from genome data without using nonnegative matrix factorization
(NMF). *K-means' computational cost is a fraction of NMF's. Using 1,389
published samples for 14 cancer types, we find that 3 cancers (liver cancer,
lung cancer and renal cell carcinoma) stand out and do not have cluster-like
structures. Two clusters have especially high within-cluster correlations with
11 other cancers indicating common underlying structures. Our approach opens a
novel avenue for studying such structures. *K-means is universal and can be
applied in other fields. We discuss some potential applications in quantitative
finance.


  Many microbial systems are known to actively reshape their proteomes in
response to changes in growth conditions induced e.g. by nutritional stress or
antibiotics. Part of the re-allocation accounts for the fact that, as the
growth rate is limited by targeting specific metabolic activities, cells simply
respond by fine-tuning their proteome to invest more resources into the
limiting activity (i.e. by synthesizing more proteins devoted to it). However,
this is often accompanied by an overall re-organization of metabolism, aimed at
improving the growth yield under limitation by re-wiring resource through
different pathways. While both effects impact proteome composition, the latter
underlies a more complex systemic response to stress. By focusing on E. coli's
`acetate switch', we use mathematical modeling and a re-analysis of empirical
data to show that the transition from a predominantly fermentative to a
predominantly respirative metabolism in carbon-limited growth results from the
trade-off between maximizing the growth yield and minimizing its costs in terms
of required the proteome share. In particular, E. coli's metabolic phenotypes
appear to be Pareto-optimal for these objective functions over a broad range of
dilutions.


  The role of proton tunneling in biological catalysis is investigated here
within the frameworks of quantum information theory and thermodynamics. We
consider the quantum correlations generated through two hydrogen bonds between
a substrate and a prototypical enzyme that first catalyzes the tautomerization
of the substrate to move on to a subsequent catalysis, and discuss how the
enzyme can derive its catalytic potency from these correlations. In particular,
we show that classical changes induced in the binding site of the enzyme
spreads the quantum correlations among all of the four hydrogen-bonded atoms
thanks to the directionality of hydrogen bonds. If the enzyme rapidly returns
to its initial state after the binding stage, the substrate ends in a new
transition state corresponding to a quantum superposition. Open quantum system
dynamics can then naturally drive the reaction in the forward direction from
the major tautomeric form to the minor tautomeric form without needing any
additional catalytic activity. We find that in this scenario the enzyme lowers
the activation energy so much that there is no energy barrier left in the
tautomerization, even if the quantum correlations quickly decay.


  In this work it is shown that scale free tails in metabolic flux
distributions inferred from realistic large scale models can be simply an
artefact due to reactions involved in thermodynamically unfeasible cycles, that
are unbounded by physical constraints and would be able to perform work without
expenditure of free energy. After correcting for thermodynamics, the metabolic
space scales meaningfully with the physical limiting factors, acquiring in turn
a richer multimodal structure potentially leading to symmetry breaking while
optimizing for objective functions.


  In many natural situations one observes a local system with many competing
species which is coupled by weak immigration to a regional species pool. The
dynamics of such a system is dominated by its stable and uninvadable (SU)
states. When the competition matrix is random, the number of SUs depends on the
average value of its entries and the variance. Here we consider the problem in
the limit of weak competition and large variance. Using a yes/no interaction
model, we show that the number of SUs corresponds to the number of maximum
cliques in a network close to its fully connected limit. The number of SUs
grows exponentially with the number of species in this limit, unless the
network is completely asymmetric. In the asymmetric limit the number of SUs is
${\cal O} (1)$. Numerical simulations suggest that these results are valid for
models with continuous distribution of competition terms.


  In the course of evolution, proteins undergo important changes in their amino
acid sequences, while their three-dimensional folded structure and their
biological function remain remarkably conserved. Thanks to modern sequencing
techniques, sequence data accumulate at unprecedented pace. This provides large
sets of so-called homologous, i.e.~evolutionarily related protein sequences, to
which methods of inverse statistical physics can be applied. Using sequence
data as the basis for the inference of Boltzmann distributions from samples of
microscopic configurations or observables, it is possible to extract
information about evolutionary constraints and thus protein function and
structure. Here we give an overview over some biologically important questions,
and how statistical-mechanics inspired modeling approaches can help to answer
them. Finally, we discuss some open questions, which we expect to be addressed
over the next years.


  Proteins have evolved to perform diverse cellular functions, from serving as
reaction catalysts to coordinating cellular propagation and development.
Frequently, proteins do not exert their full potential as monomers but rather
undergo concerted interactions as either homo-oligomers or with other proteins
as hetero-oligomers. The experimental study of such protein complexes and
interactions has been arduous. Theoretical structure prediction methods are an
attractive alternative. Here, we investigate homo-oligomeric interfaces by
tracing residue coevolution via the global statistical Direct Coupling Analysis
(DCA). DCA can accurately infer spatial adjacencies between residues. These
adjacencies can be included as constraints in structure-prediction techniques
to predict high-resolution models. By taking advantage of the on-going
exponential growth of sequence databases, we go significantly beyond anecdotal
cases of a few protein families and apply DCA to a systematic large-scale study
of nearly 2000 PFAM protein families with sufficient sequence information and
structurally resolved homo-oligomeric interfaces. We find that large interfaces
are commonly identified by DCA. We further demonstrate that DCA can
differentiate between subfamilies of different binding modes within one large
PFAM family. Sequence derived contact information for the subfamilies proves
sufficient to assemble accurate structural models of the diverse
protein-oligomers. Thus, we provide an approach to investigate oligomerization
for arbitrary protein families leading to structural models complementary to
often difficult experimental methods. Combined with ever more abundant
sequential data, we anticipate that this study will be instrumental to allow
the structural description of many hetero-protein complexes in the future.


  Understanding how desertification takes place in different ecosystems is an
important step in attempting to forecast and prevent such transitions. Dryland
ecosystems often exhibit patchy vegetation, which has been shown to be an
important factor on the possible regime shifts that occur in arid regions in
several model studies. In particular, both gradual shifts that occur by front
propagation, and abrupt shifts where patches of vegetation vanish at once, are
a possibility in dryland ecosystems due to their emergent spatial
heterogeneity. However, recent theoretical work has suggested that the final
step of desertification - the transition from spotted vegetation to bare soil -
occurs only as an abrupt shift, but the generality of this result, and its
underlying origin, remain unclear. We investigate two models that detail the
dynamics of dryland vegetation using a markedly different functional structure,
and find that in both models the final step of desertification can only be
abrupt. Using a careful numerical analysis, we show that this behavior is
associated with the disappearance of confined spot-pattern domains as
stationary states, and identify the mathematical origin of this behavior. Our
findings show that a gradual desertification to bare soil due to a front
propagation process can not occur in these and similar models, and opens the
question of whether these dynamics can take place in nature.


  Microbial communities frequently communicate via quorum sensing (QS), where
cells produce, secrete, and respond to a threshold level of an autoinducer (AI)
molecule, thereby modulating gene expression. However, the biology of QS
remains incompletely understood in heterogeneous communities, where variant
bacterial strains possess distinct QS systems that produce chemically unique
AIs. AI molecules bind to 'cognate' receptors, but also to 'non-cognate'
receptors found in other strains, resulting in inter-strain crosstalk.
Understanding these interactions is a prerequisite for deciphering the
consequences of crosstalk in real ecosystems, where multiple AIs are regularly
present in the same environment. As a step towards this goal, we map crosstalk
in a heterogeneous community of variant QS strains onto an artificial neural
network model. This formulation allows us to systematically analyze how
crosstalk regulates the community's capacity for flexible decision making, as
quantified by the Boltzmann entropy of all QS gene expression states of the
system. In a mean-field limit of complete cross-inhibition between variant
strains, the model is exactly solvable, allowing for an analytical formula for
the number of variants that maximize capacity as a function of signal kinetics
and activation parameters. An analysis of previous experimental results on the
Staphylococcus aureus two-component Agr system indicates that the observed
combination of variant numbers, gene expression rates and threshold
concentrations lies near this critical regime of parameter space where capacity
peaks. The results are suggestive of a potential evolutionary driving force for
diversification in certain QS systems.


  The Ebola virus disease is a severe viral haemorrhagic fever syndrome caused
by Ebola virus. This disease is transmitted by direct contact with the body
fluids of an infected person and objects contaminated with virus or infected
animals, with a death rate close to 90% in humans. Recently, some mathematical
models have been presented to analyse the spread of the 2014 Ebola outbreak in
West Africa. In this paper, we introduce vaccination of the susceptible
population with the aim of controlling the spread of the disease and analyse
two optimal control problems related with the transmission of Ebola disease
with vaccination. Firstly, we consider the case where the total number of
available vaccines in a fixed period of time is limited. Secondly, we analyse
the situation where there is a limited supply of vaccines at each instant of
time for a fixed interval of time. The optimal control problems have been
solved analytically. Finally, we have performed a number of numerical
simulations in order to compare the models with vaccination and the model
without vaccination, which has recently been shown to fit the real data. Three
vaccination scenarios have been considered for our numerical simulations,
namely: unlimited supply of vaccines; limited total number of vaccines; and
limited supply of vaccines at each instant of time.


  Understanding the influence of structure of dispersal network on the species
persistence and modeling a much realistic species dispersal in nature are two
central issues in spatial ecology. A realistic dispersal structure which favors
the persistence of interacting ecological systems has been studied in [Holland
\& Hastings, Nature, 456:792--795 (2008)], where it is shown that a
randomization of the structure of dispersal network in a metapopulation model
of prey and predator increases the species persistence via clustering,
prolonged transient dynamics, and amplitudes of population fluctuations. In
this paper, by contrast, we show that a deterministic network topology in a
metapopulation can also favor asynchrony and prolonged transient dynamics if
species dispersal obeys a long-range interaction governed by a
distance-dependent power-law. To explore the effects of power-law coupling, we
take a realistic ecological model, namely the Rosenzweig-MacArthur model in
each patch (node) of the network of oscillators, and show that the coupled
system is driven from synchrony to asynchrony with an increase in the power-law
exponent. Moreover, to understand the relationship between species persistence
and variations in power-law exponent, we compute correlation coefficient to
characterize cluster formation, synchrony order parameter and median predator
amplitude. We further show that smaller metapopulations with less number of
patches are more vulnerable to extinction as compared to larger metapopulations
with higher number of patches. We believe that the present work improves our
understanding of the interconnection between the random network and
deterministic network in theoretical ecology.


  Which properties of metabolic networks can be derived solely from
stoichiometric information about the network's constituent reactions?
Predictive results have been obtained by Flux Balance Analysis (FBA), by
postulating that cells set metabolic fluxes within the allowed stoichiometry so
as to maximize their growth. Here, we generalize this framework to single cell
level using maximum entropy models from statistical physics. We define and
compute, for the core metabolism of Escherichia coli, a joint distribution over
all fluxes that yields the experimentally observed growth rate. This solution,
containing FBA as a limiting case, provides a better match to the measured
fluxes in the wild type and several mutants. We find that E. coli metabolism is
close to, but not at, the optimality assumed by FBA. Moreover, our model makes
a wide range of predictions: (i) on flux variability, its regulation, and flux
correlations across individual cells; (ii) on the relative importance of
stoichiometric constraints vs. growth rate optimization; (iii) on quantitative
scaling relations for singe-cell growth rate distributions. We validate these
scaling predictions using data from individual bacterial cells grown in a
microfluidic device at different sub-inhibitory antibiotic concentrations.
Under mild dynamical assumptions, fluctuation-response relations further
predict the autocorrelation timescale in growth data and growth rate adaptation
times following an environmental perturbation.


  Whether evolution can be predicted is a key question in evolutionary biology.
Here we set out to better understand the repeatability of evolution. We
explored experimentally the effect of mutation supply and the strength of
selective pressure on the repeatability of selection from standing genetic
variation. Different sizes of mutant libraries of an antibiotic resistance
gene, TEM-1 $\beta$-lactamase in Escherichia coli, were subjected to different
antibiotic concentrations. We determined whether populations went extinct or
survived, and sequenced the TEM gene of the surviving populations. The
distribution of mutations per allele in our mutant libraries- generated by
error-prone PCR- followed a Poisson distribution. Extinction patterns could be
explained by a simple stochastic model that assumed the sampling of beneficial
mutations was key for survival. In most surviving populations, alleles
containing at least one known large-effect beneficial mutation were present.
These genotype data also support a model which only invokes sampling effects to
describe the occurrence of alleles containing large-effect driver mutations.
Hence, evolution is largely predictable given cursory knowledge of mutational
fitness effects, the mutation rate and population size. There were no clear
trends in the repeatability of selected mutants when we considered all
mutations present. However, when only known large-effect mutations were
considered, the outcome of selection is less repeatable for large libraries, in
contrast to expectations. Furthermore, we show experimentally that alleles
carrying multiple mutations selected from large libraries confer higher
resistance levels relative to alleles with only a known large-effect mutation,
suggesting that the scarcity of high-resistance alleles carrying multiple
mutations may contribute to the decrease in repeatability at large library
sizes.


  While we once thought of cancer as single monolithic diseases affecting a
specific organ site, we now understand that there are many subtypes of cancer
defined by unique patterns of gene mutations. These gene mutational data, which
can be more reliably obtained than gene expression data, help to determine how
the subtypes develop, evolve, and respond to therapies. Different from dense
continuous-value gene expression data, which most existing cancer subtype
discovery algorithms use, somatic mutational data are extremely sparse and
heterogeneous, because there are less than 0.5\% mutated genes in discrete
value 1/0 out of 20,000 human protein-coding genes, and identical mutated genes
are rarely shared by cancer patients.
  Our focus is to search for cancer subtypes from extremely sparse and high
dimensional gene mutational data in discrete 1 and 0 values using unsupervised
learning. We propose a new network-based distance metric. We project cancer
patients' mutational profile into their gene network structure and measure the
distance between two patients using the similarity between genes and between
the gene vertexes of the patients in the network. Experimental results in
synthetic data and real-world data show that our approach outperforms the top
competitors in cancer subtype discovery. Furthermore, our approach can identify
cancer subtypes that cannot be detected by other clustering algorithms in real
cancer data.


  Many of the most important processes in cells take place on and across
membranes. With the rise of an impressive array of powerful quantitative
methods for characterizing these membranes, it is an opportune time to reflect
on the structure and function of membranes from the point of view of biological
numeracy. To that end, in this article, I review the quantitative parameters
that characterize the mechanical, electrical and transport properties of
membranes and carry out a number of corresponding order of magnitude estimates
that help us understand the values of those parameters.


  In this note, we characterize the embeddability of generic Kimura 3ST Markov
matrices in terms of their eigenvalues. As a consequence, we are able to
compute the volume of such matrices relative to the volume of all Markov
matrices within the model. We also provide examples showing that, in general,
mutation rates are not identifiable from substitution probabilities. These
examples also illustrate that symmetries between mutation probabilities do not
necessarily arise from symmetries between the corresponding mutation rates.


  We study the large population limit of a stochastic individual-based model
which describes the time evolution of a diploid hermaphroditic population
reproducing according to Mendelian rules. In [Neukirch, Bovier, 2016] it is
proved that sexual reproduction allows unfit alleles to survive in individuals
with mixed genotype much longer than they would in populations reproducing
asexually. In the present paper we prove that this indeed opens the possibility
that individuals with a pure genotype can reinvade in the population after the
appearance of further mutations. We thus expose a formal description of a
mechanism by which a recessive allele can re-emerge in a population. This can
be seen as a statement of genetic robustness exhibited by diploid populations
performing sexual reproduction.


  It is well-established that many physical properties of DNA at sufficiently
long length scales can be understood by means of simple polymer models. One of
the most widely used elasticity models for DNA is the twistable worm-like chain
(TWLC), which describes the double helix as a continuous elastic rod with
bending and torsional stiffness. An extension of the TWLC, which has recently
received some attention, is the model by Marko and Siggia, who introduced an
additional twist-bend coupling, expected to arise from the groove asymmetry. By
performing computer simulations of two available versions of oxDNA, a
coarse-grained model of nucleic acids, we investigate the microscopic origin of
twist-bend coupling. We show that this interaction is negligible in the oxDNA
version with symmetric grooves, while it appears in the oxDNA version with
asymmetric grooves. Our analysis is based on the calculation of the covariance
matrix of equilibrium deformations, from which the stiffness parameters are
obtained. The estimated twist-bend coupling coefficient from oxDNA simulations
is $G=30\pm1$~nm. The groove asymmetry induces a novel twist length scale and
an associated renormalized twist stiffness $\kappa_{\rm t} \approx 80$~nm,
which is different from the intrinsic torsional stiffness $C \approx 110$~nm.
This naturally explains the large variations on experimental estimates of the
intrinsic stiffness performed in the past.


  Here we report on a set of programs developed at the ZMBH Bio-Imaging
Facility for tracking real-life images of cellular processes. These programs
perform 1) automated tracking; 2) quantitative and comparative track analyses
of different images in different groups; 3) different interactive visualization
schemes; and 4) interactive realistic simulation of different cellular
processes for validation and optimal problem-specific adjustment of image
acquisition parameters (tradeoff between speed, resolution, and quality with
feedback from the very final results). The collection of programs is primarily
developed for the common bio-image analysis software ImageJ (as a single Java
Plugin). Some programs are also available in other languages (C++ and
Javascript) and may be run simply with a web-browser; even on a low-end Tablet
or Smartphone. The programs are available at
https://github.com/nurlicht/CellulyzerDemo


  The epigenome, i.e. the whole of chromatin modifications, is transferred from
mother to daughter cells during cell differentiation. When de novo chromatin
modifications (establishment or erasure of, respectively, new or pre-existing
DNA methylations and/or histone modifications) are made in a daughter cell,
however, it has a different epigenome than its mother cell. Although de novo
chromatin modifications are an important event that comprises elementary
processes of cell differentiation, its molecular mechanism remains poorly
understood. We argue in this Letter that a key to solving this problem lies in
understanding the role of long non-coding RNAs (lncRNAs)- a type of RNA that is
becoming increasingly prominent in epigenetic studies. Many studies show that
lncRNAs form ribonucleo-protein complexes in the nucleus and are involved in
chromatin modifications. However, chromatin-modifying enzymes lack the
information about genomic positions on which they act. It is known, on the
other hand, that a single-stranded RNA in general can bind to a double-stranded
DNA to form a triple helix. If each lncRNA forms a ribonucleo-protein complex
with chromatin-modifying enzymes on one hand and, at the same time, a triple
helix with a genomic region based on its specific nucleotide sequence on the
other hand, it can induce de novo chromatin modifications at specific sites.
Thus, the great variety of lncRNAs can be explained by the requirement for the
diversity of "genomic address codes" specific to their cognate genomic regions
where de novo chromatin modifications take place.


  Liquid-liquid phase separation of intrinsically disordered proteins (IDPs) is
a major undergirding factor in the regulated formation of membraneless
organelles in the cell. The phase behavior of an IDP is sensitive to its amino
acid sequence. Here we apply a recent random-phase-approximation polymer theory
to investigate how the tendency for multiple chains of a protein to phase
separate, as characterized by the critical temperature $T^*_{\rm cr}$, is
related to the protein's single-chain average radius of gyration $\langle
R_{\rm g} \rangle$. For a set of sequences containing different permutations of
an equal number of positively and negatively charged residues, we found a
striking correlation $T^*_{\rm cr}\sim \langle R_{\rm g} \rangle^{-\gamma}$
with $\gamma$ as large as $\sim 6.0$, indicating that electrostatic effects
have similarly significant impact on promoting single-chain conformational
compactness and phase separation. Moreover, $T^*_{\rm cr}\propto -{\rm SCD}$,
where SCD is a recently proposed "sequence charge decoration" parameter
determined solely by sequence information. Ramifications of our findings for
deciphering the sequence dependence of IDP phase separation are discussed.


  Gene expression is a noisy process and several mechanisms, both
transcriptional and posttranscriptional, can stabilize protein levels in cells.
Much work has focused on the role of miRNAs, showing in particular that
miRNA-mediated regulation can buffer expression noise for lowly expressed
genes. Here, using in silico simulations and mathematical modeling, we
demonstrate that miRNAs can exert a much broader influence on protein levels by
orchestrating competition-induced crosstalk between mRNAs. Most notably, we
find that miRNA-mediated cross-talk (i) can stabilize protein levels across the
full range of gene expression rates, and (ii) modifies the correlation pattern
of co-regulated interacting proteins, changing the sign of correlations from
negative to positive. The latter feature may constitute a potentially robust
signature of the existence of RNA crosstalk induced by endogenous competition
for miRNAs in standard cellular conditions.


  A key problem in modelling the evolution dynamics of infectious diseases is
the mathematical representation of the mechanism of transmission of the
contagion. Models with a finite number of subpopulations can be described via
systems of ordinary differential equations. When dealing with populations with
space structure the relevant quantities are spatial densities, whose evolution
in time requires nonlinear partial differential equations, which are known as
reaction-diffusion systems. Here we present an (historical) outline of
mathematical epidemiology, with a particular attention to the role of spatial
heterogeneity and dispersal in the population dynamics of infectious diseases.
Two specific examples are discussed, which have been the subject of intensive
research by the authors, i.e. man-environment-man epidemics, and malaria. In
addition to the epidemiological relevance of these epidemics all over the
world, their treatment requires a large amount of different sophisticate
mathematical methods, and has even posed new non trivial mathematical problems,
as one can realize from the list of references. One of the most relevant
problems posed by the authors, i.e. regional control, has been emphasized here:
the public health concern consists of eradicating the disease in the relevant
population, as fast as possible. On the other hand, very often the entire
domain of interest for the epidemic, is either unknown, or difficult to manage
for an affordable implementation of suitable environmental programmes. For
regional control instead it might be sufficient to implement such programmes
only in a given subregion conveniently chosen so to lead to an effective
(exponentially fast) eradication of the epidemic in the whole habitat; it is
evident that this practice may have an enormous importance in real cases with
respect to both financial and practical affordability.