  When a self-localized quasiparticle excitation propagates along a discrete
one dimensional lattice, it becomes subject to a dissipation that converts the
kinetic energy into lattice vibrations. Eventually the kinetic energy does no
longer enable the excitation to cross over the minimum energy barrier between
neighboring sites, and the excitation becomes localized within a lattice cell.
In the case of a protein, the lattice structure consists of the C-alpha
backbone. The self-localized quasiparticle excitation is the elemental building
block of loops. It can be modeled by a kink which solves a variant of the
discrete non-linear Schroedinger equation (DNLS). We study the propagation of
such a kink in the case of protein G related albumin-binding domain, using the
UNRES coarse-grained molecular dynamics force field. We estimate the height of
the energy barriers the kink needs to cross over, in order to propagate along
the backbone lattice. We analyse how these barriers gives rise to both stresses
and reliefs which control the kink movement. For this, we deform a natively
folded protein structure by parallel translating the kink along the backbone
away from its native position. We release the transposed kink, and we follow
how it propagates along the backbone towards the native location. We observe
that the dissipative forces which are exerted on the kink by the various energy
barriers, have a pivotal role in determining how a protein folds towards its
native state.


  The behaviour of gene regulatory networks (GRNs) is typically analysed using
simulation-based statistical testing-like methods. In this paper, we
demonstrate that we can replace this approach by a formal verification-like
method that gives higher assurance and scalability. We focus on Wagner weighted
GRN model with varying weights, which is used in evolutionary biology. In the
model, weight parameters represent the gene interaction strength that may
change due to genetic mutations. For a property of interest, we synthesise the
constraints over the parameter space that represent the set of GRNs satisfying
the property. We experimentally show that our parameter synthesis procedure
computes the mutational robustness of GRNs -an important problem of interest in
evolutionary biology- more efficiently than the classical simulation method. We
specify the property in linear temporal logics. We employ symbolic bounded
model checking and SMT solving to compute the space of GRNs that satisfy the
property, which amounts to synthesizing a set of linear constraints on the
weights.


  Infectious diseases are one of the leading causes of morbidity and mortality
around the world; thus, forecasting their impact is crucial for planning an
effective response strategy. According to the Centers for Disease Control and
Prevention (CDC), seasonal influenza affects between 5% to 20% of the U.S.
population and causes major economic impacts resulting from hospitalization and
absenteeism. Understanding influenza dynamics and forecasting its impact is
fundamental for developing prevention and mitigation strategies.
  We combine modern data assimilation methods with Wikipedia access logs and
CDC influenza like illness (ILI) reports to create a weekly forecast for
seasonal influenza. The methods are applied to the 2013--2014 influenza season
but are sufficiently general to forecast any disease outbreak, given incidence
or case count data. We adjust the initialization and parametrization of a
disease model and show that this allows us to determine systematic model bias.
In addition, we provide a way to determine where the model diverges from
observation and evaluate forecast accuracy.
  Wikipedia article access logs are shown to be highly correlated with
historical ILI records and allow for accurate prediction of ILI data several
weeks before it becomes available. The results show that prior to the peak of
the flu season, our forecasting method projected the actual outcome with a high
probability. However, since our model does not account for re-infection or
multiple strains of influenza, the tail of the epidemic is not predicted well
after the peak of flu season has past.


  Recent studies have shown that single-stranded viral RNAs fold into more
compact structures than random RNA sequences with similar chemical composition
and identical length. Based on this comparison it has been suggested that
wild-type viral RNA may have evolved to be atypically compact so as to aid its
encapsidation and assist the viral assembly process. In order to further
explore the compactness selection hypothesis, we systematically compare the
predicted sizes of more than one hundred wild-type viral sequences with those
of their mutants, which are evolved in silico and subject to a number of known
evolutionary constraints. In particular, we enforce mutation synonynimity,
preserve the codon-bias, and leave untranslated regions intact. It is found
that progressive accumulation of these restricted mutations still suffices to
completely erase the characteristic compactness imprint of the viral RNA
genomes, making them in this respect physically indistinguishable from randomly
shuffled RNAs. This shows that maintaining the physical compactness of the
genome is indeed a primary factor among ssRNA viruses evolutionary constraints,
contributing also to the evidence that synonymous mutations in viral ssRNA
genomes are not strictly neutral.


  Liquid chromatography coupled with tandem mass spectrometry, also known as
shotgun proteomics, is a widely-used high-throughput technology for identifying
proteins in complex biological samples. Analysis of the tens of thousands of
fragmentation spectra produced by a typical shotgun proteomics experiment
begins by assigning to each observed spectrum the peptide hypothesized to be
responsible for generating the spectrum, typically done by searching each
spectrum against a database of peptides. We have recently described a machine
learning method---Dynamic Bayesian Network for Rapid Identification of Peptides
(DRIP)---that not only achieves state-of-the-art spectrum identification
performance on a variety of datasets but also provides a trainable model
capable of returning valuable auxiliary information regarding specific
peptide-spectrum matches. In this work, we present two significant improvements
to DRIP. First, we describe how to use word lattices, which are widely used in
natural language processing, to significantly speed up DRIP's computations. To
our knowledge, all existing shotgun proteomics search engines compute
independent scores between a given observed spectrum and each possible
candidate peptide from the database. The key idea of the word lattice is to
represent the set of candidate peptides in a single data structure, thereby
allowing sharing of redundant computations among the different candidates. We
demonstrate that using lattices in conjunction with DRIP leads to speedups on
the order of tens across yeast and worm data sets. Second, we introduce a
variant of DRIP that uses a discriminative training framework, performing
maximum mutual entropy estimation rather than maximum likelihood estimation.
This modification improves DRIP's statistical power, enabling us to increase
the number of identified spectrum at a 1% false discovery rate on yeast and
worm data sets.


  Something about structure of the water shell of DNA and how it can help to
design of novel biologically-active molecules and potential drugs with
sequence-specific binding to nucleic acids.


  Terraces are potentially large sets of trees with precisely the same
likelihood or parsimony score, which can be induced by missing sequences in
partitioned multi-locus phylogenetic data matrices. The set of trees on a
terrace can be characterized by enumeration algorithms or consensus methods
that exploit the pattern of partial taxon coverage in the data, independent of
the sequence data themselves. Terraces add ambiguity and complexity to
phylogenetic inference particularly in settings where inference is already
challenging: data sets with many taxa and relatively few loci. In this paper we
present five new findings about terraces and their impacts on phylogenetic
inference. First we clarify assumptions about model parameters that are
necessary for the existence of terraces. Second, we explore the dependence of
terrace size on partitioning scheme and indicate how to find the partitioning
scheme associated with the largest terrace containing a given tree. Third, we
highlight the impact of terraces on bootstrap estimates of confidence limits in
clades, and characterize the surprising result that the bootstrap proportion
for a clade can be entirely determined by the frequency of bipartitions on a
terrace, with some bipartitions receiving high support even when incorrect.
Fourth, we dissect some effects of prior distributions of edge lengths on the
computed posterior probabilities of clades on terraces, to understand an
example in which long edges "attract" each other in Bayesian inference. Fifth,
we show that even if data are not partitioned, patterns of missing data studied
in the terrace problem can lead to instances of apparent statistical
inconsistency when even a small element of heterotachy is introduced to the
model generating the sequence data. Finally, we discuss strategies for
remediation of some of these problems.


  The development and evolution of malware including computer viruses, worms,
and trojan horses, is shown to be closely analogous to the process of community
succession long recognized in ecology. In particular, both changes in the
overall environment by external disturbances, as well as, feedback effects from
malware competition and antivirus coevolution have driven community succession
and the development of different types of malware with varying modes of
transmission and adaptability.


  We propose a general parametrizable model to capture the dynamic interaction
among bacteria in the formation of micro-colonies. micro-colonies represent the
first social step towards the formation of structured multicellular communities
known as bacterial biofilms, which protect the bacteria against antimicrobials.
In our model, bacteria can form links in the form of intercellular adhesins
(such as polysaccharides) to collaborate in the production of resources that
are fundamental to protect them against antimicrobials. Since maintaining a
link can be costly, we assume that each bacterium forms and maintains a link
only if the benefit received from the link is larger than the cost, and we
formalize the interaction among bacteria as a dynamic network formation game.
We rigorously characterize some of the key properties of the network evolution
depending on the parameters of the system. In particular, we derive the
parameters under which it is guaranteed that all bacteria will join
micro-colonies and the parameters under which it is guaranteed that some
bacteria will not join micro-colonies. Importantly, our study does not only
characterize the properties of networks emerging in equilibrium, but it also
provides important insights on how the network dynamically evolves and on how
the formation history impacts the emerging networks in equilibrium. This
analysis can be used to develop methods to influence on- the-fly the evolution
of the network, and such methods can be useful to treat or prevent
biofilm-related diseases.


  West Africa is currently experiencing a severe outbreak of Ebola virus
disease (EVD). As part of the international effort to address this outbreak,
the United States has committed to building specialized Ebola treatment
facilities with 1700 beds. However, the effectiveness of this increase in the
available healthcare facilities to treat Ebola is unclear, especially in light
of the rapidly increasing number of cases. Adapting a previously validated
mathematical model of Ebola in West Africa, we examine the potential impact of
an increase in hospital capacity to mitigate the impact of Ebola under several
scenarios, ranging from the planned scenario of 1700 beds in 10 weeks to a
considerably more aggressive approach of twice the number of beds in 5 weeks.
We find that even for the most aggressive scenarios, while increasing the
availability of healthcare reduces the number of Ebola cases and slows the
outbreak, it is not sufficient to stop the epidemic within the next three
months. We find that only a combination of increased hospital beds and a
dramatic decrease in the rate of transmission within the community can bring
the epidemic under control within the near future.


  A new method is proposed to infer unobserved epidemic sub-populations by
