  DNA unzipping, the separation of its double helix into single strands, is
crucial in modulating a host of genetic processes. Although the large-scale
separation of double-stranded DNA has been studied with a variety of
theoretical and experimental techniques, the minute details of the very first
steps of unzipping are still unclear. Here, we use atomistic molecular dynamics
(MD) simulations, coarse-grained simulations and a statistical-mechanical model
to study the initiation of DNA unzipping by an external force. The calculation
of the potential of mean force profiles for the initial separation of the first
few terminal base pairs in a DNA oligomer reveal that forces ranging between
130 and 230 pN are needed to disrupt the first base pair, values of an order of
magnitude larger than those needed to disrupt base pairs in partially unzipped
DNA. The force peak has an "echo," of approximately 50 pN, at the distance that
unzips the second base pair. We show that the high peak needed to initiate
unzipping derives from a free energy basin that is distinct from the basins of
subsequent base pairs because of entropic contributions and we highlight the
microscopic origin of the peak. Our results suggest a new window of exploration
for single molecule experiments.


  Physical properties of capsids of plant and animal viruses are important
factors in capsid self-assembly, survival of viruses in the extracellular
environment, and their cell infectivity. Virus shells can have applications as
nanocontainers and delivery vehicles in biotechnology and medicine. Combined
AFM experiments and computational modeling on sub-second timescales of the
indentation nanomechanics of Cowpea Chlorotic Mottle Virus (CCMV) capsid show
that the capsid's physical properties are dynamic and local characteristics of
the structure, which depend on the magnitude and geometry of mechanical input.
Surprisingly, under large deformations the CCMV capsid transitions to the
collapsed state without substantial local structural alterations. The enthalpy
change in this deformation state dH = 11.5 - 12.8 MJ/mol is mostly due to
large-amplitude out-of-plane excitations, which contribute to the capsid
bending, and the entropy change TdS = 5.1 - 5.8 MJ/mol is mostly due to
coherent in-plane rearrangements of protein chains, which result in the capsid
stiffening. Dynamic coupling of these modes defines the extent of elasticity
and reversibility of capsid mechanical deformation. This emerging picture
illuminates how unique physico-chemical properties of protein nanoshells help
define their structure and morphology, and determine their viruses' biological
function.


  Motivation: Complex computational pipelines are becoming a staple of modern
scientific research. Often these pipelines are resource intensive and require
days of computing time. In such cases, it makes sense to run them over
distributed computer clusters where they can take advantage of the aggregated
resources of many powerful computers. In addition to this, researchers often
want to integrate their workflows into their own web servers. In these cases,
software is needed to manage the submission of jobs from the web interface to
the cluster and then return the results once the job has finished executing.
  Results: We have developed the Job Management System (JMS), a workflow
management system and interface for the Torque resource manager. The JMS
provides users with a user-friendly interface for creating complex workflows
with multiple stages. It integrates this workflow functionality with Torque, a
tool that is used to control and manage batch jobs on distributed computing
clusters. The JMS can be used by researchers to build and run complex
computational pipelines and provides functionality to include these pipelines
in external interfaces. The JMS is currently being used to house a number of
structural bioinformatics pipelines at the Research Unit in Bioinformatics
(RUBi) at Rhodes University.
  Availability: The JMS is an open-source project and is freely available at
https://github.com/RUBi-ZA/JMS


  Information theory is gaining popularity as a tool to characterize
performance of biological systems. However, information is commonly quantified
without reference to whether or how a system could extract and use it; as a
result, information-theoretic quantities are easily misinterpreted. Here we
take the example of pattern-forming developmental systems which are commonly
structured as cascades of sequential gene expression steps. Such a multi-tiered
structure appears to constitute sub-optimal use of the positional information
provided by the input morphogen because noise is added at each tier. However,
the conventional theory fails to distinguish between the total information in a
morphogen and information that can be usefully extracted and interpreted by
downstream elements. We demonstrate that quantifying the information that is
_accessible_ to the system naturally explains the prevalence of multi-tiered
network architectures as a consequence of the noise inherent to the control of
gene expression. We support our argument with empirical observations from
patterning along the major body axis of the fruit fly embryo. Our results
exhibit the limitations of the standard information-theoretic characterization
of biological signaling and illustrate how they can be resolved.


  Genome rearrangements are evolutionary events that shuffle genomic
architectures. Most frequent genome rearrangements are reversals,
translocations, fusions, and fissions. While there are some more complex genome
rearrangements such as transpositions, they are rarely observed and believed to
constitute only a small fraction of genome rearrangements happening in the
course of evolution. The analysis of transpositions is further obfuscated by
intractability of the underlying computational problems.
  We propose a computational method for estimating the rate of transpositions
in evolutionary scenarios between genomes. We applied our method to a set of
mammalian genomes and estimated the transpositions rate in mammalian evolution
to be around 0.26.


  Many single molecule experiments for molecular motors comprise not only the
motor but also large probe particles coupled to it. The theoretical analysis of
these assays, however, often takes into account only the degrees of freedom
representing the motor. We present a coarse-graining method that maps a model
comprising two coupled degrees of freedom which represent motor and probe
particle to such an effective one-particle model by eliminating the dynamics of
the probe particle in a thermodynamically and dynamically consistent way. The
coarse-grained rates obey a local detailed balance condition and reproduce the
net currents. Moreover, the average entropy production as well as the
thermodynamic efficiency is invariant under this coarse-graining procedure. Our
analysis reveals that only by assuming unrealistically fast probe particles,
the coarse-grained transition rates coincide with the transition rates of the
traditionally used one-particle motor models. Additionally, we find that for
multicyclic motors the stall force can depend on the probe size. We apply this
coarse-graining method to specific case studies of the F1-ATPase and the
kinesin motor.


  Large scale models of physical phenomena demand the development of new
