  We investigate invasions from a biological reservoir to an initially empty,
heterogeneous habitat in the presence of advection. The habitat consists of a
periodic alternation of favorable and unfavorable patches. In the latter the
population dies at fixed rate. In the former it grows either with the logistic
or with an Allee effect type dynamics, where the population has to overcome a
threshold to grow. We study the conditions for successful invasions and the
speed of the invasion process, which is numerically and analytically
investigated in several limits. Generically advection enhances the downstream
invasion speed but decreases the population size of the invading species, and
can even inhibit the invasion process. Remarkably, however, the rate of
population increase, which quantifies the invasion efficiency, is maximized by
an optimal advection velocity. In models with Allee effect, differently from
the logistic case, above a critical unfavorable patch size the population
localizes in a favorable patch, being unable to invade the habitat. However, we
show that advection, when intense enough, may activate the invasion process.


  An ubiquitous property of biological sensory systems is adaptation: a step
increase in stimulus triggers an initial change in a biochemical or
physiological response, followed by a more gradual relaxation toward a basal,
pre-stimulus level. Adaptation helps maintain essential variables within
acceptable bounds and allows organisms to readjust themselves to an optimum and
non-saturating sensitivity range when faced with a prolonged change in their
environment. Recently, it was shown theoretically and experimentally that many
adapting systems, both at the organism and single-cell level, enjoy a
remarkable additional feature: scale invariance, meaning that the initial,
transient behavior remains (approximately) the same even when the background
signal level is scaled. In this work, we set out to investigate under what
conditions a broadly used model of biochemical enzymatic networks will exhibit
scale-invariant behavior. An exhaustive computational study led us to discover
a new property of surprising simplicity and generality, uniform linearizations
with fast output (ULFO), whose validity we show is both necessary and
sufficient for scale invariance of enzymatic networks. Based on this study, we
go on to develop a mathematical explanation of how ULFO results in scale
invariance. Our work provides a surprisingly consistent, simple, and general
framework for understanding this phenomenon, and results in concrete
experimental predictions.


  A recurring motif in gene regulatory networks is transcription factors (TFs)
that regulate each other, and then bind to overlapping sites on DNA, where they
interact and synergistically control transcription of a target gene. Here, we
suggest that this motif maximizes information flow in a noisy network. Gene
expression is an inherently noisy process due to thermal fluctuations and the
small number of molecules involved. A consequence of multiple TFs interacting
at overlapping binding-sites is that their binding noise becomes correlated.
Using concepts from information theory, we show that in general a signaling
pathway transmits more information if 1) noise of one input is correlated with
that of the other, 2) input signals are not chosen independently. In the case
of TFs, the latter criterion hints at up-stream cross-regulation. We
demonstrate these ideas for competing TFs and feed-forward gene regulatory
modules, and discuss generalizations to other signaling pathways. Our results
challenge the conventional approach of treating biological noise as
uncorrelated fluctuations, and present a systematic method for understanding TF
cross-regulation networks either from direct measurements of binding noise, or
bioinformatic analysis of overlapping binding-sites.


  Drug gradients are believed to play an important role in the evolution of
bacteria resistant to antibiotics and tumors resistant to anti-cancer drugs. We
use a statistical physics model to study the evolution of a population of
malignant cells exposed to drug gradients, where drug resistance emerges via a
mutational pathway involving multiple mutations. We show that a non-uniform
drug distribution has the potential to accelerate the emergence of resistance
when the mutational pathway involves a long sequence of mutants with increasing
resistance, but if the pathway is short or crosses a fitness valley, the
evolution of resistance may actually be slowed down by drug gradients. These
predictions can be verified experimentally, and may help to improve strategies
for combatting the emergence of resistance.


  The thermal mean square displacement (MSD) of hydrogen in proteins and its
associated hydration water is measured by neutron scattering experiments and
used an indicator of protein function. The observed MSD as currently determined
depends on the energy resolution width of the neutron scattering instrument
employed. We propose a method for obtaining the intrinsic MSD of H in the
proteins, one that is independent of the instrument resolution width. The
intrinsic MSD is defined as the infinite time value of <r^2> that appears in
the Debye-Waller factor. The method consists of fitting a model to the
resolution broadened elastic incoherent structure factor or to the resolution
dependent MSD. The model contains the intrinsic MSD, the instrument resolution
width and a rate constant characterizing the motions of H in the protein. The
method is illustrated by obtaining the intrinsic MSD <r^2> of heparan sulphate
(HS-0.4), Ribonuclease A and Staphysloccal Nuclase (SNase) from data in the
literature.


  We propose a novel framework for studying causal inference of gene
interactions using a combination of compressive sensing and Granger causality
techniques. The gist of the approach is to discover sparse linear dependencies
between time series of gene expressions via a Granger-type elimination method.
The method is tested on the Gardner dataset for the SOS network in E. coli, for
which both known and unknown causal relationships are discovered.


  In various chemical systems enthalpy-entropy compensation (EEC) is a
well-known rule of behavior, although the physical roots of it are still not
completely understood. It has been frequently questioned whether EEC is a truly
physical phenomenon or a coincidence due to trivial mathematical connections
between statistical-mechanical parameters - or even simpler: A phantom effect
resulting from the misinterpretation of experimental data. Here, we review EEC
from a new standpoint using the notion of correlation which is essential for
the method of factor analysis, but is not conventional in physics and
chemistry. We conclude that the EEC may be rationalized in terms of hidden (not
directly measurable with the help of the current experimental set-up) but
physically real factors, implying a Carnot-cycle model in which a micro-phase
transition (MPT) plays a crucial role. Examples of such MPTs underlying
physically valid EEC should be typically cooperative processes in
supramolecular aggregates, like changes of structured water at hydrophobic
surfaces, conformational transitions upon ligand-biopolymer binding, and so on,
so forth. The MPT notion could help rationalize the occurrence of EEC in
connection with hydration and folding of proteins,enzymatic reactions,
functioning of molecular motors, DNA de- and rehybridization, as well as
similar phenomena.


  We consider a ${\Lambda}$-coalescent and we study the asymptotic behavior of
the total length $L^{(n)}_{ext}$ of the external branches of the associated
$n$-coalescent. For Kingman coalescent, i.e. ${\Lambda}={\delta}_0$, the result
is well known and is useful, together with the total length $L^{(n)}$, for Fu
and Li's test of neutrality of mutations% under the infinite sites model
asumption . For a large family of measures ${\Lambda}$, including
Beta$(2-{\alpha},{\alpha})$ with $0<\alpha<1$, M{\"o}hle has proved asymptotics
of $L^{(n)}_{ext}$. Here we consider the case when the measure ${\Lambda}$ is
Beta$(2-{\alpha},{\alpha})$, with $1<\alpha<2$. We prove that
$n^{{\alpha}-2}L^{(n)}_{ext}$ converges in $L^2$ to
$\alpha(\alpha-1)\Gamma(\alpha)$. As a consequence, we get that
$L^{(n)}_{ext}/L^{(n)}$ converges in probability to $2-\alpha$. To prove the
asymptotics of $L^{(n)}_{ext}$, we use a recursive construction of the
$n$-coalescent by adding individuals one by one. Asymptotics of the
distribution of $d$ normalized external branch lengths and a related moment
result are also given.


  Precopulatory courtship is a high-cost, non-well understood animal world
mystery. Drosophila's (=D.'s) precopulatory courtship not only shows marked
structural similarities with mammalian courtship, but also with human spoken
language. This suggests the study of purpose, modalities and in particular of
the power of this language and to compare it to human language. Following a
mathematical symbolic dynamics approach, we translate courtship videos of D.'s
body language into a formal language. This approach made it possible to show
that D. may use its body language to express individual information -
information that may be important for evolutionary optimization, on top of the
sexual group membership. Here, we use Chomsky's hierarchical language
classification to characterize the power of D.'s body language, and then
compare it with the power of languages spoken by humans. We find that from a
formal language point of view, D.'s body language is at least as powerful as
the languages spoken by humans. From this we conclude that human intellect
cannot be the direct consequence of the formal grammar complexity of human
language.


  We studied in silico effects of structural stress in Mycoplasma pneumoniae, a
genome-reduced model bacterial organism, by tracking the damage propagating on
its metabolic network after a deleterious perturbation. First, we analyzed
failure cascades spreading from individual reactions and pairs of reactions and
compared the results to those in Staphylococcus aureus and Escherichia coli. To
alert to the potential damage caused by the failure of individual reactions, we
propose a generic predictor based on local information that identifies target
reactions for structural vulnerability. With respect to the simultaneous
failure of pairs of reactions, we detected strong non-linear amplification
effects that can be predicted by the presence of specific motifs in the
intersection of single cascades. We further connected the metabolic and gene
co-expression networks of M. pneumoniae through enzyme activities, and studied
the consequences of knocking out individual genes and clusters of genes. Damage
caused by single gene knockouts reveals a strong correlation between
genome-scale cascades of large impact and gene essentiality. At the same time,
we found that genes controlling high-damage reactions tend to operate in
functional isolation, as a metabolic protection mechanism. We conclude that the
architecture of M. pneumoniae, both at the level of metabolism and genome,
seems to have evolved towards increased structural robustness, similarly to
other more complex model bacterial organisms, despite its reduced genome size
and its greater metabolic network linearity. Our approach, although motivated
biochemically, is generic enough to be of potential use toward analyzing and
predicting spreading of structural stress in any bipartite complex network.


  We review the statistical properties of the genealogies of a few models of
