  The high accuracy exhibited by biological information transcription processes
is due to kinetic proofreading, i.e., by a mechanism which reduces the error
rate of the information-handling process by driving it out of equilibrium. We
provide a consistent thermodynamic description of enzyme-assisted assembly
processes involving competing substrates, in a Master Equation framework. We
introduce and evaluate a measure of the efficiency based on rigorous
non-equilibrium inequalities. The performance of several proofreading models
are thus analyzed and the related time, dissipation and efficiency vs. error
trade-offs exhibited for different discrimination regimes. We finally introduce
and analyze in the same framework a simple model which takes into account
correlations between consecutive enzyme-assisted assembly steps. This work
highlights the relevance of the distinction between energetic and kinetic
discrimination regimes in enzyme-substrate interactions.


  Ecosystems are commonly conceptualized as networks of interacting species.
However, partitioning natural diversity of organisms into discrete units is
notoriously problematic, and mounting experimental evidence raises the
intriguing question whether this perspective is appropriate for the microbial
world. Here, an alternative formalism is proposed that does not require
postulating the existence of species as fundamental ecological variables, and
provides a naturally hierarchical description of community dynamics. This
formalism allows approaching the "species problem" from the opposite direction.
While the classical models treat a world of imperfectly clustered organism
types as a perturbation around well-clustered "species", the presented approach
allows gradually adding structure to a fully disordered background. The
relevance of this theoretical construct for describing highly diverse natural
ecosystems is discussed.


  Using high-throughput Tethered Particle Motion single molecule experiments,
the double-stranded DNA persistence length, $L_p$, is measured in solutions
with Na$^+$ and Mg$^{2+}$ ions of various ionic strengths, $I$. Several
theoretical equations for $L_p(I)$ are fitted to the experimental data, but no
decisive theory is found which fits all the $L_p$ values for the two ion
valencies. Properly extracted from the particle trajectory using simulations,
$L_p$ varies from 30~nm to 55~nm, and is compared to previous experimental
results. For the Na$^+$ only case, $L_p$ is an increasing concave function of
$I^{-1}$, well fitted by Manning's electrostatic stretching approach, but not
by classical Odjik-Skolnick-Fixman theories with or without counter-ion
condensation. With added Mg$^{2+}$ ions, $L_p$ shows a marked decrease at low
$I$, interpreted as an ion-ion correlation effect, with an almost linear law in
$I^{-1}$, fitted by a proposed variational approach.


  DNA condensation by multivalent cations plays a crucial role in genome
packaging in viruses and sperm heads, and has been extensively studied using
single-molecule experimental methods. In those experiments, the values of the
critical condensation forces have been used to estimate the amplitude of the
attractive DNA-DNA interactions. Here, to describe these experiments, we
developed an analytical model and a rigid body Langevin dynamics assay to
investigate the behavior of a polymer with self-interactions, in the presence
of a traction force applied at its extremities. We model self-interactions
using a pairwise attractive potential, thereby treating the counterions
implicitly. The analytical model allows to accurately predict the equilibrium
structures of toroidal and rod-like condensed structures, and the dependence of
the critical condensation force on the DNA length. We find that the critical
condensation force depends strongly on the length of the DNA, and finite-size
effects are important for molecules of length up to 10^5 {\mu}m. Our Langevin
dynamics simulations show that the force-extension behavior of the rod-like
structures is very different from the toroidal ones, so that their presence in
experiments should be easily detectable. In double-stranded DNA condensation
experiments, the signature of the presence of rod-like structures was not
unambiguously detected, suggesting that the polyamines used to condense DNA may
protect it from bending sharply as needed in the rod-like structures


  Complex biological systems have been successfully modeled by biochemical and
genetic interaction networks, typically gathered from high-throughput (HTP)
data. These networks can be used to infer functional relationships between
genes or proteins. Using the intuition that the topological role of a gene in a
network relates to its biological function, local or diffusion based
"guilt-by-association" and graph-theoretic methods have had success in
inferring gene functions. Here we seek to improve function prediction by
integrating diffusion-based methods with a novel dimensionality reduction
technique to overcome the incomplete and noisy nature of network data. In this
paper, we introduce diffusion component analysis (DCA), a framework that plugs
in a diffusion model and learns a low-dimensional vector representation of each
node to encode the topological properties of a network. As a proof of concept,
we demonstrate DCA's substantial improvement over state-of-the-art
diffusion-based approaches in predicting protein function from molecular
interaction networks. Moreover, our DCA framework can integrate multiple
networks from heterogeneous sources, consisting of genomic information,
biochemical experiments and other resources, to even further improve function
prediction. Yet another layer of performance gain is achieved by integrating
the DCA framework with support vector machines that take our node vector
representations as features. Overall, our DCA framework provides a novel
representation of nodes in a network that can be used as a plug-in architecture
to other machine learning algorithms to decipher topological properties of and
obtain novel insights into interactomes.


  Wild poliovirus type 3 (WPV3) has not been seen anywhere since the last case
of WPV3-associated paralysis in Nigeria in November 2012. At the time of
writing, the most recent case of wild poliovirus type 1 (WPV1) in Nigeria
occurred in July 2014, and WPV1 has not been seen in Africa since a case in
Somalia in August 2014. No cases associated with circulating vaccine-derived
type 2 poliovirus (cVDPV2) have been detected in Nigeria since November 2014.
Has WPV1 been eliminated from Africa? Has WPV3 been eradicated globally? Has
Nigeria interrupted cVDPV2 transmission? These questions are difficult because
polio surveillance is based on paralysis and paralysis only occurs in a small
fraction of infections.
  This report provides estimates for the probabilities of poliovirus
elimination in Nigeria given available data as of March 31, 2015. It is based
on a model of disease transmission that is built from historical polio
incidence rates and is designed to represent the uncertainties in transmission
dynamics and poliovirus detection that are fundamental to interpreting long
time periods without cases.
  The model estimates that, as of March 31, 2015, the probability of WPV1
elimination in Nigeria is 84%, and that if WPV1 has not been eliminated, a new
case will be detected with 99% probability by the end of 2015. The probability
of WPV3 elimination (and thus global eradication) is >99%. However, it is
unlikely that the ongoing transmission of cVDPV2 has been interrupted; the
probability of cVDPV2 elimination rises to 83% if no new cases are detected by
April 2016.
  Added July 10, 2015: On June 26, a cVDPV2 case was confirmed by the Global
Polio Laboratory Network. The date of paralysis was May 16. The case provides
new information about cVDPV2 prevalence that is useful for assessing the
accuracy of previous predictions and informing an updated forecast for the time
to cVDPV2 elimination.


  Bacterial vaginosis (BV) increases transmission of HIV, enhances the risk of
preterm labour, and its associated malodour impacts the quality of life for
many women. Clinical diagnosis primarily relies on microscopy to presumptively
detect a loss of lactobacilli and acquisition of anaerobes. This diagnostic
does not reflect the microbiota composition accurately as lactobacilli can
assume different morphotypes, and assigning BV associated morphotypes to
specific organisms is challenging. Using an untargeted metabolomics approach we
identify novel biomarkers for BV in a cohort of 131 Rwandan women, and
demonstrate that metabolic products in the vagina are strongly associated with
bacterial diversity. Metabolites associated with high diversity and clinical BV
include 2-hydroxyisovalerate and gamma-hydroxybutyrate (GHB), but not the
anaerobic end-product succinate. Low diversity, and high relative abundance of
lactobacilli, is characterized by lactate and amino acids. Biomarkers
associated with diversity and BV are independent of pregnancy status, and were
validated in a blinded replication cohort from Tanzania (n=45), in which we
predicted clinical BV with 91% accuracy. Correlations between the metabolome
and microbiota identified Gardnerella vaginalis as a putative producer of GHB,
and we demonstrate production by this species in vitro. This work provides a
deeper understanding of the relationship between the vaginal microbiota and
biomarkers of vaginal health and dysbiosis.


  The thermodynamical stability of a set of circular double helical molecules
is analyzed by path integral techniques. The minicircles differ only in
\textit{i)} the radius and \textit{ii)} the number of base pairs ($N$) arranged
along the molecule axis. Instead, the rise distance is kept constant. For any
molecule size, the computational method simulates a broad ensemble of possible
helicoidal configurations while the partition function is a sum over the path
trajectories describing the base pair fluctuational states. The stablest
helical repeat of every minicircle is determined by free energy minimization.
We find that, for molecules with $N$ larger than $100$, the helical repeat
grows linearly with the size and the twist number is constant. On the other
hand, by reducing the size below $100$ base pairs, the double helices sharply
unwind and the twist number drops to one for $N=\,20$. This is predicted as the
minimum size for the existence of helicoidal molecules in the closed form. The
helix unwinding appears as a strategy to release the bending stress associated
to the circularization of the molecules.


  Morphological transitions of phase separation associated with the asymmetry
of lipid composition were investigated using micrometer-sized vesicles of lipid
bilayers made from a lipid mixture. The complete macro-phase-separated
morphology undergoes a transition to a micro-phase-separation-like morphology
via a lorate morphology as a metastable state. The transition leads to the
emergence of monodisperse nanosized domains through repeated domain scission
events. Moreover, we have numerically confirmed the transitions using the
time-dependent Ginzburg-Landau model describing phase separation and the
bending elastic membrane, which is quantitatively consistent with experimental
results by fixing one free parameter. Our findings suggest that the local
spontaneous curvature due to the asymmetric composition plays an essential role
in the thermodynamic stabilization of micro-phase separation in lipid bilayers.


  For a genomically unstable cancer, a single tumour biopsy will often contain