statistical and computational tools in order to be effective. Many such models
are `sloppy', i.e., exhibit behavior controlled by a relatively small number of
parameter combinations. We review an information theoretic framework for
analyzing sloppy models. This formalism is based on the Fisher Information
Matrix, which we interpret as a Riemannian metric on a parameterized space of
models. Distance in this space is a measure of how distinguishable two models
are based on their predictions. Sloppy model manifolds are bounded with a
hierarchy of widths and extrinsic curvatures. We show how the manifold boundary
approximation can extract the simple, hidden theory from complicated sloppy
models. We attribute the success of simple effective models in physics as
likewise emerging from complicated processes exhibiting a low effective
dimensionality. We discuss the ramifications and consequences of sloppy models
for biochemistry and science more generally. We suggest that the reason our
complex world is understandable is due to the same fundamental reason: simple
theories of macroscopic behavior are hidden inside complicated microscopic
processes.


  Plant survival is a key factor in forest dynamics and survival probabilities
often vary across life stages. Studies specifically aimed at assessing tree
survival are unusual and so data initially designed for other purposes often
need to be used; such data are more likely to contain errors than data
collected for this specific purpose. We investigate the survival rates of ten
tree species in a dataset designed to monitor growth rates. As some individuals
were not included in the census at some time points we use
capture-mark-recapture methods both to allow us to account for missing
individuals, and to estimate relocation probabilities. Growth rates, size, and
light availability were included as covariates in the model predicting survival
rates. The study demonstrates that tree mortality is best described as constant
between years and size-dependent at early life stages and size independent at
later life stages for most species of UK hardwood. We have demonstrated that
even with a twenty-year dataset it is possible to discern variability both
between individuals and between species. Our work illustrates the potential
utility of the method applied here for calculating plant population dynamics
parameters in time replicated datasets with small sample sizes and missing
individuals without any loss of sample size, and including explanatory
covariates.


  Given a mathematical model quantifying the viral infection of pandemic
influenza H1N1pdm09-H275 wild type (WT) and H1N1pdm09-H275Y mutant (MUT)
strains, we describe a simple method of estimating the model's constant
parameters using Monte Carlo methods. Monte Carlo parameter estimation methods
present certain advantages over the bootstrapping methods previously used in
such studies: the result comprises actual parameter distributions (posteriors)
that can be used to compare different viral strains; the recovered parameter
distributions offer an exact method to compute credible intervals (similar to
the frequentist 95% parametric confidence intervals (CI)), that, in turn, using
a suitable analysis statistic, will be narrower than the ones obtained from
bootstrapping; given an appropriate computational parallelization, Monte Carlo
methods are also faster and less computationally intensive than bootstrapping.
We fit Gaussian distributions to the parameter posterior distributions and use
a two-sided Kolmogorov-Smirnoff test to compare the two strains from a
parametric point of view; our example result shows that the two strains are 94%
different. Furthermore, based on the obtained parameter values, we estimate the
reproductive number R0 for each strain and show that the infectivity of the
mutant strain is larger than the wild type strain.


  Polymerization of fibrin, the primary structural protein of blood clots and
thrombi, occurs through binding of knobs 'A' and 'B' in the central nodule of
fibrin monomer to complementary holes 'a' and 'b' in the beta- and
gamma-nodules, respectively, of another monomer. We characterized the A:a and
B:b knob-hole interactions under varying solution conditions using Molecular
Dynamics simulations of the structural models of fibrin(ogen) fragment D
complexed with synthetic peptides GPRP (knob 'A' mimetic) and GHRP (knob 'B'
mimetic). The strength of A:a and B:b knob-hole complexes was roughly equal,
decreasing with pulling force; yet, the dissociation kinetics were sensitive to
variations in acidity (pH=5-7) and temperature (T=25-37 C). There were similar
structural changes in holes 'a' and 'b' during forced dissociation of the
knob-hole complexes: elongation of loop I, stretching of interior region, and
translocation of the moveable flap. The disruption of the knob-hole
interactions was not an "all-or-none" transition, as it occurred through
distinct two-step or single-step pathways with or without intermediate states.
The knob-hole bonds were stronger, tighter, and more brittle at pH=7 than at
pH=5. The B:b knob-hole bonds were weaker, looser, and more compliant than the
A:a knob-hole bonds at pH=7, but stronger, tighter, and less compliant at pH=5.
Surprisingly, the knob-hole bonds were stronger, not weaker, at elevated
temperature (T=37 C) compared to T=25 C due to the helix-to-coil transition in
loop I, which helps stabilize the bonds. These results provide detailed
qualitative and quantitative characteristics underlying the most significant
non-covalent interactions involved in fibrin polymerization.


  Many viruses rely on the self-assembly of their capsids to protect and
transport their genomic material. For many viral systems, in particular for
human viruses like hepatitis B, adeno or human immunodeficiency virus, that
lead to persistent infections, capsomeres are continuously produced in the
cytoplasm of the host cell while completed capsids exit the cell for a new
round of infection. Here we use coarse-grained Brownian dynamics simulations of
a generic patchy particle model to elucidate the role of the dynamic supply of
capsomeres for the reversible self-assembly of empty T1 icosahedral virus
capsids. We find that for high rates of capsomere influx only a narrow range of
bond strengths exists for which a steady state of continuous capsid production
is possible. For bond strengths smaller and larger than this optimal value, the
reaction volume becomes crowded by small and large intermediates, respectively.
For lower rates of capsomere influx a broader range of bond strengths exists
for which a steady state of continuous capsid production is established,
although now the production rate of capsids is smaller. Thus our simulations
suggest that the importance of an optimal bond strength for viral capsid
assembly typical for in vitro conditions can be reduced by the dynamic influx
of capsomeres in a cellular environment.


  In order to replicate within their cellular host, many viruses have developed
self-assembly strategies for their capsids which are sufficiently robust as to
be reconstituted in vitro. Mathematical models for virus self-assembly usually
assume that the bonds leading to cluster formation have constant reactivity
over the time course of assembly (direct assembly). In some cases, however,
binding sites between the capsomers have been reported to be activated during
the self-assembly process (hierarchical assembly). In order to study possible
advantages of such hierarchical schemes for icosahedral virus capsid assembly,
we use Brownian dynamics simulations of a patchy particle model that allows us
to switch binding sites on and off during assembly. For T1 viruses, we
implement a hierarchical assembly scheme where inter-capsomer bonds become
active only if a complete pentamer has been assembled. We find direct assembly
to be favorable for reversible bonds allowing for repeated structural
reorganizations, while hierarchical assembly is favorable for strong bonds with
small dissociation rate, as this situation is less prone to kinetic trapping.
However, at the same time it is more vulnerable to monomer starvation during
the final phase. Increasing the number of initial monomers does have only a
weak effect on these general features. The differences between the two assembly
schemes become more pronounced for more complex virus geometries, as shown here
for T3 viruses, which assemble through homogeneous pentamers and heterogeneous
hexamers in the hierarchical scheme. In order to complement the simulations for
this more complicated case, we introduce a master equation approach that agrees
well with the simulation results.


  Was the spread of agropastoralism from the Eurasian founder regions dominated
by demic or by cultural diffusion? This study employs a mathematical model of
regional sociocultural development that includes different diffusion processes,
local innovation and societal adaptation. Simulations hindcast the emergence
and expansion of agropastoral life style in 294 regions of Eurasia and North
Africa. Different scenarios for demic and diffusive exchange processes between
adjacent regions are contrasted and the spatiotemporal pattern of diffusive
events is evaluated. This study supports from a modeling perspective the
hypothesis that there is no simple or exclusive demic or cultural diffusion,
but that in most regions of Eurasia a combination of demic and cultural
processes were important. Furthermore, we demonstrate the strong spatial and
temporal variability in the balance of spread processes. Each region shows
sometimes more demic, and at other times more cultural diffusion. Only few,
possibly environmentally marginal, areas show a dominance of demic diffusion.
This study affirms that diffusion processes should be investigated in a
diachronic fashion and not from a time-integrated perspective.


  To recover the long-time behavior and the statistics of molecular
trajectories from the large number (tens of thousands) of their short
fragments, obtained by super-resolution methods at the single molecule level,
data analysis based on a stochastic model of their non-equilibrium motion is
required. Recently, we characterized the local biophysical properties
underlying receptor motion based on coarse-grained long-range interactions,
corresponding to attracting potential wells of large sizes. The purpose of this
letter is to discuss optimal estimators and show what happens when thing goes
wrong.


  The growing probabilities of additional offspring with the beneficial
reversal allele for various population sizes, $N$, sequence lengths, $L$,
selective advantages, $s$, fitness parameters, $k$, and measuring parameters,
$C$, were calculated in an asymmetric sharply-peaked landscape using the
decoupled continuous-time mutation-selection model. The growing probability in
the stochastic region was inversely proportional to the measuring parameter
when $C<1/Ns^*$, bent when $C\approx 1/Ns^*$ and became saturated when
$C>1/Ns^*$, where $s^*$ is the effective selective advantage. The saturated
growing probability in the stochastic region was approximately the effective
selective advantage when $C\gg 1/Ns^*$ and $s^*\ll 1$. The present study
suggests that the growing probability in the stochastic region in an asymmetric
sharply-peaked landscape in the decoupled continuous-time mutation-selection
model can be described using the theoretical formula for the growing
probability in the Moran two-allele model. The selective advantage ratio, which
represents the ratio of the effective selective advantage to the selective
advantage, does not depend on the population size, selective advantage,
measuring parameter, and fitness parameter; instead the selective advantage
ratio decreases with increasing the sequence length.


  Recent studies suggest that cooperative decision-making in one-shot
interactions is a history-dependent dynamic process: promoting intuition versus
deliberation has typically a positive effect on cooperation (dynamism) among
people living in a cooperative setting and with no previous experience in
economic games on cooperation (history-dependence). Here we report on a lab
experiment exploring how these findings transfer to a non-cooperative setting.
We find two major results: (i) promoting intuition versus deliberation has no
effect on cooperative behavior among inexperienced subjects living in a
non-cooperative setting; (ii) experienced subjects cooperate more than
inexperienced subjects, but only under time pressure. These results suggest
that cooperation is a learning process, rather than an instinctive impulse or a
self-controlled choice, and that experience operates primarily via the channel
of intuition. In doing so, our findings shed further light on the cognitive
basis of human cooperative decision-making and provide further support for the
recently proposed Social Heuristics Hypothesis.


  Cells generally change their internal state to adapt to an environmental
change, and accordingly evolve in response to the new conditions. This process
involves phenotypic changes that occur over several different time scales,
ranging from faster environmental adaptation without a corresponding change in
the genomic sequence to slower evolutionary dynamics involving genetic
mutations and subsequent selection. In this regard, a question arises as to
whether there are any relationships between such phenotypic changes over the
different time scales at which adaptive evolution occurs. In this study, we
analyzed simulated adaptive evolution in a simple cell model, and found that
proportionality between concentration changes in adaptation and evolution over
all components, and the proportion coefficients were closely linked to the
change in the growth rate of a cell. Furthermore, we demonstrated that the
phenotypic variances in component concentrations due to (non-genetic) noise and
genomic alternations are proportional across all components. These global
relationships in cellular states were also supported by phenomenological theory
and transcriptome analysis of laboratory evolution in {\it Escherichia coli}.
These findings provide a basis for the development of a quantitative theory of
plasticity and robustness, and to determine the general restriction of
phenotypic changes imposed by evolution.


  FoF1-ATP synthases are membrane-embedded protein machines that catalyze the
synthesis of adenosine triphosphate. Using photoactivation-based localization
microscopy (PALM) in TIR-illumination as well as structured illumination
microscopy (SIM), we explore the spatial distribution and track single FoF1-ATP
synthases in living E. coli cells under physiological conditions at different
temperatures. For quantitative diffusion analysis by mean-squared-displacement
measurements, the limited size of the observation area in the membrane with its
significant membrane curvature has to be considered. Therefore, we applied a
'sliding observation window' approach (M. Renz et al., Proc. SPIE 8225, 2012)
and obtained the one-dimensional diffusion coefficient of FoF1-ATP synthase
diffusing on the long axis in living E. coli cells.


  Defining the extent of epistasis - the non-independence of the effects of
mutations - is essential for understanding the relationship of genotype,
phenotype, and fitness in biological systems. The applications cover many areas
of biological research, including biochemistry, genomics, protein and systems
engineering, medicine, and evolutionary biology. However, the quantitative
definitions of epistasis vary among fields, and its analysis beyond just
pairwise effects remains obscure in general. Here, we show that different
definitions of epistasis are versions of a single mathematical formalism - the
weighted Walsh-Hadamard transform. We discuss that one of the definitions, the
backgound-averaged epistasis, is the most informative when the goal is to
uncover the general epistatic structure of a biological system, a description
that can be rather different from the local epistatic structure of specific
model systems. Key issues are the choice of effective ensembles for averaging
and to practically contend with the vast combinatorial complexity of mutations.
In this regard, we discuss possible approaches for optimally learning the
epistatic structure of biological systems.


  Lifetime distributions of social entities, such as enterprises, products, and
media contents, are one of the fundamental statistics characterizing the social
dynamics. To investigate the lifetime distribution of mutually interacting
systems, simple models having a rule for additions and deletions of entities
are investigated. We found a quite universal lifetime distribution for various
kinds of inter-entity interactions, and it is well fitted by a
stretched-exponential function with an exponent close to 1/2. We propose a
"modified Red-Queen" hypothesis to explain this distribution. We also review
empirical studies on the lifetime distribution of social entities, and
discussed the applicability of the model.


  The Gompertz law of mortality quantitatively describes the mortality rate of
humans and almost all multicellular animals. However, its underlying kinetic
mechanism is unclear. The Gompertz law cannot explain the effect of temperature
on lifespan and the mortality plateau at advanced ages. In this study a
reaction kinetics model with a time dependent rate coefficient is proposed to
describe the survival and senescence processes. A temperature-dependent
mortality function was derived. The new mortality function becomes the Gompertz
mortality function with the same relationship of parameters prescribed by the
Strehler-Mildvan correlation when age is smaller than a characteristic value
{\delta}, and reaches the mortality plateau when age is greater than {\delta}.
A closed-form analytical expression for describing the relationship of average
lifespan with temperature and other equations are derived from the new
mortality function. The derived equations can be used to estimate the limit of
average lifespan, predict the maximal longevity, calculate the temperature
coefficient of lifespan, and explain the tendency of survival curve. This
prediction is consistent with the most recently reported mortality trajectories
for single-year birth cohorts. This study suggests that the senescence process
results from the imbalance between damaging energy and protecting energy for
the critical chemical substance in the body. The rate of senescence of the
organism increases while the protecting energy decreases. The mortality plateau
is reached when the protecting energy decreases to its minimal levels. The
decreasing rate of the protecting energy is temperature-dependent. This study
is exploring the connection between biochemical mechanism and demography.


  The emergence of autocatalytic sets of molecules seems to have played an
important role in the origin of life context. Although the possibility to
reproduce this emergence in laboratory has received considerable attention,
this is still far from being achieved. In order to unravel some key properties
enabling the emergence of structures potentially able to sustain their own
existence and growth, in this work we investigate the probability to observe
them in ensembles of random catalytic reaction networks characterized by
different structural properties. From the point of view of network topology, an
autocatalytic set have been defined either in term of strongly connected
components (SCCs) or as reflexively autocatalytic and food-generated sets
(RAFs). We observe that the average level of catalysis differently affects the
probability to observe a SCC or a RAF, highlighting the existence of a region
where the former can be observed, whereas the latter cannot. This parameter
also affects the composition of the RAF, which can be further characterized
into linear structures, autocatalysis or SCCs. Interestingly, we show that the
different network topology (uniform as opposed to power-law catalysis systems)
does not have a significantly divergent impact on SCCs and RAFs appearance,
whereas the proportion between cleavages and condensations seems instead to
play a role. A major factor that limits the probability of RAF appearance and
that may explain some of the difficulties encountered in laboratory seems to be
the presence of molecules which can accumulate without being substrate or
catalyst of any reaction.


  This study examines features such as the population structure, growth and
reproduction of the pearl mullet captured in the Lake Ercek. Between the dates
January 2008 and January 2010 a total number of 527 individual were sampled
through using trammel net and mesh size of 20,22 and 24 mm of fish net. The
fork length, total and gonad weight, sex and fecundity of the pearl mullet
which were captured has been identified; their age determination has been done
through their operculum and otolith; their condition factor and gonadosomatic
index value has been estimated.


  This chapter presents SMuGA, an integration of symbiogenesis with the
Multiset Genetic Algorithm (MuGA). The symbiogenetic approach used here is
based on the host-parasite model with the novelty of varying the length of
parasites along the evolutionary process. Additionally, it models
collaborations between multiple parasites and a single host. To improve
efficiency, we introduced proxy evaluation of parasites, which saves fitness
function calls and exponentially reduces the symbiotic collaborations produced.
Another novel feature consists of breaking the evolutionary cycle into two
phases: a symbiotic phase and a phase of independent evolution of both hosts
and parasites. SMuGA was tested in optimization of a variety of deceptive
functions, with results one order of magnitude better than state of the art
symbiotic algorithms. This allowed to optimize deceptive problems with large
sizes, and showed a linear scaling in the number of iterations to attain the
optimum.


  Which nodes are most vulnerable to an epidemic spreading through a network,
and which carry the highest risk of causing a major outbreak if they are the
source of the infection? Here we show how these questions can be answered to
good approximation using the cavity method. Several curious properties of node
vulnerability and risk are explored: some nodes are more vulnerable than others
to weaker infections, yet less vulnerable to stronger ones; a node is always
more likely to be caught in an outbreak than it is to start one, except when
the disease has a deterministic lifetime; the rank order of node risk depends
on the details of the distribution of infectious periods.


  Tomasetti and Vogelstein (2015a) find that the incidence of a set of cancer
types is correlated with the total number of normal stem cell divisions. Here,
we separate the effects of standing stem cell number (i.e., organ or tissue
size) and per stem cell lifetime replication rate. We show that each has a
statistically significant and independent effect on explaining variation in
cancer incidence over the 31 cases considered by Tomasetti and Vogelstein. When
considering the total number of stem cell divisions and when removing cases
associated with disease or carcinogens, we find that cancer incidence attains a
plateau of approximately 0.6% incidence for the cases considered by these
authors. We further demonstrate that grouping by anatomical site explains most
of the remaining variation in risk between cancer types. This new analysis
suggests that cancer risk depends not only on the number of stem cell divisions
but varies enormously ($\sim$10,000 times) depending on the stem cell's
environment. Future research should investigate how tissue characteristics
(anatomical site, type, size, stem cell divisions) explain cancer incidence
over a wider range of cancers, to what extent different tissues express
specific protective mechanisms, and whether any differential protection can be
attributed to natural selection.


  Epigenetics has captured the attention of scientists in the past decades, yet
its scope has been continuously changing. In this paper, we give an overview on
how and why its definition has evolved and suggest several clarification on the
concepts used in this field, in particular, on the notions of epigenetic
information, epigenetic stability and epigenetic templating. Another issue that
we address is the role of epigenetic information. Not only it is important in
allowing alternative interpretations of genetic information, but it appears to
be important in protecting the genetic information, moreover, we suggest that
this function appeared first in evolution and only later on the epigenetic
mechanisms were recruited to play a role in cell differentiation.


  Quantifying the attack ratio of disease is key to epidemiological inference
and Public Health planning. For multi-serotype pathogens, however, different
levels of serotype-specific immunity make it difficult to assess the population
at risk. In this paper we propose a Bayesian method for estimation of the
attack ratio of an epidemic and the initial fraction of susceptibles using
aggregated incidence data. We derive the probability distribution of the
effective reproductive number, R t , and use MCMC to obtain posterior
distributions of the parameters of a single-strain SIR transmission model with
time-varying force of infection. Our method is showcased in a data set
consisting of 18 years of dengue incidence in the city of Rio de Janeiro,
Brazil. We demonstrate that it is possible to learn about the initial fraction
of susceptibles and the attack ratio even in the absence of serotype specific
data. On the other hand, the information provided by this approach is limited,
stressing the need for detailed serological surveys to characterise the
distribution of serotype-specific immunity in the population.


  The Ebola virus is spreading throughout West Africa and is causing thousands
of deaths. In order to quantify the effectiveness of different strategies for
controlling the spread, we develop a mathematical model in which the
propagation of the Ebola virus through Liberia is caused by travel between
counties. For the initial months in which the Ebola virus spreads, we find that
the arrival times of the disease into the counties predicted by our model are
compatible with World Health Organization data, but we also find that reducing
mobility is insufficient to contain the epidemic because it delays the arrival
of Ebola virus in each county by only a few weeks. We study the effect of a
strategy in which safe burials are increased and effective hospitalisation
instituted under two scenarios: (i) one implemented in mid-July 2014 and (ii)
one in mid-August---which was the actual time that strong interventions began
in Liberia. We find that if scenario (i) had been pursued the lifetime of the
epidemic would have been three months shorter and the total number of infected
individuals 80\% less than in scenario (ii). Our projection under scenario (ii)
is that the spreading will stop by mid-spring 2015.


  Motivation: Predictive modelling of gene expression is a powerful framework
for the in silico exploration of transcriptional regulatory interactions
through the integration of high-throughput -omics data. A major limitation of
previous approaches is their inability to handle conditional and synergistic
interactions that emerge when collectively analysing genes subject to different
regulatory mechanisms. This limitation reduces overall predictive power and
thus the reliability of downstream biological inference.
  Results: We introduce an analytical modelling framework (TREEOME: tree of
models of expression) that integrates epigenetic and transcriptomic data by
separating genes into putative regulatory classes. Current predictive modelling
approaches have found both DNA methylation and histone modification epigenetic
data to provide little or no improvement in accuracy of prediction of
transcript abundance despite, for example, distinct anti-correlation between
mRNA levels and promoter-localised DNA methylation. To improve on this, in
TREEOME we evaluate four possible methods of formulating gene-level DNA
methylation metrics, which provide a foundation for identifying gene-level
methylation events and subsequent differential analysis, whereas most previous
techniques operate at the level of individual CpG dinucleotides. We demonstrate
TREEOME by integrating gene-level DNA methylation (bisulfite-seq) and histone
modification (ChIP-seq) data to accurately predict genome-wide mRNA transcript
abundance (RNA-seq) for H1-hESC and GM12878 cell lines.
  Availability: TREEOME is implemented using open-source software and made
available as a pre-configured bootable reference environment. All scripts and
data presented in this study are available online at
http://sourceforge.net/projects/budden2015treeome/.


  Data tables in the form of spreadsheets or delimited text files are the most
utilised data format in Systems Biology. However, they are often not
sufficiently structured and lack clear naming conventions that would be
required for modelling. We propose the SBtab format as an attempt to establish
an easy-to-use table format that is both flexible and clearly structured. It
comprises defined table types for different kinds of data; syntax rules for
usage of names, shortnames, and database identifiers used for annotation; and
standardised formulae for reaction stoichiometries. Predefined table types can
be used to define biochemical network models and the biochemical constants
therein. The user can also define own table types, adjusting SBtab to other
types of data. Software code, tools, and further information can be found at
www.sbtab.net.


  In this paper we describe a method to identify "relevant subsets" of
variables, useful to understand the organization of a dynamical system. The
variables belonging to a relevant subset should have a strong integration with
the other variables of the same relevant subset, and a much weaker interaction
with the other system variables. On this basis, extending previous works on
neural networks, an information-theoretic measure is introduced, i.e. the
Dynamical Cluster Index, in order to identify good candidate relevant subsets.
The method does not require any previous knowledge of the relationships among
the system variables, but relies on observations of their values in time. We
show its usefulness in several application domains, including: (i) random
boolean networks, where the whole network is made of different subnetworks with
different topological relationships (independent or interacting subnetworks);
(ii) leader-follower dynamics, subject to noise and fluctuations; (iii)
catalytic reaction networks in a flow reactor; (iv) the MAPK signaling pathway
in eukaryotes. The validity of the method has been tested in cases where the
data are generated by a known dynamical model and the Dynamical Cluster Index
method is applied in order to uncover significant aspects of its organization;
however it is important to stress that it can also be applied to time series
coming from field data without any reference to a model. Given that it is based
on relative frequencies of sets of values, the method could be applied also to
cases where the data are not ordered in time. Several indications to improve
the scope and effectiveness of the Dynamical Cluster Index to analyze the
organization of complex systems are finally given.


  Clathrin-mediated endocytosis (CME) is a key pathway for transporting cargo
into cells via membrane vesicles. It plays an integral role in nutrient import,
signal transduction, neurotransmission and cellular entry of pathogens and
drug-carrying nanoparticles. As CME entails substantial local remodeling of the
plasma membrane, the presence of membrane tension offers resistance to bending
and hence, vesicle formation. Experiments show that in such high tension
conditions, actin dynamics is required to carry out CME successfully. In this
study, we build upon these pioneering experimental studies to provide
fundamental mechanistic insights into the roles of two key endocytic proteins,
namely, actin and BAR proteins in driving vesicle formation in high membrane
tension environment. Our study reveals a new actin force induced `snap-through
instability' that triggers a rapid shape transition from a shallow invagination
to a highly invaginated tubular structure. We show that the association of BAR
proteins stabilizes vesicles and induces a milder instability. In addition, we
present a new counterintuitive role of BAR depolymerization in regulating the
shape evolution of vesicles. We show that the dissociation of BAR proteins,
supported by actin-BAR synergy, leads to considerable elongation and squeezing
of vesicles. Going beyond the membrane geometry, we put forth a new
stress-based perspective for the onset of vesicle scission and predict the
shapes and composition of detached vesicles. We present the snap-through
transition and the high in-plane stress as possible explanations for the
intriguing direct transformation of broad and shallow invaginations into
detached vesicles in BAR mutant yeast cells.


  Concentration gradients of signaling molecules are essential for patterning