evolution. In the asexual case, selection leads to coalescence times which grow
logarithmically with the size of the population in contrast with the linear
growth of the neutral case. Moreover for a whole class of models, the
statistics of the genealogies are those of the Bolthausen-Sznitman coalescent
rather than the Kingman coalescent in the neutral case. For sexual
reproduction, the time to reach the first common ancestors to the whole
population and the time for all individuals to have all their ancestors in
common are also logarithmic in the neutral case, as predicted by Chang []. We
discuss how these times are modified in a simple way of introducing selection.


  We propose a (time) multiscale method for the coarse-grained analysis of
self--propelled particle models of swarms comprising a mixture of `na\"{i}ve'
and `informed' individuals, used to address questions related to collective
motion and collective decision--making in animal groups. The method is based on
projecting the particle configuration onto a single `meta-particle' that
consists of the group elongation and the mean group velocity and position. The
collective states of the configuration can be associated with the transient and
asymptotic transport properties of the random walk followed by the
meta-particle. These properties can be accurately predicted by an
advection-diffusion equation with memory (ADEM) whose parameters are obtained
from a mean group velocity time series obtained from a single simulation run of
the individual-based model.


  The pervasive presence spatial and size structure in biological populations
challenges fundamental assumptions at the heart of continuum models of
population dynamics based on mean densities (local or global) only.
Individual-based models (IBM's) were introduced over the last decade in an
attempt to overcome this limitation by following explicitly each individual in
the population. Although the IBM approach has been quite insightful, the
capability to follow each individual usually comes at the expense of analytical
tractability, which limits the generality of the statements that can be made.
For the specific case of spatial structure in populations of sessile (and
identical) organisms, space-time point processes with local regulation seem to
cover the middle ground between analytical tractability and a higher degree of
biological realism. Continuum approximations of these stochastic processes
distill their fundamental properties, but they often result in infinite
hierarchies of moment equations. We use the principle of constrained maximum
entropy to derive a closure relationship for one such hierarchy truncated at
second order using normalization and the product densities of first and second
orders as constraints. The resulting `maxent' closure is similar to the
Kirkwood superposition approximation, but it is complemented with previously
unknown correction terms that depend on on the area for which third order
correlations are irreducible. This region also serves as a validation check,
since it can only be found if the assumptions of the closure are met.
Comparisons between simulations of the point process, alternative heuristic
closures, and the maxent closure show significant improvements in the ability
of the maxent closure to predict equilibrium values for mildly aggregated
spatial patterns.


  We study the synchronisation and phase-lag of fluctuations in the number of
infected individuals in a network of cities between which individuals commute.
The frequency and amplitude of these oscillations is known to be very well
captured by the van Kampen system-size expansion, and we use this approximation
to compute the complex coherence function that describes their correlation. We
find that, if the infection rate differs from city to city and the coupling
between them is not too strong, these oscillations are synchronised with a well
defined phase lag between cities. The analytic description of the effect is
shown to be in good agreement with the results of stochastic simulations for
realistic population sizes.


  Understanding the differentiation, a biological process from a multipotent
stem or progenitor state to a mature cell is critically important. We develop a
theoretical framework to quantify the underlying potential landscape and
biological paths for cell development and differentiation. We propose a new
mechanism of differentiation and development through binding/unbinding of
regulatory proteins to the gene promoters. We found indeed the differentiated
states can emerge from the slow promoter binding/unbinding processes.
Furthermore, under slow promoter binding/unbinding, we found multiple
meta-stable differentiated states. This can explain the origin of multiple
states observed in the recent experiments. In addition, the kinetic time
quantified by mean first passage transition time for the differentiation and
reprogramming strongly depends on the time scale of the promoter
binding/unbinding processes. We discovered an optimal speed for differentiation
for certain binding/unbinding rates of regulatory proteins to promoters. More
experiments in the future might be able to tell if cells differentiate at at
that optimal speed. In addition, we quantify kinetic pathways for the
differentiation and reprogramming. We found that they are irreversible. This
captures the non-equilibrium dynamics in multipotent stem or progenitor cells.
Such inherent time-asymmetry as a result of irreversibility of state transition
pathways as shown may provide the origin of time arrow for cell development.


  We give a new explanation for why some biological systems can stay quantum
coherent for long times at room temperatures, one of the fundamental puzzles of
quantum biology. We show that systems with the right level of complexity
between chaos and regularity can increase their coherence time by orders of
magnitude. Systems near Critical Quantum Chaos or Metal-Insulator Transition
(MIT) can have long coherence times and coherent transport at the same time.
The new theory tested in a realistic light harvesting system model can
reproduce the scaling of critical fluctuations reported in recent experiments.
Scaling of return probability in the FMO light harvesting complex shows the
signs of universal return probability decay observed at critical MIT. The
results may open up new possibilities to design low loss energy and information
transport systems in this Poised Realm hovering reversibly between quantum
coherence and classicality.


  Macromolecular crowding in living biological cells effects subdiffusion of
larger biomolecules such as proteins and enzymes. Mimicking this subdiffusion
in terms of random walks on a critical percolation cluster, we here present a
case study of EcoRV restriction enzymes involved in vital cellular defence. We
show that due to its so far elusive propensity to an inactive state the enzyme
avoids non-specific binding and remains well-distributed in the bulk cytoplasm
of the cell. Despite the reduced volume exploration capability of subdiffusion
processes, this mechanism guarantees a high efficiency of the enzyme. By
variation of the non-specific binding constant and the bond occupation
probability on the percolation network, we demonstrate that reduced
non-specific binding are beneficial for efficient subdiffusive enzyme activity
even in relatively small bacteria cells. Our results corroborate a more local
picture of cellular regulation.


  One of the main applications of balance indices is in tests of null models of
evolutionary processes. The knowledge of an exact formula for a statistic of a
balance index, holding for any number n of leaves, is necessary in order to use
this statistic in tests of this kind involving trees of any size. In this paper
we obtain exact formulas for the variance under the Yule model of the Sackin
index, the Colless index and the total cophenetic index of binary rooted
phylogenetic trees with n leaves. We also obtain the covariance of the Sackin
and the total cophenetic index.


  We investigate distances on binary (presence/absence) data in the context of
a Dollo process, where a trait can only arise once on a phylogenetic tree but
may be lost many times. We introduce a novel distance, the Additive Dollo
Distance (ADD), which is consistent for data generated under a Dollo model, and
show that it has some useful theoretical properties including an intriguing
link to the LogDet distance. Simulations of Dollo data are used to compare a
number of binary distances including ADD, LogDet, Nei Li and some simple, but
to our knowledge previously unstudied, variations on common binary distances.
The simulations suggest that ADD outperforms other distances on Dollo data.
Interestingly, we found that the LogDet distance performs poorly in the context
of a Dollo process, which may have implications for its use in connection with
conditioned genome reconstruction. We apply the ADD to two Diversity Arrays
Technology (DArT) datasets, one that broadly covers Eucalyptus species and one
that focuses on the Eucalyptus series Adnataria. We also reanalyse gene family
presence/absence data on bacteria from the COG database and compare the results
to previous phylogenies estimated using the conditioned genome reconstruction
approach.


  The biological activity and functional specificity of proteins depend on
their native three-dimensional structures determined by inter- and
intra-molecular interactions. In this paper, we investigate the geometrical
factor of protein conformation as a consequence of energy minimization in
protein folding. Folding simulations of 10 polypeptides with chain length
ranging from 183 to 548 residues manifest that the dimensionless ratio
(V/(A<r>)) of the van der Waals volume V to the surface area A and average
atomic radius <r> of the folded structures, calculated with atomic radii
setting used in SMMP [Eisenmenger F., et. al., Comput. Phys. Commun., 138
(2001) 192], approach 0.49 quickly during the course of energy minimization. A
large scale analysis of protein structures show that the ratio for real and
well-designed proteins is universal and equal to 0.491\pm0.005. The fractional
composition of hydrophobic and hydrophilic residues does not affect the ratio
substantially. The ratio also holds for intrinsically disordered proteins,
while it ceases to be universal for polypeptides with bad folding properties.


  We consider networks with two types of nodes. The v-nodes, called centers,
are hyper- connected and interact one to another via many u-nodes, called
satellites. This central- ized architecture, widespread in gene networks,
possesses two fundamental properties. Namely, this organization creates
feedback loops that are capable to generate practically any prescribed
patterning dynamics, chaotic or periodic, or having a number of equilib- rium
states. Moreover, this organization is robust with respect to random
perturbations of the system.


  In this paper, we provide a polynomial time algorithm to calculate the
probability of a {\it ranked} gene tree topology for a given species tree,
where a ranked tree topology is a tree topology with the internal vertices
being ordered. The probability of a gene tree topology can thus be calculated
in polynomial time if the number of orderings of the internal vertices is a
polynomial number. However, the complexity of calculating the probability of a
gene tree topology with an exponential number of rankings for a given species
tree remains unknown.


  Was the spread of agropastoralism from the Fertile Crescent throughout Europe
influenced by extreme climate events, or was it independent of climate? We here
generate idealized climate events using palaeoclimate records. In a
mathematical model of regional sociocultural development, these events disturb
the subsistence base of simulated forager and farmer societies. We evaluate the
regional simulated transition timings and durations against a published large
set of radiocarbon dates for western Eurasia; the model is able to
realistically hindcast much of the inhomogeneous space-time evolution of
regional Neolithic transitions. Our study shows that the consideration of
climate events improves the simulation of typical lags between cultural
complexes, but that the overall difference to a model without climate events is
not significant. Climate events may not have been as important for early
sociocultural dynamics as endogenous factors.


  The human organism is an integrated network where complex physiologic
