  There is a critical need for standard approaches to assess, report, and
compare the technical performance of genome-scale differential gene expression
experiments. We assess technical performance with a proposed "standard"
dashboard of metrics derived from analysis of external spike-in RNA control
ratio mixtures. These control ratio mixtures with defined abundance ratios
enable assessment of diagnostic performance of differentially expressed
transcript lists, limit of detection of ratio (LODR) estimates, and expression
ratio variability and measurement bias. The performance metrics suite is
applicable to analysis of a typical experiment, and here we also apply these
metrics to evaluate technical performance among laboratories. An
interlaboratory study using identical samples shared amongst 12 laboratories
with three different measurement processes demonstrated generally consistent
diagnostic power across 11 laboratories. Ratio measurement variability and bias
were also comparable amongst laboratories for the same measurement process.
Different biases were observed for measurement processes using different mRNA
enrichment protocols.


  Firestone (2014) asserted evidence for numerous (23) nearby (d<300 pc)
supernovae within the Middle and Late Pleistocene. If true, this would have
strong implications for the irradiation of the Earth; at this rate, mass
extinction level events due to supernovae would be more frequent than 100 Myr.
However, there are numerous errors in the application of past research. The
paper overestimates likely nitrate and 14C production from moderately nearby
supernovae by about four orders of magnitude. Moreover, the results are based
on wrongly selected (obsolete) nitrate and 14C datasets. The use of correct and
up-to-date datasets does not confirm the claimed results. The claims in the
paper are invalidated.


  We study the evolution of cooperation in the collective-risk social dilemma
game, where the risk is determined by a collective target that must be reached
with individual contributions. All players initially receive endowments from
the available amount of common resources. While cooperators contribute part of
their endowment to the collective target, defectors do not. If the target is
not reached, the endowments of all players are lost. In our model, we introduce
a feedback between the amount of common resources and the contributions of
cooperators. We show that cooperation can be sustained only if the common
resources are preserved but never excessively abound. This, however, requires a
delicate balance between the amount of common resources that initially exist,
and the amount cooperators contribute to the collective target. Exceeding
critical thresholds in either of the two amounts leads to loss of cooperation,
and consequently to the depletion of common resources.


  Thuiller et al. analyzed the consequences of anticipated climate change on
plant, bird, and mammal phylogenetic diversity (PD) across Europe. They
concluded that species loss will not be clade specific across the Tree of Life,
and that there will not be an overall decline in PD across the whole of Europe.
We applaud their attempt to integrate phylogenetic knowledge into scenarios of
future extinction but their analyses raise a series of concerns. We focus here
on their analyses of plants.


  As the beneficial effects of curcumin have often been reported to be limited
to its small concentrations, we have undertaken a study to find the aggregation
properties of curcumin in water by varying the number of monomers. Our
molecular dynamics simulation results show that the equilibrated structure is
always an aggregated state with remarkable structural rearrangements as we vary
the number of curcumin monomers from 4 to 16 monomers. We find that curcumin
monomers form clusters in a very definite pattern where they tend to aggregate
both in parallel and anti-parallel orientation of the phenyl rings, often seen
in the formation of beta-sheet in proteins. A considerable enhancement in the
population of parallel alignments is observed with increasing the system size
from 12 to 16 curcumin monomers. Due to the prevalence of such parallel
alignment for large system size, a more closely packed cluster is formed with
maximum number of hydrophobic contacts. We also follow the pathway of cluster
growth, in particular the transition from the initial segregated to the final
aggregated state. We find the existence of a metastable structural intermediate
involving a number of intermediate-sized clusters dispersed in the solution.
The course of aggregation bears similarity to nucleation and growth in highly
metastable state. The final aggregated form remains stable with total exclusion
of water from its sequestered hydrophobic core. We also investigate water
structure near the cluster surface along with their orientation. We find that
water molecules form a distorted tetrahedral geometry in the 1st solvation
layer of the cluster, interacting strongly with hydrophilic groups at the
surface of curcumin. The dynamics of such quasi-bound water molecules near the
surface of curcumin cluster is considerably slower than the bulk signifying a
restricted motion as often found in protein hydration layer.


  System-environment interactions are intrinsically nonlinear and dependent on
the interplay between many degrees of freedom. The complexity may be even more
pronounced when one aims to describe biologically motivated systems. In that
case, it is useful to resort to simplified models relying on effective
stochastic equations. A natural consideration is to assume that there is a
noisy contribution from the environment, such that the parameters which
characterize it are not constant but instead fluctuate around their
characteristic values. From this perspective, we propose a stochastic
generalization of the nonlocal Fisher-KPP equation where, as a first step,
environmental fluctuations are Gaussian white noises, both in space and time.
We apply analytical and numerical techniques to study how noise affects
stability and pattern formation in this context. Particularly, we investigate
noise induced coherence by means of the complementary information provided by
the dispersion relation and the structure function.


  We show that a chemostat community of bacteria and bacteriophage in which
bacteria compete for a single nutrient and for which the bipartite infection
network is perfectly nested is permanent, a.k.a. uniformly persistent, provided
that bacteria that are superior competitors for nutrient devote the least to
defence against infection and the virus that are the most efficient at
infecting host have the smallest host range. This confirms earlier work of
Jover et al \cite{Jover} who raised the issue of whether nested infection
networks are permanent. In addition, we provide sufficient conditions that a
bacteria-phage community of arbitrary size with nested infection network can
arise through a succession of permanent subcommunties each with a nested
infection network by the successive addition of one new population.


  We present a minimal model to describe the onset of collective motion seen
when a population of locusts are placed in an annular arena. At low densities
motion is disordered, while at high densities locusts march in a common
direction, which may reverse during the experiment. The data is well-captured
by an individual-based model, in which demographic noise leads to the observed
density-dependent effects. By fitting the model parameters to equation-free
coefficients, we give a quantitative comparison, showing time series,
stationary distributions and the mean switching times between states.


  The environment in which a population evolves can have a crucial impact on
selection. We study evolutionary dynamics in finite populations of fixed size
in a changing environment. The population dynamics are driven by birth and
death events. The rates of these events may vary in time depending on the state
of the environment, which follows an independent Markov process. We develop a
general theory for the fixation probability of a mutant in a population of
wild-types, and for mean unconditional and conditional fixation times. We apply
our theory to evolutionary games for which the payoff structure varies in time.
The mutant can exploit the environmental noise; a dynamic environment that
switches between two states can lead to a probability of fixation that is
higher than in any of the individual environmental states. We provide an
intuitive interpretation of this surprising effect. We also investigate
stationary distributions when mutations are present in the dynamics. In this
regime, we find two approximations of the stationary measure. One works well
for rapid switching, the other for slowly fluctuating environments.


  We show that a mesoscale model, with a minimal number of parameters, can well
describe the thermomechanical and mechanochemical behavior of homogeneous DNA
at thermal equilibrium under tension and torque. We predict critical
temperatures for denaturation under torque and stretch, phase diagrams for
stable DNA, probe/response profiles under mechanical loads, and the density of
dsDNA as a function of stretch and twist. We compare our predictions with
available single molecule manipulation experiments and find strong agreement.
In particular we elucidate the difference between angularly constrained and
unconstrained overstretching. We propose that the smoothness of the angularly
constrained overstreching transition is a consequence of the molecule being in
the vicinity of criticality for a broad range of values of applied tension.


  The question of biological stability (permanence) of a replicator
reaction-diffusion system is considered. Sufficient conditions of biological
stability are found. It is proved that there are situations when biologically
unstable non-distributed replicator system becomes biologically stable in the
distributed case. Numerical examples illustrate analytical findings. This
manuscript is a continuation of arXiv:1308.5631.


  Building better models of cellular pathways is one of the major challenges of
systems biology and functional genomics. There is a need for methods to build
on established expert knowledge and reconcile it with results of
high-throughput studies. Moreover, the available data sources are heterogeneous
and need to be combined in a way specific for the part of the pathway in which
they are most informative. Here, we present a compartment specific strategy to
integrate edge, node and path data for the refinement of a network hypothesis.
Specifically, we use a local-move Gibbs sampler for refining pathway hypotheses
from a compendium of heterogeneous data sources, including novel methodology
for integrating protein attributes. We demonstrate the utility of this approach
in a case study of the pheromone response MAPK pathway in the yeast S.
cerevisiae.


  Here we will use results of Cox, Durrett, and Perkins for voter model
perturbations to study spatial evolutionary games on $Z^d$, $d\ge 3$ when the
interaction kernel is finite range, symmetric, and has covariance matrix
$\sigma^2I$. The games we consider have payoff matrices of the form ${\bf 1} +
wG$ where ${\bf 1}$ is matrix of all 1's and $w$ is small and positive. Since
our population size $N=\infty$, we call our selection small rather than weak
which usually means $w =O(1/N)$.
  The key to studying these games is the fact that when the dynamics are
suitably rescaled in space and time they convergence to solutions of a reaction
diffusion equation (RDE). Inspired by work of Ohtsuki and Nowak and Tarnita et
al we show that the reaction term is the replicator equation for a modified
game matrix and the modifications of the game matrix depend on the interaction
kernel only through the values of two or three simple probabilities for an
associated coalescing random walk.
  Two strategy games lead to an RDE with a cubic nonlinearity, so we can
describe the phase diagram completely. Three strategy games lead to a pair of
coupled RDE, but using an idea from our earlier work, we are able to show that
if there is a repelling function for the replicator equation for the modified
game, then there is coexistence in the spatial game when selection is small.
This enables us to prove coexistence in the spatial model in a wide variety of
examples where the replicator equation of the modified game has an attracting
equilibrium with all components positive. Using this result we are able to
analyze the behavior of four evolutionary games that have recently been used in
cancer modeling.


  In two previous papers, I introduced SuperSpreader (SS) epidemic models,
