  The evolution of aligned DNA sequence sites is generally modeled by a Markov
process operating along the edges of a phylogenetic tree. It is well known that
the probability distribution on the site patterns at the tips of the tree
determines the tree and its branch lengths. However, the number of patterns is
typically much larger than the number of edges, suggesting considerable
redundancy in the branch length estimation. In this paper we ask whether the
probabilities of just the `edge-specific' patterns (the ones that correspond to
a change of state on a single edge) suffice to recover the branch lengths of
the tree, under a symmetric 2-state Markov process. We first show that this
holds provided the branch lengths are sufficiently short, by applying the
inverse function theorem. We then consider whether this restriction to short
branch lengths is necessary, and show that for trees with up to four leaves it
can be lifted. This leaves open the interesting question of whether this holds
in general.


  Large-scale DNA deformation is ubiquitous in transcriptional regulation in
prokaryotes and eukaryotes alike. Though much is known about how transcription
factors and constellations of binding sites dictate where and how gene
regulation will occur, less is known about the role played by the intervening
DNA. In this work we explore the effect of sequence flexibility on
transcription factor-mediated DNA looping, by drawing on sequences identified
in nucleosome formation and ligase-mediated cyclization assays as being
especially favorable for or resistant to large deformations. We examine a
poly(dA:dT)-rich, nucleosome-repelling sequence that is often thought to belong
to a class of highly inflexible DNAs; two strong nucleosome positioning
sequences that share a set of particular sequence features common to
nucleosome-preferring DNAs; and a CG-rich sequence representative of high
G+C-content genomic regions that correlate with high nucleosome occupancy in
vivo. To measure the flexibility of these sequences in the context of DNA
looping, we combine the in vitro single-molecule tethered particle motion
assay, a canonical looping protein, and a statistical mechan- ical model that
allows us to quantitatively relate the looping probability to the looping free
energy. We show that, in contrast to the case of nucleosome occupancy, G+C
content does not positively correlate with looping probability, and that
despite sharing sequence features that are thought to determine nucleosome
affinity, the two strong nucleosome positioning sequences behave markedly
dissimilarly in the context of looping. Most surprisingly, the poly(dA:dT)-rich
DNA that is often characterized as highly inflexible in fact exhibits one of
the highest propensities for looping that we have measured.


  It is becoming routine to obtain datasets on DNA sequence variation across
several thousands of chromosomes, providing unprecedented opportunity to infer
the underlying biological and demographic forces. Such data make it vital to
study summary statistics which offer enough compression to be tractable, while
preserving a great deal of information. One well-studied summary is the site
frequency spectrum---the empirical distribution, across segregating sites, of
the sample frequency of the derived allele. However, most previous theoretical
work has assumed that each site has experienced at most one mutation event in
its genealogical history, which becomes less tenable for very large sample
sizes. In this work we obtain, in closed-form, the predicted frequency spectrum
of a site that has experienced at most two mutation events, under very general
assumptions about the distribution of branch lengths in the underlying
coalescent tree. Among other applications, we obtain the frequency spectrum of
a triallelic site in a model of historically varying population size. We
demonstrate the utility of our formulas in two settings: First, we show that
triallelic sites are more sensitive to the parameters of a population that has
experienced historical growth, suggesting that they will have use if they can
be incorporated into demographic inference. Second, we investigate a recently
proposed alternative mechanism of mutation in which the two derived alleles of
a triallelic site are created simultaneously within a single individual, and we
develop a test to determine whether it is responsible for the excess of
triallelic sites in the human genome.


  Shared genealogies introduce allele dependencies in diploid genotypes, as
alleles within an individual or between different individuals will likely match
when they originate from a recent common ancestor. At a locus shared by a pair
of diploid individuals, there are nine combinatorially distinct modes of
identity-by-descent (IBD), capturing all possible combinations of coancestry
and inbreeding. A distribution over the IBD modes is described by the nine
associated probabilities, known as (Jacquard's) identity coefficients. The
genetic relatedness between two individuals can be succinctly characterized by
the identity coefficients corresponding to the joint genealogy. The identity
coefficients (together with allele frequencies) determine the distribution of
joint genotypes at a locus. At a locus with two possible alleles, identity
coefficients are not identifiable because different coefficients can generate
the same genotype distribution.
  We analyze precisely how different IBD modes combine into identical genotype
distributions at diallelic loci. In particular, we describe IBD mode mixtures
that result in identical genotype distributions at all allele frequencies,
implying the non-identifiability of the identity coefficients from independent
loci. Our analysis yields an exhaustive characterization of relatedness
statistics that are always identifiable. Importantly, we show that identifiable
relatedness statistics include the kinship coefficient (probability that a
random pair of alleles are identical by descent between individuals) and
inbreeding-related measures, which can thus be estimated from genotype
distributions at independent loci.


  Cancer forms a robust system and progresses as stages over time typically
with increasing aggressiveness and worsening prognosis. Characterizing these
stages and identifying the genes driving transitions between them is critical
to understand cancer progression and to develop effective anti-cancer
therapies. Here, we propose a novel model of the 'cancer system' as a Boolean
state space in which a Boolean network, built from protein interaction and
gene-expression data from different stages of cancer, transits between Boolean
satisfiability states by "editing" interactions and "flipping" genes. The
application of our model (called BoolSpace) on three case studies - pancreatic
and breast tumours in human and post spinal-cord injury in rats - reveals
valuable insights into the phenomenon of cancer progression. In particular, we
notice that several of the genes flipped are serine/threonine kinases which act
as natural cellular switches and that different sets of genes are flipped
during the initial and final stages indicating a pattern to tumour progression.
We hypothesize that robustness of cancer partly stems from "passing of the
baton" between genes at different stages, and therefore an effective therapy
should target a "cover set" of these genes. A C/C++ implementation of BoolSpace
is freely available at: http://www.bioinformatics.org.au/tools-data


  Thermodynamics constrains the flow of matter in a reaction network to occur
through routes along which the Gibbs energy decreases, implying that viable
steady-state flux patterns should be void of closed reaction cycles.
Identifying and removing cycles in large reaction networks can unfortunately be
a highly challenging task from a computational viewpoint. We propose here a
method that accomplishes it by combining a relaxation algorithm and a Monte
Carlo procedure to detect loops, with ad hoc rules (discussed in detail) to
eliminate them. As test cases, we tackle (a) the problem of identifying
infeasible cycles in the E. coli metabolic network and (b) the problem of
correcting thermodynamic infeasibilities in the Flux-Balance-Analysis solutions
for 15 human cell-type specific metabolic networks. Results for (a) are
compared with previous analyses of the same issue, while results for (b) are
weighed against alternative methods to retrieve thermodynamically viable flux
patterns based on minimizing specific global quantities. Our method on one hand
outperforms previous techniques and, on the other, corrects loopy solutions to
Flux Balance Analysis. As a byproduct, it also turns out to be able to reveal
possible inconsistencies in model reconstructions.


  It has recently been shown that structural conditions on the reaction
network, rather than a 'fine-tuning' of system parameters, often suffice to
impart 'absolute concentration robustness' on a wide class of biologically
relevant, deterministically modeled mass-action systems [Shinar and Feinberg,
Science, 2010]. We show here that fundamentally different conclusions about the
long-term behavior of such systems are reached if the systems are instead
modeled with stochastic dynamics and a discrete state space. Specifically, we
characterize a large class of models that exhibit convergence to a positive
robust equilibrium in the deterministic setting, whereas trajectories of the
corresponding stochastic models are necessarily absorbed by a set of states
that reside on the boundary of the state space, i.e. the system undergoes an
extinction event. If the time to extinction is large relative to the relevant
time-scales of the system, the process will appear to settle down to a
stationary distribution long before the inevitable extinction will occur. This
quasi-stationary distribution is considered for two systems taken from the
literature, and results consistent with absolute concentration robustness are
recovered by showing that the quasi-stationary distribution of the robust
species approaches a Poisson distribution.


  Demographic change of human populations is one of the central questions for
delving into the past of human beings. To identify major population expansions
related to male lineages, we sequenced 78 East Asian Y chromosomes at 3.9 Mbp
of the non-recombining region (NRY), discovered >4,000 new SNPs, and identified
many new clades. The relative divergence dates can be estimated much more
precisely using molecular clock. We found that all the Paleolithic divergences
were binary; however, three strong star-like Neolithic expansions at ~6 kya
(thousand years ago) (assuming a constant substitution rate of 1e-9/bp/year)
indicates that ~40% of modern Chinese are patrilineal descendants of only three
super-grandfathers at that time. This observation suggests that the main
patrilineal expansion in China occurred in the Neolithic Era and might be
related to the development of agriculture.


  We examine the foraging behavior of stingless bees (Tetragonula biroi Friese)
with regard to their preference of food location. The three factors considered
are distance of the food source from the hive, direction of the food source,
and height of the food source from the ground. Single-factor experiments for
each factor were performed to determine the preferred food location.
Statistical analysis shows that T. biroi species prefer the feeders located 1m
away from the hive, which is the shortest distance used in the experiment.
Results also show that this group of bees has no conclusive preferred
direction, which might be a result of non-apparent cause, such as random
search. In addition, the bees prefer the feeders located 1m above the ground,
which are directly in front of the hive opening. Using the results from the
single-factor experiments, a Three-Factor Factorial Experimental Design is used
to analyze the interaction of the three factors. Factorial Analysis of Variance
shows that the interaction among the three factors affects the preference of
the bees. We hypothesize that the preference of the bees can be explained by
(1) optimal foraging theory, (2) memory of food location, and (3) chemical
marks left by conspecifics.


  In the Susceptible-Infectious-Recovered (SIR) model of disease spreading, the
time to extinction of the epidemics happens at an intermediate value of the
per-contact transmission probability. Too contagious infections burn out fast
in the population. Infections that are not contagious enough die out before
they spread to a large fraction of people. We characterize how the maximal
extinction time in SIR simulations on networks depend on the network structure.
For example we find that the average distances in isolated components, weighted
by the component size, is a good predictor of the maximal time to extinction.
Furthermore, the transmission probability giving the longest outbreaks is
larger than, but otherwise seemingly independent of, the epidemic threshold.


  The denaturation of the double helix is a template for fundamental biological
functions such as replication and transcription involving the formation of
local fluctuational openings. The denaturation transition is studied for
heterogeneous short sequences of DNA, i.e. $\sim 100$ base pairs, in the
framework of a mesoscopic Hamiltonian model which accounts for the helicoidal
geometry of the molecule. The theoretical background for the application of the
path integral formalism to predictive analysis of the molecule thermodynamical
properties is discussed. The base pair displacements with respect to the ground
state are treated as paths whose temperature dependent amplitudes are governed
by the thermal wavelength. The ensemble of base pairs paths is selected, at any
temperature, consistently with both the model potential and the second law of
thermodynamics. The partition function incorporates the effects of the base
pair thermal fluctuations which become stronger close to the denaturation. The
transition appears as a gradual phenomenon starting from the molecule segments
rich in adenine-thymine base pairs. Computing the equilibrium thermodynamics,
we focus on the interplay between twisting of the complementary strands around
the molecule axis and nonlinear stacking potential: it is shown that the latter
affects the melting profiles only if the rotational degrees of freedom are
included in the Hamiltonian. The use of ladder Hamiltonian models for the DNA
complementary strands in the pre-melting regime is questioned.


  G-protein coupled receptors (GPCRs), a major gatekeeper of extracellular
signals on plasma membrane, are unarguably one of the most important
therapeutic targets. Given the recent discoveries of allosteric modulations, an
allosteric wiring diagram of intramolecular signal transductions would be of
great use to glean the mechanism of receptor regulation. Here, by evaluating
betweenness centrality ($C_B$) of each residue, we calculate maps of
information flow in GPCRs and identify key residues for signal transductions
and their pathways. Compared with preexisting approaches, the allosteric
hotspots that our $C_B$-based analysis detects for A$_{2A}$ adenosine receptor
(A$_{2A}$AR) and bovine rhodopsin are better correlated with biochemical data.
In particular, our analysis outperforms other methods in locating the rotameric
microswitches, which are generally deemed critical for mediating orthosteric
signaling in class A GPCRs. For A$_{2A}$AR, the inter-residue cross-correlation
map, calculated using equilibrium structural ensemble from molecular dynamics
simulations, reveals that strong signals of long-range transmembrane
communications exist only in the agonist-bound state. A seemingly subtle
variation in structure, found in different GPCR subtypes or imparted by agonist
bindings or a point mutation at an allosteric site, can lead to a drastic
difference in the map of signaling pathways and protein activity. The signaling
map of GPCRs provides valuable insights into allosteric modulations as well as
reliable identifications of orthosteric signaling pathways.


  DNA Polymerase is an enzyme that participates in the process of DNA
replication. In general, forward DNA polymerase moves along the template strand
in a 3'-5' direction, and the daughter strand is formed in a 5'-3' direction.
This paper theoretically analyzes the mechanism of forward DNA polymerase, and
gives the design proposal of a new reverse DNA polymerase (3'-5' direction).


  Recent empirical research has shown that links between groups reinforce
individuals within groups to adopt cooperative behaviour. Moreover, links
between networks may induce cascading failures, competitive percolation, or
contribute to efficient transportation. Here we show that there in fact exists
an intermediate fraction of links between groups that is optimal for the
evolution of cooperation in the prisoner's dilemma game. We consider individual
groups with regular, random, and scale-free topology, and study their different
combinations to reveal that an intermediate interdependence optimally
facilitates the spreading of cooperative behaviour between groups. Excessive
between-group links simply unify the two groups and make them act as one, while
too rare between-group links preclude a useful information flow between the two
groups. Interestingly, we find that between-group links are more likely to
connect two cooperators than in-group links, thus supporting the conclusion
that they are of paramount importance.


  We consider the problem of exact learning of parameters of a linear RNA
energy model from secondary structure data. A necessary and sufficient
condition for learnability of parameters is derived, which is based on
computing the convex hull of union of translated Newton polytopes of input
sequences. The set of learned energy parameters is characterized as the convex
cone generated by the normal vectors to those facets of the resulting polytope
that are incident to the origin. In practice, the sufficient condition may not
be satisfied by the entire training data set; hence, computing a maximal subset
of training data for which the sufficient condition is satisfied is often
desired. We show that problem is NP-hard in general for an arbitrary
dimensional feature space. Using a randomized greedy algorithm, we select a
subset of RNA STRAND v2.0 database that satisfies the sufficient condition for
separate A-U, C-G, G-U base pair counting model. The set of learned energy
parameters includes experimentally measured energies of A-U, C-G, and G-U
pairs; hence, our parameter set is in agreement with the Turner parameters.


  In order to control the growth of human population it is helpful to
understand correctly the mechanism of growth, and the first essential step is
to investigate current interpretations and reject any unscientific
explanations. One of such popular but questionable interpretations is the
concept of the Epoch of Malthusian Stagnation. We discuss its origin, narrative
and claims. We explain why this concept is scientifically unacceptable. This
investigation questions also the closely-related Demographic Transition Theory,
whose essential component is the assumed mechanism of Malthusian stagnation for
the first stage of growth.


  As the development of a dengue vaccine is ongoing, we simulate an
hypothetical vaccine as an extra protection to the population. In a first
phase, the vaccination process is studied as a new compartment in the model,
and different ways of distributing the vaccines investigated: pediatric and
random mass vaccines, with distinct levels of efficacy and durability. In a
second step, the vaccination is seen as a control variable in the
epidemiological process. In both cases, epidemic and endemic scenarios are
included in order to analyze distinct outbreak realities.


  Many intracellular processes continue to oscillate during the cell cycle.
Although it is not well-understood how they are affected by discontinuities in
the cellular environment, the general assumption is that oscillations remain
robust provided the period of cell divisions is much larger than the period of
the oscillator. Here, I will show that under these conditions a cell will in
fact have to correct for an additional quantity added to the phase of
oscillation upon every repetition of the cell cycle. The resulting phase shift
is an analogue of the geometric phase, a curious entity first discovered in
quantum mechanics. In this Letter, I will discuss the theory of the geometric
phase shift and demonstrate its relevance to biological oscillations.


  A central problem when studying adaptation to a new environment is the
interplay between genetic variation and phenotypic plasticity. Arabidopsis
thaliana has colonized a wide range of habitats across the world and it is
therefore an attractive model for studying the genetic mechanisms underlying
environmental adaptation 4,5. Here, we used publicly available data from two
collections of A. thaliana accessions, covering the native range of the
species, to identify loci associated with differences in climates at the
sampling sites. To address the confounding between geographic location, climate
and population structure, a new genome-wide association analysis method was
developed that facilitates detection of potentially adaptive loci where the
alternative alleles display different tolerable climate ranges. Sixteen novel
such loci, many of which contained candidate genes with amino acid changes,
were found including a strong association between Chromomethylase 2 (CMT2) and
variability in seasonal temperatures. The reference allele dominated in areas
with less seasonal variability in temperature, and the alternative allele,
which disrupts genome-wide CHH-methylation, existed in both stable and variable
regions. Our results link natural variation in CMT2, and differential
genome-wide CHH methylation, to the distribution of A. thaliana accessions
across habitats with different seasonal temperature variability. They also
suggest a role for genetic regulation of epigenetic modifications in natural
adaptation, potentially through differential allelic plasticity, and illustrate
the importance of re-analyses of existing data using new analytical methods to
obtain a more complete understanding of the mechanisms contributing to
adaptation.


  In previous work, RAF theory has been developed as a tool for making
theoretical progress on the origin of life question, providing insight into the
structure and occurrence of self-sustaining and collectively autocatalytic sets
within catalytic polymer networks. We present here an extension in which there
are two "independent" polymer sets, where catalysis occurs within and between
the sets, but there are no reactions combining polymers from both sets. Such an
extension reflects the interaction between nucleic acids and peptides observed
in modern cells and proposed forms of early life.


  Unusually high bending flexibility has been recently reported for DNA on
short length scales. We use atomic force microscopy (AFM) in solution to obtain
a direct estimate of DNA bending statistics for scales down to one helical
turn. It appears that DNA behaves as a Gaussian chain and is well described by
the worm-like chain model at length scales beyond 3 helical turns (10.5nm).
Below this threshold, the AFM data exhibit growing noise because of
experimental limitations. This noise may hide small deviations from the
Gaussian behavior, but they can hardly be significant.


  Networks play a prominent role in the study of complex systems of interacting
entities in biology, sociology, and economics. Despite this diversity, we
demonstrate here that a statistical model decomposing networks into matching
and centrality components provides a comprehensive and unifying quantification
of their architecture. First we show, for a diverse set of networks, that this
decomposition provides an extremely tight fit to observed networks.
Consequently, the model allows very accurate prediction of missing links in
partially known networks. Second, when node characteristics are known, we show
how the matching-centrality decomposition can be related to this external
information. Consequently, it offers a simple and versatile tool to explore how
node characteristics explain network architecture. Finally, we demonstrate the
efficiency and flexibility of the model to forecast the links that a novel node
would create if it were to join an existing network.


  Analysis of two independent data sets with increased taxonomic resolution
(genera rather than families) using the revised 2012 time scale reveals that an
extinction periodicity first detected by Raup and Sepkoski (1984) for only the
post-Paleozoic actually runs through the entire Phanerozoic. Although there is
not a local peak of extinction every 27 million years, an excess of the
fraction of genus extinction by interval follows a 27 million year timing
interval and differs from a random distribution at the p ~ 0.02 level. A 27
million year periodicity in the spectrum of interval lengths no longer appears,
removing the question of a possible artifact arising from it. Using a method
originally developed in Bambach (2006) we identify 19 intervals of marked
extinction intensity, including mass extinctions, spanning the last 470 million
years (and with another six present in the Cambrian) and find that 10 of the 19
lie within 3 Myr of the maxima in the spacing of the 27 Myr periodicity, which
differs from a random distribution at the p = 0.004 level. These 19 intervals
of marked extinction intensity also preferentially occur during decreasing
diversity phases of a well-known 62 Myr periodicity in diversity (16 of 19, p =
0.002). Both periodicities appear to enhance the likelihood of increased
severity of extinction, but the cause of neither periodicity is known.
Variations in the strength of the many suggested causes of extinction coupled
to the variation in combined effect of the two different periodicities as they
shift in and out of phase is surely one of the reasons that definitive
comparative study of the causes of major extinction events is so elusive.


  One important problem in genome science is to determine sets of co-regulated
genes based on measurements of gene expression levels across samples, where the
quantification of expression levels includes substantial technical and
biological noise. To address this problem, we developed a Bayesian sparse
latent factor model that uses a three parameter beta prior to flexibly model
shrinkage in the loading matrix. By applying three layers of shrinkage to the
loading matrix (global, factor-specific, and element-wise), this model has
non-parametric properties in that it estimates the appropriate number of
factors from the data. We added a two-component mixture to model each factor
loading as being generated from either a sparse or a dense mixture component;
this allows dense factors that capture confounding noise, and sparse factors
that capture local gene interactions. We developed two statistics to quantify
the stability of the recovered matrices for both sparse and dense matrices. We
tested our model on simulated data and found that we successfully recovered the
true latent structure as compared to related models. We applied our model to a
large gene expression study and found that we recovered known covariates and
small groups of co-regulated genes. We validated these gene subsets by testing
for associations between genotype data and these latent factors, and we found a
substantial number of biologically important genetic regulators for the
recovered gene subsets.


  The potential of synthetic biology techniques for designing complex cellular
circuits able to solve complicated computations opens a whole domain of
exploration, beyond experiments and theory. Such cellular circuits could be
used to carry out hard tasks involving decision-making, storage of information,
or signal processing. Since Gene Regulatory Networks (GRNs) are the best known
technical approach to synthetic designs, it would be desirable to know in
advance the potential of such circuits in performing tasks and how classical
approximations dealing with neural networks can be translated into GRNs. In
this paper such a potential is analyzed. Here we show that feed-forward GRNs
are capable of performing classic machine intelligence tasks. Therefore, two
important milestones in the success of Artificial Neural Networks are reached
for models of GRNs based on Hill equations, namely the back-propagation
algorithm and the proof that GRNs can approximate arbitrary positive functions.
Potential extensions and implications for synthetic designs are outlined.


  Clustering is typically measured by the ratio of triangles to all triples,
open or closed. Generating clustered networks, and how clustering affects
dynamics on networks, is reasonably well understood for certain classes of
networks \cite{vmclust, karrerclust2010}, e.g., networks composed of lines and
non-overlapping triangles. In this paper we show that it is possible to
generate networks which, despite having the same degree distribution and equal
clustering, exhibit different higher-order structure, specifically, overlapping
triangles and other order-four (a closed network motif composed of four nodes)
structures. To distinguish and quantify these additional structural features,
we develop a new network metric capable of measuring order-four structure
which, when used alongside traditional network metrics, allows us to more
accurately describe a network's topology. Three network generation algorithms
are considered: a modified configuration model and two rewiring algorithms. By
generating homogeneous networks with equal clustering we study and quantify
their structural differences, and using SIS (Susceptible-Infected-Susceptible)
and SIR (Susceptible-Infected-Recovered) dynamics we investigate
computationally how differences in higher-order structure impact on epidemic
threshold, final epidemic or prevalence levels and time evolution of epidemics.
Our results suggest that characterising and measuring higher-order network
structure is needed to advance our understanding of the impact of network
topology on dynamics unfolding on the networks.


  The influence of migration on the stochastic dynamics of subdivided
populations is still an open issue in various evolutionary models. We develop
here a self-consistent mean-field-like method in order to determine the effects
of migration on relevant nonequilibrium properties, such as the mean fixation
time. If evolution strongly favors coexistence of species (e.g., balancing
selection), the mean fixation time develops an unexpected minimum as a function
of the migration rate. Our analysis hinges only on the presence of a separation
of time scales between local and global dynamics and therefore it carries over
to other nonequilibrium processes in physics, biology, ecology, and social
sciences.


  Different models such as diffusion-collision and nucleation-condensation have
been used to unravel how secondary and tertiary structures form during protein
folding. However, a simple mechanism based on physical principles that provide
an accurate description of kinetics and thermodynamics for such phenomena has
not yet been identified. This study introduces the hypothesis that the
synchronization of the peptide plane oscillatory movements throughout the
backbone must also play a key role in the folding mechanism. Based on that, we
draw a parallel between the folding process and the dynamics for a network of
coupled oscillators described by the Kuramoto model. The amino acid coupling
may explain the mean-field character of the force that propels an amino acid
sequence into a structure through self-organization. Thus, the pattern of
synchronized cluster formation and growing helps to solve the Levinthal's
paradox.Synchronization may also help us to understand the success of homology
structural modeling, allosteric effect, and the mechanism responsible for the
recognition of odorants by olfactory receptors.


  Economic experiments reveal that humans value cooperation and fairness.
Punishing unfair behavior is therefore common, and according to the theory of
strong reciprocity, it is also directly related to rewarding cooperative
behavior. However, empirical data fail to confirm that positive and negative
reciprocity are correlated. Inspired by this disagreement, we determine whether
the combined application of reward and punishment is evolutionary advantageous.
We study a spatial public goods game, where in addition to the three elementary
strategies of defection, rewarding and punishment, a fourth strategy combining
the later two competes for space. We find rich dynamical behavior that gives
rise to intricate phase diagrams where continuous and discontinuous phase
transitions occur in succession. Indirect territorial competition, spontaneous
emergence of cyclic dominance, as well as divergent fluctuations of
oscillations that terminate in an absorbing phase are observed. Yet despite the
high complexity of solutions, the combined strategy can survive only in very
narrow and unrealistic parameter regions. Elementary strategies, either in pure
or mixed phases, are much more common and likely to prevail. Our results
highlight the importance of patterns and structure in human cooperation, which
should be considered in future experiments.


  One goal of human genetics is to understand how the information for precise
and dynamic gene expression programs is encoded in the genome. The interactions
of transcription factors (TFs) with DNA regulatory elements clearly play an
important role in determining gene expression outputs, yet the regulatory logic
underlying functional transcription factor binding is poorly understood. Many
studies have focused on characterizing the genomic locations of TF binding, yet
it is unclear to what extent TF binding at any specific locus has functional
consequences with respect to gene expression output. To evaluate the context of
functional TF binding we knocked down 59 TFs and chromatin modifiers in one
HapMap lymphoblastoid cell line. We then identified genes whose expression was
affected by the knockdowns. We intersected the gene expression data with
transcription factor binding data (based on ChIP-seq and DNase-seq) within 10
kb of the transcription start sites of expressed genes. This combination of
data allowed us to infer functional TF binding. On average, 14.7% of genes
bound by a factor were differentially expressed following the knockdown of that
factor, suggesting that most interactions between TF and chromatin do not
result in measurable changes in gene expression levels of putative target
genes. We found that functional TF binding is enriched in regulatory elements
that harbor a large number of TF binding sites, at sites with predicted higher
binding affinity, and at sites that are enriched in genomic regions annotated
as active enhancers.


  In this paper we introduce and study three new measures for efficient
discriminative comparison of phylogenetic trees. The NNI navigation
dissimilarity $d_{nav}$ counts the steps along a "combing" of the Nearest
Neighbor Interchange (NNI) graph of binary hierarchies, providing an efficient
approximation to the (NP-hard) NNI distance in terms of "edit length". At the
same time, a closed form formula for $d_{nav}$ presents it as a weighted count
of pairwise incompatibilities between clusters, lending it the character of an
edge dissimilarity measure as well. A relaxation of this formula to a simple
count yields another measure on all trees --- the crossing dissimilarity
$d_{CM}$. Both dissimilarities are symmetric and positive definite (vanish only
between identical trees) on binary hierarchies but they fail to satisfy the
triangle inequality. Nevertheless, both are bounded below by the widely used
Robinson-Foulds metric and bounded above by a closely related true metric, the
cluster-cardinality metric $d_{CC}$. We show that each of the three proposed
new dissimilarities is computable in time $O(n^2)$ in the number of leaves $n$,
and conclude the paper with a brief numerical exploration of the distribution
over tree space of these dissimilarities in comparison with the Robinson-Foulds
metric and the more recently introduced matching-split distance.


  The introduction and persistence of novel sexually antagonistic alleles can
depend upon factors that differ between males and females. Understanding the
conditions for invasion in a two-locus model can elucidate these processes. For
instance, selection can act differently upon the sexes, or sex-linkage can
facilitate the invasion of genetic variation with opposing fitness effects
between the sexes. Two factors that deserve further attention are recombination
rates and allele frequencies -- both of which can vary substantially between
the sexes. We find that sex-specific recombination rates in a two-locus diploid
model can affect the invasion outcome of sexually antagonistic alleles and that
the sex-averaged recombination rate is not necessarily sufficient to predict
invasion. We confirm that the range of permissible recombination rates is
smaller in the sex benefitting from invasion and larger in the sex harmed by
invasion. However, within the invasion space, male recombination rate can be
greater than, equal to, or less than female recombination rate in order for a
male-benefit, female-detriment allele to invade (and similarly for a
female-benefit, male-detriment allele). We further show that a novel, sexually
antagonistic allele that is also associated with a lowered recombination rate
can invade more easily when present in the double heterozygote genotype.
Finally, we find that sexual dimorphism in resident allele frequencies can
impact the invasion of new sexually antagonistic alleles at a second locus. Our
results suggest that accounting for sex-specific recombination rates and allele
frequencies can determine the difference between invasion and non-invasion of
novel sexually antagonistic alleles in a two-locus model.


  Empirical evidence questions fundamental concepts of the human population
dynamics. One of the key conclusions of this study is that positive checks
activate the efficient Malthusian regeneration mechanism, suggesting that the
Epoch of Malthusian Stagnation, the first stage of growth claimed by the
Demographic Transition Theory, did not exist.


  Short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) are
two kinds of commonly used markers in Y chromosome studies of forensic and
population genetics. There has been increasing interest in the cost saving
strategy by using the STR haplotypes to predict SNP haplogroups. However, the
convergence of Y chromosome STR haplotypes from different haplogroups might
compromise the accuracy of haplogroup prediction. Here, we compared the
worldwide Y chromosome lineages at both haplogroup level and haplotype level to
search for the possible haplotype similarities among haplogroups. The similar
haplotypes between haplogroups B and I2, C1 and E1b1b1, C2 and E1b1a1, H1 and
J, L and O3a2c1, O1a and N, O3a1c and O3a2b, and M1 and O3a2 have been found,
and those similarities reduce the accuracy of prediction.


  Sayyid Ajjal is the ancestor of many Muslims in areas all across China. And
one of his descendants is the famous Navigator of Ming Dynasty, Zheng He, who
led the largest armada in the world of 15th century. The origin of Sayyid
Ajjal's family remains unclear although many studies have been done on this
topic of Muslim history. In this paper, we studied the Y chromosomes of his
present descendants, and found they all have haplogroup L1a-M76, proving a
southern Persian origin.


  Properly coordinating cooperation is relevant for resolving public good
problems such as clean energy and environmental protection. However, little is
known about how individuals can coordinate themselves for a certain level of
cooperation in large populations of strangers. In a typical situation, a
consensus-building process hardly succeeds due to lack of face and standing.
The evolution of cooperation in this type of situation is studied using
threshold public good games in which cooperation prevails when it is initially
sufficient, or otherwise, it perishes. While punishment is a powerful tool to
shape human behaviours, institutional punishment is often too costly to start
with only a few contributors, which is another coordination problem. Here we
show that whatever the initial conditions, reward funds based on voluntary
contribution can evolve. The voluntary reward paves the way for effectively
overcoming the coordination problem and efficiently transforms freeloaders to
cooperators with a perceived small risk of collective failure.


  The rapid worldwide spread of severe viral infections, often involving novel
modifications of viruses, poses major challenges to our health care systems.
This means that tools that can efficiently and specifically diagnose viruses
are much needed. To be relevant for a broad application in local health care
centers, such tools should be relatively cheap and easy to use. Here we discuss
the biophysical potential for the macroscopic detection of viruses based on the
induction of a mechanical stress in a bundle of pre-stretched DNA molecules
upon binding of viruses to the DNA. We show that the affinity of the DNA to the
charged virus surface induces a local melting of the double-helix into two
single-stranded DNA. This process effects a mechanical stress along the DNA
chains leading to an overall contraction of the DNA. Our results suggest that
when such DNA bundles are incorporated in a supporting matrix such as a
responsive hydrogel, the presence of viruses may indeed lead to a significant,
macroscopic mechanical deformation of the matrix. We discuss the biophysical
basis for this effect and characterize the physical properties of the
associated DNA melting transition. In particular, we reveal several scaling
relations between the relevant physical parameters of the system. We promote
this DNA-based assay for efficient and specific virus screening.


  We analyze the solutions, on single network instances, of a recently
introduced class of constraint-satisfaction problems (CSPs), describing
feasible steady states of chemical reaction networks. First, we show that the
CSPs generalize the scheme known as Network Expansion, which is recovered in a
specific limit. Next, a full statistical mechanics characterization (including
the phase diagram and a discussion of physical origin of the phase transitions)
for Network Expansion is obtained. Finally, we provide a message-passing
algorithm to solve the original CSPs in the most general form.


  We consider a birth and death process in which death is due to both `natural