systems, each with its own regulatory mechanisms, continuously interact, and
where failure of one system can trigger a breakdown of the entire network.
Identifying and quantifying dynamical networks of diverse systems with
different types of interactions is a challenge. Here, we develop a framework to
probe interactions among diverse systems, and we identify a physiologic
network. We find that each physiologic state is characterized by a specific
network structure, demonstrating a robust interplay between network topology
and function. Across physiologic states the network undergoes topological
transitions associated with fast reorganization of physiologic interactions on
time scales of a few minutes, indicating high network flexibility in response
to perturbations. The proposed system-wide integrative approach may facilitate
the development of a new field, Network Physiology.


  When analysing gene expression time series data an often overlooked but
crucial aspect of the model is that the regulatory network structure may change
over time. Whilst some approaches have addressed this problem previously in the
literature, many are not well suited to the sequential nature of the data. Here
we present a method that allows us to infer regulatory network structures that
may vary between time points, utilising a set of hidden states that describe
the network structure at a given time point. To model the distribution of the
hidden states we have applied the Hierarchical Dirichlet Process Hideen Markov
Model, a nonparametric extension of the traditional Hidden Markov Model, that
does not require us to fix the number of hidden states in advance. We apply our
method to exisiting microarray expression data as well as demonstrating is
efficacy on simulated test data.


  Diversity patterns of tree species in a tropical forest community are
approached by a simple lattice model and investigated by Monte Carlo
simulations using a backtracking method. Our spatially explicit neutral model
is based on a simple statistical physics process, namely the diffusion of
seeds. The model has three parameters: the speciation rate, the size of the
meta-community in which the studied tree-community is embedded, and the average
surviving time of the seeds. By extensive computer simulations we aim the
reproduction of relevant statistical measures derived from the experimental
data of the Barro Colorado Island tree census in year 1995. The first two
parameters of the model are fixed to known values, characteristic of the
studied community, thus obtaining a model with only one freely adjustable
parameter. As a result of this, the average number of species in the considered
territory, the relative species abundance distribution, the species-area
relationship and the spatial auto-correlation function of the individuals in
abundant species are simultaneously fitted with only one parameter which is the
average surviving time of the seeds.


  Human brain contains about 10 billion neurons, each of which has about
10~10,000 nerve endings from which neurotransmitters are released in response
to incoming spikes, and the released neurotransmitters then bind to receptors
located in the postsynaptic neurons. However, individually, neurons are noisy
and synaptic release is in general unreliable. But groups of neurons that are
arranged in specialized modules can collectively perform complex information
processing tasks robustly and reliably. How functionally groups of neurons
perform behavioural related tasks crucial rely on a coherent organization of
dynamics from membrane ionic kinetics to synaptic coupling of the network and
dynamics of rhythmic oscillations that are tightly linked to behavioural state.
  To capture essential features of the biological system at multiple
spatial-temporal scales, it is important to construct a suitable computational
model that is closely or solely based on experimental data. Depending on what
one wants to understand, these models can either be very functional and
biologically realistic descriptions with thousands of coupled differential
equations (Hodgkin-Huxley type) or greatly simplified caricatures
(integrate-and-fire type) which are useful for studying large interconnected
networks.


  Co-localization of networks of genes in the nucleus is thought to play an
important role in determining gene expression patterns. Based upon experimental
data, we built a dynamical model to test whether pure diffusion could account
for the observed co-localization of genes within a defined subnuclear region. A
simple standard Brownian motion model in two and three dimensions shows that
preferential co-localization is possible for co-regulated genes without any
direct interaction, and suggests the occurrence may be due to a limitation in
the number of available transcription factors. Experimental data of chromatin
movements demonstrates that fractional rather than standard Brownian motion is
more appropriate to model gene mobilizations, and we tested our dynamical model
against recent static experimental data, using a sub-diffusion process by which
the genes tend to colocalize more easily. Moreover, in order to compare our
model with recently obtained experimental data, we studied the association
level between genes and factors, and presented data supporting the validation
of this dynamic model. As further applications of our model, we applied it to
test against more biological observations. We found that increasing
transcription factor number, rather than factory number and nucleus size, might
be the reason for decreasing gene co-localization. In the scenario of
frequency- or amplitude-modulation of transcription factors, our model
predicted that frequency-modulation may increase the co-localization between
its targeted genes.


  Next-generation sequencing techniques have facilitated a large scale analysis
of human genetic variation. Despite the advances in sequencing speeds, the
computational discovery of structural variants is not yet standard. It is
likely that many variants have remained undiscovered in most sequenced
individuals. Here we present a novel internal segment size based approach,
which organizes all, including also concordant reads into a read alignment
graph where max-cliques represent maximal contradiction-free groups of
alignments. A specifically engineered algorithm then enumerates all max-cliques
and statistically evaluates them for their potential to reflect insertions or
deletions (indels). For the first time in the literature, we compare a large
range of state-of-the-art approaches using simulated Illumina reads from a
fully annotated genome and present various relevant performance statistics. We
achieve superior performance rates in particular on indels of sizes 20--100,
which have been exposed as a current major challenge in the SV discovery
literature and where prior insert size based approaches have limitations. In
that size range, we outperform even split read aligners. We achieve good
results also on real data where we make a substantial amount of correct
predictions as the only tool, which complement the predictions of split-read
aligners. CLEVER is open source (GPL) and available from
http://clever-sv.googlecode.com.


  We measured mobility patterns that describe walking trajectories of
individual Me'Phaa peasants searching and collecting fuelwood in the forests of
"La Monta\~na de Guerrero" in Mexico. These one-day excursions typically follow
a mixed pattern of nearly-constant steps when individuals displace from their
homes towards potential collecting sites and a mixed pattern of steps of
different lengths when actually searching for fallen wood in the forest.
Displacements in the searching phase seem not to be compatible with L\'evy
flights described by power-laws with optimal scaling exponents. These findings
however can be interpreted in the light of deterministic searching on heavily
degraded landscapes where the interaction of the individuals with their scarce
environment produces alternative searching strategies than the expected L\'evy
flights. These results have important implications for future management and
restoration of degraded forests and the improvement of the ecological services
they may provide to their inhabitants.


  A Potts associative memory network has been proposed as a simplified model of
macroscopic cortical dynamics, in which each Potts unit stands for a patch of
cortex, which can be activated in one of S local attractor states. The internal
neuronal dynamics of the patch is not described by the model, rather it is
subsumed into an effective description in terms of graded Potts units, with
adaptation effects both specific to each attractor state and generic to the
patch. If each unit, or patch, receives effective (tensor) connections from C
other units, the network has been shown to be able to store a large number p of
global patterns, or network attractors, each with a fraction a of the units
active, where the critical load p_c scales roughly like p_c ~ (C S^2)/(a
ln(1/a)) (if the patterns are randomly correlated). Interestingly, after
retrieving an externally cued attractor, the network can continue jumping, or
latching, from attractor to attractor, driven by adaptation effects. The
occurrence and duration of latching dynamics is found through simulations to
depend critically on the strength of local attractor states, expressed in the
Potts model by a parameter w. Here we describe with simulations and then
analytically the boundaries between distinct phases of no latching, of
transient and sustained latching, deriving a phase diagram in the plane w-T,
where T parametrizes thermal noise effects. Implications for real cortical
dynamics are briefly reviewed in the conclusions.


  Many biological oscillators are arranged in networks composed of many
inter-coupled cellular oscillators. However, results are still lacking on the
collective oscillation period of inter-coupled gene regulatory oscillators,
which, as has been reported, may be different from the oscillation period of an
autonomous cellular oscillator. Based on the Goodwin oscillator, we analyze the
collective oscillation pattern of coupled cellular oscillator networks. First
we give a condition under which the oscillator network exhibits oscillatory and
synchronized behavior, then we estimate the collective oscillation period based
on a multivariable harmonic balance technique. Analytical results are derived
in terms of biochemical parameters, thus giving insight into the basic
mechanism of biological oscillation and providing guidance in synthetic biology
design. Simulation results are given to confirm the theoretical predictions.


  This paper is concerned with the evolution of haploid organisms that
reproduce asexually. In a seminal piece of work, Eigen and coauthors proposed
the quasispecies model in an attempt to understand such an evolutionary
process. Their work has impacted antiviral treatment and vaccine design
strategies. Yet, predictions of the quasispecies model are at best viewed as a
guideline, primarily because it assumes an infinite population size, whereas
realistic population sizes can be quite small. In this paper we consider a
population genetics-based model aimed at understanding the evolution of such
organisms with finite population sizes and present a rigorous study of the
convergence and computational issues that arise therein. Our first result is
structural and shows that, at any time during the evolution, as the population
size tends to infinity, the distribution of genomes predicted by our model
converges to that predicted by the quasispecies model. This justifies the
continued use of the quasispecies model to derive guidelines for intervention.
While the stationary state in the quasispecies model is readily obtained, due
to the explosion of the state space in our model, exact computations are
prohibitive. Our second set of results are computational in nature and address
this issue. We derive conditions on the parameters of evolution under which our
stochastic model mixes rapidly. Further, for a class of widely used fitness
landscapes we give a fast deterministic algorithm which computes the stationary
distribution of our model. These computational tools are expected to serve as a
framework for the modeling of strategies for the deployment of mutagenic drugs.


  The simple circadian oscillator found in cyanobacteria can be reconstituted
in vitro using three proteins-KaiA, KaiB and KaiC. The total phosphorylation
level of KaiC oscillates with a circadian period. Recent experiments elucidated
the structure and function of the cyanobacterial circadian oscillator, which is
driven by sunlight intensity variation and therefore by Earth's rotation. In
this paper, we derive the oscillator's equation of motion directly from
experimental data, couple it with Earth's rotation and compute its natural
periods and resonance curve. Finally, we propose the use of a quantum circuit
that is capable of describing the function of the oscillator.


  Constraints in embryonic development are thought to bias the direction of