during development and they have been observed in both unicellular and
multicellular systems. In subcellular systems, clustering of the signaling
molecule has been observed. We develop a theoretical model of cluster- mediated
concentration gradient formation based on the Becker-Doring equations of
aggregation- fragmentation processes. We show that such a mechanism produces
robust concentration gradients on realistic time and spatial scales so long as
the process of clustering does not significantly stabilize the signaling
molecule. Finally, we demonstrate that such a model is applicable to the pom1p
subcellular gradient in fission yeast.


  We propose a lightweight data structure for indexing and querying collections
of NGS reads data in main memory. The data structure supports the interface
proposed in the pioneering work by Philippe et al. for counting and locating
$k$-mers in sequencing reads. Our solution, PgSA (pseudogenome suffix array),
based on finding overlapping reads, is competitive to the existing algorithms
in the space use, query times, or both. The main applications of our index
include variant calling, error correction and analysis of reads from RNA-seq
experiments.


  Self-assembly of polypeptides into fibrillar structures can be initiated by
planar surfaces that interact favorably with certain residues. Using a coarse
grained model, we systematically studied the folding and adsorption behavior of
a $\beta$-roll forming polypeptide. We find that there are two different
folding pathways depending on the temperature: (i) at low temperature, the
polypeptide folds in solution into a $\beta$-roll before adsorbing onto the
attractive surface, (ii) at higher temperature, the polypeptide first adsorbs
in a disordered state, and folds while on the surface. The folding temperature
increases with increasing attraction, as the folded $\beta$-roll is stabilized
by the surface. Surprisingly, further increasing the attraction lowers the
folding temperature again, as strong attraction also stabilizes the adsorbed
disordered state, which competes with folding of the polypeptide. Our results
suggest that to enhance the folding, one should use a weakly attractive
surface. They also explain the recent experimental observation of the
non-monotonic effect of charge on the fibril formation on an oppositely charged
surface [C. Charbonneau, et al., ACS Nano, 8, 2328 (2014)].


  The last decade has seen an explosion in models that describe phenomena in
