  The value of a continuous character evolving on a phylogenetic tree is
commonly modelled as the location of a particle moving under one-dimensional
Brownian motion with constant rate. The Brownian motion model is best suited to
characters evolving under neutral drift or tracking an optimum that drifts
neutrally. We present a generalization of the Brownian motion model which
relaxes assumptions of neutrality and gradualism by considering increments to
evolving characters to be drawn from a heavy-tailed stable distribution (of
which the normal distribution is a specialized form). We describe Markov chain
Monte Carlo methods for fitting the model to biological data paying special
attention to ancestral state reconstruction, and study the performance of the
model in comparison with a selection of existing comparative methods, using
both simulated data and a database of body mass in 1,679 mammalian species. We
discuss hypothesis testing and model selection. The new model is well suited to
a stochastic process with a volatile rate of change in which biological
characters undergo a mixture of neutral drift and occasional evolutionary
events of large magnitude.


  The connection between chromatin nuclear organization and gene activity is
vividly illustrated by the observation that transcriptional coregulation of
certain genes appears to be directly influenced by their spatial proximity.
This fact poses the more general question of whether it is at all feasible that
the numerous genes that are coregulated on a given chromosome, especially those
at large genomic distances, might become proximate inside the nucleus. This
problem is studied here using steered molecular dynamics simulations in order
to enforce the colocalization of thousands of knowledge-based gene sequences on
a model for the gene-rich human chromosome 19. Remarkably, it is found that
most (~80%) gene pairs can be brought simultaneously into contact. This is made
possible by the low degree of intra-chromosome entanglement and the large
number of cliques in the gene coregulatory network, that is the many groups of
genes that are all mutually coregulated. The constrained conformations for the
model chromosome 19 are further shown to be organised in spatial macrodomains
that are similar to those inferred from recent HiC measurements. The findings
indicate that gene coregulation and colocalization are largely compatible and
that this relationship can be exploited to draft the overall spatial
organization of the chromosome in vivo. The more general validity and
implications of these findings could be investigated by applying to other
eukaryotic chromosomes the general and transferable computational strategy
introduced here.


  The comment is intended to answer the criticism presented on `Steady-state
fluctuations of a genetic feedback loop: an exact solution' [J. Chem. Phys.
{\bf 137}, 035104 (2012).] and provides the missing component for the complete
analytic solutions to the author's model.


  Genetic and pharmacological perturbation experiments, such as deleting a gene
and monitoring gene expression responses, are powerful tools for studying
cellular signal transduction pathways. However, it remains a challenge to
automatically derive knowledge of a cellular signaling system at a conceptual
level from systematic perturbation-response data. In this study, we explored a
framework that unifies knowledge mining and data mining approaches towards the
goal. The framework consists of the following automated processes: 1) applying
an ontology-driven knowledge mining approach to identify functional modules
among the genes responding to a perturbation in order to reveal potential
signals affected by the perturbation; 2) applying a graph-based data mining
approach to search for perturbations that affect a common signal with respect
to a functional module, and 3) revealing the architecture of a signaling system
organize signaling units into a hierarchy based on their relationships.
Applying this framework to a compendium of yeast perturbation-response data, we
have successfully recovered many well-known signal transduction pathways; in
addition, our analysis have led to many hypotheses regarding the yeast signal
transduction system; finally, our analysis automatically organized perturbed
genes as a graph reflecting the architect of the yeast signaling system.
Importantly, this framework transformed molecular findings from a gene level to
a conceptual level, which readily can be translated into computable knowledge
in the form of rules regarding the yeast signaling system, such as "if genes
involved in MAPK signaling are perturbed, genes involved in pheromone responses
will be differentially expressed".


  We develop a way of simulating disease spread in networks faster at the cost
of some accuracy. Instead of a discrete event simulation (DES) we use a
discrete time simulation. This aggregates events into time periods. We prove a
bound on the accuracy attained. We also discuss the choice of step size and do
an analytical comparison of the computational costs. Our error bound concept
comes from the theory of numerical methods for SDEs and the basic proof
structure comes from the theory of numerical methods for ODEs.


  To tackle the exponentially increasing throughput of Next-Generation
Sequencing (NGS), most of the existing short-read aligners can be configured to
favor speed in trade of accuracy and sensitivity. SOAP3-dp, through leveraging
the computational power of both CPU and GPU with optimized algorithms, delivers
high speed and sensitivity simultaneously. Compared with widely adopted
aligners including BWA, Bowtie2, SeqAlto, GEM and GPU-based aligners including
BarraCUDA and CUSHAW, SOAP3-dp is two to tens of times faster, while
maintaining the highest sensitivity and lowest false discovery rate (FDR) on
Illumina reads with different lengths. Transcending its predecessor SOAP3,
which does not allow gapped alignment, SOAP3-dp by default tolerates alignment
similarity as low as 60 percent. Real data evaluation using human genome
demonstrates SOAP3-dp's power to enable more authentic variants and longer
Indels to be discovered. Fosmid sequencing shows a 9.1 percent FDR on newly
discovered deletions. SOAP3-dp natively supports BAM file format and provides a
scoring scheme same as BWA, which enables it to be integrated into existing
analysis pipelines. SOAP3-dp has been deployed on Amazon-EC2, NIH-Biowulf and
Tianhe-1A.


  In the natural world, there are many insect species whose individual members
have a life history that takes them through two stages, immature and mature.
Moreover, the rates of survival, development, and reproduction almost always
depend on age, size, or development stage. Keeping this in mind, in this paper,
a three species crop-pest-natural enemy food chain model with two stages for
natural enemies is investigated. Using characteristic equations, a set of
sufficient conditions for local asymptotic stability of all the feasible
equilibria is obtained. Moreover, using approach as in (Beretta and Kuang,
2002), the possibility of the existence of a Hopf bifurcation for the interior
equilibrium with respect to maturation delay is explored, which shows that the
maturation delay plays an important role in the dynamical behavior of three
species system. Also obtain some threshold values of maturation delay for the
stability-switching of the particular system. In succession, using the normal
form theory and center manifold argument, we derive the explicit formulas which
determine the stability and direction of bifurcating periodic solutions.
Finally, a numerical simulation for supporting the theoretical analysis is
given.


  FoF1-ATP synthase catalyzes the synthesis of adenosine triphosphate (ATP).
The F1 portion can be stripped from the membrane-embedded Fo portion of the
enzyme. F1 acts as an ATP hydrolyzing enzyme, and ATP hydrolysis is associated
with stepwise rotation of the gamma and epsilon subunits of F1. This rotary
motion was studied in great detail for the last 15 years using single F1 parts
attached to surfaces. Subunit rotation of gamma was monitored by
videomicroscopy of bound fluorescent actin filaments, nanobeads or nanorods, or
single fluorophores. Alternatively, we applied single-molecule F\"orster
resonance energy transfer (FRET) to monitor subunit rotation in the holoenzyme
FoF1-ATP synthase which was reconstituted in liposomes. Now we aim to extend
the observation times of single FRET-labeled F1 in solution using a modified
version of the anti-Brownian electrokinetic trap (ABELtrap) invented by A. E.
Cohen and W. E. Moerner. We used Monte Carlo simulations to reveal that
stepwise FRET efficiency changes can be analyzed by Hidden Markov Models even
at the limit of a low signal-to-background ratio that was expected due to high
background count rates caused by the microfluidics of the ABELtrap.


  An ecological flow network is a weighted directed graph in which nodes are
species, edges are "who eats whom" relationships and weights are rates of
energy or nutrients transfer between species. Allometric scaling is a
ubiquitous feature for flow systems like river basins, vascular networks and
food webs. By "ecological network analysis" method, we can reveal the hidden
allometry directly on the original flow networks without cutting edges. On the
other hand, dissipation law, which is another significant scaling relationship
between the energy dissipation (respiration) and the throughflow of any species
is also discovered on the collected flow networks. Interestingly, the exponents
of allometric law ($\eta$) and the dissipation law ($\gamma$) have a strong
connection for both empirical and simulated flow networks. The dissipation law
exponent $\gamma$ rather than the topology of the network is the most important
ingredient to the allometric exponent $\eta$. By reinterpreting $\eta$ as the
inequality of species impacts (direct and indirect influences) to the whole
network along all energy flow pathways but not the energy transportation
efficiency, we found that as $\gamma$ increases, the relative energy loss of
large nodes (with high throughflow) increases, $\eta$ decreases, and the
inequality of the whole flow network as well as the relative importance of
large species decreases. Therefore, flow structure and thermodynamic constraint
are connected.


  GC-content, the ratio of guanine and cytosine bases in an entire nucleotide
sequence, and palindromic sequences are unique for every organism due to
genomic evolution. The goals of our research was to establish a correlation
between GC-content and palindromic densities in wild-type viral and
randomly-generated genomes. Forty viral genomes were downloaded from GenBank
and their GC-ratios and palindromic densities were calculated and plotted using
Mathematica. The palindromic densities-by-GC-ratios plot of randomly generated
sequences (palindromic density curve) exhibited a quadratic relationship and
was superimposed over the viral genome plot. It was observed that the viral
plots followed the curvature of the random sequences' quadratic curve,
signifying a directly proportional relationship between GC-content and
palindrome density in viral genomes. However, because viral genomes require
certain non-palindromic sequences to function, the palindromic densities of
most wild-type genomes were under the palindromic density curve. The variance
in palindrome densities of wild-type genomes in respect to the random
sequences' quadratic curve may be examined to determine evolutionary traits in
genomes. A better understanding of viral palindromic densities and GC-ratios
would help in understanding conserved secondary RNA structures in viral genomes
and future drug discovery. In addition, certain viral genomes were found to be
viable recombinant viruses, which are used in gene therapy.


  In ecology, prey switching refers to a predator's adaptive change of habitat
or diet in response to prey abundance. In this paper, we study piecewise-smooth
models of predator-prey interactions with a linear trade-off in a predator's
prey preference. We consider optimally foraging predators and derive a model
for a 1 predator-2 prey interaction with a tilted switching manifold between
the two sides of discontinuous vector fields. We show that the 1 predator-2
prey system undergoes a novel adding-sliding-like (center to two-part periodic
orbit; "C2PO") bifurcation in which the prey ratio transitions from constant to
time-dependent. Further away from the bifurcation point, the period of the
oscillating prey ratio period doubles, suggesting a possible cascade to chaos.
We compare our model predictions with data and demonstrate that we successfully
capture the periodicity in the ratio between the predator's preferred and
alternative prey types in data on freshwater plankton. Our study suggests that
it is useful to investigate prey ratio as a possible indicator of how
population dynamics can be influenced by ecosystem diversity.


  The family Nymphalidae is the largest family within the true butterflies and
has been used to develop hypotheses explaining evolutionary interactions
between plants and insects. Theories of insect and hostplant dynamics predict
accelerated diversification in some scenarios. We investigated whether
phylogenetic uncertainty affects a commonly used method (MEDUSA, modelling
evolutionary diversity using stepwise AIC) for estimating shifts in
diversification rates in lineages of the family Nymphalidae, by extending the
method to run across a random sample of phylogenetic trees from the posterior
distribution of a Bayesian run. We found that phylogenetic uncertainty greatly
affects diversification rate estimates. Different trees from the posterior
distribution can give diversification rates ranging from high values to almost
zero for the same clade, and for some clades both significant rate increase and
decrease were estimated. Only three out of 13 significant shifts found on the
maximum credibility tree were consistent across more than 95% of the trees from
the posterior: (i) accelerated diversification for Solanaceae feeders in the
tribe Ithomiini; (ii) accelerated diversification in the genus Charaxes, and
(iii) deceleration in the Danaina. By using the binary speciation and
extinction model (BISSE), we found that a hostplant shift to Solanaceae or a
codistributed character is responsible for the increase in diversification rate
in Ithomiini, and the result is congruent with the diffuse cospeciation
hypothesis. A shift to Apocynaceae is not responsible for the slowdown of
diversification in Danaina. Our results show that taking phylogenetic
uncertainty into account when estimating diversification rate shifts is of
great importance, and relying on the maximum credibility tree alone potentially
can give erroneous results.


  Photosystem II (PSII) and its associated light-harvesting complex II (LHCII)
are highly concentrated in the stacked grana regions of photosynthetic
thylakoid membranes. Within the membrane, PSII-LHCII supercomplexes can be
arranged in disordered packings, ordered arrays, or mixtures thereof. The
physical driving forces underlying array formation are unknown, complicating
attempts to determine a possible functional role for arrays in regulating light
harvesting or energy conversion efficiency. Here we introduce a coarse-grained
model of protein interactions in coupled photosynthetic membranes, focusing on
just two particle types that feature simple shapes and potential energies
motivated by structural studies. Reporting on computer simulations of the
model's equilibrium fluctuations, we demonstrate its success in reproducing
diverse structural features observed in experiments, including extended
PSII-LHCII arrays. Free energy calculations reveal that the appearance of
arrays marks a phase transition from the disordered fluid state to a
system-spanning crystal, which can easily be arrested by thermodynamic
constraints or slow dynamics. The region of fluid-crystal coexistence is broad,
encompassing much of the physiologically relevant parameter regime. Our results
suggest that grana membranes lie at or near phase coexistence, conferring
significant structural and functional flexibility to this densely packed
membrane protein system.


  Motivation: While the majority of gene histories found in a clade of
organisms are expected to be generated by a common process (e.g. the coalescent
process), it is well-known that numerous other coexisting processes (e.g.
horizontal gene transfers, gene duplication and subsequent
neofunctionalization) will cause some genes to exhibit a history quite distinct
from those of the majority of genes. Such "outlying" gene trees are considered
to be biologically interesting and identifying these genes has become an
important problem in phylogenetics.
  Results: We propose and implement KDETREES, a nonparametric method of
estimating distributions of phylogenetic trees, with the goal of identifying
trees which are significantly different from the rest of the trees in the
sample. Our method compares favorably with a similar recently-published method,
featuring an improvement of one polynomial order of computational complexity
(to quadratic in the number of trees analyzed), with simulation studies
suggesting only a small penalty to classification accuracy. Application of
KDETREES to a set of Apicomplexa genes identified several unreliable sequence
alignments which had escaped previous detection, as well as a gene
independently reported as a possible case of horizontal gene transfer. We also
analyze a set of Epichloe genes, fungi symbiotic with grasses, successfully
identifying a contrived instance of paralogy.
  Availability: Our method for estimating tree distributions and identifying
outlying trees is implemented as the R package KDETREES, and is available for
download from CRAN.


  Protein-protein interaction networks are useful for studying human diseases
and to look for possible health care through a holistic approach. Networks are
playing an increasing and important role in the understanding of physiological
processes such as homeostasis, signaling, spatial and temporal organizations,
and pathological conditions. In this article we show the complex system of
interactions determined by human Sirtuins (Sirt) largely involved in many
metabolic processes as well as in different diseases. The Sirtuin family
consists of seven homologous Sirt-s having structurally similar cores but
different terminal segments, being rather variable in length and/or
intrinsically disordered. Many studies have determined their cellular location
as well as biological functions although molecular mechanisms through which
they act are actually little known. Therefore, the aim of this work was to
define, explore and understand the Sirtuin-related human interactome. As a
first step, we have integrated the experimentally determined protein-protein
interactions of the Sirtuin-family as well as their first and second neighbors
to a Sirtuin-related sub-interactome. Our data showed that the second-neighbor
network of Sirtuins encompasses 25% of the entire human interactome, exhibits a
scale-free degree distribution and interconnectedness among top degree nodes.
Moreover, the Sirtuin sub interactome showed a modular structure around the
core comprising mixed functions. Finally, we extracted from the Sirtuin
sub-interactome subnets related to cancer, aging and post-translational
modifications for information on key nodes and topological space of the subnets
in Sirt family network.


  Spatial heterogeneity is a hallmark of living systems, even at the molecular
scale in individual cells. A key example is the partitioning of membrane-bound
proteins via lipid domain formation or cytoskeleton-induced corralling. Yet the
impact of this spatial heterogeneity on biochemical signaling processes is
poorly understood. Here we demonstrate that partitioning improves the
reliability of biochemical signaling. We exactly solve a stochastic model
describing a ubiquitous motif in membrane signaling. The solution reveals that
partitioning improves signaling reliability via two effects: it moderates the
non-linearity of the switching response, and it reduces noise in the response
by suppressing correlations between molecules. An optimal partition size arises
from a trade-off between minimizing the number of proteins per partition to
improve signaling reliability and ensuring sufficient proteins per partition to
maintain signal propagation. The predicted optimal partition size agrees
quantitatively with experimentally observed systems. These results persist in
spatial simulations with explicit diffusion barriers. Our findings suggest that
molecular partitioning is not merely a consequence of the complexity of
cellular substructures, but also plays an important functional role in cell
signaling.


  A goal in the kinetic characterization of a macromolecular system is the
description of its slow relaxation processes, involving (i) identification of
the structural changes involved in these processes, and (ii) estimation of the
rates or timescales at which these slow processes occur. Most of the approaches
to this task, including Markov models, Master-equation models, and kinetic
network models, start by discretizing the high-dimensional state space and then
characterize relaxation processes in terms of the eigenvectors and eigenvalues
of a discrete transition matrix. The practical success of such an approach
depends very much on the ability to finely discretize the slow order
parameters. How can this task be achieved in a high-dimensional configuration
space without relying on subjective guesses of the slow order parameters? In
this paper, we use the variational principle of conformation dynamics to derive
an optimal way of identifying the "slow subspace" of a large set of prior order
parameters - either generic internal coordinates (distances and dihedral
angles), or a user-defined set of parameters. It is shown that a method to
identify this slow subspace exists in statistics: the time-lagged independent
component analysis (TICA). Furthermore, optimal indicators-order parameters
indicating the progress of the slow transitions and thus may serve as reaction
coordinates-are readily identified. We demonstrate that the slow subspace is
well suited to construct accurate kinetic models of two sets of molecular
dynamics simulations, the 6-residue fluorescent peptide MR121-GSGSW and the
30-residue natively disordered peptide KID. The identified optimal indicators
reveal the structural changes associated with the slow processes of the
molecular system under analysis.


  Rewards and penalties are common practical tools that can be used to promote
cooperation in social institutions. The evolution of cooperation under reward
and punishment incentives in joint enterprises has been formalized and
investigated, mostly by using compulsory public good games. Recently, Sasaki et
al. (2012, Proc Natl Acad Sci USA 109:1165-1169) considered optional
participation as well as institutional incentives and described how the
interplay between these mechanisms affects the evolution of cooperation in
public good games. Here, we present a full classification of these evolutionary
dynamics. Specifically, whenever penalties are large enough to cause the
bi-stability of both cooperation and defection in cases in which participation
in the public good game is compulsory, these penalties will ultimately result
in cooperation if participation in the public good game is optional. The global
stability of coercion-based cooperation in this optional case contrasts
strikingly with the bi-stability that is observed in the compulsory case. We
also argue that optional participation is not so effective at improving
cooperation under rewards.


  The adaptation of large asexual populations is hampered by the competition
between independently arising beneficial mutations in different individuals,
which is known as clonal interference. Fisher and Muller proposed that
recombination provides an evolutionary advantage in large populations by
alleviating this competition. Based on recent progress in quantifying the speed
of adaptation in asexual populations undergoing clonal interference, we present
a detailed analysis of the Fisher-Muller mechanism for a model genome
consisting of two loci with an infinite number of beneficial alleles each and
multiplicative fitness effects. We solve the infinite population dynamics
exactly and show that, for a particular, natural mutation scheme, the speed of
adaptation in sexuals is twice as large as in asexuals. Guided by the infinite
population result and by previous work on asexual adaptation, we postulate an
expression for the speed of adaptation in finite sexual populations that agrees
with numerical simulations over a wide range of population sizes and
recombination rates. The ratio of the sexual to asexual adaptation speed is a
function of population size that increases in the clonal interference regime
and approaches 2 for extremely large populations. The simulations also show
that the imbalance between the numbers of accumulated mutations at the two loci
is strongly suppressed even by a small amount of recombination. The
generalization of the model to an arbitrary number $L$ of loci is briefly
discussed. If each offspring samples the alleles at each locus from the gene
pool of the whole population rather than from two parents, the ratio of the
sexual to asexual adaptation speed is approximately equal to $L$ in large
populations. A possible realization of this scenario is the reassortment of
genetic material in RNA viruses with $L$ genomic segments.


  Motivation: Human cancer is caused by the accumulation of somatic mutations
in tumor suppressors and oncogenes within the genome. In the case of oncogenes,
recent theory suggests that there are only a few key "driver" mutations
responsible for tumorigenesis. As there have been significant pharmacological
successes in developing drugs that treat cancers that carry these driver
mutations, several methods that rely on mutational clustering have been
developed to identify them. However, these methods consider proteins as a
single strand without taking their spatial structures into account. We propose
a new methodology that incorporates protein tertiary structure in order to
increase our power when identifying mutation clustering.
  Results: We have developed a novel algorithm, iPAC: identification of Protein
Amino acid Clustering, for the identification of non-random somatic mutations
in proteins that takes into account the three dimensional protein structure. By
using the tertiary information, we are able to detect both novel clusters in
proteins that are known to exhibit mutation clustering as well as identify
clusters in proteins without evidence of clustering based on existing methods.
For example, by combining the data in the Protein Data Bank (PDB) and the
Catalogue of Somatic Mutations in Cancer, our algorithm identifies new
mutational clusters in well known cancer proteins such as KRAS and PI3KCa.
Further, by utilizing the tertiary structure, our algorithm also identifies
clusters in EGFR, EIF2AK2, and other proteins that are not identified by
current methodology.


  In many biological applications, we would like to be able to computationally
predict mutational effects on affinity in protein-protein interactions.
However, many commonly used methods to predict these effects perform poorly in
important test cases. In particular, the effects of multiple mutations,
non-alanine substitutions, and flexible loops are difficult to predict with
available tools and protocols. We present here an existing method applied in a
novel way to a new test case; we interrogate affinity differences resulting
from mutations in a host-virus protein-protein interface. We use steered
molecular dynamics (SMD) to computationally pull the machupo virus (MACV) spike
glycoprotein (GP1) away from the human transferrin receptor (hTfR1). We then
approximate affinity using the maximum applied force of separation and the area
under the force-versus-distance curve. We find, even without the rigor and
planning required for free energy calculations, that these quantities can
provide novel biophysical insight into the GP1/hTfR1 interaction. First, with
no prior knowledge of the system we can differentiate among wild type and
mutant complexes. Moreover, we show that this simple SMD scheme correlates well
with relative free energy differences computed via free energy perturbation.
Second, although the static co-crystal structure shows two large
hydrogen-bonding networks in the GP1/hTfR1 interface, our simulations indicate
that one of them may not be important for tight binding. Third, one viral site
known to be critical for infection may mark an important evolutionary
suppressor site for infection-resistant hTfR1 mutants. Finally, our approach
provides a framework to compare the effects of multiple mutations, individually
and jointly, on protein-protein interactions.


  The periodic transference of nucleotide strings in bacterial and archaeal
complete genomes is investigated by using the metric representation and the
recurrence plot method. The generated periodic correlation structures exhibit
four kinds of fundamental transferring characteristics: a single increasing
period, several increasing periods, an increasing quasi-period and almost
noincreasing period. The mechanism of the periodic transference is further
analyzed by determining all long periodic nucleotide strings in the bacterial
and archaeal complete genomes and is explained as follows: both the repetition
of basic periodic nucleotide strings and the transference of non-periodic
nucleotide strings would form the periodic correlation structures with
approximately the same increasing periods.


  Transgenerational inheritance of a trait is presumably affected by both
genetic and environmental factors but remains poorly understood. We studied the
effect of genetic polymorphisms on transgenerational inheritance of yeast
segregants that were derived from a cross between a laboratory strain and a
wild strain of Saccharomyces cerevisiae. For each SNP analyzed, the parental
allele present in less than half of the segregants panel was called the minor
allele (MA). We found a nonrandom distribution of MAs in the segregants,
indicating natural selection. We compared segregants with high MA content (MAC)
relative to those with less and found a more dramatic shortening of the lag
phase length for the high MAC group in response to 14 days of ethanol training.
Also, the short lag phase as acquired and epigenetically memorized by ethanol
training was more dramatically lost after 7 days of recovery in ethanol free
medium for the high MAC group. Sodium chloride treatment produced similar
observations. Using public datasets, we found MAC linkage to mRNA expression of
hundreds of genes. Finally, we found preferential effect of MAC on traits with
high number of known additive quantitative trait loci (QTLs). These results
provide evidence for the slightly deleterious nature of most MAs and a lower
capacity to maintain inheritance of traits in individuals or cells with greater
MAC, which have implications for disease prevention and treatment and the
"missing heritability" problem in complex traits and diseases.


  This paper further develops work from a previous article which introduced a
new way of modelling evolutionary game models with an emphasis on ecological
realism, concerned with how ecological factors determine payoffs in
evolutionary games. The current paper is focused on the impact of selectively
neutral factors (i.e. those that are the same for all strategies), such as
background fitness and a neutral density dependence on game dynamics. In the
previous paper background mortality was modelled by post-reproductive mortality
describing aggregated mortality between game interactions. In this paper a
novel approach to the modelling of background fitness, based on different
occurrence rates (intensities) of game interactions and background events is
presented. It is shown that the new method is more realistic and natural for
dynamic models than the method used in the previous paper. The previous method
in effect produces biased trajectories, but it is compatible with the new
approach at the restpoints. Also in this paper a rigorous stability analysis of
the restpoints of the dynamics is presented. It is shown that the stability or
instability of the restpoints, can be explained by a mechanism of
eco-evolutionary feedback. The results obtained are illustrated by a Hawk-Dove
game example.


  Open-source software is a complex system; its development depends on the
self-coordinated action of a large number of agents. This study follows the
size of the building blocks, called "packages", of the Ubuntu Linux operating
system over its entire history. The analysis reveals a multiplicative diffusion
process, constrained by size-dependent bounds, driving the dynamics of the
package-size distribution. A formalization of this into a quantitative model is
able to match the data without relying on any adjustable parameters, and
generates definite predictions. Finally, we formulate the hypothesis that a
similar non-stationary mechanism could be shaping the distribution of mammal
body sizes.


  Understanding the mechanisms underlying ion channel function from the
atomic-scale requires accurate ab initio modelling as well as careful
experiments. Here, we present a density functional theory (DFT) study of the
ion channel gramicidin A, whose inner pore conducts only monovalent cations and
whose conductance has been shown to depend on the side chains of the amino
acids in the channel. We investigate the ground-state geometry and electronic
properties of the channel in vacuum, focusing on their dependence on the side
chains of the amino acids. We find that the side chains affect the ground state
geometry, while the electrostatic potential of the pore is independent of the
side chains. This study is also in preparation for a full, linear scaling DFT
study of gramicidin A in a lipid bilayer with surrounding water. We demonstrate
that linear scaling DFT methods can accurately model the system with reasonable
computational cost. Linear scaling DFT allows ab initio calculations with
10,000 to 100,000 atoms and beyond, and will be an important new tool for
biomolecular simulations.


  Many of the chemical reactions that take place within a living cell are
irreversible. Due to evolutionary pressures, the number of allowable reactions
within these systems are highly constrained and thus the resulting metabolic
networks display considerable asymmetry. In this paper, we explore possible
evolutionary factors pertaining to the reduced symmetry observed in these
networks, and demonstrate the important role environmental variability plays in
shaping their structural organization. Interpreting the returnability index as
an equilibrium constant for a reaction network in equilibrium with a
hypothetical reference system, enables us to quantify the extent to which a
metabolic network is in disequilibrium. Further, by introducing a new directed
centrality measure via an extension of the subgraph centrality metric to
directed networks, we are able to characterise individual metabolites by their
participation within metabolic pathways. To demonstrate these ideas, we study
116 metabolic networks of bacteria. In particular, we find that the equilibrium
constant for the metabolic networks decreases significantly in-line with
variability in bacterial habitats, supporting the view that environmental
variability promotes disequilibrium within these biochemical reaction systems.


  Competitive interactions represent one of the driving forces behind evolution
and natural selection in biological and sociological systems. For example,
animals in an ecosystem may vie for food or mates; in a market economy, firms
may compete over the same group of customers; sensory stimuli may compete for
limited neural resources in order to enter the focus of attention. Here, we
derive rules based on the spectral properties of the network governing the
competitive interactions between groups of agents organized in networks. In the
scenario studied here the winner of the competition, and the time needed to
prevail, essentially depend on the way a given network connects to its
competitors and on its internal structure. Our results allow assessing the
extent to which real networks optimize the outcome of their interaction, but
also provide strategies through which competing networks can improve on their
situation. The proposed approach is applicable to a wide range of systems that
can be modeled as networks.


  Extensive divergence of transcription factor binding in Drosophila embryos
with highly conserved gene expression


  Soliton rate equations are based on non-Kolmogorovian models of probability
and naturally include autocatalytic processes. The formalism is not widely
known but has great unexplored potential for applications to systems
interacting with environments. Beginning with links of contextuality to
non-Kolmogorovity we introduce the general formalism of soliton rate equations
and work out explicit examples of subsystems interacting with environments. Of
particular interest is the case of a soliton autocatalytic rate equation
coupled to a linear conservative environment, a formal way of expressing
seasonal changes. Depending on strength of the system-environment coupling we
observe phenomena analogous to hibernation or even complete blocking of decay
of a population.


  Dominant negative mutations in the insulin gene are the second most common
cause of permanent neonatal diabetes. However, variation in severity and
penetrance of neonatal diabetes, as in other complex genetic diseases, cannot
be accounted for by known disease mutations. In a novel approach to this
problem, we have utilized the genetic tools available in Drosophila to model
the effects of the C96Y mutation, a cysteine to tyrosine mutation in the
insulin protein that can cause permanent neonatal diabetes in humans. This
mutation, which disrupts a disulfide bridge in the proinsulin molecule, has
been shown to lead to partial protein unfolding and aggregation in the
endoplasmic reticulum. It is thought to induce beta cell death in humans and
mice through endoplasmic reticulum stress-mediated apoptosis. We employed the
UAS/GAL4 system to create a stable fly laboratory stock expressing human
insulin with the C96Y mutation (INSC96Y) in the Drosophila eye. We crossed
these transgenic INSC96Y flies to 36 isogenic lines derived from a wild
population. The F1 flies displayed a disrupted eye development phenotype, with
both reduced eye size and irregularity in ommatidia, which varied between lines
and with sex. By backcrossing F1 offspring to their wild parents, we were able
to analyze the effects of naturally occurring genetic variation on the INSC96Y
phenotype in flies both heterozygous and homozygous for wild third chromosomes.
We observed significant differences between lines in terms of eye size, eye
shape, the presence or absence of lesions, and degree of ommatidial structure
present, indicating that an approach that takes into account quantitative
variation in the genetic background can yield valuable insights into the
phenotypic severity and penetrance of a complex genetic disease.


  The classical approaches to the modeling of sex ratio evolution can be
divided into two classes. The first class contains the static strategic models
related to the Dusing Fisher Shaw Mohler fitness measure, based on the
reproductive value of the offspring of the focal female. The second class
contains the population genetic models focused on the dynamics of the allele
frequencies. The approaches are not fully compatible because the strategic
models disregard the role of the male individuals as the passive carriers of
the strategy genes. In the previous two papers in this cycle, a new model
combining the strategic analysis with more rigorous genetics was presented. The
new model shows that sex ratio self-regulation is a multistage complex process
which can be regarded as an example of multilevel selection. One of the
elements of this process is the dynamic equilibrium between male and female
gene carriers associated with convergence of the dynamics to the manifold
termed the male subpopulation equilibrium (MSE). This paper attempts to explain
this phenomenon and analyze its properties. We show that the MSE phenomenon
affects every stage of sex ratio self-regulation (Lemmas 1-4). The MSE plays a
crucial role in synchronizing two levels of selection in the double-level
selection process. We show that the MSE condition can be generalized as an
interesting synergistic property allowing for the estimation of the primary sex
ratio of the entire population according to the state of some arbitrarily
chosen subpopulation (Lemma 5). We also show that the classical Dusing Fisher
Shaw Mohler fitness measure is a biased approximation of the new approach, but
that it produces compatible strategic predictions (Lemma 6).


  One of the key phenomena in the adaptive immune response to infection and
immunization is affinity maturation, during which antibody genes are mutated
and selected, typically resulting in a substantial increase in binding affinity
to the eliciting antigen. Advances in technology on several fronts have made it
possible to clone large numbers of heavy-chain light-chain pairs from
individual B cells and thereby identify whole sets of clonally related
antibodies. These collections could provide the information necessary to
reconstruct their own history--the sequence of changes introduced into the
lineage during the development of the clone--and to study affinity maturation
in detail. But the success of such a program depends entirely on accurately
inferring the founding ancestor and the other unobserved intermediates.
  I have developed an Empirical Bayes method that allows one to compute the
posterior distribution over ancestors, thereby giving a thorough accounting of
the uncertainty inherent in the reconstruction. I demonstrate the application
of this method on heavy-chain and light-chain clones, assess the reliability of
the inference, and discuss the sources of uncertainty.


  Discovering all the genetic causes of a phenotype is an important goal in
functional genomics. In this paper we combine an experimental design for
multiple independent detections of the genetic causes of a phenotype, with a
high-throughput sequencing analysis that maximizes sensitivity for
comprehensively identifying them. Testing this approach on a set of 24 mutant
strains generated for a metabolic phenotype with many known genetic causes, we
show that this pathway-based phenotype sequencing analysis greatly improves
sensitivity of detection compared with previous methods, and reveals a wide
range of pathways that can cause this phenotype. We demonstrate our approach on
a metabolic re-engineering phenotype, the PEP/OAA metabolic node in E. coli,
which is crucial to a substantial number of metabolic pathways and under
renewed interest for biofuel research. Out of 2157 mutations in these strains,
pathway-phenoseq discriminated just five gene groups (12 genes) as
statistically significant causes of the phenotype. Experimentally, these five
gene groups, and the next two high-scoring pathway-phenoseq groups, either have
a clear connection to the PEP metabolite level or offer an alternative path of
producing oxaloacetate (OAA), and thus clearly explain the phenotype. These
high-scoring gene groups also show strong evidence of positive selection
pressure, compared with strictly neutral selection in the rest of the genome.


  Modern developments in population dynamics emphasize the role of the turnover
of individuals. In the new approaches stable population size is a dynamic
equilibrium between different mortality and fecundity factors instead of an
arbitrary fixed carrying capacity. The latest replicator dynamics models assume
that regulation of the population size acts through feedback driven by density
dependent juvenile mortality. Here, we consider a simplified model to extract
the properties of this approach. We show that at the stable population size,
the structure of the frequency dependent evolutionary game emerges. Turnover of
individuals induces a lottery mechanism where for each nest site released by a
dead adult individual a single newborn is drawn from the pool of newborn
candidates. This frequency dependent selection leads toward the strategy
maximizing the number of newborns per adult death. However, multiple strategies
can maximize this value. Among them, the strategy with the greatest mortality
(which implies greater instantaneous growth rate) is selected. This result is
important for the discussion about universal fitness measures and which
parameters are maximized by natural selection. This is related to the fitness
measures R0 and r, because the number of newborns per single dead individual
equals lifetime production of newborn R0 in models without ageing. Our model
suggests the existence of another fitness measure which is the combination of
R0 and r. According to the nest site lottery mechanism, at stable population
size, selection favours strategies with the greatest r from those with the
greatest R0.


  Hybrid-Lambda is a software package that simulates gene trees under Kingman
or two Lambda-coalescent processes within species networks or species trees. It
is written in C++, and re- leased under GNU General Public License (GPL)
version 3. Users can modify and make new dis- tribution under the terms of this
license. For details of this license, visit http://www.gnu.org/licenses/.
Hybrid Lambda is available at https://code.google.com/p/hybrid-lambda.


  Anthropometric measurements such as weight, stature (height), and body mass
index (BMI) provide reliable indicators of children's growth. The 2000 CDC
growth charts are the national standards in the United States for these
important measures. But these growth charts were generated using data from
1963-1994. To understand the growth patterns of US children since 1994, we
generate weight-for-age, stature-for-age and BMI-for-age percentile curves for
both boys and girls aged 2-20 through the methods used to generate the 2000 CDC
growth charts. Our datasets are from the National Health and Nutrition
Examination Survey (NHANES) for years 1999-2010 and and from NorthShore
University HealthSystem's Enterprise Data Warehouse (NS-EDW) for years
2006-2012. The weight and BMI percentile curves generated from NS-EDW and
NHANES data differ substantially from the CDC percentile curves, while those
for stature do not differ substantially. We conclude that the population weight
and BMI values of US children in recent years have increased significantly
since 2000 and the 2000 CDC growth charts may no longer be applicable to the
current population of US children. Our charts poignantly reveals the increasing
obesity of American children.


  Intrapatient HIV-1 evolution is dominated by selection on the protein level
in the arms race with the adaptive immune system. When cytotoxic CD8+ T-cells
or neutralizing antibodies target a new epitope, the virus often escapes via
nonsynonymous mutations that impair recognition. Synonymous mutations do not
affect this interplay and are often assumed to be neutral. We analyze
longitudinal intrapatient data from the C2-V5 part of the envelope gene (env)
and observe that synonymous derived alleles rarely fix even though they often
reach high frequencies in the viral population. We find that synonymous
mutations that disrupt base pairs in RNA stems flanking the variable loops of
gp120 are more likely to be lost than other synonymous changes, hinting at a
direct fitness effect of these stem-loop structures in the HIV-1 RNA.
Computational modeling indicates that these synonymous mutations have a
(Malthusian) selection coefficient of the order of -0.002 and that they are
brought up to high frequency by hitchhiking on neighboring beneficial
nonsynonymous alleles. The patterns of fixation of nonsynonymous mutations
estimated from the longitudinal data and comparisons with computer models
suggest that escape mutations in C2-V5 are only transiently beneficial, either
because the immune system is catching up or because of competition between
equivalent escapes.


  A reduction in number and an increase in size of inflorescences is a common
aspect of plant domestication. When maize was domesticated from teosinte, the
number and arrangement of ears changed dramatically. Teosinte has long lateral
branches that bear multiple small ears at their nodes and tassels at their
tips. Maize has much shorter lateral branches that are tipped by a single large
ear with no additional ears at the branch nodes. To investigate the genetic
basis of this difference in prolificacy (the number of ears on a plant), we
performed a genome-wide QTL scan. A large effect QTL for prolificacy (prol1.1)
was detected on the short arm of chromosome one in a location that has
previously been shown to influence multiple domestication traits. We
fine-mapped prol1.1 to a 2.7 kb interval or causative region upstream of the
grassy tillers1 gene, which encodes a homeodomain leucine zipper transcription
factor. Tissue in situ hybridizations reveal that the maize allele of prol1.1
is associated with up-regulation of gt1 expression in the nodal plexus. Given
that maize does not initiate secondary ear buds, the expression of gt1 in the
nodal plexus in maize may suppress their initiation. Population genetic
analyses indicate positive selection on the maize allele of prol1.1, causing a
partial sweep that fixed the maize allele throughout most of domesticated
maize. This work shows how a subtle cis-regulatory change in tissue specific
gene expression altered plant architecture in a way that improved the
harvestability of maize.


  Rapid adaptation has been observed in numerous organisms in response to
selective pressures, such as the application of pesticides and the presence of
pathogens. When rapid adaptation is driven by rare alleles from the standing
genetic variation or by a high population rate of de novo adaptive mutation,
positive selection should commonly generate soft rather that hard selective
sweeps. In a soft sweep, multiple adaptive haplotypes sweep through the
population simultaneously, in contrast to hard sweeps in which only a single
adaptive haplotype rises to high frequency. Current statistical methods were
not designed to detect soft sweeps, and are therefore likely to miss these
possibly numerous adaptive events. Here, we develop a statistical test (H12)
based on haplotype homozygosity that is capable of detecting both hard and soft
sweeps with similar power. We use H12 to identify multiple genomic regions that
have undergone recent and strong adaptation in a population sample of fully
sequenced Drosophila melanogaster strains from the Drosophila Genetic Reference
Panel (DGRP). Visual inspection of the top 50 peaks revealed that multiple
haplotypes are at high frequency, consistent with signatures of soft sweep. We
developed a second statistic (H2/H1) that is sensitive to signatures common to
soft sweeps but not hard sweeps, in order to determine whether sweeps detected
by H12 can be more easily generated by hard versus soft sweeps. Surprisingly,
we find that the H12 and H2/H1 values for all top 50 peaks are more easily
generated by soft sweeps than hard sweeps under several evolutionary scenarios.


  Taking precautions before or during the start of a virus outbreak can heavily
reduce the number of infected. The question which individuals should be
immunized in order to mitigate the impact of the virus on the rest of
population has received quite some attention in the literature. The dynamics of
the of a virus spread through a population is often represented as information
spread over a complex network. The strategies commonly proposed to determine
which nodes are to be selected for immunization often involve only one
centrality measure at a time, while often the topology of the network seems to
suggest that a single metric is insufficient to capture the influence of a node
entirely.
  In this work we present a generic method based on a genetic algorithm (GA)
which does not rely explicitly on any centrality measures during its search but
only exploits this type of information to narrow the search space. The fitness
of an individual is defined as the estimated expected number of infections of a
virus following SIR dynamics. The proposed method is evaluated on two contact
networks: the Goodreau's Faux Mesa high school and the US air transportation
network. The GA method manages to outperform the most common strategies based
on a single metric for the air transportation network and its performance is
comparable with the best performing strategy for the high school network.


  Genome-wide binding preferences of the key components of eukaryotic
pre-initiation complex (PIC) have been recently measured with high resolution
in Saccharomyces cerevisiae by Rhee and Pugh (Nature (2012) 483:295-301). Yet
the rules determining the PIC binding specificity remain poorly understood. In
this study we show that nonconsensus protein-DNA binding significantly
influences PIC binding preferences. We estimate that such nonconsensus binding
contribute statistically at least 2-3 kcal/mol (on average) of additional
attractive free energy per protein, per core promoter region. The predicted
attractive effect is particularly strong at repeated poly(dA:dT) and
poly(dC:dG) tracts. Overall, the computed free energy landscape of nonconsensus
protein-DNA binding shows strong correlation with the measured genome-wide PIC
occupancy. Remarkably, statistical PIC binding preferences to both
TFIID-dominated and SAGA-dominated genes correlate with the nonconsensus free
energy landscape, yet these two groups of genes are distinguishable based on
the average free energy profiles. We suggest that the predicted nonconsensus
binding mechanism provides a genome-wide background for specific promoter
elements, such as transcription factor binding sites, TATA-like elements, and
specific binding of the PIC components to nucleosomes. We also show that
nonconsensus binding influences transcriptional frequency genome-wide.


  We study the influence of a THz field on thermal properties of DNA molecules.
A Peyrard- Bishop-Dauxois model with the inclusion of a solvent interaction
term is considered. The THz field is included as a sinusoidal driven force in
the equation of mo tion. We show how under certain field and system parameters,
melting transition and bubble formation are modified.


  We consider a population subdivided into two demes connected by migration in
which selection acts in opposite direction. We explore the effects of
recombination and migration on the maintenance of multilocus polymorphism, on
local adaptation, and on differentiation by employing a deterministic model
with genic selection on two linked diallelic loci (i.e., no dominance or
epistasis). For the following cases, we characterize explicitly the possible
equilibrium configurations: weak, strong, highly asymmetric, and
super-symmetric migration, no or weak recombination, and independent or
strongly recombining loci. For independent loci (linkage equilibrium) and for
completely linked loci, we derive the possible bifurcation patterns as
functions of the total migration rate, assuming all other parameters are fixed
but arbitrary. For these and other cases, we determine analytically the maximum
migration rate below which a stable fully polymorphic equilibrium exists. In
this case, differentiation and local adaptation are maintained. Their degree is
quantified by a new multilocus version of $\Fst$ and by the migration load,
respectively. In addition, we investigate the invasion conditions of locally
beneficial mutants and show that linkage to a locus that is already in
migration-selection balance facilitates invasion. Hence, loci of much smaller
effect can invade than predicted by one-locus theory if linkage is sufficiently
tight. We study how this minimum amount of linkage admitting invasion depends
on the migration pattern. This suggests the emergence of clusters of locally
beneficial mutations, which may form `genomic islands of divergence'. Finally,
the influence of linkage and two-way migration on the effective migration rate
at a linked neutral locus is explored. Numerical work complements our
analytical results.


  In this paper we present a new modelling framework combining replicator