evolution by making some changes less likely, and others more likely, depending
on their consequences on ontogeny. Here, we characterize the constraints acting
on genome evolution in vertebrates. We used gene expression data from two
vertebrates: zebrafish, using a microarray experiment spanning 14 stages of
development, and mouse, using EST counts for 26 stages of development. We show
that, in both species, genes expressed early in development (1) have a more
dramatic effect of knock-out or mutation and (2) are more likely to revert to
single copy after whole genome duplication, relative to genes expressed late.
This supports high constraints on early stages of vertebrate development,
making them less open to innovations (gene gain or gene loss). Results are
robust to different sources of data-gene expression from microarrays, ESTs, or
in situ hybridizations; and mutants from directed KO, transgenic insertions,
point mutations, or morpholinos. We determine the pattern of these constraints,
which differs from the model used to describe vertebrate morphological
conservation ("hourglass" model). While morphological constraints reach a
maximum at mid-development (the "phylotypic" stage), genomic constraints appear
to decrease in a monotonous manner over developmental time.


  Fixed points are fundamental states in any dynamical system. In the case of
gene regulatory networks (GRNs) they correspond to stable genes profiles
associated to the various cell types. We use Kauffman's approach to model GRNs
with random Boolean networks (RBNs). We start this paper by proving that, if we
fix the values of the source nodes (nodes with in-degree 0), the expected
number of fixed points of any RBN is one (independently of the topology we
choose). For finding such fixed points we use the {\alpha}-asynchronous
dynamics (where every node is updated independently with probability 0 <
{\alpha} < 1). In fact, it is well-known that asynchrony avoids the cycle
attractors into which parallel dynamics tends to fall. We perform simulations
and we show the remarkable property that, if for a given RBN with scale-free
topology and {\alpha}-asynchronous dynamics an initial configuration reaches a
fixed point, then every configuration also reaches a fixed point. By contrast,
in the parallel regime, the percentage of initial configurations reaching a
fixed point (for the same networks) is dramatically smaller. We contrast the
results of the simulations on scale-free networks with the classical
Erdos-Renyi model of random networks. Everything indicates that scale-free
networks are extremely robust. Finally, we study the mean and maximum time/work
needed to reach a fixed point when starting from randomly chosen initial
configurations.


  In this working paper, we present a simple theoretical framework based on
network theory to study how speciation, the process by which new species
appear, shapes spatial patterns of diversity. We show that this framework can
be expanded to account for different types of networks and interactions, and
incorporates different modes of speciation.


  Direct gating of mechanoelectrical-transduction channels by mechanical force
is a basic feature of hair cells that assures fast transduction and underpins
the mechanical amplification of acoustic inputs. But the associated
nonlinearity - the gating compliance - inevitably distorts signals. Because
reducing distortion would make the ear a better detector, we sought mechanisms
with that effect. Mimicking in vivo stimulation, we used stiff probes to
displace individual hair bundles at physiological amplitudes and measured the
coherence and phase of the relative stereociliary motions with a dual-beam
differential interferometer. Although stereocilia moved coherently and in phase
at the stimulus frequencies, large phase lags at the frequencies of the
internally generated distortion products indicated dissipative relative
motions. Tip links engaged these relative modes and decreased the coherence in
both stimulated and free hair bundles. These results show that a hair bundle
breaks into a highly dissipative serial arrangement of stereocilia at
distortion frequencies, precluding their amplification.


  Recent work has shown that the resistive force arising from viscous effects
within the pore region could explain [Ghosal, S. Phys. Rev. E. vol. 71, 051904
(2006) and Phys. Rev. Lett. vol. 98, 238104 (2007)] observed translocation
times in certain experiments involving voltage driven translocations of DNA
through nanopores. The electrokinetic flow inside the pore and the accompanying
viscous effects also play a crucial role in the interpretation of experiments
where the DNA is immobilized inside a nanopore[U. Keyser, et al., Nature
Physics vol. 2, 473 (2006)]. In this paper the viscous force is explicitly
calculated for a nanopore of cylindrical geometry. It is found that the
reduction of the tether force due to viscous drag and due to charge reduction
by Manning condensation are of similar size. The result is of importance in the
interpretation of experimental data on tethered DNA.


  In this research, we consider a mixture of genome fragments of a certain
bacteria set. The problem of mixture separation is studied under the assumption
that all the genomes present in the mixture are completely sequenced or are
close to those already sequenced. Such assumption is relevant, e.g., in regular
observations of ecological or biomedical objects, where the possible set of
microorganisms is known and it is only necessary to follow their
concentrations.


  In this paper, we present a technique for parameterizing Leslie transition
matrices from simple age and sex population counts, using an implementation of
"Wood's Method" [wood]; these matrices can forecast population by age and sex
(the "cohort component" method) using simple matrix multiplication and a
starting population. Our approach improves on previous methods for creating
Leslie matrices in two respects: it eliminates the need to calculate input
demographic rates from "raw" data, and our new format for the Leslie matrix
more elegantly reveals the population's demographic components of change
(fertility, mortality, and migration). The paper is organized around three main
themes. First, we describe the underlying algorithm, "Wood's Method," which
uses quadratic optimization to fit a transition matrix to age and sex
population counts. Second, we use demographic theory to create constraint sets
that make the algorithm useable for human populations. Finally, we use the
method to forecast 3,120 US counties and show that it holds promise for
automating cohort-component forecasts. This paper describes the first published
successful application of Wood's method to human populations; it also points to
more general promise of constrained optimization techniques in demographic
modeling.


  We consider a continuous time stochastic individual based model for a
population structured only by an inherited vector trait and with logistic
interactions. We consider its limit in a context from adaptive dynamics: the
population is large, the mutations are rare and we view the process in the
timescale of mutations. Using averaging techniques due to Kurtz (1992), we give
a new proof of the convergence of the individual based model to the trait
substitution sequence of Metz et al. (1992) first worked out by Dieckman and
Law (1996) and rigorously proved by Champagnat (2006): rigging the model such
that "invasion implies substitution", we obtain in the limit a process that
jumps from one population equilibrium to another when mutations occur and
invade the population.


  In an effort to design efficient platform for siRNA delivery, we combine all
atom classical and quantum simulations to study the binding of small
interfering RNA (siRNA) by pristine single wall carbon nanotube (SWCNT). Our
results show that siRNA strongly binds to SWCNT surface via unzipping its
base-pairs and the propensity of unzipping increases with the increase in the
diameter of the SWCNTs. The unzipping and subsequent wrapping events are
initiated and driven by van der Waals interactions between the aromatic rings
of siRNA nucleobases and the SWCNT surface. However, MD simulations of double
strand DNA (dsDNA) of the same sequence show that the dsDNA undergoes much less
unzipping and wrapping on the SWCNT in the simulation time scale of 70 ns. This
interesting difference is due to smaller interaction energy of thymidine of
dsDNA with the SWCNT compared to that of uridine of siRNA, as calculated by
dispersion corrected density functional theory (DFT) methods. After the optimal
binding of siRNA to SWCNT, the complex is very stable which serves as one of
the major mechanisms of siRNA delivery for biomedical applications. Since siRNA
has to undergo unwinding process with the effect of RNA- induced silencing
complex, our proposed delivery mechanism by SWCNT possesses potential
advantages in achieving RNA interference (RNAi).


  Genome-wide protein-protein interaction (PPI) data are readily available
thanks to recent breakthroughs in biotechnology. However, PPI networks of
extant organisms are only snapshots of the network evolution. How to infer the
whole evolution history becomes a challenging problem in computational biology.
In this paper, we present a likelihood-based approach to inferring network
evolution history from the topology of PPI networks and the duplication
relationship among the paralogs. Simulations show that our approach outperforms
the existing ones in terms of the accuracy of reconstruction. Moreover, the
growth parameters of several real PPI networks estimated by our method are more
consistent with the ones predicted in literature.


  The path-difference metric is one of the oldest and most popular distances
for the comparison of phylogenetic trees, but its statistical properties are
still quite unknown. In this paper we compute the expected value under the Yule
model of evolution of its square on the space of fully resolved rooted
phylogenetic trees with n leaves. This complements previous work by Steel-Penny
and Mir-Rossell\'o, who computed this mean value for fully resolved unrooted
and rooted phylogenetic trees, respectively, under the uniform distribution.


  Molecular machines consist of either a single protein or a macromolecular
complex composed of protein and RNA molecules. Just like their macroscopic
counterparts, each of these nano-machines has an engine that "transduces" input
energy into an output form which is then utilized by its coupling to a
transmission system for appropriate operations. The theory of heat engines,
pioneered by Carnot, rests on the second law of equilibrium thermodynamics.
However, the engines of molecular machines, operate under isothermal conditions
far from thermodynamic equilibrium. Moreover, one of the possible mechanisms of
energy transduction, popularized by Feynman and called Brownian ratchet, does
not even have any macroscopic counterpart. But, {\it molecular machine is not
synonymous with Brownian ratchet}; a large number of molecular machines
actually execute a noisy power stroke, rather than operating as Brownian
ratchet. The man-machine analogy, a topic of intense philosophical debate in
which many leading philosophers like Aristotle and Descartes participated, was
extended to similar analogies at the cellular and subcellular levels after the
invention of optical microscope. The idea of molecular machine, pioneered by
Marcelo Malpighi, has been pursued vigorously in the last fifty years. It has
become a well established topic of current interdisciplinary research as
evident from the publication of a very influential paper by Alberts towards the
end of the twentieth century. Here we give a non-technical overview of the
strategies for (a) stochastic modeling of mechano-chemical kinetic processes,
and (b) model selection based on statistical inference drawn from analysis of
experimental data. It is written for non-experts and from a broad perspective,
showing overlapping concepts from several different branches of physics and
from other areas of science and technology.


  A stochastic model for the dynamics of enzymatic catalysis in explicit,