systems medicine. Such models are especially useful for studying signaling
pathways, such as the Wnt pathway. In this chapter we use the Wnt pathway to
showcase current mathematical and statistical techniques that enable modelers
to gain insight into (models of) gene regulation, and generate testable
predictions. We introduce a range of modeling frameworks, but focus on ordinary
differential equation (ODE) models since they remain the most widely used
approach in systems biology and medicine and continue to offer great potential.
We present methods for the analysis of a single model, comprising applications
of standard dynamical systems approaches such as nondimensionalization, steady
state, asymptotic and sensitivity analysis, and more recent statistical and
algebraic approaches to compare models with data. We present parameter
estimation and model comparison techniques, focusing on Bayesian analysis and
coplanarity via algebraic geometry. Our intention is that this (non exhaustive)
review may serve as a useful starting point for the analysis of models in
systems medicine.


  Humans have $23$ pairs of homologous chromosomes. The homologous pairs are
almost identical pairs of chromosomes. For the most part, differences in
homologous chromosome occur at certain documented positions called single
nucleotide polymorphisms (SNPs). A haplotype of an individual is the pair of
sequences of SNPs on the two homologous chromosomes. In this paper, we study
the problem of inferring haplotypes of individuals from mate-pair reads of
their genome. We give a simple formula for the coverage needed for haplotype
assembly, under a generative model. The analysis here leverages connections of
this problem with decoding convolutional codes.


  Brown algae are sessile macro-organisms of great ecological relevance in
coastal ecosystems. They evolved independently from land plants and other
multicellular lineages, and therefore hold several original ontogenic and
metabolic features. Most brown algae grow along the coastal zone where they
face frequent environmental changes, including exposure to toxic levels of
heavy metals such as copper (Cu). We carried out large-scale transcriptomic and
metabolomic analyses to decipher the short-term acclimation of the brown algal
model E. siliculosus to Cu stress, and compared these data to results known for
other abiotic stressors. This comparison demonstrates that Cu induces oxidative
stress in E. siliculosus as illustrated by the transcriptomic overlap between
Cu and H2O2 treatments. The common response to Cu and H2O2 consisted in the
activation of the oxylipin and the repression of inositol signaling pathways,
together with the regulation of genes coding for several
transcription-associated proteins. Concomitantly, Cu stress specifically
activated a set of genes coding for orthologs of ABC transporters, a P1B-type
ATPase, ROS detoxification systems such as a vanadium-dependent
bromoperoxidase, and induced an increase of free fatty acid contents. Finally
we observed, as a common abiotic stress mechanism, the activation of autophagic
processes on one hand and the repression of genes involved in nitrogen
assimilation on the other hand. Comparisons with data from green plants
indicate that some processes involved in Cu and oxidative stress response are
conserved across these two distant lineages. At the same time the high number
of yet uncharacterized brown alga-specific genes induced in response to copper
stress underlines the potential to discover new components and molecular
interactions unique to these organisms. Of particular interest for future
research is the potential cross-talk between reactive oxygen species (ROS)-,
myo-inositol-, and oxylipin signaling.


  Experimental evidences show that in gene transcription, RNA polymerase has
the possibility to be stalled at certain position of the transcription
template. This may be due to the template damage, or protein barriers. Once
stalled, polymerase may backtrack along the template to the previous nucleotide
to wait for the repair of the damaged site, or simply bypass the barrier or
damaged site and consequently synthesize an incorrect messenger RNA, or degrade
and detach from the template. Thus, the {\it effective} transcription rate (the
rate to synthesize correct product mRNA) and the transcription {\it
effectiveness} (the ratio of the {\it effective} transcription rate to the {\it
effective} transcription initiation rate) are both influenced by polymerase
stalling events. This study shows that, Without backtracking, detachment of
stalled polymerase can also help to increase the {\it effective} transcription
rate and transcription {\it effectiveness}. Generally, the increase of bypass
rate of the stalled polymerase will lead to the decrease of the {\it effective}
transcription rate and transcription {\it effectiveness}. But when both
detachment rate and backtracking rate of the stalled polymerase vanish, the
{\it effective} transcription rate may also be increased by the bypass
mechanism.


  Translation is an important process for prokaryotic and eukaryotic cells to
produce necessary proteins for cell growth. Numerious experiments have been
performed to explore the translational properties. Diverse models have also
been developed to determine the biochemical mechanism of translation. However,
to simplify the majority of the existing models, the frameshifting of ribosome
along the mRNA transcript is neglected, which actually occurs in real cells and
has been extensively experimentally studied. The frameshifting of ribosome
evidently influences the efficiency and speed of translation, considering that
the peptide chains synthesized by shifted ribosomes will not fold into
functional proteins and will degrade rapidly. In this study, a theoretical
model is presented to describe the translational process based on the model for
totally asymmetric simple exclusion process. In this model, the frameshifting
of the ribosome along the mRNA transcript and the attachment/detachment of the
ribosome to/from the main body of mRNA codons during translation elongation
process, are explicitly included. The results show that, with ribosome
frameshifing, the speed of correctly synthesized peptide chains may increase
first and then decrease with both the translation initiation rate $\alpha$ and
the ribosome detachment rate $\omega_d$. This results indicates that regulating
the translation process to reach maximal synthesized speed of proteins is
theoretically feasible. Traffic-related problems of ribosome motion along the
mRNA transcript are also addressed theoretically. Depending on parameter
values, shock wave (or domain wall) may exist for ribosome probabilities along
the mRNA.


  30% of the DNA in E. coli bacteria is covered by proteins. Such high degree