dynamics (which is the standard model of frequency dependent selection) with
the model of an age-structured population. The new framework allows for the
modelling of populations consisting of competing strategies carried by
individuals who change across their life cycle. Firstly the discretization of
the McKendrick von Foerster model is derived. It is shown that the Euler--Lotka
equation is satisfied when the new model reaches a steady state (i.e. stable
frequencies between the age classes). This discretization consists of the unit
age classes and the timescale is chosen that only a fraction of individuals
play single game round. This implies linear dynamics within single time unit
when individuals not killed during game round are moved from one age class to
another. Since its local linear behaviour the system is equivalent to large
Bernadelli-Lewis-Leslie matrix. Then the methodology of multipopulation games
is used for the derivation of two, mutually equivalent systems of equations.
The first contains equations describing the evolution of the strategy
frequencies in the whole population completed by subsystems of equations
describing the evolution of the age structure for each strategy. The second
system contains equations describing the changes of the general population's
age structure, completed with subsystems of equations describing the selection
of the strategies within each age class. Then the obtained system of replicator
dynamics is presented in the form of the mixed ODE-PDE system. The obtained
results are illustrated by example of the sex ratio model which shows that when
different mortalities of both sexes are assumed, the sex ratio of 0.5 is
obtained but that Fisher's mechanism driven by the reproductive value of the
different sexes is not in equilibrium.


  Secular trends of the incidence of NCDs are especially important as they
