  Approximations based on moment-closure (MA) are commonly used to obtain
estimates of the mean molecule numbers and of the variance of fluctuations in
the number of molecules of chemical systems. The advantage of this approach is
that it can be far less computationally expensive than exact stochastic
simulations of the chemical master equation. Here we numerically study the
conditions under which the MA equations yield results reflecting the true
stochastic dynamics of the system. We show that for bistable and oscillatory
chemical systems with deterministic initial conditions, the solution of the MA
equations can be interpreted as a valid approximation to the true moments of
the CME, only when the steady-state mean molecule numbers obtained from the
chemical master equation fall within a certain finite range. The same validity
criterion for monostable systems implies that the steady-state mean molecule
numbers obtained from the chemical master equation must be above a certain
threshold. For mean molecule numbers outside of this range of validity, the MA
equations lead to either qualitatively wrong oscillatory dynamics or to
unphysical predictions such as negative variances in the molecule numbers or
multiple steady-state moments of the stationary distribution as the initial
conditions are varied. Our results clarify the range of validity of the MA
approach and show that pitfalls in the interpretation of the results can only
be overcome through the systematic comparison of the solutions of the MA
equations of a certain order with those of higher orders.


  Genome-wide association studies (GWAS) have identified many genetic factors
underlying complex human traits. However, these factors have explained only a
small fraction of these traits' genetic heritability. It is argued that many
more genetic factors remain undiscovered. These genetic factors likely are
weakly associated at the population level and sparsely distributed across the
genome. In this paper, we adapt the recent innovations on Tukey's Higher
Criticism (Tukey [The Higher Criticism (1976) Princeton Univ.]; Donoho and Jin
[Ann. Statist. 32 (2004) 962-994]) to SNP-set analysis of GWAS, and develop a
new theoretical framework in large-scale inference to assess the joint
significance of such rare and weak effects for a quantitative trait. In the
core of our theory is the so-called detection boundary, a curve in the
two-dimensional phase space that quantifies the rarity and strength of genetic
effects. Above the detection boundary, the overall effects of genetic factors
are strong enough for reliable detection. Below the detection boundary, the
genetic factors are simply too rare and too weak for reliable detection. We
show that the HC-type methods are optimal in that they reliably yield detection
once the parameters of the genetic effects fall above the detection boundary
and that many commonly used SNP-set methods are suboptimal. The superior
performance of the HC-type approach is demonstrated through simulations and the
analysis of a GWAS data set of Crohn's disease.


  Large case/control Genome-Wide Association Studies (GWAS) often include
groups of related individuals with known relationships. When testing for
associations at a given locus, current methods incorporate only the familial
relationships between individuals. Here, we introduce the chromosome-based
Quasi Likelihood Score (cQLS) statistic that incorporates local
Identity-By-Descent (IBD) to increase the power to detect associations. In
studies robust to population stratification, such as those with case/control
sibling pairs, simulations show that the study power can be increased by over
50%. In our example, a GWAS examining late-onset Alzheimer's disease, the
$p$-values among the most strongly associated SNPs in the APOE gene tend to
decrease, with the smallest $p$-value decreasing from $1.23\times10^{-8}$ to
$7.70\times 10^{-9}$. Furthermore, as a part of our simulations, we reevaluate
our expectations about the use of families in GWAS. We show that, although
adding only half as many unique chromosomes, genotyping affected siblings is
more efficient than genotyping randomly ascertained cases. We also show that
genotyping cases with a family history of disease will be less beneficial when
searching for SNPs with smaller effect sizes.


  The relationship between the Moran model and stochastic Lotka-Volterra
competition (SLVC) model is explored via timescale separation arguments. For
neutral systems the two are found to be equivalent at long times. For systems
with selective pressure, their behavior differs. It is argued that the SLVC is
preferable to the Moran model since in the SLVC population size is regulated by
competition, rather than arbitrarily fixed as in the Moran model. As a
consequence, ambiguities found in the Moran model associated with the
introduction of more complex processes, such as selection, are avoided.


  Protocells are supposed to have played a key role in the self-organizing
processes leading to the emergence of life. Existing models either (i) describe
protocell architecture and dynamics, given the existence of sets of
collectively self-replicating molecules for granted, or (ii) describe the
emergence of the aforementioned sets from an ensemble of random molecules in a
simple experimental setting (e.g. a closed system or a steady-state flow
reactor) that does not properly describe a protocell. In this paper we present
a model that goes beyond these limitations by describing the dynamics of sets
of replicating molecules within a lipid vesicle. We adopt the simplest possible
protocell architecture, by considering a semi-permeable membrane that selects
the molecular types that are allowed to enter or exit the protocell and by
assuming that the reactions take place in the aqueous phase in the internal
compartment. As a first approximation, we ignore the protocell growth and
division dynamics. The behavior of catalytic reaction networks is then
simulated by means of a stochastic model that accounts for the creation and the
extinction of species and reactions. While this is not yet an exhaustive
protocell model, it already provides clues regarding some processes that are
relevant for understanding the conditions that can enable a population of
protocells to undergo evolution and selection.


  We suggest a novel approach to treating symbiotic relations between
biological species or social entities. The main idea is the characterisation of
symbiotic relations of coexisting species through their mutual influence on
their respective carrying capacities, taking into account that this influence
can be quite strong and requires a nonlinear functional framework. We
distinguish three variants of mutual influence, representing the main types of
relations between species: (i) passive symbiosis, when the mutual carrying
capacities are influenced by other species without their direct interactions;
(ii) active symbiosis, when the carrying capacities are transformed by
interacting species; and (iii) mixed symbiosis, when the carrying capacity of
one species is influenced by direct interactions, while that of the other
species is not. The approach allows us to describe all kinds of symbiosis,
mutualism, commensalism, and parasitism within a unified scheme. The case of
two symbiotic species is analysed in detail, demonstrating several dynamical
regimes of coexistence, unbounded growth of both populations, growth of one and
elimination of the other population, convergence to evolutionary stable states,
and everlasting population oscillations. The change of the dynamical regimes
occurs by varying the system parameters characterising the destruction or
creation of the mutual carrying capacities. The regime changes are associated
with several dynamical system bifurcations.


  Here, we provide revised gene models for D. ananassae, D. yakuba, and D.
simulans, which include UTRs and empirically verified intron-exon boundaries,
as well as ortholog groups identified using a fuzzy reciprocal-best-hit blast
comparison. Using these revised annotations, we perform differential expression
testing using the cufflinks suite to provide a broad overview of differential
expression between reproductive tissues and the carcass. We identify thousands
of genes that are differentially expressed across tissues in D. yakuba and D.
simulans, with roughly 60% agreement in expression patterns of orthologs in D.
yakuba and D. simulans. We identify several cases of putative polycistronic
transcripts, pointing to a combination of transcriptional read-through in the
genome as well as putative gene fusion and fission events across taxa. We
furthermore identify hundreds of lineage specific genes in each species with no
blast hits among transcripts of any other Drosophila species, which are
candidates for neofunctionalized proteins and a potential source of genetic
novelty.


  It is becoming clear that traditional, single-structure models of proteins
are insufficient for understanding their biological function. Here, we outline
one method for inferring, from experiments, not only the most common structure
a protein adopts (native state), but the entire ensemble of conformations the
system can adopt. Such ensemble mod- els are necessary to understand
intrinsically disordered proteins, enzyme catalysis, and signaling. We suggest
that the most difficult aspect of generating such a model will be finding a
small set of configurations to accurately model structural heterogeneity and
present one way to overcome this challenge.


  We develop a new perturbation method for studying quasi-neutral competition
in a broad class of stochastic competition models, and apply it to the analysis
of fixation of competing strains in two epidemic models. The first model is a
two-strain generalization of the stochastic Susceptible-Infected-Susceptible
(SIS) model. Here we extend previous results due to Parsons and Quince (2007),
Parsons et al (2008) and Lin, Kim and Doering (2012). The second model, a
two-strain generalization of the stochastic Susceptible-Infected-Recovered
(SIR) model with population turnover, has not been studied previously. In each
of the two models, when the basic reproduction numbers of the two strains are
identical, a system with an infinite population size approaches a point on the
deterministic coexistence line (CL): a straight line of fixed points in the
phase space of sub-population sizes. Shot noise drives one of the strain
populations to fixation, and the other to extinction, on a time scale
proportional to the total population size. Our perturbation method explicitly
tracks the dynamics of the probability distribution of the sub-populations in
the vicinity of the CL. We argue that, whereas the slow strain has a
competitive advantage for mathematically "typical" initial conditions, it is
the fast strain that is more likely to win in the important situation when a
few infectives of both strains are introduced into a susceptible population.


  Several animal species are considered to exhibit what is called negligible
senescence, i.e. they do not show signs of functional decline or any increase
of mortality with age, and do not have measurable reductions in reproductive
capacity with age. Recent studies in Naked Mole Rat (NMR) and long- lived sea
urchin showed that the level of gene expression changes with age is lower than
in other organisms. These phenotypic observations correlate well with
exceptional endurance of NMR tissues to various genotoxic stresses. Therefore,
the lifelong transcriptional stability of an organism may be a key determinant
of longevity. However, the exact relation between genetic network stability,
stress-resistance and aging has not been defined. We analyze the stability of a
simple genetic- network model of a living organism under the influence of
external and endogenous factors. We demonstrate that under most common
circumstances a gene network is inherently unstable and suffers from
exponential accumulation of gene-regulation deviations leading to death.
However, should the repair systems be sufficiently effective, the gene network
can stabilize so that gene damage remains constrained along with mortality of
the organism, which may then enjoy a remarkable degree of stability over very
long times. We clarify the relation between stress-resistance and aging and
suggest that stabilization of the genetic network may provide a mathematical
explanation of the Gompertz equation describing the relationship between age
and mortality in many species, and of the apparently negligible senescence
observed in exceptionally long-lived animals. The model may support a range of
applications, such as systematic searches for therapeutics to extend lifespan
and healthspan.


  Module network inference is an established statistical method to reconstruct