of crowding affect the dynamics of generic biological processes (e.g. gene
regulation, DNA repair, protein diffusion etc.) in ways that are not yet fully
understood. In this paper, we theoretically address the diffusion constant of a
tracer particle in a one dimensional system surrounded by impenetrable crowder
particles. While the tracer particle always stays on the lattice, crowder
particles may unbind to a surrounding bulk and rebind at another or the same
location. In this scenario we determine how the long time diffusion constant
${\cal D}$ (after many unbinding events) depends on (i) the unbinding rate of
crowder particles $k_{\rm off}$, and (ii) crowder particle line density $\rho$,
from simulations (Gillespie algorithm) and analytical calculations. For small
$k_{\rm off}$, we find ${\cal D}\sim k_{\rm off}/\rho^2$ when crowder particles
are immobile on the line, and ${\cal D}\sim \sqrt{D k_{\rm off}}/\rho$ when
they are diffusing; $D$ is the free particle diffusion constant. For large
$k_{\rm off}$, we find agreement with mean-field results which do not depend on
$k_{\rm off}$. From literature values of $k_{\rm off}$ and $D$, we show that
the small $k_{\rm off}$-limit is relevant for in vivo protein diffusion on a
crowded DNA. Our results applies to single-molecule tracking experiments.


  A wide variety of large-scale data has been produced in bioinformatics. In
response, the need for efficient handling of biomedical big data has been
partly met by parallel computing. However, the time demand of many
bioinformatics programs still remains high for large-scale practical uses due
to factors that hinder acceleration by parallelization. Recently, new
generations of storage devices have emerged, such as NAND flash-based
solid-state drives (SSDs), and with the renewed interest in near-data
processing, they are increasingly becoming acceleration methods that can
accompany parallel processing. In certain cases, a simple drop-in replacement
of hard disk drives (HDDs) by SSDs results in dramatic speedup. Despite the
various advantages and continuous cost reduction of SSDs, there has been little
review of SSD-based profiling and performance exploration of important but
time-consuming bioinformatics programs. For an informative review, we perform
in-depth profiling and analysis of 23 key bioinformatics programs using
multiple types of devices. Based on the insight we obtain from this research,
we further discuss issues related to design and optimize bioinformatics
algorithms and pipelines to fully exploit SSDs. The programs we profile cover
traditional and emerging areas of importance, such as alignment, assembly,
mapping, expression analysis, variant calling, and metagenomics. We explain how
acceleration by parallelization can be combined with SSDs for improved
performance and also how using SSDs can expedite important bioinformatics
pipelines, such as variant calling by the Genome Analysis Toolkit (GATK) and
transcriptome analysis using RNA sequencing (RNA-seq). We hope that this review
can provide useful directions and tips to accompany future bioinformatics
algorithm design procedures that properly consider new generations of powerful
storage devices.


  We present an algorithm to estimate the configurational entropy $S$ of a
polymer. The algorithm uses the statistics of coincidences among random samples
of configurations and is related to the catch-tag-release method for estimation
of population sizes, and to the classic "birthday paradox". Bias in the entropy
estimation is decreased by grouping configurations in nearly equiprobable
partitions based on their energies, and estimating entropies separately within
each partition. Whereas most entropy estimation algorithms require $N\sim
2^{S}$ samples to achieve small bias, our approach typically needs only $N\sim
\sqrt{2^{S}}$. Thus the algorithm can be applied to estimate protein free
energies with increased accuracy and decreased computational cost.


  For a paradigmatic model of chemotaxis, we analyze the effect how a nonzero
affinity driving receptors out of equilibrium affects sensitivity. This
affinity arises whenever changes in receptor activity involve ATP hydrolysis.
The sensitivity integrated over a ligand concentration range is shown to be
enhanced by the affinity, providing a measure of how much energy consumption
improves sensing. With this integrated sensitivity we can establish an
intriguing analogy between sensing with nonequilibrium receptors and kinetic
proofreading: the increase in integrated sensitivity is equivalent to the
decrease of the error in kinetic proofreading. The influence of the occupancy
of the receptor on the phosphorylation and dephosphorylation reaction rates is
shown to be crucial for the relation between integrated sensitivity and
affinity. This influence can even lead to a regime where a nonzero affinity
decreases the integrated sensitivity, which corresponds to anti-proofreading.


  Epistasis is a key concept in the theory of adaptation. Indicators of
epistasis are of interest for large system where systematic fitness
measurements may not be possible. Some recent approaches depend on information
theory. We show that considering shared entropy for pairs of loci can be
misleading. The reason is that shared entropy does not imply epistasis for the
pair. This observation holds true also in the absence of higher order
epistasis. We discuss a refined approach for identifying pairwise interactions
using entropy.


  Ecological patterns arise from the interplay of many different processes, and
yet the emergence of consistent phenomena across a diverse range of ecological
systems suggests that many patterns may in part be determined by statistical or
numerical constraints. Differentiating the extent to which patterns in a given
system are determined statistically, and where it requires explicit ecological
processes, has been difficult. We tackled this challenge by directly comparing
models from a constraint-based theory, the Maximum Entropy Theory of Ecology
(METE) and models from a process-based theory, the size-structured neutral
theory (SSNT). Models from both theories were capable of characterizing the
distribution of individuals among species and the distribution of body size
among individuals across 76 forest communities. However, the SSNT models
consistently yielded higher overall likelihood, as well as more realistic
characterizations of the relationship between species abundance and average
body size of conspecific individuals. This suggests that the details of the
biological processes contain additional information for understanding community
structure that are not fully captured by the METE constraints in these systems.
Our approach provides a first step towards differentiating between process- and
constraint-based models of ecological systems and a general methodology for
comparing ecological models that make predictions for multiple patterns.


  We generated genome-wide data from 69 Europeans who lived between 8,000-3,000
years ago by enriching ancient DNA libraries for a target set of almost four
hundred thousand polymorphisms. Enrichment of these positions decreases the
sequencing required for genome-wide ancient DNA analysis by a median of around
250-fold, allowing us to study an order of magnitude more individuals than
previous studies and to obtain new insights about the past. We show that the
populations of western and far eastern Europe followed opposite trajectories
between 8,000-5,000 years ago. At the beginning of the Neolithic period in
Europe, ~8,000-7,000 years ago, closely related groups of early farmers
appeared in Germany, Hungary, and Spain, different from indigenous
hunter-gatherers, whereas Russia was inhabited by a distinctive population of
hunter-gatherers with high affinity to a ~24,000 year old Siberian6 . By
~6,000-5,000 years ago, a resurgence of hunter-gatherer ancestry had occurred
throughout much of Europe, but in Russia, the Yamnaya steppe herders of this
time were descended not only from the preceding eastern European
hunter-gatherers, but from a population of Near Eastern ancestry. Western and
Eastern Europe came into contact ~4,500 years ago, as the Late Neolithic Corded
Ware people from Germany traced ~3/4 of their ancestry to the Yamnaya,
documenting a massive migration into the heartland of Europe from its eastern
periphery. This steppe ancestry persisted in all sampled central Europeans
until at least ~3,000 years ago, and is ubiquitous in present-day Europeans.
These results provide support for the theory of a steppe origin of at least
some of the Indo-European languages of Europe.


  We present a computational modeling framework for data-driven simulations and
analysis of infectious disease spread in large populations. For the purpose of
efficient simulations, we devise a parallel solution algorithm targeting
multi-socket shared memory architectures. The model integrates infectious
dynamics as continuous-time Markov chains and available data such as animal
movements or aging are incorporated as externally defined events. To bring out
parallelism and accelerate the computations, we decompose the spatial domain
and optimize cross-boundary communication using dependency-aware task
scheduling. Using registered livestock data at a high spatio-temporal
resolution, we demonstrate that our approach not only is resilient to varying
model configurations, but also scales on all physical cores at realistic work
loads. Finally, we show that these very features enable the solution of inverse
problems on national scales.


  Transcriptional repression may cause transcriptional noise by a competition
between repressor and RNA polymerase binding. Although promoter activity is
often governed by a single limiting step, we argue here that the size of the
noise strongly depends on whether this step is the initial equilibrium binding
or one of the subsequent unidirectional steps. Overall, we show that
nonequilibrium steps of transcription initiation systematically increase the
cell-to-cell heterogeneity in bacterial populations. In particular, this allows
also weak promoters to give substantial transcriptional noise.


  We quantify the VDJ recombination and somatic hypermutation processes in
human B-cells using probabilistic inference methods on high-throughput DNA
sequence repertoires of human B-cell receptor heavy chains. Our analysis
captures the statistical properties of the naive repertoire, first after its
initial generation via VDJ recombination and then after selection for
functionality. We also infer statistical properties of the somatic
hypermutation machinery (exclusive of subsequent effects of selection). Our
main results are the following: the B-cell repertoire is substantially more
diverse than T-cell repertoires, due to longer junctional insertions; sequences
that pass initial selection are distinguished by having a higher probability of
being generated in a VDJ recombination event; somatic hypermutations have a
non-uniform distribution along the V gene that is well explained by an
independent site model for the sequence context around the hypermutation site.


  1. Accurate estimates of demographic parameters are required to infer