a mixture of competing tumour clones. These tumour clones frequently differ
with respect to their genomic content (copy number of each gene) and structure
(order of genes on each chromosome). Modern bulk genome sequencing mixes the
signals of tumour clones and contaminating normal cells, complicating inference
of genomic content and structure. We propose a method to unmix tumour and
contaminating normal signals and jointly predict genomic structure and content
of each tumour clone. We use genome graphs to represent tumour clones, and
model the likelihood of the observed reads given clones and mixing proportions.
Our use of haplotype blocks allows us to accurately measure allele specific
read counts, and infer allele specific copy number for each clone. The proposed
method is a heuristic local search based on applying incremental, locally
optimal modifications of the genome graphs. Using simulated data, we show that
our method predicts copy counts and gene adjacencies with reasonable accuracy.


  When DNA molecules are heated they denature. This occurs locally so that
loops of molten single DNA strands form, connected by intact double-stranded
DNA pieces. The properties of this "melting" transition have been intensively
investigated. Recently there has been a surge of interest in this question,
caused by experiments determining the properties of partially bound DNA
confined to nanochannels. But how does such confinement affect the melting
transition? To answer this question we introduce, and solve a model predicting
how confinement affects the melting transition for a simple model system by
first disregarding the effect of self-avoidance. We find that the transition is
smoother for narrower channels. By means of Monte-Carlo simulations we then
show that a model incorporating self-avoidance shows qualitatively the same
behaviour and that the effect of confinement is stronger than in the ideal
case.


  Conventional studies of biomolecular behaviors rely largely on the
construction of kinetic schemes. Since the selection of these networks is not
unique, a concern is raised whether and under which conditions hierarchical
schemes can reveal the same experimentally measured fluctuating behaviors and
unique fluctuation related physical properties. To clarify these questions, we
introduce stochasticity into the traditional lumping analysis, generalize it
from rate equations to chemical master equations and stochastic differential
equations, and extract the fluctuation relations between kinetically and
thermodynamically equivalent networks under intrinsic and extrinsic noises. The
results provide a theoretical basis for the legitimate use of low-dimensional
models in the studies of macromolecular fluctuations and, more generally, for
exploring stochastic features in different levels of contracted networks in
chemical and biological kinetic systems.


  The behavior of long DNA molecules in a cell-sized confined space was
investigated. We prepared water-in-oil droplets covered by phospholipids, which
mimic the inner space of a cell, following the encapsulation of DNA molecules
with unfolded coil and folded globule conformations. Microscopic observation
revealed that the adsorption of coiled DNA onto the membrane surface depended
on the size of the vesicular space. Globular DNA showed a cell-size-dependent
unfolding transition after adsorption on the membrane. Furthermore, when DNA
interacted with a two-phase membrane surface, DNA selectively adsorbed on the
membrane phase, such as an ordered or disordered phase, depending on its
conformation. We discuss the mechanism of these trends by considering the free
energy of DNA together with a polyamine in the solution. The free energy of our
model was consistent with the present experimental data. The cooperative
interaction of DNA and polyamines with a membrane surface leads to the
size-dependent behavior of molecular systems in a small space. These findings
may contribute to a better understanding of the physical mechanism of molecular
events and reactions inside a cell.


  We study a mechanistic mathematical model of extinction and coexistence in a
generic hunter-prey ecosystem. The model represents typical scenarios of human
invasion and environmental change, characteristic of the late Pleistocene,
concomitant with the extinction of fauna in many regions of the world. As a
first approach we focus on a small trophic web of three species, including two
herbivores in asymmetric competition, in order to characterize the generic
behaviors. Specifically, we use a stochastic dynamical system, allowing the
study of the role of fluctuations and spatial correlations. We show that the
presence of hunters drives the superior herbivore to extinction even in
habitats that would allow coexistence, and even when the pressure of hunting is
lower than on the inferior one. The role of system size and fluctuating
populations is addressed, showing an ecological meltdown in small systems in
the presence of humans. The time to extinction as a function of the system
size, as calculated with the model, shows a good agreement with paleontological
data. Other findings show the intricate play of the anthropic and environmental
factors that may have caused the extinction of megafauna.


  In this paper we consider a strategy based on the values of the left
one-sided maximal Hardy-Littlewood operator to regulate the pulsating seeding
of young fish in the case where two or more competitive species (eg perch and
trout in lakes) interact; this strategy may occur as a natural from the point
of view of those who propose seeding and also allows controlling the magnitudes
resulting from competing species.


  Recent technological advances in cutting-edge ultrasensitive fluorescence
microscopy have allowed single-molecule imaging experiments in living cells
across all three domains of life to become commonplace. Single-molecule
live-cell data is typically obtained in a low signal-to-noise ratio (SNR)
regime sometimes only marginally in excess of 1, in which a combination of
detector shot noise, sub-optimal probe photophysics, native cell
autofluorescence and intrinsically underlying stochastic of molecules result in
highly noisy datasets for which underlying true molecular behaviour is
non-trivial to discern. The ability to elucidate real molecular phenomena is
essential in relating experimental single-molecule observations to both the
biological system under study as well as offering insight into the fine details
of the physical and chemical environments of the living cell. To confront this
problem of faithful signal extraction and analysis in a noise-dominated regime,
the needle in a haystack challenge, such experiments benefit enormously from a
suite of objective, automated, high-throughput analysis tools that can home in
on the underlying molecular signature and generate meaningful statistics across
a large population of individual cells and molecules. Here, I discuss the
development and application of several analytical methods applied to real case
studies, including objective methods of segmenting cellular images from light
microscopy data, tools to robustly localize and track single
fluorescently-labelled molecules, algorithms to objectively interpret molecular
mobility, analysis protocols to reliably estimate molecular stoichiometry and
turnover, and methods to objectively render distributions of molecular
parameters


  As proof-of-principle for generating superresolution structural information
from DNA we applied a method of localization microscopy utilizing photoblinking
comparing intercalating dye YOYO-1 against minor groove binding dye SYTO-13,
using a bespoke multicolor single-molecule fluorescence microscope. We used a
full-length ~49 kbp {\lambda} DNA construct possessing oligo inserts at either
terminus allowing conjugation of digoxigenin and biotin at opposite ends for
tethering to a glass coverslip surface and paramagnetic microsphere
respectively. We observed stochastic DNA-bound dye photoactivity consistent
with dye photoblinking as opposed to binding/unbinding events, evidenced
through both discrete simulations and continuum kinetics analysis. We analyzed
dye photoblinking images of immobilized DNA molecules using superresolution
reconstruction software from two existing packages, rainSTORM and QuickPALM,
and compared the results against our own novel home-written software called
ADEMS code. ADEMS code generated lateral localization precision values of 30-40
nm and 60-70 nm for YOYO-1 and SYTO-13 respectively at video-rate sampling,
similar to rainSTORM, running more slowly than rainSTORM and QuickPALM
algorithms but having a complementary capability over both in generating
automated centroid distribution and cluster analyses. Our imaging system allows
us to observe dynamic topological changes to single molecules of DNA in
real-time, such as rapid molecular snapping events. This will facilitate
visualization of fluorescently-labeled DNA molecules conjugated to a magnetic
bead in future experiments involving newly developed magneto-optical tweezers
combined with superresolution microscopy.


  DNA-interacting proteins have roles multiple processes, many operating as
molecular machines which undergo dynamic metastable transitions to bring about
their biological function. To fully understand this molecular heterogeneity,
DNA and the proteins that bind to it must ideally be interrogated at a single
molecule level in their native in vivo environments, in a time-resolved manner
fast to sample the molecular transitions across the free energy landscape.
Progress has been made over the past decade in utilising cutting-edge tools of
the physical sciences to address challenging biological questions concerning
the function and modes of action of several different proteins which bind to
DNA. These physiologically relevant assays are technically challenging, but can
be complemented by powerful and often more tractable in vitro experiments which
confer advantages of the chemical environment with enhanced detection
single-to-noise of molecular signatures and transition events. Here, we discuss
a range of techniques we have developed to monitor DNA-protein interactions in
vivo, in vitro and in silico. These include bespoke single-molecule
fluorescence microscopy techniques to elucidate the architecture and dynamics
of the bacterial replisome and the structural maintenance of bacterial
chromosomes, as well as new computational tools to extract single-molecule
molecular signatures from live cells to monitor stoichiometry, spatial
localization and mobility in living cells. We also discuss recent developments
from our lab made in vitro, complementing these in vivo studies, which combine
optical and magnetic tweezers to manipulate and image single molecules of DNA,
with and without bound protein, in a new superresolution fluorescence
microscope.


  In this paper, we use a finite-state continuous-time Markov chain with one
absorbing state to model an individual's lifetime. Under this model, the time
of death follows a phase-type distribution, and the transient states of the
Markov chain are known as phases. We then attempt to provide an answer to the
simple question "What is the conditional age distribution of the individual,
given its current phase"? We show that the answer depends on how we interpret
the question, and in particular, on the phase observation scheme under
consideration. We then apply our results to the computation of the age pyramid
for the endangered Chatham Island black robin Petroica traversi during years of
intensive conservation efforts in 1980-1989.


  In molecular mechanics, current generation potential energy functions provide
a reasonably good compromise between accuracy and effectiveness. This paper
firstly reviewed several most commonly used classical potential energy
functions and their optimization methods used for energy minimization. To
minimize a potential energy function, about 95\% efforts are spent on the
Lennard-Jones potential of van der Waals interactions; we also give a detailed
review on some effective computational optimization methods in the Cambridge
Cluster Database to solve the problem of Lennard-Jones clusters. From the
reviews, we found the hybrid idea of optimization methods is effective,
necessary and efficient for solving the potential energy minimization problem
and the Lennard-Jones clusters problem. An application to prion protein
structures is then done by the hybrid idea. We focus on the $\beta$2-$\alpha$2
loop of prion protein structures, and we found (i) the species that has the
clearly and highly ordered $\beta$2-$\alpha$2 loop usually owns a
3$_{10}$-helix in this loop, (ii) a "$\pi$-circle"
Y128--F175--Y218--Y163--F175--Y169--R164--Y128(--Y162) is around the
$\beta$2-$\alpha$2 loop.


  In order to transmit biochemical signals, biological regulatory systems