offered some theoretical discussion of prevention issues, and fitted some
models to data derived from published accounts of the ongoing MERS epidemic
(concluding that a pandemic is likely). Continuing on this theme, here I
discuss prediction: whether, in a disease outbreak driven by superspreader
events, a rigorous decision point---meaning a declaration that a pandemic is
imminent---can be defined. I show that all sources of prediction bias
contribute to generating false negatives (i.e., discounting the chance of a
pandemic when it is looming or has already started). Nevertheless, the
statistical difficulties can be overcome by improved data gathering and use of
known techniques that decrease bias. One peculiarity of the SS epidemic is that
the prediction can sometimes be made long before the actual pandemic onset,
generating lead time to alert the medical community and the public. Thus
modeling is useful to overcome a false sense of security arising from the long
"kindling times" characteristic of SS epidemics and certain
political/psychological factors, as well as improve the public health response.


  The large scale behaviour of a population of cells that grow and interact
through the concentration field of the chemicals they secrete is studied using
dynamical renormalization group methods. The combination of the effective
long-range chemotactic interaction and lack of number conservation leads to a
rich variety of phase behaviour in the system, which includes a sharp
transition from a phase that has moderate (or controlled) growth and regulated
chemical interactions to a phase with strong (or uncontrolled) growth and no
chemical interactions. The transition point has nontrivial critical exponents.
Our results might help shed light on the interplay between chemical signalling
and growth in tissues and colonies, and in particular on the challenging
problem of cancer metastasis.


  Approximate Bayesian computation (ABC) methods provide an elaborate approach
to Bayesian inference on complex models, including model choice. Both
theoretical arguments and simulation experiments indicate, however, that model
posterior probabilities may be poorly evaluated by standard ABC techniques. We
propose a novel approach based on a machine learning tool named random forests
to conduct selection among the highly complex models covered by ABC algorithms.
We thus modify the way Bayesian model selection is both understood and
operated, in that we rephrase the inferential goal as a classification problem,
first predicting the model that best fits the data with random forests and
postponing the approximation of the posterior probability of the predicted MAP
for a second stage also relying on random forests. Compared with earlier
implementations of ABC model choice, the ABC random forest approach offers
several potential improvements: (i) it often has a larger discriminative power
among the competing models, (ii) it is more robust against the number and
choice of statistics summarizing the data, (iii) the computing effort is
drastically reduced (with a gain in computation efficiency of at least fifty),
and (iv) it includes an approximation of the posterior probability of the
selected model. The call to random forests will undoubtedly extend the range of
size of datasets and complexity of models that ABC can handle. We illustrate
the power of this novel methodology by analyzing controlled experiments as well
as genuine population genetics datasets. The proposed methodologies are
implemented in the R package abcrf available on the CRAN.


  The Letter from Rosenfeld (2014, Biological Conservation) in response to
Jamieson and Allendorf (2012, Trends in Ecology and Evolution) and Frankham et
al. (2014, Biological Conservation) and related papers is misleading in places
and requires clarification and correction. We provide those here.


  Synthetic lethal reaction/gene-sets are sets of reactions/genes where only
the simultaneous removal of all reactions/genes in the set abolishes growth of
an organism. In silico, synthetic lethal sets can be identified by simulating
the effect of removal of gene sets from the reconstructed genome-scale
metabolic network of an organism. Flux balance analysis (FBA), based on linear
programming, has emerged as a powerful tool for the in silico analyses of
metabolic networks. To identify all possible synthetic lethal reactions
combinations, an exhaustive sampling of all possible combinations is
computationally expensive. We surmount the computational complexity of
exhaustive search by iteratively restricting the sample space of reaction
combinations for search, resulting in a substantial reduction in the running
time. We here propose an algorithm, Fast-SL, which provides an efficient way to
analyse metabolic networks for higher order lethal reaction sets. Fast-SL
offers a substantial speed-up through a massive reduction in the search space
for synthetic lethals; in the case of E. coli, Fast-SL reduces the search space
for synthetic lethal triplets by over 4000-fold. Fast-SL also compares
favourably with SL Finder, an algorithm for identifying synthetic lethal sets,
by Suthers et al (2009), which involves the solution of a bi-level Mixed
Integer Linear Programming problem. We have implemented the Fast-SL algorithm
in MATLAB, building upon COBRA toolbox v2.0.


  Dynamical systems arising from chemical reaction networks with mass action
kinetics are the subject of chemical reaction network theory (CRNT). In
particular, this theory provides statements about uniqueness, existence, and
stability of positive steady states for all rate constants and initial
conditions. In terms of the corresponding polynomial equations, the results
guarantee uniqueness and existence of positive solutions for all positive
parameters.
  We address a recent extension of CRNT, called generalized mass-action
systems, where reaction rates are allowed to be power-laws in the
concentrations. In particular, the (real) kinetic orders can differ from the
(integer) stoichiometric coefficients. As with mass-action kinetics, complex
balancing equilibria are determined by the graph Laplacian of the underlying
network and can be characterized by binomial equations and parametrized by
monomials. In algebraic terms, we focus on a constructive characterization of
positive solutions of polynomial equations with real and symbolic exponents.
  Uniqueness and existence for all rate constants and initial conditions
additionally depend on sign vectors of the stoichiometric and kinetic-order
subspaces. This leads to a generalization of Birch's theorem, which is robust
with respect to certain perturbations in the exponents. In this context, we
discuss the occurrence of multiple complex balancing equilibria. We illustrate
our results by a running example and provide a MAPLE worksheet with
implementations of all algorithmic methods.


  Cluster-based descriptions of biological networks have received much
attention in recent years fostered by accumulated evidence of the existence of
meaningful correlations between topological network clusters and biological
functional modules. Several well-performing clustering algorithms exist to
infer topological network partitions. However, due to respective technical
idiosyncrasies they might produce dissimilar modular decompositions of a given
network. In this contribution, we aimed to analyze how alternative modular
descriptions could condition the outcome of follow-up network biology analysis.
We considered a human protein interaction network and two paradigmatic cluster
recognition algorithms, namely: the Clauset-Newman-Moore and the infomap
procedures. We analyzed at what extent both methodologies yielded different
results in terms of granularity and biological congruency. In addition, taking
into account Guimera cartographic role characterization of network nodes, we
explored how the adoption of a given clustering methodology impinged on the
ability to highlight relevant network meso-scale connectivity patterns. As a
case study we considered a set of aging related proteins, and showed that only
the high-resolution modular description provided by infomap, could unveil
statistically significant associations between them and inter-intra modular
cartographic features. Besides reporting novel biological insights that could
be gained from the discovered associations, our contribution warns against
possible technical concerns that might affect the tools used to mine for
interaction patterns in network biology studies. In particular our results
suggested that sub-optimal partitions from the strict point of view of their
modularity levels might still be worth being analyzed when meso-scale features
were to be explored in connection with external source of biological knowledge.


  We suggest a possible correlation between the ionization events caused by the
background neutron radiation and the experimental data on mutations with damage
in the DNA repair mechanism, coming from the Long Term Evolution Experiment in
E. Coli populations.


  Collective behavior in cellular populations is coordinated by biochemical
signaling networks within individual cells. Connecting the dynamics of these
intracellular networks to the population phenomena they control poses a
considerable challenge because of network complexity and our limited knowledge
of kinetic parameters. However, from physical systems we know that behavioral
changes in the individual constituents of a collectively-behaving system occur
in a limited number of well-defined classes, and these can be described using
simple models. Here we apply such an approach to the emergence of collective
oscillations in cellular populations of the social amoeba Dictyostelium
discoideum. Through direct tests of our model with quantitative in vivo
measurements of single-cell and population signaling dynamics, we show how a
simple model can effectively describe a complex molecular signaling network and
its effects at multiple size and temporal scales. The model predicts novel
noise-driven single-cell and population-level signaling phenomena that we then
experimentally observe. Our results suggest that like physical systems,
collective behavior in biology may be universal and described using simple
mathematical models.


  We consider the adaptation dynamics of an asexual population that walks
uphill on a rugged fitness landscape which is endowed with large number of
local fitness peaks. We work in a parameter regime where only those mutants
that are single mutation away are accessible, as a result of which the
population eventually gets trapped at a local fitness maximum and the adaptive
walk terminates. We study how the number of adaptive steps taken by the
population before reaching a local fitness peak depends on the initial fitness
of the population, the extreme value distribution of the beneficial mutations
and correlations amongst the fitnesses. Assuming that the relative fitness
difference between successive steps is small, we analytically calculate the
average walk length for both uncorrelated and correlated fitnesses in all
extreme value domains for a given initial fitness. We present numerical results
for the model where the fitness differences can be large, and find that the
walk length behavior differs from that in the former model in the Fr\'echet
domain of extreme value theory. We also discuss the relevance of our results to
microbial experiments.


  The United Nations Joint Programme on HIV/AIDS (UNAIDS) has set a target to
ensure that 15 million HIV-positive people in the world would be receiving
combination anti-retroviral treatment (ART) by 2015. This target is likely to
be reached and new targets for 2020 and 2030 are needed. Eastern and Southern
Africa (ESA) account for approximately half of all people living with HIV in
the world and it will be especially important to set reachable and affordable
targets for this region. In order to make future projections of HIV and TB
prevalence, incidence and mortality assuming different levels of ART scale-up
and coverage, it is first necessary to assess the current state of the
epidemic. Here we review national data on the prevalence of HIV, the coverage
of ART and the notification rates of TB to provide a firm basis for making
future projections. We use the data to assess the extent to which behaviour
change and ART have reduced the number of people living with HIV who remain
infectious.


  Populations can evolve in order to adapt to external changes. The capacity to