appropriate ecological relationships and inform management actions. Recently
developed N-mixture models use count data from unmarked individuals to estimate
demographic parameters, but a joint approach combining the strengths of both
analytical tools has not been developed. 2. We present an integrated model
combining known-fate and open N-mixture models, allowing the estimation of
detection probability, recruitment, and the joint estimation of survival. We
first use a simulation study to evaluate the performance of the model relative
to known values. We then provide an applied example using 4 years of wolf
survival data consisting of relocations of radio-collared wolves within packs
and counts of associated pack-mates. The model is implemented in both
maximum-likelihood and Bayesian frameworks using a new R package kfdnm and the
BUGS language. 3. The simulation results indicated that the integrated model
was able to reliably recover parameters with no evidence of bias, and estimates
were more precise under the joint model as expected. Results from the applied
example indicated that the marked sample of wolves was biased towards
individuals with higher apparent survival rates (including losses due to
mortality and emigration) than the unmarked pack-mates, suggesting estimates of
apparent survival based on joint estimation could be more representative of the
overall population. Estimates of recruitment were similar to direct
observations of pup production, and overlap of the credible intervals suggested
no clear differences in recruitment rates. 4. Our integrated model is a
practical approach for increasing the amount of information gained from future
and existing radio-telemetry and other similar mark-resight datasets.


  Steady state analysis of dynamical systems for biological networks give rise
to algebraic varieties in high-dimensional spaces whose study is of interest in
their own right. We demonstrate this for the shuttle model of the Wnt signaling
pathway. Here the variety is described by a polynomial system in 19 unknowns
and 36 parameters. Current methods from computational algebraic geometry and
combinatorics are applied to analyze this model.


  Random Boolean networks (RBNs) are frequently employed for modelling complex
systems driven by information processing, e.g. for gene regulatory networks
(GRNs). Here we propose a hierarchical adaptive RBN (HARBN) as a system
consisting of distinct adaptive RBNs - subnetworks - connected by a set of
permanent interlinks. Information measures and internal subnetworks topology of
HARBN coevolve and reach steady-states that are specific for a given network
structure. We investigate mean node information, mean edge information as well
as a mean node degree as functions of model parameters and demonstrate HARBN's
ability to describe complex hierarchical systems.


  A fundamental problem in the study of phylogenetic networks is to determine
whether or not a given phylogenetic network contains a given phylogenetic tree.
We develop a quadratic-time algorithm for this problem for binary nearly-stable
phylogenetic networks. We also show that the number of reticulations in a
reticulation visible or nearly stable phylogenetic network is bounded from
above by a function linear in the number of taxa.


  Recently, a growing interest has been addressed to the electrical properties
of bacteriorhodopsin (bR), a protein belonging to the transmembrane protein
family. To take into account the structure-dependent nature of the current, in
a previous set of papers we suggested a mechanism of sequential tunneling among
neighbouring amino acids. As a matter of fact, it is well accepted that, when
irradiated with green light, bR undergoes a conformational change at a
molecular level. Thus, the role played by the protein tertiary-structure in
modeling the charge transfer cannot be neglected. The aim of this paper is to
go beyond previous models, in the framework of a new branch of electronics, we
called proteotronics, which exploits the ability to use proteins as reliable,
well understood materials, for the development of novel bioelectronic devices.
In particular, the present approach assumes that the conformational change is
not the unique transformation that the protein undergoes when irradiated by
light. Instead, the light can also promote a free-energy increase of the
protein state that, in turn, should modify its internal degree of connectivity,
here described by the change in the value of an interaction radius associated
with the physical interactions among amino acids. The implemented model enables
us to achieve a better agreement between theory and experiments in the region
of a low applied bias by preserving the level of agreement at high values of
applied bias. Furthermore, results provide new insights on the mechanisms
responsible for bR photoresponse.


  An approach for approximately calculating the number of genes in a genome is
presented, which takes into account the average protein length expected for the
species. A number of virus, bacterial and eukaryotic genomes are scrutinized.
Genome figures are presented, which support the average protein size of a
species as a criterion for assessing life complexity. The human gene
distribution in the 23 chromosomes is investigated emphasizing the genomic
rate, the mean 'exon' length, and the mean 'exons per gene'. It is shown that
storing all genes of a single human definitely requires less than 12 MB.


  Summary: We present a new tool to correct sequencing errors in Illumina data
produced from high-coverage whole-genome shotgun resequencing. It uses a
non-greedy algorithm and shows comparable performance and higher accuracy in an
evaluation on real human data. This evaluation has the most complete collection
of high-performance error correctors so far.
  Availability and implementation: https://github.com/lh3/bfc
  Contact: hengli@broadinstitute.org


  Cooperation or defection and participation or withdrawal are well-known
options of behavior in game-like activities in free societies, yet the
co-evolutionary dynamics of these behavioral traits in the individual level are
not well understood. Here we investigate the continuous voluntary public good
game, in which individuals have two types of continuous-valued options: a
probability of joining the public good game and a level of cooperative
investment in the game. Our numerical results reveal hitherto unreported
phenomena: (i) The evolutionary dynamics are initially characterized by
oscillations in individual cooperation and participation levels, in contrast to
the population-level oscillations that have previously been reported. (ii)
Eventually, the population's average cooperation and participation levels
converge to and stabilize at a center. (iii) Then, a most peculiar phenomenon
unfolds: The strategies present in the population diversify and give rise to a
"cloud" of tinkering individuals who each tries out a different strategy, and
this process continues unchecked as long as the population's cooperation and
participation levels remain balanced. Over time, however, imbalances build up
as a consequence of random drift and there is a sudden and abrupt collapse of
the strategy-diversity cloud. The process then repeats again in a cyclic
manner. To understand the three aforementioned phenomena, we investigate the
system analytically using adaptive-dynamics techniques. Our analysis casts
light on the mechanisms which underpin the unexpected and surprising
evolutionary dynamics.


  We explore the hypothesis that the relative abundance of feedback loops in
many empirical complex networks is severely reduced owing to the presence of an
inherent global directionality. Aimed at quantifying this idea, we propose a
simple probabilistic model in which a free parameter $\gamma$ controls the
degree of inherent directionality. Upon strengthening such directionality, the
model predicts a drastic reduction in the fraction of loops which are also
feedback loops. To test this prediction, we extensively enumerated loops and
feedback loops in many empirical biological, ecological and socio-
technological directed networks. We show that, in almost all cases, empirical
networks have a much smaller fraction of feedback loops than network
randomizations. Quite remarkably, this empirical finding is quantitatively
reproduced, for all loop lengths, by our model by fitting its only parameter
$\gamma$. Moreover, the fitted value of $\gamma$ correlates quite well with
another direct measurement of network directionality, performed by means of a
novel algorithm. We conclude that the existence of an inherent network
directionality provides a parsimonious quantitative explanation for the
observed lack of feedback loops in empirical networks.


  We give a criterion for unlimited growth with positive probability for a
large class of multidimensional stochastic models. As a by-product, we recover
the necessary and sufficient conditions for recurrence and transience for
critical multitype Galton-Watson with immigration processes. We also
significantly improve some results on multitype size-dependent Galton-Watson
processes.


  To monitor conformational changes of individual membrane transporters in
liposomes in real time, we attach two fluorophores to selected domains of a
protein. Sequential distance changes between the dyes are recorded and analyzed
by Forster resonance energy transfer (FRET). Using freely diffusing membrane
proteins reconstituted in liposomes, observation times are limited by Brownian
motion through the confocal detection volume. A. E. Cohen and W. E. Moerner
have invented and built microfluidic devices to actively counteract Brownian
motion of single nanoparticles in electrokinetic traps (ABELtrap). Here we
present a version of an ABELtrap with a laser focus pattern generated by
electro-optical beam deflectors and controlled by a programmable FPGA. This
ABELtrap could hold single fluorescent nanobeads for more than 100 seconds,
increasing the observation times of a single particle by more than a factor of
1000. Conformational changes of single FRET-labeled membrane enzymes FoF1-ATP
synthase can be detected in the ABELtrap.


  We extend the theory of quasipotentials in dynamical systems by calculating,
within a broad class of period-doubling maps, an exact potential for the
critical fluctuations of pitchfork bifurcations in the weak noise limit. These
far-from-equilibrium fluctuations are described by finite-size mean field
theory, placing their static properties in the same universality class as the
Ising model on a complete graph. We demonstrate that the effective system size
of noisy period-doubling bifurcations exhibits universal scaling behavior along
period-doubling routes to chaos.


  We present MLEHaplo, a maximum likelihood de novo assembly algorithm for
reconstructing viral haplotypes in a virus population from paired-end next
generation sequencing (NGS) data. Using the pairing information of reads in our
proposed Viral Path Reconstruction Algorithm (ViPRA), we generate a small
subset of paths from a De Bruijn graph of reads that serve as candidate paths
for true viral haplotypes. Our proposed method MLEHaplo then generates a
maximum likelihood estimate of the viral population using the paths
reconstructed by ViPRA. We evaluate and compare MLEHaplo on simulated datasets
of 1200 base pairs at different sequence coverage, on HCV strains with
sequencing errors, and on a lab mixture of five HIV-1 strains. MLEHaplo
reconstructs full length viral haplotypes having a 100% sequence identity to
the true viral haplotypes in most of the small genome simulated viral
populations at 250x sequencing coverage. While reference based methods either
under-estimate or over-estimate the viral haplotypes, MLEHaplo limits the
over-estimation to 3 times the size of true viral haplotypes, reconstructs the
full phylogeny in the HCV to greater than 99% sequencing identity and captures
more sequencing variation for the HIV-1 strains dataset compared to their known
consensus sequences.


  A population protocol *stably elects a leader* if, for all $n$, starting from
an initial configuration with $n$ agents each in an identical state, with
probability 1 it reaches a configuration $\mathbf{y}$ that is correct (exactly
one agent is in a special leader state $\ell$) and stable (every configuration
reachable from $\mathbf{y}$ also has a single agent in state $\ell$). We show
that any population protocol that stably elects a leader requires $\Omega(n)$
expected "parallel time" --- $\Omega(n^2)$ expected total pairwise interactions
--- to reach such a stable configuration. Our result also informs the
understanding of the time complexity of chemical self-organization by showing
an essential difficulty in generating exact quantities of molecular species
quickly.


  Although accumulation of molecular damage is suggested to be an important