co-expression modules and their upstream regulatory programs from integrated
multi-omics datasets measuring the activity levels of various cellular
components across different individuals, experimental conditions or time points
of a dynamic process. We have developed Lemon-Tree, an open-source,
platform-independent, modular, extensible software package implementing
state-of-the-art ensemble methods for module network inference. We benchmarked
Lemon-Tree using large-scale tumor datasets and showed that Lemon-Tree
algorithms compare favorably with state-of-the-art module network inference
software. We also analyzed a large dataset of somatic copy-number alterations
and gene expression levels measured in glioblastoma samples from The Cancer
Genome Atlas and found that Lemon-Tree correctly identifies known glioblastoma
oncogenes and tumor suppressors as master regulators in the inferred module
network. Novel candidate driver genes predicted by Lemon-Tree were validated
using tumor pathway and survival analyses. Lemon-Tree is available from
http://lemon-tree.googlecode.com under the GNU General Public License version
2.0.


  A self-regulatory genetic circuit, where a protein acts as a positive
regulator of its own production, is known to be a simplest form of biological
network with a positive feedback loop. Although at least three components, DNA,
RNA, and the protein, are required to form such a circuit, the stability
analysis of fixed points of the self-regulatory circuit has been performed only
after reducing the system into to a two-component system consisting of RNA and
protein only, assuming a fast equilibration of the DNA component. Here, the
stability of fixed points of the three-component positive feedback loop is
analyzed by obtaining eigenvalues of full three dimensional Hessian matrix. In
addition to rigorously identifying the stable fixed points and the saddle
points, detailed information can be obtained, such as the number of positive
eigenvalues near a saddle point. In particular, complex eigenvalues is shown to
exist for sufficiently slow binding and unbinding of the auto-regulatory
transcription factor to DNA, leading to oscillatory convergence to the steady
states, a novel feature unseen in the two-dimensional analysis.


  Transmission rates are key in understanding the spread of infectious
diseases. Using the framework of compartmental models, we introduce a simple
method that enables us to reconstruct time series of transmission rates
directly from incidence or disease-related mortality data. The reconstruction
exploits differential equations, which model the time evolution of infective
stages and strains. Being sensitive to initial values, the method produces
asymptotically correct solutions. The computations are fast, with time
complexity being quadratic. We apply the reconstruction to data of measles
(England and Wales, 1948-67), dengue (Thailand, 1982-99), and influenza (U.S.,
1910-27). The Measles example offers comparison with earlier work. Here we
re-investigate reporting corrections, include and exclude demographic
information. The dengue example deals with the failure of vector-control
measures in reducing dengue hemorrhagic fever (DHF) in Thailand. Two competing
mechanisms have been held responsible: strain interaction and demographic
transitions. Our reconstruction reveals that both explanations are possible,
showing that the increase in DHF cases is consistent with decreasing
transmission rates resulting from reduced vector counts. The flu example
focuses on the 1918/19 pandemic, examining the transmission rate evolution for
an invading strain. Our analysis indicates that the pandemic strain could have
circulated in the population for many months before the pandemic was initiated
by an event of highly increased transmission.


  The rubeola virus, commonly known as measles, is one of the major causes of
vaccine-preventable deaths among children worldwide. This is the case despite
the fact that an effective vaccine is widely available. Even in developed
countries elimination efforts have fallen short as seen by recent outbreaks in
Europe, which had over 30,000 cases reported in 2010. The string of measles
outbreaks in France from 2008-2011 is of particular interest due to the
documented disparity in regional vaccination coverage. The impact of
heterogeneous vaccine coverage on disease transmission is a broad interest and
the focus of this study. A Susceptible-Exposed-Infectious-Recovered (SEIR)
multi-patch epidemiological model capturing the regional differences in
vaccination rates and mixing is introduced. The mathematical analysis of a
two-patch system is carried out to help our understanding of the behavior of
multi-patch systems. Numerical simulations are generated to aid the study of
the system's qualitative dynamics. Data from the recent French outbreaks were
used to generate parameter values and to help connect theory with application.
Our findings show that heterogeneous vaccination coverage increases controlled
reproduction number compared to comparable homogeneous coverage.


  It has been recently proposed that the two "emerging" hallmarks of cancer,
namely altered glucose metabolism and immune evasion, may in fact be
fundamentally linked (Kareva and Hahnfeldt, 2013). This connection comes from
up-regulation of glycolysis by tumor cells, which can lead to active
competition for resources in the tumor microenvironment between tumor and
immune cells. Here it is further proposed that cancer stem cells (CSCs) can
circumvent the anti-tumor immune response by creating a 'protective shield' of
non-stem cancer cells around them. This shield can protect the CSCs both by
creating a physical barrier between them and cytotoxic lymphocytes (CTLs), and
by promoting competition for the common resources, such as glucose, between
non-stem cancer cells and CTLs. The implications of this hypothesis are
investigate using an agent-based model, leading to a prediction that relative
CSC to non-CSC ratio will vary with the strength of the host immune response,
with the highest occurring at an intermediate state of immune activation. A
discussion of possible therapeutic approaches concludes the paper, suggesting
that a chemotherapeutic regimen consisting of regular pulsed doses, i.e.,
metronomic chemotherapy, would yield the best clinical outcome by allowing CTLs
to most effectively reach and eliminate CSCs.


  We present a thermodynamic theory for a generic population of $M$ individuals
distributed into $N$ groups (clusters). We construct the ensemble of all
distributions with fixed $M$ and $N$, introduce a selection functional that
embodies the physics that governs the population, and obtain the distribution
that emerges in the scaling limit as the most probable among all distributions
consistent with the given physics. We develop the thermodynamics of the
ensemble and establish a rigorous mapping to thermodynamics. We treat the
emergence of a so-called "giant component" as a formal phase transition and
show that the criteria for its emergence are entirely analogous to the
equilibrium conditions in molecular systems. We demonstrate the theory by an
analytic model and confirm the predictions by Monte Carlo simulation.


  A Belief Propagation approach has been recently proposed for the zero-patient
problem in a SIR epidemics. The zero-patient problem consists in finding the
initial source of an epidemic outbreak given observations at a later time. In
this work, we study a harder but related inference problem, in which
observations are noisy and there is confusion between observed states. In
addition to studying the zero-patient problem, we also tackle the problem of
completing and correcting the observations possibly finding undiscovered
infected individuals and false test results. Moreover, we devise a set of
equations, based on the variational expression of the Bethe free energy, to
find the zero patient along with maximum-likelihood epidemic parameters. We
show, by means of simulated epidemics, how this method is able to infer details
on the past history of an epidemic outbreak based solely on the topology of the
contact network and a single snapshot of partial and noisy observations.


  Myosin-V is a highly processive dimeric protein that walks with 36nm steps
along actin tracks, powered by coordinated ATP hydrolysis reactions in the two
myosin heads. No previous theoretical models of the myosin-V walk reproduce all
the observed trends of velocity and run-length with [ADP], [ATP] and external
forcing. In particular, a result that has eluded all theoretical studies based
upon rigorous physical chemistry is that run length decreases with both
increasing [ADP] and [ATP]. We systematically analyse which mechanisms in
existing models reproduce which experimental trends and use this information to
guide the development of models that can reproduce them all. We formulate
models as reaction networks between distinct mechanochemical states with
energetically determined transition rates. For each network architecture, we
compare predictions for velocity and run length to a subset of experimentally
measured values, and fit unknown parameters using a bespoke MCSA optimization
routine. Finally we determine which experimental trends are replicated by the
best-fit model for each architecture. Only two models capture them all: one
involving [ADP]-dependent mechanical detachment, and another including
[ADP]-dependent futile cycling and nucleotide pocket collapse. Comparing
model-predicted and experimentally observed kinetic transition rates favors the
latter.


  A key challenge in systems biology is the elucidation of the underlying
principles, or fundamental laws, which determine the cellular phenotype.
Understanding how these fundamental principles are altered in diseases like
cancer is important for translating basic scientific knowledge into clinical
advances. While significant progress is being made, with the identification of
novel drug targets and treatments by means of systems biological methods, our
fundamental systems level understanding of why certain treatments succeed and
others fail is still lacking. We here advocate a novel methodological framework
for systems analysis and interpretation of molecular omic data, which is based
on statistical mechanical principles. Specifically, we propose the notion of
cellular signalling entropy (or uncertainty), as a novel means of analysing and
interpreting omic data, and more fundamentally, as a means of elucidating
systems-level principles underlying basic biology and disease. We describe the
power of signalling entropy to discriminate cells according to differentiation
potential and cancer status. We further argue the case for an empirical
cellular entropy-robustness correlation theorem and demonstrate its existence
in cancer cell line drug sensitivity data. Specifically, we find that high
signalling entropy correlates with drug resistance and further describe how
entropy could be used to identify the achilles heels of cancer cells. In
summary, signalling entropy is a deep and powerful concept, based on rigorous
statistical mechanical principles, which, with improved data quality and
coverage, will allow a much deeper understanding of the systems biological
principles underlying normal and disease physiology.


  Knowledge of the mechanisms of interaction between self-aggregating peptides
