  Apoptosis is an important area of research because of its role in keeping a
mature multicellular organism's number of cells constant hence, ensuring that
the organism does not have cell accumulation that may transform into cancer
with additional hallmarks. Firstly, we have carried out sensitivity analysis on
an existing mitochondria-dependent mathematical apoptosis model to find out
which parameters have a role in causing monostable cell survival i.e.,
malfunction in apoptosis. We have then generated three healthy cell models by
changing these sensitive parameters while preserving bistability i.e., healthy
functioning. For each healthy cell, we varied the proapoptotic production
rates, which were found to be among the most sensitive parameters, to yield
cells that have malfunctioning apoptosis. We simulated caspase-3 activation, by
numerically integrating the governing ordinary differential equations of a
mitochondria-dependent apoptosis model, in a hypothetical malfunctioning cell
which is treated by four potential treatments, namely: (i) proteasome inhibitor
treatment, (ii) Bcl-2 inhibitor treatment, (iii) IAP inhibitor treatment, (iv)
Bid-like synthetic peptides treatment. The simulations of the present model
suggest that proteasome inhibitor treatment is the most effective treatment
though it may have severe side effects. For this treatment, we observed that
the amount of proteasome inhibitor needed for caspase-3 activation may be
different for cells in individuals with a different proapoptotic protein
deficiency. We also observed that caspase-3 can be activated by Bcl-2 inhibitor
treatment only in those hypothetical malfunctioning cells with Bax deficiency
but not in others. These support the view that molecular heterogeneity in
individuals may be an important factor in determining the individuals' positive
or negative responses to treatments.


  Non-genetic perturbations, such as environmental change or developmental
noise, can induce novel phenotypes. If an induced phenotype confers a fitness
advantage, selection may promote its genetic stabilization. Non-genetic
perturbations can thus initiate evolutionary innovation. Genetic variation that
is not usually phenotypically visible may play an important role in this
process. Populations under stabilizing selection on a phenotype that is robust
to mutations can accumulate such variation. After non-genetic perturbations,
this variation can become a source of new phenotypes. We here study the
relationship between a phenotype's robustness to mutations and a population's
potential to generate novel phenotypic variation. To this end, we use a
well-studied model of transcriptional regulation circuits. Such circuits are
important in many evolutionary innovations. We find that phenotypic robustness
promotes phenotypic variability in response to non-genetic perturbations, but
not in response to mutation. Our work suggests that non-genetic perturbations
may initiate innovation more frequently in mutationally robust gene expression
traits.


  The immune system recognizes a myriad of invading pathogens and their toxic
products. It does so with a finite repertoire of antibodies and T cell
receptors. We here describe theories that quantify the immune system dynamics.
We describe how the immune system recognizes antigens by searching the large
space of receptor molecules. We consider in some detail the theories that
quantify the immune response to influenza and dengue fever. We review
theoretical descriptions of the complementary evolution of pathogens that
occurs in response to immune system pressure. Methods including bioinformatics,
molecular simulation, random energy models, and quantum field theory contribute
to a theoretical understanding of aspects of immunity.


  The normalized maximum likelihood (NML) is a recent penalized likelihood that
has properties that justify defining the amount of discrimination information
(DI) in the data supporting an alternative hypothesis over a null hypothesis as
the logarithm of an NML ratio, namely, the alternative hypothesis NML divided
by the null hypothesis NML. The resulting DI, like the Bayes factor but unlike
the p-value, measures the strength of evidence for an alternative hypothesis
over a null hypothesis such that the probability of misleading evidence
vanishes asymptotically under weak regularity conditions and such that evidence
can support a simple null hypothesis. Unlike the Bayes factor, the DI does not
require a prior distribution and is minimax optimal in a sense that does not
involve averaging over outcomes that did not occur. Replacing a (possibly
pseudo-) likelihood function with its weighted counterpart extends the scope of
the DI to models for which the unweighted NML is undefined. The likelihood
weights leverage side information, either in data associated with comparisons
other than the comparison at hand or in the parameter value of a simple null
hypothesis. Two case studies, one involving multiple populations and the other
involving multiple biological features, indicate that the DI is robust to the
type of side information used when that information is assigned the weight of a
single observation. Such robustness suggests that very little adjustment for
multiple comparisons is warranted if the sample size is at least moderate.


  Transcription factors (TF) play an essential role in the cell as locus- and
condition-specific recruiters of transcriptional machinery or
chromatin-modifying complexes. However, predicting the in vivo profile of TF
occupancy along the genome, which depends on complex interactions with other
chromatin-associated proteins, from the DNA sequence remains a major challenge.
Through careful reanalysis of ChIP-chip data for 138 TFs obtained in rich
media, we were able to classify the upstream promoter regions of S. cerevisiae
into 15 distinct chromatin types. One of these encompasses 5% of all promoters
and is unique in that it is highly occupied by (essentially) all TFs expressed
in rich media. These "hotspots" of TF occupancy are strongly
nucleosome-depleted and preferentially targeted by chromatin-remodeling
complexes and the origin-of-replication complex (ORC). They are also the only
chromatin type enriched for predicted Rap1p and Pdr1p binding sites, which we
found to work cooperatively with AAA/TTT motifs, known to affect local DNA
structure, to reduce nucleosome occupancy. Taken together, our results reveal
and characterize a new type of local chromatin structure in yeast.


  NF-{\kappa}B is a pleiotropic protein whose nucleo-cytoplasmic trafficking is
tightly regulated by multiple negative feedback loops embedded in the
NF-{\kappa}B signaling network and contributes to diverse gene expression
profiles important in immune cell differentiation, cell apoptosis, and innate
immunity. The intracellular signaling processes and their control mechanisms,
however, are susceptible to both extrinsic and intrinsic noise. In this
article, we present numerical evidence for a universal dynamic behavior of
NF-{\kappa}B, namely oscillatory nucleo-cytoplasmic shuttling, due to the
fundamentally stochastic nature of the NF-{\kappa}B signaling network. We
simulated the effect of extrinsic noise with a deterministic ODE model, using a
statistical ensemble approach, generating many copies of the signaling network
with different kinetic rates sampled from a biologically feasible parameter
space. We modeled the effect of intrinsic noise by simulating the same networks
stochastically using the Gillespie algorithm. The results demonstrate that
extrinsic noise diversifies the shuttling patterns of NF-{\kappa}B response,
whereas intrinsic noise induces oscillatory behavior in many of the otherwise
non-oscillatory patterns. We identify two key model parameters which
significantly affect the NF-{\kappa}B dynamic response and deduce a
two-dimensional phase-diagram of the NF-{\kappa}B response as a function of
these parameters. We conclude that if single-cell experiments are performed, a
rich variety of NF-{\kappa}B response will be observed, even if
population-level experiments, which average response over large numbers of
cells, do not evidence oscillatory behavior.


  Cellular responses in the single cells are known to be highly heterogeneous
and individualistic due to the strong influence by extrinsic and intrinsic
noise. Here, we are concerned about how to model the extrinsic noise-induced
heterogeneous response in the single cells under the constraints of
experimentally obtained population-averaged response, but without much detailed
kinetic information. We propose a novel statistical ensemble scheme where
extrinsic noise is regarded as fluctuations in the values of kinetic parameters
and such fluctuations are modeled by randomly sampling the kinetic rate
constants from a uniform distribution. We consider a large number of signaling
system replicates, each of which has the same network topology, but a uniquely
different set of kinetic rate constants. A protein dynamic response from each
replicate should represent the dynamics in a single cell and the statistical
ensemble average should be regarded as a population-level response averaged
over a population of the cells. We devise an optimization algorithm to find the
correct uniform distribution of the network parameters, which produces the
correct statistical distribution of the response whose ensemble average and
distribution agree well with the population-level experimental data and the
experimentally observed heterogeneity. We apply this statistical ensemble
analysis to a NF-{\kappa}B signaling system and (1) predict the distributions
of the heterogeneous NF-{\kappa}B (either oscillatory or non-oscillatory)
dynamic patterns and of the dynamic features (e.g., period), (2) predict that
both the distribution and the statistical ensemble average of the NF-{\kappa}B
dynamic response depends sensitively on the dosage of stimulant, and lastly (3)
demonstrate the sigmoidally shaped dose-response from the statistical ensemble
average and the individual replicates.


  Essential oils have antifungal and antipathogenic effects and therefore are
targets in plant pathology research for their potential uses as natural
substitutes for inorganic plant pesticides. Beauveria bassiana, an
entomopathogenic fungus, can endophytically colonize a vast number of plant
species and trigger induced systemic resistance against plant pathogens. Spore
germination is the most vulnerable in the fungal life cycle and is therefore a
good candidate for monitoring the effect of essential oils on the growth of B.
bassiana. Percentage germination of fungal spores and length of germination
tubes were recorded from experiments. A mathematical model that was able to
capture the effects of cymene, an essential oil produced by Monarda, on the
germination and growth was developed. This is the first report of a model for
the impact of essential oils on B. bassiana spore germination.


  Metabolic networks are complex systems that comprise hundreds of chemical