molecular mechanism of aging, a quantitative link between the dynamics of
damage accumulation and mortality of species has so far remained elusive. To
address this question, we examine stability properties of a generic gene
regulatory network (GRN) and demonstrate that many characteristics of aging and
the associated population mortality rate emerge as inherent properties of the
critical dynamics of gene regulation and metabolic levels. Based on the
analysis of age-dependent changes in gene-expression and metabolic profiles in
Drosophila melanogaster, we explicitly show that the underlying GRNs are nearly
critical and inherently unstable. This instability manifests itself as aging in
the form of distortion of gene expression and metabolic profiles with age, and
causes the characteristic increase in mortality rate with age as described by a
form of the Gompertz law. In addition, we explain late-life mortality
deceleration observed at very late ages for large populations. We show that
aging contains a stochastic component, related to accumulation of regulatory
errors in transcription/translation/metabolic pathways due to imperfection of
signaling cascades in the network and of responses to environmental factors. We
also establish that there is a strong deterministic component, suggesting
genetic control. Since mortality in humans, where it is characterized best, is
strongly associated with the incidence of age-related diseases, our findings
support the idea that aging is the driving force behind the development of
chronic human diseases.


  Classification is one of the most popular and widely used supervised learning
tasks, which categorizes objects into predefined classes based on known
knowledge. Classification has been an important research topic in machine
learning and data mining. Different classification methods have been proposed
and applied to deal with various real-world problems. Unlike unsupervised
learning such as clustering, a classifier is typically trained with labeled
data before being used to make prediction, and usually achieves higher accuracy
than unsupervised one.
  In this paper, we first define classification and then review several
representative methods. After that, we study in details the application of
classification to a critical problem in drug discovery, i.e., drug-target
prediction, due to the challenges in predicting possible interactions between
drugs and targets.


  The importance of studying properties of networks is manifest in diverse
fields ranging from biology, engineering, physics, chemistry, neuroscience, and
medicine. The functionality of networks with regard to performance, throughput,
reliability and robustness is strongly linked to the underlying geometric and
topological properties of the network and this is the focus of this paper,
especially as applied to certain biological networks. The fundamental
mathematical abstraction of a network as a weighted graph brings to bear the
tools of graph theory--a highly developed subject of mathematical research. But
more importantly, recently proposed geometric notions of curvature on very
general metric measure spaces allow us to utilize a whole new set of tools and
ideas that help quantify functionality and robustness of graphs. In particular,
robustness is closely connected to network entropy which, in turn, is very
closely related to curvature. We will see that there are a number of
alternative notions of discrete curvature that are compatible with the
classical Riemannian definition, each having its own advantages and
disadvantages, and are relevant to networks of interest. We will concentrate on
the role of curvature for certain key cancer networks in order to
quantitatively indicate their apparent functional robustness relative to their
normal counterparts.


  Mitochondrial research is important to ageing, apoptosis, and mitochondrial
diseases. In previous works, mitochondria are usually stimulated indirectly by
proapoptotic drugs to study mitochondrial development, which is in lack of
controllability, or spatial and temporal resolution. These chemicals or even
gene techniques regulating mitochondrial dynamics may also activate other
inter- or intra-cellular processes simultaneously. Here we demonstrate a
photostimulation method on single-mitochondrion level by tightly-focused
femtosecond laser that can precisely activate restorable fragmentation of
mitochondria which soon recover their original tubular structure after tens of
seconds. In this process, series of mitochondrial reactive oxygen species
(mROS) flashes are observed and found very critical to mitochondrial
fragmentation. Meanwhile, transient openings of mitochondrial permeability
transition pores (mPTP), suggested by oscillations of mitochondrial membrane
potential, contribute to the scavenging of redundant mROS and recovery of
fragmented mitochondria. Those results demonstrate photostimulation as an
active, precise and controllable method for the study of mitochondrial
oxidative and morphological dynamics or related fields.


  Infectious diseases spread through human networks.
Susceptible-Infected-Removed (SIR) model is one of the epidemic models to
describe infection dynamics on a complex network connecting individuals. In the
metapopulation SIR model, each node represents a population (group) which has
many individuals. In this paper, we propose a modified metapopulation SIR model
in which a latent period is taken into account. We call it SIIR model. We
divide the infection period into two stages: an infected stage, which is the
same as the previous model, and a seriously ill stage, in which individuals are
infected and cannot move to the other populations. The two infectious stages in
our modified metapopulation SIR model produce a discontinuous final size
distribution. Individuals in the infected stage spread the disease like
individuals in the seriously ill stage and never recover directly, which makes
an effective recovery rate smaller than the given recovery rate.


  Binomial ideals are special polynomial ideals with many algorithmically and
theoretically nice properties. We discuss the problem of deciding if a given
polynomial ideal is binomial. While the methods are general, our main
motivation and source of examples is the simplification of steady state
equations of chemical reaction networks. For homogeneous ideals we give an
efficient, Gr\"obner-free algorithm for binomiality detection, based on linear
algebra only. On inhomogeneous input the algorithm can only give a sufficient
condition for binomiality. As a remedy we construct a heuristic toolbox that
can lead to simplifications even if the given ideal is not binomial.


  During the lifecycle of many single-stranded RNA viruses, including many
human pathogens, a protein shell called the capsid spontaneously assembles
around the viral genome. Understanding the mechanisms by which capsid proteins
selectively assemble around the viral RNA amidst diverse host RNAs is a key
question in virology. In one proposed mechanism, sequence elements (packaging
sites) within the genomic RNA promote rapid and efficient assembly through
specific interactions with the capsid proteins. In this work we develop a
coarse-grained particle-based computational model for capsid proteins and RNA
which represents protein-RNA interactions arising both from non-specific
electrostatics and specific packaging sites interactions. Using Brownian
dynamics simulations, we explore how the efficiency and specificity of assembly
depend on solution conditions (which control protein-protein and nonspecific
protein-RNA interactions) as well as the strength and number of packaging
sites. We identify distinct regions in parameter space in which packaging sites
lead to highly specific assembly via different mechanisms, and others in which
packaging sites lead to kinetic traps. We relate these computational
predictions to in vitro assays for specificity in which cognate viral RNAs are
compete against non-cognate RNAs for assembly by capsid proteins.


  Background: Identifying all possible mapping locations of next-generation
sequencing (NGS) reads is highly essential in several applications such as
prediction of genomic variants or protein binding motifs located in repeat
regions, isoform expression quantification, metagenomics analysis, etc.
However, this task is very time-consuming and majority of mapping tools only
focus on one or a few best mapping locations. Results: We propose AMAS, an
alignment tool specialized in identifying all possible mapping locations of NGS
reads in a reference sequence. AMAS features an effective use of adaptive seeds
to speed up read mapping while preserving sensitivity. Specifically, an index
is designed to pre-store the locations of adaptive seeds in the reference
sequence, efficiently reducing the time for seed matching and partitioning. An
accurate filtration of adaptive seeds is further applied to substantially
tighten the candidate alignment space. As a result, AMAS runs several times
faster than other state-of-the-art read mappers while achieving similar
accuracy. Conclusions: AMAS provides a valuable resource to speed up the
important yet time-consuming task of identifying all mapping locations of NGS
reads. AMAS is implemented in C++ based on the SeqAn library and is freely
available at https://sourceforge.net/projects/ngsamas/. Keywords:
next-generation sequencing, read mapping, sequence alignment, adaptive seeds,
seed partition, filtration


  The Moran model with recombination is considered, which describes the
evolution of the genetic composition of a population under recombination and
resampling. There are $n$ sites (or loci), a finite number of letters (or
alleles) at every site, and we do not make any scaling assumptions. In
particular, we do not assume a diffusion limit. We consider the following
marginal ancestral recombination process. Let $S = \{1,...,n\}$ and $\mathcal
A=\{A_1, ..., A_m\}$ be a partition of $S$. We concentrate on the joint
probability of the letters at the sites in $A_1$ in individual $1$, $...$, and
at the sites in $A_m$ in individual $m$, where the individuals are sampled from
the current population without replacement. Following the ancestry of these
sites backwards in time yields a process on the set of partitions of $S$,
which, in the diffusion limit, turns into a marginalised version of the
$n$-locus ancestral recombination graph. With the help of an
inclusion-exclusion principle, we show that the type distribution corresponding
to a given partition may be represented in a systematic way, with the help of
so-called recombinators and sampling functions. The same is true of correlation
functions (known as linkage disequilibria in genetics) of all orders.
  We prove that the partitioning process (backward in time) is dual to the
Moran population process (forward in time), where the sampling function plays
the role of the duality function. This sheds new light on the work of
Bobrowski, Wojdyla, and Kimmel (2010). The result also leads to a closed system
of ordinary differential equations for the expectations of the sampling
functions, which can be translated into expected type distributions and
expected linkage disequilibria.


  Dynamic interspecific interactions are thought to contribute to fundamental
structures of ecological communities. Using the Yodzis framework of dominance
control versus founder control versus niche control, the importance of
intransitive loops in sub-communities in a source/sink or metacommunity
competition context is examined. It is found that a three-species sub-community
sampled from a species pool composed of nothing but dominance controlled
species pairs can result in a stable three species coexistence. It is shown
that if effect and response competition are balanced, the overall sub-community
will be persistent, anchored by one or more intransitive loops. It is proposed
that such an arrangement be referred to as an intransitive structure. It is
concluded that persistence of a sub-community (when species pairs in the
species pool are all dominance controlled) of more than three species is
possible only if it has an intransitive structure.


  Long-range and highly accurate de novo assembly from short-read data is one