and nucleic acids or other polyanions is key to the understanding of many
aggregation processes underlying several human diseases (e.g. Alzheimer's and
Parkinson's diseases). Determining the affinity and kinetic steps of such
interactions is challenging due to the competition between hydrophobic
self-aggregating forces and electrostatic binding forces. Kahalalide F (KF) is
an anticancer hydrophobic peptide which contains a single positive charge that
confers strong aggregative properties with polyanions. This makes KF an ideal
model to elucidate the mechanisms by which self-aggregation competes with
binding to a strongly charged polyelectrolyte such as DNA. We use optical
tweezers to apply mechanical forces to single DNA molecules and show that KF
and DNA interact in a two-step kinetic process promoted by the electrostatic
binding of DNA to the aggregate surface followed by the stabilization of the
complex due to hydrophobic interactions. From the measured pulling curves we
determine the spectrum of binding affinities, kinetic barriers and lengths of
DNA segments sequestered within the KF-DNA complex. We find there is a capture
distance beyond which the complex collapses into compact aggregates stabilized
by strong hydrophobic forces, and discuss how the bending rigidity of the
nucleic acid affects such process. We hypothesize that within an in vivo
context, the enhanced electrostatic interaction of KF due to its aggregation
might mediate the binding to other polyanions. The proposed methodology should
be useful to quantitatively characterize other compounds or proteins in which
the formation of aggregates is relevant.


  A comprehensive approach to cellular differentiation in cyanobacteria is
developed. To this aim, the process of heterocyst cell formation is studied
under a systems biology point of view. By relying on statistical physics
techniques, we translate the essential ingredients and mechanisms of the
genetic circuit into a set of differential equations that describes the
continuous time evolution of combined nitrogen, PatS, HetR and NtcA
concentrations. The detailed analysis of these equations gives insight into the
single cell dynamics. On the other hand, the inclusion of diffusion and noisy
conditions allows simulating the formation of heterocysts patterns in
cyanobacteria strains. The time evolution of relevant component concentrations
are calculated allowing for a comparison with experiments. Finally, we discuss
the validity and the possible improvements of the model.


  Challenges of assessing complexity and clonality in populations of mixed
species arise in diverse areas of modern biology, including estimating
diversity and clonality in microbiome populations, measuring patterns of T and
B cell clonality, and determining the underlying tumor cell population
structure in cancer. Here we address the problem of quantifying populations,
with our analysis directed toward systems for which previously defined
algorithms allow the sequence-based identification of clonal subpopulations.
Data come from replicate sequencing libraries generated from a sample,
potentially with very different depths. While certain properties of the
underlying clonal distribution (most notably the total number of clones) are
difficult to estimate accurately from data representing a small fraction of the
total population, the population-level "clonality" metric that is the sum of
squared probabilities of the respective species can be calculated. (This is the
sum of squared entries of a high-dimensional vector $p$ of relative
frequencies.) The clonality score is the probability of a clonal relationship
between two randomly chosen members of the population of interest. A principal
takeaway message is that knowing a functional of $p$ well may not depend on
knowing $p$ itself very well.
  Our work has led to software, which we call {\it lymphclon}; it has been
deposited in the CRAN library.


  Background: Differential expression analysis of (individual) genes is often
used to study their roles in diseases. However, diseases such as cancer are a
result of the combined effect of multiple genes. Gene products such as proteins
seldom act in isolation, but instead constitute stable multi-protein complexes
performing dedicated functions. Therefore, complexes aggregate the effect of
individual genes (proteins) and can be used to gain a better understanding of
cancer mechanisms. Here, we observe that complexes show considerable changes in
their expression, in turn directed by the concerted action of transcription
factors (TFs), across cancer conditions. We seek to gain novel insights into
cancer mechanisms through a systematic analysis of complexes and their
transcriptional regulation.
  Results: We integrated large-scale protein-interaction (PPI) and
gene-expression datasets to identify complexes that exhibit significant changes
in their expression across different conditions in cancer. We devised a
log-linear model to relate these changes to the differential regulation of
complexes by TFs. The application of our model on two case studies involving
pancreatic and familial breast tumour conditions revealed: (i) complexes in
core cellular processes, especially those responsible for maintaining genome
stability and cell proliferation (e.g. DNA damage repair and cell cycle) show
considerable changes in expression; (ii) these changes include decrease and
countering increase for different sets of complexes indicative of compensatory
mechanisms coming into play in tumours; and (iii) TFs work in cooperative and
counteractive ways to regulate these mechanisms. Such aberrant complexes and
their regulating TFs play vital roles in the initiation and progression of
cancer.


  Taking the opportunity of the 20th anniversary of the word "proteomics", this
young adult age is a good time to remember how proteomics came from enormous
progress in protein separation and protein microanalysis techniques, and from
the conjugation of these advances into a high performance and streamlined
working setup. However, in the history of the almost three decades that
encompass the first attempts to perform large scale analysis of proteins to the
current high throughput proteomics that we can enjoy now, it is also
interesting to underline and to recall how difficult the first decade was.
Indeed when the word was cast, the battle was already won. This recollection is
mostly devoted to the almost forgotten period where proteomics was being
conceived and put to birth, as this collective scientific work will never
appear when searched through the keyword "proteomics". BIOLOGICAL SIGNIFICANCE:
The significance of this manuscript is to recall and review the two decades
that separated the first attempts of performing large scale analysis of
proteins from the solid technical corpus that existed when the word
"proteomics" was coined twenty years ago. This recollection is made within the
scientific historical context of this decade, which also saw the blossoming of
DNA cloning and sequencing. This article is part of a Special Issue entitled:
20 years of Proteomics in memory of Viatliano Pallini. Guest Editors: Luca Bini
, Juan J. Calvete, Natacha Turck, Denis Hochstrasser and Jean-Charles Sanchez.


  The radical pair model has been successful in explaining behavioral
characteristics of the geomagnetic compass believed to underlie the navigation
capability of certain avian species. In this study, the spin dynamics of the
radical pair model and decoherence therein are interpreted from a microscopic
state transition point of view. This helps to elucidate the interplay between
the hyperfine and Zeeman interactions that enables the avian compass, and the
distinctive effects of nuclear and environmental decoherence on it. Using a
quantum information theoretic quantifier of coherence, we find that nuclear
decoherence induces new structure in the spin dynamics without materially
affecting the compass action; environmental decoherence, on the other hand,
completely disrupts it.


  Self propelled particle (SPP) models are often compared with animal swarms.
However, the collective behaviour observed in experiments usually leaves
considerable unconstrained freedom in the structure of these models. To tackle
this degeneracy, and better distinguish between candidate models, we study
swarms of SPPs circulating in channels (like spins) where we permit information
to pass through windows between neighbouring channels. Co-alignment between
particles then couples the channels (antiferromagnetically) so that they tend
to counter-rotate. We study channels arranged to mimic a geometrically
frustrated antiferromagnet and show how the effects of this frustration allow
us to better distinguish between SPP models. Similar experiments could
therefore improve our understanding of collective motion in animals. Finally we
discuss how the spin analogy can be exploited to construct universal logic
gates and therefore swarming systems that can function as Turing machines.


  We categorise coherent band (aka nullification) pathways between knots and
2-component links. Additionally, we characterise the minimal coherent band
pathways (with intermediates) between any two knots or 2-component links with
small crossing number. We demonstrate these band surgeries for knots and links
with small crossing number. We apply these results to place lower bounds on the
minimum number of recombinant events separating DNA configurations, restrict
the recombination pathways and determine chirality and/or orientation of the
resulting recombinant DNA molecules.


  This paper compares the auto-repressed gene to a simple one (a gene without
auto-regulation) in terms of response time and output noise under the
assumption of fixed metabolic cost. The analysis shows that, in the case of
non-vanishing leak expression rate, the negative feedback reduces both the
switching on and switching off times of a gene. The noise of the auto-repressed
gene will be lower than the one of the simple gene only for low leak expression
rates. Summing up, for low, but non-vanishing leak expression rates, the
auto-repressed gene is both faster and less noisier compared to the simple one.


  Despite the fact that grouping behavior has been actively studied for over a
century, the relative importance of the numerous proposed fitness benefits of
grouping remain unclear. We use a digital model of evolving prey under
simulated predation to directly explore the evolution of gregarious foraging
behavior according to one such benefit, the "many eyes" hypothesis. According
to this hypothesis, collective vigilance allows prey in large groups to detect
predators more efficiently by making alarm signals or behavioral cues to each
other, thereby allowing individuals within the group to spend more time
foraging. Here, we find that collective vigilance is sufficient to select for
gregarious foraging behavior as long there is not a direct cost for grouping
(e.g., competition for limited food resources), even when controlling for
confounding factors such as the dilution effect. Further, we explore the role
of the genetic relatedness and reproductive strategy of the prey, and find that
highly related groups of prey with a semelparous reproductive strategy are the
most likely to evolve gregarious foraging behavior mediated by the benefit of
vigilance. These findings, combined with earlier studies with evolving digital
organisms, further sharpen our understanding of the factors favoring grouping
behavior.


  Cooperators that refuse to participate in sanctioning defectors create the
second-order free-rider problem. Such cooperators will not be punished because
they contribute to the public good, but they also eschew the costs associated
with punishing defectors. Altruistic punishers - those that cooperate and
punish - are at a disadvantage, and it is puzzling how such behaviour has
evolved. We show that sharing the responsibility to sanction defectors rather
than relying on certain individuals to do so permanently can solve the problem
of costly punishment. Inspired by the fact that humans have strong but also
emotional tendencies for fair play, we consider probabilistic sanctioning as
the simplest way of distributing the duty. In well-mixed populations the public
goods game is transformed into a coordination game with full cooperation and
defection as the two stable equilibria, while in structured populations pattern
formation supports additional counterintuitive solutions that are reminiscent
of Parrondo's paradox.


  We discuss a cancer hallmark network framework for modelling
genome-sequencing data to predict cancer clonal evolution and associated
clinical phenotypes. Strategies of using this framework in conjunction with
genome sequencing data in an attempt to predict personalized drug targets, drug
resistance, and metastasis for a cancer patient, as well as cancer risks for a
healthy individual are discussed. Accurate prediction of cancer clonal
evolution and clinical phenotypes will have substantial impact on timely
diagnosis, personalized management and prevention of cancer.


  Nanotechnology and synthetic biology currently constitute one of the most
innovative, interdisciplinary fields of research, poised to radically transform
society in the 21st century. This paper concerns the synthetic design of
ribonucleic acid molecules, using our recent algorithm, RNAiFold, which can
determine all RNA sequences whose minimum free energy secondary structure is a
user-specified target structure. Using RNAiFold, we design ten cis-cleaving
hammerhead ribozymes, all of which are shown to be functional by a cleavage
assay. We additionally use RNAiFold to design a functional cis-cleaving
hammerhead as a modular unit of a synthetic larger RNA. Analysis of kinetics on
this small set of hammerheads suggests that cleavage rate of computationally
designed ribozymes may be correlated with positional entropy, ensemble defect,
structural flexibility/rigidity and related measures. Artificial ribozymes have
been designed in the past either manually or by SELEX (Systematic Evolution of
Ligands by Exponential Enrichment); however, this appears to be the first
purely computational design and experimental validation of novel functional
ribozymes. RNAiFold is available at
http://bioinformatics.bc.edu/clotelab/RNAiFold/.


  All life, including cells and artificial protocells, must integrate diverse
molecules into a single unit in order to reproduce. Despite expected pressure
to evolve a simple system with the fastest replication speed, the mechanism by
which the use of the great variety of components, and the coexistence of
diverse cell-types with different compositions are achieved is as yet unknown.
Here we show that coexistence of such diverse compositions and cell-types is
the result of competitions for a variety of limited resources. We find that a
transition to diversity occurs both in chemical compositions and in protocell
types, as the resource supply is decreased, when the maximum inflow and
consumption of resources are balanced.


  It has been suggested that the origins of cognitive modernity in the
Middle/Upper Paleolithic following the appearance of anatomically modern humans
was due to the onset of dual processing or contextual focus (CF), the ability
to shift between different modes of thought: an explicit mode conducive to
logical problem solving, and an implicit mode conducive to free-association and
breaking out of a rut. Mathematical and computational models of CF supported
this hypothesis, showing that CF is conducive to making creative connections by
placing concepts in new contexts. This paper proposes that CF was made possible
by mutation of the FOXP2 gene in the Paleolithic. FOXP2, once thought to be the
'language gene', turned out not to be uniquely associated with language. In its
modern form FOXP2 enabled fine-tuning of the neurological mechanisms underlying
the capacity to shift between processing modes by varying the size of the
activated region of memory.


  Recent experimental results suggest that the E. coli chromosome feels a
self-attracting interaction of osmotic origin, and is condensed in foci by
bridging interactions. Motivated by these findings, we explore a generic
modeling framework combining solely these two ingredients, in order to
characterize their joint effects. Specifically, we study a simple polymer
physics computational model with weak ubiquitous short-ranged self attraction
and stronger sparse bridging interactions. Combining theoretical arguments and
simulations, we study the general phenomenology of polymer collapse induced by
these dual contributions, in the case of regularly-spaced bridging. Our results
distinguish a regime of classical Flory-like coil-globule collapse dictated by
the interplay of excluded volume and attractive energy and a switch-like
collapse where bridging interaction compete with entropy loss terms from the
looped arms of a star-like rosette. Additionally, we show that bridging can
induce stable compartmentalized domains. In these configurations, different
"cores" of bridging proteins are kept separated by star-like polymer loops in
an entropically favorable multi-domain configuration, with a mechanism that
parallels micellar polysoaps. Such compartmentalized domains are stable, and do
not need any intra-specific interactions driving their segregation. Domains can
be stable also in presence of uniform attraction, as long as the uniform
collapse is above its theta point.


  We study networks of biochemical reactions modelled by continuous-time Markov
processes. Such networks typically contain many molecular species and reactions
and are hard to study analytically as well as by simulation. Particularly, we
are interested in reaction networks with intermediate species such as the
substrate-enzyme complex in the Michaelis-Menten mechanism. These species are
virtually in all real-world networks, they are typically short-lived, degraded
at a fast rate and hard to observe experimentally.
  We provide conditions under which the Markov process of a multiscale reaction
network with intermediate species is approximated in finite dimensional
distribution by the Markov process of a simpler reduced reaction network
without intermediate species. We do so by embedding the Markov processes into a
one-parameter family of processes, where reaction rates and species abundances
are scaled in the parameter. Further, we show that there are close links
between these stochastic models and deterministic ODE models of the same
networks.


  Statistical machine learning methods, especially nonparametric Bayesian
methods, have become increasingly popular to infer clonal population structure
of tumors. Here we describe the treeCRP, an extension of the Chinese restaurant
process (CRP), a popular construction used in nonparametric mixture models, to
infer the phylogeny and genotype of major subclonal lineages represented in the
population of cancer cells. We also propose new split-merge updates tailored to
the subclonal reconstruction problem that improve the mixing time of Markov
chains. In comparisons with the tree-structured stick breaking prior used in
PhyloSub, we demonstrate superior mixing and running time using the treeCRP
with our new split-merge procedures. We also show that given the same number of
samples, TSSB and treeCRP have similar ability to recover the subclonal
structure of a tumor.


  Background
  Non-small cell lung cancer (NSCLC) is the predominant histological type of
lung cancer, accounting for up to 85% of cases. Disease stage is commonly used
to determine adjuvant treatment eligibility of NSCLC patients, however, it is
an imprecise predictor of the prognosis of an individual patient. Currently,
many researchers resort to microarray technology for identifying relevant
genetic prognostic markers, with particular attention on trimming or extending
a Cox regression model.
  Among NSCLC, adenocarcinoma (AC) and squamous cell carcinoma (SCC) are two
major histology subtypes. It has been demonstrated that there exist fundamental
differences in the underlying mechanisms between them, which motivated us to
postulate there might exist specific genes relevant to prognosis of each
histology subtype.
  Results
  In this article, we propose a simple filterer feature selection algorithm
with a Cox regression model as the base. Applying this method to a real-world
microarray data, no evidence has been found to support the existence of
histology-specific prognostic gene signature. Nevertheless, a 31-gene
prognostic gene signature for the early-stage AC and SCC samples is obtained,
which provides comparable performance when compared with other relevant
signatures.
  Conclusions
  Our proposal is conceptually simple and straightforward to implement.
Therefore, it is expected that other researchers, especially those with less
statistical knowledge and experience, can adapt this method readily to test
their own research hypotheses.


  In recent years the research community has accumulated overwhelming evidence
for the emergence of complex and heterogeneous connectivity patterns in a wide
range of biological and sociotechnical systems. The complex properties of
real-world networks have a profound impact on the behavior of equilibrium and
nonequilibrium phenomena occurring in various systems, and the study of
epidemic spreading is central to our understanding of the unfolding of
dynamical processes in complex networks. The theoretical analysis of epidemic
spreading in heterogeneous networks requires the development of novel
analytical frameworks, and it has produced results of conceptual and practical
relevance. A coherent and comprehensive review of the vast research activity
concerning epidemic processes is presented, detailing the successful
theoretical approaches as well as making their limits and assumptions clear.
Physicists, mathematicians, epidemiologists, computer, and social scientists
share a common interest in studying epidemic spreading and rely on similar
models for the description of the diffusion of pathogens, knowledge, and
innovation. For this reason, while focusing on the main results and the
paradigmatic models in infectious disease modeling, the major results
concerning generalized social contagion processes are also presented. Finally,
the research activity at the forefront in the study of epidemic spreading in
coevolving, coupled, and time-varying networks is reported.


  Recent experimental breakthroughs have finally allowed to implement in-vitro
reaction kinetics (the so called {\em enzyme based logic}) which code for
two-inputs logic gates and mimic the stochastic AND (and NAND) as well as the
stochastic OR (and NOR). This accomplishment, together with the already-known
single-input gates (performing as YES and NOT), provides a logic base and paves
the way to the development of powerful biotechnological devices. The
investigation of this field would enormously benefit from a self-consistent,
predictive, theoretical framework. Here we formulate a complete statistical
mechanical description of the Monod-Wyman-Changeaux allosteric model for both
single and double ligand systems, with the purpose of exploring their practical
capabilities to express logical operators and/or perform logical operations.
Mixing statistical mechanics with logics, and quantitatively our findings with
the available biochemical data, we successfully revise the concept of
cooperativity (and anti-cooperativity) for allosteric systems, with particular
emphasis on its computational capabilities, the related ranges and scaling of
the involved parameters and its differences with classical cooperativity (and
anti-cooperativity).


  Epistatic interactions between residues determine a protein's adaptability
and shape its evolutionary trajectory. When a protein experiences a changed
environment, it is under strong selection to find a peak in the new fitness
landscape. It has been shown that strong selection increases epistatic
interactions as well as the ruggedness of the fitness landscape, but little is
known about how the epistatic interactions change under selection in the
long-term evolution of a protein. Here we analyze the evolution of epistasis in
the protease of the human immunodeficiency virus type 1 (HIV-1) using protease
sequences collected for almost a decade from both treated and untreated
patients, to understand how epistasis changes and how those changes impact the
long-term evolvability of a protein. We use an information-theoretic proxy for
epistasis that quantifies the co-variation between sites, and show that
positive information is a necessary (but not sufficient) condition that detects
epistasis in most cases. We analyze the "fossils" of the evolutionary
trajectories of the protein contained in the sequence data, and show that
epistasis continues to enrich under strong selection, but not for proteins
whose environment is unchanged. The increase in epistasis compensates for the
information loss due to sequence variability brought about by treatment, and
facilitates adaptation in the increasingly rugged fitness landscape of
treatment. While epistasis is thought to enhance evolvability via
valley-crossing early-on in adaptation, it can hinder adaptation later when the
landscape has turned rugged. However, we find no evidence that the HIV-1
protease has reached its potential for evolution after 9 years of adapting to a
drug environment that itself is constantly changing.


  Entropy measures of probability distributions are widely used measures in
ecology, biology, genetics, and in other fields, to quantify species diversity
of a community. Unfortunately, entropy-based diversity indices, or diversity
indices for short, suffer from three problems. First, when computing the
diversity for samples withdrawn from communities with different structures,
diversity indices can easily yield non-comparable and hard to interpret
results. Second, diversity indices impose weighting schemes on the species
distributions that unnecessarily emphasize low abundant rare species, or
erroneously identified ones. Third, diversity indices do not allow for
comparing distributions against each other, which is necessary when a community
has a well-known species' distribution.
  In this paper we propose a new general methodology based on information
theoretic principles to quantify the species diversity of a community. Our
methodology, comprised of two steps, naturally overcomes the previous mentioned
problems, and yields comparable and easy to interpret diversity values. We show
that our methodology retains all the functional properties of any diversity
index, and yet is far more flexible than entropy--based diversity indices. Our
methodology is easy to implement and is applicable to any community of
interest.


  In cells and in vitro assays the number of motor proteins involved in
biological transport processes is far from being unlimited. The cytoskeletal
binding sites are in contact with the same finite reservoir of motors (either
the cytosol or the flow chamber) and hence compete for recruiting the available
motors, potentially depleting the reservoir and affecting cytoskeletal
transport. In this work we provide a theoretical framework to study,
analytically and numerically, how motor density profiles and crowding along
cytoskeletal filaments depend on the competition of motors for their binding
sites. We propose two models in which finite processive motor proteins actively
advance along cytoskeletal filaments and are continuously exchanged with the
motor pool. We first look at homogeneous reservoirs and then examine the
effects of free motor diffusion in the surrounding medium. We consider as a
reference situation recent in vitro experimental setups of kinesin-8 motors
binding and moving along microtubule filaments in a flow chamber. We
investigate how the crowding of linear motor proteins moving on a filament can
be regulated by the balance between supply (concentration of motor proteins in
the flow chamber) and demand (total number of polymerised tubulin
heterodimers). We present analytical results for the density profiles of bound
motors, the reservoir depletion, and propose novel phase diagrams that present
the formation of jams of motor proteins on the filament as a function of two
tuneable experimental parameters: the motor protein concentration and the
concentration of tubulins polymerized into cytoskeletal filaments. Extensive
numerical simulations corroborate the analytical results for parameters in the
experimental range and also address the effects of diffusion of motor proteins
in the reservoir.


  What can we learn from the collective dynamics of a complex network about its
interaction topology? Taking the perspective from nonlinear dynamics, we
briefly review recent progress on how to infer structural connectivity (direct
interactions) from accessing the dynamics of the units. Potential applications
range from interaction networks in physics, to chemical and metabolic
reactions, protein and gene regulatory networks as well as neural circuits in
biology and electric power grids or wireless sensor networks in engineering.
Moreover, we briefly mention some standard ways of inferring effective or
functional connectivity.


  It is a known fact that, given two rooted binary phylogenetic trees, the
concept of maximum acyclic agreement forests is sufficient to compute
hybridization networks with minimum hybridization number. In this work, we
demonstrate by first presenting an algorithm and then showing its correctness,
that this concept is also sufficient in the case of multiple input trees. More
precisely, we show that for computing minimum hybridization networks for
multiple rooted binary phylogenetic trees on the same set of taxa it suffices
to take only maximum acyclic agreement forests into account. Moreover, this
article contains a proof showing that the minimum hybridization number for a
set of rooted binary phylogenetic trees on the same set of taxa can be also
computed by solving subproblems referring to common clusters of the input
trees.


  Ion channels are of major interest and form an area of intensive research in
the fields of biophysics and medicine since they control many vital
physiological functions. The aim of this work is on one hand to propose a fully
stochastic and discrete model describing the main characteristics of a multiple
channel system. The movement of the ions is coupled, as usual, with a Poisson
equation for the electrical field; we have considered, in addition, the
influence of exclusion forces. On the other hand, we have discussed about the
nondimensionalization of the stochastic system by using real physical
parameters, all supported by numerical simulations. The specific features of
both cases of micro- and nanochannels have been taken in due consideration with
particular attention to the latter case in order to show that it is necessary
to consider a discrete and stochastic model for ions movement inside the
channels.


  Theoretical reasoning suggests that human cancer may result from knocking
down the genetic constraints evolved for maintenance of the metazoan
multicellularity, which, however, requires a critical test. Using
xenograft-based experimental evolution we characterized for the first time the
full life history from initiation to metastasis of a tumor at the genomic and
transcriptomic levels, and observed metastasis-driving positive selection for
generally loss-of-function mutations on a set of multicellularity-related
genes, which is further supported by large-scale exome data of clinical tumor
samples. Subsequent expression analysis revealed mainly expression
down-regulation of multicellularity-related genes, which form an evolving
expression profile approaching that of embryonic stem cells, the cell type with
the most characteristics of unicellular life. The theoretical conjecture
predicts that genes born at the emergence of metazoan multicellularity tend to
be cancer drivers, which we validated using a rigorous phylostratigraphy
analysis on the birth rate of genes annotated by Cancer Gene Census. Also, the
number of loss-of-function tumor suppressors often predominates over activated
oncogenes in a typical tumor of human patients. These data collectively suggest
that, different from typical organismal evolution in which gain of new genes is
the mainstream, cancer represents a loss-of-function-driven degenerative
evolution back to the unicellular ground state. This cancer evolution model may
explain the enormous tumoral genetic heterogeneity in the clinic, underlie how
distant-organ metastases originate in primary tumors despite distinct
environmental requirements, and hold implications for designing effective
cancer therapy.


  The sequence of a protein is not only constrained by its physical and
biochemical properties under current selection, but also by features of its
past evolutionary history. Understanding the extent and the form that these
evolutionary constraints may take is important to interpret the information in
protein sequences. To study this problem, we introduce a simple but physical
model of protein evolution where selection targets allostery, the functional
coupling of distal sites on protein surfaces. This model shows how the
geometrical organization of couplings between amino acids within a protein
structure can depend crucially on its evolutionary history. In particular, two
scenarios are found to generate a spatial concentration of functional
constraints: high mutation rates and fluctuating selective pressures. This
second scenario offers a plausible explanation for the high tolerance of
natural proteins to mutations and for the spatial organization of their least
tolerant amino acids, as revealed by sequence analyses and mutagenesis
experiments. It also implies a faculty to adapt to new selective pressures that
is consistent with observations. Besides, the model illustrates how several
independent functional modules may emerge within a same protein structure,
depending on the nature of past environmental fluctuations. Our model thus
relates the evolutionary history and evolutionary potential of proteins to the
geometry of their functional constraints, with implications for decoding and
engineering protein sequences.


  This paper is an attempt to formalize analytically the question raised in
"World Population Explained: Do Dead People Outnumber Living, Or Vice Versa?"
Huffington Post, \cite{HJ}. We start developing simple deterministic Malthusian
growth models of the problem (with birth and death rates either constant or
time-dependent) before running into both linear birth and death Markov chain
models and age-structured models.


  How do genes affect cognitive ability or other human quantitative traits such
as height or disease risk? Progress on this challenging question is likely to
be significant in the near future. I begin with a brief review of psychometric
measurements of intelligence, introducing the idea of a "general factor" or g
score. The main results concern the stability, validity (predictive power), and
heritability of adult g. The largest component of genetic variance for both
height and intelligence is additive (linear), leading to important
simplifications in predictive modeling and statistical estimation. Due mainly
to the rapidly decreasing cost of genotyping, it is possible that within the
coming decade researchers will identify loci which account for a significant
fraction of total g variation. In the case of height analogous efforts are well
under way. I describe some unpublished results concerning the genetic
architecture of height and cognitive ability, which suggest that roughly 10k
moderately rare causal variants of mostly negative effect are responsible for
normal population variation. Using results from Compressed Sensing
(L1-penalized regression), I estimate the statistical power required to
characterize both linear and nonlinear models for quantitative traits. The main
unknown parameter s (sparsity) is the number of loci which account for the bulk
of the genetic variation. The required sample size is of order 100s, or roughly
a million in the case of cognitive ability.


  Background: Influenza A/H3N2 has been circulating in humans since 1968,
causing considerable morbidity and mortality. Although H3N2 incidence is highly
seasonal, how such seasonality contributes to global phylogeographic migration
dynamics has not yet been established.
  Results: Incorporating seasonally varying migration rates improves the
modeling of migration. In our global model, windows of increased immigration
map to the seasonal timing of epidemic spread, while windows of increased
emigration map to epidemic decline. Seasonal patterns also correlate with the
probability that local lineages go extinct and fail to contribute to long term
viral evolution, as measured through the trunk of the phylogeny. However, the
fraction of the trunk in each community was found to be better determined by
its overall human population size
  Conclusions: Seasonal migration and rapid turnover within regions is
sustained by the invasion of 'fertile epidemic grounds' at the end of older
epidemics. Thus, the current emphasis on connectivity, including air-travel,
should be complemented with a better understanding of the conditions and timing
required for successful establishment.Models which account for migration
seasonality will improve our understanding of the seasonal drivers of
influenza,enhance epidemiological predictions, and ameliorate vaccine updating
by identifying strains that not only escape immunity but also have the seasonal
opportunity to establish and spread. Further work is also needed on additional
conditions that contribute to the persistence and long term evolution of
influenza within the human population,such as spatial heterogeneity with
respect to climate and seasonality


  The 2014 Ebola virus (EBOV) outbreak in West Africa is the largest outbreak
of the genus Ebolavirus to date. To better understand the spread of infection
in the affected countries, it is crucial to know the number of secondary cases
generated by an infected index case in the absence and presence of control
measures, i.e., the basic and effective reproduction number. In this study, I
describe the EBOV epidemic using an SEIR
(susceptible-exposed-infectious-recovered) model and fit the model to the most
recent reported data of infected cases and deaths in Guinea, Sierra Leone and
Liberia. The maximum likelihood estimates of the basic reproduction number are
1.51 (95% confidence interval [CI]: 1.50-1.52) for Guinea, 2.53 (95% CI:
2.41-2.67) for Sierra Leone and 1.59 (95% CI: 1.57-1.60) for Liberia. The model
indicates that in Guinea and Sierra Leone the effective reproduction number
might have dropped to around unity by the end of May and July 2014,
respectively. In Liberia, however, the model estimates no decline in the
effective reproduction number by end-August 2014. This suggests that control
efforts in Liberia need to be improved substantially in order to stop the
current outbreak.


  There is little doubt in scientific circles that--counting from the origin of
life towards today--evolution has led to an increase in the amount of
information stored within the genomes of the biosphere. This trend of
increasing information on average likely holds for every successful line of
descent, but it is not clear whether this increase is due to a general law, or
whether it is a secondary effect linked to an overall increase in fitness.
Here, we use "digital life" evolution experiments to study whether information
is under selection if treated as an organismal trait, using the Price equation.
By measuring both sides of the equation individually in an adapting population,
the strength of selection on a trait appears as a "gap" between the two terms
of the right-hand-side of the Price equation. We find that information is
strongly selected (as it encodes all fitness-producing traits) by comparing the
strength of selection on information to a weakly selected trait (sequence
length), as well as to a neutral marker. We observe that while strength of
selection on arbitrary traits can vary during an experiment (including
reversing sign), information is a selectable trait that must increase in a
fixed environment.


  Binding interactions between proteins and other molecules mediate numerous
cellular processes, including metabolism, signaling, and regulation of gene
expression. These interactions evolve in response to changes in the protein's
chemical or physical environment (such as the addition of an antibiotic), or
when genes duplicate and diverge. Several recent studies have shown the
importance of folding stability in constraining protein evolution. Here we
investigate how structural coupling between protein folding and binding -- the
fact that most proteins can only bind their targets when folded -- gives rise
to evolutionary coupling between the traits of folding stability and binding
strength. Using biophysical and evolutionary modeling, we show how these
protein traits can emerge as evolutionary "spandrels" even if they do not
confer an intrinsic fitness advantage. In particular, proteins can evolve
strong binding interactions that have no functional role but merely serve to
stabilize the protein if misfolding is deleterious. Furthermore, such proteins
may have divergent fates, evolving to bind or not bind their targets depending
on random mutation events. These observations may explain the abundance of
apparently nonfunctional interactions among proteins observed in
high-throughput assays. In contrast, for proteins with both functional binding
and deleterious misfolding, evolution may be highly predictable at the level of
biophysical traits: adaptive paths are tightly constrained to first gain extra
folding stability and then partially lose it as the new binding function is
developed. These findings have important consequences for our understanding of
fundamental evolutionary principles of both natural and engineered proteins.


  Intracellular fluctuations, mainly triggered by gene expression, are an
inevitable phenomenon observed in living cells. It influences generation of
phenotypic diversity in genetically identical cells. Such variation of cellular
components is beneficial in some contexts but detrimental in others. To
quantify the fluctuations in a gene product, we undertake an analytical scheme
for studying few naturally abundant linear as well as branched chain network
motifs. We solve the Langevin equations associated with each motif under the
purview of linear noise approximation and quantify Fano factor and mutual
information. Both quantifiable expressions exclusively depend on the relaxation
time (decay rate constant) and steady state population of the network
components. We investigate the effect of relaxation time constraints on Fano
factor and mutual information to indentify a time scale domain where a network
can recognize the fluctuations associated with the input signal more reliably.
We also show how input population affects both quantities. We extend our
calculation to long chain linear motif and show that with increasing chain
length, the Fano factor value increases but the mutual information processing
capability decreases. In this type of motif, the intermediate components are
shown to act as a noise filter that tune up input fluctuations and maintain
optimum fluctuations in the output. For branched chain motifs, both quantities
vary within a large scale due to their network architecture and facilitate
survival of living system in diverse environmental conditions.


  In this paper, we analyze the mean number $E(n,d)$ of internal equilibria in
a general $d$-player $n$-strategy evolutionary game where the agents' payoffs
are normally distributed. First, we give a computationally implementable
formula for the general case. Next we characterize the asymptotic behavior of
$E(2,d)$, estimating its lower and upper bounds as $d$ increases. Two important
consequences are obtained from this analysis. On the one hand, we show that in
both cases the probability of seeing the maximal possible number of equilibria
tends to zero when $d$ or $n$ respectively goes to infinity. On the other hand,
we demonstrate that the expected number of stable equilibria is bounded within
a certain interval. Finally, for larger $n$ and $d$, numerical results are
provided and discussed.


  We evaluate a version of the recently-proposed classification system named
Optimized Dissimilarity Space Embedding (ODSE) that operates in the input space
of sequences of generic objects. The ODSE system has been originally presented
as a classification system for patterns represented as labeled graphs. However,
since ODSE is founded on the dissimilarity space representation of the input
data, the classifier can be easily adapted to any input domain where it is
possible to define a meaningful dissimilarity measure. Here we demonstrate the
effectiveness of the ODSE classifier for sequences by considering an
application dealing with the recognition of the solubility degree of the
Escherichia coli proteome. Solubility, or analogously aggregation propensity,
is an important property of protein molecules, which is intimately related to
the mechanisms underlying the chemico-physical process of folding. Each protein
of our dataset is initially associated with a solubility degree and it is
represented as a sequence of symbols, denoting the 20 amino acid residues. The
herein obtained computational results, which we stress that have been achieved
with no context-dependent tuning of the ODSE system, confirm the validity and
generality of the ODSE-based approach for structured data classification.


  The dynamics of a local community of competing species with weak immigration
from a static regional pool is studied. Implementing the generalized
competitive Lotka-Volterra model with demographic noise, a rich dynamics
structure with four qualitatively distinct phases is unfolded. When the overall
interspecies competition is weak, the island species are a sample of the
mainland species. For higher values of the competition parameter the system
still admit an equilibrium community, but now some of the mainland species are
absent on the island. Further increase in competition leads to an intermittent
"chaotic" phase, where the dynamics is controlled by invadable combinations of
species and the turnover rate is governed by the migration. Finally, the strong
competition phase is glassy, dominated by uninvadable state and noise-induced
transitions. Our model contains, as a spatial case, the celebrated neutral
island theories of Wilson-MacArthur and Hubbell. Moreover, we show that slight
deviations from perfect neutrality may lead to each of the phases, as the
Hubbell point appears to be quadracritical.


  We present CloudSTRUCTURE, an application for running parallel analyses with
the population genetics program STRUCTURE. The HPC ready application, powered
by iPlant cyber-infrastructure, provides a fast (by parallelization) and
convenient (through a user friendly GUI) way to calculate like-lihood values
across multiple values of K (number of genetic groups) and numbers of
iterations. The results are automati-cally summarized for easier determination
of the K value that best fit the data. In addition, CloudSTRUCTURE will
reformat STRUCTURE output for use in downstream programs, such as TASSEL for
association analysis with population structure ef-fects stratified.


  Genetic interactions can strongly influence the fitness effects of individual
mutations, yet the impact of these epistatic interactions on evolutionary
dynamics remains poorly understood. Here we investigate the evolutionary role
of epistasis over 50,000 generations in a well-studied laboratory evolution
experiment in E. coli. The extensive duration of this experiment provides a
unique window into the effects of epistasis during long-term adaptation to a
constant environment. Guided by analytical results in the weak-mutation limit,
we develop a computational framework to assess the compatibility of a given
epistatic model with the observed patterns of fitness gain and mutation
accumulation through time. We find that a decelerating fitness trajectory alone
provides little power to distinguish between competing models, including those
that lack any direct epistatic interactions between mutations. However, when
combined with the mutation trajectory, these observables place strong
constraints on the set of possible models of epistasis, ruling out many
existing explanations of the data. Instead, we find that the data are
consistent with "two-epoch" model of adaptation, in which an initial burst of
diminishing returns epistasis is followed by a steady accumulation of mutations
under a constant distribution of fitness effects. Our results highlight the
need for additional DNA sequencing of these populations, as well as for more
sophisticated models of epistasis that are compatible with all of the
experimental data.


  An integro-differential model for evolutionary dynamics with mutations is
investigated by improving the understanding of its behavior using numerical
simulations. The proposed numerical approach can handle also density dependent
fitness, and gives new insights about the role of mutation in the preservation
of cooperation.


  16S-ribosomal sequencing and other metagonomic techniques provide snapshots
of microbial communities, revealing phylogeny and the abundances of microbial
populations across diverse ecosystems. While changes in microbial community
structure are demonstrably associated with certain environmental conditions,
identification of underlying mechanisms requires new statistical tools, as
these datasets present several technical challenges. First, the abundances of
microbial operational taxonomic units (OTUs) from 16S datasets are
compositional, and thus, microbial abundances are not independent. Secondly,
microbial sequencing-based studies typically measure hundreds of OTUs on only
tens to hundreds of samples; thus, inference of OTU-OTU interaction networks is
severely under-powered, and additional assumptions are required for accurate
inference. Here, we present SPIEC-EASI (SParse InversE Covariance Estimation
for Ecological Association Inference), a statistical method for the inference
of microbial ecological interactions from metagenomic datasets that addresses
both of these issues. SPIEC-EASI combines data transformations developed for
compositional data analysis with a graphical model inference framework that
assumes the underlying ecological interaction network is sparse. To reconstruct
the interaction network, SPIEC-EASI relies on algorithms for sparse
neighborhood and inverse covariance selection. Because no large-scale microbial
ecological networks have been experimentally validated, SPIEC-EASI comprises
computational tools to generate realistic OTU count data from a set of diverse
underlying network topologies. SPIEC-EASI outperforms state-of-the-art methods
in terms of edge recovery and network properties on realistic synthetic data
under a variety of scenarios. SPIEC-EASI also reproducibly predicts previously
unknown microbial interactions using data from the American Gut project.


  Diverse biological networks exhibit universal features distinguished from
those of random networks, calling much attention to their origins and
implications. Here we propose a minimal evolution model of Boolean regulatory
networks, which evolve by selectively rewiring links towards enhancing
adaptability to a changing environment and stability against dynamical
perturbations. We find that sparse and heterogeneous connectivity patterns
emerge, which show qualitative agreement with real transcriptional regulatory
networks and metabolic networks. The characteristic scaling behavior of
stability reflects the balance between robustness and flexibility. The scaling
of fluctuation in the perturbation spread shows a dynamic crossover, which is
analyzed by investigating separately the stochasticity of internal dynamics and
the network structures different depending on the evolution pathways. Our study
delineates how the ambivalent pressure of evolution shapes biological networks,
which can be helpful for studying general complex systems interacting with
environments.


  We study theoretically the replication of Kinetoplast DNA consisting of
several thousands separate mini-circles found in organisms of the class
Kinetoplastida. When the cell is not actively dividing these are topologically
connected in a marginally linked network of rings with only one connected
component. During cell division each mini-circle is removed from the network,
duplicated and then re-attached, along with its progeny. We study this process
under the hypothesis that there is a coupling between the topological state of
the mini-circles and the expression of genetic information encoded on them,
leading to the production of Topoisomerase. This model describes a
self-regulating system capable of full replication that reproduces several
previous experimental findings. We find that the fixed point of the system
depends on a primary free parameter of the model: the ratio between the rate of
removal of mini-circles from the network (R) and their (re)attachment rate (A).
The final topological state is found to be that of a marginally linked network
structure in which the fraction of mini-circles linked to the largest connected
component approaches unity as R/A decreases. Finally we discuss how this may
suggest an evolutionary trade-off between the speed of replication and the
accuracy with which a fully topologically linked state is produced.


  The effect of secondary structure on DNA duplex formation is poorly
understood. We use a coarse-grained model of DNA to show that specific 3- and
4-base pair hairpins reduce hybridization rates by factors of 2 and 10
respectively, in good agreement with experiment. By contrast, melting rates are
accelerated by factors of ~100 and ~2000. This surprisingly large speed-up
occurs because hairpins form during the melting process, stabilizing partially
melted states, and facilitating dissociation. These results may help guide the
design of DNA devices that use hairpins to modulate hybridization and
dissociation pathways and rates.


  We reformulate the eigenvalue problem for the selection--mutation equilibrium
distribution in the case of a haploid asexually reproduced population in the
form of an equation for an unknown probability generating function of this
distribution. The special form of this equation in the infinite sequence limit
allows us to obtain analytically the steady state distributions for a number of
particular cases of the fitness landscape. The general approach is illustrated
by examples; theoretical findings are compared with numerical calculations.


  A tandem cell is proposed for DNA sequencing in which an exonuclease enzyme
cleaves bases (mononucleotides) from a strand of DNA for identification inside
a nanopore. It has two nanopores and three compartments with the structure
[cis1, upstream nanopore (UNP), trans1=cis2, downstream nanopore (DNP),
trans2]. The exonuclease is attached to the exit side of UNP in trans1/cis2. A
cleaved base cannot regress into cis1 because of the remaining DNA strand in
UNP. A profiled electric field over DNP with positive and negative components
slows down base translocation through DNP. The proposed structure is modeled
with a Fokker-Planck equation and a piecewise solution presented. Results from
the model indicate that with probability approaching 1 bases enter DNP in their
natural order, are detected without any loss, and do not regress into DNP after
progressing into trans2. Sequencing efficiency with a tandem cell would then be
determined solely by the level of discrimination among the base types inside
DNP.


  Applications of first passage times in stochastic processes arise across a
wide range of length and time scales in biological settings. After an initial
technical overview, we survey representative applications and their
corresponding models. Within models that are effectively Markovian, we discuss
canonical examples of first passage problems spanning applications to molecular
dissociation and self-assembly, molecular search, transcription and
translation, neuronal spiking, cellular mutation and disease, and organismic
evolution and population dynamics. In this last application, a simple model for
stem-cell ageing is presented and some results derived. Various approximation
methods and the physical and mathematical subtleties that arise in the chosen
applications are also discussed.


  System-level properties of metabolic networks may be the direct product of
natural selection or arise as a by-product of selection on other properties.
Here we study the effect of direct selective pressure for growth or viability
in particular environments on two properties of metabolic networks: latent
versatility to function in additional environments and carbon usage efficiency.
Using a Markov Chain Monte Carlo (MCMC) sampling based on Flux Balance Analysis
(FBA), we sample from a known biochemical universe random viable metabolic
networks that differ in the number of directly constrained environments. We
find that the latent versatility of sampled metabolic networks increases with
the number of directly constrained environments and with the size of the
networks. We then show that the average carbon wastage of sampled metabolic
networks across the constrained environments decreases with the number of
directly constrained environments and with the size of the networks. Our work
expands the growing body of evidence about nonadaptive origins of key
functional properties of biological networks.


  We propose here to interpret and model peculiar plant morphologies (cushions,
tussocks) observed in the Andean altiplano as localized structures. Such
structures resulting in a patchy, aperiodic aspect of the vegetation cover are
hypothesized to self-organize thanks to the interplay between facilitation and
competition processes occurring at the scale of basic plant components
biologically referred to as 'ramets'. (Ramets are often of clonal origin.) To
verify this interpretation, we applied a simple, fairly generic model (one
integro-differential equation) emphasizing via Gaussian kernels non-local
facilitative and competitive feedbacks of the vegetation biomass density on its
own dynamics. We show that under realistic assumptions and parameter values
relating to ramet scale, the model can reproduce some macroscopic features of
the observed systems of patches and predict values for the inter-patch distance
that match the distances encountered in the reference area (Sajama National
Park in Bolivia). Prediction of the model can be confronted in the future to
data on vegetation patterns along environmental gradients as to anticipate the
possible effect of global change on those vegetation systems experiencing
constraining environmental conditions.


  Within the framework of population genetics we consider the evolution of an
asexual haploid population under the effect of a rapidly varying natural
selection (microevolution). We focus on the case in which the environment
exerting selection changes stochastically. We derive the effective genotype and
fitness dynamics on the slower time-scales at which the relevant genetic
modifications take place. We find that, despite the fast environmental
switches, the population manages to adapt on the fast time-scales yielding a
finite positive contribution to the fitness. However, such contribution is
balanced by the continuous loss in fitness due to the varying selection so that
the statistics of the global fitness can be described neglecting the details of
the fast environmental process. The occurrence of adaptation on fast
time-scales would be undetectable if one were to consider only the effective
genotype and fitness dynamics on the slow time-scales. We therefore propose an
experimental observable to detect it.


  According to the chemical reaction network theory, the topology of a certain
class of chemical reaction networks, regardless of the kinetic details, sets a
limit on the dynamical properties that a particular network can potentially
admit; the structure of a network predetermines the dynamic capacity of the
network. We note that stochastic fluctuations can possibly confer a new
dynamical capability to a network. Thus, it is of tremendous value to
understand and be able to control the landscape of stochastic dynamical
behaviors of a biochemical reaction network as a function of network
architecture. Here we investigate such a case where stochastic fluctuations can
give rise to the new capability of noise-induced oscillation in a subset of
biochemical reaction networks, the networks with only three biochemical species
whose reactions are governed by mass action kinetics and with the coupling of
positive and negative feedback loops. We model the networks with the master
equations and approximate them, using the linear noise approximation. For each
network, we read the signal-to-noise ratio value, an indicator of amplified and
coherent noise-induced oscillation, off from the analytically derived power
spectra. We classify the networks into three performance groups based on the
average values of the signal-to-noise ratio and the robustness. We identify the
common network architecture among the networks belonging to the same
performance group, from which we learn that the coupling of negative and
positive feedback loops generally enhance the noise-induced oscillation
performance better than the negative feedback loops alone. The performance of
networks also depends on the relative size of the positive and negative
feedback loops; the networks with the bigger positive and smaller negative
feedbacks are much worse oscillators than the networks with only negative
feedback loops.


  The twin-arginine translocation (Tat) system transports folded proteins of
various sizes across both bacterial and plant thylakoid membranes. The
membrane-associated TatA protein is an essential component of the Tat
translocon, and a broad distribution of different sized TatA-clusters is
observed in bacterial membranes. We assume that the size dynamics of TatA
clusters are affected by substrate binding, unbinding, and translocation to
associated TatBC clusters, where clusters with bound translocation substrates
favour growth and those without associated substrates favour shrinkage. With a
stochastic model of substrate binding and cluster dynamics, we numerically
determine the TatA cluster size distribution. We include a proportion of
targeted but non-translocatable (NT) substrates, with the simplifying
hypothesis that the substrate translocatability does not directly affect
cluster dynamical rate constants or substrate binding or unbinding rates. This
amounts to a translocation model without specific quality control.
Nevertheless, NT substrates will remain associated with TatA clusters until
unbound and so will affect cluster sizes and translocation rates. We find that
the number of larger TatA clusters depends on the NT fraction $f$. The
translocation rate can be optimized by tuning the rate of spontaneous substrate
unbinding, $\Gamma_U$. We present an analytically solvable three-state model of
substrate translocation without cluster size dynamics that follows our computed
translocation rates, and that is consistent with {\em in vitro}
Tat-translocation data in the presence of NT substrates.


  We study biological evolution on a random fitness landscape where
correlations are introduced through a linear fitness gradient of strength $c$.
When selection is strong and mutations rare the dynamics is a directed uphill
walk that terminates at a local fitness maximum. We analytically calculate the
dependence of the walk length on the genome size $L$. When the distribution of
the random fitness component has an exponential tail we find a phase transition
of the walk length $D$ between a phase at small $c$ where walks are short $(D
\sim \ln L)$ and a phase at large $c$ where walks are long $(D \sim L)$. For
all other distributions only a single phase exists for any $c > 0$. The
considered process is equivalent to a zero temperature Metropolis dynamics for
the random energy model in an external magnetic field, thus also providing
insight into the aging dynamics of spin glasses.


  Phenotypic heterogeneity is a strategy commonly used by bacteria to rapidly
adapt to changing environmental conditions. Here, we study the interplay
between phenotypic heterogeneity and genetic diversity in spatially extended
populations. By analyzing the spatio-temporal dynamics, we show that the level
of mobility and the type of competition qualitatively influence the persistence
of phenotypic heterogeneity. While direct competition generally promotes
persistence of phenotypic heterogeneity, specialization dominates in models
with indirect competition irrespective of the degree of mobility.


  Cytoplasmic dynein exhibits a directional processive movement on microtubule
filaments and is known to move in steps of varying length based on the number
of ATP molecules bound to it and the load that it carries. It is experimentally
observed that dynein takes occasional backward steps and the frequency of such
backward steps increases as the load approaches the stall force. Using a
stochastic process model, we investigate the bidirectional movement of single
head of a dynein motor. The probability for backward step is implemented based
on Crook's fluctuation theorem of non-equilibrium statistical mechanics. We
find that the movement of dynein motor is characterized with negative velocity
implying backward motion beyond stall force. We observe that the motor moves
backward for super stall forces by hydrolyzing the ATP exactly the same way as
it does while moving forward for sub stall forces.


  There is increased interest in the identification and analysis of gene
fusions and chimeric RNA transcripts. While most recent efforts focused on the
analysis of genomic and transcriptomic data, identi-fication of novel peptides
corresponding to such events in mass spectrometry-based proteomic datasets
would provide complemen-tary, protein-level evidence. The process of
identifying fusion pro-teins from mass spectrometry data is inherently
difficult because such events are rare. It is also complicated due to large
amount of spectra collected and the explosion in the number of candidate
peptide sequences that need to be considered, which makes ex-haustive search
for all possible fusion partner proteins impractical. In this work, we present
a sequence tag based fusion protein identi-fication algorithm, FIT, that
combines the virtue of both de novo sequence tag retrieval and peptide-spectrum
matching for identifi-cation of fusion proteins. Results on simulated datasets
show high sensitivity and low false positive rates for fusion protein
identifica-tion by the FIT algorithm.


  We investigate the statistics of the time taken for a system driven by
recruitment to reach fixation. Our model describes a series of experiments
where a population is confronted with two identical options, resulting in the
system fixating on one of the options. For a specific population size, we show
that the time distribution behaves like an inverse Gaussian with an exponential
decay. Varying the population size reveals that the timescale of the decay
depends on the population size and allows the critical population number, below
which fixation occurs, to be estimated from experimental data.


  We consider a branching particle system where particles reproduce according
to the pure birth Yule process with the birth rate L, conditioned on the
observed number of particles to be equal n. Particles are assumed to move
independently on the real line according to the Brownian motion with the local
variance s2. In this paper we treat $n$ particles as a sample of related
species. The spatial Brownian motion of a particle describes the development of
a trait value of interest (e.g. log-body-size). We propose an unbiased
estimator Rn2 of the evolutionary rate r2=s2/L. The estimator Rn2 is
proportional to the sample variance Sn2 computed from n trait values. We find
an approximate formula for the standard error of Rn2 based on a neat asymptotic
relation for the variance of Sn2.


  We propose a new model of permanent monogamous pair formation in zoological
populations with multiple types of females and males. According to this model,
animals randomly encounter members of the opposite sex at their so-called
firing times to form temporary pairs which then become permanent if mating
happens. Given the distributions of the firing times and the mating preferences
upon encounter, we analyze the contingency table of permanent pair types in
three cases: (i) definite mating upon encounter; (ii) Poisson firing times; and
(iii) Bernoulli firing times. In the first case, the contingency table has a
multiple hypergeometric distribution which implies panmixia. The other two
cases generalize the encounter-mating models of Gimelfarb (1988) who gives
conditions that he conjectures to be sufficient for panmixia. We formulate
adaptations of his conditions and prove that they not only characterize
panmixia but also allow us to reduce the model to the first case by changing
its underlying parameters. Finally, when there are only two types of females
and males, we provide a full characterization of panmixia, homogamy and
heterogamy.


  Many methods have been developed for finding the commonalities between
different organisms to study their phylogeny. The structure of metabolic
networks also reveal valuable insights into metabolic capacity of species as
well as into the habitats where they have evolved. We constructed metabolic
networks of 79 fully sequenced organisms and compared their architectures. We
used spectral density of normalized Laplacian matrix for comparing the
structure of networks. The eigenvalues of this matrix reflect not only the
global architecture of a network but also the local topologies that are
produced by different graph evolutionary processes like motif duplication or
joining. A divergence measure on spectral densities is used to quantify the
distances between various metabolic networks, and a split network is
constructed to analyze the phylogeny from these distances. In our analysis, we
focus on the species, which belong to different classes, but appear more
related to each other in the phylogeny. We tried to explore whether they have
evolved under similar environmental conditions or have similar life histories.
With this focus, we have obtained interesting insights into the phylogenetic
commonality between different organisms.


  The Matthew effect describes the phenomenon that in societies the rich tend
to get richer and the potent even more powerful. It is closely related to the
concept of preferential attachment in network science, where the more connected
nodes are destined to acquire many more links in the future than the auxiliary
nodes. Cumulative advantage and success-breads-success also both describe the
fact that advantage tends to beget further advantage. The concept is behind the
many power laws and scaling behaviour in empirical data, and it is at the heart
of self-organization across social and natural sciences. Here we review the
methodology for measuring preferential attachment in empirical data, as well as
the observations of the Matthew effect in patterns of scientific collaboration,
socio-technical and biological networks, the propagation of citations, the
emergence of scientific progress and impact, career longevity, the evolution of
common English words and phrases, as well as in education and brain
development. We also discuss whether the Matthew effect is due to chance or
optimisation, for example related to homophily in social systems or efficacy in
technological systems, and we outline possible directions for future research.


  Repeated game theory has been one of the most prevailing tools for
understanding the long-run relationships, which are footstones in building
human society. Recent works have revealed a new set of "zero-determinant (ZD)"
strategies, which is an important advance in repeated games. A ZD strategy
player can exert a unilaterally control on two players' payoffs. In particular
he can deterministically set the opponent's payoff, or enforce an unfair linear
relationship between the players' payoffs, thereby always seizing an
advantageous share of payoffs. One of the limitations of the original ZD
strategy, however, is that it does not capture the notion of robustness when
the game is subjected to stochastic errors. In this paper, we propose a general
model of ZD strategies for noisy repeated games, and find that ZD strategies
have high robustness against errors. We further derive the pinning strategy
under noise, by which the ZD strategy player coercively set the opponent's
expected payoff to his desired level, although his payoff control ability
declines with the increase of noise strength. Due to the uncertainty caused by
noise, the ZD strategy player cannot secure his payoff to be higher than the
opponent's, which implies strong extortions do not exist even under low noise.
While we show that the ZD strategy player can still establish a novel kind of
extortions, named weak extortions, where any increase of his own payoff always
exceeds that of the opponent's by a fixed percentage, and the conditions under
which the weak extortions can be realized are more stringent as the noise
becomes stronger.