reactions which synthesize biomass molecules from chemicals in an organism's
environment. The metabolic network of any one organism is encoded by a
metabolic genotype, defined by a set of enzyme-coding genes whose products
catalyze the network's reactions. Each metabolic genotype has a metabolic
phenotype, such as the ability to synthesize biomass on a spectrum of different
sources of chemical elements and energy. We here focus on sulfur metabolism,
which is attractive to study the evolution of metabolic networks, because it
involves many fewer reactions than carbon metabolism. Specifically, we study
properties of the space of all possible metabolic genotypes, and analyze
properties of random metabolic genotypes that are viable on different numbers
of sulfur sources. We show that metabolic genotypes with the same phenotype
form large connected genotype networks that extend far through metabolic
genotype space. How far they reach through this space is a linear function of
the number of super-essential reactions in such networks, the number of
reactions that occur in all networks with the same phenotype. We show that
different neighborhoods of any genotype network harbor very different novel
phenotypes, metabolic innovations that can sustain life on novel sulfur
sources. We also analyze the ability of evolving populations of metabolic
networks to explore novel metabolic phenotypes. This ability is facilitated by
the existence of genotype networks, because different neighborhoods of these
networks contain very different novel phenotypes. In contrast to
macromolecules, where phenotypic robustness may facilitate phenotypic
innovation, we show that here the ability to access novel phenotypes does not
monotonically increase with robustness.


  This paper introduces two new probabilistic graphical models for
reconstruction of genetic regulatory networks using DNA microarray data. One is
an Independence Graph (IG) model with either a forward or a backward search
algorithm and the other one is a Gaussian Network (GN) model with a novel
greedy search method. The performances of both models were evaluated on four
MAPK pathways in yeast and three simulated data sets. Generally, an IG model
provides a sparse graph but a GN model produces a dense graph where more
information about gene-gene interactions is preserved. Additionally, we found
two key limitations in the prediction of genetic regulatory networks using DNA
microarray data, the first is the sufficiency of sample size and the second is
the complexity of network structures may not be captured without additional
data at the protein level. Those limitations are present in all prediction
methods which used only DNA microarray data.


  By integrating heterogeneous functional genomic datasets, we have developed a
new framework for detecting combinatorial control of gene expression, which
includes estimating transcription factor activities using a singular value
decomposition method and reducing high-dimensional input gene space by
considering genomic properties of gene clusters. The prediction of cooperative
gene regulation is accomplished by either Gaussian Graphical Models or Pairwise
Mixed Graphical Models. The proposed framework was tested on yeast cell cycle
datasets: (1) 54 known yeast cell cycle genes with 9 cell cycle regulators and
(2) 676 putative yeast cell cycle genes with 9 cell cycle regulators. The new
framework gave promising results on inferring TF-TF and TF-gene interactions.
It also revealed several interesting mechanisms such as negatively correlated
protein-protein interactions and low affinity protein-DNA interactions that may
be important during the yeast cell cycle. The new framework may easily be
extended to study other higher eukaryotes.


  Traditional approaches to ecosystem modelling have relied on spatially
homogeneous approximations to interaction, growth and death. More recently,
spatial interaction and dispersal have also been considered. While these leads
to certain changes in community dynamics, their effect is sometimes fairly
minimal, and demographic scenarios in which this difference is important have
not been systematically investigated.
  We take a simple mean-field model which simulates birth, growth and death
processes, and rewrite it with spatially distributed discrete individuals. Each
individual's growth and mortality is determined by a competition measure which
captures the effects of neighbours in a way which retains the conceptual
simplicity of a generic, analytically-solvable model. Although the model is
generic, we here parameterise it using data from Caledonian Scots Pine stands.
The dynamics of simulated populations, starting from a plantation lattice
configuration, mirror those of well-established qualitative descriptions of
natural forest stand behaviour; an analogy which assists in understanding the
transition from artificial to old-growth structure.
  When parameterised for Scots Pine populations, the signature of spatial
processes is evident, but they do not have a large effect on first-order
statistics such as density and biomass. The sensitivity of this result to
variation in each individual rate parameter is investigated; distinct
differences between spatial and mean-field models are seen only upon alteration
of the interaction strength parameters, and in low density populations. Under
the Scots Pine parameterisation, dispersal also has an effect of spatial
structure, but not first-order properties. Only in more intense competitive
scenarios does altering the relative scales of dispersal and interaction lead
to a clear signal in first order behaviour.


  The influence of DNA cis-regulatory elements on a gene's expression has been
intensively studied. However, little is known about expressions driven by
trans-acting DNA hotspots. DNA hotspots harboring copy number aberrations are
recognized to be important in cancer as they influence multiple genes on a
global scale. The challenge in detecting trans-effects is mainly due to the
computational difficulty in detecting weak and sparse trans-acting signals
amidst co-occuring passenger events. We propose an integrative approach to
learn a sparse interaction network of DNA copy-number regions with their
downstream targets in a breast cancer dataset. Information from this network
helps distinguish copy-number driven from copy-number independent expression
changes on a global scale. Our result further delineates cis- and trans-effects
in a breast cancer dataset, for which important oncogenes such as ESR1 and
ERBB2 appear to be highly copy-number dependent. Further, our model is shown to
be efficient and in terms of goodness of fit no worse than other state-of the
art predictors and network reconstruction models using both simulated and real
data.


  From genetic regulatory networks to nervous systems, the interactions between
elements in biological networks often take a sigmoidal or S-shaped form. This
paper develops a probabilistic characterization of the parameter space of
continuous-time sigmoidal networks (CTSNs), a simple but dynamically-universal
model of such interactions. We describe an efficient and accurate method for
calculating the probability of observing effectively M-dimensional dynamics in
an N-element CTSN, as well as a closed-form but approximate method. We then
study the dependence of this probability on N, M, and the parameter ranges over
which sampling occurs. This analysis provides insight into the overall
structure of CTSN parameter space.


  NF-{\kappa}B, one of key regulators of inflammation, apoptosis, and
differentiation, was found to have noisy oscillatory shuttling between the
nucleus and the cytoplasm in single cells when cells are stimulated by cytokine
TNF{\alpha}. We present the analytical analysis which uncovers the underlying
physical mechanisms of this spectacular noise-induced transition in biological
networks. Starting with the master equation describing both signaling and
transcription events in NF-{\kappa}B signaling network, we derived the
macroscopic and the Fokker-Planck equations by using van Kampen's sysem size
expansion. Using the noise-induced oscillatory signatures present in the power
spectrum, we constructed the two-dimensional phase diagram where the
noise-induced oscillation emerges in the dynamically stable parameter space.


  The NF-{\kappa}B signaling network plays an important role in many different
compartments of the immune system during immune activation. Using a
computational model of the NF-{\kappa}B signaling network involving two
negative regulators, I{\kappa}B{\alpha} and A20, we performed sensitivity
analyses with three different sampling methods and present a ranking of the
kinetic rate variables by the strength of their influence on the NF-{\kappa}B
signaling response. We also present a classification of temporal response
profiles of nuclear NF-{\kappa}B concentration into six clusters, which can be
regrouped to three biologically relevant clusters. Lastly, based upon the
ranking, we constructed a reduced network of the
IKK-NF-{\kappa}B-I{\kappa}B{\alpha}-A20 signal transduction.


  We review models and approaches for error-control in order to prevent the
buildup of noise when gates for digital chemical and biomolecular computing
based on (bio)chemical reaction processes are utilized to realize stable,
scalable networks for information processing. Solvable rate-equation models
illustrate several recently developed methodologies for gate-function
optimization. We also survey future challenges and possible new research
avenues.


  We present an algorithm for phylogenetic reconstruction using quartets that
returns the correct topology for $n$ taxa in $O(n \log n)$ time with high
probability, in a probabilistic model where a quartet is not consistent with
the true topology of the tree with constant probability, independent of other
quartets. Our incremental algorithm relies upon a search tree structure for the
phylogeny that is balanced, with high probability, no matter what the true
topology is. Our experimental results show that our method is comparable in
runtime to the fastest heuristics, while still offering consistency guarantees.


  Weak selection, which means a phenotype is slightly advantageous over
another, is an important limiting case in evolutionary biology. Recently it has
been introduced into evolutionary game theory. In evolutionary game dynamics,
the probability to be imitated or to reproduce depends on the performance in a
game. The influence of the game on the stochastic dynamics in finite
populations is governed by the intensity of selection. In many models of both
unstructured and structured populations, a key assumption allowing analytical
calculations is weak selection, which means that all individuals perform
approximately equally well. In the weak selection limit many different
microscopic evolutionary models have the same or similar properties. How
universal is weak selection for those microscopic evolutionary processes? We
answer this question by investigating the fixation probability and the average
fixation time not only up to linear, but also up to higher orders in selection
intensity. We find universal higher order expansions, which allow a rescaling
of the selection intensity. With this, we can identify specific models which
violate (linear) weak selection results, such as the one--third rule of
coordination games in finite but large populations.


This paper has been withdrawn by the author.


  Twisting DNA under a constant applied force reveals a thermally activated
transition into a state with a supercoiled structure known as a plectoneme.
Using transition state theory, we predict the rate of this plectoneme
nucleation to be of order 10^4 Hz. We reconcile this with experiments that have
measured hopping rates of order 10 Hz by noting that the viscosity of the bead
used to manipulate the DNA limits the measured rate. We find that the intrinsic
bending caused by disorder in the base-pair sequence is important for
understanding the free energy barrier that governs the transition. Both
analytic and numerical methods are used in the calculations. We provide
extensive details on the numerical methods for simulating the elastic rod model
with and without disorder.


  Over the past years several authors have used the approach of generalized
modeling to study the dynamics of food chains and food webs. Generalized models
come close to the efficiency of random matrix models, while being as directly
interpretable as conventional differential-equation-based models. Here we
present a pedagogical introduction to the approach of generalized modeling.
This introduction places more emphasis on the underlying concepts of
generalized modeling than previous publications. Moreover, we propose a
shortcut that can significantly accelerate the formulation of generalized
models and introduce an iterative procedure that can be used to refine existing
generalized models by integrating new biological insights.


  In the absence of selection, the structure of allelic diversity is described