exploiting the synchronization properties of multistrain epidemic models. A
model for dengue fever is driven by simulated data from secondary infective
populations. Primary infective populations in the driven system synchronize to
the correct values from the driver system. Most hospital cases of dengue are
secondary infections, so this method provides a way to deduce unobserved
primary infection levels. We derive center manifold equations that relate the
driven system to the driver system and thus motivate the use of synchronization
to predict unobserved primary infectives. Synchronization stability between
primary and secondary infections is demonstrated through numerical measurements
of conditional Lyapunov exponents and through time series simulations.


  Random mixing in host populations has been a convenient simplifying
assumption in the study of epidemics, but neglects important differences in
contact rates within and between population groups. For HIV/AIDS, the
assumption of random mixing is inappropriate for epidemics that are
concentrated in groups of people at high risk, including female sex workers
(FSW) and their male clients (MCF), injection drug users (IDU) and men who have
sex with men (MSM). To find out who transmits infection to whom and how that
affects the spread and containment of infection remains a major empirical
challenge in the epidemiology of HIV/AIDS. Here we develop a technique, based
on the routine sampling of infection in linked population groups, which shows
how an Asian HIV/AIDS epidemic began in FSW, was propagated mainly by IDU, and
ultimately generated most cases among the female partners of MCF (FPM).
Calculation of the case reproduction numbers within and between groups, and for
the whole network, provides insights into control that cannot be deduced simply
from observations on the prevalence of infection. Specifically, the per capita
rate of HIV transmission was highest from FSW to MCF, and most HIV infections
occurred in FPM, but the number of infections in the whole network is best
reduced by interrupting transmission to and from IDU. This network analysis can
be used to guide HIV/AIDS interventions based on needle exchange, condom
distribution and antiretroviral therapy. The method requires only routine data
and could be applied to infections in other populations.


  In 2014, a major epidemic of human Ebola virus disease emerged in West
Africa, where human-to-human transmission has now been been sustained for
greater than 10 months. In the summer of 2014, there was great uncertainty
about the answers to several key policy questions concerning the path to
containment. In recent years, epidemic models have been used to guide public
health interventions. But, model-based policy relies on high quality causal
understanding of transmission, including the availability of appropriate
dynamic transmission models and reliable reporting about the sequence of case
incidence for model fitting, which were lacking for this epidemic. To
investigate the range of potential transmission scenarios, we developed a
multi-type branching process model that incorporates key heterogeneities and
time-varying parameters to reflect changing human behavior and deliberate
interventions. Ensembles of this model were evaluated at a set of parameters
that were both epidemiologically plausible and capable of reproducing the
observed trajectory. Results suggest that epidemic outcome depends on both
hospital capacity and individual behavior. The model predicts that if hospital
capacity is not increased soon, then transmission may outpace the rate of
isolation and the ability to provide care for the ill, infectious, and dying.
Similarly, containment will probably require individuals to adopt behaviors
that increase the rates of case identification and isolation and secure burial
of the deceased. Given current knowledge, it is uncertain that this epidemic
will be contained even with 99% hospitalization rate at the currently projected
hospital capacity.


  Correlations between two variables of a high-dimensional system can be
indicative of an underlying interaction, but can also result from indirect
effects. Inverse Ising inference is a method to distinguish one from the other.
Essentially, the parameters of the least constrained statistical model are
learned from the observed correlations such that direct interactions can be
separated from indirect correlations. Among many other applications, this
approach has been helpful for protein structure prediction, because residues
which interact in the 3D structure often show correlated substitutions in a
multiple sequence alignment. In this context, samples used for inference are
not independent but share an evolutionary history on a phylogenetic tree. Here,
we discuss the effects of correlations between samples on global inference.
Such correlations could arise due to phylogeny but also via other slow
dynamical processes. We present a simple analytical model to address the
resulting inference biases, and develop an exact method accounting for
background correlations in alignment data by combining phylogenetic modeling
with an adaptive cluster expansion algorithm. We find that popular reweighting
schemes are only marginally effective at removing phylogenetic bias, suggest a
rescaling strategy that yields better results, and provide evidence that our
conclusions carry over to the frequently used mean-field approach to the
inverse Ising problem.


  Our planet is experiencing an accelerated process of change associated to a
variety of anthropogenic phenomena. The future of this transformation is
uncertain, but there is general agreement about its negative unfolding that
might threaten our own survival. Furthermore, the pace of the expected changes
is likely to be abrupt: catastrophic shifts might be the most likely outcome of
this ongoing, apparently slow process. Although different strategies for
geo-engineering the planet have been advanced, none seem likely to safely
revert the large-scale problems associated to carbon dioxide accumulation or
ecosystem degradation. An alternative possibility considered here is inspired
in the rapidly growing potential for engineering living systems. It would
involve designing synthetic organisms capable of reproducing and expanding to
large geographic scales with the goal of achieving a long-term or a transient
restoration of ecosystem-level homeostasis. Such a regional or even
planetary-scale engineering would have to deal with the complexity of our
biosphere. It will require not only a proper design of organisms but also
understanding their place within ecological networks and their evolvability.
This is a likely future scenario that will require integration of ideas coming
from currently weakly connected domains, including synthetic biology,
ecological and genome engineering, evolutionary theory, climate science,
biogeography and invasion ecology, among others.


  Genes are connected in regulatory networks, often modelled by ordinary
differential equations. Changes in expression of a gene propagate to other
genes along paths in the network. At a stable state, the system's Jacobian
matrix confers information about network connectivity. To disclose the
functional properties of genes, knowledge of network connections is essential.
We present a new method to reconstruct the Jacobian matrix of models for gene
regulatory systems from equilibrium protein concentrations. In a recent paper
we defined propagation and feedback functions describing how genetic variation
at one gene propagates to the other genes in the network and possibly also back
to itself. Here we show how propagation and feedback functions provide
relations between equilibrium protein levels which are in principle observable,
and Jacobi elements which are not directly observable. We establish exact
formulae expressing the Jacobian in terms of derivatives of propagation and
feedback functions. Approximating these derivatives from perturbed and
unperturbed protein levels, we derive formulae for estimating the Jacobian. We
apply the method to models of the Drosophila segment polarity network and
randomly generated gene networks. Genes could be perturbed in two ways: by
modifying mRNA degradation rates, or by allele knockout in diploid models.
Comparison with the true Jacobians shows that for noiseless data we obtain hit
rates close to 100% in the former case and in the range 80-90% in the latter.
Our method adds to the network interference toolbox and provides a sign
estimate of the Jacobian from steady state data, and a value estimate of the
Jacobian if protein degradation rates are known. Also the approach identifies
some predicted connections as much more reliable than others, and could point
to further experiments for resolving uncertainties in the less dependable
Jacobian elements.


  The Gibbs free energy of transferring a solute at infinite dilution between
two solvents quantifies differences in solute-solvent interactions $-$ if the
transfer takes place at constant molarity of the solute. Yet, many calculation
formulae and measuring instructions that are commonly used to quantify
solute-solvent interactions correspond to transfer processes in which not the
molarity of the solute but its concentration measured in another concentration
scale is constant. Here, we demonstrate that in this case, not only the change
in solute-solvent interactions is quantified but also the entropic effect of a
volume change during the transfer. Consequently, the "phenomenon" which is
known as "concentration-scale dependence" of transfer free energies is simply
explained by a volume-entropy effect. Our explanations are of high importance
for the study of cosolvent effects on protein stability.


  We consider the problem of extinction processes on random networks with a
given structure. For sufficiently large well-mixed populations, the process of
extinction of one or more state variable components occurs in the tail of the
quasi-stationary probability distribution, thereby making it a rare event. Here
we show how to extend the theory of large deviations to random networks to
predict extinction times. In particular, we use the theory to find the most
probable path leading to extinction. We apply the methodology to epidemic
models and discover how mean extinction times scale with epidemiological and
network parameters in Erdos-Renyi networks. The results are shown to compare
quite well with Monte Carlo simulations of the network in predicting both the
most probable paths to extinction and mean extinction times.


  Continuous-time birth-death-shift (BDS) processes are frequently used in
stochastic modeling, with many applications in ecology and epidemiology. In
particular, such processes can model evolutionary dynamics of transposable
elements - important genetic markers in molecular epidemiology. Estimation of
the effects of individual covariates on the birth, death, and shift rates of
the process can be accomplished by analyzing patient data, but inferring these
rates in a discretely and unevenly observed setting presents computational
challenges. We propose a mutli-type branching process approximation to BDS
processes and develop a corresponding expectation maximization (EM) algorithm,
where we use spectral techniques to reduce calculation of expected sufficient
statistics to low dimensional integration. These techniques yield an efficient
and robust optimization routine for inferring the rates of the BDS process, and
apply more broadly to multi-type branching processes where rates can depend on
many covariates. After rigorously testing our methodology in simulation
studies, we apply our method to study intrapatient time evolution of IS6110
transposable element, a frequently used element during estimation of
epidemiological clusters of Mycobacterium tuberculosis infections.


  A general procedure to formulate asexual (unstructured, deterministic)
population dynamical models resulting from individual pairwise interactions is
proposed. Individuals are characterized by a continuous strategy that
represents all their behavioral, morphological, and functional traits.
Populations group conspecific individuals with identical strategy and are
measured by densities in space. Species can be monomorphic, if only one
strategy is present, or polymorphic otherwise. The procedure highlights the
structural properties fulfilled by the population per-capita growth rates. In
particular, the effect of perturbing a set of similar strategies is
proportional to the product of the corresponding densities, with a
proportionality coefficient that is density-dependent only through the total
density. This generalizes the law of mass action, which traditionally refers to
the case in which the per-capita growth rates are linearly density-dependent
and insensitive to joint strategy perturbations. Being underpinned with
individual strategies, the proposed procedure is most useful for evolutionary
considerations, in the case strategies are inheritable. The developed body of
theory is exemplified on a Holling-type-II many-prey-one-predator system and on
a model of a cannibalistic community.


  Background: Diversity estimates in cultivated plants provide a rationale for