indicate changes of the risk profile of a population. The article describes a
method for detecting secular trends in the incidence from a series of
prevalence data - without requiring costly follow-up studies or running a
register. After describing the theory, the method is applied to the incidence
of diabetes in Denmark.


  We present a study of the DNA translocation of the bacteriophage \phi 29
packaging molecular motor. From the experimental available information we
present a model system based in an stochastic fashing potential, which
reproduces the experimental observations such as: detailed trajectories, steps
and substeps, spatial correlation, and velocity. Moreover the model allows the
evaluation of power and efficiency of this motor. We have found that the
maximum power regime does not correspond with that of the maximum efficiency.
These informations can stimulate further experiments.


  Stereoisomers have the same molecular formula and the same atom connectivity
and their existence can be related to the presence of different
three-dimensional arrangements. Stereoisomerism is of great importance in many
different fields since the molecular properties and biological effects of the
stereoisomers are often significantly different. Most drugs for example, are
often composed of a single stereoisomer of a compound, and while one of them
may have therapeutic effects on the body, another may be toxic. A challenging
task is the automatic detection of stereoisomers using line input
specifications such as SMILES or InChI since it requires information about
group theory (to distinguish stereoisomers using mathematical information about
its symmetry), topology and geometry of the molecule. There are several
software packages that include modules to handle stereochemistry, especially
the ones to name a chemical structure and/or view, edit and generate chemical
structure diagrams. However, there is a lack of software capable of
automatically analyzing a molecule represented as a graph and generate a
classification of the type of isomerism present in a given atom or bond.
Considering the importance of stereoisomerism when comparing chemical
structures, this report describes a computer program for analyzing and
processing steric information contained in a chemical structure represented as
a molecular graph and providing as output a binary classification of the isomer
type based on the recommended conventions. Due to the complexity of the
underlying issue, specification of stereochemical information is currently
limited to explicit stereochemistry and to the two most common types of
stereochemistry caused by asymmetry around carbon atoms: chiral atom and double
bond. A Webtool to automatically identify and classify stereochemistry is
available at http://nams.lasige.di.fc.ul.pt/tools.php


  High-confidence prediction of complex traits such as disease risk or drug