by the elegant sampling formula of Ewens. This formula has helped shape our
expectations of empirical patterns of molecular variation. Along with
coalescent theory, it provides statistical techniques for rejecting the null
model of neutrality. However, we still do not fully understand the statistics
of the allelic diversity we expect to see in the presence of natural selection.
Earlier work has described the effects of strongly deleterious mutations linked
to many neutral sites, and allelic variation in models where offspring fitness
is unrelated to parental fitness, but it has proven difficult to understand
allelic diversity in the presence of purifying selection at many linked sites.
Here, we study the population genetics of infinitely many perfectly linked
sites, some neutral and some deleterious. Our approach is based on studying the
lineage structure within each class of individuals of similar fitness in the
deleterious mutation-selection balance. Analogous to the Ewens sampling
formula, we derive expressions for the likelihoods of any configuration of
allelic types in a sample. We find that for moderate and weak selection
pressures the patterns of allelic diversity cannot be described by a neutral
model for any choice of the effective population size, indicating that there is
power to detect selection from patterns of sampled allelic diversity.


  Compared to a neutral model, purifying selection distorts the structure of
genealogies and hence alters the patterns of sampled genetic variation.
Although these distortions may be common in nature, our understanding of how we
expect purifying selection to affect patterns of molecular variation remains
incomplete. Genealogical approaches such as coalescent theory have proven
difficult to generalize to situations involving selection at many linked sites,
unless selection pressures are extremely strong. Here, we introduce an
effective coalescent theory (a "fitness-class coalescent") to describe the
structure of genealogies in the presence of purifying selection at many linked
sites. We use this effective theory to calculate several simple statistics
describing the expected patterns of variation in sequence data, both at the
sites under selection and at linked neutral sites. Our analysis combines our
earlier description of the allele frequency spectrum in the presence of
purifying selection (Desai et al. 2010) with the structured coalescent approach
of Nordborg (1997), to trace the ancestry of individuals through the
distribution of fitnesses within the population. Alternatively, we can derive
our results using an extension of the coalescent approach of Hudson and Kaplan
(1994). We find that purifying selection leads to patterns of genetic variation
that are related but not identical to a neutrally evolving population in which
population size has varied in a specific way in the past.


  It is shown, by means of numerical simulations, that intercellular spiral
waves of calcium can be initiated in a network of coupled cells as a result of
a de-synchronization between Ca2+ oscillations in two domains. No artificial
heterogeneities need to be imposed to the system for spontaneous formation of
spiral waves. The de-synchronization occurs near the interface of the
stimulated region (which acts as a pacemaker) and propagates over the entire
network. We also find the outcome of the collision of two spiral waves.


  In the present paper we study a lattice model of two species competing for
the same resources. Monte Carlo simulations for d=1, 2, and 3 show that when
resources are easily available both species coexist. However, when the supply
of resources is on an intermediate level, the species with slower metabolism
becomes extinct. On the other hand, when resources are scarce it is the species
with faster metabolism that becomes extinct. The range of coexistence of the
two species increases with dimension. We suggest that our model might describe
some aspects of the competition between normal and tumor cells. With such an
interpretation, examples of tumor remission, recurrence and of different
morphologies are presented. In the d=1 and d=2 models, we analyse the nature of
phase transitions: they are either discontinuous or belong to the
directed-percolation universality class, and in some cases they have an active
subcritical phase. In the d=2 case, one of the transitions seems to be
characterized by critical exponents different than directed-percolation ones,
but this transition could be also weakly discontinuous. In the d=3 version,
Monte Carlo simulations are in a good agreement with the solution of the
mean-field approximation. This approximation predicts that oscillatory
behaviour occurs in the present model, but only for d>2. For d>=2, a steady
state depends on the initial configuration in some cases.


  We introduce the first analytical model of asymmetric community dynamics to
yield Hubbell's neutral theory in the limit of functional equivalence among all
species. Our focus centers on an asymmetric extension of Hubbell's local
community dynamics, while an analogous extension of Hubbell's metacommunity
dynamics is deferred to an appendix. We find that mass-effects may facilitate
coexistence in asymmetric local communities and generate unimodal species
abundance distributions indistinguishable from those of symmetric communities.
Multiple modes, however, only arise from asymmetric processes and provide a
strong indication of non-neutral dynamics. Although the exact stationary
distributions of fully asymmetric communities must be calculated numerically,
we derive approximate sampling distributions for the general case and for
nearly neutral communities where symmetry is broken by a single species
distinct from all others in ecological fitness and dispersal ability. In the
latter case, our approximate distributions are fully normalized, and novel
asymptotic expansions of the required hypergeometric functions are provided to
make evaluations tractable for large communities. Employing these results in a
Bayesian analysis may provide a novel statistical test to assess the
consistency of species abundance data with the neutral hypothesis.


  By exploiting an analogy between population genetics and statistical
mechanics, we study the evolution of a polygenic trait under stabilizing
selection, mutation, and genetic drift. This requires us to track only four
macroscopic variables, instead of the distribution of all the allele
frequencies that influence the trait. These macroscopic variables are the
expectations of: the trait mean and its square, the genetic variance, and of a
measure of heterozygosity, and are derived from a generating function that is
in turn derived by maximizing an entropy measure. These four macroscopics are
enough to accurately describe the dynamics of the trait mean and of its genetic
variance (and in principle of any other quantity). Unlike previous approaches
that were based on an infinite series of moments or cumulants, which had to be
truncated arbitrarily, our calculations provide a well-defined approximation
procedure. We apply the framework to abrupt and gradual changes in the optimum,
as well as to changes in the strength of stabilizing selection. Our
approximations are surprisingly accurate, even for systems with as few as 5
loci. We find that when the effects of drift are included, the expected genetic
variance is hardly altered by directional selection, even though it fluctuates
in any particular instance. We also find hysteresis, showing that even after
averaging over the microscopic variables, the macroscopic trajectories retain a
memory of the underlying genetic states.


  We consider several one-species population dynamics model with finite and
infinite carrying capacity, time dependent growth and effort rates and solve
them analytically. We show that defining suitable scaling functions for a given
time, one is able to demonstrate that their ratio with respect to its initial
value is universal. This ratio is independent from the initial condition and
from the model parameters. Although the effort rate does not break the model
universality it produces a transition between the species extinction and
survival. A general formula is furnished to obtain the scaling functions.


  Conformational changes upon protein-protein association are the key element
of the binding mechanism. The study presents a systematic large-scale analysis
of such conformational changes in the side chains. The results indicate that
short and long side chains have different propensities for the conformational
changes. Long side chains with three or more dihedral angles are often subject
to large conformational transition. Shorter residues with one or two dihedral
angles typically undergo local conformational changes not leading to a
conformational transition. The relationship between the local readjustments and
the equilibrium fluctuations of a side chain around its unbound conformation is
suggested. Most of the side chains undergo larger changes in the dihedral angle
most distant from the backbone. The amino acids with symmetric aromatic (Phe
and Tyr) and charged (Asp and Glu) groups show the opposite trend where the
near-backbone dihedral angles change the most. The frequencies of the
core-to-surface interface transitions of six nonpolar residues and Tyr exceed
the frequencies of the opposite, surface-to-core transitions. The binding
increases both polar and nonpolar interface areas. However, the increase of the
nonpolar area is larger for all considered classes of protein complexes. The
results suggest that the protein association perturbs the unbound interfaces to
increase the hydrophobic forces. The results facilitate better understanding of
the conformational changes in proteins and suggest directions for efficient
conformational sampling in docking protocols.


  Disentangling the processes leading populations to extinction is a major
topic in ecology and conservation biology. The difficulty to find a mate in
many species is one of these processes. Here, we investigate the impact of
self-incompatibility in flowering plants, where several inter-compatible
classes of individuals exist but individuals of the same class cannot mate. We
model pollen limitation through different relationships between mate
availability and fertilization success. After deriving a general stochastic
model, we focus on the simple case of distylous plant species where only two
classes of individuals exist. We first study the dynamics of such a species in
a large population limit and then, we look for an approximation of the
extinction probability in small populations. This leads us to consider
inhomogeneous random walks on the positive quadrant. We compare the dynamics of
distylous species to self-fertile species with and without inbreeding
depression, to obtain the conditions under which self-incompatible species
could be less sensitive to extinction while they can suffer more pollen
limitation.


  Transcription factors (TFs) are key regulators of gene expression. Based on
the classical scenario in which the TF search process switches between
one-dimensional motion along the DNA molecule and free Brownian motion in the
nucleus, we study the arrival time of several TFs to multiple binding sites and
derive, in the presence of competitive binding ligands, the probability that
several target sites are bound. We then apply our results to the hunchback
regulation by bicoid in the fly embryo and we propose a general mechanism that
allows cells to read a morphogenetic gradient and specialize according to their
position in the embryo.


  Repulsive point processes arise in models where competition forces entities
to be more spread apart than if placed independently. Simulation of these types
of processes can be accomplished using dominated coupling from the past with a
running time that varies as the intensity of the number of points. These
algorithms usually exhibit what is called an artificial phase transition, where
below a critical intensity the algorithm runs in finite expected time, but
above the critical intensity the expected number of steps is infinite. Here the
artificial phase transition is examined. In particular, an earlier lower bound
on this artificial phase transition is improved by including a new type of term
in the analysis. In addition, the results of computer experiments to locate the
transition are presented.


  Using Monte-Carlo simulations, we deconvolved the sliding and hopping