effective solvents under physiological conditions is presented.
Analytically-computed first passage time densities of a diffusing particle in a
spherical shell with absorbing boundaries are combined with densities obtained
from explicit simulation to obtain the overall probability density for the
total reaction cycle time of the enzymatic system. The method is used to
investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate
kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The
direct simulation of the enzyme-substrate binding and reaction is carried out
using an elastic network model for the protein, and the solvent motions are
described by multiparticle collision dynamics, which incorporates hydrodynamic
flow effects. Systems where solvent-enzyme coupling occurs through explicit
intermolecular interactions, as well as systems where this coupling is taken
into account by including the protein and substrate in the multiparticle
collision step, are investigated and compared with simulations where
hydrodynamic coupling is absent. It is demonstrated that the flow of solvent
particles around the enzyme facilitates the large-scale hinge motion of the
enzyme with bound substrates, and has a significant impact on the shape of the
probability densities and average time scales of substrate binding for
substrates near the enzyme, the closure of the enzyme after binding, and the
overall time of completion of the cycle.


  Atomistic molecular dynamics simulations have been carried out to reveal the
characteristic features of ethylenediamine (EDA) cored protonated poly amido
amine (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) that are
functionalized with single stranded DNAs (ssDNAs). The four ssDNA strands that
are attached via alkythiolate [-S (CH2)6-] linker molecule to the free amine
groups on the surface of the PAMAM dendrimers observed to undergo a rapid
conformational change during the 25 ns long simulation period. From the RMSD
values of ssDNAs, we find relative stability in the case of purine rich ssDNA
strands than pyrimidine rich ssDNA strands. The degree of wrapping of ssDNA
strands on the dendrimer molecule was found to be influenced by the charge
ratio of DNA and the dendrimer. As G4 dendrimer contains relatively more
positive charge than G3 dendrimer, we observe extensive wrapping of ssDNAs on
the G4 dendrimer. The ssDNA strands along with the linkers are seen to
penetrate the surface of the dendrimer molecule and approach closer to the
center of the dendrimer indicating the soft sphere nature of the dendrimer
molecule. The effective radius of DNA-functionalized dendrimer nanoparticle was
found to be independent of base composition of ssDNAs and was observed to be
around 19.5 {\AA} and 22.4 {\AA} when we used G3 and G4 PAMAM dendrimer as the
core of the nanoparticle respectively. The observed effective radius of
DNA-functionalized dendrimer molecule apparently indicates the significant
shrinkage in the structure that has taken place in dendrimer, linker and DNA
strands. As a whole our results describe the characteristic features of
DNA-functionalized dendrimer nanoparticle and can be used as strong inputs to
design effectively the DNA-dendrimer nanoparticle self-assembly for their
active biological applications.


  Background: Identification of causal SNPs in most genome wide association
studies relies on approaches that consider each SNP individually. However,
there is a strong correlation structure among SNPs that need to be taken into
account. Hence, increasingly modern computationally expensive regression
methods are employed for SNP selection that consider all markers simultaneously
and thus incorporate dependencies among SNPs.
  Results: We develop a novel multivariate algorithm for large scale SNP
selection using CAR score regression, a promising new approach for prioritizing
biomarkers. Specifically, we propose a computationally efficient procedure for
shrinkage estimation of CAR scores from high-dimensional data. Subsequently, we
conduct a comprehensive comparison study including five advanced regression
approaches (boosting, lasso, NEG, MCP, and CAR score) and a univariate approach
(marginal correlation) to determine the effectiveness in finding true causal
SNPs.
  Conclusions: Simultaneous SNP selection is a challenging task. We demonstrate
that our CAR score-based algorithm consistently outperforms all competing
approaches, both uni- and multivariate, in terms of correctly recovered causal
SNPs and SNP ranking. An R package implementing the approach as well as R code
to reproduce the complete study presented here is available from
http://strimmerlab.org/software/care/ .


  One of the important questions in biological evolution is to know if certain
changes along protein coding genes have contributed to the adaptation of
species. This problem is known to be biologically complex and computationally
very expensive. It, therefore, requires efficient Grid or cluster solutions to
overcome the computational challenge. We have developed a Grid-enabled tool
(gcodeml) that relies on the PAML (codeml) package to help analyse large
phylogenetic datasets on both Grids and computational clusters. Although we
report on results for gcodeml, our approach is applicable and customisable to
related problems in biology or other scientific domains.


  We propose a new toy model of a heteropolymer chain capable of forming planar
secondary structures typical for RNA molecules. In this model the sequential
intervals between neighboring monomers along a chain are considered as quenched
random variables. Using the optimization procedure for a special class of
concave--type potentials, borrowed from optimal transport analysis, we derive
the local difference equation for the ground state free energy of the chain
with the planar (RNA--like) architecture of paired links. We consider various
distribution functions of intervals between neighboring monomers (truncated
Gaussian and scale--free) and demonstrate the existence of a topological
crossover from sequential to essentially embedded (nested) configurations of
paired links.


  We apply optimal control theory to a tuberculosis model given by a system of
ordinary differential equations. Optimal control strategies are proposed to
minimize the cost of interventions. Numerical simulations are given using data
from Angola.


  We study a character-based phylogeny reconstruction problem when an
incomplete set of data is given. More specifically, we consider the situation
under the directed perfect phylogeny assumption with binary characters in which
for some species the states of some characters are missing. Our main object is
to give an efficient algorithm to enumerate (or list) all perfect phylogenies
that can be obtained when the missing entries are completed. While a simple
branch-and-bound algorithm (B&B) shows a theoretically good performance, we
propose another approach based on a zero-suppressed binary decision diagram
(ZDD). Experimental results on randomly generated data exhibit that the ZDD
approach outperforms B&B. We also prove that counting the number of
phylogenetic trees consistent with a given data is #P-complete, thus providing
an evidence that an efficient random sampling seems hard.


  Learning how proteins fold will hardly have any impact in the way
conventional -- active site centered -- drugs are designed. On the other hand,
this knowledge is proving instrumental in defining a new paradigm for the
identification of drugs against any target protein: folding inhibition.
Targeting folding renders drugs less prone to elicit spontaneous genetic
mutations which in many cases, notably in connection with viruses like the
Human Immunodeficiency Virus (HIV), can block therapeutic action. From the
progress which has taken place during the last years in the understanding of
the becoming of a protein, and how to read from the corresponding sequences the
associated three-dimensional, biologically active, native structure, the idea
of non-conventional (folding) inhibitors and thus of leads to eventual drugs to
fight disease, arguably, without creating resistance, emerges as a distinct
possibility.


  Contact tracing data collected from disease outbreaks has received relatively
little attention in the epidemic modelling literature because it is thought to
be unreliable: infection sources might be wrongly attributed, or data might be
missing due to resource contraints in the questionnaire exercise. Nevertheless,
these data might provide a rich source of information on disease transmission
rate. This paper presents novel methodology for combining contact tracing data
with rate-based contact network data to improve posterior precision, and
therefore predictive accuracy. We present an advancement in Bayesian inference
for epidemics that assimilates these data, and is robust to partial contact
tracing. Using a simulation study based on the British poultry industry, we
show how the presence of contact tracing data improves posterior predictive
accuracy, and can directly inform a more effective control strategy.


  We present a general framework to describe the evolutionary dynamics of an
arbitrary number of types in finite populations based on stochastic
differential equations (SDE). For large, but finite populations this allows to
include demographic noise without requiring explicit simulations. Instead, the
population size only rescales the amplitude of the noise. Moreover, this
framework admits the inclusion of mutations between different types, provided
that mutation rates, $\mu$, are not too small compared to the inverse
population size 1/N. This ensures that all types are almost always represented
in the population and that the occasional extinction of one type does not
result in an extended absence of that type. For $\mu N\ll1$ this limits the use
of SDE's, but in this case there are well established alternative
approximations based on time scale separation. We illustrate our approach by a
Rock-Scissors-Paper game with mutations, where we demonstrate excellent
agreement with simulation based results for sufficiently large populations. In
the absence of mutations the excellent agreement extends to small population
sizes.


  Biochemical reaction systems may be viewed as discrete event processes
characterized by a number of states and state transitions. These systems may be
modeled as state transition systems with transitions representing individual
reaction events. Since they often involve a large number of interactions, it
can be difficult to construct such a model for a system, and since the
resulting state-level model can involve a huge number of states, model analysis
can be difficult or impossible. Here, we describe methods for the high-level
specification of a system using hypergraphs, for the automated generation of a
state-level model from a high-level model, and for the exact reduction of a
state-level model using information from the high-level model. Exact reduction
is achieved through the automated application of symmetry reduction and
invariant manifold reduction techniques to the high-level model, allowing
potentially significant reductions without the need to generate a full model.
The application of the method to biochemical reaction systems is illustrated by
models describing a hypothetical ion-channel at several levels of complexity.
The method allows for the reduction of the otherwise intractable example models
to a manageable size.


  The Luria-Delbr\"uck distribution is a classical model of mutations in cell
kinetics. It is obtained as a limit when the probability of mutation tends to
zero and the number of divisions to infinity. It can be interpreted as a
compound Poisson distribution (for the number of mutations) of exponential
mixtures (for the developing time of mutant clones) of geometric distributions
(for the number of cells produced by a mutant clone in a given time). The
probabilistic interpretation, and a rigourous proof of convergence in the
general case, are deduced from classical results on Bellman-Harris branching
processes. The two parameters of the Luria-Delbr\"uck distribution are the
expected number of mutations, which is the parameter of interest, and the
relative fitness of normal cells compared to mutants, which is the heavy tail
exponent. Both can be simultaneously estimated by the maximum likehood method.
However, the computation becomes numerically unstable as soon as the maximal
value of the sample is large, which occurs frequently due to the heavy tail
property. Based on the empirical generating function, robust estimators are
proposed and their asymptotic variance is given. They are comparable in
precision to maximum likelihood estimators, with a much broader range of
calculability, a better numerical stability, and a negligible computing time.


  The study of biological cells in terms of mesoscopic, nonequilibrium,