response is an ultimate goal of personalized medicine. Although genome-wide
association studies have discovered thousands of well-replicated polymorphisms
associated with a broad spectrum of complex traits, the combined predictive
power of these associations for any given trait is generally too low to be of
clinical relevance. We propose a novel systems approach to complex trait
prediction, which leverages and integrates similarity in genetic,
transcriptomic or other omics-level data. We translate the omic similarity into
phenotypic similarity using a method called Kriging, commonly used in
geostatistics and machine learning. Our method called OmicKriging emphasizes
the use of a wide variety of systems-level data, such as those increasingly
made available by comprehensive surveys of the genome, transcriptome and
epigenome, for complex trait prediction. Furthermore, our OmicKriging framework
allows easy integration of prior information on the function of subsets of
omics-level data from heterogeneous sources without the sometimes heavy
computational burden of Bayesian approaches. Using seven disease datasets from
the Wellcome Trust Case Control Consortium (WTCCC), we show that OmicKriging
allows simple integration of sparse and highly polygenic components yielding
comparable performance at a fraction of the computing time of a recently
published Bayesian sparse linear mixed model method. Using a cellular growth
phenotype, we show that integrating mRNA and microRNA expression data
substantially increases performance over either dataset alone. We also
integrate genotype and expression data to predict change in LDL cholesterol
levels after statin treatment and show that OmicKriging performs better than
the polygenic score method. We provide an R package to implement OmicKriging.


  The entropy rates of the Wright-Fisher process, the Moran process, and
