  Actin networks in certain single-celled organisms exhibit a complex
pattern-forming dynamics that starts with the appearance of static spots of
actin on the cell cortex. Spots soon become mobile, executing persistent random
walks, and eventually give rise to traveling waves of actin. Here we describe a
possible physical mechanism for this distinctive set of dynamic
transformations, by equipping an excitable reaction-diffusion model with a
field describing the spatial orientation of its chief constituent (which we
consider to be actin). The interplay of anisotropic actin growth and spatial
inhibition drives a transformation at fixed parameter values from static spots
to moving spots to waves.


  We study the robustness and stability of the yeast cell regulatory network by
using a general inhomogeneous discrete model. We find that inhomogeneity, on
average, enhances the stability of the biggest attractor of the dynamics and
that the large size of the basin of attraction is robust against changes in the
parameters of inhomogeneity. We find that the most frequent orbit, which
represents the cell-cycle pathway, has a better biological meaning than the one
exhibited by the homogeneous model.


  Cellular function is widely believed to be organized in a modular fashion. On
all scales and at all levels of complexity, relatively independent sub-units
perform relatively independent sub-tasks of biological function. This
functional modularity must be reflected in the topology of molecular networks.
But how a functional module should be represented in an interaction network is
an open question. In protein-interaction networks (PIN), one can identify a
protein-complex as a module on a small scale, i.e. modules are understood as
densely linked, resp. interacting, groups of proteins, that are only sparsely
interacting with the rest of the network.
  In this contribution, we show that extrapolating this concept of cohesively
linked clusters of proteins as modules to the scale of the entire PIN
inevitable misses important and functionally relevant structure inherent in the
network. As an alternative, we introduce a novel way of decomposing a network
into functional roles and show that this represents network structure and
function more efficiently. This finding should have a profound impact on all
module assisted methods of protein function prediction and should shed new
light on how functional modules can be represented in molecular interaction
networks in general.


  Quantifying human group dynamics represents a unique challenge. Unlike