evolve and adapt makes successful treatment of infectious diseases and cancer
difficult. Indeed, therapy resistance has quickly become a key challenge for
global health. Therefore, ideas of how to control evolving populations in order
to overcome this threat are valuable. Here we use the mathematical concepts of
stochastic optimal control to study what is needed to control evolving
populations. Following established routes to calculate control strategies, we
first study how a polymorphism can be maintained in a finite population by
adaptively tuning selection. We then introduce a minimal model of drug
resistance in a stochastically evolving cancer cell population and compute
adaptive therapies, where decisions are based on monitoring the response of the
tumor, which can outperform established therapy paradigms. For both case
studies, we demonstrate the importance of high-resolution monitoring of the
target population in order to achieve a given control objective: to control one
must monitor.


  We study the evolution of cooperation in the spatial prisoner's dilemma game,
where besides unconditional cooperation and defection, tit-for-tat,
win-stay-lose-shift and extortion are the five competing strategies. While
pairwise imitation fails to sustain unconditional cooperation and extortion
regardless of game parametrization, myopic updating gives rise to the
coexistence of all five strategies if the temptation to defect is sufficiently
large or if the degree distribution of the interaction network is
heterogeneous. This counterintuitive evolutionary outcome emerges as a result
of an unexpected chain of strategy invasions. Firstly, defectors emerge and
coarsen spontaneously among players adopting win-stay-lose-shift. Secondly,
extortioners and players adopting tit-for-tat emerge and spread via neutral
drift among the emerged defectors. And lastly, among the extortioners,
cooperators become viable too. These recurrent evolutionary invasions yield a
five-strategy phase that is stable irrespective of the system size and the
structure of the interaction network, and they reveal the most unexpected
mechanism that stabilizes extortion and cooperation in an evolutionary setting.


  Breast cancer has been reported to account for the maximum cases among all
female cancers till date. In order to gain a deeper insight into the
complexities of the disease, we analyze the breast cancer network and its
normal counterpart at the proteomic level. While the short range correlations
in the eigenvalues exhibiting universality provide an evidence towards the
importance of random connections in the underlying networks, the long range
correlations along with the localization properties reveal insightful
structural patterns involving functionally important proteins. The analysis
provides a benchmark for designing drugs which can target a subgraph instead of
individual proteins.


  At the present time reliably established that probability density functions
of gene expression of microarray experiments possess a number of universal
properties. First of all these distributions have power asymptotic and secondly
the shape of these distributions are inherent for all organisms and tissues.
This fact led to appearance of a number works where authors are investigating
various probability distributions for approximation of empirical distributions
of gene expression. Nevertheless all these distributions aren't limit
distribution and aren't solution of any equations. These facts by our opinion
are essential shortcoming of these probability laws. Besides, expression of
individual gene aren't accidental event and it depends from expression other
genes. This allows to talk about existence of genic regulatory net in the cell.
In the work the class of fractional stable distributions (FSD) are described.
This class of distributions is limit distribution of sum independent identical
distributed random variables. These distributions have power-law asymptotic and
this fact allow us to apply their for approximation of experimental densities
gene expression of microarray experiments. The parameters of FSDs are
statistically estimated by experimental dates and empirical density is compared
whit theoretical density. In the work the algorithms of parameters estimation
and simulating of FSD variables are presented. The results of such comparison
allow to make conclusion that empirical densities of gene expression can be
approximate by FSD.


  Noncoding RNAs are integral to a wide range of biological processes,
including translation, gene regulation, host-pathogen interactions and
environmental sensing. While genomics is now a mature field, our capacity to
identify noncoding RNA elements in bacterial and archaeal genomes is hampered
by the difficulty of de novo identification. The emergence of new technologies
for characterizing transcriptome outputs, notably RNA-seq, are improving
noncoding RNA identification and expression quantification. However, a major
challenge is to robustly distinguish functional outputs from transcriptional
noise. To establish whether annotation of existing transcriptome data has
effectively captured all functional outputs, we analysed over 400 publicly
available RNA-seq datasets spanning 37 different Archaea and Bacteria. Using
comparative tools, we identify close to a thousand highly-expressed candidate
noncoding RNAs. However, our analyses reveal that capacity to identify
noncoding RNA outputs is strongly dependent on phylogenetic sampling.
Surprisingly, and in stark contrast to protein-coding genes, the phylogenetic
window for effective use of comparative methods is perversely narrow:
aggregating public datasets only produced one phylogenetic cluster where these
tools could be used to robustly separate unannotated noncoding RNAs from a null
hypothesis of transcriptional noise. Our results show that for the full
potential of transcriptomics data to be realized, a change in experimental
design is paramount: effective transcriptomics requires phylogeny-aware
sampling.


  Here we present the results of a large-scale bioinformatic annotation of
non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models
of hand-curated families from the Rfam database to infer conserved RNA families
within each avian genome. We supplement these annotations with predictions from
the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We show that
a number of lncRNA-associated loci are conserved between birds and mammals,
including several intriguing cases where the reported mammalian lncRNA function
is not conserved in birds. We also demonstrate extensive conservation of
classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g.,
snoRNAs and miRNAs) in birds. Furthermore, we describe numerous "losses" of
several RNA families, and attribute these to genuine loss, divergence or
missing data. In particular, we show that many of these losses are due to the
challenges associated with assembling Avian microchromosomes. These combined
results illustrate the utility of applying homology-based methods for
annotating novel vertebrate genomes.


  Tumors often contain multiple subpopulations of cancerous cells defined by
distinct somatic mutations. We describe a new method, PhyloWGS, that can be
applied to WGS data from one or more tumor samples to reconstruct complete
genotypes of these subpopulations based on variant allele frequencies (VAFs) of
point mutations and population frequencies of structural variations. We
introduce a principled phylogenic correction for VAFs in loci affected by copy
number alterations and we show that this correction greatly improves subclonal
reconstruction compared to existing methods.


  We prove identifiability of the tree parameters of the 3-class Jukes-Cantor
mixture model. The proof uses ideas from algebraic statistics, in particular:
finding phylogenetic invariants that separate the varieties associated to
different triples of trees; computing dimensions of the resulting phylogenetic
varieties; and using the disentangling number to reduce to trees with a small
number of leaves. Symbolic computation also plays a key role in handling the
many different cases and finding relevant phylogenetic invariants.


  We define a measure of coherent activity for gene regulatory networks, a
property that reflects the unity of purpose between the regulatory agents with
a common target. We propose that such harmonious regulatory action is desirable
under a demand for energy efficiency and may be selected for under evolutionary
pressures. We consider two recent models of the cell-cycle regulatory network
of the budding yeast, Saccharomyces cerevisiae, as a case study and calculate
their degree of coherence. A comparison with random networks of similar size
and composition reveals that the yeast's cell-cycle regulation is wired to
yield and exceptionally high level of coherent regulatory activity. We also
investigate the mean degree of coherence as a function of the network size,
connectivity and the fraction of repressory/activatory interactions.


  It's been said that "mathematics is biology's next microscope, only better;
biology is mathematics' next physics, only better". Here we aim for something
even better. We try to combine mathematical physics and biology into a
picoscope of life. For this we merge techniques which have been introduced and
developed in modern mathematical physics, largely by Ludvig Faddeev to describe
objects such as solitons and Higgs and to explain phenomena such as anomalies
in gauge fields. We propose a synthesis that can help to resolve the protein
folding problem, one of the most important conundrums in all of science. We
apply the concept of gauge invariance to scrutinize the extrinsic geometry of
strings in three dimensional space. We evoke general principles of symmetry in
combination with Wilsonian universality and derive an essentially unique
Landau-Ginzburg energy that describes the dynamics of a generic string-like
configuration in the far infrared. We observe that the energy supports
topological solitons, that pertain to an anomaly in the manner how a string is
framed around its inflection points. We explain how the solitons operate as
modular building blocks from which folded proteins are composed. We describe
crystallographic protein structures by multi-solitons with experimental
precision, and investigate the non-equilibrium dynamics of proteins under
varying temperature. We simulate the folding process of a protein at in vivo
speed and with close to pico-scale accuracy using a standard laptop computer:
With pico-biology as mathematical physics' next pursuit, things can only get
better.


  Modeling and optimization of metabolic networks has been one of the hottest
topics in computational systems biology within recent years. However, the
complexity and uncertainty of these networks in addition to the lack of
necessary data has resulted in more efforts to design and usage of more capable
models which fit to realistic conditions. In this paper, instead of optimizing
networks in equilibrium condition, the optimization of dynamic networks in
non-equilibrium states including low number of molecules has been studied using
a 2-D lattice simulation. A prototyped network has been simulated with such
approach, and has been optimized using Swarm Particle Algorithm the results of
which are presented in addition to the relevant plots.


  Concepts rooted in physics are becoming increasingly important in biology as
we transition to an era in which quantitative descriptions of all processes
from molecular to cellular level are needed. In this essay I discuss two
unexpected findings of universal behavior, uncommon in biology, in the
self-assembly of proteins and RNA. These findings, which are surprising, reveal
that physics ideas applied to biological problems ranging from folding to gene
expression to cellular movement and communication between cells might lead to
discovery of universal principles operating in adoptable living systems.


  Analysis of high-density superresolution imaging of receptors reveal the
organization of dendrites at the nano-scale resolution. We present here
simulations in empirical live cell images, which allows converting local
information extracted from short range trajectories into simulations of long
range trajectories. Based on these empirical simulations, we compute the
residence time of an AMPA receptor (AMPAR) in dendritic spines that accounts
for receptors local interactions and geometrical organization. We report here
that depending on the type of the spine, the residence time varies from one to
five minutes. Moreover, we show that there exists transient organized
structures, previously described as potential wells that can regulate the
trafficking of AMPARs to dendritic spines.


  A simple force-probe setup is employed to study the mechanical properties of