kinetics of GFP-LacI proteins on elongated DNA from their experimentally
observed seconds-long diffusion trajectories. Our simulations suggest the
following results: (1) in each diffusion trajectory, a protein makes on average
hundreds of alternating slides and hops with a mean sliding time of several
tens of ms; (2) sliding dominates the root mean square displacement of fast
diffusion trajectories, whereas hopping dominates slow ones; (3) flow and
variations in salt concentration have limited effects on hopping kinetics,
while in vivo DNA configuration is not expected to influence sliding kinetics;
furthermore, (4) the rate of occurrence for hops longer than 200 nm agrees with
experimental data for EcoRV proteins.


  Embryonic stem cells (ESC) have the potential to self-renew indefinitely and
to differentiate into any of the three germ layers. The molecular mechanisms
for self-renewal, maintenance of pluripotency and lineage specification are
poorly understood, but recent results point to a key role for epigenetic
mechanisms. In this study, we focus on quantifying the impact of histone 3
acetylation (H3K9,14ac) on gene expression in murine embryonic stem cells. We
analyze genome-wide histone acetylation patterns and gene expression profiles
measured over the first five days of cell differentiation triggered by
silencing Nanog, a key transcription factor in ESC regulation. We explore the
temporal and spatial dynamics of histone acetylation data and its correlation
with gene expression using supervised and unsupervised statistical models. On a
genome-wide scale, changes in acetylation are significantly correlated to
changes in mRNA expression and, surprisingly, this coherence increases over
time. We quantify the predictive power of histone acetylation for gene
expression changes in a balanced cross-validation procedure. In an in-depth
study we focus on genes central to the regulatory network of Mouse ESC,
including those identified in a recent genome-wide RNAi screen and in the
PluriNet, a computationally derived stem cell signature. We find that compared
to the rest of the genome, ESC-specific genes show significantly more
acetylation signal and a much stronger decrease in acetylation over time, which
is often not reflected in an concordant expression change. These results shed
light on the complexity of the relationship between histone acetylation and
gene expression and are a step forward to dissect the multilayer regulatory
mechanisms that determine stem cell fate.


  The process of protein folding from an unfolded state to a biologically
active, folded conformation is governed by many parameters e.g the sequence of
amino acids, intermolecular interactions, the solvent, temperature and chaperon
molecules. Our study, based on random matrix modeling of the interactions,
shows however that the evolution of the statistical measures e.g Gibbs free
energy, heat capacity, entropy is single parametric. The information can
explain the selection of specific folding pathways from an infinite number of
possible ways as well as other folding characteristics observed in computer
simulation studies.


  The two-species population dynamics model is the simplest paradigm of
interspecies interaction. Here, we include intraspecific competition to the
Lotka-Volterra model and solve it analytically. Despite being simple and
thoroughly studied, this model presents a very rich behavior and some
characteristics not so well explored, which are unveiled. The forbidden region
in the mutualism regime and the dependence on initial conditions in the
competition regime are some examples of these characteristics. From the
stability of the steady state solutions, three phases are obtained: (i)
extinction of one species (Gause transition), (ii) their coexistence and (iii)
a forbidden region. Full analytical solutions have been obtained for the
considered ecological regimes. The time transient allows one to defined time
scales for the system evolution, which can be relevant for the study of tumor
growth by theoretical or computer simulation models.


  The capacity to resist attacks from the environment is crucial to the
survival of all organisms. We quantitatively analyze the susceptibility of
protein interaction networks of numerous organisms to random and malicious
attacks. We find for all organisms studied that random rewiring improves
protein network robustness, so that actual networks are more fragile than
rewired surrogates. This unexpected fragility contrasts with the behavior of
networks such as the Internet, whose robustness decreases with random rewiring.
We trace this surprising effect to the modular structure of protein networks.


  Several pathogens use evolvability as a survival strategy against acquired
immunity of the host. Despite their high variability in time, some of them
exhibit quite low variability within the population at any given time, a
somehow paradoxical behavior often called the evolving quasispecies. In this
paper we introduce a simplified model of an evolving viral population in which
the effects of the acquired immunity of the host are represented by the
decrease of the fitness of the corresponding viral strains, depending on the
frequency of the strain in the viral population. The model exhibits evolving
quasispecies behavior in a certain range of its parameters, ans suggests how
punctuated evolution can be induced by a simple feedback mechanism.


  Dynamical networks with time delays can pose a considerable challenge for
mathematical analysis. Here, we extend the approach of generalized modeling to
investigate the stability of large networks of delay-coupled delay oscillators.
When the local dynamical stability of the network is plotted as a function of
the two delays then a pattern of tongues is revealed. Exploiting a link between
structure and dynamics, we identify conditions under which perturbations of the
topology have a strong impact on the stability. If these critical regions are
avoided the local stability of large random networks can be well approximated
analytically.


  We study a circular order on labelled, m-edge-coloured trees with k vertices,
and show that the set of such trees with a fixed circular order is in bijection
with the set of RNA m-diagrams of degree k, combinatorial objects which can be
regarded as RNA secondary structures of a certain kind. We enumerate these sets
and show that the set of trees with a fixed circular order can be characterized
as an equivalence class for the transitive closure of an operation which, in
the case m=3, arises as an induction in the context of interval exchange
transformations.


  Existing theoretical models of evolution focus on the relative fitness
advantages of different mutants in a population while the dynamic behavior of
the population size is mostly left unconsidered. We here present a generic
stochastic model which combines the growth dynamics of the population and its
internal evolution. Our model thereby accounts for the fact that both
evolutionary and growth dynamics are based on individual reproduction events
and hence are highly coupled and stochastic in nature. We exemplify our
approach by studying the dilemma of cooperation in growing populations and show
that genuinely stochastic events can ease the dilemma by leading to a transient
but robust increase in cooperation


  Microphase separation of membrane components is thought to play an important
role in many physiological processes, from cell signaling to endocytosis and
cellular trafficking. Here, we study how variations in the membrane composition
can be driven by fluctuating forces. We show that the membrane steady state is
not only controlled by the strength of the forces and how they couple to the
membrane, but also by their dynamics: In a simple class of models this is
captured by a single a correlation time. We conclude that the coupling of
membrane composition to normal mechanical forces, such as might be exerted by
polymerizing cytoskeleton filaments, could play an important role in
controlling the steady state of a cell membrane that exhibits transient
microphase separation on lengthscales in the 10-100 nm regime.


  We have performed Molecular Dynamics simulations of ectoine, hydroxyectoine
and urea in explicit solvent. Special attention has been spent on the local
surrounding structure of water molecules. Our results indicate that ectoine and
hydroxyectoine are able to accumulate more water molecules than urea by a
pronounced ordering due to hydrogen bonds. We have validated that the charging
of the molecules is of main importance resulting in a well defined hydration
sphere. The influence of a varying salt concentration is also investigated.
Finally we present experimental results of a DPPC monolayer phase transition
that validate our numerical findings.


  We propose an exactly solvable simplified statistical mechanical model for
the thermodynamics of beta-amyloid aggregation, generalizing a well-studied
model for protein folding. The monomer concentration is explicitly taken into
account as well as a non trivial dependence on the microscopic degrees of
freedom of the single peptide chain, both in the alpha-helix folded isolated
state and in the fibrillar one. The phase diagram of the model is studied and
compared to the outcome of fibril formation experiments which is qualitatively
reproduced.


  Reversibility, weak reversibility and deficiency, detailed and complex
balancing are generally not "encoded" in the kinetic differential equations but
they are realization properties that may imply local or even global asymptotic
stability of the underlying reaction kinetic system when further conditions are
also fulfilled. In this paper, efficient numerical procedures are given for
finding complex balanced or detailed balanced realizations of mass action type
chemical reaction networks or kinetic dynamical systems in the framework of
linear programming. The procedures are illustrated on numerical examples.


  We present TT2NE, a new algorithm to predict RNA secondary structures with
pseudoknots. The method is based on a classification of RNA structures
according to their topological genus. TT2NE guarantees to find the minimum free
energy structure irrespectively of pseudoknot topology. This unique proficiency
is obtained at the expense of the maximum length of sequence that can be
treated but comparison with state-of-the-art algorithms shows that TT2NE is a
very powerful tool within its limits. Analysis of TT2NE's wrong predictions
sheds light on the need to study how sterical constraints limit the range of
pseudoknotted structures that can be formed from a given sequence. An
implementation of TT2NE on a public server can be found at
http://ipht.cea.fr/rna/tt2ne.php.


  I investigate a disordered version of a simplified model of protein folding,
with binary degrees of freedom, applied to an ideal beta-hairpin structure.
Disorder is introduced by assuming that the contact energies are independent
and identically distributed random variables. The equilibrium free-energy of
the model is studied, performing the exact calculation of its quenched value
and proving the self-averaging feature.


  We study the growth of a population of bacteria in a dynamical hostile
environment corresponding to the immune system of the colonised organism. The
immune cells evolve as subcritical open clusters of oriented percolation and
are perpetually reinforced by an immigration process, while the bacteria try to
grow as a supercritical oriented percolation in the remaining empty space. For
appropriate values of the parameters, we prove that the population of bacteria
grows linearly. In this perspective, we build general tools to study dependent
percolation models issued from renormalization processes.


  We consider the unwinding of two lattice polymer strands of length N that are