of the most pressing challenges in genomics. Recently, it has been shown that
read pairs generated by proximity ligation of DNA in chromatin of living tissue
can address this problem. These data dramatically increase the scaffold
contiguity of assemblies and provide haplotype phasing information. Here, we
describe a simpler approach ("Chicago") based on in vitro reconstituted
chromatin. We generated two Chicago datasets with human DNA and used a new
software pipeline ("HiRise") to construct a highly accurate de novo assembly
and scaffolding of a human genome with scaffold N50 of 30 Mb. We also
demonstrated the utility of Chicago for improving existing assemblies by
re-assembling and scaffolding the genome of the American alligator. With a
single library and one lane of Illumina HiSeq sequencing, we increased the
scaffold N50 of the American alligator from 508 kb to 10 Mb. Our method uses
established molecular biology procedures and can be used to analyze any genome,
as it requires only about 5 micrograms of DNA as the starting material.


  Understanding the architectural subtleties of ecological networks, believed
to confer them enhanced stability and robustness, is a subject of outmost
relevance. Mutualistic interactions have been profusely studied and their
corresponding bipartite networks, such as plant-pollinator networks, have been
reported to exhibit a characteristic "nested" structure. Assessing the
importance of any given species in mutualistic networks is a key task when
evaluating extinction risks and possible cascade effects. Inspired in a
recently introduced algorithm --similar in spirit to Google's PageRank but with
a built-in non-linearity-- here we propose a method which --by exploiting their
nested architecture-- allows us to derive a sound ranking of species importance
in mutualistic networks. This method clearly outperforms other existing ranking
schemes and can become very useful for ecosystem management and biodiversity
preservation, where decisions on what aspects of ecosystems to explicitly
protect need to be made.


  Many life-history traits, like the age at maturity or adult longevity, are
important determinants of the generation time. For instance, semelparous
species whose adults reproduce once and die have shorter generation times than
iteroparous species that reproduce on several occasions. A shorter generation
time ensures a higher growth rate in stable environments where resources are in
excess, and is therefore a positively selected feature in this (rarely met)
situation. In a stable and limiting environment, all combination of traits (or
strategies) that produce the same number of viable offspring on average are
strictly neutral even when their generation times differ. We first study the
evolution of life-history strategies with different generation times in this
context, and show that those with the longest generation time represent the
most likely evolutionary outcomes. Indeed, strategies with longer generation
times generate fewer mutants per time unit, which makes them less likely to be
replaced within a given time period. This `turnover bias' inevitably exists and
favors the evolution of strategies with long generation times. Its real impact,
however, should depend on the strength and direction of other evolutionary
forces; selection for short generation times, for instance, may oppose turnover
bias. Likewise, the evolutionary outcome depends on the strength of such
selection and population size, comparably to other biases acting on the
occurrence of mutations.


  Evolutionary and Interpretive Archaeologies, edited by Ethan E. Cochrane and
Andrew Gardner, grew out of a seminar at the Institute for Archaeology at
University College London in 2007. It consists of 15 chapters by archaeologists
who self-identify themselves as practitioners who emphasize the benefits of
evolutionary or interpretive approaches to the study of the archaeological
record. While the authors' theoretical views are dichotomous, the editors' aim
for the book as a whole is not to expound on the differences between these two
kinds of archaeology but to bring forward a richer understanding of the
discipline and to highlight areas of mutual concern. Some chapters come across
as a bit of a sales pitch, but the majority of the contributions emphasize how
each approach can be productively used to address the goals of the other. The
book seeks to contribute to a mutually beneficial and more productive
discipline, and overall, it succeeds in this effort.


  We present a dual optimization concept of predicting optimal sequences as
well as optimal folds of off-lattice protein models in the context of
multi-scale modeling. We validate the utility of the recently introduced
hidden-force Monte Carlo optimization algorithm by finding significantly lower
energy folds for minimalist and detailed protein models than previously
reported. Further, we also find the protein sequence that yields the lowest
energy fold amongst all sequences for a given chain length and residue mixture.
In particular, for protein models with a binary sequence, we show that the
sequence-optimized folds form more compact cores than the lowest energy folds
of the historically fixed, Fibonacci-series sequences of chain lengths of 13,
21, 34, 55, and 89. We then extend our search algorithm to use UNRES, one of
the leading united-residue protein force fields. Our combined fold and sequence
optimization on three test proteins reveal an inherent bias in UNRES favoring
alpha helical structures even when secondary structure prediction clearly
suggests only beta sheets besides random coil, and virtually no helices. One
test in particular, a triple-stranded antiparallel beta-sheet protein domain,
demonstrates that by permutations of its sequence UNRES re-folds this structure
into a perfect alpha helix but, in fact, the helix is just an artefact of the
force field, the structure quickly unfolds in all-atom state-of-the-art
molecular dynamics simulation.


  RNA is a fundamental class of biomolecules that mediate a large variety of
molecular processes within the cell. Computational algorithms can be of great
help in the understanding of RNA structure-function relationship. One of the
main challenges in this field is the development of structure-prediction
algorithms, which aim at the prediction of the three-dimensional (3D) native
fold from the sole knowledge of the sequence. In a recent paper, we have
introduced a scoring function for RNA structure prediction. Here, we analyze in
detail the performance of the method, we underline strengths and shortcomings,
and we discuss the results with respect to state-of-the-art techniques. These
observations provide a starting point for improving current methodologies, thus
paving the way to the advances of more accurate approaches for RNA 3D structure
prediction.


  Allometry and growth rates of 8 forest species in the UK. The data were
collected from two United Kingdom woodlands - Wytham Woods and Alice Holt. Here
we present data from 582 individual trees of eight taxa in the form of summary
variables. In addition the raw data files containing the variables from which
the summary data were obtained. Large sample sizes with longitudinal data
spanning 22 years make these datasets useful for future studies concerned with
the way trees change in size and shape over their life-span. The allometric
relationships include (1) trunk diameter, (2) height, (3) crown height, (4)
crown radius and (5) trunk radial growth rate to (A) the light environment of
each tree and (B) diameter at breast height.


  Metagenomics is a powerful approach to study genetic content of environmental
samples that has been strongly promoted by NGS technologies. To cope with
massive data involved in modern metagenomic projects, recent tools [4, 39] rely
on the analysis of k-mers shared between the read to be classified and sampled
reference genomes. Within this general framework, we show in this work that
spaced seeds provide a significant improvement of classification accuracy as
opposed to traditional contiguous k-mers. We support this thesis through a
series a different computational experiments, including simulations of
large-scale metagenomic projects. Scripts and programs used in this study, as
well as supplementary material, are available from
http://github.com/gregorykucherov/spaced-seeds-for-metagenomics.


  Gene expression levels are important molecular quantitative traits that link
genotypes to molecular functions and fitness. In Drosophila, population-genetic
studies in recent years have revealed substantial adaptive evolution at the
genomic level. However, the evolutionary modes of gene expression have remained
controversial. Here we present evidence that adaptation dominates the evolution
of gene expression levels in flies. We show that 63% of the observed expression
divergence across seven Drosophila species are adaptive changes driven by
directional selection. Our results are derived from the variation of expression
within species and the time-resolved divergence across a family of related
species, using a new inference method for selection. We identify functional
classes of adaptively regulated genes, as well as sex-specific adaptation
occurring predominantly in males. Our analysis opens a new avenue to map
system-wide selection on molecular quantitative traits independently of their
genetic basis.


  The rise of trait-based ecology has led to an increased focus on the
distribution and dynamics of traits in communities. However, a general theory
of trait-based ecology, that can apply across different scales (e.g., species
that differ in size) and gradients (e.g., temperature), has yet to be
formulated. While research focused on metabolic and allometric scaling theory
provides the basis for such a theory it does not explicitly account for
differences traits within and across taxa, such as variation in the optimal
temperature for growth. Here we synthesize trait-based and metabolic scaling
approaches into a framework that we term Trait Drivers Theory or TDT. It shows
that the shape and dynamics of trait distributions can be uniquely linked to
fundamental drivers of community assembly and how the community will respond to
future drivers. To assess predictions and assumptions of TDT, we review several
theoretical studies, recent empirical studies spanning local and biogeographic
gradients. Further, we analyze how the shift in trait distributions influences
ecosystem productivity across an elevational gradient and a 140-year long
ecological experiment. We argue that our general TDT provides a baseline for
(i) recasting the predictions of ecological theories based on species richness
in terms of the shape of trait distributions; and (ii) integrating how specific
traits, including body size, and functional diversity scale up to influence the
dynamics of species assemblages across climatic gradients and how shifts in
functional composition influences ecosystem functioning. Further, it offers a
novel framework to integrate trait, metabolic/allometric, and species-richness
based approaches in order to build a more predictive functional biogeography to
show how assemblages of species have and will respond to climate change.


  We discuss the implementation of a multiscale biophysico-chemical model able
to cope with the main mechanisms underlying cumulative exocytosis in cells. The
model is based on a diffusion equation in the presence of external forces that
links calcium signaling and the biochemistry associated to the activity of
cytoskeletal-based protein motors. This multiscale model offers an excellent
quantitative spatio-temporal description of the cumulative exocytosis measured
by means of fluorescence experiments. We also review pre-existing models
reported in the literature on calcium waves, protein motor activation and
dynamics, and intracellular directed transport of vesicles. As an example of
the proposed model, we analyze the formation of the shield against polyspermy
in the early events of fertilization in sea urchin eggs.


  Why life persists at the edge of chaos is a question at the very heart of
evolution. Here we show that molecules taking part in biochemical processes
from small molecules to proteins are critical quantum mechanically. Electronic
Hamiltonians of biomolecules are tuned exactly to the critical point of the
metal-insulator transition separating the Anderson localized insulator phase
from the conducting disordered metal phase. Using tools from Random Matrix
Theory we confirm that the energy level statistics of these biomolecules show
the universal transitional distribution of the metal-insulator critical point
and the wave functions are multifractals in accordance with the theory of
Anderson transitions. The findings point to the existence of a universal
mechanism of charge transport in living matter. The revealed bio-conductor
material is neither a metal nor an insulator but a new quantum critical
material which can exist only in highly evolved systems and has unique material
properties.