dissipate energy with concomitant entropy production. Additionally, signaling
often takes place in challenging environmental conditions. In a simple model
regulatory circuit given by an input and a delayed output, we explore the
trade-offs between information transmission and the system's energetic
efficiency. We determine the maximally informative network, given a fixed
amount of entropy production and delayed response, exploring both the case with
and without feedback. We find that feedback allows the circuit to overcome
energy constraints and transmit close to the maximum available information even
in the dissipationless limit. Negative feedback loops, characteristic of shock
responses, are optimal at high dissipation. Close to equilibrium positive
feedback loops, known for their stability, become more informative. Asking how
the signaling network should be constructed to best function in the worst
possible environment, rather than an optimally tuned one or in steady state, we
discover that at large dissipation the same universal motif is optimal in all
of these conditions.


  We present a structural data set of the 20 proteinogenic amino acids and
their amino-methylated and acetylated (capped) dipeptides. Different
protonation states of the backbone (uncharged and zwitterionic) were considered
for the amino acids as well as varied side chain protonation states.
Furthermore, we studied amino acids and dipeptides in complex with divalent
cations (Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+). The database covers the
conformational hierarchies of 280 systems in a wide relative energy range of up
to 4 eV (390 kJ/mol), summing up to an overall of 45,892 stationary points on
the respective potential-energy surfaces. All systems were calculated on equal
first-principles footing, applying density-functional theory in the generalized
gradient approximation corrected for long-range van der Waals interactions. We
show good agreement to available experimental data for gas-phase ion
affinities. Our curated data can be utilized, for example, for a wide
comparison across chemical space of the building blocks of life, for the
parametrization of protein force fields, and for the calculation of reference
spectra for biophysical applications.


  Beyond GWAS: alternatives to localize QTLs in farm animals. Two methods that
could be used for QTL mapping as alternatives to standard GWAS are presented.
The first relies on the differential frequency of runs of homozygosity (ROH) in
groups of animals (e.g. cases and controls), while the second stems from
resampling techniques used for the prediction of carriers of a mutation, and is
based on the frequency of inclusion of polymorphisms (SNP) in the predictive
model. ROH were applied to the detection of reproductive diseases in
Holstein-Friesian cattle, while resampling was applied to the detection of
carriers of the BH2 haplotype in Brown Swiss cattle. These alternative
approaches may complement GWAS analyses in localizing more accurately QTLs for
traits of interest in livestock.


  We analyze evolutionary dynamics on graphs, where the nodes represent
individuals of a population. The links of a node describe which other
individuals can be displaced by the offspring of the individual on that node.
Amplifiers of selection are graphs for which the fixation probability is
increased for advantageous mutants and decreased for disadvantageous mutants. A
few examples of such amplifiers have been developed, but so far it is unclear
how many such structures exist and how to construct them. Here, we show that
almost any undirected random graph is an amplifier of selection for Birth-death
updating, where an individual is selected to reproduce with probability
proportional to its fitness and one of its neighbors is replaced by that
offspring at random. If we instead focus on death-Birth updating, in which a
random individual is removed and its neighbors compete for the empty spot, then
the same ensemble of graphs consists of almost only suppressors of selection
for which the fixation probability is decreased for advantageous mutants and
increased for disadvantageous mutants. Thus, the impact of population structure
on evolutionary dynamics is a subtle issue that will depend on seemingly minor
details of the underlying evolutionary process.


  A method based on mapping a symbolic sequence into a set of patterns (strings
resulting from the sequence parsing) is proposed as a tool for the
reconstruction of ancestral sequences. The set union of patterns comprises all
the patterns present in the family of related proteins sequences of an extant
species. The set of most frequent patterns among protein sequences is selected
and concatenated. The resulting sequence of amino acids is supposed to be the
ancestral protein of the family. No sequences alignment and phylogenetic tree
of the species family are necessary. The method is used for inferring the
ancestral amino acid sequences of thirteen mitochondria-encoded protein
families of mammal species. Statistical distribution of the similarity between
extant and ancestral sequences exhibits some structures related to
environmental changes in the past.


  We discuss a model of protein conformations where the conformations are
combinations of short fragments from some small set. For these fragments we
consider a distribution of frequencies of occurrence of pairs (sequence of
amino acids, conformation), averaged over some balls in the spaces of sequences
and conformations. These frequencies can be estimated due to smallness of
epsilon-entropy of the set of conformations of protein fragments.
  We consider statistical potentials for protein fragments which describe the
mentioned frequencies of occurrence and discuss model of free energy of a
protein where the free energy is equal to a sum of statistical potentials of
the fragments.
  We also consider contribution of contacts of fragments to the energy of
protein conformation, and contribution from statistical potentials of some
hierarchical set of larger protein fragments. This set of fragments is
constructed using the distribution of frequencies of occurrence of short
fragments.
  We discuss applications of this model to problem of prediction of the native
conformation of a protein from its primary structure and to description of
dynamics of a protein. Modification of structural alignment taking into account
statistical potentials for protein fragments is considered and application to
threading procedure for proteins is discussed.


  Although not laying claim to being the inventor of the light microscope,
Antonj van Leeuwenhoek, (1632-1723) was arguably the first person to bring this
new technological wonder of the age properly to the attention of natural
scientists interested in the study of living things (people we might now term
biologists). He was a Dutch draper with no formal scientific training. From
using magnifying glasses to observe threads in cloth, he went on to develop
over 500 simple single lens microscopes with which he used to observe many
different biological samples. He communicated his finding to the Royal Society
in a series of letters including the one republished in this edition of Open
Biology. Our review here begins with the work of van Leeuwenhoek before
summarising the key developments over the last ca. 300 years which has seen the
light microscope evolve from a simple single lens device of van Leeuwenhoek's
day into an instrument capable of observing the dynamics of single biological
molecules inside living cells, and to tracking every cell nucleus in the
development of whole embryos and plants.


  Studies of fixation dynamics in Markov processes predominantly focus on the
mean time to absorption. This may be inadequate if the distribution is broad
and skewed. We compute the distribution of fixation times in one-step
birth-death processes with two absorbing states. These are expressed in terms
of the spectrum of the process, and we provide different representations as
forward-only processes in eigenspace. These allow efficient sampling of
fixation time distributions. As an application we study evolutionary game
dynamics, where invading mutants can reach fixation or go extinct. We also
highlight the median fixation time as a possible analog of mixing times in
systems with small mutation rates and no absorbing states, whereas the mean
fixation time has no such interpretation.


  Networks form the backbone of many complex systems, ranging from the Internet
to human societies. Accordingly, not only is the range of our interactions
limited and thus best described and modeled by networks, it is also a fact that
the networks that are an integral part of such models are often interdependent
or even interconnected. Networks of networks or multilayer networks are
therefore a more apt description of social systems. This colloquium is devoted
to evolutionary games on multilayer networks, and in particular to the
evolution of cooperation as one of the main pillars of modern human societies.
We first give an overview of the most significant conceptual differences
between single-layer and multilayer networks, and we provide basic definitions
and a classification of the most commonly used terms. Subsequently, we review
fascinating and counterintuitive evolutionary outcomes that emerge due to
different types of interdependencies between otherwise independent populations.
The focus is on coupling through the utilities of players, through the flow of
information, as well as through the popularity of different strategies on
different network layers. The colloquium highlights the importance of pattern
formation and collective behavior for the promotion of cooperation under
adverse conditions, as well as the synergies between network science and
evolutionary game theory.


  It is pointed out that the mystery of how biological systems measure their
lengths vanishes away if one premises that they have discovered a way to
generate linear waves analogous to compressional sound. These can be used to
detect length at either large or small scales using echo timing and fringe
counting. It is shown that suitable linear chemical potential waves can, in
fact, be manufactured by tuning to criticality conventional reaction-diffusion
with a small number substances. Min oscillations in E. coli are cited as
precedent resonant length measurement using chemical potential waves analogous
to laser detection. Mitotic structures in eucaryotes are identified as
candidates for such an effect at higher frequency. The engineering principle is
shown to be very general and functionally the same as that used by hearing
organs.
  PNAS Significance Statement: This paper invokes physical principles to
address the question of how living things might use reaction-diffusion to
measure out and regulate the many thousands of lengths required to make their
body parts and internal organs. It argues that two ideas have been missing. One
is that oscillation is necessary to achieve the necessary design stability and
plasticity. The other is that the system must be tuned to criticality to
stabilize the propagation velocity, thus enabling clocks to function as meter
sticks. The broader significance is twofold: First, a fundamental piece of the
machinery of life is probably invisible to present-day biochemical methods
because they are too slow. Second, the simplicity of growth and form identified
a century ago by D'Arcy Thompson is probably a symptom of biological
engineering strategies, not primitive law.


  We study an alternative approach to model the dynamical behaviors of
biological feedback loop, that is, a type-dependent spin system, this class of
stochastic models was introduced by Fern\'andez et. al (2009), and are useful
since take account to inherent variability of gene expression. We analyze a
non-symmetric feedback module being an extension for the repressilator, the
first synthetic biological oscillator, invented by Elowitz and Leibler (2000).
We consider a mean-field dynamics for a type-dependent Ising model, and then
study the empirical-magnetization vector representing concentration of
molecules. We apply a convergence result from stochastic jump processes to
deterministic trajectories and present a bifurcation analysis for the
associated dynamical system. We show that non-symmetric module under study can
exhibit very rich behaviours, including the empirical oscillations described by
repressilator.


  To characterize natural selection, various analytical methods for detecting