conservation strategies and support the selection of starting material for
breeding programs. Diversity measures applied to crops usually have been
limited to the assessment of genome polymorphism at the DNA level.
Occasionally, selected morphological features are recorded and the content of
key chemical constituents determined, but unbiased and comprehensive chemical
phenotypes have not been included systematically in diversity surveys. Our
objective in this study was to assess metabolic diversity in sesame by
nontargeted metabolic profiling and elucidate the relationship between
metabolic and genome diversity in this crop. Results: Ten sesame accessions
were selected that represent most of the genome diversity of sesame grown in
India, Western Asia, Sudan and Venezuela based on previous AFLP studies.
Ethanolic seed extracts were separated by HPLC, metabolites were ionized by
positive and negative electrospray and ions were detected with an ion trap mass
spectrometer in full-scan mode for m/z from 50 to 1000. Genome diversity was
determined by Amplified Fragment Length Polymorphism (AFLP) using eight primer
pair combinations. The relationship between biodiversity at the genome and at
the metabolome levels was assessed by correlation analysis and multivariate
statistics. Conclusion: Patterns of diversity at the genomic and metabolic
levels differed, indicating that selection played a significant role in the
evolution of metabolic diversity in sesame. This result implies that when used
for the selection of genotypes in breeding and conservation, diversity
assessment based on neutral DNA markers should be complemented with metabolic
profiles. We hypothesize that this applies to all crops with a long history of
domestication that possess commercially relevant traits affected by chemical
phenotypes.


  Consider the network of all secondary structures of a given RNA sequence,
where nodes are connected when the corresponding structures have base pair
distance one. The expected degree of the network is the average number of
neighbors, where average may be computed with respect to the either the uniform
or Boltzmann probability. Here we describe the first algorithm, RNAexpNumNbors,
that can compute the expected number of neighbors, or expected network degree,
of an input sequence. For RNA sequences from the Rfam database, the expected
degree is significantly less than the CMFE structure, defined to have minimum
free energy over all structures consistent with the Rfam consensus structure.
The expected degree of structural RNAs, such as purine riboswitches,
paradoxically appears to be smaller than that of random RNA, yet the difference
between the degree of the MFE structure and the expected degree is larger than
that of random RNA. Expected degree does not seem to correlate with standard
structural diversity measures of RNA, such as positional entropy, ensemble
defect, etc. The program {\tt RNAexpNumNbors} is written in C, runs in cubic
time and quadratic space, and is publicly available at
http://bioinformatics.bc.edu/clotelab/RNAexpNumNbors.


  Among all kinds of apes, only gibbons have the slim body as human. Gibbons
can move in the forest by cross arm swing, what was the locomotion mode of our
arboreal ancestor. Since our ancestor had much heavier body but weaker arms
than gibbons, we suppose they had to move with two arm brachiation. Such mode
of locomotion can account reasonably for the transition to bipedalism. Firstly,
it needed our ancestor to straighten knee and hip joints and flex their lumbar
spine. secondly, it evolved the feet of our ancestor with longitudinal arche.
And most importantly, it made the ratio of the length of the upper limbs to
that of the lower limbs unsuitable for quadruped walking.


  Plasmids are major players in Horizontal Gene Transfer mechanisms, hence they
are highly variable in their gene content and length. We propose a model for
the fitness of a plasmid as a function of its length, which predicts
diminishing returns. We infer the distribution of plasmid lengths by a simple
evolutionary model and we show that there is a positive correlation between the
presence and efficiency of the conjugation machinery and the length of the
plasmid. The model predicts an "insertion load" on plasmids, which could
explain also the fact that plasmids are widespread in bacterial populations but
rarely established. Finally we discuss how the typical length of plasmids
increases with the amount of stress in the environment, focusing on the recent
human-driven increase in antibiotic concentrations.


  High read depth can be used to assemble short sequence repeats. The existing
genome assemblers fail in repetitive regions of longer than average read.
  I propose a new algorithm for a DNA assembly which uses the relative
frequency of reads to properly reconstruct repetitive sequences. The
mathematical model shows the upper limits of accuracy of the results as a
function of read coverage. For high coverage, the estimation error depends
linearly on repetitive sequence length and inversely proportional to the
sequencing coverage. The algorithm requires high read depth, provided by the
next-generation sequencers and could use the existing data. The tests on
errorless reads, generated in silico from several model genomes, pointed the
properly reconstructed repetitive sequences, where existing assemblers fail.


  The histone-DNA interaction in the nucleosome is a fundamental mechanism of
genomic compaction and regulation, which remains largely unkown despite a
growing structural knowledge of the complex. Here, we propose a framework for
the extraction of a nanoscale histone-DNA force-field from a collection of
high-resolution structures, which may be adapted to a larger class of
protein-DNA complexes. We apply the procedure on a large crystallographic
database extended by snapshots from molecular dynamics simulations. The
comparison of the structural models first shows that, at the sites of
histone-DNA contact, the DNA base-pairs are locally shifted outwards,
consistent with locally repulsive forces exerted by the histones. In a second
step, we show that the various force profiles of the analyzed structures derive
locally from a unique, sequence-independent, quadratic repulsive force field,
while the sequence preferences are entirely due to the internal DNA mechanics.
We thus obtain the first knowledge-derived nanoscale potential for the
histone-DNA interaction in the nucleosome. The conformations obtained by
relaxation of nucleosomal DNA with high-affinity sequences in this potential
accurately reproduce experimental values of binding preferences. We finally
address the more generic binding mechanisms relevant to the 80% genomic
sequences incorporated in nucleosomes, by computing the conformation of
nucleosomal DNA with sequence-averaged properties. This conformation is found
to differ from those found in crystals, and the analysis suggests that
repulsive histone forces are related to a local stretch tension in nucleosomal
DNA, mostly between successive contact points. This tension could play a role
in the stability of the complex.


  The thermodynamic properties of DNA circular molecules are investigated by a
new path integral computational method which treats in the real space the
fundamental forces stabilizing the molecule. The base pair and stacking
contributions to the classical action are evaluated separately by simulating a
broad ensemble of twisted conformations. We obtain, for two short sequences, a
free energy landscape with multiple wells corresponding to the most convenient
values of helical repeat. Our results point to a intrinsic flexibility of the
circular structures in which the base pair fluctuations move the system from
one well to the next thus causing the local unwinding of the helix. The latter
is more pronounced in the shorter sequence whose cyclization causes a higher
bending stress. The entropic reductions associated to the formation of the
ordered helicoidal structure are estimated.


  Book review of Stochastic Chemical Kinetics. Theory and (Mostly) Systems
Biological Applications, P. Erdi, G. Lente. Springer (2014)


  Omics technologies enable unbiased investigation of biological systems
through massively parallel sequence acquisition or molecular measurements,
bringing the life sciences into the era of Big Data. A central challenge posed
by such omics datasets is how to transform this data into biological knowledge.
For example, how to use this data to answer questions such as: which functional
pathways are involved in cell differentiation? Which genes should we target to
stop cancer? Network analysis is a powerful and general approach to solve this
problem consisting of two fundamental stages, network reconstruction and
network interrogation. Herein, we provide an overview of network analysis
including a step by step guide on how to perform and use this approach to
investigate a biological question. In this guide, we also include the software
packages that we and others employ for each of the steps of a network analysis
workflow.


  The mesoscopic reaction-diffusion master equation (RDME) is a popular
modeling framework, frequently applied to stochastic reaction-diffusion
kinetics in systems biology. The RDME is derived from assumptions about the
underlying physical properties of the system, and it may produce unphysical
results for models where those assumptions fail. In that case, other more
comprehensive models are better suited, such as hard-sphere Brownian dynamics
(BD). Although the RDME is a model in its own right, and not inferred from any
specific microscale model, it proves useful to attempt to approximate a
microscale model by a specific choice of mesoscopic reaction rates. In this
paper we derive mesoscopic reaction rates by matching certain statistics of the
RDME solution to statistics of the solution of a widely used microscopic BD
model: the Smoluchowski model with a mixed boundary condition at the reaction
radius of two molecules. We also establish fundamental limits for the range of
mesh resolutions for which this approach yields accurate results, and show both
theoretically and in numerical examples that as we approach the lower
fundamental limit, the mesoscopic dynamics approach the microscopic dynamics.


  Here we focus on the challenge of verifying the correctness of molecular
implementations of abstract chemical reaction networks, where operation in a
well-mixed "soup" of molecules is stochastic, asynchronous, concurrent, and
often involves multiple intermediate steps in the implementation, parallel
pathways, and side reactions. This problem relates to the verification of Petri
nets, but existing approaches are not sufficient for providing a single
guarantee covering an infinite set of possible initial states (molecule counts)
and an infinite state space potentially explored by the system given any
initial state. We address these issues by formulating a new theory of pathway
decomposition that provides an elegant formal basis for comparing chemical
reaction network implementations, and we present an algorithm that computes
this basis. Our theory naturally handles certain situations that commonly arise
in molecular implementations, such as what we call "delayed choice," that are
not easily accommodated by other approaches. We further show how pathway
decomposition can be combined with weak bisimulation to handle a wider class
that includes most currently known enzyme-free DNA implementation techniques.
We anticipate that our notion of logical equivalence between chemical reaction
network implementations will be valuable for other molecular implementations
such as biochemical enzyme systems, and perhaps even more broadly in
concurrency theory.


  High-throughput molecular screening is required to investigate the origin and