initially wound around each other in a double-helical conformation and evolve
through Rouse dynamics. The problem relates to quickly bringing a
double-stranded polymer well above its melting temperature, i.e., the binding
interactions between the strands are neglected, and the strands separate from
each other as it is entropically favorable for them to do so. The strands
unwind by rotating around each other until they separate. We find that the
process proceeds from the ends inward; intermediate conformations can be
characterized by a tightly wound inner part, from which loose strands are
sticking out, with length l~t^0.39. The total time needed for the two strands
to unwind scales as a power of N as tu~N^(2.57+-0.03). We present a theoretical
argument, which suggests that during this unwinding process, these loose
strands are far out of equilibrium.


  Because of the high cost of commercial genotyping chip technologies, many
investigations have used a two-stage design for genome-wide association
studies, using part of the sample for an initial discovery of ``promising''
SNPs at a less stringent significance level and the remainder in a joint
analysis of just these SNPs using custom genotyping. Typical cost savings of
about 50% are possible with this design to obtain comparable levels of overall
type I error and power by using about half the sample for stage I and carrying
about 0.1% of SNPs forward to the second stage, the optimal design depending
primarily upon the ratio of costs per genotype for stages I and II. However,
with the rapidly declining costs of the commercial panels, the generally low
observed ORs of current studies, and many studies aiming to test multiple
hypotheses and multiple endpoints, many investigators are abandoning the
two-stage design in favor of simply genotyping all available subjects using a
standard high-density panel. Concern is sometimes raised about the absence of a
``replication'' panel in this approach, as required by some high-profile
journals, but it must be appreciated that the two-stage design is not a
discovery/replication design but simply a more efficient design for discovery
using a joint analysis of the data from both stages. Once a subset of
highly-significant associations has been discovered, a truly independent
``exact replication'' study is needed in a similar population of the same
promising SNPs using similar methods.


  The standard paradigm for the analysis of genome-wide association studies
involves carrying out association tests at both typed and imputed SNPs. These
methods will not be optimal for detecting the signal of association at SNPs
that are not currently known or in regions where allelic heterogeneity occurs.
We propose a novel association test, complementary to the SNP-based approaches,
that attempts to extract further signals of association by explicitly modeling
and estimating both unknown SNPs and allelic heterogeneity at a locus. At each
site we estimate the genealogy of the case-control sample by taking advantage
of the HapMap haplotypes across the genome. Allelic heterogeneity is modeled by
allowing more than one mutation on the branches of the genealogy. Our use of
Bayesian methods allows us to assess directly the evidence for a causative SNP
not well correlated with known SNPs and for allelic heterogeneity at each
locus. Using simulated data and real data from the WTCCC project, we show that
our method (i) produces a significant boost in signal and accurately identifies
the form of the allelic heterogeneity in regions where it is known to exist,
(ii) can suggest new signals that are not found by testing typed or imputed
SNPs and (iii) can provide more accurate estimates of effect sizes in regions
of association.


  We review the problem of confounding in genetic association studies, which
arises principally because of population structure and cryptic relatedness.
Many treatments of the problem consider only a simple ``island'' model of
population structure. We take a broader approach, which views population
structure and cryptic relatedness as different aspects of a single confounder:
the unobserved pedigree defining the (often distant) relationships among the
study subjects. Kinship is therefore a central concept, and we review methods
of defining and estimating kinship coefficients, both pedigree-based and
marker-based. In this unified framework we review solutions to the problem of
population structure, including family-based study designs, genomic control,
structured association, regression control, principal components adjustment and
linear mixed models. The last solution makes the most explicit use of the
kinships among the study subjects, and has an established role in the analysis
of animal and plant breeding studies. Recent computational developments mean
that analyses of human genetic association data are beginning to benefit from
its powerful tests for association, which protect against population structure
and cryptic kinship, as well as intermediate levels of confounding by the
pedigree.


  Genome-wide association studies, in which as many as a million single
nucleotide polymorphisms (SNP) are measured on several thousand samples, are
quickly becoming a common type of study for identifying genetic factors
associated with many phenotypes. There is a strong assumption that interactions
between SNPs or genes and interactions between genes and environmental factors
substantially contribute to the genetic risk of a disease. Identification of
such interactions could potentially lead to increased understanding about
disease mechanisms; drug $\times$ gene interactions could have profound
applications for personalized medicine; strong interaction effects could be
beneficial for risk prediction models. In this paper we provide an overview of
different approaches to model interactions, emphasizing approaches that make
specific use of the structure of genetic data, and those that make specific
modeling assumptions that may (or may not) be reasonable to make. We conclude
that to identify interactions it is often necessary to do some selection of
SNPs, for example, based on prior hypothesis or marginal significance, but that
to identify SNPs that are marginally associated with a disease it may also be
useful to consider larger numbers of interactions.


  In genome-wide association studies (GWAS), hundreds of thousands of genetic
markers (SNPs) are tested for association with a trait or phenotype. Reported
effects tend to be larger in magnitude than the true effects of these markers,
the so-called ``winner's curse.'' We argue that the classical definition of
unbiasedness is not useful in this context and propose to use a different
definition of unbiasedness that is a property of the estimator we advocate. We
suggest an integrated approach to the estimation of the SNP effects and to the
prediction of trait values, treating SNP effects as random instead of fixed
effects. Statistical methods traditionally used in the prediction of trait
values in the genetics of livestock, which predates the availability of SNP
data, can be applied to analysis of GWAS, giving better estimates of the SNP
effects and predictions of phenotypic and genetic values in individuals.


  Nested sampling is a Bayesian sampling technique developed to explore
probability distributions lo- calised in an exponentially small area of the
parameter space. The algorithm provides both posterior samples and an estimate
of the evidence (marginal likelihood) of the model. The nested sampling algo-
rithm also provides an efficient way to calculate free energies and the
expectation value of thermodynamic observables at any temperature, through a
simple post-processing of the output. Previous applications of the algorithm
have yielded large efficiency gains over other sampling techniques, including
parallel tempering (replica exchange). In this paper we describe a parallel
implementation of the nested sampling algorithm and its application to the
problem of protein folding in a Go-type force field of empirical potentials
that were designed to stabilize secondary structure elements in
room-temperature simulations. We demonstrate the method by conducting folding
simulations on a number of small proteins which are commonly used for testing
protein folding procedures: protein G, the SH3 domain of Src tyrosine kinase
and chymotrypsin inhibitor 2. A topological analysis of the posterior samples
is performed to produce energy landscape charts, which give a high level
description of the potential energy surface for the protein folding
simulations. These charts provide qualitative insights into both the folding
process and the nature of the model and force field used.


  Plant-pollinator mutualistic networks are asymmetric in their interactions:
specialist plants are pollinated by generalist animals, while generalist plants
are pollinated by a broad involving specialists and generalists. It has been
suggested that this asymmetric ---or disassortative--- assemblage could play an
important role in determining the equal susceptibility of specialist and
generalist plants under habitat destruction. At the core of the argument lies
the observation that specialist plants, otherwise candidates to extinction,
could cope with the disruption thanks to their interaction with generalist
pollinators. We present a theoretical framework that supports this thesis. We
analyze a dynamical model of a system of mutualistic plants and pollinators,
subject to the destruction of their habitat. We analyze and compare two
families of interaction topologies, ranging from highly assortative to highly
disassortative ones, as well as real pollination networks. We found that
several features observed in natural systems are predicted by the mathematical
model. First, there is a tendency to increase the asymmetry of the network as a
result of the extinctions. Second, an entropy measure of the differential
susceptibility to extinction of specialist and generalist species show that
they tend to balance when the network is disassortative. Finally, the
disappearance of links in the network, as a result of extinctions, shows that
specialist plants preserve more connections than the corresponding plants in an
assortative system, enabling them to resist the disruption.


  We show that autocatalytic networks of epsilon-machines and their population
dynamics differ substantially between spatial (geographically distributed) and
nonspatial (panmixia) populations. Generally, regions of spacetime-invariant
autocatalytic networks---or domains---emerge in geographically distributed
populations. These are separated by functional membranes of complementary
epsilon-machines that actively translate between the domains and are
responsible for their growth and stability. We analyze both spatial and
nonspatial populations, determining the algebraic properties of the
autocatalytic networks that allow for space to affect the dynamics and so
generate autocatalytic domains and membranes. In addition, we analyze
populations of intermediate spatial architecture, delineating the thresholds at
which spatial memory (information storage) begins to determine the character of
the emergent auto-catalytic organization.


  Copy number variants (CNVs) account for more polymorphic base pairs in the
human genome than do single nucleotide polymorphisms (SNPs). CNVs encompass
genes as well as noncoding DNA, making these polymorphisms good candidates for
functional variation. Consequently, most modern genome-wide association studies
test CNVs along with SNPs, after inferring copy number status from the data
generated by high-throughput genotyping platforms. Here we give an overview of
CNV genomics in humans, highlighting patterns that inform methods for
identifying CNVs. We describe how genotyping signals are used to identify CNVs
and provide an overview of existing statistical models and methods used to
infer location and carrier status from such data, especially the most commonly
used methods exploring hybridization intensity. We compare the power of such
methods with the alternative method of using tag SNPs to identify CNV carriers.
As such methods are only powerful when applied to common CNVs, we describe two
alternative approaches that can be informative for identifying rare CNVs
contributing to disease risk. We focus particularly on methods identifying de
novo CNVs and show that such methods can be more powerful than case-control
designs. Finally we present some recommendations for identifying CNVs
contributing to common complex disorders.


  The notion of information pervades informal descriptions of biological