transcription activator-like effector (TALE) proteins in computer experiments.
It is shown that their spring-like arrangement benefits superelastic behaviour
which is manifested by large-scale global conformational changes along the
helical axis, thus linking structure and dynamics in TALE proteins. As
evidenced from the measured force-extension curves the dHax3 and PthXo1 TALEs
behave like linear springs, obeying Hooke's law, for moderate global structural
changes. For larger deformations, however, the proteins exhibit nonlinearities
and the structures become stiffer the more they are stretched. Flexibility is
not homogeneously distributed over TALE structure, but instead soft spots which
correspond to the RVD loop residues and present key agents in the transmission
of conformational motions are identified.


  Biological cells are often found to sense their chemical environment near the
single-molecule detection limit. Surprisingly, this precision is higher than
simple estimates of the fundamental physical limit, hinting towards active
sensing strategies. In this work, we analyse the effect of cell memory, e.g.
from slow biochemical processes, on the precision of sensing by cell-surface
receptors. We derive analytical formulas, which show that memory significantly
improves sensing in weakly fluctuating environments. However, surprisingly when
memory is adjusted dynamically, the precision is always improved, even in
strongly fluctuating environments. In support of this prediction we quantify
the directional biases in chemotactic Dictyostelium discoideum cells in a flow
chamber with alternating chemical gradients. The strong similarities between
cell sensing and control engineering suggest universal problem-solving
strategies of living matter.


  We investigated global patterns of variation in 157 whole genome sequences of
Vibrio parahaemolyticus, a free-living and seafood associated marine bacterium.
Pandemic clones, responsible for recent outbreaks of gastroenteritis in humans
have spread globally. However, there are oceanic gene pools, one located in the
oceans surrounding Asia and another in the Mexican Gulf. Frequent recombination
means that most isolates have acquired the genetic profile of their current
location. We investigated the genetic structure in the Asian gene pool by
calculating the effective population size in two different ways. Under standard
neutral models, the two estimates should give similar answers but we found a
thirty fold difference. We propose that this discrepancy is caused by the
subdivision of the species into a hundred or more ecotypes which are maintained
stably in the population. To investigate the genetic factors involved, we used
51 unrelated isolates to conduct a genome-wide scan for epistatically
interacting loci. We found a single example of strong epistasis between distant
genome regions. A majority of strains had a type VI secretion system associated
with bacterial killing. The remaining strains had genes associated with biofilm
formation and regulated by c-di-GMP signaling. All strains had one or other of
the two systems and none of isolate had complete complements of both systems,
although several strains had remnants. Further top-down analysis of patterns of
linkage disequilibrium within frequently recombining species will allow a
detailed understanding of how selection acts to structure the pattern of
variation within natural bacterial populations.


  Stochastic modelling provides an indispensable tool for understanding how
random events at the molecular level influence cellular functions. In practice,
the common challenge is to calibrate a large number of model parameters against
the experimental data. A related problem is to efficiently study how the
behaviour of a stochastic model depends on its parameters, i.e. whether a
change in model parameters can lead to a significant qualitative change in
model behaviour (bifurcation). In this paper, tensor-structured parametric
analysis (TPA) is presented. It is based on recently proposed low-parametric
tensor-structured representations of classical matrices and vectors. This
approach enables simultaneous computation of the model properties for all
parameter values within a parameter space. This methodology is exemplified to
study the parameter estimation, robustness, sensitivity and bifurcation
structure in stochastic models of biochemical networks. The TPA has been
implemented in Matlab and the codes are available at http://www.stobifan.org .


  Admixture models are a ubiquitous approach to capture latent population
structure in genetic samples. Despite the widespread application of admixture
models, little thought has been devoted to the quality of the model fit or the
accuracy of the estimates of parameters of interest for a particular study.
Here we develop methods for validating admixture models based on posterior
predictive checks (PPCs), a Bayesian method for assessing the quality of a
statistical model. We develop PPCs for five population-level statistics of
interest: within-population genetic variation, background linkage
disequilibrium, number of ancestral populations, between-population genetic
variation, and the downstream use of admixture parameters to correct for
population structure in association studies. Using PPCs, we evaluate the
quality of the model estimates for four qualitatively different population
genetic data sets: the POPRES European individuals, the HapMap phase 3
individuals, continental Indians, and African American individuals. We found
that the same model fitted to different genomic studies resulted in highly
study-specific results when evaluated using PPCs, illustrating the utility of
PPCs for model-based analyses in large genomic studies.


  Gene translation is the process in which intracellular macro-molecules,
called ribosomes, decode genetic information in the mRNA chain into the
corresponding proteins. Gene translation includes several steps. During the
elongation step, ribosomes move along the mRNA in a sequential manner and link
amino-acids together in the corresponding order to produce the proteins.
  The homogeneous ribosome flow model(HRFM) is a deterministic computational
model for translation-elongation under the assumption of constant elongation
rates along the mRNA chain. The HRFM is described by a set of n first-order
nonlinear ordinary differential equations, where n represents the number of
sites along the mRNA chain. The HRFM also includes two positive parameters:
ribosomal initiation rate and the (constant) elongation rate. In this paper, we
show that the steady-state translation rate in the HRFM is a concave function
of its parameters. This means that the problem of determining the parameter
values that maximize the translation rate is relatively simple. Our results may
contribute to a better understanding of the mechanisms and evolution of
translation-elongation. We demonstrate this by using the theoretical results to
estimate the initiation rate in M. musculus embryonic stem cell. The underlying
assumption is that evolution optimized the translation mechanism.
  For the infinite-dimensional HRFM, we derive a closed-form solution to the
problem of determining the initiation and transition rates that maximize the
protein translation rate. We show that these expressions provide good
approximations for the optimal values in the n-dimensional HRFM already for
relatively small values of n. These results may have applications for synthetic
biology where an important problem is to re-engineer genomic systems in order
to maximize the protein production rate.


  Social insect colonies can be seen as a distinct form of biological