animals and other biological systems, humans form groups in both real (offline)
and virtual (online) spaces -- from potentially dangerous street gangs
populated mostly by disaffected male youths, through to the massive global
guilds in online role-playing games for which membership currently exceeds tens
of millions of people from all possible backgrounds, age-groups and genders. We
have compiled and analyzed data for these two seemingly unrelated offline and
online human activities, and have uncovered an unexpected quantitative link
between them. Although their overall dynamics differ visibly, we find that a
common team-based model can accurately reproduce the quantitative features of
each simply by adjusting the average tolerance level and attribute range for
each population. By contrast, we find no evidence to support a version of the
model based on like-seeking-like (i.e. kinship or `homophily').


  Natural selection and random drift are competing phenomena for explaining the
evolution of populations. Combining a highly fit mutant with a population
structure that improves the odds that the mutant spreads through the whole
population tips the balance in favor of natural selection. The probability that
the spread occurs, known as the fixation probability, depends heavily on how
the population is structured. Certain topologies, albeit highly artificially
contrived, have been shown to exist that favor fixation. We introduce a
randomized mechanism for network growth that is loosely inspired in some of
these topologies' key properties and demonstrate, through simulations, that it
is capable of giving rise to structured populations for which the fixation
probability significantly surpasses that of an unstructured population. This
discovery provides important support to the notion that natural selection can
be enhanced over random drift in naturally occurring population structures.


  Competition between random genetic drift and natural selection plays a
central role in evolution: Whereas non-beneficial mutations often prevail in
small populations by chance, mutations that sweep through large populations
typically confer a selective advantage. Here, however, we observe chance
effects during range expansions that dramatically alter the gene pool even in
large microbial populations. Initially well-mixed populations of two
fluorescently labeled strains of Escherichia coli develop well-defined,
sector-like regions with fractal boundaries in expanding colonies. The
formation of these regions is driven by random fluctuations that originate in a
thin band of pioneers at the expanding frontier. A comparison of bacterial and
yeast colonies (Saccharomyces cerevisiae) suggests that this large-scale
genetic sectoring is a generic phenomenon that may provide a detectable
footprint of past range expansions.


  The protective effect exerted by polyamines (Put and Spd) against cadmium
(Cd) stress was investigated in Brassica juncea plants. Treatment with CdCl2
(75 micro-Mole) resulted in a rise of Cd accumulation, a decrease of fresh and
dry weights in every plant organ, an increase of free polyamine content at limb
and stem levels as well as a decrease at root level. On the other hand, the
total conjugated polyamine levels in the stem tissues were unaffected by Cd. In
the leaf tissues, this metal caused a reduction of chlorophyll a content, a
rise of guaiacol peroxidase (GPOX) activity and an increase of malondialdehyde
(MDA), soluble glucide, proline and amino acid contents. Exogenous application,
by spraying, of putrescine (Put) and spermidine (Spd) to leaf tissues reduced
CdCl2-induced stress. These polyamines proved to exert a partial, though
significant, protection of the foliar fresh weight and to alleviate the
oxidative stress generated by Cd through reductions of MDA amounts and GPOX
(E.C.1.11.1.7) activity. The enhancement of chlorophyll a content in plants by
Put and those of Chl a and Chl b by Spd both constitute evidences of their
efficacy against the Cd2+-induced loss of pigments. Conversely to Put, Spd
caused a decrease of Cd content in leave tissues and a rise in the stems and
roots; these findings are in favour of a stimulation of Cd uptake by Spd. The
proline stimulation observed with Cd was reduced further to the spraying of Put
onto tissues, but the decrease induced by Spd was more limited. In the plants
treated with Cd, the amino acid contents in the leaves were unaffected by Put
and Spd spraying; on the other hand, Cd2+ disturbed polyamine levels (free and
acido-soluble conjugated-forms); we notice the rise of total free PAs and the
decrease of their conjugated-ones.


  We normalize the combinatorial Laplacian of a graph by the degree sum, look
at its eigenvalues as a probability distribution and then study its Shannon
entropy. Equivalently, we represent a graph with a quantum mechanical state and
study its von Neumann entropy. At the graph-theoretic level, this quantity may
be interpreted as a measure of regularity; it tends to be larger in relation to
the number of connected components, long paths and nontrivial symmetries. When
the set of vertices is asymptotically large, we prove that regular graphs and
the complete graph have equal entropy, and specifically it turns out to be
maximum. On the other hand, when the number of edges is fixed, graphs with
large cliques appear to minimize the entropy.


  BACKGROUND: The uncoupling protein (UCP) genes belong to the superfamily of
electron transport carriers of the mitochondrial inner membrane. Members of the
uncoupling protein family are involved in thermogenesis and determining the
functional evolution of UCP genes is important to understand the evolution of
thermo-regulation in vertebrates. RESULTS: Sequence similarity searches of
genome and scaffold data identified homologues of UCP in eutherians, teleosts
and the first squamates uncoupling proteins. Phylogenetic analysis was used to
characterize the family evolutionary history by identifying two duplications
early in vertebrate evolution and two losses in the avian lineage (excluding
duplications within a species, excluding the losses due to incompletely
sequenced taxa and excluding the losses and duplications inferred through
mismatch of species and gene trees). Estimates of synonymous and nonsynonymous
substitution rates (dN/dS) and more complex branch and site models suggest that
the duplication events were not associated with positive Darwinian selection
and that the UCP is constrained by strong purifying selection except for a
single site which has undergone positive Darwinian selection, demonstrating
that the UCP gene family must be highly conserved. CONCLUSION: We present a
phylogeny describing the evolutionary history of the UCP gene family and show
that the genes have evolved through duplications followed by purifying
selection except for a single site in the mitochondrial matrix between the 5th
and 6th alpha-helices which has undergone positive selection.


  The relaxation of initially straight semiflexible polymers has been discussed
mainly with respect to the longest relaxation time. The biologically relevant
non-equilibrium dynamics on shorter times is comparatively poorly understood,
partly because "initially straight" can be realized in manifold ways. Combining
Brownian dynamics simulations and systematic theory, we demonstrate how
different experimental preparations give rise to specific short-time and
universal long-time dynamics. We also discuss boundary effects and the onset of
the stretch--coil transition.


  The large-scale properties of chemical reaction systems, such as the
metabolism, can be studied with graph-based methods. To do this, one needs to
reduce the information -- lists of chemical reactions -- available in
databases. Even for the simplest type of graph representation, this reduction
can be done in several ways. We investigate different simple network
representations by testing how well they encode information about one
biologically important network structure -- network modularity (the propensity
for edges to be cluster into dense groups that are sparsely connected between
each other). To reach this goal, we design a model of reaction-systems where
network modularity can be controlled and measure how well the reduction to
simple graphs capture the modular structure of the model reaction system. We
find that the network types that best capture the modular structure of the
reaction system are substrate-product networks (where substrates are linked to
products of a reaction) and substance networks (with edges between all
substances participating in a reaction). Furthermore, we argue that the
proposed model for reaction systems with tunable clustering is a general
framework for studies of how reaction-systems are affected by modularity. To
this end, we investigate statistical properties of the model and find, among
other things, that it recreate correlations between degree and mass of the
molecules.


  We study the impact of DNA coiling on the search rate of proteins moving
along the DNA contour interspersed by three-dimensional (3D) bulk excursions.
When the DNA is coiled proteins performing short 3D hops along a DNA segment
can be captured by foreign DNA segments that have looped back close to the
original segment. These intersegmental jumps are shown to enhance the search
rate for a specific site on the DNA by lowering the tendency to resample
previously visited sites on the DNA. The model developed here offers a
quantitative description of recent single molecule experiments on facilitated
diffusion of restriction enzymes EcoRV.


  Rhythmic and sequential subdivision of the elongating vertebrate embryonic
body axis into morphological somites is controlled by an oscillating
multicellular genetic network termed the segmentation clock. This clock
operates in the presomitic mesoderm (PSM), generating dynamic stripe patterns
of oscillatory gene-expression across the field of PSM cells. How these spatial
patterns, the clock's collective period, and the underlying cellular-level
interactions are related is not understood. A theory encompassing temporal and
spatial domains of local and collective aspects of the system is essential to
tackle these questions. Our delayed coupling theory achieves this by
representing the PSM as an array of phase oscillators, combining four key
elements: a frequency profile of oscillators slowing across the PSM; coupling
between neighboring oscillators; delay in coupling; and a moving boundary
describing embryonic axis elongation. This theory predicts that the
segmentation clock's collective period depends on delayed coupling. We derive
an expression for pattern wavelength across the PSM and show how this can be
used to fit dynamic wildtype gene-expression patterns, revealing the
quantitative values of parameters controlling spatial and temporal organization
of the oscillators in the system. Our theory can be used to analyze
experimental perturbations, thereby identifying roles of genes involved in
segmentation.


  Gene expression is significantly stochastic making modeling of genetic
networks challenging. We present an approximation that allows the calculation
of not only the mean and variance but also the distribution of protein numbers.
We assume that proteins decay substantially slower than their mRNA and confirm
that many genes satisfy this relation using high-throughput data from budding
yeast. For a two-stage model of gene expression, with transcription and
translation as first-order reactions, we calculate the protein distribution for
all times greater than several mRNA lifetimes and thus qualitatively predict
the distribution of times for protein levels to first cross an arbitrary
threshold. If in addition the promoter fluctuates between inactive and active
states, we can find the steady-state protein distribution, which can be bimodal
if promoter fluctuations are slow. We show that our assumptions imply that
protein synthesis occurs in geometrically distributed bursts and allows mRNA to
be eliminated from a master equation description. In general, we find that
protein distributions are asymmetric and may be poorly characterized by their
mean and variance. Through maximum likelihood methods, our expressions should
therefore allow more quantitative comparisons with experimental data. More
generally, we introduce a technique to derive a simpler, effective dynamics for
a stochastic system by eliminating a fast variable.


  Chloroplast microsatellites are becoming increasingly popular markers for
population genetic studies in plants, but there has been little focus on their
potential for demographic inference. In this work the utility of chloroplast
microsatellites for the study of population expansions was explored. First, we
investigated the power of mismatch distribution analysis and the F(S) test with
coalescent simulations of different demographic scenarios. We then applied
these methods to empirical data obtained for the Canary Island pine (Pinus
canariensis). The results of the simulations showed that chloroplast
microsatellites are sensitive to sudden population growth. The power of the
F(S) test and accuracy of demographic parameter estimates, such as the time of
expansion, were reduced proportionally to the level of homoplasy within the
data. The analysis of Canary Island pine chloroplast microsatellite data
indicated population expansions for almost all sample localities. Demographic
expansions at the island level can be explained by the colonization of the
archipelago by the pine, while population expansions of different ages in
different localities within an island could be the result of local extinctions
and recolonization dynamics. Comparable mitochondrial DNA sequence data from a
parasite of P. canariensis, the weevil Brachyderes rugatus, supports this
scenario, suggesting a key role for volcanism in the evolution of pine forest
communities in the Canary Islands.


  Characterization of transient intermediate or transition states is crucial
for the description of biomolecular folding pathways, which is however
difficult in both experiments and computer simulations. Such transient states
are typically of low population in simulation samples. Even for simple systems
such as RNA hairpins, recently there are mounting debates over the existence of
multiple intermediate states. In this paper, we develop a computational
approach to explore the relatively low populated transition or intermediate
states in biomolecular folding pathways, based on a topological data analysis
tool, Mapper, with simulation data from large-scale distributed computing. The
method is inspired by the classical Morse theory in mathematics which
characterizes the topology of high dimensional shapes via some functional level
sets. In this paper we exploit a conditional density filter which enables us to
focus on the structures on pathways, followed by clustering analysis on its
level sets, which helps separate low populated intermediates from high
populated uninteresting structures. A successful application of this method is
given on a motivating example, a RNA hairpin with GCAA tetraloop, where we are
able to provide structural evidence from computer simulations on the multiple
intermediate states and exhibit different pictures about unfolding and
refolding pathways. The method is effective in dealing with high degree of
heterogeneity in distribution, capturing structural features in multiple
pathways, and being less sensitive to the distance metric than nonlinear
dimensionality reduction or geometric embedding methods. It provides us a
systematic tool to explore the low density intermediate states in complex
biomolecular folding systems.


  Epidemic models are always simplifications of real world epidemics. Which
real world features to include, and which simplifications to make, depend both
on the disease of interest and on the purpose of the modelling. In the present
paper we discuss some such purposes for which a \emph{stochastic} model is
preferable to a \emph{deterministic} counterpart. The two main examples
illustrate the importance of allowing the infectious and latent periods to be
random when focus lies on the \emph{probability} of a large epidemic outbreak
and/or on the initial \emph{speed}, or growth rate, of the epidemic. A
consequence of the latter is that estimation of the basic reproduction number
$R_0$ is sensitive to assumptions about the distributions of the infectious and
latent periods when using the data from the early stages of an outbreak, which
we illustrate with data from the SARS outbreak. Some further examples are also
discussed as are some practical consequences related to these stochastic
aspects.


  Quantitative scaling relationships among body mass, temperature and metabolic
rate of organisms are still controversial, while resolution may be further
complicated through the use of different and possibly inappropriate approaches
to statistical analysis. We propose the application of a modelling strategy
based on Akaike's information criteria and non-linear model fitting (nlm).
Accordingly, we collated and modelled available data at intraspecific level on
the individual standard metabolic rate of Antarctic microarthropods as a
function of body mass (M), temperature (T), species identity (S) and high rank
taxa to which species belong (G) and tested predictions from Metabolic Scaling
Theory. We also performed allometric analysis based on logarithmic
transformations (lm). Conclusions from lm and nlm approaches were different.
Best-supported models from lm incorporated T, M and S. The estimates of the
allometric scaling exponent b linking body mass and metabolic rate indicated no
interspecific difference and resulted in a value of 0.696 +/- 0.105 (mean +/-
95% CI). In contrast, the four best-supported nlm models suggested that both
the scaling exponent and activation energy significantly vary across the high
rank taxa to which species belong, with mean values of b ranging from about 0.6
to 0.8. We therefore reached two conclusions: 1) published analyses of
arthropod metabolism based on logarithmic data may be biased by data
transformation; 2) non-linear models applied to Antarctic microarthropod
metabolic rate suggest that intraspecific scaling of standard metabolic rate in
Antarctic microarthropods is highly variable and can be characterised by
scaling exponents that greatly vary within taxa, which may have biased previous
interspecific comparisons that neglected intraspecific variability.


  This paper examines the interplay of the effect of cross immunity and
antibody-dependent enhancement (ADE) in mutistrain diseases. Motivated by
dengue fever, we study a model for the spreading of epidemics in a population
with multistrain interactions mediated by both partial temporary cross immunity
and ADE. Although ADE models have previously been observed to cause chaotic
outbreaks, we show analytically that weak cross immunity has a stabilizing
effect on the system. That is, the onset of disease fluctuations requires a
larger value of ADE with small cross immunity than without. However, strong
cross immunity is shown numerically to cause oscillations and chaotic outbreaks
even for low values of ADE.


  We study the statistical properties of a simple genetic regulatory network
that provides heterogeneity within a population of cells. This network consists
of a binary genetic switch in which stochastic flipping between the two switch
states is mediated by a "flipping" enzyme. Feedback between the switch state
and the flipping rate is provided by a linear feedback mechanism: the flipping
enzyme is only produced in the on switch state and the switching rate depends
linearly on the copy number of the enzyme. This work generalises the model of
[Phys. Rev. Lett., 101, 118104] to a broader class of linear feedback systems.
We present a complete analytical solution for the steady-state statistics of
the number of enzyme molecules in the on and off states, for the general case
where the enzyme can mediate flipping in either direction. For this general
case we also solve for the flip time distribution, making a connection to first
passage and persistence problems in statistical physics. We show that the
statistics of the model are non-Poissonian, leading to a peak in the flip time
distribution. The occurrence of such a peak is analysed as a function of the
parameter space. We present a new relation between the flip time distributions
measured for two relevant choices of initial condition. We also introduce a new
correlation measure to show that this model can exhibit long-lived temporal
correlations, thus providing a primitive form of cellular memory. Motivated by
DNA replication as well as by evolutionary mechanisms involving gene
duplication, we study the case of two switches in the same cell. This results
in correlations between the two switches; these can either positive or negative
depending on the parameter regime.


  Proteins from the kinesin-8 family promote microtubule (MT) depolymerization,
a process thought to be important for the control of microtubule length in
living cells. In addition to this MT shortening activity, kinesin 8s are motors
that show plus-end directed motility on MTs. Here we describe a simple model
that incorporates directional motion and destabilization of the MT plus end by
kinesin 8. Our model quantitatively reproduces the key features of
length-vs-time traces for stabilized MTs in the presence of purified kinesin 8,
including length-dependent depolymerization. Comparison of model predictions
with experiments suggests that kinesin 8 depolymerizes processively, i.e., one
motor can remove multiple tubulin dimers from a stabilized MT. Fluctuations in
MT length as a function of time are related to depolymerization processivity.
We have also determined the parameter regime in which the rate of MT
depolymerization is length dependent: length-dependent depolymerization occurs
only when MTs are sufficiently short; this crossover is sensitive to the bulk
motor concentration.


  To regulate a particular gene, a transcription factor (TF) needs to bind a
specific genome location. How is this genome address specified amid the
presence of ~10^6-10^9 decoy sites? Our analysis of 319 known TF binding motifs
clearly demonstrates that prokaryotes and eukaryotes use strikingly different
strategies to target TFs to specific genome locations; eukaryotic TFs exhibit
widespread nonfunctional binding and require clustering of sites in regulatory
regions for specificity.


  Protein detection methods after electrophoresis have to be sensitive,
homogeneous, and not to impair downstream analysis of proteins by MS. Speed,
low cost, and user friendliness are also favored features. Silver staining
combines many of these features, but its compatibility with MS is limited. We
describe here, a new variant of silver staining that is completely
formaldehyde-free. Reducing sugars in alkaline borate buffer are used as
developers. While keeping the benefits of silver staining, this method is shown
to afford a much better performance in terms of compatibility with MS, both in
PMF by MALDI and in LC/ESI/MS/MS.


  This paper considers a stochastic SIR (susceptible$\to$infective$\to$removed)
epidemic model in which individuals may make infectious contacts in two ways,
both within `households' (which for ease of exposition are assumed to have
equal size) and along the edges of a random graph describing additional social
contacts. Heuristically-motivated branching process approximations are
described, which lead to a threshold parameter for the model and methods for
calculating the probability of a major outbreak, given few initial infectives,
and the expected proportion of the population who are ultimately infected by
such a major outbreak. These approximate results are shown to be exact as the
number of households tends to infinity by proving associated limit theorems.
Moreover, simulation studies indicate that these asymptotic results provide
good approximations for modestly sized finite populations. The extension to
unequal sized households is discussed briefly.


  In this paper we establish a relation between the spread of infectious