generalizations are computed and used to compare these processes and their
dependence on standard evolutionary parameters. Entropy rates are measures of
the variation dependent on both short-run and long-run behavior, and allow the
relationships between mutation, selection, and population size to be examined.
Bounds for the entropy rate are given for the Moran process (independent of
population size) and for the Wright-Fisher process (bounded for fixed
population size). A generational Moran process is also presented for comparison
to the Wright-Fisher Process. Results include analytic results and
computational extensions.


  GC-biased gene conversion (gBGC) is a recombination-associated process that
favors the fixation of G/C alleles over A/T alleles. In mammals, gBGC is
hypothesized to contribute to variation in GC content, rapidly evolving
sequences, and the fixation of deleterious mutations, but its prevalence and
general functional consequences remain poorly understood. gBGC is difficult to
incorporate into models of molecular evolution and so far has primarily been
studied using summary statistics from genomic comparisons. Here, we introduce a
new probabilistic model that captures the joint effects of natural selection
and gBGC on nucleotide substitution patterns, while allowing for correlations
along the genome in these effects. We implemented our model in a computer
program, called phastBias, that can accurately detect gBGC tracts ~1 kilobase
or longer in simulated sequence alignments. When applied to real primate genome
sequences, phastBias predicts gBGC tracts that cover roughly 0.3% of the human
and chimpanzee genomes and account for 1.2% of human-chimpanzee nucleotide
differences. These tracts fall in clusters, particularly in subtelomeric
regions; they are enriched for recombination hotspots and fast-evolving
sequences; and they display an ongoing fixation preference for G and C alleles.
We also find some evidence that they contribute to the fixation of deleterious
alleles, including an enrichment for disease-associated polymorphisms. These
tracts provide a unique window into historical recombination processes along
the human and chimpanzee lineages; they supply additional evidence of long-term
conservation of megabase-scale recombination rates accompanied by rapid
turnover of hotspots. Together, these findings shed new light on the
evolutionary, functional, and disease implications of gBGC. The phastBias
program and our predicted tracts are freely available.


  We present and study a minimal structure-based model for the self-assembly of