organization because they function as superorganisms. Understanding how natural
selection acts on the emergence and maintenance of these colonies remains a
major question in evolutionary biology and ecology. Here, we explore this by
using multi-type branching processes to calculate the basic reproductive ratios
and the extinction probabilities for solitary versus eusocial reproductive
strategies. In order to derive precise mathematical results, we use a simple
haploid, asexual model. In general, we show that eusocial reproductive
strategies are unlikely to materialize unless large fitness advantages are
gained by the production of only a few workers. These fitness advantages are
maximized through obligate rather than facultative eusocial strategies.
Furthermore, we find that solitary reproduction is `unbeatable' as long as the
solitary reproductive ratio exceeds a critical value. In these cases, no
eusocial parameters exist that would reduce their probability of extinction.
Our results help to explain why the number of solitary species exceeds that of
eusocial ones: eusociality is a high risk, high reward strategy, while solitary
reproduction is better at risk management.


  Differences in transcriptional regulatory networks underlie much of the
phenotypic variation observed across organisms. Changes to cis-regulatory
elements are widely believed to be the predominant means by which regulatory
networks evolve, yet examples of regulatory network divergence due to
transcription factor (TF) variation have also been observed. To systematically
ascertain the extent to which TFs contribute to regulatory divergence, we
analyzed the evolution of the largest class of metazoan TFs, Cys2-His2 zinc
finger (C2H2-ZF) TFs, across 12 Drosophila species spanning ~45 million years
of evolution. Remarkably, we uncovered that a significant fraction of all
C2H2-ZF 1-to-1 orthologs in flies exhibit variations that can affect their
DNA-binding specificities. In addition to loss and recruitment of C2H2-ZF
domains, we found diverging DNA-contacting residues in ~47% of domains shared
between D. melanogaster and the other fly species. These diverging
DNA-contacting residues, found in ~66% of the D. melanogaster C2H2-ZF genes in
our analysis and corresponding to ~24% of all annotated D. melanogaster TFs,
show evidence of functional constraint: they tend to be conserved across
phylogenetic clades and evolve slower than other diverging residues. These same
variations were rarely found as polymorphisms within a population of D.
melanogaster flies, indicating their rapid fixation. The predicted
specificities of these dynamic domains gradually change across phylogenetic
distances, suggesting stepwise evolutionary trajectories for TF divergence.
Further, whereas proteins with conserved C2H2-ZF domains are enriched in
developmental functions, those with varying domains exhibit no functional
enrichments. Our work suggests that a subset of highly dynamic and largely
unstudied TFs are a likely source of regulatory variation in Drosophila and
other metazoans.


  Motivation: Estimation of bacterial community composition from a
high-throughput sequenced sample is an important task in metagenomics
applications. Since the sample sequence data typically harbors reads of
variable lengths and different levels of biological and technical noise,
accurate statistical analysis of such data is challenging. Currently popular
estimation methods are typically very time consuming in a desktop computing
environment.
  Results: Using sparsity enforcing methods from the general sparse signal
processing field (such as compressed sensing), we derive a solution to the
community composition estimation problem by a simultaneous assignment of all
sample reads to a pre-processed reference database. A general statistical model
based on kernel density estimation techniques is introduced for the assignment
task and the model solution is obtained using convex optimization tools.
Further, we design a greedy algorithm solution for a fast solution. Our
approach offers a reasonably fast community composition estimation method which
is shown to be more robust to input data variation than a recently introduced
related method.
  Availability: A platform-independent Matlab implementation of the method is
freely available at http://www.ee.kth.se/ctsoftware; source code that does not
require access to Matlab is currently being tested and will be made available
later through the above website.


  We investigate the emergence and persistence of communities through a
recently proposed mechanism of adaptive rewiring in coevolutionary networks. We
characterize the topological structures arising in a coevolutionary network
subject to an adaptive rewiring process and a node dynamics given by a simple
voterlike rule. We find that, for some values of the parameters describing the
adaptive rewiring process, a community structure emerges on a connected
network. We show that the emergence of communities is associated to a decrease
in the number of active links in the system, i.e. links that connect two nodes
in different states. The lifetime of the community structure state scales
exponentially with the size of the system. Additionally, we find that a small
noise in the node dynamics can sustain a diversity of states and a community
structure in time in a finite size system. Thus, large system size and/or local
noise can explain the persistence of communities and diversity in many real
systems.


  We perform extensive coarse-grained (CG) Langevin dynamics simulations of
intrinsically disordered proteins (IDPs), which possess fluctuating
conformational statistics between that for excluded volume random walks and
collapsed globules. Our CG model includes repulsive steric, attractive
hydrophobic, and electrostatic interactions between residues and is calibrated
to a large collection of single-molecule fluorescence resonance energy transfer
data on the inter-residue separations for 36 pairs of residues in five IDPs:
$\alpha$-, $\beta$-, and $\gamma$-synuclein, the microtubule-associated protein
$\tau$, and prothymosin $\alpha$. We find that our CG model is able to
recapitulate the average inter-residue separations regardless of the choice of
the hydrophobicity scale, which shows that our calibrated model can robustly
capture the conformational dynamics of IDPs. We then employ our model to study
the scaling of the radius of gyration with chemical distance in 11 known IDPs.
We identify a strong correlation between the distance to the dividing line
between folded proteins and IDPs in the mean charge and hydrophobicity space
and the scaling exponent of the radius of gyration with chemical distance along
the protein.


  The spatio-temporal arrangement of interacting populations often influences
the maintenance of species diversity and is a subject of intense research.
Here, we study the spatio-temporal patterns arising from the cyclic competition
between three species in two dimensions. Inspired by recent experiments, we
consider a generic metapopulation model comprising "rock-paper-scissors"
interactions via dominance removal and replacement, reproduction, mutations,
pair-exchange and hopping of individuals. By combining analytical and numerical
methods, we obtain the model's phase diagram near its Hopf bifurcation and
quantitatively characterize the properties of the spiraling patterns arising in
each phase. The phases characterizing the cyclic competition away far from the
Hopf bifurcation (at low mutation rate) are also investigated. Our analytical
approach relies on the careful analysis of the properties of the complex
Ginzburg-Landau equation derived through a controlled (perturbative) multiscale
expansion around the model's Hopf bifurcation. Our results allows us to clarify
when spatial "rock-paper-scissors" competition leads to stable spiral waves and
under which circumstances they are influenced by nonlinear mobility.


  Piecewise-Linear in Rates (PWLR) Lyapunov functions are introduced for a
class of Chemical Reaction Networks (CRNs). In addition to their simple
structure, these functions are robust with respect to arbitrary monotone
reaction rates, of which mass-action is a special case. The existence of such
functions ensures the convergence of trajectories towards equilibria, and
guarantee their asymptotic stability with respect to the corresponding
stoichiometric compatibility class. We give the definition of these Lyapunov
functions, prove their basic properties, and provide algorithms for
constructing them. Examples are provided, relationship with consensus dynamics
are discussed, and future directions are elaborated.


  We develop revised theoretical ideas on the mechanism by which the
transcription factor proteins locate their specific binding sites on DNA faster
than the three-dimensional (3D) diffusion controlled rate limit. We demonstrate
that the 3D-diffusion controlled rate limit can be enhanced when the protein
molecule reads several possible binding stretches of the template DNA via
one-dimensional (1D) diffusion upon each 3D-diffusion mediated collision or
nonspecific binding event. The overall enhancement of site-specific association
rate is directly proportional to the maximum possible sliding length (LA,
square root of (6Do/kr) where Do is the 1D-diffusion coefficient and kr is the
dissociation rate constant associated with the nonspecific DNA-protein complex)
associated with the 1D-diffusion of protein molecule along DNA. Upon
considering several possible mechanisms we find that the DNA binding proteins
can efficiently locate their cognate sites on DNA by switching across
fast-moving, slow-moving and reading states of their DNA binding domains in a
cyclic manner. Irrespective of the type of mechanism the overall rate
enhancement factor asymptotically approaches a limiting value which is directly
proportional to LA as the total length of DNA that contains the cognate site
increases. These results are consistent with the in vitro experimental
observations.


  The diffusive arrival of transcription factors at the promoter sites on the
DNA sets a lower bound on how accurately a cell can regulate its protein
levels. Using results from the literature on diffusion-influenced reactions, we
derive an analytical expression for the lower bound on the precision of
transcriptional regulation. In our theory, transcription factors can perform
multiple rounds of 1D diffusion along the DNA and 3D diffusion in the cytoplasm
before binding to the promoter. Comparing our expression for the lower bound on
the precision against results from Green's Function Reaction Dynamics
simulations shows that the theory is highly accurate under biologically
relevant conditions. Our results demonstrate that, to an excellent
approximation, the promoter switches between the transcription-factor bound and
unbound state in a Markovian fashion. This remains true even in the presence of
sliding, i.e. with 1D diffusion along the DNA. This has two important
implications: (1) minimizing the noise in the promoter state is equivalent to
minimizing the search time of transcription factors for their promoters; (2)
the complicated dynamics of 3D diffusion in the cytoplasm and 1D diffusion
along the DNA can be captured in a well-stirred model by renormalizing the
promoter association and dissociation rates, making it possible to efficiently
simulate the promoter dynamics using Gillespie simulations. Based on the recent
experimental observation that sliding can speed up the promoter search by a
factor of 4, our theory predicts that sliding can enhance the precision of
transcriptional regulation by a factor of 2.


  Complex social behaviors lie at the heart of many of the challenges facing
evolutionary biology, sociology, economics, and beyond. For evolutionary
biologists in particular the question is often how such behaviors can arise
\textit{de novo} in a simple evolving system. How can group behaviors such as
collective action, or decision making that accounts for memories of past
experience, emerge and persist? Evolutionary game theory provides a framework
for formalizing these questions and admitting them to rigorous study. Here we
develop such a framework to study the evolution of sustained collective action
in multi-player public-goods games, in which players have arbitrarily long
memories of prior rounds of play and can react to their experience in an
arbitrary way. To study this problem we construct a coordinate system for
memory-$m$ strategies in iterated $n$-player games that permits us to
characterize all the cooperative strategies that resist invasion by any mutant
strategy, and thus stabilize cooperative behavior. We show that while larger
games inevitably make cooperation harder to evolve, there nevertheless always
exists a positive volume of strategies that stabilize cooperation provided the
population size is large enough. We also show that, when games are small,
longer-memory strategies make cooperation easier to evolve, by increasing the
number of ways to stabilize cooperation. Finally we explore the co-evolution of
behavior and memory capacity, and we find that longer-memory strategies tend to
evolve in small games, which in turn drives the evolution of cooperation even
when the benefits for cooperation are low.


  Competitions can occur on an absolute scale, to be faster or more efficient,
or they can occur on a relative scale, to "beat" one's competitor in a zero-sum
game. Ecological models have focused on absolute competitions, in which optima
exist. Classic evolutionary models such as the Wright-Fisher model, as well as
more recent models of travelling waves, have focused on purely relative
competitions, in which fitness continues to increase indefinitely, without
actually progressing anywhere. This manuscript proposes a new way to describe
both at the same time. It begins with a revised version of r/K-selection
theory. r continues to describe maximum reproductive speed, but the new version
of K, with a different subscript, now describes parsimoniousness in territory
use, a group-selected, anti-tragedy-of-the-commons trait. A third dimension c
of fitness is then added to this novel system, one which is unitless and
normalized, and hence capable of capturing the population genetics concept w of
a strictly relative, genetically-limited competitive race. MacArthur's original
version of r/K-selection theory is shown to confound parsimoniousness K with
competitive ability c, despite the fact that available data suggests a negative
correlation between the two; here they are disentangled. A rotation of the
resulting three-dimensional system provides a population genetic underpinning
for Grime's universal adaptive strategy theory of ruderals (selected for high
r), stress tolerators (selected for a combination of high r and high K), and
competitors (selected for a combination of high r and high c).


  Due to the scarcity of quantitative details about biological phenomena,
quantitative modeling in systems biology can be compromised, especially at the
subcellular scale. One way to get around this is qualitative modeling because
it requires few to no quantitative information. One of the most popular
qualitative modeling approaches is the Boolean network formalism. However,
Boolean models allow variables to take only two values, which can be too
simplistic in some cases. The present work proposes a modeling approach derived
from Boolean networks where continuous logical operators are used and where
edges can be tuned. Using continuous logical operators allows variables to be
more finely valued while remaining qualitative. To consider that some
biological interactions can be slower or weaker than other ones, edge states
are also computed in order to modulate in speed and strength the signal they
convey. The proposed formalism is illustrated on a toy network coming from the
epidermal growth factor receptor signaling pathway. The obtained simulations
show that continuous results are produced, thus allowing finer analysis. The
simulations also show that modulating the signal conveyed by the edges allows
to incorporate knowledge about the interactions they model. The goal is to
provide enhancements in the ability of qualitative models to simulate the
dynamics of biological networks while limiting the need of quantitative
information.


  Double phosphorylation of protein kinases is a common feature of signalling
cascades. This motif may reduce cross-talk between signalling pathways, as the
second phosphorylation site allows for proofreading, especially when
phosphorylation is distributive rather than processive. Recent studies suggest
that phosphorylation can be `pseudo-processive' in the crowded cellular
environment, as rebinding after the first phosphorylation is enhanced by slow
diffusion. Here, we use a simple model with unsaturated reactants to show that
specificity for one substrate over another drops as rebinding increases and
pseudo-processive behavior becomes possible. However, this loss of specificity
with increased rebinding is typically also observed if two distinct enzyme
species are required for phosphorylation, i.e. when the system is necessarily
distributive. Thus the loss of specificity is due to an intrinsic reduction in
selectivity with increased rebinding, which benefits inefficient reactions,
rather than pseudo-processivity itself. We also show that proofreading can
remain effective when the intended signalling pathway exhibits high levels of
rebinding-induced pseudo-processivity, unlike other proposed advantages of the
dual phosphorylation motif.


  Motivation: Building the histogram of occurrences of every $k$-symbol long