candidate genomic regions have been developed. We propose to perform
genome-wide scans of natural selection using principal component analysis. We
show that the common Fst index of genetic differentiation between populations
can be viewed as a proportion of variance explained by the principal
components. Considering the correlations between genetic variants and each
principal component provides a conceptual framework to detect genetic variants
involved in local adaptation without any prior definition of populations. To
validate the PCA-based approach, we consider the 1000 Genomes data (phase 1)
after removal of recently admixed individuals resulting in 850 individuals
coming from Africa, Asia, and Europe. The number of genetic variants is of the
order of 36 millions obtained with a low-coverage sequencing depth (3X). The
correlations between genetic variation and each principal component provide
well-known targets for positive selection (EDAR, SLC24A5, SLC45A2, DARC), and
also new candidate genes (APPBPP2, TP1A1, RTTN, KCNMA, MYO5C) and non-coding
RNAs. In addition to identifying genes involved in biological adaptation, we
identify two biological pathways involved in polygenic adaptation that are
related to the innate immune system (beta defensins) and to lipid metabolism
(fatty acid omega oxidation). An additional analysis of European data shows
that a genome scan based on PCA retrieves classical examples of local
adaptation even when there are no well-defined populations. PCA-based
statistics, implemented in the PCAdapt R package and the PCAdapt open-source
software, retrieve well-known signals of human adaptation, which is encouraging
for future whole-genome sequencing project, especially when defining
populations is difficult.


  We introduce the graph parameter readability and study it as a function of
the number of vertices in a graph. Given a digraph D, an injective overlap
labeling assigns a unique string to each vertex such that there is an arc from
x to y if and only if x properly overlaps y. The readability of D is the
minimum string length for which an injective overlap labeling exists. In
applications that utilize overlap digraphs (e.g., in bioinformatics),
readability reflects the length of the strings from which the overlap digraph
is constructed. We study the asymptotic behaviour of readability by casting it
in purely graph theoretic terms (without any reference to strings). We prove
upper and lower bounds on readability for certain graph families and general
graphs


  In this letter, a generalization of pairwise models to non-Markovian
epidemics on networks is presented. For the case of infectious periods of fixed
length, the resulting pairwise model is a system of delay differential
equations (DDEs), which shows excellent agreement with results based on
explicit stochastic simulations of non-Markovian epidemics on networks.
Furthermore, we analytically compute a new R0-like threshold quantity and an
implicit analytical relation between this and the final epidemic size. In
addition we show that the pairwise model and the analytic calculations can be
generalized in terms of integro-differential equations to any distribution of
the infectious period, and we illustrate this by presenting a closed form
expression for the final epidemic size. By showing the rigorous mathematical
link between non-Markovian network epidemics and pairwise DDEs, we provide the
framework for a deeper and more rigorous understanding of the impact of
non-Markovian dynamics with explicit results for final epidemic size and
threshold quantities.


  Cooperation is a very common, yet not fully-understood phenomenon in natural
and human systems. The introduction of a network within the population is known
to affect the outcome of cooperative dynamics, allowing for the survival of
cooperation in adverse scenarios. Recently, the introduction of multiplex
networks has yet again modified the expectations for the outcome of the
Prisoner's Dilemma game, compared to the monoplex case. However, much remains
unstudied regarding other social dilemmas on multiplex, as well as the
unexplored microscopic underpinnings of it. In this paper, we systematically
study the evolution of cooperation in all four games in the $T-S$ plane on
multiplex. More importantly, we find some remarkable and previously unknown
features in the microscopic organization of the strategies, that are
responsible for the important differences between cooperative dynamics in
monoplex and multiplex. Specifically, we find that in the stationary state,
there are individuals that play the same strategy in all layers (coherent), and
others that don't (incoherent). This second group of players is responsible for
the surprising fact of a non full-cooperation in the Harmony Game on multiplex,
never observed before, as well as a higher-than-expected cooperation rates in
some regions of the other three social dilemmas.


  The Triander program is an interactive software package for nucleotide
sequence visualization. The program was developed using the freeware Pascal RAD
IDE Lazarus, and its source code and binaries compiled for Windows are freely
accessible at http://www.icbge.org.ua/eng/Triander . Triander can produce four
types of plots. It is possible to build three DNA walks independently for each
nucleotide position in triplets. The use of nucleotide vectors with unequal
modulus leads to a significant reduction in the visual information lost in DNA
walks. The program can be used in the investigation of the fine structure of
sequences, to find standard patterns in them and to locate nontrivial regions
for further detailed analysis.


  Many fields use search algorithms, which automatically explore a search space
to find high-performing solutions: chemists search through the space of
molecules to discover new drugs; engineers search for stronger, cheaper, safer
designs, scientists search for models that best explain data, etc. The goal of
search algorithms has traditionally been to return the single
highest-performing solution in a search space. Here we describe a new,
fundamentally different type of algorithm that is more useful because it
provides a holistic view of how high-performing solutions are distributed
throughout a search space. It creates a map of high-performing solutions at
each point in a space defined by dimensions of variation that a user gets to
choose. This Multi-dimensional Archive of Phenotypic Elites (MAP-Elites)
algorithm illuminates search spaces, allowing researchers to understand how
interesting attributes of solutions combine to affect performance, either
positively or, equally of interest, negatively. For example, a drug company may
wish to understand how performance changes as the size of molecules and their
cost-to-produce vary. MAP-Elites produces a large diversity of high-performing,
yet qualitatively different solutions, which can be more helpful than a single,
high-performing solution. Interestingly, because MAP-Elites explores more of
the search space, it also tends to find a better overall solution than
state-of-the-art search algorithms. We demonstrate the benefits of this new
algorithm in three different problem domains ranging from producing modular
neural networks to designing simulated and real soft robots. Because MAP-
Elites (1) illuminates the relationship between performance and dimensions of
interest in solutions, (2) returns a set of high-performing, yet diverse
solutions, and (3) improves finding a single, best solution, it will advance
science and engineering.


  If a deuterated molecule containing strong intramolecular hydrogen bonds is
placed in a hydrogenated solvent it may preferentially exchange deuterium for
hydrogen. This preference is due to the difference between the vibrational
zero-point energy for hydrogen and deuterium. It is found that the associated
fractionation factor $\Phi$ is correlated with the strength of the
intramolecular hydrogen bonds. This correlation has been used to determine the
length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes
and has been argued to support the existence of short low-barrier H-bonds.
Starting with a potential energy surface based on a simple diabatic state model
for H-bonds we calculate $\Phi$ as a function of the proton donor-acceptor
distance $R$. For numerical results, we use a parameterization of the model for
symmetric O-H.... O bonds. We consider the relative contributions of the O-H
stretch vibration, O-H bend vibrations (both in plane and out of plane),
tunnelling splitting effects at finite temperature, and the secondary geometric
isotope effect. We compare our total $\Phi$ as a function of $R$ with NMR
experimental results for enzymes, and in particular with an empirical
parametrisation $\Phi(R)$, used previously to determine bond lengths.


  The risk-sensitive foraging theory formulated in terms of the (daily) energy