peptides into ordered beta-sheet-rich fibrils. The peptides are represented by
unit-length sticks on a cubic lattice and interact by hydrogen bonding and
hydrophobicity forces. By Monte Carlo simulations with >100,000 peptides, we
show that fibril formation occurs with sigmoidal kinetics in the model. To
determine the mechanism of fibril nucleation, we compute the joint distribution
in length and width of the aggregates at equilibrium, using an efficient
cluster move and flat-histogram techniques. This analysis, based on simulations
with 256 peptides in which aggregates form and dissolve reversibly, shows that
the main free-energy barriers that a nascent fibril has to overcome are
associated with changes in width.


We introduce design and optimization considerations for the 'khmer' package.


  The evolution of leadership in migratory populations depends not only on
costs and benefits of leadership investments but also on the opportunities for
individuals to rely on cues from others through social interactions. We derive
an analytically tractable adaptive dynamic network model of collective
migration with fast timescale migration dynamics and slow timescale adaptive
dynamics of individual leadership investment and social interaction. For large
populations, our analysis of bifurcations with respect to investment cost
explains the observed hysteretic effect associated with recovery of migration
in fragmented environments. Further, we show a minimum connectivity threshold
above which there is evolutionary branching into leader and follower
populations. For small populations, we show how the topology of the underlying
social interaction network influences the emergence and location of leaders in
the adaptive system. Our model and analysis can describe other adaptive network
dynamics involving collective tracking or collective learning of a noisy,
unknown signal, and likewise can inform the design of robotic networks where
agents use decentralized strategies that balance direct environmental
measurements with agent interactions.


  We present a data structure called a history graph that offers a practical
basis for the analysis of genome evolution. It conceptually simplifies the
study of parsimonious evolutionary histories by representing both substitutions
and double cut and join (DCJ) rearrangements in the presence of duplications.
The problem of constructing parsimonious history graphs thus subsumes related
maximum parsimony problems in the fields of phylogenetic reconstruction and
genome rearrangement. We show that tractable functions can be used to define
upper and lower bounds on the minimum number of substitutions and DCJ
rearrangements needed to explain any history graph. These bounds become tight
for a special type of unambiguous history graph called an ancestral variation
graph (AVG), which constrains in its combinatorial structure the number of
operations required. We finally demonstrate that for a given history graph $G$,
a finite set of AVGs describe all parsimonious interpretations of $G$, and this
set can be explored with a few sampling moves.


  How to compare whole genome sequences at large scale has not been achieved
via conventional methods based on pair-wisely base-to-base comparison;
nevertheless, no attention was paid to handle in-one-sitting a number of
genomes crossing genetic category (chromosome, plasmid, and phage) with farther
divergences (much less or no homologous) over large size ranges (from Kbp to
Mbp). We created a new method, GenomeFingerprinter, to unambiguously produce
three-dimensional coordinates from a sequence, followed by one
three-dimensional plot and six two-dimensional trajectory projections to
illustrate whole genome fingerprints. We further developed a set of concepts
and tools and thereby established a new method, universal genome fingerprint
analysis. We demonstrated their applications through case studies on over a
hundred of genome sequences. Particularly, we defined the total genetic
component configuration (TGCC) (i.e., chromosome, plasmid, and phage) for
describing a strain as a system, and the universal genome fingerprint map
(UGFM) of TGCC for differentiating a strain as a universal system, as well as
the systematic comparative genomics (SCG) for comparing in-one-sitting a number
of genomes crossing genetic category in diverse strains. By using UGFM,
UGFM-TGCC, and UGFM-TGCC-SCG, we compared a number of genome sequences with
farther divergences (chromosome, plasmid, and phage; bacterium, archaeal
bacterium, and virus) over large size ranges (6Kbp~5Mbp), giving new insights
into critical problematic issues in microbial genomics in the post-genomic era.
This paper provided a new method for rapidly computing, geometrically
visualizing, and intuitively comparing genome sequences at fingerprint level,
and hence established a new method of universal genome fingerprint analysis for
systematic comparative genomics.


  Cancer cells maintain a high level of aerobic glycolysis (the Warburg
effect), which is associated with their rapid proliferation. Many studies have
reported that the suppression of glycolysis and activation of oxidative
phosphorylation can repress the growth of cancer cells through regulation of
key regulators. Whether Warburg effect of cancer cells could be switched by
some other environmental stimulus? Herein, we report an interesting phenomenon
in which cells alternated between glycolysis and mitochondrial respiration
depending on the type of radiation they were exposed to. We observed enhanced
glycolysis and mitochondrial respiration in HeLa cells exposed to 2-Gy X-ray
and 2-Gy carbon ion radiation, respectively. This discovery may provide novel
insights for tumor therapy.


  The physiology and genetics underlying juvenility is poorly understood. Here
we exploit Arabidopsis as a system to understand the mechanisms that regulate
floral incompetence during juvenility. Using an experimental assay that allows
the length of juvenility to be estimated, and mutants impaired in different
pathways, we show that multiple inputs influence juvenility. Juvenile phase
lengths of wild type (WT) accessions Col-0, Ler-0 and Ws-4 are shown to differ,
with Col-0 having the shortest and Ws-4 the longest length. Plants defective in
sugar signalling [gin1-1, gin2-1, gin6 (abi4)] and floral repressor mutants
[hst1, tfl1, tfl2 (lhp1)] showed shortened juvenile phase lengths, compared to
their respective WTs. Mutants defective in starch anabolism (adg1-1, pgm1) and
catabolism (sex1, sex4, bam3) showed prolonged juvenile phase lengths compared
to Col-0. Examination of diurnal metabolite changes in adg1-1 and sex1 mutants
indicates that their altered juvenile phase length may be due to lack of starch
turnover, which influences carbohydrate availability. In this article we
propose a model in which a variety of signals including floral activators and
repressors modulate the juvenile-to-adult phase transition. The role of
carbohydrates may be in their capacity as nutrients, osmotic regulators,
signalling molecules and/ or through their interaction with phytohormonal
networks.


  Whole transcriptome sequencing is increasingly being used as a functional