systems, but formal treatments face the problem of defining a quantitative
measure of information rooted in a concept of fitness, which is itself an
elusive notion. Here, we present a model of population dynamics where this
problem is amenable to a mathematical analysis. In the limit where any
information about future environmental variations is common to the members of
the population, our model is equivalent to known models of financial
investment. In this case, the population can be interpreted as a portfolio of
financial assets and previous analyses have shown that a key quantity of
Shannon's communication theory, the mutual information, sets a fundamental
limit on the value of information. We show that this bound can be violated when
accounting for features that are irrelevant in finance but inherent to
biological systems, such as the stochasticity present at the individual level.
This leads us to generalize the measures of uncertainty and information usually
encountered in information theory.


  Objections to Darwinian evolution are often based on the time required to
carry out the necessary mutations. Seemingly, exponential numbers of mutations
are needed. We show that such estimates ignore the effects of natural
selection, and that the numbers of necessary mutations are thereby reduced to
about $K\log{L}$, rather than $K^L$, where $L$ is the length of the genomic
"word", and $K$ is the number of possible "letters" that can occupy any
position in the word. The required theory makes contact with the theory of
radix-exchange sorting in theoretical computer science, and the asymptotic
analysis of certain sums that occur there.


  Background and Aims: Prediction of phenotypic traits from new genotypes under
untested environmental conditions is crucial to build simulations of breeding
strategies to improve target traits. Although the plant response to
environmental stresses is characterized by both architectural and functional
plasticity, recent attempts to integrate biological knowledge into genetics
models have mainly concerned specific physiological processes or crop models
without architecture, and thus may prove limited when studying genotype x
environment interactions. Consequently, this paper presents a simulation study
introducing genetics into a functional-structural growth model, which gives
access to more fundamental traits for quantitative trait loci (QTL) detection
and thus to promising tools for yield optimization. Methods: The GreenLab model
was selected as a reasonable choice to link growth model parameters to QTL.
Virtual genes and virtual chromosomes were defined to build a simple genetic
model that drove the settings of the species-specific parameters of the model.
The QTL Cartographer software was used to study QTL detection of simulated
plant traits. A genetic algorithm was implemented to define the ideotype for
yield maximization based on the model parameters and the associated allelic
combination. Key Results and Conclusions: By keeping the environmental factors
constant and using a virtual population with a large number of individuals
generated by a Mendelian genetic model, results for an ideal case could be
simulated. Virtual QTL detection was compared in the case of phenotypic traits
- such as cob weight - and when traits were model parameters, and was found to
be more accurate in the latter case. The practical interest of this approach is
illustrated by calculating the parameters (and the corresponding genotype)
associated with yield optimization of a GreenLab maize model. The paper
discusses the potentials of GreenLab to represent environment x genotype
interactions, in particular through its main state variable, the ratio of
biomass supply over demand.


  Studies of cooperation have traditionally focused on discrete games such as
the well-known prisoner's dilemma, in which players choose between two pure
strategies: cooperation and defection. Increasingly, however, cooperation is
being studied in continuous games that feature a continuum of strategies
determining the level of cooperative investment. For the continuous snowdrift
game, it has been shown that a gradually evolving monomorphic population may
undergo evolutionary branching, resulting in the emergence of a defector
strategy that coexists with a cooperator strategy. This phenomenon has been
dubbed the 'tragedy of the commune'. Here we study the effects of fluctuating
group size on the tragedy of the commune and derive analytical conditions for
evolutionary branching. Our results show that the effects of fluctuating group
size on evolutionary dynamics critically depend on the structure of payoff
functions. For games with additively separable benefits and costs, fluctuations
in group size make evolutionary branching less likely, and sufficiently large
fluctuations in group size can always turn an evolutionary branching point into
a locally evolutionarily stable strategy. For games with multiplicatively
separable benefits and costs, fluctuations in group size can either prevent or
induce the tragedy of the commune. For games with general interactions between
benefits and costs, we derive a general classification scheme based on second
derivatives of the payoff function, to elucidate when fluctuations in group
size help or hinder cooperation.


  The Wnt signalling pathway transducing the stabilization of beta-catenin is
essential for metazoan embryo development and is misregulated in many diseases
such as cancers. In recent years models have been proposed for the Wnt
signalling pathway during the segmentation process in developing embryos. Many
of these include negative feedback loops build around Axin2 regulation. In this
article we propose a new negative feedback model for the Wnt pathway with
Dickkopf1 (Dkk1) at its core. Dkk1 is a negative regulator of Wnt signalling.
In chicken and mouse embryos there is a gradient of Wnt in the presomitic
mesoderm (PSM) decreasing from the posterior to the anterior end. The formation
of somites and the oscillations of Wnt target genes are controlled by this
gradient. Here we incorporate a Wnt gradient and show that synchronization of
neighbouring cells in the PSM is important in accordance with experimental
observations.


  A random null model termed the Blind Watchmaker network (BW) has been shown
to reproduce the degree distribution found in metabolic networks. This might
suggest that natural selection has had little influence on this particular
network property. We here investigate to what extent other structural network
properties have evolved under selective pressure from the corresponding ones of
the random null model: The clustering coefficient and the assortativity
measures are chosen and it is found that these measures for the metabolic
network structure are close enough to the BW-network so as to fit inside its
reachable random phase space. It is furthermore shown that the use of this null
model indicates an evolutionary pressure towards low assortativity and that
this pressure is stronger for larger networks. It is also shown that selecting
for BW networks with low assortativity causes the BW degree distribution to
slightly deviate from its power-law shape in the same way as the metabolic
networks. This implies that an equilibrium model with fluctuating degree
distribution is more suitable as a null model, when identifying selective
pressures, than a randomized counterpart with fixed degree sequence, since the
overall degree sequence itself can change under selective pressure on other
global network properties.


  A major challenge in the analysis of population genomics data consists of
isolating signatures of natural selection from background noise caused by
random drift and gene flow. Analyses of massive amounts of data from many
related populations require high-performance algorithms to determine the
likelihood of different demographic scenarios that could have shaped the
observed neutral single nucleotide polymorphism (SNP) allele frequency
spectrum. In many areas of applied mathematics, Fourier Transforms and Spectral
Methods are firmly established tools to analyze spectra of signals and model
their dynamics as solutions of certain Partial Differential Equations (PDEs).
When spectral methods are applicable, they have excellent error properties and
are the fastest possible in high dimension. In this paper we present an
explicit numerical solution, using spectral methods, to the forward Kolmogorov
equations for a Wright-Fisher process with migration of K populations, influx
of mutations, and multiple population splitting events.


  In the first part of this note, we show (following Hofbauer and Sigmund) that
Dawkins' "Battle of the Sexes" defines an interesting map from a space of
economic parameters to psychosocial coordinates. The second part discusses an
even more elementary game, but one which is not completely trivial.


  An important question in genome evolution is whether there exist fragile
regions (rearrangement hotspots) where chromosomal rearrangements are happening
over and over again. Although nearly all recent studies supported the existence
of fragile regions in mammalian genomes, the most comprehensive phylogenomic
study of mammals (Ma et al. (2006) Genome Research 16, 1557-1565) raised some
doubts about their existence. We demonstrate that fragile regions are subject
to a "birth and death" process, implying that fragility has limited
evolutionary lifespan. This finding implies that fragile regions migrate to
different locations in different mammals, explaining why there exist only a few
chromosomal breakpoints shared between different lineages. The birth and death
of fragile regions phenomenon reinforces the hypothesis that rearrangements are
promoted by matching segmental duplications and suggests putative locations of
the currently active fragile regions in the human genome.


  To reduce global biodiversity loss, there is an urgent need to determine the
most efficient allocation of conservation resources. Recently, there has been a
growing trend for many governments to supplement public ownership and
management of reserves with incentive programs for conservation on private
land. At the same time, policies to promote conservation on private land are
rarely evaluated in terms of their ecological consequences. This raises
important questions, such as the extent to which private land conservation can
improve conservation outcomes, and how it should be mixed with more traditional
public land conservation. We address these questions, using a general framework
for modelling environmental policies and a case study examining the
conservation of endangered native grasslands to the west of Melbourne,
Australia. Specifically, we examine three policies that involve: i) spending
all resources on creating public conservation areas; ii) spending all resources
on an ongoing incentive program where private landholders are paid to manage
vegetation on their property with 5-year contracts; and iii) splitting
resources between these two approaches. The performance of each strategy is
quantified with a vegetation condition change model that predicts future
changes in grassland quality. Of the policies tested, no one policy was always
best and policy performance depended on the objectives of those enacting the
policy. This work demonstrates a general method for evaluating environmental
policies and highlights the utility of a model which combines ecological and
socioeconomic processes.


  The mean fixation time of a deleterious mutant allele is studied beyond the
diffusion approximation. As in Kimura's classical work [M. Kimura, Proc. Natl.
Acad. Sci. U.S.A. Vol.77, 522 (1980)], that was motivated by the problem of
fixation in the presence of amorphic or hypermorphic mutations, we consider a
diallelic model at a single locus comprising a wild-type A and a mutant allele
A' produced irreversibly from A at small uniform rate v. The relative fitnesses
of the mutant homozygotes A'A', mutant heterozygotes A'A and wild-type
homozygotes AA are 1-s, 1-h and 1, respectively, where it is assumed that v<<
s. Here, we adopt an approach based on the direct treatment of the underlying
Markov chain (birth-death process) obeyed by the allele frequency (whose
dynamics is prescribed by the Moran model), which allows to accurately account
for the effects of large fluctuations. After a general description of the
theory, we focus on the case of a deleterious mutant allele (i.e. s>0) and
discuss three situations: when the mutant is (i) completely dominant (s=h);
(ii) completely recessive (h=0), and (iii) semi-dominant (h=s/2). Our
theoretical predictions for the mean fixation time and the quasi-stationary
distribution of the mutant population in the coexistence state, are shown to be
in excellent agreement with numerical simulations. Furthermore, when s is
finite, we demonstrate that our results are superior to those of the diffusion
theory that is shown to be an accurate approximation only when N_e s^2 << 1,
where N_e is the effective population size.


  We present a mean field model for the phase diagram of a community of