death' and to competition between individuals, modelled as a quadratic function
of population size. The resulting `logistic branching process' has been
proposed as a model for numbers of individuals in populations competing for
some resource, or for numbers of species. However, because of the quadratic
death rate, even if the intrinsic growth rate is positive, the population will,
with probability one, die out in finite time. There is considerable interest in
understanding the process conditioned on non-extinction.
  In this paper, we exploit a connection with the ancestral selection graph of
population genetics to find expressions for the transition rates in the
logistic branching process conditioned on survival until some fixed time $T$,
in terms of the distribution of a certain one-dimensional diffusion process at
time $T$. We also find the probability generating function of the Yaglom
distribution of the process and rather explicit expressions for the transition
rates for the so-called Q-process, that is the logistic branching process
conditioned to stay alive into the indefinite future. For this process, one can
write down the joint generator of the (time-reversed) total population size and
what in population genetics would be called the `genealogy' and in
phylogenetics would be called the `reconstructed tree' of a sample from the
population.
  We explore some ramifications of these calculations numerically.


  RNA protein interactions control the fate of cellular RNAs and play an
important role in gene regulation. An interdependency between such interactions
allows for the implementation of logic functions in gene regulation. We
investigate the interplay between RNA binding partners in the context of the
statistical physics of RNA secondary structure, and define a linear correlation
function between the two partners as a measurement of the interdependency of
their binding events. We demonstrate the emergence of a long-range power-law
behavior of this linear correlation function. This suggests RNA secondary
structure driven interdependency between binding sites as a general mechanism
for combinatorial post-transcriptional gene regulation.


  A recent study (Fish. Bull. 109:394-401 (2011)) purportedly tests two
hypotheses: 1. that the capture of elasmobranchs would be reduced with hooks
containing magnets in comparison with control hooks in hook-and-line and
longline studies. 2. that the presence of permanent magnets on hooks would not
alter teleost capture because teleosts lack the ampullary organ. Review of this
paper shows some inconsistencies in the data supporting the first hypothesis
and insufficient data and poor experimental design to adequately test the
second hypothesis. Further, since several orders of teleosts are known to
possess ampullary organs and demonstrate electroreception, grouping all
teleosts in a study design or data analysis of magnetic hook catch rates is not
warranted. Adequate tests of the hypothesis that permanent magnets or
magnetized hooks do not alter teleost capture requires a more careful study
design and much larger sample sizes than O'Connell et al. (Fish. Bull.
109:394-401 (2011)).


  The paternally inherited Y chromosome has been widely used in population
genetic studies to understand relationships among human populations. Our
interpretation of Y chromosomal evidence about population history and genetics
has rested on the assumption that all the Y chromosomal markers in the
male-specific region (MSY) are selectively neutral. However, the very low
diversity of Y chromosome has drawn a long debate about whether natural
selection has affected this chromosome or not. In recent several years, the
progress in Y chromosome sequencing has helped to address this dispute.
Purifying selection has been detected in the X-degenerate genes of human Y
chromosomes and positive selection might also have an influence in the
evolution of testis-related genes in the ampliconic regions. Those new findings
remind us to take the effect of natural selection into account when we use Y
chromosome in population genetic studies.


  Animal grouping behaviors have been widely studied due to their implications
for understanding social intelligence, collective cognition, and potential
applications in engineering, artificial intelligence, and robotics. An
important biological aspect of these studies is discerning which selection
pressures favor the evolution of grouping behavior. In the past decade,
researchers have begun using evolutionary computation to study the evolutionary
effects of these selection pressures in predator-prey models. The selfish herd
hypothesis states that concentrated groups arise because prey selfishly attempt
to place their conspecifics between themselves and the predator, thus causing
an endless cycle of movement toward the center of the group. Using an
evolutionary model of a predator-prey system, we show that how predators attack
is critical to the evolution of the selfish herd. Following this discovery, we
show that density-dependent predation provides an abstraction of Hamilton's
original formulation of ``domains of danger.'' Finally, we verify that
density-dependent predation provides a sufficient selective advantage for prey
to evolve the selfish herd in response to predation by coevolving predators.
Thus, our work corroborates Hamilton's selfish herd hypothesis in a digital
evolutionary model, refines the assumptions of the selfish herd hypothesis, and
generalizes the domain of danger concept to density-dependent predation.


  The existence of complex (multiple-step) genetic adaptations that are
"irreducible" (i.e., all partial combinations are less fit than the original
genotype) is one of the longest standing problems in evolutionary biology. In
standard genetics parlance, these adaptations require the crossing of a wide
adaptive valley of deleterious intermediate stages. Here we demonstrate, using
a simple model, that evolution can cross wide valleys to produce "irreducibly
complex" adaptations by making use of previously cryptic mutations. When
revealed by an evolutionary capacitor, previously cryptic mutants have higher
initial frequencies than do new mutations, bringing them closer to a
valley-crossing saddle in allele frequency space. Moreover, simple
combinatorics imply an enormous number of candidate combinations exist within
available cryptic genetic variation. We model the dynamics of crossing of a
wide adaptive valley after a capacitance event using both numerical simulations
and analytical approximations. Although individual valley crossing events
become less likely as valleys widen, by taking the combinatorics of genotype
space into account, we see that revealing cryptic variation can cause the
frequent evolution of complex adaptations. This finding also effectively
dismantles "irreducible complexity" as an argument against evolution by
providing a general mechanism for crossing wide adaptive valleys.


  How the neutral diversity is affected by selection and adaptation is
investigated in an eco-evolutionary framework. In our model, we study a finite
population in continuous time, where each individual is characterized by a
trait under selection and a completely linked neutral marker. Population
dynamics are driven by births and deaths, mutations at birth, and competition
between individuals. Trait values influence ecological processes (demographic
events, competition), and competition generates selection on trait variation,
thus closing the eco-evolutionary feedback loop. The demographic effects of the
trait are also expected to influence the generation and maintenance of neutral
variation. We consider a large population limit with rare mutation, under the
assumption that the neutral marker mutates faster than the trait under
selection. We prove the convergence of the stochastic individual-based process
to a new measure-valued diffusive process with jumps that we call Substitution
Fleming-Viot Process (SFVP). When restricted to the trait space this process is
the Trait Substitution Sequence first introduced by Metz et al. (1996). During
the invasion of a favorable mutation, a genetical bottleneck occurs and the
marker associated with this favorable mutant is hitchhiked. By rigorously
analysing the hitchhiking effect and how the neutral diversity is restored
afterwards, we obtain the condition for a time-scale separation; under this
condition, we show that the marker distribution is approximated by a
Fleming-Viot distribution between two trait substitutions. We discuss the
implications of the SFVP for our understanding of the dynamics of neutral
variation under eco-evolutionary feedbacks and illustrate the main phenomena
with simulations. Our results highlight the joint importance of mutations,
ecological parameters, and trait values in the restoration of neutral diversity
after a selective sweep.


  We propose a new approach to model composition, based on reducing several
models to the same level of complexity and subsequent combining them together.
Firstly, we suggest a set of model reduction tools that can be systematically
applied to a given model. Secondly, we suggest a notion of a minimal complexity
model. This model is the simplest one that can be obtained from the original
model using these tools and still able to approximate experimental data.
Thirdly, we propose a strategy for composing the reduced models together.
Connection with the detailed model is preserved, which can be advantageous in
some applications. A toolbox for model reduction and composition has been
implemented as part of the BioUML software and tested on the example of
integrating two previously published models of the CD95 (APO-1/Fas) signaling
pathways. We show that the reduced models lead to the same dynamical behavior
of observable species and the same predictions as in the precursor models. The
composite model is able to recapitulate several experimental datasets which
were used by the authors of the original models to calibrate them separately,
but also has new dynamical properties.


  Market trade-routes can support infectious-disease transmission, impacting
biological populations and even disrupting causal trade. Epidemiological models
increasingly account for reductions in infectious contact, such as
risk-aversion behaviour in response to pathogen outbreaks. However, market
dynamics clearly differ from simple risk-aversion, as are driven by different
motivation and conditioned by trade constraints, known in economics as
friction, that arise because exchanges are costly. Here we develop a novel
economic-market model where transient and long-term market dynamics are
determined by trade friction and agent adaptation, and can influence disease
transmission. We specify the participants, frequency, volume, and price in
trade transactions, and investigate, using analytical insights and simulation,
how trade friction affects joint market and epidemiological dynamics. The
friction values explored encompass estimates from French cattle and pig
markets. We show that, when trade is the dominant route of transmission, market
friction can be a significantly stronger determinant of epidemics than
risk-aversion behaviour. In particular, there is a critical friction level
above which epidemics do not occur. For a given level of friction, open
unregulated markets can boost epidemics compared with closed or tightly
regulated markets. Our results are robust to model specificities and can hold
in the presence of non-trade disease-transmission routes. In particular, we try
to explain why outbreaks in French livestock markets appear more frequently in
cattle than swine despite swine trade-flow being larger. To minimize contagion
in markets, safety policies could generate incentives for larger-volume,
less-frequent transactions, increasing trade friction without necessarily
affecting overall trade flow.


  Risk aversion is a common behavior universal to humans and animals alike.
Economists have traditionally defined risk preferences by the curvature of the
utility function. Psychologists and behavioral economists also make use of
concepts such as loss aversion and probability weighting to model risk
aversion. Neurophysiological evidence suggests that loss aversion has its
origins in relatively ancient neural circuitries (e.g., ventral striatum).
Could there thus be an evolutionary origin to risk avoidance? We study this
question by evolving strategies that adapt to play the equivalent mean payoff
gamble. We hypothesize that risk aversion in the equivalent mean payoff gamble
is beneficial as an adaptation to living in small groups, and find that a
preference for risk averse strategies only evolves in small populations of less
than 1,000 individuals, while agents exhibit no such strategy preference in
larger populations. Further, we discover that risk aversion can also evolve in
larger populations, but only when the population is segmented into small groups
of around 150 individuals. Finally, we observe that risk aversion only evolves
when the gamble is a rare event that has a large impact on the individual's
fitness. These findings align with earlier reports that humans lived in small
groups for a large portion of their evolutionary history. As such, we suggest
that rare, high-risk, high-payoff events such as mating and mate competition
could have driven the evolution of risk averse behavior in humans living in
small groups.


  Acquiring genomes at single-cell resolution has many applications such as in
the study of microbiota. However, deep sequencing and assembly of all of
millions of cells in a sample is prohibitively costly. A property that can come
to rescue is that deep sequencing of every cell should not be necessary to
capture all distinct genomes, as the majority of cells are biological
replicates. Biologically important samples are often sparse in that sense. In
this paper, we propose an adaptive compressed method, also known as distilled
sensing, to capture all distinct genomes in a sparse microbial community with
reduced sequencing effort. As opposed to group testing in which the number of
distinct events is often constant and sparsity is equivalent to rarity of an
event, sparsity in our case means scarcity of distinct events in comparison to
the data size. Previously, we introduced the problem and proposed a distilled
sensing solution based on the breadth first search strategy. We simulated the
whole process which constrained our ability to study the behavior of the
algorithm for the entire ensemble due to its computational intensity. In this
paper, we modify our previous breadth first search strategy and introduce the
depth first search strategy. Instead of simulating the entire process, which is
intractable for a large number of experiments, we provide a dynamic programming
algorithm to analyze the behavior of the method for the entire ensemble. The
ensemble analysis algorithm recursively calculates the probability of capturing
every distinct genome and also the expected total sequenced nucleotides for a
given population profile. Our results suggest that the expected total sequenced
nucleotides grows proportional to $\log$ of the number of cells and
proportional linearly with the number of distinct genomes.


  We investigate the potential of numerical algorithms to decipher the kinetic
parameters involved in multi-step chemical reactions. To this end we study a
dimerization kinetics of protein as a model system. We follow the dimerization
kinetics using a stochastic simulation algorithm and combine it with three
different optimization techniques (Genetic Algorithm, Simulated Annealing and
Parallel Tempering) to obtain the rate constants involved in each reaction
step. We find good convergence of the numerical scheme to the rate constants of
the process. We also perform a sensitivity test on the reaction kinetic
parameters to see the relative effects of the parameters for the associated
profile of the monomer/dimer distribution.


  The energetics of cerebral activity critically relies on the functional and
metabolic interactions between neurons and astrocytes. Important open questions
include the relation between neuronal versus astrocytic energy demand, glucose
uptake and intercellular lactate transfer, as well as their dependence on the
level of activity. We have developed a large-scale, constraint-based network
model of the metabolic partnership between astrocytes and glutamatergic neurons
that allows for a quantitative appraisal of the extent to which stoichiometry
alone drives the energetics of the system. We find that the velocity of the
glutamate-glutamine cycle ($V_{cyc}$) explains part of the uncoupling between
glucose and oxygen utilization at increasing $V_{cyc}$ levels. Thus, we are
able to characterize different activation states in terms of the tissue
oxygen-glucose index (OGI). Calculations show that glucose is taken up and
metabolized according to cellular energy requirements, and that partitioning of
the sugar between different cell types is not significantly affected by
$V_{cyc}$. Furthermore, both the direction and magnitude of the lactate shuttle
between neurons and astrocytes turn out to depend on the relative cell glucose
uptake while being roughly independent of $V_{cyc}$. These findings suggest
that, in absence of ad hoc activity-related constraints on neuronal and
astrocytic metabolism, the glutamate-glutamine cycle does not control the
relative energy demand of neurons and astrocytes, and hence their glucose
uptake and lactate exchange.


  Current Y chromosome research is limited in the poor resolution of Y
chromosome phylogenetic tree. Entirely sequenced Y chromosomes in numerous
human individuals have only recently become available by the advent of
next-generation sequencing technology. The 1000 Genomes Project has sequenced Y
chromosomes from more than 1000 males. Here, we analyzed 1000 Genomes Project Y
chromosome data of 1269 individuals and discovered about 25,000 phylogenetic
relevant SNPs. Those new markers are useful in the phylogeny of Y chromosome
and will lead to an increased phylogenetic resolution for many Y chromosome
studies.


  As proposed by Alan Turing in 1952 as a ubiquitous mechanism for
nonequilibrium pattern formation, diffusional effects may destabilize uniform
distributions of reacting chemical species and lead to both spatially and
temporally heterogeneous patterns. While stationary Turing patterns are broadly
known, the oscillatory instability, leading to traveling waves in continuous
media and also called the wave bifurcation, is rare for chemical systems. Here,
we extend the analysis by Turing to general networks and apply it to ecological
metapopulations of biological species with dispersal connections between
habitats. Remarkably, the oscillatory Turing instability does not lead to wave
patterns in networks, but to spontaneous development of heterogeneous
oscillations and possible extinction of some species, even though they are
absent for isolated populations. Furthermore, our theoretical analysis reveals
that this instability is more common in ecological metapopulations than in
chemical reactions. Indeed, we find the instabilities for all possible food
webs with three predator or prey species, under various assumptions about the
mobility of individual species and nonlinear interactions between them.
Therefore, we suggest that the oscillatory Turing instability is generic and
must play a fundamental role in metapopulation dynamics, providing a common
mechanism for dispersal-induced destabilization of ecosystems.


  The molecular motor myosin V exhibits a wide repertoire of pathways during
the stepping process, which is intimately connected to its biological function.
The best understood of these is hand-over-hand stepping by a swinging lever arm
movement toward the plus-end of actin filaments, essential to its role as a
cellular transporter. However, single-molecule experiments have also shown that
the motor "foot stomps", with one hand detaching and rebinding to the same
site, and backsteps under sufficient load. Explaining the complete taxonomy of
myosin V's load-dependent stepping pathways, and the extent to which these are
constrained by motor structure and mechanochemistry, are still open questions.
Starting from a polymer model, we develop an analytical theory to understand
the minimal physical properties that govern motor dynamics. In particular, we
solve the first-passage problem of the head reaching the target binding site,
investigating the competing effects of load pulling back at the motor, strain
in the leading head that biases the diffusion in the direction of the target,
and the possibility of preferential binding to the forward site due to the
recovery stroke. The theory reproduces a variety of experimental data,
including the power stroke and slow diffusive search regimes in the mean
trajectory of the detached head, and the force dependence of the
forward-to-backward step ratio, run length, and velocity. The analytical
approach yields a formula for the stall force, identifying the relative
contributions of the chemical cycle rates and mechanical features like the
bending rigidities of the lever arms. Most importantly, by fully exploring the
design space of the motor, we predict that myosin V is a robust motor whose
dynamical behavior is not compromised by reasonable perturbations to the
reaction cycle, and changes in the architecture of the lever arm.


  The search space of rooted phylogenetic trees is vast and a major research
focus of recent decades has been the development of algorithms to effectively
navigate this space. However this space is tiny when compared with the space of
rooted phylogenetic networks, and navigating this enlarged space remains a
poorly understood problem. This, and the difficulty of biologically
interpreting such networks, obstructs adoption of networks as tools for
modelling reticulation. Here, we argue that the superimposition of biologically
motivated constraints, via an SQL-style language, can both stimulate use of
network software by biologists and potentially significantly prune the search
space.


  Tissue heterogeneity is a major confounding factor in studying individual
populations that cannot be resolved directly by global profiling. Experimental
solutions to mitigate tissue heterogeneity are expensive, time consuming,
inapplicable to existing data, and may alter the original gene expression
patterns. Here we ask whether it is possible to deconvolute two-source mixed
expressions (estimating both proportions and cell-specific profiles) from two
or more heterogeneous samples without requiring any prior knowledge. Supported
by a well-grounded mathematical framework, we argue that both constituent
proportions and cell-specific expressions can be estimated in a completely
unsupervised mode when cell-specific marker genes exist, which do not have to
be known a priori, for each of constituent cell types. We demonstrate the
performance of unsupervised deconvolution on both simulation and real gene
expression data, together with perspective discussions.


  Differentiation is a key cellular process in normal tissue development that
is significantly altered in cancer. Although molecular signatures
characterising pluripotency and multipotency exist, there is, as yet, no single
quantitative mark of a cellular sample's position in the global differentiation
hierarchy. Here we adopt a systems view and consider the sample's network
entropy, a measure of signaling pathway promiscuity, computable from a sample's
genome-wide expression profile. We demonstrate that network entropy provides a
quantitative, in-silico, readout of the average undifferentiated state of the
profiled cells, recapitulating the known hierarchy of pluripotent, multipotent
and differentiated cell types. Network entropy further exhibits dynamic changes
in time course differentiation data, and in line with a sample's
differentiation stage. In disease, network entropy predicts a higher level of
cellular plasticity in cancer stem cell populations compared to ordinary cancer
cells. Importantly, network entropy also allows identification of key
differentiation pathways. Our results are consistent with the view that
pluripotency is a statistical property defined at the cellular population
level, correlating with intra-sample heterogeneity, and driven by the degree of
signaling promiscuity in cells. In summary, network entropy provides a
quantitative measure of a cell's undifferentiated state, defining its elevation
in Waddington's landscape.


  Mimulus guttatus and M. nasutus are an evolutionary and ecological model
sister species pair differentiated by ecology, mating system, and partial
reproductive isolation. Despite extensive research on this system, the history
of divergence and differentiation in this sister pair is unclear. We present
and analyze a novel population genomic data set which shows that M. nasutus
"budded" off of a central Californian M. guttatus population within the last
200 to 500 thousand years. In this time, the M. nasutus genome has accrued
numerous genomic signatures of the transition to predominant selfing. Despite
clear biological differentiation, we document ongoing, bidirectional
introgression. We observe a negative relationship between the recombination
rate and divergence between M. nasutus and sympatric M. guttatus samples,
suggesting that selection acts against M. nasutus ancestry in M. guttatus.


  We propose a method for obtaining joint probabilistic projections of
migration rates for all countries, broken down by age and sex. Joint
trajectories for all countries are constrained to satisfy the requirement of
zero global net migration. We evaluate our model using out-of-sample validation
and compare point projections to the projected migration rates from a
persistence model similar to the method used in the United Nations' World
Population Prospects, and also to a state of the art gravity model. We also
resolve an apparently paradoxical discrepancy between growth trends in the
proportion of the world population migrating and the average absolute migration
rate across countries.


  The hydration thermodynamics of the GXG tripeptide relative to the reference
GGG is often used to define the conditional hydration contribution of X. This
quantity or the hydration thermodynamics of a small molecule analog of the
side-chain or some combination of such estimates, have anchored the
interpretation of seminal experiments on protein stability and folding. We show
that such procedures to model protein hydration have significant limitations.
We study the conditional hydration thermodynamics of the isoleucine side-chain
in an extended pentapeptide and in helical deca-peptides, using as appropriate
an extended penta-glycine or appropriate helical deca-peptides as reference.
Hydration of butane in the gauche conformation provides a small molecule
reference for the side-chain. We use the quasichemical theory to parse the
hydration thermodynamics into chemical, packing, and long-range interaction
contributions. The chemical contribution reflects the contribution of solvent
clustering within the defined inner-shell of the solute; the chemical
contribution of g-butane is substantially more negative than the conditional
chemical contribution of isoleucine. The packing contribution gives the work
required to create a cavity in the solvent, a quantity of interest in
understanding hydrophobic hydration. The packing contribution for g-butane
substantially overestimates the conditional packing of isoleucine. The net of
such compensating contributions still disagrees with the conditional free
energy of isoleucine but by a lesser magnitude. The excess enthalpy and entropy
of hydration of g-butane model are also more negative than the corresponding
conditional quantities for the side-chain. The conditional solvation of
isoleucine in GGIGG also proves unsatisfactory in describing the conditional
solvation of isoleucine in the helical peptides.


  We present ProbMetab, an R package which promotes substantial improvement in
automatic probabilistic LC-MS based metabolome annotation. The inference engine
core is based on a Bayesian model implemented to: (i) allow diverse source of
experimental data and metadata to be systematically incorporated into the model
with alternative ways to calculate the likelihood function and; (ii) allow
sensitive selection of biologically meaningful biochemical reactions databases
as Dirichlet-categorical prior distribution. Additionally, to ensure result
interpretation by system biologists, we display the annotation in a network
where observed mass peaks are connected if their candidate metabolites are
substrate/product of known biochemical reactions. This graph can be overlaid
with other graph-based analysis, such as partial correlation networks, in a
visualization scheme exported to Cytoscape, with web and stand alone versions.
ProbMetab was implemented in a modular fashion to fit together with established
upstream (xcms, CAMERA, AStream, mzMatch.R, etc) and downstream R package tools
(GeneNet, RCytoscape, DiffCorr, etc). ProbMetab, along with extensive
documentation and case studies, is freely available under GNU license at:
http://labpib.fmrp.usp.br/methods/probmetab/.


  Background: Current research suggests that a small set of "driver" mutations
are responsible for tumorigenesis while a larger body of "passenger" mutations
occurs in the tumor but does not progress the disease. Due to recent
pharmacological successes in treating cancers caused by driver mutations, a
variety of of methodologies that attempt to identify such mutations have been
developed. Based on the hypothesis that driver mutations tend to cluster in key
regions of the protein, the development of cluster identification algorithms
has become critical.
  Results: We have developed a novel methodology, SpacePAC (Spatial Protein
Amino acid Clustering), that identifies mutational clustering by considering
the protein tertiary structure directly in 3D space. By combining the
mutational data in the Catalogue of Somatic Mutations in Cancer (COSMIC) and
the spatial information in the Protein Data Bank (PDB), SpacePAC is able to
identify novel mutation clusters in many proteins such as FGFR3 and CHRM2. In
addition, SpacePAC is better able to localize the most significant mutational
hotspots as demonstrated in the cases of BRAF and ALK. The R package is
available on Bioconductor at:
http://www.bioconductor.org/packages/release/bioc/html/SpacePAC.html
  Conclusion: SpacePAC adds a valuable tool to the identification of mutational
clusters while considering protein tertiary structure


  Cardiac development is complex, multiscale process encompassing cell fate
adoption, differentiation and morphogenesis. To elucidate pathways underlying
this process, a recently developed algorithm to reverse engineer gene
regulatory networks was applied to time-course microarray data obtained from
the developing mouse heart. The algorithm generates many different putative
network topologies that are capable of explaining the experimental data via
model simulation. To cull specious network interactions, thousands of
topologies are merged and filtered to generate a scale-free, hierarchical
network. The network is validated with known gene interactions and used to
identify regulatory pathways critical to the developing mammalian heart. The
predicted gene interactions are prioritized using semantic similarity and gene
profile uniqueness metrics. Using these metrics, the network is expanded to
include all known mouse genes to form the most likely cardiogenic gene
regulatory network. The method outlined herein provides an informative approach
to network inference and leads to clear testable hypotheses related to gene
regulation.


  Modeling biological networks serves as both a major goal and an effective
tool of systems biology in studying mechanisms that orchestrate the activities
of gene products in cells. Biological networks are context specific and dynamic
in nature. To systematically characterize the selectively activated regulatory
components and mechanisms, the modeling tools must be able to effectively
distinguish significant rewiring from random background fluctuations. We
formulated the inference of differential dependency networks that incorporates
both conditional data and prior knowledge as a convex optimization problem, and
developed an efficient learning algorithm to jointly infer the conserved
biological network and the significant rewiring across different conditions. We
used a novel sampling scheme to estimate the expected error rate due to random
knowledge and based on which, developed a strategy that fully exploits the
benefit of this data-knowledge integrated approach. We demonstrated and
validated the principle and performance of our method using synthetic datasets.
We then applied our method to yeast cell line and breast cancer microarray data
and obtained biologically plausible results.


  The study of human mobility is both of fundamental importance and of great
potential value. For example, it can be leveraged to facilitate efficient city
planning and improve prevention strategies when faced with epidemics. The
newfound wealth of rich sources of data---including banknote flows, mobile
phone records, and transportation data---has led to an explosion of attempts to
characterize modern human mobility. Unfortunately, the dearth of comparable
historical data makes it much more difficult to study human mobility patterns
from the past. In this paper, we present an analysis of long-term human
migration, which is important for processes such as urbanization and the spread
of ideas. We demonstrate that the data record from Korean family books (called
"jokbo") can be used to estimate migration patterns via marriages from the past
750 years. We apply two generative models of long-term human mobility to
quantify the relevance of geographical information to human marriage records in
the data, and we find that the wide variety in the geographical distributions
of the clans poses interesting challenges for the direct application of these
models. Using the different geographical distributions of clans, we quantify
the "ergodicity" of clans in terms of how widely and uniformly they have spread
across Korea, and we compare these results to those obtained using surname data
from the Czech Republic. To examine population flow in more detail, we also
construct and examine a population-flow network between regions. Based on the
correlation between ergodicity and migration in Korea, we identify two
different types of migration patterns: diffusive and convective. We expect the
analysis of diffusive versus convective effects in population flows to be
widely applicable to the study of mobility and migration patterns across
different cultures.


  Intercellular heterogeneity serves as both a confounding factor in studying
individual clones and an information source in characterizing any heterogeneous
tissues, such as blood, tumor systems. Due to inevitable sequencing errors and
other sample preparation artifacts such as PCR errors, systematic efforts to
characterize intercellular genomic heterogeneity must effectively distinguish
genuine clonal sequences from fake derivatives. We developed a novel approach
(SIGH) for identifying significant genuine clonal sequences directly from mixed
sequencing reads that can improve genomic analyses in many biological contexts.
This method offers several attractive features: (1) it automatically estimates
the error rate from raw sequence reads and identifies genuine clonal sequences;
(2) it is robust to the large variety of error rate due to the various
experimental conditions; (3) it is supported by a well grounded statistical
framework that exploits probabilistic characteristics of sequencing errors; (4)
its unbiased strategy allows detecting rare clone(s) despite that clone
relative abundance; and (5) it estimates constituent proportions in each
sample. Extensive realistic simulation studies show that our method can
reliably estimate the error rates and faithfully distinguish the genuine clones
from fake derivatives, paving the way for follow up analysis that is otherwise
ruined by the often dominant fake clones.


  Language universals have long been attributed to an innate Universal Grammar.
An alternative explanation states that linguistic universals emerged
independently in every language in response to shared cognitive or perceptual
biases. A computational model has recently shown how this could be the case,
focusing on the paradigmatic example of the universal properties of colour
naming patterns, and producing results in quantitative agreement with the
experimental data. Here we investigate the role of an individual perceptual
bias in the framework of the model. We study how, and to what extent, the
structure of the bias influences the corresponding linguistic universal
patterns. We show that the cultural history of a group of speakers introduces
population-specific constraints that act against the pressure for uniformity
arising from the individual bias, and we clarify the interplay between these
two forces.


  It has been recently claimed that it is possible to predict the rate of de
novo mutation of each site in the human genome with almost perfect accuracy
(Michaelson et al. (2012) Cell, 151, 1431-1442). We show that this claim is
unwarranted. By considering the correlation between the rate of de novo
mutation and the predictions from the model of Michaelson et al., we show that
there could be substantial unexplained variance in the mutation rate. We also
demonstrate that the model of Michaelson et al. fails to capture a major
component of the variation in the mutation rate, that which is local but not
associated with simple context.


  BACOM is a statistically principled and unsupervised method that detects copy
number deletion types (homozygous versus heterozygous), estimates normal cell
fraction, and recovers cancer specific copy number profiles, using allele
specific copy number signals. In a subsequent analysis of TCGA ovarian cancer
dataset, the average normal cell fraction estimated by BACOM was found higher
than expected. In this letter, we first discuss the advantages of the BACOM in
relation to alternative approaches. Then, we show that this elevated estimate
of normal cell fraction is the combined result of inaccurate signal modeling
and normalization. Lastly, we describe an allele specific signal modeling and
normalization scheme that can enhance BACOM applications in many biological
contexts. An open source MATLAB program was developed to implement our extended
method and it is publically available.


  A series of studies have revealed the among-population components of genetic
variation are higher for the paternal Y chromosome than for the maternal
mitochondrial DNA (mtDNA), which indicates sex-biased migrations in human
populations. However, this phenomenon might be also an ascertainment bias due
to nonrandom sampling of SNPs. To eliminate the possible bias, we used the
whole Y chromosome and mtDNA sequence data of 491 individuals from the 1000
Genomes Project Phase I to address the sex-biased migration dispute. We found
that genetic differentiation between populations was higher for Y chromosome
than for the mtDNA at global scales. The migration rate of female might be
three times higher than that of male, assuming the effective population size is
the same for male and female.


  The evolution of autocatalytic sets (ACS) is a widespread process in
biological, chemical and ecological systems which is of great significance in
many applications, such as the evolution of new species or complex chemical
organizations. In this paper, we propose a competitive model with a m-selection
rule in which an abrupt emergence of a macroscopic independent ACS is observed.
By numerical simulations, we find that the maximal increase of the size grows
linearly with the system size. We analytically derive the threshold t{\alpha}
where the abrupt jump happens and verify it by simulations. Moreover, our
analysis explains how this giant independent ACS grows and reveals that, as the
selection rule becomes more strict, the phase transition is dramatically
postponed, and the number of the largest independent ACSs coexisting in the
system increases accordingly. Our research work deepens the understanding of
the evolution of ACS and should provide useful information for designing
strategies to control the emergence of ACS in corresponding applications.


  Several recent studies show that companion cells in flowering plant
gametophytes relax epigenetic control of transposable elements (TEs) to promote
production of small RNA that presumably assist nearby reproductive cells in
management of TEs. In light of this possibility, a closer look at the timing of
cell division in relation to angiosperm double fertilization is warranted. From
such an analysis, it is conceivable that double fertilization can drive
angiosperm evolution by facilitating crosses between genetically diverse
parents. A key feature of this ability is the order of cell division following
double fertilization, since division of the primary endosperm nucleus prior to
the zygote would produce small RNA capable of identifying TEs and defining
chromatin states in the zygote prior to its entry into S-phase of the cell
cycle. Consequently, crosses leading to increased ploidy or between genetically
diverse parents would yield offspring better capable of managing a diverse
complement of TEs and repetitive DNA. Considering double fertilization in this
regard challenges previous notions that the primary purpose of endosperm is for
improved seed reserve storage and utilization.


  One outcome of interspecific hybridization and subsequent effects of
evolutionary forces is introgression, which is the integration of genetic
material from one species into the genome of an individual in another species.
The evolution of several groups of eukaryotic species has involved
hybridization, and cases of adaptation through introgression have been already
established. In this work, we report on a new comparative genomic framework for
detecting introgression in genomes, called PhyloNet-HMM, which combines
phylogenetic networks, that capture reticulate evolutionary relationships among
genomes, with hidden Markov models (HMMs), that capture dependencies within
genomes. A novel aspect of our work is that it also accounts for incomplete
lineage sorting and dependence across loci.
  Application of our model to variation data from chromosome 7 in the mouse
(Mus musculus domesticus) genome detects a recently reported adaptive
introgression event involving the rodent poison resistance gene Vkorc1, in
addition to other newly detected introgression regions. Based on our analysis,
it is estimated that about 12% of all sites withinchromosome 7 are of
introgressive origin (these cover about 18 Mbp of chromosome 7, and over 300
genes). Further, our model detects no introgression in two negative control
data sets. Our work provides a powerful framework for systematic analysis of
introgression while simultaneously accounting for dependence across sites,
point mutations, recombination, and ancestral polymorphism.


  One of the central themes in modern ecology is the enduring debate on whether
there is a relationship between the complexity of a biological community and
its stability. In this paper, we focus on the role of detritus and spatial
dispersion on the stability of ecosystems. Using Monte Carlo simulations we
analyze two three level models of food webs: a grazing one with the basal
species (i.e. primary producers) having unlimited food resources and a detrital
one in which the basal species uses detritus as a food resource. While the vast
majority of theoretical studies neglects detritus, from our results it follows
that the detrital food web is more stable than its grazing counterpart, because
the interactions mediated by detritus damp out fluctuations in species'
densities. Since the detritus model is the more complex one in terms of
interaction patterns, our results provide new evidence for the advocates of the
complexity as one of the factors enhancing stability of ecosystems.


  We demonstrate how a genetic ring oscillator network with quorum sensing
feedback can operate as a robust logic gate. Specifically we show how a range
of logic functions, namely AND/NAND, OR/NOR and XOR/XNOR, can be realized by
the system, thus yielding a versatile unit that can morph between different
logic operations. We further demonstrate the capacity of this system to yield
complementary logic operations in parallel. Our results then indicate the
computing potential of this biological system, and may lead to bio-inspired
computing devices.


  Gene covariation networks are commonly used to study biological processes.
The inference of gene covariation networks from observational data can be
challenging, especially considering the large number of players involved and
the small number of biological replicates available for analysis. We propose a
new statistical method for estimating the number of erroneous edges in
reconstructed networks that strongly enhances commonly used inference
approaches. This method is based on a special relationship between sign of
correlation (positive/negative) and directionality (up/down) of gene
regulation, and allows for the identification and removal of approximately half
of all erroneous edges. Using the mathematical model of Bayesian networks and
positive correlation inequalities we establish a mathematical foundation for
our method. Analyzing existing biological datasets, we find a strong
correlation between the results of our method and false discovery rate (FDR).
Furthermore, simulation analysis demonstrates that our method provides a more
accurate estimate of network error than FDR.


  Motivation: Computational methods are essential to extract actionable
information from raw sequencing data, and to thus fulfill the promise of
next-generation sequencing technology. Unfortunately, computational tools
developed to call variants from human sequencing data disagree on many of their
predictions, and current methods to evaluate accuracy and computational
performance are ad-hoc and incomplete. Agreement on benchmarking variant
calling methods would stimulate development of genomic processing tools and
facilitate communication among researchers.
  Results: We propose SMaSH, a benchmarking methodology for evaluating human
genome variant calling algorithms. We generate synthetic datasets, organize and
interpret a wide range of existing benchmarking data for real genomes, and
propose a set of accuracy and computational performance metrics for evaluating
variant calling methods on this benchmarking data. Moreover, we illustrate the
utility of SMaSH to evaluate the performance of some leading single nucleotide
polymorphism (SNP), indel, and structural variant calling algorithms.
  Availability: We provide free and open access online to the SMaSH toolkit,
along with detailed documentation, at smash.cs.berkeley.edu.


  One of the most challenging problems in microbiology is to understand how a
small fraction of microbes that resists killing by antibiotics can emerge in a
population of genetically identical cells, the phenomenon known as persistence
or drug tolerance. Its characteristic signature is the biphasic kill curve,
whereby microbes exposed to a bactericidal agent are initially killed very
rapidly but then much more slowly. Here we relate this problem to the more
general problem of understanding the emergence of distinct growth phenotypes in
clonal populations. We address the problem mathematically by adopting the
framework of the phenomenon of so-called weak ergodicity breaking, well known
in dynamical physical systems, which we extend to the biological context. We
show analytically and by direct stochastic simulations that distinct growth
phenotypes can emerge as a consequence of slow-down of stochastic fluctuations
in the expression of a gene controlling growth rate. In the regime of fast gene
transcription, the system is ergodic, the growth rate distribution is unimodal,
and accounts for one phenotype only. In contrast, at slow transcription and
fast translation, weakly non-ergodic components emerge, the population
distribution of growth rates becomes bimodal, and two distinct growth
phenotypes are identified. When coupled to the well-established growth rate
dependence of antibiotic killing, this model describes the observed fast and
slow killing phases, and reproduces much of the phenomenology of bacterial
persistence. The model has major implications for efforts to develop control
strategies for persistent infections.


  Relevance of mode coupling to energy/information transfer during protein
function, particularly in the context of allosteric interactions is widely
accepted. However, existing evidence in favor of this hypothesis comes
essentially from model systems. We here report a novel formal analysis of the
near-native dynamics of myosin II, which allows us to explore the impact of the
interaction between possibly non-Gaussian vibrational modes on fluctutational
dynamics. We show that, an information-theoretic measure based on mode coupling
{\it alone} yields a ranking of residues with a statistically significant bias
favoring the functionally critical locations identified by experiments on
myosin II.


  Speciation is driven by many different factors. Among those are trade-offs
between different ways an organism utilizes resources, and these trade-offs can
constrain the manner in which selection can optimize traits. Limited migration
among allopatric populations and species interactions can also drive
speciation, but here we ask if trade-offs alone are sufficient to drive
speciation in the absence of other factors. We present a model to study the
effects of trade-offs on specialization and adaptive radiation in asexual
organisms based solely on competition for limiting resources, where trade-offs
are stronger the greater an organism's ability to utilize resources. In this
model resources are perfectly substitutable, and fitness is derived from the
consumption of these resources. The model contains no spatial parameters, and
is therefore strictly sympatric. We quantify the degree of specialization by
the number of ecotypes formed and the niche breadth of the population, and
observe that these are sensitive to resource influx and trade-offs. Resource
influx has a strong effect on the degree of specialization, with a clear
transition between minimal diversification at high influx and multiple species
evolving at low resource influx. At low resource influx the degree of
specialization further depends on the strength of the trade-offs, with more
ecotypes evolving the stronger trade-offs are. The specialized organisms
persist through negative frequency-dependent selection. In addition, by
analyzing one of the evolutionary radiations in greater detail we demonstrate
that a single mutation alone is not enough to establish a new ecotype, even
though phylogenetic reconstruction identifies that mutation as the branching
point. Instead, it takes a series of additional mutations to ensure the stable
coexistence of the new ecotype in the background of the existing ones,
reminiscent of a recent observa


  Transcription of genes is the focus of most forms of regulation of gene
expression. Even though careful biochemical experimentation has revealed the
molecular mechanisms of transcription initiation for a number of different
promoters in vitro, the dynamics of this process in cells is still poorly
understood. One approach has been to measure the transcriptional output
(fluorescently labeled messenger RNAs or proteins) from single cells in a
genetically identical population, which could then be compared to predictions
from models that incorporate different molecular mechanisms of transcription
initiation. However, this approach suffers from the problem, that processes
downstream from transcription can affect the measured output and therefore mask
the signature of stochastic transcription initiation on the cell-to-cell
variability of the transcriptional outputs. Here we show theoretically that
measurements of the cell-to-cell variability in the number of nascent RNAs
provide a more direct test of the mechanism of transcription initiation. We
derive exact expressions for the first two moments of the distribution of
nascent RNA molecules and apply our theory to published data for a collection
of constitutively expressed yeast genes. We find that the measured nascent RNA
distributions are inconsistent with transcription initiation proceeding via one
rate-limiting step, which has been generally inferred from measurements of
cytoplasmic messenger RNA. Instead, we propose a two-step mechanism of
initiation, which is consistent with the available data. These findings for the
yeast promoters highlight the utility of our theory for deciphering
transcriptional dynamics in vivo from experiments that count nascent RNA
molecules in single cells.


  The mapping between biological genotypes and phenotypes is central to the
study of biological evolution. Here we introduce a rich, intuitive, and
biologically realistic genotype-phenotype (GP) map, that serves as a model of
self-assembling biological structures, such as protein complexes, and remains
computationally and analytically tractable. Our GP map arises naturally from
the self-assembly of polyomino structures on a 2D lattice and exhibits a number
of properties: $\textit{redundancy}$ (genotypes vastly outnumber phenotypes),
$\textit{phenotype bias}$ (genotypic redundancy varies greatly between
phenotypes), $\textit{genotype component disconnectivity}$ (phenotypes consist
of disconnected mutational networks) and $\textit{shape space covering}$ (most
phenotypes can be reached in a small number of mutations). We also show that
the mutational robustness of phenotypes scales very roughly logarithmically
with phenotype redundancy and is positively correlated with phenotypic
evolvability. Although our GP map describes the assembly of disconnected
objects, it shares many properties with other popular GP maps for connected
units, such as models for RNA secondary structure or the HP lattice model for
protein tertiary structure. The remarkable fact that these important properties
similarly emerge from such different models suggests the possibility that
universal features underlie a much wider class of biologically realistic GP
maps.


  A biomolecular system representing the first realization of associative
memory based on enzymatic reactions in vitro has been designed. The system
demonstrated "training" and "forgetting" features characteristic of memory in
biological systems, but presently realized in simple biocatalytic cascades.


  We report the first systematic study of designed two-input biochemical
systems as information processing gates with favorable noise-transmission
properties accomplished by modifying the gate's response from convex shape to
sigmoid in both inputs. This is realized by an added chemical "filter" process
which recycles some of the output back into one of the inputs. We study a
system involving the biocatalytic function of the enzyme horseradish
peroxidase, functioning as an AND gate. We consider modularity properties, such
as the use of three different input chromogens that, when oxidized yield
signal-detection outputs for various ranges of the primary input, hydrogen
peroxide. We also examine possible uses of different filter-effect chemicals
(reducing agents) to induce the sigmoid-response. A modeling approach is
developed and applied to our data, allowing us to describe the enzymatic
kinetics in the framework of a formulation suitable for evaluating the
noise-handling properties of the studied systems as logic gates for information
processing steps.