genomics tool to study non- model organisms. However, when the reference
transcriptome used to calculate differential expression is incomplete,
significant error in the inferred expression levels can result. In this study,
we use simulated reads generated from real transcriptomes to determine the
accuracy of read mapping, and measure the error resulting from using an
incomplete transcriptome. We show that the two primary sources of count- ing
error are 1) alternative splice variants that share reads and 2) missing
transcripts from the reference. Alternative splice variants increase the false
positive rate of mapping while incomplete reference tran- scriptomes decrease
the true positive rate, leading to inaccurate transcript expression levels.
Grouping transcripts by gene or read sharing (similar to mapping to a reference
genome) significantly decreases false positives, but only by improving the
reference transcriptome itself can the missing transcript problem be addressed.
We also demonstrate that employing different mapping software does not yield
substantial increases in accuracy on simulated data. Finally, we show that read
lengths or insert sizes must increase past 1kb to resolve mapping ambiguity.


  The paradox of pesticides was observed experimentally, which says that
pesticides may dramatically increase the population of a pest when the pest has
a natural predator. Here we use a mathematical model to study the paradox. We
find that the timing for the application of pesticides is crucial for the
resurgence or non-resurgence of the pests. In particular, regularly applying
pesticides is not a good idea as also observed in experiments. In fact, the
best time to apply pesticides is when the pest population is reasonably high.


  Cybersecurity attacks are a major and increasing burden to economic and
social systems globally. Here we analyze the principles of security in
different domains and demonstrate an architectural flaw in current
cybersecurity. Cybersecurity is inherently weak because it is missing the
ability to defend the overall system instead of individual computers. The
current architecture enables all nodes in the computer network to communicate
transparently with one another, so security would require protecting every
computer in the network from all possible attacks. In contrast, other systems
depend on system-wide protections. In providing conventional security, police
patrol neighborhoods and the military secures borders, rather than defending
each individual household. Likewise, in biology, the immune system provides
security against viruses and bacteria using primarily action at the skin,
membranes, and blood, rather than requiring each cell to defend itself. We
propose applying these same principles to address the cybersecurity challenge.
This will require: (a) Enabling pervasive distribution of self-propagating
securityware and creating a developer community for such securityware, and (b)
Modifying the protocols of internet routers to accommodate adaptive security
software that would regulate internet traffic. The analysis of the immune
system architecture provides many other principles that should be applied to
cybersecurity. Among these principles is a careful interplay of detection and
action that includes evolutionary improvement. However, achieving significant
security gains by applying these principles depends strongly on remedying the
underlying architectural limitations.


  Stochasticity plays important roles in molecular networks when molecular
concentrations are in the range of $0.1 \mu$M to $10 n$M (about 100 to 10
copies in a cell). The chemical master equation provides a fundamental
framework for studying these networks, and the time-varying landscape
probability distribution over the full microstates provide a full
characterization of the network dynamics. A complete characterization of the
space of the microstates is a prerequisite for obtaining the full landscape
probability distribution of a network. However, there are neither closed-form
solutions nor algorithms fully describing all microstates for a given molecular
network.
  We have developed an algorithm that can exhaustively enumerate the
microstates of a molecular network of small copy numbers under the finite
buffer condition that the net gain in newly synthesized molecules is smaller
than a predefined limit. We also describe a simple method for computing the
exact mean or steady state landscape probability distribution over microstates.
We show how the full landscape probability for the gene networks of the
self-regulating gene and the toggle-switch in the steady state can be fully
characterized. We also give an example using the MAPK cascade network.
  Our algorithm works for networks of small copy numbers buffered with a finite
copy number of net molecules that can be synthesized, regardless of the
reaction stoichiometry, and is optimal in both storage and time complexity. The
buffer size is limited by the available memory or disk storage. Our algorithm
is applicable to a class of biological networks when the copy numbers of
molecules are small and the network is closed, or the network is open but the
net gain in newly synthesized molecules does not exceed a predefined buffer
capacity.


  The local mechanical properties of DNA are believed to play an important role
in their biological functions and DNA-based nanomechanical devices. Using a
simple sphere-tip compression system, the local radial mechanical properties of
DNA are systematically studied by changing the tip size. The compression
simulation results for the 16 nm diameter sphere tip are well consistent with
the experimental results. With the diameter of the tip decreasing, the radial
compressive elastic properties under external loads become sensitive to the tip
size and the local DNA conformation. There appears a suddenly force break in
the compression-force curve when the sphere size is less than or equal to 12 nm
diameter. The analysis of the hydrogen bonds and base stacking interaction
shows there is a local unwinding process occurs. During the local unwinding
process, first the hydrogen bonds between complement base pairs are broken.
With the compression aggregating, the local backbones in the compression center
are unwound from the double helix conformation to a kind of parallel
conformation. This local unwinding behavior deducing by external loads is
helpful to understand the biological process, and important to DNA-based
nanomechanical devices.


  The ability to represent intracellular biochemical dynamics via deterministic
and stochastic modelling is one of the crucial components to move biological
sciences in the observe-predict-control-design knowledge ladder. Compared to
the engineering or physics problems, dynamical models in quantitative biology
typically dependent on a relatively large number of parameters. Therefore, the
relationship between model parameters and dynamics is often prohibitively
difficult to determine. We developed a method to depict the input-output
relationship for multi-parametric stochastic and deterministic models via
information-theoretic quantification of similarity between model parameters and
modules. Identification of most information-theoretically orthogonal biological
components, provided mathematical language to precisely communicate and
visualise compensation like phenomena such as biological robustness, sloppiness
and statistical non-identifiability. A comprehensive analysis of the
multi-parameter NF-$\kappa$B signalling pathway demonstrates that the
information-theoretic similarity reflects a topological structure of the
network. Examination of the currently available experimental data on this
system reveals the number of identifiable parameters and suggests informative
experimental protocols.


  In order to better understand the interplay of partnership and competition in
population dynamics we study a family of generalized May-Leonard models with
$N$ species. These models have a very rich structure, characterized by
different types of space-time patterns. Interesting partnership formations
emerge following the maxim that "the enemy of my enemy is my friend." In
specific cases cyclic dominance within coarsening clusters yields a peculiar
coarsening behavior with intriguing pattern formation. We classify the
different types of dynamics through the analysis of the square of the adjacency
matrix. The dependence of the population densities on emerging pattern and
propagating wave fronts is elucidated through a Fourier analysis. Finally,
after having identified collaborating teams, we study interface fluctuations
where we initially populate different parts of the system with different teams.


  In this paper, we consider the problem of optimal exogenous control of gene
regulatory networks. Our approach consists in adapting an established
reinforcement learning algorithm called the fitted Q iteration. This algorithm
infers the control law directly from the measurements of the system's response
to external control inputs without the use of a mathematical model of the
system. The measurement data set can either be collected from wet-lab
experiments or artificially created by computer simulations of dynamical models
of the system. The algorithm is applicable to a wide range of biological
systems due to its ability to deal with nonlinear and stochastic system
dynamics. To illustrate the application of the algorithm to a gene regulatory
network, the regulation of the toggle switch system is considered. The control
objective of this problem is to drive the concentrations of two specific
proteins to a target region in the state space.


  The biological significance of genomic features is often context-dependent.
We present CruzDB, a fast and intuitive programmatic interface to the UCSC
genome browser that facilitates integrative analyses of diverse local and
remotely hosted datasets. We showcase the syntax of CruzDB using miRNA-binding
sites as examples, and further demonstrate its utility with 3 novel biological
discoveries. First, we find that while exons replicate early, introns tend to
replicate late, suggesting a complex replication pattern in gene regions.
Second, variants associated with cognitive functions map to lincRNA transcripts
of relevant function. Third, lamina-associated domains are highly enriched in
olfaction-related genes. CruzDB is available at
https://github.com/brentp/cruzdb


  While the thermodynamics of DNA hybridization is well understood, much less
is known about the kinetics of this classic system. Filling this gap in our
understanding has new urgency because DNA nanotechnology often depends
critically on binding rates. Here we use a coarse-grained model to explore the
hybridization kinetics of DNA oligomers, finding that strand association
proceeds through a complex set of intermediate states. Successful binding
events start with the formation of a few metastable base-pairing interactions,
followed by zippering of the remaining bonds. However, despite reasonably
strong interstrand interactions, initial contacts frequently fail to lead to
zippering because the typical configurations in which they form differ from
typical states of similar enthalpy in the double-stranded equilibrium ensemble.
Therefore, if the association process is analyzed on the base-pair (secondary
structure) level, it shows non-Markovian behavior. Initial contacts must be
stabilized by two or three base pairs before full zippering is likely,
resulting in negative effective activation enthalpies. Non-Arrhenius behavior
is observed as the number of base pairs in the effective transition state
increases with temperature. In addition, we find that alternative pathways
involving misbonds can increase association rates. For repetitive sequences,
misaligned duplexes frequently rearrange to form fully paired duplexes by two
distinct processes which we label `pseudoknot' and `inchworm' internal
displacement. We show how the above processes can explain why experimentally
observed association rates of GC-rich oligomers are higher than rates of
AT-rich equivalents. More generally, we argue that association rates can be
modulated by sequence choice.


  By monitoring multiple molecular transitions, force-clamp and
trap-position-clamp methods have led to precise determinations of the free
energies and free energy landscapes for molecular states populated in
equilibrium at the same or similar forces. Here, we present a powerful new
elaboration of the force-clamp and force-jump methods, applicable to
transitions far from equilibrium. Specifically, we have implemented a live jump
detection and force-clamp algorithm that intelligently adjusts and maintains
the force on a single molecule in response to the measured state of that
molecule. We are able to collect hundreds of individual molecular transitions
at different forces, many times faster than previously, permitting us to
accurately determine force-dependent lifetime distributions and reaction rates.
Application of our method to unwinding and rewinding the nucleosome inner turn,
using optical tweezers reveals experimental lifetime distributions that
comprise a statistically-meaningful number of transitions, and that are
accurately single exponential. These measurements significantly reduce the
error in the previously measured rates, and demonstrate the existence of a
single, dominant free energy barrier at each force studied. A key benefit of
the molecular yoyo method for nucleosomes is that it reduces as far as possible
the time spent in the tangentially-bound state, which minimizes the loss of
nucleosomes by dissociation.


  In this work we analyze the evolution of voluntary vaccination in networked