substring of nucleotide data is a standard step in many bioinformatics
applications, known under the name of $k$-mer counting. Its applications
include developing de Bruijn graph genome assemblers, fast multiple sequence
alignment and repeat detection. The tremendous amounts of NGS data require fast
algorithms for $k$-mer counting, preferably using moderate amounts of memory.
  Results: We present a novel method for $k$-mer counting, on large datasets at
least twice faster than the strongest competitors (Jellyfish~2, KMC~1), using
about 12\,GB (or less) of RAM memory. Our disk-based method bears some
resemblance to MSPKmerCounter, yet replacing the original minimizers with
signatures (a carefully selected subset of all minimizers) and using $(k,
x)$-mers allows to significantly reduce the I/O, and a highly parallel overall
architecture allows to achieve unprecedented processing speeds. For example,
KMC~2 allows to count the 28-mers of a human reads collection with 44-fold
coverage (106\,GB of compressed size) in about 20 minutes, on a 6-core Intel i7
PC with an SSD.
  Availability: KMC~2 is freely available at http://sun.aei.polsl.pl/kmc.
  Contact: sebastian.deorowicz@polsl.pl


  Some species exhibit very high levels of DNA sequence variability; there is
also evidence for the existence of heritable epigenetic variants that
experience state changes at a much higher rate than sequence variants. In both
cases, the resulting high diversity levels within a population (hyperdiversity)
mean that standard population genetics methods are not trustworthy. We analyze
a population genetics model that incorporates purifying selection, reversible
mutations and genetic drift, assuming a stationary population size. We derive
analytical results for both population parameters and sample statistics, and
discuss their implications for studies of natural genetic and epigenetic
variation. In particular, we find that (1) many more intermediate frequency
variants are expected than under standard models, even with moderately strong
purifying selection (2) rates of evolution under purifying selection may be
close to, or even exceed, neutral rates. These findings are related to
empirical studies of sequence and epigenetic variation.


  We study the effect of membrane proteins on the shape, composition and
thermodynamic stability of the surrounding membrane. When the coupling between
membrane composition and curvature is strong enough the nearby composition and
shape both undergo a transition from over-damped to under-damped spatial
variation, well before the membrane becomes unstable in the bulk. This
transition is associated with a change in the sign of the thermodynamic energy
and hence has the unusual features that it can favour the early stages of coat
assembly necessary for vesiculation (budding), while suppressing the activity
of mechanosensitive membrane channels and transporters. Our results also
suggest an approach to obtain physical parameters that are otherwise difficult
to measure.


  In this perspective we describe the critical role membranes play in
modulating the structures of the Amyloid Precursor Proteins to produce the
peptides involved in the Alzheimer's disease. Some of the key concepts related
to protein aggregation including the potential role of the excited states of
monomers in initiating protein aggregation are described.


  Estimation of epidemiological and population parameters from molecular
sequence data has become central to the understanding of infectious disease
dynamics. Various models have been proposed to infer details of the dynamics
that describe epidemic progression. These include inference approaches derived
from Kingman's coalescent theory. Here, we use recently described coalescent
theory for epidemic dynamics to develop stochastic and deterministic coalescent
SIR tree priors. We implement these in a Bayesian phylogenetic inference
framework to permit joint estimation of SIR epidemic parameters and the sample
genealogy. We assess the performance of the two coalescent models and also
juxtapose results obtained with BDSIR, a recently published
birth-death-sampling model for epidemic inference. Comparisons are made by
analyzing sets of genealogies simulated under precisely known epidemiological
parameters. Additionally, we analyze influenza A (H1N1) sequence data sampled
in the Canterbury region of New Zealand and HIV-1 sequence data obtained from
known UK infection clusters. We show that both coalescent SIR models are
effective at estimating epidemiological parameters from data with large
fundamental reproductive number $R_0$ and large population size $S_0$.
Furthermore, we find that the stochastic variant generally outperforms its
deterministic counterpart in terms of error, bias, and highest posterior
density coverage, particularly for smaller $R_0$ and $S_0$. However, each of
these inference models are shown to have undesirable properties in certain
circumstances, especially for epidemic outbreaks with $R_0$ close to one or
with small effective susceptible populations.


  Protein aggregation, linked to many of diseases, is initiated when monomers
access rogue conformations that are poised to form amyloid fibrils. We show,
using simulations of src SH3 domain, that mechanical force enhances the
population of the aggregation prone ($N^*$) states, which are rarely populated
under force free native conditions, but are encoded in the spectrum of native
fluctuations. The folding phase diagrams of SH3 as a function of denaturant
concentration ($[C]$), mechanical force ($f$), and temperature exhibit an
apparent two-state behavior, without revealing the presence of the elusive
$N^*$ states. Interestingly, the phase boundaries separating the folded and
unfolded states at all [C] and $f$ fall on a master curve, which can can be
quantitatively described using an analogy to superconductors in a magnetic
field. The free energy profiles as a function of the molecular extension ($R$),
which are accessible in pulling experiments, ($R$), reveal the presence of a
native-like $N^*$ with a disordered solvent-exposed amino terminal
$\beta$-strand. The structure of the $N^*$ state is identical to that found in
Fyn SH3 by NMR dispersion experiments. We show that the time scale for fibril
formation can be estimated from the population of the $N^*$ state, determined
by the free energy gap separating the native structure and the $N^*$ state, a
finding that can be used to assess fibril forming tendencies of proteins. The
structures of the $N^*$ state are used to show that oligomer formation and
likely route to fibrils occur by a domain-swap mechanism in SH3 domain.


  The human microbiome is the ensemble of genes in the microbes that live
inside and on the surface of humans. Because microbial sequencing information
is now much easier to come by than phenotypic information, there has been an
explosion of sequencing and genetic analysis of microbiome samples. Much of the
analytical work for these sequences involves phylogenetics, at least
indirectly, but methodology has developed in a somewhat different direction
than for other applications of phylogenetics. In this paper I review the field
and its methods from the perspective of a phylogeneticist, as well as
describing current challenges for phylogenetics coming from this type of work.


  We study the two-species symbiotic contact process (2SCP), recently proposed
in [de Oliveira, Santos and Dickman, Phys. Rev. E {\bf 86}, 011121 (2012)] . In
this model, each site of a lattice may be vacant or host single individuals of
species A and/or B. Individuals at sites with both species present interact in
a symbiotic manner, having a reduced death rate, $\mu < 1$. Otherwise, the
dynamics follows the rules of the basic CP, with individuals reproducing to
vacant neighbor sites at rate $\lambda$ and dying at a rate of unity. We
determine the full phase diagram in the $\lambda-\mu$ plane in one and two
dimensions by means of exact numerical quasistationary distributions, cluster
approximations, and Monte Carlo simulations. We also study the effects of
asymmetric creation rates and diffusion of individuals. In two dimensions, for
sufficiently strong symbiosis (i.e., small $\mu$), the absorbing-state phase
transition becomes discontinuous for diffusion rates $D$ within a certain
range. We report preliminary results on the critical surface and tricritical
line in the $\lambda-\mu-D$ space. Our results raise the possibility that
strongly symbiotic associations of mobile species may be vulnerable to sudden
extinction under increasingly adverse conditions.


  Fractal structure of shortest paths depends strongly on interresidue
interaction cutoff distance. The dimensionality of shortest paths is calculated
as a function of interaction cutoff distance. Shortest paths are self similar
with a fractal dimension of 1.12 when calculated with step lengths larger than
6.8 {\AA}. Paths are multifractal below 6.8 {\AA}. The number of steps to
traverse a shortest path is a discontinuous function of cutoff size at short
cutoff values, showing abrupt decreases to smaller values as cutoff distance
increases. As information progresses along the direction of a shortest path a
large set of residues are affected because they are interacting neighbors to
the residues of the shortest path. Thus, several residues are involved
diffusively in information transport which may be identified with the present
model. An algorithm is introduced to determine the residues of a given shortest
path. The shortest path residues are the highly visited residues during
information transport. These paths are shown to lie on the high entropy
landscape of the protein where entropy is taken to increase with abundance of
visits to nodes during signal transport.


  Background: To understand individual genomes it is necessary to look at the
variations that lead to changes in phenotype and possibly to disease. However,
genotype information alone is often not sufficient and additional knowledge
regarding the phase of the variation is needed to make correct interpretations.
Interactive visualizations, that allow the user to explore the data in various
ways, can be of great assistance in the process of making well informed
decisions. But, currently there is a lack for visualizations that are able to
deal with phased haplotype data. Results: We present inPHAP, an interactive
visualization tool for genotype and phased haplotype data. inPHAP features a
variety of interaction possibilities such as zooming, sorting, filtering and
aggregation of rows in order to explore patterns hidden in large genetic data
sets. As a proof of concept, we apply inPHAP to the phased haplotype data set
of Phase 1 of the 1000 Genomes Project. Thereby, inPHAP's ability to show
genetic variations on the population as well as on the individuals level is
demonstrated for several disease related loci. Conclusions: As of today, inPHAP
is the only visual analytical tool that allows the user to explore unphased and
phased haplotype data interactively. Due to its highly scalable design, inPHAP
can be applied to large datasets with up to 100 GB of data, enabling users to
visualize even large scale input data. inPHAP closes the gap between common
visualization tools for unphased genotype data and introduces several new
features, such as the visualization of phased data.


  The last two decades have seen a surge in kinetic and macroscopic models