diffusion of antimicrobial resistance in pathogen outbreaks. The most frequent
cause of human infection is Escherichia coli, which is dominated by sequence
type 131 (ST131), a set of rapidly radiating pandemic clones. The highly
infectious clades of ST131 originated firstly by a mutation enhancing virulence
and adhesion. Secondly, single-nucleotide polymorphisms occurred enabling
fluoroquinolone-resistance, which is near-fixed in all ST131. Thirdly, broader
resistance through beta-lactamases has been gained and lost frequently,
symptomatic of conflicting environmental selective effects. This flexible
approach to gene exchange is worrying and supports the proposition that ST131
will develop an even wider range of plasmid and chromosomal elements promoting
antimicrobial resistance. To stymie ST131, deep genome sequencing is required
to understand the origin, evolution and spread of antimicrobial resistance
genes. Phylogenetic methods that decipher past events can predict future
patterns of virulence and transmission based on genetic signatures of
adaptation and gene exchange. Both the effect of partial antimicrobial exposure
and cell dormancy caused by variation in gene expression may accelerate the
development of resistance. High-throughput sequencing can decode measurable
evolution of cell populations within patients associated with systems-wide
changes in gene expression during treatments. A multi-faceted approach can
enhance assessment of antimicrobial resistance in E. coli ST131 by examining
transmission dynamics between hosts to achieve a goal of pre-empting resistance
before it emerges by optimising antimicrobial treatment protocols.


  Classifications and phylogenetic inferences of organismal groups change in
light of new insights. Over time these changes can result in an imperfect
tracking of taxonomic perspectives through the re-/use of Code-compliant or
informal names. To mitigate these limitations, we introduce a novel approach
for aligning taxonomies through the interaction of human experts and logic
reasoners. We explore the performance of this approach with the Perelleschus
use case of Franz & Cardona-Duque (2013). The use case includes six taxonomies
published from 1936 to 2013, 54 taxonomic concepts (i.e., circumscriptions of
names individuated according to their respective source publications), and 75
expert-asserted Region Connection Calculus articulations (e.g., congruence,
proper inclusion, overlap, or exclusion). An Open Source reasoning toolkit is
used to analyze 13 paired Perelleschus taxonomy alignments under heterogeneous
constraints and interpretations. The reasoning workflow optimizes the logical
consistency and expressiveness of the input and infers the set of maximally
informative relations among the entailed taxonomic concepts. The latter are
then used to produce merge visualizations that represent all congruent and
non-congruent taxonomic elements among the aligned input trees. In this small
use case with 6-53 input concepts per alignment, the information gained through
the reasoning process is on average one order of magnitude greater than in the
input. The approach offers scalable solutions for tracking provenance among
succeeding taxonomic perspectives that may have differential biases in naming
conventions, phylogenetic resolution, ingroup and outgroup sampling, or
ostensive (member-referencing) versus intensional (property-referencing)
concepts and articulations.


  A number of studies have suggested using comparisons between DNA sequences of
closely related bacterial isolates to estimate the relative rate of
recombination to mutation for that bacterial species. We consider such an
approach which uses single-locus variants: pairs of isolates whose DNA differ
at a single gene locus. One way of deriving point estimates for the relative
rate of recombination to mutation from such data is to use composite likelihood
methods. We extend recent work in this area so as to be able to construct
confidence intervals for our estimates, without needing to resort to
computationally-intensive bootstrap procedures, and to develop a test for
whether the relative rate varies across loci. Both our test and method for
constructing confidence intervals are obtained by modeling the dependence
structure in the data, and then applying asymptotic theory regarding the
distribution of estimators obtained using a composite likelihood. We applied
these methods to multi-locus sequence typing (MLST) data from eight bacteria,
finding strong evidence for considerable rate variation in three of these:
Bacillus cereus, Enterococcus faecium and Klebsiella pneumoniae.


  US public health authorities claim imposing quarantines on healthcare workers
returning from West Africa is incorrect according to science. Their positions
rely upon a set of studies and experience about outbreaks and transmission
mechanisms in Africa as well as assumptions about what those studies imply
about outbreaks in the US. According to this view the probability of a single
infection is low and that of a major outbreak is non-existent. In a series of
brief reports we will provide insight into why properties of networks of
contagion that are not considered in traditional statistics suggest that risks
are higher than those assumptions suggest. We begin with the difference between
thin and fat tailed distributions applied to the number of infected individuals
that can arise from a single one. Traditional epidemiological models consider
the contagion process as described by $R_0$, the average number of new infected
individuals arising from a single case. However, in a complex interdependent
society it is possible for the actual number due to a single individual to
dramatically differ from the average number, with severe consequences for the
ability to contain an outbreak when it is just beginning. Our analysis raises
doubts about the scientific validity of policy recommendations of public health
authorities. We also point out that existing CDC public health policies and
actions are inconsistent with their claims.


  Characterizing interactions between drugs is important to avoid potentially
harmful combinations, to reduce off-target effects of treatments and to fight
antibiotic resistant pathogens, among others. Here we present a network
inference algorithm to predict uncharacterized drug-drug interactions. Our
algorithm takes, as its only input, sets of previously reported interactions,
and does not require any pharmacological or biochemical information about the
drugs, their targets or their mechanisms of action. Because the models we use
are abstract, our approach can deal with adverse interactions,
synergistic/antagonistic/suppressing interactions, or any other type of drug
interaction. We show that our method is able to accurately predict
interactions, both in exhaustive pairwise interaction data between small sets
of drugs, and in large-scale databases. We also demonstrate that our algorithm
can be used efficiently to discover interactions of new drugs as part of the
drug discovery process.


  Ion transport is the major determining factor of salinity tolerance in
plants. A simple scheme of a plant cell with ion fluxes provides basic
understanding of ion transport and the corresponding changes of ion
concentrations under salinity. The review describes in detail basic principles
of ion transport for a plant cell, introduces set of transporters essential for
sodium and potassium uptake and efflux, analyses driving forces of ion
transport and compares ion fluxes measured by several techniques. Study of
differences in ion transport between salt tolerant halophytes and
salt-sensitive plants with an emphasis on transport of potassium and sodium via
plasma membranes offers knowledge for increasing salinity tolerance. Effects of
salt stress on ion transport properties of membranes show huge opportunities
for manipulating ion transport. Several attempts to overexpress or knockout ion
transporters for changing salinity tolerance are described. Future perspectives
are questioned with more attention given to potential candidate ion channels
and transporters for altered expression. The potential direction of increasing
salinity tolerance by modifying ion channels and transporters is discussed and
questioned. An alternative approach from synthetic biology is to modify the
existing membrane transport proteins or create new ones with desired properties
for transforming agricultural crops. The approach had not been widely used
earlier and leads also to theoretical and pure scientific aspects of protein
chemistry, structure-function relations of membrane proteins, systems biology
and physiology of stress and ion homeostasis.


  Pathogenic enterotoxigenic Escherichia coli (ETEC) are the major bacterial
cause of diarrhea in young children in developing countries and in travelers,
causing significant mortality in children. Adhesive fimbriae are a prime
virulence factor for ETEC, initiating colonization of the small intestinal
epithelium. Similar to other Gram-negative bacteria, ETEC express one or more
diverse fimbriae, some assembled by the chaperone-usher pathway and others by
the alternate chaperone pathway. Here we elucidate structural and biophysical
aspects and adaptations of each fimbrial type to its respective host niche.
CS20 fimbriae are compared to CFA/I fimbriae, which are two ETEC fimbriae
assembled via different pathways, and to Pfimbriae from uropathogenic E. coli.
Many fimbriae unwind from their native helical filament to an extended linear
conformation under force, thereby sustaining adhesion by reducing load at the
point of contact between the bacterium and the target cell. CFA/I fimbriae
require the least force to unwind, followed by CS20 fimbriae and then
P-fimbriae, which require the highest unwinding force. We conclude from our
electron microscopy reconstructions, modeling, and force spectroscopy data that
the target niche plays a central role in the biophysical properties of fimbriae
that are critical for bacterial pathophysiology.


  We model biochemical signal transduction, based on a ligand-receptor binding
mechanism, as a discrete-time finite-state Markov channel, which we call the
BIND channel. We show how to obtain the capacity of this channel, for the case
of binary output, binary channel state, and arbitrary finite input alphabets.
We show that the capacity-achieving input distribution is IID. Further, we show
that feedback does not increase the capacity of this channel. We show how the
capacity of the discrete-time channel approaches the capacity of Kabanov's
Poisson channel, in the limit of short time steps and rapid ligand release.


  Toxic fibrillar aggregates of Islet Amyloid PolyPeptide (IAPP) appear as the
physical outcome of a peptidic phase-transition signaling the onset of type-2
diabetes mellitus in different mammalian species. In particular, experimentally
verified mutations on the amyloidogenic segment 20-29 in humans, cats and rats
are highly correlated with the molecular aggregation propensities. Through a
microcanonical analysis of the aggregation of IAPP_{20-29} isoforms, we show
that a minimalist one-bead hydrophobic-polar continuum model for protein
interactions properly quantifies those propensities from free-energy barriers.
Our results highlight the central role of sequence-dependent hydrophobic
mutations on hot spots for stabilization, and so for the engineering, of such
biological peptides.


  The 2014 epidemic of the Ebola virus is governed by a genetically diverse