populations by entangling the spreading dynamics of an influenza-like disease
with an evolutionary framework taking place at the end of each influenza season
so that individuals take or not the vaccine upon their previous experience. Our
framework thus put in competition two well-known dynamical properties of
scale-free networks: the fast propagation of diseases and the promotion of
cooperative behaviors. Our results show that when vaccine is perfect scale-free
networks enhance the vaccination behavior with respect to random graphs with
homogeneous connectivity patterns. However, when imperfection appears we find a
cross-over effect so that the number of infected (vaccinated) individuals
increases (decreases) with respect to homogeneous networks, thus showing up the
competition between the aforementioned properties of scale-free graphs.


  The tremdendous advances in high-throughput sequencing technologies have made
population-scale sequencing as performed in the 1000 Genomes project and the
Genome of the Netherlands project possible. Next-generation sequencing has
allowed genom-wide discovery of variations beyond single-nucleotide
polymorphisms (SNPs), in particular of structural variations (SVs) like
deletions, insertions, duplications, translocations, inversions, and even more
complex rearrangements. Here, we design a read aligner with special emphasis on
the following properties: (1) high sensitivity, i.e. find all (reasonable)
alignments; (2) ability to find (long) indels; (3) statistically sound
alignment scores; and (4) runtime fast enough to be applied to whole genome
data. We compare performance to BWA, bowtie2, stampy and find that our methods
is especially advantageous on reads containing larger indels.


  Modelling the Calvin cycle of photosynthesis leads to various systems of
ordinary differential equations and reaction-diffusion equations. They differ
by the choice of chemical substances included in the model, the choices of
stoichiometric coefficients and chemical kinetics and whether or not diffusion
is taken into account. This paper studies the long-time behaviour of solutions
of several of these systems, concentrating on the ODE case. In some examples it
is shown that there exist two positive stationary solutions. In several cases
it is shown that there exist solutions where the concentrations of all
substrates tend to zero at late times and others (runaway solutions) where the
concentrations of all substrates increase without limit. In another case, where
the concentration of ATP is explicitly included, runaway solutions are ruled
out.


  The evolution of drug resistance in HIV occurs by the fixation of specific,
well-known, drug-resistance mutations, but the underlying population genetic
processes are not well understood. By analyzing within-patient longitudinal
sequence data, we make four observations that shed a light on the underlying
processes and allow us infer the short-term effective population size of the
viral population in a patient. Our first observation is that the evolution of
drug resistance usually occurs by the fixation of one drug-resistance mutation
at a time, as opposed to several changes simultaneously. Second, we find that
these fixation events are accompanied by a reduction in genetic diversity in
the region surrounding the fixed drug-resistance mutation, due to the
hitchhiking effect. Third, we observe that the fixation of drug-resistance
mutations involves both hard and soft selective sweeps. In a hard sweep, a
resistance mutation arises in a single viral particle and drives all linked
mutations with it when it spreads in the viral population, which dramatically
reduces genetic diversity. On the other hand, in a soft sweep, a resistance
mutation occurs multiple times on different genetic backgrounds, and the
reduction of diversity is weak. Using the frequency of occurrence of hard and
soft sweeps we estimate the effective population size of HIV to be 1.5 * 10^5
(95% confidence interval 0.8 * 10^5, 4.8 * 10^5). This number is much lower
than the actual number of infected cells, but much larger than previous
population size estimates based on synonymous diversity. We propose several
explanations for the observed discrepancies. Finally, our fourth observation is
that genetic diversity at non-synonymous sites recovers to its pre-fixation
value within 18 months, whereas diversity at synonymous sites remains depressed
after this time period.


  Traditional epidemic models consider that individual processes occur at
constant rates. That is, an infected individual has a constant probability per
unit time of recovering from infection after contagion. This assumption
certainly fails for almost all infectious diseases, in which the infection time
usually follows a probability distribution more or less spread around a mean
value. We show a general treatment for an SIRS model in which both the infected
and the immune phases admit such a description. The general behavior of the
system shows transitions between endemic and oscillating situations that could
be relevant in many real scenarios. The interaction with the other main source
of oscillations, seasonality, is also discussed.


  Fitness landscapes are central in analyzing evolution, in particular for drug
resistance mutations for bacteria and virus. We show that the fitness
landscapes associated with antibiotic resistance are not compatible with any of
the classical models; additive, uncorrelated and block fitness landscapes. The
NK model is also discussed. It is frequently stated that virtually nothing is
known about fitness landscapes in nature. We demonstrate that available records
of antimicrobial drug mutations can reveal interesting properties of fitness
landscapes in general. We apply the methods to analyze the TEM family of
$\beta$-lactamases associated with antibiotic resistance. Laboratory results
agree with our observations. The qualitative tools we suggest are well suited
for comparisons of empirical fitness landscapes. Fitness landscapes are central
in the theory of recombination and there is a potential for finding relations
between the tools and recombination strategies.


  Characterizing the spatial extent of epidemics at the outbreak stage is key
to controlling the evolution of the disease. At the outbreak, the number of
infected individuals is typically small, so that fluctuations around their
average are important: then, it is commonly assumed that the
susceptible-infected-recovered (SIR) mechanism can be described by a stochastic
birth-death process of Galton-Watson type. The displacements of the infected
individuals can be modelled by resorting to Brownian motion, which is
applicable when long-range movements and complex network interactions can be
safely neglected, as in case of animal epidemics. In this context, the spatial
extent of an epidemic can be assessed by computing the convex hull enclosing
the infected individuals at a given time. We derive the exact evolution
equations for the mean perimeter and the mean area of the convex hull, and
compare them with Monte Carlo simulations.


  We propose a new framework for the study of continuous time dynamical systems
on networks. We view such dynamical systems as collections of interacting
control systems. We show that a class of maps between graphs called graph
fibrations give rise to maps between dynamical systems on networks. This allows
us to produce conjugacy between dynamical systems out of combinatorial data. In
particular we show that surjective graph fibrations lead to synchrony subspaces
in networks. The injective graph fibrations, on the other hand, give rise to
surjective maps from large dynamical systems to smaller ones. One can view
these surjections as a kind of "fast/slow" variable decompositions or as
"abstractions" in the computer science sense of the word.


  Demographic stochasticity is important in determining extinction risks of
small populations, but it is largely unknown how its effect depends on the life
histories of species. We modeled effects of demographic stochasticity on
extinction risk in a broad range of generalized life histories, using matrix
models and branching processes. Extinction risks of life histories varied
greatly in their sensitivity to demographic stochasticity. Comparing life
histories, extinction risk generally increased with increasing fecundity and
decreased with higher ages of maturation. Effects of adult survival depended on
age of maturation. At lower ages of maturation, extinction risk peaked at
intermediate levels of adult survival, but it increased along with adult
survival at higher ages of maturation. These differences were largely explained
by differences in sensitivities of population growth to perturbations of
life-history traits. Juvenile survival rate contributed most to total
demographic variance in the majority of life histories. Our general results
confirmed earlier findings, suggesting that empirical patterns can be explained
by a relatively simple model. Thus, basic life history information can be used
to assign life-history-specific sensitivity to demographic stochasticity. This
is of great value when assessing the vulnerability of small populations.


  Systems biology uses large networks of biochemical reactions to model the
functioning of biological cells from the molecular to the cellular scale. The
dynamics of dissipative reaction networks with many well separated time scales
can be described as a sequence of successive equilibrations of different
subsets of variables of the system. Polynomial systems with separation are
equilibrated when at least two monomials, of opposite signs, have the same
order of magnitude and dominate the others. These equilibrations and the
corresponding truncated dynamics, obtained by eliminating the dominated terms,
find a natural formulation in tropical analysis and can be used for model
reduction.


  Summary: BWA-MEM is a new alignment algorithm for aligning sequence reads or
long query sequences against a large reference genome such as human. It
automatically chooses between local and end-to-end alignments, supports
paired-end reads and performs chimeric alignment. The algorithm is robust to
sequencing errors and applicable to a wide range of sequence lengths from 70bp
to a few megabases. For mapping 100bp sequences, BWA-MEM shows better
performance than several state-of-art read aligners to date.
  Availability and implementation: BWA-MEM is implemented as a component of
BWA, which is available at http://github.com/lh3/bwa.
  Contact: hengli@broadinstitute.org


  Okasha (2006) proposed distinguishing aspects of selection: those based in
particle level traits (MSL1), and those based in group level traits (MSL2). It
is proposed here that MSL1 can usefully be further split into two aspects, one
(MLS1E) representing selection of particles based in their individual
interaction with environmental properties, and one (MLS1G) representing Multi
Level Selection of particles based in their relation to group properties.
Similarly MSL2 can be split into two parts based in this distinction. This
splitting enables a characterisation of how emergent group properties can
affect particle selection, and thus affect group traits that are important for
survival. This proposal is illustrated by considering a key aspect of animal
and human life, namely the formation of social groups, which greatly enhances
survival prospects. The biological mechanism that underlies such group
formation is the existence of innate primary emotional systems studied by
Panksepp (1998), effective through the ascending systems in the brain that
diffuse neurotransmitters such as dopamine and norepinephrine through the
cortex. Evolutionary emergence of such brain mechanisms, and hence the
emergence of social groups, can only result from multi-level selection
characterised by the combination of MSL2 and MLS1G. The distinctions proposed
here should be useful in other contexts.