nonlinear, stochastic dynamics of open chemical systems provides a paradigm for
other complex, self-organizing systems with ultra-fast stochastic fluctuations,
short-time deterministic nonlinear dynamics, and long-time evolutionary
behavior with exponentially distributed rare events, discrete jumps among
punctuated equilibria, and catastrophe.


  Speciation is the "elephant in the room" of community ecology. As the
ultimate source of biodiversity, its integration in ecology's theoretical
corpus is necessary to understand community assembly. Yet, speciation is often
completely ignored or stripped of its spatial dimension. Recent approaches
based on network theory have allowed ecologists to effectively model complex
landscapes. In this study, we use this framework to model allopatric and
parapatric speciation in networks of communities and focus on the relationship
between speciation, richness, and the spatial structure of communities. We find
a strong opposition between speciation and local richness, with speciation
being more common in isolated communities and local richness being higher in
more connected communities. Unlike previous models, we also find a transition
to a positive relationship between speciation and local richness when dispersal
is low and the number of communities is small. Also, we use several measures of
centrality to characterize the effect of network structure on diversity. The
degree, the simplest measure of centrality, is found to be the best predictor
of local richness and speciation, although it loses some of its predictive
power as connectivity grows. Our framework shows how a simple neutral model can
be combined with network theory to reveal complex relationships between
speciation, richness, and the spatial organization of populations.


  We describe the distribution of Guanine and Cytosine (GC) content in the
third codon position (GC3) distributions in different species, analyze
evolutionary trends and discuss differences between genes and organisms with
distinct GC3 levels. We scrutinize previously published theoretical frameworks
and construct a unified view of GC3 biology in eukaryotes and prokaryotes.


  It has recently been shown that the NP-hard problem of calculating the
minimum number of hybridization events that is needed to explain a set of
rooted binary phylogenetic trees by means of a hybridization network is
fixed-parameter tractable if an instance of the problem consists of precisely
two such trees. In this paper, we show that this problem remains
fixed-parameter tractable for an arbitrarily large set of rooted binary
phylogenetic trees. In particular, we present a quadratic kernel.


  The prion-forming C-terminal domain of the fungal prion HET-s forms
infectious amyloid fibrils at physiological pH. The conformational switch from
the non-prion soluble form to the prion fibrillar form is believed to have a
functional role, since HET-s in its prion form participates in a recognition
process of different fungal strains. Based on the knowledge of the
high-resolution structure of HET-s(218-289) (the prion forming-domain) in its
fibrillar form, we here present a numerical simulation of the fibril growth
process which emphasizes the role of the topological properties of the
fibrillar structure. An accurate thermodynamic analysis of the way an
intervening HET-s chain is recruited to the tip of the growing fibril suggests
that elongation proceeds through a dock and lock mechanism. First, the chain
docks onto the fibril by forming the longest $\beta$-strands. Then, the
re-arrangement in the fibrillar form of all the rest of molecule takes place.
Interestingly, we predict also that one side of the HET-s fibril is more
suitable for substaining its growth with respect to the other. The resulting
strong polarity of fibril growth is a consequence of the complex topology of
HET-s fibrillar structure, since the central loop of the intervening chain
plays a crucially different role in favouring or not the attachment of the
C-terminus tail to the fibril, depending on the growth side.


  Multistep denaturation in a short circular DNA molecule is analyzed by a
mesoscopic Hamiltonian model which accounts for the helicoidal geometry.
Computation of melting profiles by the path integral method suggests that
stacking anharmonicity stabilizes the double helix against thermal disruption
of the hydrogen bonds. Twisting is essential in the model to capture the
importance of nonlinear effects on the thermodynamical properties. In a ladder
model with zero twist, anharmonic stacking scarcely affects the thermodynamics.
Moderately untwisted helices, with respect to the equilibrium conformation,
show an energetic advantage against the overtwisted ones. Accordingly
moderately untwisted helices better sustain local fluctuational openings and
make more unlikely the thermally driven complete strand separation.


  An open-access program allowing three-item statement matrices to be generated
from data such as molecular sequences does not exist so far. The recently
developed LisBeth package (ver. 1.0) allows representing hypotheses of homology
among taxa or areas directly as rooted trees or as hierarchies; however, this
is not a standard matrix-based platform. Here we present "TAXON version 1.1"
(TAXON), a program designed for building three-item statement-matrices from
binary, additive (ordered) and non-additive (unordered) multistate characters,
with both fractional and uniform weighting of the resulted statements.


  Stochastic simulations of coarse-grained protein models are used to
investigate the propensity to form knots in early stages of protein folding.
The study is carried out comparatively for two homologous
carbamoyltransferases, a natively-knotted N-acetylornithine
carbamoyltransferase (AOTCase) and an unknotted ornithine carbamoyltransferase
(OTCase). In addition, two different sets of pairwise amino acid interactions
are considered: one promoting exclusively native interactions, and the other
additionally including non-native quasi-chemical and electrostatic
interactions. With the former model neither protein show a propensity to form
knots. With the additional non-native interactions, knotting propensity remains
negligible for the natively-unknotted OTCase while for AOTCase it is much
enhanced. Analysis of the trajectories suggests that the different entanglement
of the two transcarbamylases follows from the tendency of the C-terminal to
point away from (for OTCase) or approach and eventually thread (for AOTCase)
other regions of partly-folded protein. The analysis of the OTCase/AOTCase pair
clarifies that natively-knotted proteins can spontaneously knot during early
folding stages and that non-native sequence-dependent interactions are
important for promoting and disfavoring early knotting events.


  We present new theoretical and empirical results on the probability
distributions of species persistence times in natural ecosystems. Persistence
times, defined as the timespans occurring between species' colonization and
local extinction in a given geographic region, are empirically estimated from
local observations of species' presence/absence. A connected sampling problem
is presented, generalized and solved analytically. Species persistence is shown
to provide a direct connection with key spatial macroecological patterns like
species-area and endemics-area relationships. Our empirical analysis pertains
to two different ecosystems and taxa: a herbaceous plant community and a
estuarine fish database. Despite the substantial differences in ecological
interactions and spatial scales, we confirm earlier evidence on the general
properties of the scaling of persistence times, including the predicted effects
of the structure of the spatial interaction network. The framework tested here
allows to investigate directly nature and extent of spatial effects in the
context of ecosystem dynamics. The notable coherence between spatial and
temporal macroecological patterns, theoretically derived and empirically
verified, is suggested to underlie general features of the dynamic evolution of
ecosystems.


  One strategy for winning a coevolutionary struggle is to evolve rapidly. Most
of the literature on host-pathogen coevolution focuses on this phenomenon, and
looks for consequent evidence of coevolutionary arms races. An alternative
strategy, less often considered in the literature, is to deter rapid
evolutionary change by the opponent. To study how this can be done, we
construct an evolutionary game between a controller that must process
information, and an adversary that can tamper with this information processing.
In this game, a species can foil its antagonist by processing information in a
way that is hard for the antagonist to manipulate. We show that the structure
of the information processing system induces a fitness landscape on which the
adversary population evolves, and that complex processing logic is required to
make that landscape rugged. Drawing on the rich literature concerning rates of
evolution on rugged landscapes, we show how a species can slow adaptive
evolution in the adversary population. We suggest that this type of defensive
complexity on the part of the vertebrate adaptive immune system may be an
important element of coevolutionary dynamics between pathogens and their
vertebrate hosts.


  Many models for the origin of life have focused on understanding how
evolution can drive the refinement of a preexisting enzyme, such as the
evolution of efficient replicase activity. Here we present a model for what
was, arguably, an even earlier stage of chemical evolution, when polymer
sequence diversity was generated and sustained before, and during, the onset of
functional selection. The model includes regular environmental cycles (e.g.
hydration-dehydration cycles) that drive polymers between times of replication
and functional activity, which coincide with times of different monomer and
polymer diffusivity. Kinetic Monte Carlo simulations demonstrate that this
proposed prebiotic scenario provides a robust mechanism for the exploration of
sequence space. Introduction of a polymer sequence with monomer synthetase
activity illustrates that functional sequences can become established in a
preexisting pool of otherwise non-functional sequences. Functional selection
does not dominate system dynamics and sequence diversity remains high,
permitting the emergence and spread of more than one functional sequence. It is
also observed that polymers spontaneously form clusters in simulations where
polymers diffuse more slowly than monomers, a feature that is reminiscent of a
previous proposal that the earliest stages of life could have been defined by
the collective evolution of a system-wide cooperation of polymer aggregates.
Overall, the results presented demonstrate the merits of considering plausible
prebiotic polymer chemistries and environments that would have allowed for the
rapid turnover of monomer resources and for regularly varying monomer/polymer
diffusivities.


  The identification of cellular objectives is one of the central topics in the
research of microbial metabolic networks. In particular, the information about
a cellular objective is needed in flux balance analysis which is a commonly
used constrained-based metabolic network analysis method for the prediction of
cellular phenotypes. The cellular objective may vary depending on the organism
and its growth conditions. It is probable that nutritionally scarce conditions
are very common in the nature and, in order to survive in those conditions,
cells exhibit various highly efficient nutrient processing systems like
enzymes. In this study, we explore the efficiency of a metabolic network in
transformation of substrates to new biomass, and we introduce a new objective
function simulating growth efficiency.
  We examined the properties of growth efficiency using a metabolic model for