viral population. In the early Sierra Leone outbreak, a recent study has
identified new mutations that generate genetically distinct sequence clades.
Here we find evidence that major Sierra Leone clades have systematic
differences in growth rate and reproduction number. If this growth
heterogeneity remains stable, it will generate major shifts in clade
frequencies and influence the overall epidemic dynamics on time scales within
the current outbreak. Our method is based on simple summary statistics of clade
growth, which can be inferred from genealogical trees with an underlying
clade-specific birth-death model of the infection dynamics. This method can be
used to perform realtime tracking of an evolving epidemic and identify emerging
clades of epidemiological or evolutionary significance.


  Containing the spreading of crime in urban societies remains a major
challenge. Empirical evidence suggests that, left unchecked, crimes may be
recurrent and proliferate. On the other hand, eradicating a culture of crime
may be difficult, especially under extreme social circumstances that impair the
creation of a shared sense of social responsibility. Although our understanding
of the mechanisms that drive the emergence and diffusion of crime is still
incomplete, recent research highlights applied mathematics and methods of
statistical physics as valuable theoretical resources that may help us better
understand criminal activity. We review different approaches aimed at modeling
and improving our understanding of crime, focusing on the nucleation of crime
hotspots using partial differential equations, self-exciting point process and
agent-based modeling, adversarial evolutionary games, and the network science
behind the formation of gangs and large-scale organized crime. We emphasize
that statistical physics of crime can relevantly inform the design of
successful crime prevention strategies, as well as improve the accuracy of
expectations about how different policing interventions should impact malicious
human activity deviating from social norms. We also outline possible directions
for future research, related to the effects of social and coevolving networks
and to the hierarchical growth of criminal structures due to self-organization.


  Many signs indicate that the graphene could widely occur on the early Earth.
Here, we report a new theory that graphene might be an embryo of protocell
membrane, and found several evidences. Firstly, the graphene oxide and
phospholipid-graphene oxide composite would curl into capsules in strongly
acidic saturated solution of Pb(NO3)2 at low temperature, providing a
protective space for biochemical reactions. Secondly, L-animi acids exhibit
higher reactivity than D-animi acids for graphene oxides in favor of the
formation of left-handed proteins. Thirdly, monolayer graphene with nanopores
prepared by unfocused 84Kr25+ has high selectivity for permeation of the
monovalent metal ions (Rb+ > K+ > Cs+ > Na+ > Li+), but does not allow Cl-
through, which could be attributed to the ion exchange of oxygen-containing
groups on the rim of nanopores. It is similar to K+ channels, which would cause
efflux of some ions from capsule of graphene oxides with the decrease of pH in
the primitive ocean, creating a suitable inner condition for the origin of
life. Consequently, the strongly acidic, high salinity and strong radiation as
well as temperature changes in the early Earth, regarded as negative factors,
would be indispensable for the origin of protocell. In short, graphene bred
life, but digested gradually by the evolution.


  We consider possible designs and experimental realiza-tions in synthesized
rather than naturally occurring bio-chemical systems of a selection of basic
bio-inspired information processing steps. These include feed-forward loops,
which have been identified as the most common information processing motifs in
many natural pathways in cellular functioning, and memory-involving processes,
specifically, associative memory. Such systems should not be designed to
literally mimic nature. Rather, we can be guided by nature's mechanisms for
experimenting with new information/signal processing steps which are based on
coupled biochemical reactions, but are vastly simpler than natural processes,
and which will provide tools for the long-term goal of understanding and
harnessing nature's information processing paradigm. Our biochemical processes
of choice are enzymatic cascades because of their compatibility with
physiological processes in vivo and with electronics (e.g., electrodes) in
vitro allowing for networking and interfacing of enzyme-catalyzed processes
with other chemical and biochemical reactions. In addition to designing and
realizing feed-forward loops and other processes, one has to develop approaches
to probe their response to external control of the time-dependence of the
input(s), by measuring the resulting time-dependence of the output. The goal
will be to demonstrate the expected features, for example, the delayed response
and stabilizing effect of the feed-forward loops.


  DNA nanotechnology promises to provide controllable self-assembly on the