diseases and the dynamics of so called M/G/1 queues with processor sharing. The
in epidemiology well known relation between the spread of epidemics and
branching processes and the in queueing theory well known relation between
M/G/1 queues and birth death processes will be combined to provide a framework
in which results from queueing theory can be used in epidemiology and vice
versa.
  In particular, we consider the number of infectious individuals in a standard
SIR epidemic model at the moment of the first detection of the epidemic, where
infectious individuals are detected at a constant per capita rate. We use a
result from the literature on queueing processes to show that this number of
infectious individuals is geometrically distributed.


  Coexistence of competing species is, due to unavoidable fluctuations, always
transient. In this Letter, we investigate the ultimate survival probabilities
characterizing different species in cyclic competition. We show that they often
obey a surprisingly simple, though non-trivial behavior. Within a model where
coexistence is neutrally stable, we demonstrate a robust zero-one law: When the
interactions between the three species are (generically) asymmetric, the
`weakest' species survives at a probability that tends to one for large
population sizes, while the other two are guaranteed to extinct. We rationalize
our findings from stochastic simulations by an analytic approach.


  Eukaryotic cells are often exposed to fluctuations in growth conditions as
well as endogenous and exogenous stress-related agents. In addition, during
development global patterns of gene transcription change dramatically, and
these changes are associated with altered patterns of DNA replication. In
metazoan embryos, for example, transcription is repressed globally and any
sequence in the genome can serve as a site for the start of DNA synthesis. As
transcription is activated and a G1 phase imposed, the pattern of replication
adapts to these changes by restricting the sites where DNA synthesis begins.
Recent evidence indicates that each unit of replication, or replicon, is
specified by two or more potential replication origins, but only one is
selected to initiate replication of the replicon. How the cell distinguishes
between potential origins, and how it selects a given origin of replication
remain unclear. This raises important questions concerning the nature and
definition of the eukaryotic replicon. In the following we will review emerging
evidence concerning the mechanisms involved in regulating replication origins
during both the normal and perturbed eukaryotic cell cycle.


  Mathematical theory of selection is developed within the frameworks of
general models of inhomogeneous populations with continuous time. Methods that
allow us to study the distribution dynamics under natural selection and to
construct explicit solutions of the models are developed. All statistical
characteristics of interest, such as the mean values of the fitness or any
trait can be computed effectively, and the results depend in a crucial way on
the initial distribution. The developed theory provides an effective method for
solving selection systems; it reduces the initial complex model to a special
system of ordinary differential equations (the escort system). Applications of
the method to the Price equations are given; the solutions of some particular
inhomogeneous Malthusian, Ricker and logistic-like models used but not solved
in the literature are derived in explicit form.


  Replicator equations (RE) are among the basic tools in mathematical theory of
selection and evolution. We develop a method for reducing a wide class of the
RE, which in general are systems of differential equations in Banach space to
escort systems of ODEs that in many cases can be explored analytically. The
method has potential for different applications; some examples are given.


  The recent years have seen the emergence of diseases which have spread very
quickly all around the world either through human travels like SARS or animal
migration like avian flu. Among the biggest challenges raised by infectious
emerging diseases, one is related to the constant mutation of the viruses which
turns them into continuously moving targets for drug and vaccine discovery.
Another challenge is related to the early detection and surveillance of the
diseases as new cases can appear just anywhere due to the globalization of
exchanges and the circulation of people and animals around the earth, as
recently demonstrated by the avian flu epidemics. For 3 years now, a
collaboration of teams in Europe and Asia has been exploring some innovative in
silico approaches to better tackle avian flu taking advantage of the very large
computing resources available on international grid infrastructures. Grids were
used to study the impact of mutations on the effectiveness of existing drugs
against H5N1 and to find potentially new leads active on mutated strains. Grids
allow also the integration of distributed data in a completely secured way. The
paper presents how we are currently exploring how to integrate the existing
data sources towards a global surveillance network for molecular epidemiology.


  We introduce a new percolation model to describe and analyze the spread of an
epidemic on a general directed and locally finite graph. We assign a
two-dimensional random weight vector to each vertex of the graph in such a way
that the weights of different vertices are i.i.d., but the two entries of the
vector assigned to a vertex need not be independent. The probability for an
edge to be open depends on the weights of its end vertices, but conditionally
on the weights, the states of the edges are independent of each other. In an
epidemiological setting, the vertices of a graph represent the individuals in a
(social) network and the edges represent the connections in the network. The
weights assigned to an individual denote its (random) infectivity and
susceptibility, respectively. We show that one can bound the percolation
probability and the expected size of the cluster of vertices that can be
reached by an open path starting at a given vertex from above and below by the
corresponding quantities for respectively independent bond and site percolation
with certain densities; this generalizes a result of Kuulasmaa. Many models in
the literature are special cases of our general model.


  The symmetrical network theory is a framework for understanding the immune
system, that dates back to the mid 1970s. The symmetrical network theory is
based on symmetrical stimulatory, inhibitory and killing interactions between
clones that are specific for each other. Previous papers described roles for
helper and suppressor T cells in regulating immune responses and a model for
HIV pathogenesis. This paper extends the theory to account for regulatory T
cells that include three types of suppressor cells called Ts1, Ts2 and Ts3, and
two types of helper cells called Th1 and Th2. The theory leads to a concept for
an HIV vaccine, namely a reagent commonly known as IVIG, to be administered in
small amounts in an immunogenic form via an immunogenic route. Predictions are
made for experiments in mice and macaque monkeys.


  The bacterium E.Coli swims in a zig-zag manner, in a series of straight runs
and tumbles occurring alternately, with the run-durations dependent on the
local spatial gradient of chemo-attractants/repellants. This enables the
organism to move towards nutrient sources and move away from toxins. The signal
transduction network of E.Coli has been well-characterized, and theoretical
modeling has been used, with some success, in understanding its many remarkable
features, including the near-perfect adaptation to spatially uniform stimulus.
We study a reduced form of this network, with 3 methylation states for the
receptor instead of 5. We derive an analytical form of the response function of
the tumbling rate and use it to compute the drift velocity of the bacterium in
the presence of a weak spatial attractant gradient.


  The system-level dynamics of multivalent biomolecular interactions can be
simulated using a rule-based kinetic Monte Carlo method in which a rejection
sampling strategy is used to generate reaction events. This method becomes
inefficient when simulating aggregation processes with large biomolecular
complexes. Here, we present a rejection-free method for determining the
kinetics of multivalent biomolecular interactions, and we apply the method to
simulate simple models for ligand-receptor interactions. Simulation results
show that performance of the rejection-free method is equal to or better than
that of the rejection method over wide parameter ranges, and the rejection-free
method is more efficient for simulating systems in which aggregation is
extensive. The rejection-free method reported here should be useful for
simulating a variety of systems in which multisite molecular interactions yield
large molecular aggregates.


  The perennial problem of "how many clusters?" remains an issue of substantial
interest in data mining and machine learning communities, and becomes
particularly salient in large data sets such as populational genomic data where
the number of clusters needs to be relatively large and open-ended. This
problem gets further complicated in a co-clustering scenario in which one needs
to solve multiple clustering problems simultaneously because of the presence of
common centroids (e.g., ancestors) shared by clusters (e.g., possible descents
from a certain ancestor) from different multiple-cluster samples (e.g.,
different human subpopulations). In this paper we present a hierarchical
nonparametric Bayesian model to address this problem in the context of
multi-population haplotype inference. Uncovering the haplotypes of single
nucleotide polymorphisms is essential for many biological and medical
applications. While it is uncommon for the genotype data to be pooled from
multiple ethnically distinct populations, few existing programs have explicitly
leveraged the individual ethnic information for haplotype inference. In this
paper we present a new haplotype inference program, Haploi, which makes use of
such information and is readily applicable to genotype sequences with thousands
of SNPs from heterogeneous populations, with competent and sometimes superior
speed and accuracy comparing to the state-of-the-art programs. Underlying
Haploi is a new haplotype distribution model based on a nonparametric Bayesian
formalism known as the hierarchical Dirichlet process, which represents a
tractable surrogate to the coalescent process. The proposed model is
exchangeable, unbounded, and capable of coupling demographic information of
different populations.


  RNA polymerase (RNAP) is a mobile molecular workshop that polymerizes a RNA
molecule by adding monomeric subunits one by one, while moving step by step on
the DNA template itself. Here we develop a theoretical model by incorporating
the steric interactions of the RNAPs and their mechanochemical cycles which
explicitly captures the cyclical shape changes of each motor. Using this model,
we explain not only the dependence of the average velocity of a RNAP on the
externally applied load force, but also predict a {\it nonmotonic} variation of
the average velocity on external torque. We also show the effect of steric
interactions of the motors on the total rate of RNA synthesis. In principle,
our predictions can be tested by carrying out {\it in-vitro} experiments which
we suggest here.


  The quality and ease of proteomics analysis depends on the performance of the
analytical tools used, and thus of the performances of the protein separation
tools used to deconvolute complex protein samples. Among protein samples,
membrane proteins are one of the most difficult sample classes, because of
their hydrophobicity and embedment in the lipid bilayers. This review deals
with the recent progresses and advances made in the separation of membrane
proteins by 2-DE separating only denatured proteins. Traditional 2-D methods,
i.e., methods using IEF in the first dimension are compared to methods using
only zone electrophoresis in both dimensions, i.e., electrophoresis in the
presence of cationic or anionic detergents. The overall performances and fields
of application of both types of method is critically examined, as are future
prospects for this field.


  This paper investigates the dynamics of biomass in a marine ecosystem. A
stochastic process is defined in which organisms undergo jumps in body size as
they catch and eat smaller organisms. Using a systematic expansion of the
master equation, we derive a deterministic equation for the macroscopic
dynamics, which we call the deterministic jump-growth equation, and a linear
Fokker-Planck equation for the stochastic fluctuations. The McKendrick--von
Foerster equation, used in previous studies, is shown to be a first-order
approximation, appropriate in equilibrium systems where predators are much
larger than their prey. The model has a power-law steady state consistent with
the approximate constancy of mass density in logarithmic intervals of body mass
often observed in marine ecosystems. The behaviours of the stochastic process,
the deterministic jump-growth equation and the McKendrick--von Foerster
equation are compared using numerical methods. The numerical analysis shows two
classes of attractors: steady states and travelling waves.


  Testing fit of data to model is fundamentally important to any science, but
publications in the field of phylogenetics rarely do this. Such analyses
discard fundamental aspects of science as prescribed by Karl Popper. Indeed,
not without cause, Popper (1978) once argued that evolutionary biology was
unscientific as its hypotheses were untestable. Here we trace developments in
assessing fit from Penny et al. (1982) to the present. We compare the general
log-likelihood ratio (the G or G2 statistic) statistic between the evolutionary
tree model and the multinomial model with that of marginalized tests applied to
an alignment (using placental mammal coding sequence data). It is seen that the
most general test does not reject the fit of data to model (p~0.5), but the
marginalized tests do. Tests on pair-wise frequency (F) matrices, strongly (p <
0.001) reject the most general phylogenetic (GTR) models commonly in use. It is
also clear (p < 0.01) that the sequences are not stationary in their nucleotide
composition. Deviations from stationarity and homogeneity seem to be unevenly
distributed amongst taxa; not necessarily those expected from examining other
regions of the genome. By marginalizing the 4t patterns of the i.i.d. model to
observed and expected parsimony counts, that is, from constant sites, to
singletons, to parsimony informative characters of a minimum possible length,
then the likelihood ratio test regains power, and it too rejects the
evolutionary model with p << 0.001. Given such behavior over relatively recent
evolutionary time, readers in general should maintain a healthy skepticism of
results, as the scale of the systematic errors in published analyses may really
be far larger than the analytical methods (e.g., bootstrap) report.


  Stochastic networks are a plausible representation of the relational
information among entities in dynamic systems such as living cells or social
communities. While there is a rich literature in estimating a static or
temporally invariant network from observation data, little has been done toward
estimating time-varying networks from time series of entity attributes. In this
paper we present two new machine learning methods for estimating time-varying
networks, which both build on a temporally smoothed $l_1$-regularized logistic
regression formalism that can be cast as a standard convex-optimization problem
and solved efficiently using generic solvers scalable to large networks. We
report promising results on recovering simulated time-varying networks. For
real data sets, we reverse engineer the latent sequence of temporally rewiring
political networks between Senators from the US Senate voting records and the
latent evolving regulatory networks underlying 588 genes across the life cycle
of Drosophila melanogaster from the microarray time course.


  In this thesis, we have studied the large scale structure and system level
dynamics of certain biological networks using tools from graph theory,
computational biology and dynamical systems. We study the structure and
dynamics of large scale metabolic networks inside three organisms, Escherichia
coli, Saccharomyces cerevisiae and Staphylococcus aureus. We also study the
dynamics of the large scale genetic network controlling E. coli metabolism. We
have tried to explain the observed system level dynamical properties of these
networks in terms of their underlying structure. Our studies of the system
level dynamics of these large scale biological networks provide a different
perspective on their functioning compared to that obtained from purely
structural studies. Our study also leads to some new insights on features such
as robustness, fragility and modularity of these large scale biological
networks. We also shed light on how different networks inside the cell such as
metabolic networks and genetic networks are interrelated to each other.


  Divergence time estimation requires the reconciliation of two major sources
of data. These are fossil and/or biogeographic evidence that give estimates of
the absolute age of nodes (ancestors) and molecular estimates that give us
estimates of the relative ages of nodes in a molecular evolutionary tree. Both
forms of data are often best characterized as yielding continuous probability
distributions on nodes. Here, the distributions modeling older fossil
calibrations within the tree of placental (eutherian) mammals are reconsidered.
In particular the Horse/Rhino, Human/Tarsier, Whale/ Hippo, Rabbit/Pika and
Rodentia calibrations are reexamined and adjusted. Inferring the relative ages
of nodes in a phylogeny also requires the assumption of a model of evolutionary
rate change across the tree. Here nine models of evolutionary rate change, are
combined with various continuous distributions modeling fossil calibrations.
Fit of model is measured both relative to a normalized fit, which assumes that
all models fit well in the absence of multiple fossil calibrations, and also by
the linearity of their residuals. The normalized fit used attempts to track
twice the log likelihood difference from the best expected model. The results
suggest there is a very large difference in the age of the root proposed by
calibrations in Supraprimates (informally Euarchontoglires) versus
Laurasiatheria. Combining both sets of calibrations results in the penalty
function vastly increasing in all cases. These issues remain irrespective of
the model used or whether the newer calibrations are used.


  This supplementary information contains detailed derivations, comparison to
experiment, and discussion of other miscellaneous issues omitted from the main
text.


  Three coarse-grained models of the double-stranded DNA are proposed and
compared in the context of mechanical manipulation such as twisting and various
schemes of stretching. The models differ in the number of effective beads
(between two and five) representing each nucleotide. They all show similar
behavior and, in particular, lead to a torque-force phase diagrams
qualitatively consistent with experiments and all-atom simulations.


  We incorporate hydrodynamic interactions (HI) in a coarse-grained and
structure-based model of proteins by employing the Rotne-Prager hydrodynamic
tensor. We study several small proteins and demonstrate that HI facilitate
folding. We also study HIV-1 protease and show that HI make the flap closing
dynamics faster. The HI are found to affect time correlation functions in the
vicinity of the native state even though they have no impact on same time
characteristics of the structure fluctuations around the native state.


  Synthesis of tri functional electrically conducting, optical and magnetic
nano-chain of Nicore-Aushell has been discussed here. Our Investigation
indicates that such material attached with biomolecule DNA in chain form will
have great potentiality in medical instrument and bio computer device.


  We present a new approach to nicotinic receptor kinetics and a new model
explaining random variabilities in the duration of open events. The model gives
new interpretation on brief and long receptor openings and predicts (for two
identical binding sites) the presence of three components in the open time
distribution: two brief and a long. We also present the physical model of the
receptor block. This picture naturally and universally explains receptor
desensitization, the phenomenon of central importance in cellular signaling.
The model is based on single-channel experiments concerning the effects of
hydrocortisone (HC) on the kinetics of control wild-type (WT) and mutated
alphaD200Q mouse nicotinic acetylcholine receptors (nAChRs), expressed in HEK
293 cells. The appendix contains an original result from probability renewal
theory: a derivation of the probability distribution function for the duration
of a process performed by two independent servers.


  In a dynamic social or biological environment, the interactions between the
actors can undergo large and systematic changes. In this paper we propose a
model-based approach to analyze what we will refer to as the dynamic tomography
of such time-evolving networks. Our approach offers an intuitive but powerful
tool to infer the semantic underpinnings of each actor, such as its social
roles or biological functions, underlying the observed network topologies. Our
model builds on earlier work on a mixed membership stochastic blockmodel for
static networks, and the state-space model for tracking object trajectory. It
overcomes a major limitation of many current network inference techniques,
which assume that each actor plays a unique and invariant role that accounts
for all its interactions with other actors; instead, our method models the role
of each actor as a time-evolving mixed membership vector that allows actors to
behave differently over time and carry out different roles/functions when
interacting with different peers, which is closer to reality. We present an
efficient algorithm for approximate inference and learning using our model; and
we applied our model to analyze a social network between monks (i.e., the
Sampson's network), a dynamic email communication network between the Enron
employees, and a rewiring gene interaction network of fruit fly collected
during its full life cycle. In all cases, our model reveals interesting
patterns of the dynamic roles of the actors.


  Due to the dynamic nature of biological systems, biological networks
underlying temporal process such as the development of {\it Drosophila
melanogaster} can exhibit significant topological changes to facilitate dynamic
regulatory functions. Thus it is essential to develop methodologies that
capture the temporal evolution of networks, which make it possible to study the
driving forces underlying dynamic rewiring of gene regulation circuity, and to
predict future network structures. Using a new machine learning method called
Tesla, which builds on a novel temporal logistic regression technique, we
report the first successful genome-wide reverse-engineering of the latent
sequence of temporally rewiring gene networks over more than 4000 genes during
the life cycle of \textit{Drosophila melanogaster}, given longitudinal gene
expression measurements and even when a single snapshot of such measurement
resulted from each (time-specific) network is available. Our methods offer the
first glimpse of time-specific snapshots and temporal evolution patterns of
gene networks in a living organism during its full developmental course. The
recovered networks with this unprecedented resolution chart the onset and
duration of many gene interactions which are missed by typical static network
analysis, and are suggestive of a wide array of other temporal behaviors of the
gene network over time not noticed before.


  We propose a physical model to describe the mechanisms of two major scenarios
of the genetic code evolution, the codon capture and ambiguous intermediate
scenarios, in a consistent manner. We sketch the lowest dimensional version of
our model, a minimal model, by introducing a physical quantity, codon level. On
the basis of the hierarchical structure of the codon levels two scenarios are
classified into two different routes of the evolutional process. In the case of
the ambiguous intermediate scenario we perform a simulation implemented cost
selection of amino acids and confirm a rapid transition of the code change.
Such rapidness reduces uncomfortableness of the non-unique translation of the
code at intermediate state that is the weakness of the scenario.


  We introduce a stochastic lattice gas model including two particle species
and two parallel lanes. One lane with exclusion interaction and directed motion
and the other lane without exclusion and unbiased diffusion, mimicking a
micotubule filament and the surrounding solution. For a high binding affinity
to the filament, jam-like situations dominate the system's behaviour. The
fundamental process of position exchange of two particles is approximated. In
the case of a many-particle system, we were able to identify a regime in which
the system is rather homogenous presenting only small accumulations of
particles and a regime in which an important fraction of all particles
accumulates in the same cluster. Numerical data proposes that this cluster
formation will occur at all densities for large system sizes. Coupling of
several filaments leads to an enhanced cluster formation compared to the
uncoupled system, suggesting that efficient bidirectional transport on
one-dimensional filaments relies on long-ranged interactions and track
formation.


  The distribution of species body size within taxonomic groups exhibits a
heavy right-tail extending over many orders of magnitude, where most species
are significantly larger than the smallest species. We provide a simple model
of cladogenetic diffusion over evolutionary time that omits explicit mechanisms
for inter-specific competition and other microevolutionary processes yet fully
explains the shape of this distribution. We estimate the model's parameters
from fossil data and find that it robustly reproduces the distribution of 4002
mammal species from the late Quaternary. The observed fit suggests that the
asymmetric distribution arises from a fundamental tradeoff between the
short-term selective advantages (Cope's rule) and long-term selective risks of
increased species body size, in the presence of a taxon-specific lower limit on
body size.


  Biochemical networks can respond to temporal characteristics of time-varying
signals. To understand how reliably biochemical networks can transmit
information we must consider how an input signal as a function of time--the
input trajectory--can be mapped onto an output trajectory. Here we estimate the
mutual information between in- and output trajectories using a Gaussian model.
We study how reliably the chemotaxis network of E. coli can transmit
information on the ligand concentration to the flagellar motor, and find the
input power spectrum that maximizes the information transmission rate.


  In our previous publication, a framework for information flow in interaction
networks based on random walks with damping was formulated with two fundamental
modes: emitting and absorbing. While many other network analysis methods based
on random walks or equivalent notions have been developed before and after our
earlier work, one can show that they can all be mapped to one of the two modes.
In addition to these two fundamental modes, a major strength of our earlier
formalism was its accommodation of context-specific directed information flow
that yielded plausible and meaningful biological interpretation of protein
functions and pathways. However, the directed flow from origins to destinations
was induced via a potential function that was heuristic. Here, with a
theoretically sound approach called the channel mode, we extend our earlier
work for directed information flow. This is achieved by constructing a
potential function facilitating a purely probabilistic interpretation of the
channel mode. For each network node, the channel mode combines the solutions of
emitting and absorbing modes in the same context, producing what we call a
channel tensor. The entries of the channel tensor at each node can be
interpreted as the amount of flow passing through that node from an origin to a
destination. Similarly to our earlier model, the channel mode encompasses
damping as a free parameter that controls the locality of information flow.
Through examples involving the yeast pheromone response pathway, we illustrate
the versatility and stability of our new framework.


  Molecular motors are essential components for the biophysical functions of
the cell. Our current quantitative understanding of how multiple motors move
along a single track is not complete; even though models and theories for
single motor chemomechanics abound. Recently, M.J.I. M$\ddot{\rm u}$ller {\em
et al.} have developed a tug-of-war model to describe the bidirectional
movement of the cargo (PNAS(2008) 105(12) P4609-4614). Through Monte Carlo
simulations, they discovered that the tug-of-war model exhibits several
qualitative different motility regimes, which depend on the precise value of
single motor parameters, and they suggested the sensitivity can be used by a
cell to regulate its cargo traffic. In the present paper, we carry out a
thorough analysis of the tug-of-war model. All the stable, i.e., biophysically
observable, steady states are obtained. Depending on several parameters, the
system exhibits either uni-, bi- or tristability. Based on the separating
boundary of the different stable states and the initial numbers of the
different motor species that are bound to the track, the steady state of the
cargo movement can be predicted, and consequently the steady state velocity can
be obtained. It is found that, the velocity, even the direction, of the cargo
movement change with the initial numbers of the motors which are bound to the
track and several other parameters.


  We address the issue of the effects of considering a network of contacts on
the emergence of cooperation on social dilemmas under myopic best response
dynamics. We begin by summarizing the main features observed under less
intellectually demanding dynamics, pointing out their most relevant general
characteristics. Subsequently we focus on the new framework of best response.
By means of an extensive numerical simulation program we show that, contrary to
the rest of dynamics considered so far, best response is largely unaffected by
the underlying network, which implies that, in most cases, no promotion of
cooperation is found with this dynamics. We do find, however, nontrivial
results differing from the well-mixed population in the case of coordination
games on lattices, which we explain in terms of the formation of spatial
clusters and the conditions for their advancement, subsequently discussing
their relevance to other networks.


  The interaction between natural selection and random mutation is frequently
debated in recent years. Does similar dilemma also exist in the evolution of
real networks such as biological networks? In this paper, we try to discuss
this issue by a simple model system, in which the topological structure of
networks is repeatedly modified and selected in order to make them have better
performance in dynamical processes. Interestingly, when the networks with
optimal performance deviate from the steady state networks under pure
mutations, we find the evolution behaves as a balance between mutation and
selection. Furthermore, when the timescales of mutations and dynamical
processes are comparable with each other, the steady state of evolution is
mainly determined by mutation. On the opposite side, when the timescale of
mutations is much longer than that of dynamical processes, selection dominates
the evolution and the steady-state networks turn to have much improved
performance and highly heterogeneous structures. Despite the simplicity of our
model system, this finding could give useful indication to detect the
underlying mechanisms that rein the evolution of real systems.


  We test the claim by Yu et al. -- presented in Science 322, 104 (2008) --
that the degree distribution of the yeast (Saccharomyces cerevisiae)
protein-interaction network is best approximated by a power law. Yu et al.
consider three versions of this network. In all three cases, however, we find
the most likely power-law model of the data is distinct from and incompatible
with the one given by Yu et al. Only one network admits good statistical
support for any power law, and in that case, the power law explains only the
distribution of the upper 10% of node degrees. These results imply that there
is considerably more structure present in the yeast interactome than suggested
by Yu et al., and that these networks should probably not be called "scale
free."


  Enzymes are on the front lines of evolution. All living organisms rely on
highly efficient, specific enzymes for growth, sustenance, and reproduction;
and many diseases are a consequence of a mutation on an enzyme that affects its
catalytic function. It follows that the function of an enzyme affects the
fitness of an organism, but just as rightfully true, the function of an enzyme
affects the fitness of itself. Understanding how the complexity of enzyme
structure relates to its essential function will unveil the fundamental
mechanisms of evolution, and, perhaps, shed light on strategies used by ancient
replicators. This paper presents evidence that supports the hypothesis that
enzymes, and proteins in general, are the manifestation of the coevolution of
two opposing forces. The synthesis of enzyme architecture, stability, function,
evolutionary relationships, and evolvability shows that the complexity of
macromolecules is a consequence of the function it provides.


  The cooperative action of many molecular motors is essential for dynamic
processes such as cell motility and mitosis. This action can be studied by
using motility assays in which the motion of cytoskeletal filaments over a
surface coated with motor proteins is tracked. In previous studies of
actin-myosin II systems, fast directional motion was observed, reflecting the
tendency of myosin II motors to propagate unidirectionally along actin
filaments. Here, we present a motility assay with actin bundles consisting of
short filamentous segments with randomly alternating polarities. These actin
tracks exhibit bidirectional motion with macroscopically large time intervals
(of the order of several seconds) between direction reversals. Analysis of this
bidirectional motion reveals that the characteristic reversal time,
$\tau_{rev}$, does not depend on the size of the moving bundle or on the number
of motors, $N$. This observation contradicts previous theoretical calculations
based on a two-state ratchet model [Badoual et al., Proc. Natl. Acad. Sci. USA,
vol. 99, p. 6696 (2002)], predicting an exponential increase of $\tau_{rev}$
with $N$. We present a modified version of this model that takes into account
the elastic energy due to the stretching of the actin track by the myosin II
motors. The new model yields a very good quantitative agreement with the
experimental results.


  We investigate the fitness advantage associated with the robustness of a
phenotype against deleterious mutations using deterministic mutation-selection
models of quasispecies type equipped with a mesa shaped fitness landscape. We
obtain analytic results for the robustness effect which become exact in the
limit of infinite sequence length. Thereby, we are able to clarify a seeming
contradiction between recent rigorous work and an earlier heuristic treatment
based on a mapping to a Schr\"odinger equation. We exploit the quantum
mechanical analogy to calculate a correction term for finite sequence lengths
and verify our analytic results by numerical studies. In addition, we
investigate the occurrence of an error threshold for a general class of
epistatic landscape and show that diminishing epistasis is a necessary but not
sufficient condition for error threshold behavior.


  The bacterial flagellar motor drives the rotation of flagellar filaments and
enables many species of bacteria to swim. Torque is generated by interaction of
stator units, anchored to the peptidoglycan cell wall, with the rotor. Recent
experiments [Yuan, J. & Berg, H. C. (2008) PNAS 105, 1182-1185] show that near
zero load the speed of the motor is independent of the number of stators. Here,
we introduce a mathematical model of the motor dynamics that explains this
behavior based on a general assumption that the stepping rate of a stator
depends on the torque exerted by the stator on the rotor. We find that the
motor dynamics can be characterized by two time scales: the moving-time
interval for the mechanical rotation of the rotor and the waiting-time interval
determined by the chemical transitions of the stators. We show that these two
time scales depend differently on the load, and that their crossover provides
the microscopic explanation for the existence of two regimes in the
torque-speed curves observed experimentally. We also analyze the speed
fluctuation for a single motor using our model. We show that the motion is
smoothed by having more stator units. However, the mechanism for such
fluctuation reduction is different depending on the load. We predict that the
speed fluctuation is determined by the number of steps per revolution only at
low load and is controlled by external noise for high load. Our model can be
generalized to study other molecular motor systems with multiple
power-generating units.


  We study the cyclic dominance of three species in two-dimensional constrained
Newman-Watts networks with a four-state variant of the rock-paper-scissors
game. By limiting the maximal connection distance $R_{max}$ in Newman-Watts
networks with the long-rang connection probability $p$, we depict more
realistically the stochastic interactions among species within ecosystems. When
we fix mobility and vary the value of $p$ or $R_{max}$, the Monte Carlo
simulations show that the spiral waves grow in size, and the system becomes
unstable and biodiversity is lost with increasing $p$ or $R_{max}$. These
results are similar to recent results of Reichenbach \textit{et al.} [Nature
(London) \textbf{448}, 1046 (2007)], in which they increase the mobility only
without including long-range interactions. We compared extinctions with or
without long-range connections and computed spatial correlation functions and
correlation length. We conclude that long-range connections could improve the
mobility of species, drastically changing their crossover to extinction and
making the system more unstable.


  Degree distributions have been widely used to characterize biological
networks including food webs, and play a vital role in recent models of food
web structure. While food webs degree distributions have been suggested to
follow various functional forms, to date there has been no mechanistic or
statistical explanation for these forms. Here I introduce models for the degree
distributions of food webs based on the principle of maximum entropy (MaxEnt)
constrained by the number of species, number of links and the number of top or
basal species. The MaxEnt predictions are compared to observed distributions in
51 food webs. The distributions of the number of consumers and resources in 23
(45%) and 35 (69%) of the food webs respectively are not significantly
different at a 95% confidence level from the MaxEnt distribution. While the
resource distributions of niche model webs are well-described by the MaxEnt
model, the consumer distributions are more narrowly distributed than predicted
by the MaxEnt model. These findings offer a new null model for the most
probable degree distributions in food webs. Having an appropriate null
hypothesis in place allows informative study of the deviations from it; for
example these results suggest that there is relatively little pressure favoring
generalist versus specialist consumption strategies but that there is more
pressure driving the consumer distribution away from the MaxEnt form. Given the
methodological idiosyncrasies of current food web data, further study of such
deviations will need to consider both biological drivers and methodological
bias.


  We have used the Monte Carlo based computer models to show that selection
pressure could affect the distribution of recombination hotspots along the
chromosome. Close to critical crossover rate, where genomes may switch between
the Darwinian purifying selection or complementation of haplotypes, the
distribution of recombination events and the force of selection exerted on
genes affect the structure of chromosomes. The order of expression of gene s
and their location on chromosome may decide about the extinction or survival of
competing populations.


  Fights-to-the-death occur in many natural, medical and commercial settings.
Standard mass action theory and conventional wisdom imply that the minority
(i.e. smaller) group's survival time decreases as its relative initial size
decreases, in the absence of replenishment. Here we show that the opposite
actually happens, if the minority group features internal network dynamics. Our
analytic theory provides a unified quantitative explanation for a range of
previously unexplained data, and predicts how losing battles in a medical or
social context might be extended or shortened using third-party intervention.


  Understanding the molecular determinants of specificity in protein-protein
interaction is an outstanding challenge of postgenome biology. The availability
of large protein databases generated from sequences of hundreds of bacterial
genomes enables various statistical approaches to this problem. In this context
covariance-based methods have been used to identify correlation between amino
acid positions in interacting proteins. However, these methods have an
important shortcoming, in that they cannot distinguish between directly and
indirectly correlated residues. We developed a method that combines covariance
analysis with global inference analysis, adopted from use in statistical
physics. Applied to a set of >2,500 representatives of the bacterial
two-component signal transduction system, the combination of covariance with
global inference successfully and robustly identified residue pairs that are
proximal in space without resorting to ad hoc tuning parameters, both for
heterointeractions between sensor kinase (SK) and response regulator (RR)
proteins and for homointeractions between RR proteins. The spectacular success
of this approach illustrates the effectiveness of the global inference approach
in identifying direct interaction based on sequence information alone. We
expect this method to be applicable soon to interaction surfaces between
proteins present in only 1 copy per genome as the number of sequenced genomes
continues to expand. Use of this method could significantly increase the
potential targets for therapeutic intervention, shed light on the mechanism of
protein-protein interaction, and establish the foundation for the accurate
prediction of interacting protein partners.


  This paper develops mathematical models describing the evolutionary dynamics
of both asexually and sexually reproducing populations of diploid unicellular
organisms. We consider two forms of genome organization. In one case, we assume
that the genome consists of two multi-gene chromosomes, while in the second
case we assume that each gene defines a separate chromosome. If the organism
has $ l $ homologous pairs that lack a functional copy of the given gene, then
the fitness of the organism is $ \kappa_l $. The $ \kappa_l $ are assumed to be
monotonically decreasing, so that $ \kappa_0 = 1 > \kappa_1 > \kappa_2 > ... >
\kappa_{\infty} = 0 $. For nearly all of the reproduction strategies we
consider, we find, in the limit of large $ N $, that the mean fitness at
mutation-selection balance is $ \max\{2 e^{-\mu} - 1, 0\} $, where $ N $ is the
number of genes in the haploid set of the genome, $ \epsilon $ is the
probability that a given DNA template strand of a given gene produces a mutated
daughter during replication, and $ \mu = N \epsilon $. The only exception is
the sexual reproduction pathway for the multi-chromosomed genome. Assuming a
multiplicative fitness landscape where $ \kappa_l = \alpha^{l} $ for $ \alpha
\in (0, 1) $, this strategy is found to have a mean fitness that exceeds the
mean fitness of all of the other strategies. Furthermore, while the other
reproduction strategies experience a total loss of viability due to the steady
accumulation of deleterious mutations once $ \mu $ exceeds $ \ln 2 $, no such
transition occurs in the sexual pathway. The results of this paper demonstrate
a selective advantage for sexual reproduction with fewer and much less
restrictive assumptions than previous work.


  Much is now known about the consistency of Bayesian updating on
infinite-dimensional parameter spaces with independent or Markovian data.
Necessary conditions for consistency include the prior putting enough weight on
the correct neighborhoods of the data-generating distribution; various
sufficient conditions further restrict the prior in ways analogous to capacity
control in frequentist nonparametrics. The asymptotics of Bayesian updating
with mis-specified models or priors, or non-Markovian data, are far less well
explored. Here I establish sufficient conditions for posterior convergence when
all hypotheses are wrong, and the data have complex dependencies. The main
dynamical assumption is the asymptotic equipartition (Shannon-McMillan-Breiman)
property of information theory. This, along with Egorov's Theorem on uniform
convergence, lets me build a sieve-like structure for the prior. The main
statistical assumption, also a form of capacity control, concerns the
compatibility of the prior and the data-generating process, controlling the
fluctuations in the log-likelihood when averaged over the sieve-like sets. In
addition to posterior convergence, I derive a kind of large deviations
principle for the posterior measure, extending in some cases to rates of
convergence, and discuss the advantages of predicting using a combination of
models known to be wrong. An appendix sketches connections between these
results and the replicator dynamics of evolutionary theory.


  Intragenomic recombination (crossover) is a very important evolutionary
mechanism. The crossover events are not evenly distributed along the natural
chromosomes. Monte Carlo simulations revealed that frequency of recombinations
decides about the strategy of chromosomes' and genomes' evolution. In large
panmictic populations, under high recombination rate the Darwinian purifying
selection operates keeping the fraction of defective genes at the relatively
low level. In small populations and under low recombination rate the strategy
of complementing haplotypes seems to be more advantageous. Switching between
the two strategies has a character of phase transition - it depends on
inbreeding coefficient and crossover rate. The critical recombination rate
depends also on the size of chromosome. It is also possible, that in one genome
some chromosomes could be under complementing while some other under purifying
selection. Such situation stabilizes genome evolution and reproduction
strategy. It seems that this phenomenon can be responsible for the positive
correlation between kinship and fecundity, recently found in the Islander
population. When large population is forced to enter the complementing
strategy, the phenomenon of sympatric speciation is observed.


  We describe the simulation method of modelling the population evolution using
Monte Carlo based on the Penna model. Individuals in the populations are
represented by their diploid genomes. Genes expressed after the minimum
reproduction age are under a weaker selection pressure and accumulate more
mutations than those expressed before the minimum reproduction age. The
generated gradient of defective genes determines the ageing of individuals and
age-structured populations are very similar to the natural, sexually
reproducing populations. The genetic structure of a population depends on the
way how the random death affects the population. The improvement of the medical
care and healthier life styles are responsible for the increasing of the life
expectancy of humans during the last century. Introducing a noise into the
relations between the genotype, phenotype, and environment, it is possible to
simulate some other effects, like the role of immunological systems and a
mother care. One of the most interesting results was the evolution of sex
chromosomes. Placing the male sex determinants on one chromosome of a pair of
sex chromosomes is enough to condemn it for shrinking if the population is
panmictic (random-mating is assumed). If males are indispensable for taking
care of their offspring and have to be faithful to their females, the male sex
chromosome does not shrink.


  We have presented the basic knowledge on the structure of molecules coding
the genetic information, mechanisms of transfer of this information from DNA to
proteins and phenomena connected with replication of DNA. In particular, we
have described the differences of mutational pressure connected with
replication of the leading and lagging DNA strands. We have shown how the
asymmetric replication of DNA affects the structure of genomes, positions of
genes, their function and amino acid composition. Results of Monte Carlo
simulations of evolution of protein coding sequences have shown a specific role
of genetic code in minimizing the effect of nucleotide substitutions on the
amino acid composition of proteins. The results of simulations were compared
with the results of analyses of genomic and proteomic data bases. This chapter
is considered as an introduction to further chapters where chromosomes with
genes represented by nucleotide sequences were replaced by bitstrings with
single bits representing genes.


  The solution of high-dimensional inference and prediction problems in
computational biology is almost always a compromise between mathematical theory
and practical constraints such as limited computational resources. As time
progresses, computational power increases but well-established inference
methods often remain locked in their initial suboptimal solution. We revisit
the approach of Segal et al. (2003) to infer regulatory modules and their
condition-specific regulators from gene expression data. In contrast to their
direct optimization-based solution we use a more representative centroid-like
solution extracted from an ensemble of possible statistical models to explain
the data. The ensemble method automatically selects a subset of most
informative genes and builds a quantitatively better model for them. Genes
which cluster together in the majority of models produce functionally more
coherent modules. Regulators which are consistently assigned to a module are
more often supported by literature, but a single model always contains many
regulator assignments not supported by the ensemble. Reliably detecting
condition-specific or combinatorial regulation is particularly hard in a single
optimum but can be achieved using ensemble averaging.


  Mixed-species associations are temporary associations between individuals of
different species that are often observed in birds, primates and cetaceans.
They have been interpreted as a strategy to reduce predation risk, enhance
foraging success and/or provide a social advantage. In the archipelago of the
Azores, four species of dolphins are commonly involved in mixed-species
associations: the common dolphin, Delphinus delphis, the bottlenose dolphin,
Tursiops truncatus, the striped dolphin, Stenella coeruleoalba, and the spotted
dolphin, Stenella frontalis. In order to understand the reasons why dolphins
associate, we analysed field data collected since 1999 by research scientists
and trained observers placed onboard fishing vessels. In total, 113
mixed-species groups were observed out of 5720 sightings. The temporal
distribution, habitat (water depth, distance to the coast), behaviour (i.e.
feeding, travelling, socializing), size and composition of mixed-species groups
were compared with those of single-species groups. Results did not support the
predation avoidance hypothesis and gave little support to the social advantage
hypothesis. The foraging advantage hypothesis was the most convincing. However,
the benefits of mixed-species associations appeared to depend on the species.
Associations were likely to be opportunistic in the larger bottlenose dolphin,
while there seemed to be some evolutionary constraints favouring associations
in the rarer striped dolphin. Comparison with previous studies suggests that
the formation of mixed-species groups depends on several environmental factors,
and therefore may constitute an adaptive response.


  This paper introduces constNJ, the first algorithm for phylogenetic
reconstruction of sets of trees with constrained pairwise rooted subtree-prune
regraft (rSPR) distance. We are motivated by the problem of constructing sets
of trees which must fit into a recombination, hybridization, or similar
network. Rather than first finding a set of trees which are optimal according
to a phylogenetic criterion (e.g. likelihood or parsimony) and then attempting
to fit them into a network, constNJ estimates the trees while enforcing
specified rSPR distance constraints. The primary input for constNJ is a
collection of distance matrices derived from sequence blocks which are assumed
to have evolved in a tree-like manner, such as blocks of an alignment which do
not contain any recombination breakpoints. The other input is a set of rSPR
constraints for any set of pairs of trees. ConstNJ is consistent and a strict
generalization of the neighbor-joining algorithm; it uses the new notion of
"maximum agreement partitions" to assure that the resulting trees satisfy the
given rSPR distance constraints.


  Given a dense triplet set $\mathcal{T}$, there arise two interesting
questions: Does there exists any phylogenetic network consistent with
$\mathcal{T}$? And if so, can we find an effective algorithm to construct one?
For cases of networks of levels $k=0$ or 1 or 2, these questions were answered
with effective polynomial algorithms. For higher levels $k$, partial answers
were recently obtained with an $O(|\mathcal{T}|^{k+1})$ time algorithm for
simple networks. In this paper we give a complete answer to the general case.
The main idea is to use a special property of SN-sets in a level-k network. As
a consequence, we can also find the level-k network with the minimum number of
reticulations in polynomial time.


  Peptides and proteins exhibit a common tendency to assemble into highly
ordered fibrillar aggregates, whose formation proceeds in a
nucleation-dependent manner that is often preceded by the formation of
disordered oligomeric assemblies. This process has received much attention
because disordered oligomeric aggregates have been associated with
neurodegenerative disorders such as Alzheimer's and Parkinson's diseases. Here
we describe a self-templated nucleation mechanism that determines the
transition between the initial condensation of polypeptide chains into
disordered assemblies and their reordering into fibrillar structures. The
results that we present show that at the molecular level this transition is due
to the ability of polypeptide chains to reorder within oligomers into fibrillar
assemblies whose surfaces act as templates that stabilise the disordered
assemblies.


  The presence of oligomeric aggregates, which is often observed during the
process of amyloid formation, has recently attracted much attention since it
has been associated with neurodegenerative conditions such as Alzheimer's and
Parkinson's diseases. We provide a description of a sequence-indepedent
mechanism by which polypeptide chains aggregate by forming metastable
oligomeric intermediate states prior to converting into fibrillar structures.
Our results illustrate how the formation of ordered arrays of hydrogen bonds
drives the formation of beta-sheets within the disordered oligomeric aggregates
that form early under the effect of hydrophobic forces. Initially individual
beta-sheets form with random orientations, which subsequently tend to align
into protofilaments as their lengths increases. Our results suggest that
amyloid aggregation represents an example of the Ostwald step rule of first
order phase transitions by showing that ordered cross-beta structures emerge
preferentially from disordered compact dynamical intermediate assemblies.


  In this supporting information we briefly describe the torque-speed
measurement procedure. We show the hook spring compliance used in the
simulations. We analyze the distribution functions of the moving and waiting
time intervals. We study the dependence of the torque plateau region on the
stator jumping rate and the cutoff angle, and the robustness of the results
against different rotor-stator- potential and load-rotor forces.


  The work presented in this paper essentially focuses at providing a
scientific theory to explain the growth of information bearing molecules (size
and information contents) without the need of any enzymatic system. It infers
the footprints of molecular evolution in the cell interior for a property
common to all life forms. It is deducted that molecular crowding is a vital
cellular trait common to the all types of cells (primitive or highly evolved).
It is argued that this trait is pervasive and must have been incorporated at
some stage as a common vital feature of life. If this feature has central
importance it must have been part of the pre-biotic information growth of
information bearing molecules. The thermodynamic consequences of molecular
crowding on the growth of RNA (50-100bp long) in the absence of enzyme system
were calculated.


  The problem of \emph{statistical recognition} is considered, as it arises in
immunobiology, namely, the discrimination of foreign antigens against a
background of the body's own molecules. The precise mechanism of this
foreign-self-distinction, though one of the major tasks of the immune system,
continues to be a fundamental puzzle. Recent progress has been made by van den
Berg, Rand, and Burroughs (2001), who modelled the \emph{probabilistic} nature
of the interaction between the relevant cell types, namely, T-cells and
antigen-presenting cells (APCs). Here, the stochasticity is due to the random
sample of antigens present on the surface of every APC, and to the random
receptor type that characterises individual T-cells. It has been shown
previously that this model, though highly idealised, is capable of reproducing
important aspects of the recognition phenomenon, and of explaining them on the
basis of stochastic rare events. These results were obtained with the help of a
refined large deviation theorem and were thus asymptotic in nature. Simulations
have, so far, been restricted to the straightforward simple sampling approach,
which does not allow for sample sizes large enough to address more detailed
questions. Building on the available large deviation results, we develop an
importance sampling technique that allows for a convenient exploration of the
relevant tail events by means of simulation. With its help, we investigate the
mechanism of statistical recognition in some depth. In particular, we
illustrate how a foreign antigen can stand out against the self background if
it is present in sufficiently many copies, although no \emph{a priori}
difference between self and nonself is built into the model.


  Until recently, the use of Bayesian inference in population genetics was
limited to a few cases because for many realistic population genetic models the
likelihood function cannot be calculated analytically . The situation changed
with the advent of likelihood-free inference algorithms, often subsumed under
the term Approximate Bayesian Computation (ABC). A key innovation was the use
of a post-sampling regression adjustment, allowing larger tolerance values and
as such shifting computation time to realistic orders of magnitude (see
Beaumont et al., 2002). Here we propose a reformulation of the regression
adjustment in terms of a General Linear Model (GLM). This allows the
integration into the framework of Bayesian statistics and the use of its
methods, including model selection via Bayes factors. We then apply the
proposed methodology to the question of population subdivision among western
chimpanzees Pan troglodytes verus.


  We present a numerical study of the effect of knotting on the ejection of
flexible and semiflexible polymers from a spherical, virus-like capsid. The
polymer ejection rate is primarily controlled by the knot, which moves to the
hole in the capsid and then acts as a ratchet. Polymers with more complex knots
eject more slowly and, for large knots, the knot type, and not the flexibility
of the polymer, determines the rate of ejection. We discuss the relation of our
results to the ejection of DNA from viral capsids and conjecture that this
process has the biological advantage of unknotting the DNA before it enters a
cell.


  Many complex systems in mathematical biology and other areas can be described
by the replicator equation. We show that solutions of a wide class of
replicator equations minimize the production of information under
time-dependent constraints, which, in their turn, can be computed explicitly at
every instant due to the system dynamics. Therefore, the Kullback principle of
minimum discrimination information, as well as the maximum entropy principle,
for systems governed by the replicator equations can be derived from the system
dynamics rather than postulated. Applications to the Malthusian inhomogeneous
models, global demography, and the Eigen quasispecies equation are given.


  The results of a long-established investigation into pupal transpiration are
used as a rudimentary data set. These data are then generalised to all
temperatures and humidities by invoking the property of multiplicative
separability, as well as by converting established relationships in terms of
constant humidity at fixed temperature, to alternatives in terms of a
calculated water loss. In this way a formulation which is a series of very
simple, first order, ordinary differential equations is devised. The model is
extended to include a variety of Glossina species using their relative surface
areas, their relative pupal and puparial loss rates and their different 4th
instar excretions. The resulting computational model calculates total, pupal
water loss, consequent mortality and emergence. Remaining fat reserves are a
more tenuous result.
  The model suggests that, while conventional wisdom is often correct in
dismissing variability in transpiration-related pupal mortality as
insignificant, the effects of transpiration can be profound under adverse
conditions and for some species, in general. The model demonstrates how two
gender effects, the more significant one at the drier extremes of tsetse fly
habitat, might arise. The agreement between calculated and measured critical
water losses suggests very little difference in the behaviour of the different
species.


  A two-species spatially extended system of hosts and parasitoids is studied.
There are two distinct kinds of coexistence; one with populations distributed
homogeneously in space and another one with spatiotemporal patterns. In the
latter case, there are noise-sustained oscillations in the population
densities, whereas in the former one the densities are essentially constants in
time with small fluctuations. We introduce several metrics to characterize the
patterns and onset thereof. We also build a consistent sequence of corrections
to the mean-field equations using a posteriori knowledge from simulations.
These corrections both lead to better description of the dynamics and connect
the patterns to it. The analysis is readily applicable to realistic systems,
which we demonstrate by an example using an empirical metapopulation landscape.


  The origin of stochastic fluctuations in gene expression has received
considerable attention recently. Fluctuations in gene expression are
particularly pronounced in cellular systems because of the small copy number of
species undergoing transitions between discrete chemical states and the small
size of biological compartments. In this paper, we propose a stochastic model
for gene expression regulation including several binding sites, considering
elementary reactions only. The model is used to investigate the role of
cooperativity on the intrinsic fluctuations of gene expression, by means of
master equation formalism. We found that the Hill coefficient and the level of
noise increases as the interaction energy between activators increases.
Additionally, we show that the model allows to distinguish between two
cooperative binding mechanisms.


  Metabolic reactions of single-cell organisms are routinely observed to become
dispensable or even incapable of carrying activity under certain circumstances.
Yet, the mechanisms as well as the range of conditions and phenotypes
associated with this behavior remain very poorly understood. Here we predict
computationally and analytically that any organism evolving to maximize growth
rate, ATP production, or any other linear function of metabolic fluxes tends to
significantly reduce the number of active metabolic reactions compared to
typical non-optimal states. The reduced number appears to be constant across
the microbial species studied and just slightly larger than the minimum number
required for the organism to grow at all. We show that this massive spontaneous
reaction silencing is triggered by the irreversibility of a large fraction of
the metabolic reactions and propagates through the network as a cascade of
inactivity. Our results help explain existing experimental data on
intracellular flux measurements and the usage of latent pathways, shedding new
light on microbial evolution, robustness, and versatility for the execution of
specific biochemical tasks. In particular, the identification of optimal
reaction activity provides rigorous ground for an intriguing knockout-based
method recently proposed for the synthetic recovery of metabolic function.


  We consider the relation between the microscopic and effective descriptions
of the unfolding experiment on a model polypeptide. We evaluate the probability
distribution function of the performed work by Monte Carlo simulations and
compare it with that obtained by evaluating the work distribution generating
function on an effective Brownian motion model tailored to reproduce exactly
the equilibrium properties. The agreement is satisfactory for fast protocols,
but deteriorates for slower ones, hinting at the existence of processes on
several time scales even in such a simple system.


  Multiple Sequence Alignment (MSA) is one of the most computationally
intensive tasks in Computational Biology. Existing best known solutions for
multiple sequence alignment take several hours (in some cases days) of
computation time to align, for example, 2000 homologous sequences of average
length 300. Inspired by the Sample Sort approach in parallel processing, in
this paper we propose a highly scalable multiprocessor solution for the MSA
problem in phylogenetically diverse sequences. Our method employs an
intelligent scheme to partition the set of sequences into smaller subsets using
kmer count based similarity index, referred to as k-mer rank. Each subset is
then independently aligned in parallel using any sequential approach. Further
fine tuning of the local alignments is achieved using constraints derived from
a global ancestor of the entire set. The proposed Sample-Align-D Algorithm has
been implemented on a cluster of workstations using MPI message passing
library. The accuracy of the proposed solution has been tested on standard
benchmarks such as PREFAB. The accuracy of the alignment produced by our
methods is comparable to that of well known sequential MSA techniques. We were
able to align 2000 randomly selected sequences from the Methanosarcina
acetivorans genome in less than 10 minutes using Sample-Align-D on a 16 node
cluster, compared to over 23 hours on sequential MUSCLE system running on a
single cluster node.


  An overview of current multiple alignment systems to date are described.The
useful algorithms, the procedures adopted and their limitations are
presented.We also present the quality of the alignments obtained and in which
cases(kind of alignments, kind of sequences etc) the particular systems are
useful.


  Pyro-Align is a multiple alignment program specifically designed for
pyrosequencing reads of huge number. Multiple sequence alignment is shown to be
NP-hard and heuristics are designed for approximate solutions. Multiple
sequence alignment of pyrosequenceing reads is complex mainly because of 2
factors. One being the huge number of reads, making the use of traditional
heuristics,that scale very poorly for large number, unsuitable. The second
reason is that the alignment cannot be performed arbitrarily, because the
position of the reads with respect to the original genome is important and has
to be taken into account.In this report we present a short description of the
multiple alignment system for pyrosequencing reads.


  Pyrosequencing is among the emerging sequencing techniques, capable of
generating upto 100,000 overlapping reads in a single run. This technique is
much faster and cheaper than the existing state of the art sequencing technique
such as Sanger. However, the reads generated by pyrosequencing are short in
size and contain numerous errors. Furthermore, each read has a specific
position in the reference genome. In order to use these reads for any
subsequent analysis, the reads must be aligned . Existing multiple sequence
alignment methods cannot be used as they do not take into account the specific
positions of the sequences with respect to the genome, and are highly
inefficient for large number of sequences. Therefore, the common practice has
been to use either simple pairwise alignment despite its poor accuracy for
error prone pyroreads, or use computationally expensive techniques based on
sequential gap propagation. In this paper, we develop a computationally
efficient method based on domain decomposition, referred to as pyro-align, to
align such large number of reads. The proposed alignment algorithm accurately
aligns the erroneous reads in a short period of time, which is orders of
magnitude faster than any existing method. The accuracy of the alignment is
confirmed from the consensus obtained from the multiple alignments.


  We analyse sequential Markov coalescent algorithms for populations with
demographic structure: for a bottleneck model, a population-divergence model,
and for a two-island model with migration. The sequential Markov coalescent
method is an approximation to the coalescent suggested by McVean and Cardin,
and Marjoram and Wall. Within this algorithm we compute, for two individuals
randomly sampled from the population, the correlation between times to the most
recent common ancestor and the linkage probability corresponding to two
different loci with recombination rate R between them. We find that the
sequential Markov coalescent method approximates the coalescent well in general
in models with demographic structure. An exception is the case where
individuals are sampled from populations separated by reduced gene flow. In
this situation, the gene-history correlations may be significantly
underestimated. We explain why this is the case.