micro-organisms, interacting through their metabolism so that they are, in
effect, engaging in a cooperative social game. We show that as a function of
the concentration of the nutrients glucose and histidine, the community
undergoes a phase transition separating a state in which one strain is dominant
to a state which is characterized by coexisting populations. Our results are in
good agreement with recent experimental results, correctly predicting
quantitative trends and the phase diagram.


  Numerous experiments demonstrate a high level of promiscuity and structural
disorder in organismal proteomes. Here we ask the question what makes a protein
promiscuous, i.e., prone to non-specific interactions, and structurally
disordered. We predict that multi-scale correlations of amino acid positions
within protein sequences statistically enhance the propensity for promiscuous
intra- and inter-protein binding. We show that sequence correlations between
amino acids of the same type are statistically enhanced in structurally
disordered proteins and in hubs of organismal proteomes. We also show that
structurally disordered proteins possess a significantly higher degree of
sequence order than structurally ordered proteins. We develop an analytical
theory for this effect and predict the robustness of our conclusions with
respect to the amino acid composition and the form of the microscopic potential
between the interacting sequences. Our findings have implications for
understanding molecular mechanisms of protein aggregation diseases induced by
the extension of sequence repeats.


  Background: Speciation corresponds to the progressive establishment of
reproductive barriers between groups of individuals derived from an ancestral
stock. Since Darwin did not believe that reproductive barriers could be
selected for, he proposed that most events of speciation would occur through a
process of separation and divergence, and this point of view is still shared by
most evolutionary biologists today. Results: I do, however, contend that, if so
much speciation occurs, the most likely explanation is that there must be
conditions where reproductive barriers can be directly selected for. In other
words, situations where it is advantageous for individuals to reproduce
preferentially within a small group and reduce their breeding with the rest of
the ancestral population. This leads me to propose a model whereby new species
arise not by populations splitting into separate branches, but by small
inbreeding groups "budding" from an ancestral stock. This would be driven by
several advantages of inbreeding, and mainly by advantageous recessive
phenotypes, which could only be retained in the context of inbreeding.
Reproductive barriers would thus not arise as secondary consequences of
divergent evolution in populations isolated from one another, but under the
direct selective pressure of ancestral stocks. Many documented cases of
speciation in natural populations appear to fit the model proposed, with more
speciation occurring in populations with high inbreeding coefficients, and many
recessive characters identified as central to the phenomenon of speciation,
with these recessive mutations expected to be surrounded by patterns of limited
genomic diversity. Conclusions: Whilst adaptive evolution would correspond to
gains of function that would, most of the time, be dominant, this type of
speciation by budding would thus be driven by mutations resulting in the
advantageous loss of certain functions since recessive mutations very often
correspond to the inactivation of a gene. A very important further advantage of
inbreeding is that it reduces the accumulation of recessive mutations in
genomes. A consequence of the model proposed is that the existence of species
would correspond to a metastable equilibrium between inbreeding and
outbreeding, with excessive inbreeding promoting speciation, and excessive
outbreeding resulting in irreversible accumulation of recessive mutations that
could ultimately only lead to extinction.


  In a genetic algorithm, fluctuations of the entropy of a genome over time are
interpreted as fluctuations of the information that the genome's organism is
storing about its environment, being this reflected in more complex organisms.
The computation of this entropy presents technical problems due to the small
population sizes used in practice. In this work we propose and test an
alternative way of measuring the entropy variation in a population by means of
algorithmic information theory, where the entropy variation between two
generational steps is the Kolmogorov complexity of the first step conditioned
to the second one. As an example application of this technique, we report
experimental differences in entropy evolution between systems in which sexual
reproduction is present or absent.


  Our micromechanical experiments show that at physiological temperatures type
I collagen fibril has several basic features of the glassy state. The
transition out of this state [softening transition] essentially depends on the
speed of heating v, e.g., for v=1 C/min it occurs around 70 C and is displayed
by a peak of the internal friction and decreasing Young's modulus. The
softening transition decreases by 45 C upon decreasing the heating speed to
v=0.1 C/min. For temperatures 20-30 C the native collagen fibril demonstrates
features of mechanical glassines at oscillation frequencies 0.1-3 kHz; in
particular, the internal friction has a sharp maximum as a function of the
frequency. This is the first example of biopolymer glassines at physiological
temperatures, because well-known glassy features of DNA and globular proteins
are seen only for much lower temperatures (around 200 K).


  Metapopulation theory for a long time has assumed dispersal to be symmetric,
i.e. patches are connected through migrants dispersing bi-directionally without
a preferred direction. However, for natural populations symmetry is often
broken, e.g. for species in the marine environment dispersing through the
transport of pelagic larvae with ocean currents. The few recent studies of
asymmetric dispersal concluded, that asymmetry has a distinct negative impact
on the persistence of metapopulations. Detailed analysis however revealed, that
these previous studies might have been unable to properly disentangle the
effect of symmetry from other potentially confounding properties of dispersal
patterns. We resolve this issue by systematically investigating the symmetry of
dispersal patterns and its impact on metapopulation persistence. Our main
analysis based on a metapopulation model equivalent to previous studies but now
applied on regular dispersal patterns aims to isolate the effect of dispersal
symmetry on metapopulation persistence. Our results suggest, that asymmetry in
itself does not imply negative effects on metapopulation persistence. For this
reason we recommend to investigate it in connection with other properties of
dispersal instead of in isolation.


  Summary: Macromolecular assembly vertebrates essential cellular processes,
such as gene regulation and signal transduction. A major challenge for
conventional computational methods to study these processes is tackling the
exponential increase of the number of configurational states with the number of
components. CplexA is a Mathematica package that uses functional programming to
efficiently compute probabilities and average properties over such
exponentially large number of states from the energetics of the interactions.
The package is particularly suited to study gene expression at complex
promoters controlled by multiple, local and distal, DNA binding sites for
transcription factors. Availability: CplexA is freely available together with
documentation at http://sourceforge.net/projects/cplexa/.


  Focusing on the DNA-bridging nucleoid proteins Fis and H-NS, and integrating
several independent experimental and bioinformatic data sources, we investigate
the links between chromosomal spatial organization and global transcriptional
regulation. By means of a novel multi-scale spatial aggregation analysis, we
uncover the existence of contiguous clusters of nucleoid-perturbation sensitive
genes along the genome, whose expression is affected by a combination of
topological DNA state and nucleoid-shaping protein occupancy. The clusters
correlate well with the macrodomain structure of the genome. The most
significant of them lay symmetrically at the edges of the ter macrodomain and
involve all of the flagellar and chemotaxis machinery, in addition to key
regulators of biofilm formation, suggesting that the regulation of the physical
state of the chromosome by the nucleoid proteins plays an important role in
coordinating the transcriptional response leading to the switch between a
motile and a biofilm lifestyle.


  Using lattice models we explore the factors that determine the tendencies of
polypeptide chains to aggregate by exhaustively sampling the sequence and
conformational space. The morphologies of the fibril-like structures and the
time scales ($\tau_{fib}$) for their formation depend on a balance between
hydrophobic and coulomb interactions. The extent of population of an ensemble
of \textbf{N$^*$} structures, which are fibril-prone structures in the spectrum
of conformations of an isolated protein, is the major determinant of
$\tau_{fib}$. This observation is used to determine the aggregating sequences
by exhaustively exploring the sequence space, thus providing a basis for genome
wide search of fragments that are aggregation prone.


  One of the main necessities for population geneticists is the availability of