nanoscale, allowing for the design of static structures, dynamic machines and
computational architectures. In this article I review the state-of-the art of
DNA nanotechnology, highlighting the need for a more detailed understanding of
the key processes, both in terms of theoretical modelling and experimental
characterisation. I then consider coarse-grained models of DNA, mesoscale
descriptions that have the potential to provide great insight into the
operation of DNA nanotechnology if they are well designed. In particular, I
discuss a number of nanotechnological systems that have been studied with
oxDNA, a recently developed coarse-grained model, highlighting the subtle
interplay of kinetic, thermodynamic and mechanical factors that can determine
behaviour. Finally, new results highlighting the importance of mechanical
tension in the operation of a two-footed walker are presented, demonstrating
that recovery from an unintended `overstepped' configuration can be accelerated
by three to four orders of magnitude by application of a moderate tension to
the walker's track. More generally, the walker illustrates the possibility of
biasing strand-displacement processes to affect the overall rate.


  Rapid advances in high-throughput technologies have led to considerable
interest in analyzing genome-scale data in the context of biological pathways,
with the goal of identifying functional systems that are involved in a given
phenotype. In the most common approaches, biological pathways are modeled as
simple sets of genes, neglecting the network of interactions comprising the
pathway and treating all genes as equally important to the pathway's function.
Recently, a number of new methods have been proposed to integrate pathway
topology in the analyses, harnessing existing knowledge and enabling more
nuanced models of complex biological systems. However, there is little guidance
available to researches choosing between these methods. In this review, we
discuss eight topology-based methods, comparing their methodological approaches
and appropriate use cases. In addition, we present the results of the
application of these methods to a curated set of ten gene expression profiling
studies using a common set of pathway annotations. We report the computational
efficiency of the methods and the consistency of the results across methods and
studies to help guide users in choosing a method. We also discuss the
challenges and future outlook for improved network analysis methodologies.


  Identifying latent structure in large data matrices is essential for
exploring biological processes. Here, we consider recovering gene co-expression
networks from gene expression data, where each network encodes relationships
between genes that are locally co-regulated by shared biological mechanisms. To
do this, we develop a Bayesian statistical model for biclustering to infer
subsets of co-regulated genes whose covariation may be observed in only a
subset of the samples. Our biclustering method, BicMix, has desirable
properties, including allowing overcomplete representations of the data,
computational tractability, and jointly modeling unknown confounders and
biological signals. Compared with related biclustering methods, BicMix recovers
latent structure with higher precision across diverse simulation scenarios.
Further, we develop a method to recover gene co-expression networks from the
estimated sparse biclustering matrices. We apply BicMix to breast cancer gene
expression data and recover a gene co-expression network that is differential
across ER+ and ER- samples.


  We find an advantage of recombination for a category of complex fitness
landscapes. Recent studies of empirical fitness landscapes reveal complex gene
interactions and multiple peaks, and recombination can be a powerful mechanism
for escaping suboptimal peaks. However classical work on recombination largely
ignores the effect of complex gene interactions. The advantage we find has no
correspondence for 2-locus systems or for smooth landscapes. The effect is
sometimes extreme, in the sense that shutting off recombination could result in
that the organism fails to adapt. A standard question about recombination is if
the mechanism tends to accelerate or decelerate adaptation. However, we argue
that extreme effects may be more important than how the majority falls.


  Understanding the stabilities of proteins in nanopores requires a
quantitative description of confinement induced interactions between amino acid
side chains. We use molecular dynamics simulations to study the nature of
interactions between the side chain pairs ALA-PHE, SER-ASN and LYS-GLU in bulk
water and in water-filled nanopores. The temperature dependence of the bulk
solvent potentials of mean force and the interaction free energies in
cylindrical and spherical nanopores is used to identify the corresponding
entropic and enthalpic components. The entropically stabilized hydrophobic
interaction between ALA and PHE in bulk water is enthalpically dominated upon
confinement depending on the relative orientations between the side chains. In
the case of SER-ASN, hydrogen bonded configurations that are similar in bulk
water are thermodynamically distinct in a cylindrical pore, thus making rotamer
distributions different from those in the bulk. Remarkably, salt bridge
formation between LYS-GLU is stabilized by entropy in contrast to the bulk.
Implications of our findings for confinement-induced alterations in protein
stability are briefly outlined.


  This article proposes an age-structured compartment model for irreversible
diseases with a pre-clinical state of undiagnosed cases that precedes the
diagnosis. The model is able to cope with mortality rates differing between the
pre-clinical and the clinical state (differential mortality). Applicability is
tested in a hypothetical disease with realistic incidence and mortality rates.


  A recently experimentally observed biochemical "threshold filtering"
mechanism by processes catalyzed by the enzyme malate dehydrogenase is
explained in terms of a model that incorporates an unusual mechanism of
inhibition of this enzyme that has a reversible mechanism of action.
Experimental data for a system in which the output signal is produced by
biocatalytic processes of the enzyme glucose dehydrogenase are analyzed to
verify the model's validity. We also establish that fast reversible conversion
of the output product to another compound, without the additional inhibition,
cannot on its own result in filtering.


  The spatial arrangement of trees in a tropical forest reflects the interplay
between aggregating processes, like dispersal limitation, and negative feedback
that induces effective repulsion among individuals. Monitoring the
variance-mean ratio for conspecific individuals along length-scales, we show
that the effect of negative feedback is dominant at short scales, while
aggregation characterizes the large-scale patterns. A comparison of different
species indicates, surprisingly, that both aggregation and negative feedback
scales are related to the overall abundance of the species. This suggests a
bottom-up control mechanism, in which the negative feedback dictates the
dispersal kernel and the overall abundance.


  Various fossil evidences show that hunting is one of major means of ancient
human to get foods. But the running speed of our ancestors was much slower than
quadruped animals, and they did not have sharp claws and canines. So, they have
to rely heavily on stone and wooden tools when they hunting or fighting against
other predators, which are very different from the hunting behaviors of other
carnivores. There are mainly two types of attack and defense action during
human hunting, front or side hit with a wooden stick in hands and stone or
wooden spears throwing, and throwing had play an important role in human
evolution process. But there is almost no work to study the why only human
chose to hunting by this way. Here we suppose that ancient human chose two-arm
brachiation as main arboreal locomotion mode because of their suitable body
weight. Human body traits include slim body, parallel arranged scapulas, long
thumb and powerful grip ability are all evolved as results of two arm
brachiation. The relevant adaptive evolution of the shoulder bone structure
make human arms with a large range of movement and the long thumb makes human
activities to be more accurate and controllable. These are two important body
structure advantages of ancient human which makes them could get from arboreal
life into a whole new hunting and fighting stage.


  Background: Sphingolipids - as suggested by the prefix in their name - are
mysterious molecules, which play surprisingly various roles in opposable
cellular processes, like autophagy, apoptosis, proliferation and
differentiation. Recently they have been also recognized as important
messengers in cellular signalling pathways. More importantly, sphingolipid
metabolism disorders were observed in various pathological conditions such as
cancer and neurodegeneration. Results: Existing formal models of sphingolipids
metabolism concentrates mostly on de novo ceramide synthesis or restrict their
focus to biochemical transformations of a particular subspecies. We propose
first comprehensive computational model of sphingolipid metabolism in human
tissue. In contrast to previous approaches we explicitly model
compartmentalization what allows emphasizing the differences among individual
organelles. Conclusions: Presented here model was validated by means of
recently proposed model analysis technics allowing for detection of most
sensitive and experimentally non-identifiable parameters and determination of
main sources of model variance. Moreover, we demonstrate the utility of the
model for the study of molecular processes underlying Alzheimer's disease.


  Sequencing by synthesis is used in many next-generation DNA sequencing
technologies. Some of the technologies, especially those exploring the
principle of single-molecule sequencing, allow incomplete nucleotide
incorporation in each cycle. We derive statistical distributions for sequencing
by synthesis by taking into account the possibility that nucleotide
incorporation may not be complete in each flow cycle. The statistical
distributions are expressed in terms of nucleotide probabilities of the target
sequences and the nucleotide incorporation probabilities for each nucleotide.
We give exact distributions both for fixed number of flow cycles and for fixed
sequence length. Explicit formulas are derived for the mean and variance of
these distributions. The results are generalizations of our previous work for
pyrosequencing. Incomplete nucleotide incorporation leads to significant change
in the mean and variance of the distributions, but still they can be
approximated by normal distributions with the same mean and variance. The
results are also generalized to handle sequence context dependent
incorporation. The statistical distributions will be useful for instrument and
software development for sequencing by synthesis platforms.


  Sequencing by synthesis is the underlying technology for many next-generation
DNA sequencing platforms. We developed a new model, the fixed flow cycle model,
to derive the distributions of sequence length for a given number of flow
cycles under the general conditions where the nucleotide incorporation is
probabilistic and may be incomplete, as in some single-molecule sequencing
technologies. Unlike the previous model, the new model yields the probability
distribution for the sequence length. Explicit closed form formulas are derived
for the mean and variance of the distribution.


  Pyrosequencing is one of the important next-generation sequencing
technologies. We derive the distribution of the number of positive signals in
pyrograms of this sequencing technology as a function of flow cycle numbers and
nucleotide probabilities of the target sequences. As for the distribution of
sequence length, we also derive the distribution of positive signals for the
fixed flow cycle model. Explicit formulas are derived for the mean and variance
of the distributions. A simple result for the mean of the distribution is that
the mean number of positive signals in a pyrogram is approximately twice the
number of flow cycles, regardless of nucleotide probabilities. The statistical
distributions will be useful for instrument and software development for
pyrosequencing and other related platforms.


  Pyrosequencing is emerging as one of the important next-generation sequencing
technologies. We derive the statistical distributions of this technique in
terms of nucleotide probabilities of the target sequences. We give exact
distributions both for fixed number of flow cycles and for fixed sequence
length. Explicit formulas are derived for the mean and variance of these
distributions. In both cases, the distributions can be approximated accurately
by normal distributions with the same mean and variance. The statistical
distributions will be useful for instrument and software development for
pyrosequencing platforms.


  The de Bruijn graph $G_K$ of a set of strings $S$ is a key data structure in
genome assembly that represents overlaps between all the $K$-length substrings
of $S$. Construction and navigation of the graph is a space and time bottleneck
in practice and the main hurdle for assembling large, eukaryote genomes. This
problem is compounded by the fact that state-of-the-art assemblers do not build
the de Bruijn graph for a single order (value of $K$) but for multiple values
of $K$. More precisely, they build $d$ de Bruijn graphs, each with a specific
order, i.e., $G_{K_1}, G_{K_2}, ..., G_{K_d}$. Although, this paradigm
increases the quality of the assembly produced, it increases the memory by a
factor of $d$ in most cases. In this paper, we show how to augment a succinct
de Bruijn graph representation by Bowe et al. (Proc. WABI, 2012) to support new
operations that let us change order on the fly, effectively representing all de
Bruijn graphs of order up to some maximum $K$ in a single data structure. Our
experiments show our variable-order de Bruijn graph only modestly increases
space usage, construction time, and navigation time compared to a single order
graph.


  We formulate and characterize a model to describe the dynamics of
semiflexible polymers in the presence of activity due to motor proteins
attached irreversibly to a substrate, and a transverse pulling force acting on
one end of the filament. The stochastic binding-unbinding of the motor proteins
and their ability to move along the polymer, generates active forces. As the
pulling force reaches a threshold value, the polymer eventually desorbs from
the substrate. Performing molecular dynamics simulations of the polymer in
presence of a Langevin heat bath, and stochastic motor activity, we obtain
desorption phase diagrams. The correlation time for fluctuations in desorbed
fraction increases as one approaches complete desorption, captured
quantitatively by a power law spectral density. We present theoretical analysis
of the phase diagram using mean field approximations in the weakly bending
limit of the polymer and performing linear stability analysis. This predicts
increase in the desorption force with the polymer bending rigidity, active
velocity and processivity of the motor proteins to capture the main features of
the simulation results.


  Allosteric regulation at distant sites is central to many cellular processes.
In particular, allosteric sites in proteins are a major target to increase the
range and selectivity of new drugs, and there is a need for methods capable of
identifying intra-molecular signalling pathways leading to allosteric effects.
Here, we use an atomistic graph-theoretical approach that exploits Markov
transients to extract such pathways and exemplify our results in an important
allosteric protein, caspase-1. Firstly, we use Markov Stability community
detection to perform a multiscale analysis of the structure of caspase-1 which
reveals that the active conformation has a weaker, less compartmentalised
large-scale structure as compared to the inactive conformation, resulting in
greater intra-protein coherence and signal propagation. We also carry out a
full computational point mutagenesis and identify that only a few residues are
critical to such structural coherence. Secondly, we characterise explicitly the
transients of random walks originating at the active site and predict the
location of a known allosteric site in this protein quantifying the
contribution of individual bonds to the communication pathway between the
active and allosteric sites. Several of the bonds we find have been shown
experimentally to be functionally critical, but we also predict a number of as
yet unidentified bonds which may contribute to the pathway. Our approach offers
a computationally inexpensive method for the identification of allosteric sites
and communication pathways in proteins using a fully atomistic description.


  Biochemical mechanisms with mass action kinetics are often modeled by systems
of polynomial differential equations (DE). Determining directly if the DE
system has multiple equilibria (multistationarity) is difficult for realistic
systems, since they are large, nonlinear and contain many unknown parameters.
Mass action biochemical mechanisms can be represented by a directed bipartite
graph with species and reaction nodes. Graph-theoretic methods can then be used
to assess the potential of a given biochemical mechanism for multistationarity
by identifying structures in the bipartite graph referred to as critical
fragments. In this article we present a graph-theoretic method for conservative
biochemical mechanisms characterized by bounded species concentrations, which
makes the use of degree theory arguments possible. We illustrate the results
with an example of a MAPK network.


  In this paper we first introduce the general stochastic epidemic model for
the spread of infectious diseases. Then we give methods for inferring model
parameters such as the basic reproduction number $R_0$ and vaccination coverage
$v_c$ assuming different types of data from an outbreak such as final outbreak
details and temporal data or observations from an ongoing outbreak. Both
individual heterogeneities and heterogeneous mixing are discussed. We also
provide an overview of statistical methods to perform parameter estimation for
stochastic epidemic models. In the last section we describe the problem of
early outbreak detection in infectious disease surveillance and statistical
models used for this purpose.


  We study the thermodynamics and kinetics of an RNA toehold-mediated strand
displacement reaction with a recently developed coarse-grained model of RNA.
Strand displacement, during which a single strand displaces a different strand
previously bound to a complementary substrate strand, is an essential mechanism
in active nucleic acid nanotechnology and has also been hypothesized to occur
in vivo. We study the rate of displacement reactions as a function of the
length of the toehold and temperature and make two experimentally testable
predictions: that the displacement is faster if the toehold is placed at the 5'
end of the substrate and that the displacement slows down with increasing
temperature for longer toeholds.


  In contrast with ordinary materials, living matter drives its own motion by
generating active, out-of-equilibrium internal stresses. These stresses
typically originate from localized active elements embedded in an elastic
medium, such as molecular motors inside the cell or contractile cells in a
tissue. While many large-scale phenomenological theories of such active media
have been developed, a systematic understanding of the emergence of stress from
the local force-generating elements is lacking. In this paper, we present a
rigorous theoretical framework to study this relationship. We show that the
medium's macroscopic active stress tensor is equal to the active elements'
force dipole tensor per unit volume in both continuum and discrete linear
homogeneous media of arbitrary geometries. This relationship is conserved on
average in the presence of disorder, but can be violated in nonlinear elastic
media. Such effects can lead to either a reinforcement or an attenuation of the
active stresses, giving us a glimpse of the ways in which nature might harness
microscopic forces to create active materials.


  {\alpha}-synuclein ({\alpha}-syn) is an intrinsically disordered protein
which is considered to be one of the causes of Parkinson's disease. This
protein forms amyloid fibrils when in a highly concentrated solution. The
fibril formation of {\alpha}-syn is induced not only by increases in
{\alpha}-syn concentration but also by macromolecular crowding. In order to
investigate the coupled effect of the intrinsic disorder of {\alpha}-syn and
macromolecular crowding, we construct a lattice gas model of {\alpha}-syn in
contact with a crowding agent reservoir based on statistical mechanics. The
main assumption is that {\alpha}-syn can be expressed as coarse-grained
particles with internal states coupled with effective volume; and disordered
states are modeled by larger particles with larger internal entropy than other
states. Thanks to the simplicity of the model, we can exactly calculate the
number of conformations of crowding agents, and this enables us to prove that
the original grand canonical ensemble with a crowding agent reservoir is
mathematically equivalent to a canonical ensemble without crowding agents. In
this expression, the effect of macromolecular crowding is absorbed in the
internal entropy of disordered states; it is clearly shown that the crowding
effect reduces the internal entropy. Based on Monte Carlo simulation, we
provide scenarios of crowding-induced fibril formation. We also discuss the
recent controversy over the existence of helically folded tetramers of
{\alpha}-syn, and suggest that macromolecular crowding is the key to resolving
the controversy.


  Multiple epidemiological models have been proposed to predict the spread of
Ebola in West Africa. These models include consideration of counter-measures
meant to slow and, eventually, stop the spread of the disease. Here, we examine
one component of Ebola dynamics that is of growing concern -- the transmission
of Ebola from the dead to the living. We do so by applying the toolkit of
mathematical epidemiology to analyze the consequences of post-death
transmission. We show that underlying disease parameters cannot be inferred
with confidence from early-stage incidence data (that is, they are not
"identifiable") because different parameter combinations can produce virtually
the same epidemic trajectory. Despite this identifiability problem, we find
robustly that inferences that don't account for post-death transmission tend to
underestimate the basic reproductive number -- thus, given the observed rate of
epidemic growth, larger amounts of post-death transmission imply larger
reproductive numbers. From a control perspective, we explain how improvements
in reducing post-death transmission of Ebola may reduce the overall epidemic
spread and scope substantially. Increased attention to the proportion of
post-death transmission has the potential to aid both in projecting the course
of the epidemic and in evaluating a portfolio of control strategies.


  One approach in phylogenomics to infer the tree of life is based on
concatenated multiple sequence alignments from many genes. Unfortunately, the
resulting so-called supermatrix is usually sparse, that is, not every gene
sequence is available for all species in the supermatrix. Due to the missing
sequence information a phylogenetic inference, assuming that each gene evolves
with its own substitution model, suffers from phylogenetic terraces on which
many phylogenetic trees show the same likelihood. Here, we propose a
phylogenetic terrace aware (PTA) data structure for efficient supermatrix based
tree inference under partition models. PTA avoids likelihood computations for
trees belonging to the same terrace. PTA is implemented in the IQ-TREE
software, and leads to an 1.7 to 6-fold speedup for real data sets compared
with a na\"ive implementation. Speedups are independent on terrace sizes but
correlate with the amount of missing data. Thus, the PTA data structure is well
suited for phylogenomic analyses. IQ-TREE source codes, binaries and
documentation are freely available at http://www.cibiv.at/software/iqtree .


  DNA nanotubes are tubular structures composed of DNA crossover molecules. We
present a bottom up approach for construction and characterization of these
structures. Various possible topologies of nanotubes are constructed such as
6-helix, 8-helix and tri-tubes with different sequences and lengths. We have
used fully atomistic molecular dynamics simulations to study the structure,
stability and elasticity of these structures. Several nanosecond long MD
simulations give the microscopic details about DNA nanotubes. Based on the
structural analysis of simulation data, we show that 6-helix nanotubes are
stable and maintain their tubular structure; while 8-helix nanotubes are
flattened to stabilize themselves. We also comment on the sequence dependence
and effect of overhangs. These structures are approximately four times more
rigid having stretch modulus of ~4000 pN compared to the stretch modulus of
1000 pN of DNA double helix molecule of same length and sequence. The stretch
moduli of these nanotubes are also three times larger than those of PX/JX
crossover DNA molecules which have stretch modulus in the range of 1500-2000
pN. The calculated persistence length is in the range of few microns which is
close to the reported experimental results on certain class of the DNA
nanotubes.


  Modeling long-range epidemic spreading in a random environment, we consider a
quenched disordered, $d$-dimensional contact process with infection rates
decaying with the distance as $1/r^{d+\sigma}$. We study the dynamical behavior
of the model at and below the epidemic threshold by a variant of the
strong-disorder renormalization group method and by Monte Carlo simulations in
one and two spatial dimensions. Starting from a single infected site, the
average survival probability is found to decay as $P(t) \sim t^{-d/z}$ up to
multiplicative logarithmic corrections. Below the epidemic threshold, a
Griffiths phase emerges, where the dynamical exponent $z$ varies continuously
with the control parameter and tends to $z_c=d+\sigma$ as the threshold is
approached. At the threshold, the spatial extension of the infected cluster (in
surviving trials) is found to grow as $R(t) \sim t^{1/z_c}$ with a
multiplicative logarithmic correction, and the average number of infected sites
in surviving trials is found to increase as $N_s(t) \sim (\ln t)^{\chi}$ with
$\chi=2$ in one dimension.


  Motivation: Models of discrete concurrent systems often lead to huge and
complex state transition graphs that represent their dynamics. This makes
difficult to analyse dynamical properties. In particular, for logical models of
biological regulatory networks, it is of real interest to study attractors and
their reachability from specific initial conditions, i.e. to assess the
potential asymptotical behaviours of the system. Beyond the identification of
the reachable attractors, we propose to quantify this reachability.
  Results: Relying on the structure of the state transition graph, we estimate
the probability of each attractor reachable from a given initial condition or
from a portion of the state space. First, we present a quasi-exact solution
with an original algorithm called Firefront, based on the exhaustive
exploration of the reachable state space. Then, we introduce an adapted version
of Monte Carlo simulation algorithm, termed Avatar, better suited to larger
models. Firefront and Avatar methods are validated and compared to other
related approaches, using as test cases logical models of synthetic and
biological networks.
  Availability: Both algorithms are implemented as Perl scripts that can be
freely downloaded from http://compbio.igc.gulbenkian.pt/nmd/node/59 along with
Supplementary Material.


  The spread of tick-borne pathogens represents an important threat to human
and animal health in many parts of Eurasia. Here, we analysed a 9-year time
series of Ixodes ricinus ticks feeding on Apodemus flavicollis mice (main
reservoir-competent host for tick-borne encephalitis, TBE) sampled in Trentino
(Northern Italy). The tail of the distribution of the number of ticks per host
was fitted by three theoretical distributions: Negative Binomial (NB),
Poisson-LogNormal (PoiLN), and Power-Law (PL). The fit with theoretical
distributions indicated that the tail of the tick infestation pattern on mice
is better described by the PL distribution. Moreover, we found that the tail of
the distribution significantly changes with seasonal variations in host
abundance. In order to investigate the effect of different tails of tick
distribution on the invasion of a non-systemically transmitted pathogen, we
simulated the transmission of a TBE-like virus between susceptible and
infective ticks using a stochastic model. Model simulations indicated different
outcomes of disease spreading when considering different distribution laws of
ticks among hosts. Specifically, we found that the epidemic threshold and the
prevalence equilibria obtained in epidemiological simulations with PL
distribution are a good approximation of those observed in simulations feed by
the empirical distribution. Moreover, we also found that the epidemic threshold
for disease invasion was lower when considering the seasonal variation of tick
aggregation.


  In the collective-risk social dilemma, players lose their personal endowments
if contributions to the common pool are too small. This fact alone, however,
does not always deter selfish individuals from defecting. The temptations to
free-ride on the prosocial efforts of others are strong because we are
hardwired to maximize our own fitness regardless of the consequences this might
have for the public good. Here we show that the addition of risky assets to the
personal endowments, both of which are lost if the collective target is not
reached, can contribute to solving the collective-risk social dilemma. In
infinite well-mixed populations risky assets introduce new stable and unstable
mixed steady states, whereby the stable mixed steady state converges to full
cooperation as either the risk of collective failure or the amount of risky
assets increases. Similarly, in finite well-mixed populations the introduction
of risky assets enforces configurations where cooperative behavior thrives. In
structured populations cooperation is promoted as well, but the distribution of
assets amongst the groups is crucial. Surprisingly, we find that the completely
rational allocation of assets only to the most successful groups is not
optimal, and this regardless of whether the risk of collective failure is high
or low. Instead, in low-risk situations bounded rational allocation of assets
works best, while in high-risk situations the simplest uniform distribution of
assets among all the groups is optimal. These results indicate that prosocial
behavior depends sensitively on the potential losses individuals are likely to
endure if they fail to cooperate.


  Both evolution and ecology have long been concerned with the impact of
variable environmental conditions on observed levels of genetic diversity
within and between species. We model the evolution of a quantitative trait
under selection that fluctuates in space and time, and derive an analytical
condition for when these fluctuations promote genetic diversification. As
ecological scenario we use a generalized island model with soft selection
within patches in which we incorporate generation overlap. We allow for
arbitrary fluctuations in the environment including spatio-temporal
correlations and any functional form of selection on the trait. Using the
concepts of invasion fitness and evolutionary branching, we derive a simple and
transparent condition for the adaptive evolution and maintenance of genetic
diversity. This condition relates the strength of selection within patches to
expectations and variances in the environmental conditions across space and
time. Our results unify, clarify, and extend a number of previous results on
the evolution and maintenance of genetic variation under fluctuating selection.
Individual-based simulations show that our results are independent of the
details of the genetic architecture and on whether reproduction is clonal or
sexual. The onset of increased genetic variance is predicted accurately also in
small populations in which alleles can go extinct due to environmental
stochasticity.


  Cellular differentiation and evolution are stochastic processes that can
involve multiple types (or states) of particles moving on a complex,
high-dimensional state-space or "fitness" landscape. Cells of each specific
type can thus be quantified by their population at a corresponding node within
a network of states. Their dynamics across the state-space network involve
genotypic or phenotypic transitions that can occur upon cell division, such as
during symmetric or asymmetric cell differentiation, or upon spontaneous
mutation. Waiting times between transitions can be nonexponentially distributed
and reflect e.g., the cell cycle. Here, we use a multi-type branching processes
to study first passage time statistics for a single cell to appear in a
specific state. We present results for a sequential evolutionary process in
which $L$ successive transitions propel a population from a "wild-type" state
to a given "terminally differentiated," "resistant," or "cancerous" state.
Analytic and numeric results are also found for first passage times across an
evolutionary chain containing a node with increased death or proliferation
rate, representing a desert/bottleneck or an oasis. Processes involving cell
proliferation are shown to be "nonlinear" (even though mean-field equations for
the expected particle numbers are linear) resulting in first passage time
statistics that depend on the position of the bottleneck or oasis. Our results
highlight the sensitivity of stochastic measures to cell division fate and
quantify the limitations of using certain approximations and assumptions (such
as fixed-population and mean-field assumptions) in evaluating fixation times.


  We analyze the performance of a Brownian ratchet in the presence of
geometrical constraints. A two-state model that describes the kinetics of
molecular motors is used to characterize the energetic cost when the motor
proceeds under confinement, in the presence of an external force. We show that
the presence of geometrical constraints has a strong effect on the performance
of the motor. In particular, we show that it is possible to enhance the ratchet
performance by a proper tuning of the parameters characterizing the
environment. These results open the possibility of engineering
entropically-optimized transport devices.


  This paper provides a collection of mathematical and computational tools for
the study of robustness in nonlinear gene regulatory networks, represented by
time- and state-discrete dynamical systems taking on multiple states. The focus
is on networks governed by nested canalizing functions (NCFs), first introduced
in the Boolean context by S. Kauffman. After giving a general definition of
NCFs we analyze the class of such functions. We derive a formula for the
normalized average $c$-sensitivities of multistate NCFs, which enables the
calculation of the Derrida plot, a popular measure of network stability. We
also provide a unique canonical parametrized polynomial form of NCFs. This form
has several consequences. We can easily generate NCFs for varying parameter
choices, and derive a closed form formula for the number of such functions in a
given number of variables, as well as an asymptotic formula. Finally, we
compute the number of equivalence classes of NCFs under permutation of
variables. Together, the results of the paper represent a useful mathematical
framework for the study of NCFs and their dynamic networks.


  A variety of algorithms have been proposed for reconstructing trees that show
the evolutionary relationships between species by comparing differences in
genetic data across present-day taxa. If the leaf-to-leaf distances in a tree
can be accurately estimated, then it is possible to reconstruct this tree from
these estimated distances, using polynomial-time methods such as the popular
`Neighbor-Joining' algorithm. There is a precise combinatorial condition under
which distance-based methods are guaranteed to return a correct tree (in full
or in part) based on the requirement that the input distances all lie within
some `safety radius' of the true distances. Here, we explore a stochastic
analogue of this condition, and mathematically establish upper and lower bounds
on this `stochastic safety radius' for distance-based tree reconstruction
methods. Using simulations, we show how this notion provides a new way to
compare the performance of distance-based tree reconstruction methods. This may
help explain why Neighbor-Joining performs so well, as its stochastic safety
radius appears close to optimal (while its more classical safety radius is the
same as many other less accurate methods).


  High-throughput RNA sequencing (RNA-seq) has emerged as a revolutionary and