derived to investigate the multi-scale aspects of self-organised biological
aggregations. Because the individual-level details incorporated into the
kinetic models (e.g., individual speeds and turning rates) make them somewhat
difficult to investigate, one is interested in transforming these models into
simpler macroscopic models, by using various scaling techniques that are
imposed by the biological assumptions of the models. Here, we consider three
scaling approaches (parabolic, hydrodynamic and grazing collision limits) that
can be used to reduce a class of non-local 1D and 2D models for biological
aggregations to simpler models existent in the literature. Next, we investigate
how some of the spatio-temporal patterns exhibited by the original kinetic
models are preserved via these scalings. To this end, we focus on the parabolic
scaling for non-local 1D models and apply asymptotic preserving numerical
methods, which allow us to analyse changes in the patterns as the scaling
coefficient $\epsilon$ is varied from $\epsilon=1$ (for 1D transport models) to
$\epsilon=0$ (for 1D parabolic models). We show that some patterns (describing
stationary aggregations) are preserved in the limit $\epsilon\to 0$, while
other patterns (describing moving aggregations) are lost in this limit. To
understand the loss of these patterns, we construct bifurcation diagrams.


  The Fenna Mathews Olson (FMO) complex of green sulphur bacteria is an example
of a photosynthetic pigment protein complex, in which the electronic properties
of the pigments are modified by the protein environment to promote efficient
excitonic energy transfer from antenna complexes to the reaction centres. Many
of the electronic properties of the FMO complex can be extracted from knowledge
of the static crystal structure. However, the recent observation and analysis
of long lasting quantum dynamics in the FMO complex point to protein dynamics
as a key factor in protecting and generating quantum coherence under laboratory
conditions. While fast inter and intra molecular vibrations have been
investigated extensively, the slow dynamics which effectively determine the
optical inhomogeneous broadening of experimental ensembles has received less
attention. Our study employs constrained geometric dynamics to study the
flexibility in the protein network by efficiently generating the accessible
conformational states from the published crystal structure. Statistical and
principle component analysis reveal highly correlated low frequency motions
between functionally relevant elements, including strong correlations between
pigments that are excitonically coupled. Our analysis reveals a hierarchy of
structural interactions which enforce these correlated motions, from the level
of monomer monomer interfaces right down to the alpha helices, beta sheets and
pigments. In addition to inducing strong spatial correlations across the
conformational ensemble, we find that the overall rigidity of the FMO complex
is exceptionally high. We suggest that these observations support the idea of
highly correlated inhomogeneous disorder of the electronic excited states,
which is further supported by the remarkably low variance of the excitonic
couplings of the conformational ensemble.


  Group-level cognitive states are widely observed in human social systems, but
their discussion is often ruled out a priori in quantitative approaches. In
this paper, we show how reference to the irreducible mental states and
psychological dynamics of a group is necessary to make sense of large scale
social phenomena. We introduce the problem of mental boundaries by reference to
a classic problem in the evolution of cooperation. We then provide an explicit
quantitative example drawn from ongoing work on cooperation and conflict among
Wikipedia editors, showing how some, but not all, effects of individual
experience persist in the aggregate. We show the limitations of methodological
individualism, and the substantial benefits that come from being able to refer
to collective intentions, and attributions of cognitive states of the form
"what the group believes" and "what the group values".


  Substantial research on structured sparsity has contributed to analysis of
many different applications. However, there have been few Bayesian procedures
among this work. Here, we develop a Bayesian model for structured sparsity that
uses a Gaussian process (GP) to share parameters of the sparsity-inducing prior
in proportion to feature similarity as defined by an arbitrary positive
definite kernel. For linear regression, this sparsity-inducing prior on
regression coefficients is a relaxation of the canonical spike-and-slab prior
that flattens the mixture model into a scale mixture of normals. This prior
retains the explicit posterior probability on inclusion parameters---now with
GP probit prior distributions---but enables tractable computation via
elliptical slice sampling for the latent Gaussian field. We motivate
development of this prior using the genomic application of association mapping,
or identifying genetic variants associated with a continuous trait. Our
Bayesian structured sparsity model produced sparse results with substantially
improved sensitivity and precision relative to comparable methods. Through
simulations, we show that three properties are key to this improvement: i)
modeling structure in the covariates, ii) significance testing using the
posterior probabilities of inclusion, and iii) model averaging. We present
results from applying this model to a large genomic dataset to demonstrate
computational tractability.


  In our paper selected linguistic features of genomes to study the statistics
of the gene codes are considered. We present the information theory from which
it follows that if the system is described by distributions of hyperbolic type
it leads to the possibility of entropy loss and stability. We show that the
histograms of gene lengths are similar to that of language words. We show the
correspondence between presented theory and results for the number of
replicated genes and replicated fragments of genes in genomes for Borelia
burgdorferi, Escherichia coli and Saccharomyces cerevisiae S288c.


  Chikungunya fever (CHIKV), a viral disease transmitted by mosquitoes, is
currently affecting several territories in the Caribbean. The vector is found
in the Americas from South Florida to Brazil and the Caribbean is a highly
connected region. There is therefore a significant risk for the epidemic to
quickly expand to a wide area in the Americas, and of introduction to Europe
where competent vectors are present. Here, we describe the spread of CHIKV in
the first three territories to report cases and between territories of the
region. Local transmission of CHIKV in the Caribbean is found to be very
effective (mean number of cases generated by a human case in the range 2 to 4).
There is a strong spatial signature in the regional epidemic, with the risk of
transmission between territories estimated to be inversely proportional to the
distance rather than being driven by air transportation. So far, this simple
model has successfully predicted observed patterns of spread. The spatial
structure allows ranking territories according to their current risk of
invasion. This characterization may help national and international agencies to
optimize resource allocation for monitoring and control and encourage
territories with elevated risks to act. This is important to reduce the risk of
dissemination to other parts of the world, including Europe.


  Understanding the assembly of ecosystems to estimate the number of species at
different spatial scales is a challenging problem. Until now, maximum entropy
approaches have lacked the important feature of considering space in an
explicit manner. We propose a spatially explicit maximum entropy model suitable
to describe spatial patterns such as the species area relationship and the
endemic area relationship. Starting from the minimal information extracted from
presence/absence data, we compare the behavior of two models considering the
occurrence or lack thereof of each species and information on spatial
correlations. Our approach uses the information at shorter spatial scales to
infer the spatial organization at larger ones. We also hypothesize a possible
ecological interpretation of the effective interaction we use to characterize
spatial clustering.


  The differences between DNA-sequences within a population are the basis to
infer the ancestral relationship of the individuals. Within the classical
infinitely many sites model, it is possible to estimate the mutation rate based
on the site frequency spectrum, which is comprised by the numbers
$C_1,...,C_{n-1}$, where n is the sample size and $C_s$ is the number of site
mutations (Single Nucleotide Polymorphisms, SNPs) which are seen in $s$
genomes. Classical results can be used to compare the observed site frequency
spectrum with its neutral expectation, $E[C_s]= \theta_2/s$, where $\theta_2$
is the scaled site mutation rate. In this paper, we will relax the assumption
of the infinitely many sites model that all individuals only carry homologous
genetic material. Especially, it is today well-known that bacterial genomes
have the ability to gain and lose genes, such that every single genome is a
mosaic of genes, and genes are present and absent in a random fashion, giving
rise to the dispensable genome. While this presence and absence has been
modeled under neutral evolution within the infinitely many genes model in
previous papers, we link presence and absence of genes with the numbers of site
mutations seen within each gene. In this work we derive a formula for the
expectation of the joint gene and site frequency spectrum, denotes $G_{k,s}$
the number of mutated sites occurring in exactly $s$ gene sequences, while the
corresponding gene is present in exactly $k$ individuals. We show that standard
estimators of $\theta_2$ for dispensable genes are biased and that the site
frequency spectrum for dispensable genes differs from the classical result.


  Markov chains are a common framework for individual-based state and time
discrete models in ecology and evolution. Their use, however, is largely
limited to systems with a low number of states, since the transition matrices
involved pose considerable challenges as their size and their density increase.
Big, dense transition matrices may easily defy both the computer's memory and
the scientists' ability to interpret them, due to the very high amount of
information they contain; yet approximations using other types of models are
not always the best solution.
  We propose a set of methods to overcome the difficulties associated with big,
dense Markov chain transition matrices. Using a population genetic model as an
example, we demonstrate how big matrices can be transformed into clear and
easily interpretable graphs with the help of network analysis. Moreover, we
describe an algorithm to save computer memory by substituting the original
matrix with a sparse approximate while preserving all its mathematically
important properties. In the same model example, we manage to store about 90%
less data while keeping more than 99% of the information contained in the
matrix and a closely corresponding dominant eigenvector.
  Our approach is an example how numerical limitations for the number of states
in a Markov chain can be overcome. By facilitating the use of state-rich Markov
chain models, they may become a valuable supplement to the diversity of models
currently employed in biology.


  The stage of evolution is the population of reproducing individuals. The