budget rule has been influential in behavioural ecology as well as other
disciplines. Predicting risk-aversion on positive budgets and risk-proneness on
negative budgets, however, the budget rule has recently been challenged both
empirically and theoretically. In this paper, we critically review these
challenges as well as the original derivation of the budget rule and propose a
`gradual' budget rule, which is normatively derived from a gradual nature of
risk sensitivity and encompasses the conventional budget rule as a special
case. The gradual budget rule shows that the conventional budget rule holds
when the expected reserve is close enough to a threshold for overnight
survival, selection pressure being significant. The gradual view also reveals
that the conventional budget rule does not need to hold when the expected
reserve is not close enough to the threshold, selection pressure being
insignificant. The proposed gradual budget rule better fits the empirical
findings including those that used to challenge the conventional budget rule.


  Recent studies suggest that PML (Promyelocytic Leukemia) suppresses p53
protein degradation by inhibiting MDM2 protein in the nucleus [1], and
regulates a number of biological functions [2]. We modeled a PML-MDM2-p53
regulatory network by integrating p53-MDM2 autoregulatory model [3] and related
PML pathways. We found that p53 dynamics switched at various dynamical states
induced by PML which corresponds to different cellular states. Our results show
clear transitions among the stabilized, damped and sustain states under
different stress conditions induced by PML in the system. These states in p53
dynamics is the signature of existence of various cellular states and the phase
transition like behaviour of these states involve various biological
significance.


  The amount of completely sequenced chloroplast genomes increases rapidly
every day, leading to the possibility to build large scale phylogenetic trees
of plant species. Considering a subset of close plant species defined according
to their chloroplasts, the phylogenetic tree that can be inferred by their core
genes is not necessarily well supported, due to the possible occurrence of
"problematic" genes (i.e., homoplasy, incomplete lineage sorting, horizontal
gene transfers, etc.) which may blur phylogenetic signal. However, a
trustworthy phylogenetic tree can still be obtained if the number of
problematic genes is low, the problem being to determine the largest subset of
core genes that produces the best supported tree. To discard problematic genes
and due to the overwhelming number of possible combinations, we propose an
hybrid approach that embeds both genetic algorithms and statistical tests.
Given a set of organisms, the result is a pipeline of many stages for the
production of well supported phylogenetic trees. The proposal has been applied
to different cases of plant families, leading to encouraging results for these
families.


  In this paper, dynamical systems theory and bifurcation theory are applied to
investi- gate the rich dynamical behaviours observed in three simple disease
models. The 2- and 3-dimensional models we investigate have arisen in previous
investigations of epidemiol- ogy, in-host disease, and autoimmunity. These
closely related models display interesting dynamical behaviors including
bistability, recurrence, and regular oscillations, each of which has possible
clinical or public health implications. In this contribution we elucidate the
key role of backward bifurcation in the parameter regimes leading to the
behaviors of interest. We demonstrate that backward bifurcation facilitates the
appearance of Hopf bifurcations, and the varied dynamical behaviors are then
determined by the properties of the Hopf bifurcation(s), including their
location and direction. A Maple program devel- oped earlier is implemented to
determine the stability of limit cycles bifurcating from the Hopf bifurcation.
Numerical simulations are presented to illustrate phenomena of interest such as
bistability, recurrence and oscillation. We also discuss the physical
motivations for the models and the clinical implications of the resulting
dynamics.


  We consider the reconstruction of a phylogeny from multiple genes under the
multispecies coalescent. We establish a connection with the sparse signal
detection problem, where one seeks to distinguish between a distribution and a
mixture of the distribution and a sparse signal. Using this connection, we
derive an information-theoretic trade-off between the number of genes, $m$,
needed for an accurate reconstruction and the sequence length, $k$, of the
genes. Specifically, we show that to detect a branch of length $f$, one needs
$m = \Theta(1/[f^{2} \sqrt{k}])$.


  One of the fundamental problems in biology concerns the method by which a
cluster of organisms can regulate the proportion of individuals that perform
various roles or modes as if each individual knows a whole situation without a
leader. A specific ratio exists in various species at multiple levels from the
process of cell differentiation in multicellular organisms to the situation of
social dilemma in a group of human beings. This study found a common basis of
regulating a collective behavior which is realized by a series of local
contacts between individuals. The most essential behavior of individuals in
this theory is to change its internal mode through sharing information in
contact with others. Our numerical simulations with cellular automata model
realize to regulate the ratio of population of individuals who has either two
kinds of modes. From the theoretical analysis and numerical calculations, we
found that asymmetric properties in local contacts, are essential for adaptive
regulation in response to the global information such a group size and whole
density. Furthermore, a discrete system is crucial in no-leader groups to
realize the flexible regulation, and the critical condition which eliminates
overlap with one another (excluded volume effect) also affects the resulting
proportion in high density. The foremost advantage of this strategy is that no
global information is required for each individual, and only a couple of mode
switching can achieve the whole proportion regulation. The simple mechanism say
that proportion regulation in well-organized groups in nature can be realized
through and limited to local contacts, and has a potential to solve various
phenomena that microscopic individuals behaviors connect to macroscopic orderly
behaviors.


  Human papillomavirus, or HPV, is a sexually transmittable virus infection,
which is necessary risk factor for developing cervical cancer, first most
common type of cancer in working age women in Moldova. We observe both
behavioral change (sexuality increase) and demographical change (population
ageing). We used data since 1998 (Moldovan peace treaty) to adjust model
parameter and we project till around 2030 (for vaccination till 2050).
According to provided information, interdisciplinary model was proposed. It iss
set of deterministic differential equations. Stochasticity was introduced in
sexual partner change rates. The model has aggregated the most important paths
of infection, cancer development and prevention scenarios (more than 100
equations and 200 parameters). Moldovan cervical cancer perspective looks much
better, than in central western Europe countries, because of relatively young
society. In our setup, obligatory vaccination seems to not be so crucial (for
none of realistic scenarios increase of cancer cases is possible) for public
health, as in most countries in European Union. However, screening practice
could be verified in terms of efficiency, when cost benefit calculation would
be done. We propose more optimal screening guidelines (with prevention cost 5
-10k EUR per QALY), which could provide saving perspective in 10-15 year in
range 150-300k EUR yearly. Targeted vaccination could be also consider, because
costs are similar to high frequencies screening schema with the same cancer
cases projection. However, some positive side effects of vaccination as
reduction of pathogen circulation in society, will cause decrease of other
pathologies related to HPV like genital warts and other cancer.


  Protein side-chain packing is a critical component in obtaining the 3D
coordinates of a structure and drug discovery. Single-domain protein side-chain
packing has been thoroughly studied. A major challenge in generalizing these
methods to protein complexes is that they, unlike monomers, often have very
large treewidth, and thus algorithms such as TreePack cannot be directly
applied. To address this issue, SCWRL4 treats the complex effectively as a
monomer, heuristically excluding weak interactions to decrease treewidth; as a
result, SCWRL4 generates poor packings on protein interfaces. To date, few
side-chain packing methods exist that are specifically designed for protein
complexes. In this paper, we introduce a method, iTreePack, which solves the
side-chain packing problem for complexes by using a novel combination of dual
decomposition and tree decomposition. In particular, iTreePack overcomes the
problem of large treewidth by decomposing a protein complex into smaller
subgraphs and novelly reformulating the complex side-chain packing problem as a
dual relaxation problem; this allows us to solve the side-chain packing of each
small subgraph separately using tree-decomposition. A projected subgradient
algorithm is applied to enforcing the consistency among the side-chain packings
of all the small subgraphs. Computational results demonstrate that our
iTreePack program outperforms SCWRL4 on protein complexes. In particular,
iTreePack places side-chain atoms much more accurately on very large complexes,
which constitute a significant portion of protein-protein interactions.
Moreover, the advantage of iTreePack over SCWRL4 increases with respect to the
treewidth of a complex. Even for monomeric proteins, iTreePack is much more
efficient than SCWRL and slightly more accurate.


  The early steps of photosynthesis involve the photo-excitation of reaction
centres (RCs) and light-harvesting (LH) units. Here, we show that the
--historically overlooked-- excitonic delocalisation across RC and LH pigments
results in a redistribution of dipole strengths that benefits the absorption
cross section of the optical bands associated with the RC of several species.
While we prove that this redistribution is robust to the microscopic details of
the dephasing between these units in the purple bacterium Rhodospirillum
rubrum, we are able to show that the redistribution witnesses a more fragile,
but persistent, coherent population dynamics which directs excitations from the
LH towards the RC units under incoherent illumination and physiological
conditions. Stochastic optimisation allows us to delineate clear guidelines and
develop simple analytic expressions, in order to achieve directed coherent
population dynamics in artificial nano-structures.


  Genetic variation in the tomato clade was explored by sequencing a selection
of 84 tomato accessions and related wild species representative for the
Lycopersicon, Arcanum, Eriopersicon, and Neolycopersicon groups. We present a
reconstruction of three new reference genomes in support of our comparative
genome analyses. Sequence diversity in commercial breeding lines appears
extremely low, indicating the dramatic genetic erosion of crop tomatoes. This
is reflected by the SNP count in wild species which can exceed 10 million i.e.
20 fold higher than in crop accessions. Comparative sequence alignment reveals
group, species, and accession specific polymorphisms, which explain
characteristic fruit traits and growth habits in tomato accessions. Using gene
models from the annotated Heinz reference genome, we observe a bias in dN/dS
ratio in fruit and growth diversification genes compared to a random set of
genes, which probably is the result of a positive selection. We detected highly
divergent segments in wild S. lycopersicum species, and footprints of
introgressions in crop accessions originating from a common donor accession.
Phylogenetic relationships of fruit diversification and growth specific genes
from crop accessions show incomplete resolution and are dependent on the
introgression donor. In contrast, whole genome SNP information has sufficient
power to resolve the phylogenetic placement of each accession in the four main
groups in the Lycopersicon clade using Maximum Likelihood analyses.
Phylogenetic relationships appear correlated with habitat and mating type and
point to the occurrence of geographical races within these groups and thus are
of practical importance for introgressive hybridization breeding. Our study
illustrates the need for multiple reference genomes in support of tomato
comparative genomics and Solanum genome evolution studies.


  Breeding by introgressive hybridization is a pivotal strategy to broaden the
genetic basis of crops. Usually, the desired traits are monitored in
consecutive crossing generations by marker-assisted selection, but their
analyses fail in chromosome regions where crossover recombinants are rare or
not viable. Here, we present the Introgression Browser (IBROWSER), a novel
bioinformatics tool aimed at visualizing introgressions at nucleotide or SNP
accuracy. The software selects homozygous SNPs from Variant Call Format (VCF)
information and filters out heterozygous SNPs, Multi-Nucleotide Polymorphisms
(MNPs) and insertion-deletions (InDels). For data analysis IBROWSER makes use
of sliding windows, but if needed it can generate any desired fragmentation
pattern through General Feature Format (GFF) information. In an example of
tomato (Solanum lycopersicum) accessions we visualize SNP patterns and
elucidate both position and boundaries of the introgressions. We also show that
our tool is capable of identifying alien DNA in a panel of the closely related
S. pimpinellifolium by examining phylogenetic relationships of the introgressed
segments in tomato. In a third example, we demonstrate the power of the
IBROWSER in a panel of 600 Arabidopsis accessions, detecting the boundaries of
a SNP-free region around a polymorphic 1.17 Mbp inverted segment on the short
arm of chromosome 4. The architecture and functionality of IBROWSER makes the
software appropriate for a broad set of analyses including SNP mining, genome
structure analysis, and pedigree analysis. Its functionality, together with the
capability to process large data sets and efficient visualization of sequence
variation, makes IBROWSER a valuable breeding tool.


  We study Eigen's quasispecies model in the asymptotic regime where the length
of the genotypes goes to infinity and the mutation probability goes to 0. We
give several explicit formulas for the stationary solutions of the limiting
system of differential equations.


  Gene expression is controlled primarily by interactions between transcription
factor proteins (TFs) and the regulatory DNA sequence, a process that can be
captured well by thermodynamic models of regulation. These models, however,
neglect regulatory crosstalk: the possibility that non-cognate TFs could
initiate transcription, with potentially disastrous effects for the cell. Here
we estimate the importance of crosstalk, suggest that its avoidance strongly
constrains equilibrium models of TF binding, and propose an alternative
non-equilibrium scheme that implements kinetic proofreading to suppress
erroneous initiation. This proposal is consistent with the observed covalent
modifications of the transcriptional apparatus and would predict increased
noise in gene expression as a tradeoff for improved specificity. Using
information theory, we quantify this tradeoff to find when optimal proofreading
architectures are favored over their equilibrium counterparts.


  We investigate the role of migration patterns on the spread of epidemics in
complex networks. We enhance the SIS-diffusion model on metapopulations to a
nonlinear diffusion. Specifically, individuals move randomly over the network
but at a rate depending on the population of the departure patch. In the
absence of epidemics, the migration-driven equilibrium is described by
quantifying the total number of individuals living in heavily/lightly populated
areas. Our analytical approach reveals that strengthening the migration from
populous areas contains the infection at the early stage of the epidemic.
Moreover, depending on the exponent of the nonlinear diffusion rate, epidemic
outbreaks do not always occur in the most populated areas as one might expect.


  This paper proposes a control theoretic framework to model and analyze the
self-organized pattern formation of molecular concentrations in biomolecular
communication networks, emerging applications in synthetic biology. In
biomolecular communication networks, bionanomachines, or biological cells,
communicate with each other using a cell-to-cell communication mechanism
mediated by a diffusible signaling molecule, thereby the dynamics of molecular
concentrations are approximately modeled as a reaction-diffusion system with a
single diffuser. We first introduce a feedback model representation of the
reaction-diffusion system and provide a systematic local stability/instability
analysis tool using the root locus of the feedback system. The instability
analysis then allows us to analytically derive the conditions for the
self-organized spatial pattern formation, or Turing pattern formation, of the
bionanomachines. We propose a novel synthetic biocircuit motif called
activator-repressor-diffuser system and show that it is one of the minimum
biomolecular circuits that admit self-organized patterns over cell population.


  The way to infer well-supported phylogenetic trees that precisely reflect the
evolutionary process is a challenging task that completely depends on the way
the related core genes have been found. In previous computational biology
studies, many similarity based algorithms, mainly dependent on calculating
sequence alignment matrices, have been proposed to find them. In these kinds of
approaches, a significantly high similarity score between two coding sequences
extracted from a given annotation tool means that one has the same genes. In a
previous work article, we presented a quality test approach (QTA) that improves
the core genes quality by combining two annotation tools (namely NCBI, a
partially human-curated database, and DOGMA, an efficient annotation algorithm
for chloroplasts). This method takes the advantages from both sequence
similarity and gene features to guarantee that the core genome contains correct
and well-clustered coding sequences (\emph{i.e.}, genes). We then show in this
article how useful are such well-defined core genes for biomolecular
phylogenetic reconstructions, by investigating various subsets of core genes at
various family or genus levels, leading to subtrees with strong bootstraps that
are finally merged in a well-supported supertree.


  This work investigated the stability and asymptotic behavior of some Lotka
Volterra type models. We used the Liapunov method which consists in analyzing
the stability of systems of ordinary differential equations (ODEs) around the
equilibrium when they submitted to perturbations in the initial conditions


  To cope with the complexity of large networks, a number of dimensionality
reduction techniques for graphs have been developed. However, the extent to
which information is lost or preserved when these techniques are employed has
not yet been clear. Here we develop a framework, based on algorithmic
information theory, to quantify the extent to which information is preserved
when network motif analysis, graph spectra and spectral sparsification methods
are applied to over twenty different biological and artificial networks. We
find that the spectral sparsification is highly sensitive to high number of
edge deletion, leading to significant inconsistencies, and that graph spectral
methods are the most irregular, capturing algebraic information in a condensed
fashion but largely losing most of the information content of the original
networks. However, the approach shows that network motif analysis excels at
preserving the relative algorithmic information content of a network, hence
validating and generalizing the remarkable fact that despite their inherent
combinatorial possibilities, local regularities preserve information to such an
extent that essential properties are fully recoverable across different
networks to determine their family group to which they belong to (eg genetic vs
social network). Our algorithmic information methodology thus provides a
rigorous framework enabling a fundamental assessment and comparison between
different data dimensionality reduction methods thereby facilitating the
identification and evaluation of the capabilities of old and new methods.


  Genetic interaction can be defined as a deviation of the phenotypic
quantitative effect of a double gene mutation from the effect predicted from
single mutations using a simple (e.g., multiplicative or linear additive)
statistical model. Experimentally characterized genetic interaction networks in
model organisms provide important insights into relationships between different
biological functions. We describe a computational methodology allowing to
systematically and quantitatively characterize a Boolean mathematical model of
a biological network in terms of genetic interactions between all loss of
function and gain of function mutations with respect to all model phenotypes or
outputs. We use the probabilistic framework defined in MaBoSS software, based
on continuous time Markov chains and stochastic simulations. In addition, we
suggest several computational tools for studying the distribution of double
mutants in the space of model phenotype probabilities. We demonstrate this
methodology on three published models for each of which we derive the genetic
interaction networks and analyze their properties. We classify the obtained
interactions according to their class of epistasis, dependence on the chosen
initial conditions and phenotype. The use of this methodology for validating
mathematical models from experimental data and designing new experiments is
discussed.


  Information processing at the molecular scale is limited by thermal
fluctuations. This can cause undesired consequences in copying information
since thermal noise can lead to errors that can compromise the functionality of
the copy. For example, a high error rate during DNA duplication can lead to
cell death. Given the importance of accurate copying at the molecular scale, it
is fundamental to understand its thermodynamic features. In this paper, we
derive a universal expression for the copy error as a function of entropy
production and {\cred work dissipated by the system during wrong
incorporations}. Its derivation is based on the second law of thermodynamics,
hence its validity is independent of the details of the molecular machinery, be
it any polymerase or artificial copying device. Using this expression, we find
that information can be copied in three different regimes. In two of them, work
is dissipated to either increase or decrease the error. In the third regime,
the protocol extracts work while correcting errors, reminiscent of a Maxwell
demon. As a case study, we apply our framework to study a copy protocol
assisted by kinetic proofreading, and show that it can operate in any of these
three regimes. We finally show that, for any effective proofreading scheme,
error reduction is limited by the chemical driving of the proofreading
reaction.


  We present an analytically solvable model for self-assembly of a molecular
complex on a filament. The process is driven by a seed molecule that undergoes
facilitated diffusion, which is a search strategy that combines diffusion in
three-dimensions and one-dimension. Our study is motivated by single molecule
level observations revealing the dynamics of transcription factors that bind to
the DNA at early stages of transcription. We calculate the probability that a
complex made up of a given number of molecules is completely formed, as well as
the distribution of completion times, upon the binding of a seed molecule at a
target site on the filament (without explicitly modeling the three-dimensional
diffusion that precedes binding). We compare two different mechanisms of
assembly where molecules bind in sequential and random order. Our results
indicate that while the probability of completion is greater for random
binding, the completion time scales exponentially with the size of the complex;
in contrast, it scales as a power-law or slower for sequential binding,
asymptotically. Furthermore, we provide model predictions for the dissociation
and residence times of the seed molecule, which are observables accessible in
single molecule tracking experiments.


  Haplogroup N-M231 of human Y chromosome is a common clade from Eastern Asia
to Northern Europe, being one of the most frequent haplogroups in Altaic and
Uralic-speaking populations. Using newly discovered bi-allelic markers from
high-throughput DNA sequencing, we largely improved the phylogeny of Haplogroup
N, in which 16 subclades could be identified by 33 SNPs. More than 400 males
belonging to Haplogroup N in 34 populations in China were successfully
genotyped, and populations in Northern Asia and Eastern Europe were also
compared together. We found that all the N samples were typed as inside either
clade N1-F1206 (including former N1a-M128, N1b-P43 and N1c-M46 clades), most of
which were found in Altaic, Uralic, Russian and Chinese-speaking populations,
or N2-F2930, common in Tibeto-Burman and Chinese-speaking populations. Our
detailed results suggest that Haplogroup N developed in the region of China
since the final stage of late Paleolithic Era.


  Summary: FermiKit is a variant calling pipeline for Illumina data. It de novo
assembles short reads and then maps the assembly against a reference genome to
call SNPs, short insertions/deletions (INDELs) and structural variations (SVs).
FermiKit takes about one day to assemble 30-fold human whole-genome data on a
modern 16-core server with 85GB RAM at the peak, and calls variants in half an
hour to an accuracy comparable to the current practice. FermiKit assembly is a
reduced representation of raw data while retaining most of the original
information.
  Availability and implementation: https://github.com/lh3/fermikit
  Contact: hengli@broadinstitute.org


  How do mammalian cells that share the same genome exist in notably distinct
phenotypes, exhibiting differences in morphology, gene expression patterns, and
epigenetic chromatin statuses? Furthermore how do cells of different phenotypes
differentiate reproducibly from a single fertilized egg? These are fundamental
problems in developmental biology. Epigenetic histone modifications play an
important role in the maintenance of different cell phenotypes. The exact
molecular mechanism for inheritance of the modification patterns over cell
generations remains elusive. The complexity comes partly from the number of
molecular species and the broad time scales involved. In recent years
mathematical modeling has made significant contributions on elucidating the
molecular mechanisms of DNA methylation and histone covalent modification
inheritance. We will pedagogically introduce the typical procedure and some
technical details of performing a mathematical modeling study, and discuss
future developments.


  In the last decade a number of algorithms and associated software have been
developed to align next generation sequencing (NGS) reads with relevant
reference genomes. The accuracy of these programs may vary significantly,
especially when the NGS reads are quite different from the available reference
genome. We propose a benchmark to assess accuracy of short reads mapping based
on the pre-computed global alignment of related genome sequences.
  In this paper we propose a benchmark to assess accuracy of the short reads
mapping based on the pre-computed global alignment of closely related genome
sequences. We outline the method and also present a short report of an
experiment performed on five popular alignment tools based on the pairwise
alignments of Escherichia coli O6 CFT073 genome with genomes of seven other
bacteria.


  The stability of ecosystems as well as the relation between topology and
dynamics on multilayer networks are important questions that are usually
discussed in separate communities. Here, we combine these two topics by
investigating the influence of the topology of the migration network on the
stability of a four-species foodweb module on six patches. The parameters are
chosen such that the dynamics on an isolated patch have a periodic attractor
with all four species present as well as an attractor where the prey that is
preferred by the top predator dies out. The stability measure used here is
robustness, which is the average proportion of surviving species in the system,
and which shows a complex dependence on the migration rate. We use principal
component analysis to quantify the migration network structure in terms of the
most relevant network measures, and we evaluate correlations between these
measures and characteristics of the robustness curves. Our most important
findings are that higher connectivity of the migration network leads to a
larger maximum robustness, that a broad distribution of connectivities favors
extinction of the preferred prey at intermediate migration rates, and that
migration topologies with a larger betweenness centrality are more prone to
extinction of the preferred prey at the onset of synchronization. Our study
thus demonstrates a strong correlation between dynamical robustness and spatial
topology and can serve as an example for similar studies in other types of
multilayer networks.


  Sexual reproductive behavior has a necessary social coordination component as
willing and capable partners must both be in the right place at the right time.
While there are many known social behavioral adaptations to support solutions
to this problem, we explore the possibility and likelihood of solutions that
rely only on non-social mechanisms. We find three kinds of social organization
that help solve this social coordination problem (herding, assortative mating,
and natal philopatry) emerge in populations of simulated agents with no social
mechanisms available to support these organizations. We conclude that the
non-social origins of these social organizations around sexual reproduction may
provide the environment for the development of social solutions to the same and
different problems.


  Mechanistic multi-stage models are used to analyze lung-cancer mortality
after Plutonium exposure in the Mayak-workers cohort, with follow-up until
2008. Besides the established two-stage model with clonal expansion, models
with three mutation stages as well as a model with two distinct pathways to
cancer are studied. The results suggest that three-stage models offer an
improved description of the data. The best-fitting models point to a mechanism
where radiation increases the rate of clonal expansion. This is interpreted in
terms of changes in cell-cycle control mediated by bystander signaling or
repopulation following cell killing. No statistical evidence for a two-pathway
model is found. To elucidate the implications of the different models for
radiation risk, several exposure scenarios are studied. Models with a radiation
effect at an early stage show a delayed response and a pronounced drop-off with
older ages at exposure. Moreover, the dose-response relationship is strongly
nonlinear for all three-stage models, revealing a marked increase above a
critical dose.


  Statistically consistent estimation of phylogenetic trees or gene trees is
possible if pairwise sequence dissimilarities can be converted to a set of
distances that are proportional to the true evolutionary distances. Susko et
al. (2004) reported some strikingly broad results about the forms of
inconsistency in tree estimation that can arise if corrected distances are not
proportional to the true distances. They showed that if the corrected distance
is a concave function of the true distance, then inconsistency due to long
branch attraction will occur. If these functions are convex, then two "long
branch repulsion" trees will be preferred over the true tree -- though these
two incorrect trees are expected to be tied as the preferred true. Here we
extend their results, and demonstrate the existence of a tree shape (which we
refer to as a "twisted Farris-zone" tree) for which a single incorrect tree
topology will be guaranteed to be preferred if the corrected distance function
is convex. We also report that the standard practice of treating gaps in
sequence alignments as missing data is sufficient to produce non-linear
corrected distance functions if the substitution process is not independent of
the insertion/deletion process. Taken together, these results imply
inconsistent tree inference under mild conditions. For example, if some
positions in a sequence are constrained to be free of substitutions and
insertion/deletion events while the remaining sites evolve with independent
substitutions and insertion/deletion events, then the distances obtained by
treating gaps as missing data can support an incorrect tree topology even given
an unlimited amount of data.


  To adapt their behaviour in changing environments, cells sense concentrations
by binding external ligands to their receptors. However, incorrect ligands may
bind nonspecifically to receptors, and when their concentration is large, this
binding activity may interfere with the sensing of the ligand of interest.
Here, I derive analytically the physical limit to the accuracy of concentration
sensing amid a large number of interfering ligands. A scaling transition is
found when the mean bound time of correct ligands is twice that of incorrect
ligands. I discuss how the physical bound can be approached by a cascade of
receptor states generalizing kinetic proof-reading schemes.


  Protein data bank entries obtain distinct, reproducible flexibility
characteristics determined by normal mode analyses of their three dimensional
coordinate files. We study the effectiveness and sensitivity of this technique
by analyzing the results on one class of glycosidases: family 10 xylanases. A
conserved tryptophan that appears to affect access to the active site can be in
one of two conformations according to X-ray crystallographic electron density
data. The two alternate orientations of this active site tryptophan lead to
distinct flexibility spectra, with one orientation thwarting the oscillations
seen in the other. The particular orientation of this sidechain furthermore
affects the appearance of the motility of a distant, C terminal region we term
the mallet. The mallet region is known to separate members of this family of
enzymes into two classes.


  This paper describes a technique for analyzing the stochastic structure of
the vectorial capacity using moment--generating functions. In such formulation,
for an infectious disease transmitted by a vector, we obtain the generating
function for the distribution of the number of infectious contacts (e.g.,
infectious bites from mosquitoes) between vectors and humans after a contact to
a single infected individual. This approach permits us to generally derive the
moments of the distribution and, under some conditions, derive the distribution
function of the vectorial capacity. A stochastic modeling framework is helpful
for analyzing the dynamics of disease spreading, such as performing sensitivity
analysis.


  It is generally considered that animal life was triggered by the rise of
oxygen levels. Based on experiments evaluating the minimum range of oxygen
levels at which sponges can survive, Mills and coauthors
(doi:10.1073/pnas.1400547111) defend the opposite view. However, the authors do
not demonstrate that "animal life was not triggered by the oxygen rise" is the
only possible and unique conclusion from their observation. In this reply, it
is suggested that a mechanism to explain the ability of sponges to survive at
low oxygen biota is Ocean Anoxic Events. These lead to oxygen depletion and a
series of them would selectively favor animals able to survive at low oxygen
levels. Thus, the origin of the ability of marine animals to survive in low
oxygen biota remains to be clarified.


  Riboswitches are cis-acting regulatory RNA elements prevalently located in
the leader sequences of bacterial mRNA. An adenine sensing riboswitch
cis-regulates adeninosine deaminase gene (add) in Vibrio vulnificus. The
structural mechanism regulating its conformational changes upon ligand binding
mostly remains to be elucidated. In this open framework it has been suggested
that the ligand stabilizes the interaction of the distal "kissing loop"
complex. Using accurate full-atom molecular dynamics with explicit solvent in
combination with enhanced sampling techniques and advanced analysis methods it
could be possible to provide a more detailed perspective on the formation of
these tertiary contacts. In this work, we used umbrella sampling simulations to
study the thermodynamics of the kissing loop complex in the presence and in the
absence of the cognate ligand. We enforced the breaking/formation of the
loop-loop interaction restraining the distance between the two loops. We also
assessed the convergence of the results by using two alternative initialization
protocols. A structural analysis was performed using a novel approach to
analyze base contacts. Our simulations qualitatively indicated that the ligand
could stabilize the kissing loop complex. We also compared with previously
published simulation studies. Kissing complex stabilization given by the ligand
was compatible with available experimental data. However, the dependence of its
value on the initialization protocol of the umbrella sampling simulations posed
some questions on the quantitative interpretation of the results and called for
better converged enhanced sampling simulations.


  Rapid rebinding of molecular interaction partners that are in close proximity
after dissociation leads to a dissociation and association kinetics that can
profoundly differ from predictions based on bulk reaction models. The cause of
this effect can be traced back to the non-Markovian character of the ligand's
rebinding time probability density function, reflecting the fact that, for a
certain time span, the ligand still 'remembers' the receptor it was bound to
previously. In this manuscript, we explore the consequences of the hypothesis
that initial binding and consecutive rebinding give rise to a bond lifetime
density that is non-Markovian as well. We study the combined effect of the two
non-Markovian waiting time probability densities and show that even for very
short times the decay of the fraction of occupied receptors deviates from an
exponential. For long times, dissociation is slower than an exponential and the
fate of the the steady-state bound receptor fraction critically depends on the
extent of the deviation, relative to the rebinding time density, from the
Markovian limit: The population of occupied receptors may either decay
completely or assume a non-vanishing value for small and strong deviations,
respectively. Furthermore, we point out the important role played by fractional
calculus and demonstrate that the short- and long-time dynamics of the occupied
receptors can be naturally expressed as well as easily obtained in terms of
fractional differential equations involving the Riemann-Liouville derivative.
Our analysis shows that cells may exploit receptor memory as mechanism to
dynamically widen the range of potential response-patterns to a given signal.


  Biological organisms adapt to changes by processing informations from
different sources, most notably from their ancestors and from their
environment. We review an approach to quantify these informations by analyzing
mathematical models of evolutionary dynamics, and show how explicit results are
obtained for a solvable subclass of these models. In several limits, the
results coincide with those obtained in studies of information processing for
communication, gambling or thermodynamics. In the most general case, however,
information processing by biological populations shows unique features that
motivate the analysis of specific models.


  Condensation phenomena arise through a collective behaviour of particles.
They are observed in both classical and quantum systems, ranging from the
formation of traffic jams in mass transport models to the macroscopic
occupation of the energetic ground state in ultra-cold bosonic gases
(Bose-Einstein condensation). Recently, it has been shown that a driven and
dissipative system of bosons may form multiple condensates. Which states become
the condensates has, however, remained elusive thus far. The dynamics of this
condensation are described by coupled birth-death processes, which also occur
in evolutionary game theory. Here, we apply concepts from evolutionary game
theory to explain the formation of multiple condensates in such
driven-dissipative bosonic systems. We show that vanishing of relative entropy
production determines their selection. The condensation proceeds exponentially
fast, but the system never comes to rest. Instead, the occupation numbers of
condensates may oscillate, as we demonstrate for a rock-paper-scissors game of
condensates.


  We discuss the symbolic dynamics of biochemical networks with separate
timescales. We show that symbolic dynamics of monomolecular reaction networks
with separated rate constants can be described by deterministic, acyclic
automata with a number of states that is inferior to the number of biochemical
species. For nonlinear pathways, we propose a general approach to approximate
their dynamics by finite state machines working on the metastable states of the
network (long life states where the system has slow dynamics). For networks
with polynomial rate functions we propose to compute metastable states as
solutions of the tropical equilibration problem. Tropical equilibrations are
defined by the equality of at least two dominant monomials of opposite signs in
the differential equations of each dynamic variable. In algebraic geometry,
tropical equilibrations are tantamount to tropical prevarieties, that are
finite intersections of tropical hypersurfaces.


  Positioning of nucleosomes along eukaryotic genomes plays an important role
in their organization and regulation. There are many different factors
affecting the location of nucleosomes. Some can be viewed as preferential
binding of a single nucleosome to different locations along the DNA and some as
interactions between neighboring nucleosomes. In this study we analyzed how
well nucleosomes are positioned along the DNA as a function of strength of the
preferential binding, correlation length of the binding energy landscape,
interactions between neighboring nucleosomes and others relevant system
properties. We analyze different scenarios: designed energy landscapes and
generically disordered ones and derive conditions for good positioning. Using
analytic and numerical approaches we find that, even if the binding preferences
are very weak, synergistic interplay between the interactions and the binding
preferences is essential for a good positioning of nucleosomes, especially on
correlated energy landscapes. Analyzing empirical energy landscape, we discuss
relevance of our theoretical results to positioning of nucleosomes on DNA
\emph{in vivo.}


  Being capable of characterizing DNA local bending is essential to understand
thoroughly many biological processes because they involve a local bending of
the double helix axis, either intrinsic to the sequence or induced by the
binding of proteins. Developing a method to measure DNA bend angles that does
not perturb the conformation of the DNA itself or the DNA-protein complex is a
challenging task. Here, we propose a joint theory-experiment high throughput
approach to rigorously measure such bend angles using the Tethered Particle
Motion (TPM) technique. By carefully modeling the TPM geometry, we propose a
simple formula based on a kinked Worm-Like Chain model to extract the bend
angle from TPM measurements. Using constructs made of 575 base-pair DNAs with
in-phase assemblies of 1 to 7 6A-tracts, we find that the sequence CA6CGG
induces a bend angle of 19 [4] {\deg}. Our method is successfully compared to
more theoretically complex or experimentally invasive ones such as cyclization,
NMR, FRET or AFM. We further apply our procedure to TPM measurements from the
literature and demonstrate that the angles of bends induced by proteins, such
as Integration Host Factor (IHF) can be reliably evaluated as well.


  Nature presents multiple intriguing examples of processes which proceed at
high precision and regularity. This remarkable stability is frequently counter
to modelers' experience with the inherent stochasticity of chemical reactions
in the regime of low copy numbers. Moreover, the effects of noise and
nonlinearities can lead to "counter-intuitive" behavior, as demonstrated for a
basic enzymatic reaction scheme that can display stochastic focusing (SF).
Under the assumption of rapid signal fluctuations, SF has been shown to convert
a graded response into a threshold mechanism, thus attenuating the detrimental
effects of signal noise. However, when the rapid fluctuation assumption is
violated, this gain in sensitivity is generally obtained at the cost of very
large product variance, and this unpredictable behavior may be one possible
explanation of why, more than a decade after its introduction, SF has still not
been observed in real biochemical systems.
  In this work we explore the noise properties of a simple enzymatic reaction
mechanism with a small and fluctuating number of active enzymes that behaves as
a high-gain, noisy amplifier due to SF caused by slow enzyme fluctuations. We
then show that the inclusion of a plausible negative feedback mechanism turns
the system from a noisy signal detector to a strong homeostatic mechanism by
exchanging high gain with strong attenuation in output noise and robustness to
parameter variations. Moreover, we observe that the discrepancy between
deterministic and stochastic descriptions of stochastically focused systems in
the evolution of the means almost completely disappears, despite very low
molecule counts and the additional nonlinearity due to feedback.
  The reaction mechanism considered here can provide a possible resolution to
the apparent conflict between intrinsic noise and high precision in critical
intracellular processes.


  Many existing studies on pattern formation in the reaction-diffusion systems
rely on deterministic models. However, environmental noise is often a major
factor which leads to significant changes in the spatiotemporal dynamics. In
this paper, we focus on the spatiotemporal patterns produced by the
predator-prey model with ratio-dependent functional response and density
dependent death rate of predator. We get the reaction-diffusion equations
incorporating the self-diffusion terms, corresponding to random movement of the
individuals within two dimensional habitats, into the growth equations for the
prey and predator population. In order to have to have the noise added model,
small amplitude heterogeneous perturbations to the linear intrinsic growth
rates are introduced using uncorrelated Gaussian white noise terms. For the
noise added system, we then observe spatial patterns for the parameter values
lying outside the Turing instability region. With thorough numerical
simulations we characterize the patterns corresponding to Turing and
Turing-Hopf domain and study their dependence on different system parameters
like noise-intensity, etc.


  It is has been assumed for over half a century that the earliest animals were
obligate aerobes with relatively high oxygen requirements. However, the
conserved biochemistry and widespread phylogenetic distribution of anaerobic
energy metabolism in animals suggests a deep ancestral possession of the genes
and enzymes necessary for a facultative anaerobic lifestyle. Additionally,
non-bilaterian bodyplans are not expected to require particularly high
environmental oxygen levels. This is consistent with experimental evidence
demonstrating the low-oxygen tolerance of the sponge Halichondria panicea.
While it is conceivable that low-oxygen-adapted animals evolved only sometime
during the past 541 million years, perhaps in response to oceanic anoxic
events, they most reasonably date back to the first animals themselves, as the
last common ancestor of animals likely emerged in a relatively low-oxygen
world, possessed the genetic means for anaerobiosis, and exhibited a bodyplan
conducive to aerobic growth under oxygen levels less than 4% of modern
atmospheric saturation.


  Ethnocentrism is a behavioral strategy seen on every scale of social
interaction. Game-theory models demonstrate that evolution selects
ethnocentrism because it boosts cooperation, which increases reproductive
fitness. However, some believe that interethnic unions have the potential to
foster universal cooperation and overcome in-group biases in humans. Here, we
use agent-based computer simulations to test this hypothesis. Cooperative
intergroup mating does lend an advantage to a universal cooperation strategy
when the cost/benefit ratio of cooperation is low and local population
diversity is high.


  Knowing the strategy of an opponent in a competitive environment conveys
obvious evolutionary advantages. But this information is costly, and the
benefit of being informed may not necessarily offset the additional cost. Here
we introduce social dilemmas with informed strategies, and we show that this
gives rise to two cyclically dominant triplets that form defensive alliances.
The stability of these two alliances is determined by the rotation velocity of
the strategies within each triplet. A weaker strategy in a faster rotating
triplet can thus overcome an individually stronger competitor. Fascinating
spatial patterns favor the dominance of a single defensive alliance, but enable
also the stable coexistence of both defensive alliances in very narrow regions
of the parameter space. A continuous reentrant phase transition reveals before
unseen complexity behind the stability of strategic alliances in evolutionary
social dilemmas.


  The growth of a population divided among spatial sites, with migration
between the sites, is sometimes modelled by a product of random matrices, with
each diagonal elements representing the growth rate in a given time period, and
off-diagonal elements the migration rate. If the sites are reinterpreted as age
classes, the same model may apply to a single population with age-dependent
mortality and reproduction.
  We consider the case where the off-diagonal elements are small, representing
a situation where there is little migration or, alternatively, where a
deterministic life-history has been slightly disrupted, for example by
introducing a rare delay in development. We examine the asymptotic behaviour of
the long-term growth rate. We show that when the highest growth rate is
attained at two different sites in the absence of migration (which is always
the case when modelling a single age-structured population) the increase in
stochastic growth rate due to a migration rate $\epsilon$ is like $(\log
\epsilon^{-1})^{-1}$ as $\epsilon\downarrow 0$, under fairly generic
conditions. When there is a single site with the highest growth rate the
behavior is more delicate, depending on the tails of the growth rates. For the
case when the log growth rates have Gaussian-like tails we show that the
behavior near zero is like a power of $\epsilon$, and derive upper and lower
bounds for the power in terms of the difference in the growth rates and the
distance between the sites.


  Culinary systems, the practice of preparing a refined combination of
ingredients that is palatable as well as socially acceptable, are examples of
complex dynamical systems. They evolve over time and are affected by a large
number of factors. Modeling the dynamic nature of evolution of regional
cuisines may provide us a quantitative basis and exhibit underlying processes
that have driven them into the present day status. This is especially important
given that the potential culinary space is practically infinite because of
possible number of ingredient combinations as recipes. Such studies also
provide a means to compare and contrast cuisines and to unearth their
therapeutic value. Herein we provide rigorous analysis of modeling eight
diverse Indian regional cuisines, while also highlighting their uniqueness, and
a comparison among those models at the level of flavor compounds which opens up
molecular level studies associating them especially with non-communicable
diseases such as diabetes.