powerful technology for expression profiling. Most proposed methods for
detecting differentially expressed (DE) genes from RNA-seq are based on
statistics that compare normalized read counts between conditions. However,
there are few methods considering the expression measurement uncertainty into
DE detection. Moreover, most methods are only capable of detecting DE genes,
and few methods are available for detecting DE isoforms. In this paper, a
Bayesian framework (BDSeq) is proposed to detect DE genes and isoforms with
consideration of expression measurement uncertainty. This expression
measurement uncertainty provides useful information which can help to improve
the performance of DE detection. Three real RAN-seq data sets are used to
evaluate the performance of BDSeq and results show that the inclusion of
expression measurement uncertainty improves accuracy in detection of DE genes
and isoforms. Finally, we develop a GamSeq-BDSeq RNA-seq analysis pipeline to
facilitate users, which is freely available at the website
http://parnec.nuaa.edu.cn/liux/GSBD/GamSeq-BDSeq.html.


  Actin cables are linear cytoskeletal structures that serve as tracks for
myosin-based intracellular transport of vesicles and organelles in both yeast
and mammalian cells. In a yeast cell undergoing budding, cables are in constant
dynamic turnover yet some cables grow from the bud neck toward the back of the
mother cell until their length roughly equals the diameter of the mother cell.
This raises the question: how is the length of these cables controlled? Here we
describe a novel molecular mechanism for cable length control inspired by
recent experimental observations in cells. This antenna mechanism involves
three key proteins: formins, which polymerize actin, Smy1 proteins, which bind
formins and inhibit actin polymerization, and myosin motors, which deliver Smy1
to formins, leading to a length-dependent actin polymerization rate. We compute
the probability distribution of cable lengths as a function of several
experimentally tuneable parameters such as the formin-binding affinity of Smy1
and the concentration of myosin motors delivering Smy1. These results provide
testable predictions of the antenna mechanism of actin-cable length control.


  Transfer of nucleotide strings in the Synechocystis sp. PCC6803 genome is