structure of the population is know to affect the dynamics and outcome of
evolutionary processes, but analytical results for generic random structures
have been lacking. The most general result so far, the isothermal theorem,
assumes the propensity for change in each position is exactly the same, but
realistic biological structures are always subject to variation and noise. We
consider a population of finite size $n$ under constant selection whose
structure is given by a wide variety of weighted, directed, random graphs;
vertices represent individuals and edges interactions between individuals. By
establishing a robustness result for the isothermal theorem and using large
deviation estimates to understand the typical structure of random graphs, we
prove that for a generalization of the Erd\H{o}s-R\'{e}nyi model the fixation
probability of an invading mutant is approximately the same as that of a mutant
of equal fitness in a well-mixed population with high probability. Simulations
of perturbed lattices, small-world networks, and scale-free networks behave
similarly. We conjecture that the fixation probability in a well-mixed
population, $(1-r^{-1})/(1-r^{-n})$, is universal: for many random graph
models, the fixation probability approaches the above function uniformly as the
graphs become large.


  Two different zinc oxide nanoparticles, as well as zinc ions, are used to
study the cellular responses of the RAW 264 macrophage cell line. A proteomic
screen is used to provide a wide view of the molecular effects of zinc, and the
most prominent results are cross-validated by targeted studies. Furthermore,
the alteration of important macrophage functions (e.g. phagocytosis) by zinc is
also investigated. The intracellular dissolution/uptake of zinc is also studied
to further characterize zinc toxicity. Zinc oxide nanoparticles dissolve
readily in the cells, leading to high intracellular zinc concentrations, mostly
as protein-bound zinc. The proteomic screen reveals a rather weak response in
the oxidative stress response pathway, but a strong response both in the
central metabolism and in the proteasomal protein degradation pathway. Targeted
experiments confirm that carbohydrate catabolism and proteasome are critical
determinants of sensitivity to zinc, which also induces DNA damage. Conversely,
glutathione levels and phagocytosis appear unaffected at moderately toxic zinc
concentrations.


  Depletion forces play a role in the compaction and de-compation of
chromosomal material in simple cells but it remains debatable whether they are
sufficient to account for chromosomal collapse. We present coarse-grained
molecular dynamics simulations, which reveal that depletion-induced attraction
is sufficient to cause the collapse of a flexible chain of large structural
monomers immersed in a bath of smaller depletants. These simulations use an
explicit coarse-grained computational model that treats both the supercoiled
DNA structural monomers and the smaller protein crowding agents as
combinatorial, truncated Lennard-Jones spheres. By presenting a simple
theoretical model, we quantitatively cast the action of depletants on
supercoiled bacterial DNA as an effective solvent quality. The rapid collapse
of the simulated flexible chromosome at the predicted volume fraction of
depletants is a continous phase transition. Additional physical effects to such
simple chromosome models, such as enthalpic interactions between structural
monomers or chain rigidity, are required if the collapse is to be a first-order
phase transition.


  Polynomial dynamical systems are widely used to model and study real
phenomena. In biochemistry, they are the preferred choice for modelling the
concentration of chemical species in reaction networks with mass-action
kinetics. These systems are typically parameterised by many (unknown)
parameters. A goal is to understand how properties of the dynamical systems
depend on the parameters. Qualitative properties relating to the behaviour of a
dynamical system are locally inferred from the system at steady state. Here we
focus on steady states that are the positive solutions to a parameterised
system of generalised polynomial equations. In recent years, methods from
computational algebra have been developed to understand these solutions, but
our knowledge is limited: for example, we cannot efficiently decide how many
positive solutions the system has as a function of the parameters. Even
deciding whether there is one or more solutions is non-trivial. We present a
new method, based on so-called injectivity, to preclude or assert that multiple
positive solutions exist. The results apply to generalised polynomials and
variables can be restricted to the linear, parameter-independent first
integrals of the dynamical system. The method has been tested in a wide range
of systems.


  Analysis of probability distributions conditional on species trees has
demonstrated the existence of anomalous ranked gene trees (ARGTs), ranked gene
trees that are more probable than the ranked gene tree that accords with the
ranked species tree. Here, to improve the characterization of ARGTs, we study
enumerative and probabilistic properties of two classes of ranked labeled
species trees, focusing on the presence or avoidance of certain subtree
patterns associated with the production of ARGTs. We provide exact enumerations
and asymptotic estimates for cardinalities of these sets of trees, showing that
as the number of species increases without bound, the fraction of all ranked
labeled species trees that are ARGT-producing approaches 1. This result extends
beyond earlier existence results to provide a probabilistic claim about the
frequency of ARGTs.


  We are interested in the impact of natural selection in a prey-predator
community. We introduce an individual-based model of the community that takes
into account both prey and predator phenotypes. Our aim is to understand the
phenotypic coevolution of prey and predators. The community evolves as a
multi-type birth and death process with mutations. We first consider the
infinite particle approximation of the process without mutation. In this limit,
the process can be approximated by a system of differential equations. We prove
the existence of a unique globally asymptotically stable equilibrium under
specific conditions on the interaction among prey individuals. When mutations
are rare, the community evolves on the mutational scale according to a
Markovian jump process. This process describes the successive equilibria of the
prey-predator community and extends the Polymorphic Evolutionary Sequence to a
coevolutionary framework. We then assume that mutations have a small impact on
phenotypes and consider the evolution of monomorphic prey and predator
populations. The limit of small mutation steps leads to a system of two
differential equations which is a version of the canonical equation of adaptive
dynamics for the prey-predator coevolution. We illustrate these different
limits with an example of prey-predator community that takes into account
different prey defense mechanisms. We observe through simulations how these
various prey strategies impact the community.


  Competition between individuals drives the evolution of whole species.
Although the fittest individuals survive the longest and produce the most
offspring, in some circumstances the resulting species may not be optimally
fit. Here, using theoretical analysis and stochastic simulations of a simple
model ecology, we show how the mode of competition can profoundly affect the
fitness of evolved species. When individuals compete directly with one another,
the adaptive dynamics framework provides accurate predictions for the number
and distribution of species, which occupy positions of maximal fitness. By
contrast, if competition is mediated by the consumption of a common resource
then demographic noise leads to the stabilization of species with near minimal
fitness.


  We propose a mechanism in which two molecular knots pass through each other
and swap positions along a polymer strand. Associated free energy barriers in
our simulations only amount to a few $k_{B}T$, which may enable the interchange
of knots on a single DNA strand.


  Developing robust, quantitative methods to optimize resource allocations in
response to epidemics has the potential to save lives and minimize health care
costs. In this paper, we develop and apply a computationally efficient
algorithm that enables us to calculate the complete probability distribution
for the final epidemic size in a stochastic Susceptible-Infected-Recovered
(SIR) model. Based on these results, we determine the optimal allocations of a
limited quantity of vaccine between two non-interacting populations. We compare
the stochastic solution to results obtained for the traditional, deterministic
SIR model. For intermediate quantities of vaccine, the deterministic model is a
poor estimate of the optimal strategy for the more realistic, stochastic case.


  Analytic and computational methods developed within statistical physics have
found applications in numerous disciplines. In this letter, we use such methods
to solve a long-standing problem in statistical genetics. The problem, posed by
Haldane and Waddington [J.B.S. Haldane and C.H. Waddington, Genetics 16,
357-374 (1931)], concerns so-called recombinant inbred lines (RILs) produced by
repeated inbreeding. Haldane and Waddington derived the probabilities of RILs
when considering 2 and 3 genes but the case of 4 or more genes has remained
elusive. Our solution uses two probabilistic frameworks relatively unknown
outside of physics: Glauber's formula and self-consistent equations of the
Schwinger-Dyson type. Surprisingly, this combination of statistical formalisms
unveils the exact probabilities of RILs for any number of genes. Extensions of
the framework may have applications in population genetics and beyond.


  DNA barcoding has emerged as a cost-effective approach for species
identification. However, the scarcity of tools used for searching the booming
reference database becomes an obstacle, currently with BLAST as the only
practical choice. Here, we propose a program - LV Barcoding - based on both the
random hyperplane projection-based locality sensitive hashing method and the
composition vector-based VIP Barcoding for fast species identification. The
performance of LV Barcoding is assessed on the data release of BOLD. LV
Barcoding has higher accuracy than BLAST, and is able to match a single query
against ~114,000 reference barcodes within 10 seconds on a desktop computer.
This program is available at http://msl.sls.cuhk.edu.hk/vipbarcoding/.


  Because mutations are mostly deleterious, mutation rates should be reduced by
natural selection. However, mutations also provide the raw material for
adaptation. Therefore, evolutionary theory suggests that the mutation rate must
balance between adaptability - the ability to adapt - and adaptedness - the
ability to remain adapted. We model an asexual population crossing a fitness
valley and analyze the rate of complex adaptation with and without
stress-induced mutagenesis - the increase of mutation rates in response to
stress or maladaptation. We show that stress-induced mutagenesis increases the
rate of complex adaptation without reducing the population mean fitness, thus
breaking the evolutionary trade-off between adaptability and adaptedness. Our
theoretical results support the hypothesis that stress-induced mutagenesis
promotes adaptation and provide quantitative predictions of the rate of complex
adaptation with different mutational strategies.


  Cells adapt to different conditions by altering a vast number of components,
which is measurable using transcriptome analysis. Given that a cell undergoing
steady growth is constrained to sustain each of its internal components, the
abundance of all the components in the cell has to be roughly doubled during
each cell division event. From this steady-growth constraint, expression of all
genes is shown to change along a one-parameter curve in the state space in
response to the environmental stress. This leads to a global relationship that
governs the cellular state: By considering a relatively moderate change around
a steady state, logarithmic changes in expression are shown to be proportional
across all genes, upon alteration of stress strength, with the proportionality
coefficient given by the change in the growth rate of the cell. This theory is
confirmed by transcriptome analysis of Escherichia Coli in response to several
stresses.