statistical tools that enable to accept or reject the neutral Wright-Fisher
model with high power. A number of statistical tests have been developed to
detect specific deviations from the null frequency spectrum in different
directions (i.e., Tajima's D, Fu and Li's F and D test, Fay and Wu's H).
Recently, a general framework was proposed to generate all neutrality tests
that are linear functions of the frequency spectrum. In this framework, a
family of optimal tests was developed to have almost maximum power against a
specific alternative evolutionary scenario. Following these developments, in
this paper we provide a thorough discussion of linear and nonlinear neutrality
tests. First, we present the general framework for linear tests and emphasize
the importance of the property of scalability with the sample size (that is,
the results of the tests should not depend on the sample size), which, if
missing, can guide to errors in data interpretation. The motivation and
structure of linear optimal tests are discussed. In a further generalization,
we develop a general framework for nonlinear neutrality tests and we derive
nonlinear optimal tests for polynomials of any degree in the frequency
spectrum.


  We study the evolution of sexual and asexual populations in general fitness
landscapes. We find deep relations between the mathematics of biological
evolution and the formalism of quantum mechanics. We give the general structure
of the evolution of populations which is in general an off-equilibrium process
that can be expressed by path integrals over phylogenies. These phylogenies are
sums of linear lineages for asexual populations. For sexual populations
instead, each lineage is a tree of branching ratio two and the path integral
describing the evolving population is given by a sum over these trees. Finally,
we show that the Bose-Einstein and the Fermi-Dirac distributions describe the
stationary state of biological populations in simple cases.


  We derive a Poisson random field model for population site polymorphisms
differences within and between two species that share a relatively recent
common ancestor. The model can be either equilibrium or time inhomogeneous. We
first consider a random field of Markov chains that describes the fate of a set
of individual mutations. This field is approximated by a Poisson random field
from which we can make inferences about the amounts of mutation and selection
that have occurred in the history of observed aligned DNA sequences.


  Cooperation and self-organized criticality are two main keywords in current
studies of evolution. We propose a generalized Bak-Sneppen model and provide a
natural mechanism which accounts for both phenomena simultaneously. We use the
prisoner's dilemma games to mimic the interactions among the members of the
population. Each member is identified by its cooperation probability, and its
fitness is given by the payoffs from neighbors. The least fit member with the
minimum payoff is replaced by a new member with a random cooperation
probability. When the neighbors of the least fit member are also replaced with
a non-zero probability, a strong cooperation emerges. The Bak-Sneppen process
builds a self-organized structure so that the cooperation can emerge even in
the parameter region where a uniform or random population decreases the number
of cooperators. The emergence of cooperation is due to the same dynamical
correlation which leads to self-organized criticality in replacement
activities.


  We develop Bayesian inference methods for a recently-emerging type of
epigenetic data to study the transmission fidelity of DNA methylation patterns
over cell divisions. The data consist of parent-daughter double-stranded DNA
methylation patterns with each pattern coming from a single cell and
represented as an unordered pair of binary strings. The data are technically
difficult and time-consuming to collect, putting a premium on an efficient
inference method. Our aim is to estimate rates for the maintenance and de novo
methylation events that gave rise to the observed patterns, while accounting
for measurement error. We model data at multiple sites jointly, thus using
whole-strand information, and considerably reduce confounding between
parameters. We also adopt a hierarchical structure that allows for variation in
rates across sites without an explosion in the effective number of parameters.
Our context-specific priors capture the expected stationarity, or
near-stationarity, of the stochastic process that generated the data analyzed
here. This expected stationarity is shown to greatly increase the precision of
the estimation. Applying our model to a data set collected at the human FMR1
locus, we find that measurement errors, generally ignored in similar studies,
occur at a nontrivial rate (inappropriate bisulfite conversion error: 1.6$%$
with 80$%$ CI: 0.9--2.3$%$). Accounting for these errors has a substantial
impact on estimates of key biological parameters. The estimated average failure
of maintenance rate and daughter de novo rate decline from 0.04 to 0.024 and
from 0.14 to 0.07, respectively, when errors are accounted for. Our results
also provide evidence that de novo events may occur on both parent and daughter
strands: the median parent and daughter de novo rates are 0.08 (80$%$ CI:
0.04--0.13) and 0.07 (80$%$ CI: 0.04--0.11), respectively.


  The nucleotide composition of human genes with a special emphasis on
transcription-related strand asymmetries is analyzed. Such asymmetries may be
associated with different mutational rates in two principal factors. The first
one is transcription-coupled repair and the second one is the selective
pressure related to optimization of the translation efficiency. The former
factor affects both coding and noncoding regions of a gene, while the latter
factor is applicable only to the coding regions. Compositional asymmetries
calculated at the third position of a codon in coding (exons) and noncoding
(introns, UTR, upstream and downstream) regions of human genes are compared. It
is shown that the keto-skew (excess of the frequencies of G and T nucleotides
over the frequencies of A and C nucleotides in the same strand) is most
pronounced in intronic regions, less pronounced in coding regions, and has near
zero values in untranscribed regions. The keto-skew correlates with the level
of gene expression in germ-line cells in both introns and exons. We propose to
use the results of our analysis to estimate the contribution of different
evolutionary factors to the transcription-related compositional biases.


  Eigen's quasi-species model describes viruses as ensembles of different
mutants of a high fitness "master" genotype. Mutants are assumed to have lower
fitness than the master type, yet they coexist with it forming the
quasi-species. When the mutation rate is sufficiently high, the master type no
longer survives and gets replaced by a wide range of mutant types, thus
destroying the quasi-species. It is the so-called "error catastrophe". But
natural selection acts on phenotypes, not genotypes, and huge amounts of
genotypes yield the same phenotype. An important consequence of this is the
appearance of beneficial mutations which increase the fitness of mutants. A
model has been recently proposed to describe quasi-species in the presence of
beneficial mutations. This model lacks the error catastrophe of Eigen's model
and predicts a steady state in which the viral population grows exponentially.
Extinction can only occur if the infectivity of the quasi-species is so low
that this exponential is negative. In this work I investigate the transient of
this model when infection is started from a small amount of low fitness
virions. I prove that, beyond an initial regime where viral population
decreases (and can go extinct), the growth of the population is
super-exponential. Hence this population quickly becomes so huge that selection
due to lack of host cells to be infected begins to act before the steady state
is reached. This result suggests that viral infection may widespread before the
virus has developed its optimal form.


  Mass spectrometry provides a high-throughput way to identify proteins in
biological samples. In a typical experiment, proteins in a sample are first
broken into their constituent peptides. The resulting mixture of peptides is
then subjected to mass spectrometry, which generates thousands of spectra, each
characteristic of its generating peptide. Here we consider the problem of
inferring, from these spectra, which proteins and peptides are present in the
sample. We develop a statistical approach to the problem, based on a nested
mixture model. In contrast to commonly used two-stage approaches, this model
provides a one-stage solution that simultaneously identifies which proteins are
present, and which peptides are correctly identified. In this way our model
incorporates the evidence feedback between proteins and their constituent
peptides. Using simulated data and a yeast data set, we compare and contrast
our method with existing widely used approaches (PeptideProphet/ProteinProphet)
and with a recently published new approach, HSM. For peptide identification,
our single-stage approach yields consistently more accurate results. For
protein identification the methods have similar accuracy in most settings,
although we exhibit some scenarios in which the existing methods perform
poorly.


  The principles underlying protein folding remains one of Nature's puzzles
with important practical consequences for Life. An approach that has gathered
momentum since the late 1990's, looks at protein hetero-polymers and their
folding process through the lens of complex network analysis. Consequently,
there is now a body of empirical studies describing topological characteristics
of protein macro-molecules through their contact networks and linking these
topological characteristics to protein folding. The present paper is primarily
a review of this rich area. But it delves deeper into certain aspects by
emphasizing short-range and long-range links, and suggests unconventional
places where "power-laws" may be lurking within protein contact networks.
Further, it considers the dynamical view of protein contact networks. This
closer scrutiny of protein contact networks raises new questions for further
research, and identifies new regularities which may be useful to parameterize a
network approach to protein folding. Preliminary experiments with such a model
confirm that the regularities we identified cannot be easily reproduced through
random effects. Indeed, the grand challenge of protein folding is to elucidate
the process(es) which not only generates the specific and diverse linkage
patterns of protein contact networks, but also reproduces the dynamic behavior
of proteins as they fold. Keywords: network analysis, protein contact networks,
protein folding


  We present a network model in which words over a specific alphabet, called
{\it structures}, are associated to each node and undirected edges are added
depending on some distance between different structures. It is shown that this
model can generate, without the use of preferential attachment or any other
heuristic, networks with topological features similar to biological networks:
power law degree distribution, clustering coefficient independent from the
network size, etc. Specific biological networks ({\it C. Elegans} neural
network and {\it E. Coli} protein-protein interaction network) are replicated
using this model.


  We study the dynamics of the voter and Moran processes running on top of
complex network substrates where each edge has a weight depending on the degree
of the nodes it connects. For each elementary dynamical step the first node is
chosen at random and the second is selected with probability proportional to
the weight of the connecting edge. We present a heterogeneous mean-field
approach allowing to identify conservation laws and to calculate exit
probabilities along with consensus times. In the specific case when the weight
is given by the product of nodes' degree raised to a power theta, we derive a
rich phase-diagram, with the consensus time exhibiting various scaling laws
depending on theta and on the exponent of the degree distribution gamma.
Numerical simulations give very good agreement for small values of |theta|. An
additional analytical treatment (heterogeneous pair approximation) improves the
agreement with numerics, but the theoretical understanding of the behavior in
the limit of large |theta| remains an open challenge.


  We consider nano-sized artificial or biological machines working in steady
state enforced by imposing non-equilibrium concentrations of solutes or by
applying external forces, torques or electric fields. For unicyclic and
strongly coupled multicyclic machines, efficiency at maximum power is not
bounded by the linear response value 1/2. For strong driving, it can even
approach the thermodynamic limit 1. Quite generally, such machines fall in
three different classes characterized, respectively, as "strong and efficient",
"strong and inefficient", and "balanced". For weakly coupled multicyclic
machines, efficiency at maximum power has lost any universality even in the
linear response regime.


  Aging is a fundamental aspect of living systems that undergo a progressive
deterioration of physiological function with age and an increase of
vulnerability to disease and death. Living systems, known as complex systems,
require complexity in interactions among molecules, cells, organs, and
individuals or regulatory mechanisms to perform a variety of activities for
survival. On this basis, aging can be understood in terms of a progressive loss
of complexity with age; this suggests that complexity in living systems would
evolve with age. In general, aging dynamics is mathematically depicted by a
survival function, which monotonically changes from 1 to 0 with age. It would
be then useful to find an adequate survival function to link aging dynamics and
complexity evolution. Here we describe a flexible survival function, which is
derived from the stretched exponential function by adopting an age-dependent
exponent. We note that the exponent is associated with evolving complexity,
i.e., a fractal-like scaling in cumulative mortality. The survival function
well depicts a general feature in survival curves; healthy populations show a
tendency to evolve towards rectangular-like survival curves, as examples in
humans or laboratory animals. This tendency suggests that both aging and
complexity would evolve towards healthy survival in living systems. Our
function to link aging with complexity may contribute to better understanding
of biological aging in terms of complexity evolution.