Eschericia coli. We found that the maximal growth efficiency is obtained at a
finite nutrient uptake rate. The rate is substrate-dependent and it typically
does not exceed 20 mmol/h/gDW. We further examined whether the maximal growth
efficiency could serve as a cellular objective function in metabolic network
analysis, and found that cellular growth in batch cultivation can be predicted
reasonably well under this assumption. The fit to experimental data was found
slightly better than with the commonly used objective function of maximal
growth rate.
  Based on our results, we suggest that the maximal growth efficiency can be
considered as a plausible optimization criterion in metabolic modeling for E.
coli. In the future, it would be interesting to study growth efficiency as a
cellular objective also in other cellular systems and under different
cultivation conditions.


  Deep shotgun sequencing and analysis of genomes, transcriptomes, amplified
single-cell genomes, and metagenomes has enabled investigation of a wide range
of organisms and ecosystems. However, sampling variation in short-read data
sets and high sequencing error rates of modern sequencers present many new
computational challenges in data interpretation. These challenges have led to
the development of new classes of mapping tools and {\em de novo} assemblers.
These algorithms are challenged by the continued improvement in sequencing
throughput. We here describe digital normalization, a single-pass computational
algorithm that systematizes coverage in shotgun sequencing data sets, thereby
decreasing sampling variation, discarding redundant data, and removing the
majority of errors. Digital normalization substantially reduces the size of
shotgun data sets and decreases the memory and time requirements for {\em de
novo} sequence assembly, all without significantly impacting content of the
generated contigs. We apply digital normalization to the assembly of microbial
genomic data, amplified single-cell genomic data, and transcriptomic data. Our
implementation is freely available for use and modification.


  The adaptive evolution of a population under the influence of mutation and
selection is strongly influenced by the structure of the underlying fitness
landscape, which encodes the interactions between mutations at different
genetic loci. Theoretical studies of such landscapes have been carried out for
several decades, but only recently experimental fitness measurements
encompassing all possible combinations of small sets of mutations have become
available. The empirical studies have spawned new questions about the
accessibility of optimal genotypes under natural selection. Depending on
population dynamic parameters such as mutation rate and population size,
evolutionary accessibility can be quantified through the statistics of
accessible mutational pathways (along which fitness increases monotoncially),
or through the study of the basin of attraction of the optimal genotype under
greedy (steepest ascent) dynamics. Here we investigate these two measures of
accessibility in the framework of Kauffman's $LK$-model, a paradigmatic family
of random fitness landscapes with tunable ruggedness. The key parameter
governing the strength of genetic interactions is the number $K$ of interaction
partners of each of the $L$ sites in the genotype sequence. In general,
accessibility increases with increasing genotype dimensionality $L$ and
decreases with increasing number of interactions $K$. Remarkably, however, we
find that some measures of accessibility behave non-monotonically as a function
of $K$, indicating a special role of the most sparsely connected, non-trivial
cases K=1 and 2. The relation between models for fitness landscapes and spin
glasses is also addressed.


  A large offspring number diploid biparental multilocus population model of