investigated to exhibit periodic and non-periodic correlation structures by
using the recurrence plot method and the phase space reconstruction technique.
The periodic correlation structures are generated by periodic transfer of
several substrings in long periodic or non-periodic nucleotide strings embedded
in the coding regions of genes. The non-periodic correlation structures are
generated by non-periodic transfer of several substrings covering or
overlapping with the coding regions of genes. In the periodic and non-periodic
transfer, some gaps divide the long nucleotide strings into the substrings and
prevent their global transfer. Most of the gaps are either the replacement of
one base or the insertion/reduction of one base. In the reconstructed phase
space, the points generated from two or three steps for the continuous
iterative transfer via the second maximal distance can be fitted by two lines.
It partly reveals an intrinsic dynamics in the transfer of nucleotide strings.
Due to the comparison of the relative positions and lengths, the substrings
concerned with the non-periodic correlation structures are almost identical to
the mobile elements annotated in the genome. The mobile elements are thus
endowed with the basic results on the correlation structures.


  We present a coarse-grained model of clathrin that is simple enough to be
computationally tractable yet includes key observed qualitative features: a
triskelion structure with excluded volume between legs; assembly of polymorphic
cages in the bulk; formation of buds on a membrane. We investigate the assembly
of our model using both Monte Carlo simulations and molecular dynamics with
hydrodynamic interactions, in the latter employing a new membrane boundary
condition. In the bulk, a range of known clathrin structures are assembled. A
membrane budding pathway involving the coalescence of multiple small clusters
is identified.


  Cyclic, nonhierarchical interactions among biological species represent a
general mechanism by which ecosystems are able to maintain high levels of
biodiversity. However, species coexistence is often possible only in spatially
extended systems with a limited range of dispersal, whereas in well-mixed
environments models for cyclic competition often lead to a loss of
biodiversity. Here we consider the dispersal of biological species in a fluid
environment, where mixing is achieved by a combination of advection and
diffusion. In particular, we perform a detailed numerical analysis of a model
composed of turbulent advection, diffusive transport, and cyclic interactions
among biological species in two spatial dimensions and discuss the
circumstances under which biodiversity is maintained when external
environmental conditions, such as resource supply, are uniform in space. Cyclic
interactions are represented by a model with three competitors, resembling the
children's game of rock-paper-scissors, whereas the flow field is obtained from
a direct numerical simulation of two-dimensional turbulence with
hyperviscosity. It is shown that the space-averaged dynamics undergoes
bifurcations as the relative strengths of advection and diffusion compared to
biological interactions are varied.


  We study a classical multitype Galton--Watson process with mutation and
selection. The individuals are sequences of fixed length over a finite
alphabet. On the sharp peak fitness landscape together with independent
mutations per locus, we show that, as the length of the sequences goes to
$\infty$ and the mutation probability goes to 0, the asymptotic relative
frequency of the sequences differing on $k$ digits from the master sequence
approaches $$(\sigma
e^{-a}-1)\frac{a^k}{k!}\sum_{i\geq1}\frac{i^k}{\sigma^i}\,,$$ where $\sigma$ is
the selective advantage of the master sequence and $a$ is the product of the
length of the chains with the mutation probability. The probability
distribution $\mathcal{Q}(\sigma,a)$ on the non negative integers given by the
above formula is the quasispecies distribution with parameters $\sigma$ and
$a$.