Moran type is our object of study. At each timestep, a pair of diploid
individuals drawn uniformly at random contribute offspring to the population.
The number of offspring can be large relative to the total population size.
Similar `heavily skewed' reproduction mechanisms have been considered by
various authors recently. Each diploid parental individual contributes exactly
one chromosome to each diploid offspring, and hence ancestral lineages can only
coalesce when in distinct individuals. A separation of timescales phenomenon is
thus observed. A result of M\"{o}hle (1998) is extended to obtain convergence
of the ancestral process to an ancestral recombination graph necessarily
admitting simultaneous multiple mergers of ancestral lineages. The usual
ancestral recombination graph is obtained as a special case of our model when
the parents contribute only one offspring to the population each time. Due to
diploidy and large offspring numbers, novel effects appear. For example, the
marginal genealogy at each locus admits simultaneous multiple mergers in up to
four groups, and different loci remain substantially correlated even as the
recombination rate grows large. Thus, genealogies for loci far apart on the
same chromosome remain correlated. Correlation in coalescence times for two
loci is derived and shown to be a function of the coalescence parameters of our
model. Extending the observations by Eldon and Wakeley (2008), predictions of
linkage disequilibrium are shown to be functions of the reproduction parameters
of our model, in addition to the recombination rate. Correlations in ratios of
coalescence times between loci can be high, even when the recombination rate is
high and sample size is large.


  Modern developments in microscopy and image processing are revolutionizing
areas of physics, chemistry and biology as nanoscale objects can be tracked
with unprecedented accuracy. The goal of single particle tracking is to
determine the interaction between the particle and its environment. The price
paid for having a direct visualization of a single particle is a consequent
lack of statistics. Here we address the optimal way of extracting diffusion
constants from single trajectories for pure Brownian motion. It is shown that
the maximum likelihood estimator is much more efficient than the commonly used
least squares estimate. Furthermore we investigate the effect of disorder on
the distribution of estimated diffusion constants and show that it increases
the probability of observing estimates much smaller than the true (average)
value.


  Using analytical techniques and Langevin dynamics simulations, we investigate
the dynamics of polymer translocation through a nanochannel embedded in two
dimensions under an applied external field. We examine the translocation time
for various ratio of the channel length $L$ to the polymer length $N$. For
short channels $L\ll N$, the translocation time $\tau \sim N^{1+\nu}$ under
weak driving force $F$, while $\tau\sim F^{-1}L$ for long channels $L\gg N$,
independent of the chain length $N$. Moreover, we observe a minimum of
translocation time as a function of $L/N$ for different driving forces and
channel widths. These results are interpreted by the waiting time of a single
segment.


  Mass-action kinetics is frequently used in systems biology to model the
behaviour of interacting chemical species. Many important dynamical properties
are known to hold for such systems if they are weakly reversible and have a low
deficiency. In particular, the Deficiency Zero and Deficiency One Theorems
guarantee strong regularity with regards to the number and stability of
positive equilibrium states. It is also known that chemical reaction networks
with disparate reaction structure can exhibit the same qualitative dynamics.
The theory of linear conjugacy encapsulates the cases where this relationship
is captured by a linear transformation. In this paper, we propose a
mixed-integer linear programming algorithm capable of determining weakly
reversible reaction networks with a minimal deficiency which are linearly
conjugate to a given reaction network.


  We present a new model for seed banks, where direct ancestors of individuals
may have lived in the near as well as the very far past. The classical
Wright-Fisher model, as well as a seed bank model with bounded age distribution
considered by Kaj, Krone and Lascoux (2001) are special cases of our model. We
discern three parameter regimes of the seed bank age distribution, which lead
to substantially different behaviour in terms of genetic variability, in
particular with respect to fixation of types and time to the most recent common
ancestor. We prove that for age distributions with finite mean, the ancestral
process converges to a time-changed Kingman coalescent, while in the case of
infinite mean, ancestral lineages might not merge at all with positive
probability. Further, we present a construction of the forward in time process
in equilibrium. The mathematical methods are based on renewal theory, the urn
process introduced by Kaj et al., as well as on a paper by Hammond and
Sheffield (2011).


  Cells often perform computations in response to environmental cues. A simple
example is the classic problem, first considered by Berg and Purcell, of
determining the concentration of a chemical ligand in the surrounding media. On
general theoretical grounds (Landuer's Principle), it is expected that such
computations require cells to consume energy. Here, we explicitly calculate the
energetic costs of computing ligand concentration for a simple two-component
cellular network that implements a noisy version of the Berg-Purcell strategy.
We show that learning about external concentrations necessitates the breaking
of detailed balance and consumption of energy, with greater learning requiring
more energy. Our calculations suggest that the energetic costs of cellular
computation may be an important constraint on networks designed to function in
resource poor environments such as the spore germination networks of bacteria.


  The branching process (BP) approach has been successful in explaining the
avalanche dynamics in complex networks. However, its applications are mainly
focused on unipartite networks, in which all nodes are of the same type. Here,
motivated by a need to understand avalanche dynamics in metabolic networks, we
extend the BP approach to a particular bipartite network composed of Boolean
AND and OR logic gates. We reduce the bipartite network into a unipartite
network by integrating out OR gates, and obtain the effective branching ratio
for the remaining AND gates. Then the standard BP approach is applied to the
reduced network, and the avalanche size distribution is obtained. We test the
BP results with simulations on the model networks and two microbial metabolic
networks, demonstrating the usefulness of the BP approach.


  Motivated by geometric origami and the straight skeleton construction, we
outline a map between spaces of phylogenetic trees and spaces of planar
polygons. The limitations of this map is studied through explicit examples,
culminating in proving a structural rigidity result.


  Neutral macroevolutionary models, such as the Yule model, give rise to a
probability distribution on the set of discrete rooted binary trees over a
given leaf set. Such models can provide a signal as to the approximate location
of the root when only the unrooted phylogenetic tree is known, and this signal
becomes relatively more significant as the number of leaves grows. In this
short note, we show that among models that treat all taxa equally, and are
sampling consistent (i.e. the distribution on trees is not affected by taxa yet
to be included), all such models, except one, convey some information as to the
location of the ancestral root in an unrooted tree.


  Microbes providing public goods are widespread in nature despite running the
risk of being exploited by free-riders. However, the precise ecological factors
supporting cooperation are still puzzling. Following recent experiments, we
consider the role of population growth and the repetitive fragmentation of
populations into new colonies mimicking simple microbial life-cycles.
Individual-based modeling reveals that demographic fluctuations, which lead to
a large variance in the composition of colonies, promote cooperation. Biased by
population dynamics these fluctuations result in two qualitatively distinct
regimes of robust cooperation under repetitive fragmentation into groups.
First, if the level of cooperation exceeds a threshold, cooperators will take
over the whole population. Second, cooperators can also emerge from a single
mutant leading to a robust coexistence between cooperators and free-riders. We
find frequency and size of population bottlenecks, and growth dynamics to be
the major ecological factors determining the regimes and thereby the
evolutionary pathway towards cooperation.


  The forming and melting of complementary base pairs in RNA duplexes are
conformational transitions required to accomplish a plethora of biological
functions. Yet the dynamic steps of these transitions have not been
quantitatively characterized at the molecular level. In this work, the base
opening process was first enforced by atomistic pulling simulations and then
analyzed with a novel reweighting scheme which allowed the free-energy profile
along any suitable reaction coordinate, e.g. solvation, to be reconstructed.
The systematic application of such approach to different base-pair combinations
provides a molecular motion picture of helix opening which is validated by
comparison with an extensive set of experimental observations and links them to
the enzyme-dependent unwinding mechanism. The RNA intrinsic dynamics disclosed
in this work could rationalize the directionality observed in RNA-processing
molecular machineries.


  Summary: MALDIquant is an R package providing a complete and modular analysis
pipeline for quantitative analysis of mass spectrometry data. MALDIquant is
specifically designed with application in clinical diagnostics in mind and
implements sophisticated routines for importing raw data, preprocessing,
non-linear peak alignment, and calibration. It also handles technical
replicates as well as spectra with unequal resolution.
  Availability: MALDIquant and its associated R packages readBrukerFlexData and
readMzXmlData are freely available from the R archive CRAN
(http://cran.r-project.org). The software is distributed under the GNU General
Public License (version 3 or later) and is accompanied by example files and
data. Additional documentation is available from
http://strimmerlab.org/software/maldiquant/.


  Robustness of biochemical systems has become one of the central questions in
systems biology although it is notoriously difficult to formally capture its
multifaceted nature. Maintenance of normal system function depends not only on
the stoichiometry of the underlying interrelated components, but also on a
multitude of kinetic parameters. Invariant flux ratios, obtained within flux
coupling analysis, as well as invariant complex ratios, derived within chemical
reaction network theory, can characterize robust properties of a system at
steady state. However, the existing formalisms for the description of these
invariants do not provide full characterization as they either only focus on
the flux-centric or the concentration-centric view. Here we develop a novel
mathematical framework which combines both views and thereby overcomes the
limitations of the classical methodologies. Our unified framework will be
helpful in analyzing biologically important system properties.


  Electrical signaling via voltage-gated ion channels depends upon the function
of a voltage sensor (VS), identified with the S1-S4 domain in voltage-gated K+
channels. Here we investigate some energetic aspects of the sliding-helix model
of the VS using simulations based on VS charges, linear dielectrics and
whole-body motion. Model electrostatics in voltage-clamped boundary conditions
are solved using a boundary element method. The statistical mechanical
consequences of the electrostatic configurational energy are computed to gain
insight into the sliding-helix mechanism and to predict experimentally measured
ensemble properties such as gating charge displaced by an applied voltage.
Those consequences and ensemble properties are investigated for two alternate
S4 configurations, \alpha- and 3(10)-helical. Both forms of VS are found to
have an inherent electrostatic stability. Maximal charge displacement is
limited by geometry, specifically the range of movement where S4 charges and
counter-charges overlap in the region of weak dielectric. Charge displacement
responds more steeply to voltage in the \alpha-helical than the 3(10)-helical
sensor. This difference is due to differences on the order of 0.1 eV in the
landscapes of electrostatic energy. As a step toward integrating these VS
models into a full-channel model, we include a hypothetical external load in
the Hamiltonian of the system and analyze the energetic in/output relation of
the VS.


  The Time To the Most Recent Common Ancestor (TMRCA) based on human
mitochondrial DNA (mtDNA) is estimated to be twice that based on the
non-recombining part of the Y chromosome (NRY). These TMRCAs have special
demographic implications because mtDNA is transmitted only from mother to
child, and NRY from father to son. Therefore, mtDNA reflects female history,
and NRY, male history. To investigate what caused the two-to-one female-male
TMRCA ratio in humans, we develop a forward-looking agent-based model (ABM)
with overlapping generations and individual life cycles. We implement two main
mating systems: polygynandry and polygyny with different degrees in between. In
each mating system, the male population can be either homogeneous or
heterogeneous. In the latter case, some males are `alphas' and others are
`betas', which reflects the extent to which they are favored by female mates. A
heterogeneous male population implies a competition among males with the
purpose of signaling as alphas. The introduction of a heterogeneous male
population is found to reduce by a factor 2 the probability of finding equal
female and male TMRCAs and shifts the distribution of the TMRCA ratio to higher
values. We find that high male-male competition is necessary to reproduce a
TMRCA ratio of 2: less than half the males can be alphas and betas can have at
most half the fitness of alphas. In addition, in the modes that maximize the
probability of having a TMRCA ratio between 1.5 and 2.5, the present generation
has 1.4 times as many female as male ancestors. We also tested the effect of
sex-biased migration and sex-specific death rates and found that these are
unlikely to explain alone the sex-biased TMRCA ratio observed in humans. Our
results support the view that we are descended from males who were successful
in a highly competitive context, while females were facing a much smaller
female-female competition.


  DNA sequencing is the basic workhorse of modern day biology and medicine.
Shotgun sequencing is the dominant technique used: many randomly located short
fragments called reads are extracted from the DNA sequence, and these reads are
assembled to reconstruct the original sequence. A basic question is: given a
sequencing technology and the statistics of the DNA sequence, what is the
minimum number of reads required for reliable reconstruction? This number
provides a fundamental limit to the performance of {\em any} assembly
algorithm. For a simple statistical model of the DNA sequence and the read
process, we show that the answer admits a critical phenomena in the asymptotic
limit of long DNA sequences: if the read length is below a threshold,
reconstruction is impossible no matter how many reads are observed, and if the
read length is above the threshold, having enough reads to cover the DNA
sequence is sufficient to reconstruct. The threshold is computed in terms of
the Renyi entropy rate of the DNA sequence. We also study the impact of noise
in the read process on the performance.


  Organisms often grow, migrate and compete in liquid environments, as well as
on solid surfaces. However, relatively little is known about what happens when
competing species are mixed and compressed by fluid turbulence. In these
lectures we review our recent work on population dynamics and population
genetics in compressible velocity fields of one and two dimensions. We discuss
why compressible turbulence is relevant for population dynamics in the ocean
and we consider cases both where the velocity field is turbulent and when it is
static. Furthermore, we investigate populations in terms of a continuos density
field and when the populations are treated via discrete particles. In the last
case we focus on the competition and fixation of one species compared to
another


  Motivation: Eugene Myers in his string graph paper (Myers, 2005) suggested
that in a string graph or equivalently a unitig graph, any path spells a valid
assembly. As a string/unitig graph also encodes every valid assembly of reads,
such a graph, provided that it can be constructed correctly, is in fact a
lossless representation of reads. In principle, every analysis based on
whole-genome shotgun sequencing (WGS) data, such as SNP and insertion/deletion
(INDEL) calling, can also be achieved with unitigs.
  Results: To explore the feasibility of using de novo assembly in the context
of resequencing, we developed a de novo assembler, fermi, that assembles
Illumina short reads into unitigs while preserving most of information of the
input reads. SNPs and INDELs can be called by mapping the unitigs against a
reference genome. By applying the method on 35-fold human resequencing data, we
showed that in comparison to the standard pipeline, our approach yields similar
accuracy for SNP calling and better results for INDEL calling. It has higher
sensitivity than other de novo assembly based methods for variant calling. Our
work suggests that variant calling with de novo assembly be a beneficial
complement to the standard variant calling pipeline for whole-genome
resequencing. In the methodological aspects, we proposed FMD-index for
forward-backward extension of DNA sequences, a fast algorithm for finding all
super-maximal exact matches and one-pass construction of unitigs from an
FMD-index.
  Availability: http://github.com/lh3/fermi
  Contact: hengli@broadinstitute.org


  Motivation: Most existing methods for DNA sequence analysis rely on accurate
sequences or genotypes. However, in applications of the next-generation
sequencing (NGS), accurate genotypes may not be easily obtained (e.g.
multi-sample low-coverage sequencing or somatic mutation discovery). These
applications press for the development of new methods for analyzing sequence
data with uncertainty.
  Results: We present a statistical framework for calling SNPs, discovering
somatic mutations, inferring population genetical parameters and performing
association tests directly based on sequencing data without explicit genotyping
or linkage-based imputation. On real data, we demonstrate that our method
achieves comparable accuracy to alternative methods for estimating site allele
count, for inferring allele frequency spectrum and for association mapping. We
also highlight the necessity of using symmetric datasets for finding somatic
mutations and confirm that for discovering rare events, mismapping is
frequently the leading source of errors.
  Availability: http://samtools.sourceforge.net.
  Contact: hengli@broadinstitute.org.


  Mechanical properties of DNA, in particular their stretch dependent extension
and their loop formation characteristics, have been recognized as an effective
probe for understanding the possible biochemical role played by them in a
living cell. Single stranded DNA (ssDNA), which, till recently was presumed to
be an simple flexible polymer continues to spring surprises. Synthetic ssDNA,
like polydA (polydeoxyadenosines) has revealed an intriguing force-extension
(FX) behavior exhibiting two plateaus, absent in polydT (polydeoxythymidines)
for example. Loop closing time in polydA had also been found to scale
exponentially with inverse temperature, unexpected from generic models of
homopolymers. Here we present a new model for polydA which incorporates both a
helix-coil transition and a over-stretching transition, accounting for the two
plateaus. Using transfer matrix calculation and Monte-Carlo simulation we show
that the model reproduces different sets of experimental observations,
quantitatively. It also predicts interesting reentrant behavior in the
temperature-extension characteristics of polydA, which is yet to be verified
experimentally.