  A central theme in phylogenetics is the reconstruction and analysis of
evolutionary trees from a given set of data. To determine the optimal search
methods for reconstructing trees, it is crucial to understand the size and
structure of the neighbourhoods of trees under tree rearrangement operations.
The diameter and size of the immediate neighbourhood of a tree has been
well-studied, however little is known about the number of trees at distance
two, three or (more generally) $k$ from a given tree. In this paper we provide
a number of exact and asymptotic results concerning these quantities, and
identify some key aspects of tree shape that play a role in determining these
quantities. We obtain several new results for two of the main tree
rearrangement operations - Nearest Neighbour Interchange and Subtree Prune and
Regraft -- as well as for the Robinson-Foulds metric on trees.


  Given an RNA sequence a, consider the network G = (V;E), where the set V of
nodes consists of all secondary structures of a, and whose edge set E consists
of all edges connecting two secondary structures whose base pair distance is 1.
Define the network connectivity, or expected network degree, as the average
number of edges incident to vertices of G. Using algebraic combinatorial
methods, we prove that the asymptotic connectivity of length n homopolymer
sequences is 0:473418 ? n. This raises the question of what other network
properties are characteristic of the network of RNA secondary structures.
Programs in Python, C and Mathematica are available at the web site
http://bioinformatics.bc.edu/clotelab/ RNAexpNumNbors.


  Ant raiding, the process of identifying and returning food to the nest or
bivouac, is a fascinating example of collective motion in nature. During such
raids ants lay pheromones to form trails for others to find a food source. In
this work a coupled PDE/ODE model is introduced to study ant dynamics and
pheromone concentration. The key idea is the introduction of two forms of ant
dynamics: foraging and returning, each governed by different environmental and
social cues. The model accounts for all aspects of the raiding cycle including
local collisional interactions, the laying of pheromone along a trail, and the
transition from one class of ants to another. Through analysis of an order
parameter measuring the orientational order in the system, the model shows
self-organization into a collective state consisting of lanes of ants moving in
opposite directions as well as the transition back to the individual state once
the food source is depleted matching prior experimental results. This indicates
that in the absence of direct communication ants naturally form an efficient
method for transporting food to the nest/bivouac. The model exhibits a
continuous kinetic phase transition in the order parameter as a function of
certain system parameters. The associated critical exponents are found,
shedding light on the behavior of the system near the transition.


  The temperature dependence of DNA flexibility is studied in the presence of
stretching and unzipping forces. Two classes of models are considered. In one
case the origin of elasticity is entropic due to the polymeric correlations,
and in the other the double-stranded DNA is taken to have an intrinsic rigidity
for bending. In both cases single strands are completely flexible. The change
in the elastic constant for the flexible case due to thermally generated
bubbles is obtained exactly. For the case of intrinsic rigidity, the elastic
constant is found to be proportional to the square root of the bubble number
fluctuation.


  Living systems are often described utilizing informational analogies. An
important open question is whether information is merely a useful conceptual
metaphor, or intrinsic to the operation of biological systems. To address this
question, we provide a rigorous case study of the informational architecture of
two representative biological networks: the Boolean network model for the
cell-cycle regulatory network of the fission yeast S. pombe and that of the
budding yeast S. cerevisiae. We compare our results for these biological
networks to the same analysis performed on ensembles of two different types of
random networks. We show that both biological networks share features in common
that are not shared by either ensemble. In particular, the biological networks
in our study, on average, process more information than the random networks.
They also exhibit a scaling relation in information transferred between nodes
that distinguishes them from either ensemble: even when compared to the
ensemble of random networks that shares important topological properties, such
as a scale-free structure. We show that the most biologically distinct regime
of this scaling relation is associated with the dynamics and function of the
biological networks. Information processing in biological networks is therefore
interpreted as an emergent property of topology (causal structure) and dynamics
(function). These results demonstrate quantitatively how the informational
architecture of biologically evolved networks can distinguish them from other
classes of network architecture that do not share the same informational
properties.


  Solid-state nanopores are single molecule sensors that measure changes in
ionic current as charged polymers such as DNA pass through. Here, we present
comprehensive experiments on the length, voltage and salt dependence of the
frequency of double-stranded DNA translocations through conical quartz
nanopores with mean opening diameter 15 nm. We observe an entropic barrier
limited, length dependent translocation frequency at 4M LiCl salt concentration
and a drift-dominated, length independent translocation frequency at 1M KCl
salt concentration. These observations are described by a unifying
convection-diffusion equation which includes the contribution of an entropic
barrier for polymer entry.


  Spreading phenomena are ubiquitous in nature and society. For example,
disease, rumor, and information spread over underlying social and information
networks. It is well known that there is no threshold for epidemic models on
scale-free networks; this suggests that disease can spread on such networks,
regardless of how low the contact rate may be. In this paper, I consider six
models with different contact and propagation mechanisms. Each model is
analyzed by degree-based mean-field theory. I show that the presence or absence
of an outbreak threshold depends on the contact and propagation mechanism.


  We consider a model of asexually reproducing individuals with random
mutations and selection. The rate of mutations is proportional to the
population size, $N$. The mutations may be either beneficial or deleterious. In
a paper by Yu, Etheridge and Cuthbertson (2009) it was conjectured that the
average rate at which the mean fitness increases in this model is $O(\log
N/(\log\log N)^2)$. In this paper we show that for any time $t > 0$ there exist
values $\epsilon_N \rightarrow 0$ and a fixed $c > 0$ such that the maximum
fitness of the population is greater than $cs\log N/(\log\log N)^2$ for all
times $s \in [\epsilon_N,t]$ with probability tending to 1 as $N$ tends to
infinity.


  We introduce an algorithm for the alignment of protein- coding sequences
accounting for frameshifts. The main specificity of this algorithm as compared
to previously published protein-coding sequence alignment methods is the
introduction of a penalty cost for frameshift ex- tensions. Previous algorithms
have only used constant frameshift penal- ties. This is similar to the use of
scoring schemes with affine gap penalties in classical sequence alignment
algorithms. However, the overall penalty of a frameshift portion in an
alignment cannot be formulated as an affine function, because it should also
incorporate varying codon substitution scores. The second specificity of the
algorithm is its search space being the set of all possible alignments between
two coding sequences, under the classical definition of an alignment between
two DNA sequences. Previous algorithms have introduced constraints on the
length of the alignments, and additional symbols for the representation of
frameshift openings in an alignment. The algorithm has the same asymptotic
space and time complexity as the classical Needleman-Wunsch algorithm.


  The distribution of human linguistic groups presents a number of interesting
and non-trivial patterns. The distributions of the number of speakers per
language and the area each group covers follow log-normal distributions, while
population and area fulfill an allometric relationship. The topology of
networks of spatial contacts between different linguistic groups has been
recently characterized, showing atypical properties of the degree distribution
and clustering, among others. Human demography, spatial conflicts, and the
construction of networks of contacts between linguistic groups are mutually
dependent processes. Here we introduce an adaptive network model that takes all
of them into account and successfully reproduces, using only four model
parameters, not only those features of linguistic groups already described in
the literature, but also correlations between demographic and topological
properties uncovered in this work. Besides their relevance when modeling and
understanding processes related to human biogeography, our adaptive network
model admits a number of generalizations that broaden its scope and make it
suitable to represent interactions between agents based on population dynamics
and competition for space.


  The general deterministic recombination equation in continuous time is
analysed for various lattices, with special emphasis on the lattice of interval
(or ordered) partitions. Based on the recently constructed general solution for
the lattice of all partitions, the corresponding solution for interval
partitions is derived and analysed in detail. We focus our attention on the
recursive structure of the solution and its decay rates, and also discuss the
solution in the degenerate cases, where it comprises products of monomials with
exponentially decaying factors. This can be understood via the Markov generator
of the underlying partitioning process that was recently identified. We use
interval partitions to gain insight into the structure of the solution, while
our general framework works for arbitrary lattices.


  The possible control of competitive invasion by infection of the invader and
multiplicative noise is studied. The basic model is the Lotka-Volterra
competition system with emergent carrying capacities. Several stationary
solutions of the non-infected and infected system are identified as well as
parameter ranges of bistability. The latter are used for the numerical study of
invasion phenomena. The diffusivities, the infection but in particular the
white and colored multiplicative noise are the control parameters. It is shown
that not only competition, possible infection and mobilities are important
drivers of the invasive dynamics but also the noise and especially its color
and the functional response of populations to the emergence of noise.


  Non structural protein 3 (NS3) helicase from hepatitis C virus is an enzyme
that unwinds and translocates along nucleic acids with an ATP-dependent
mechanism and has a key role in the replication of the viral RNA. An
inchworm-like mechanism for translocation has been proposed based on crystal
structures and single molecule experiments. We here perform atomistic molecular
dynamics in explicit solvent on the microsecond time scale of the available
experimental structures. We also construct and simulate putative intermediates
for the translocation process, and we perform non-equilibrium targeted
simulations to estimate their relative stability. For each of the simulated
structures we carefully characterize the available conformational space, the
ligand binding pocket, and the RNA binding cleft. The analysis of the hydrogen
bond network and of the non-equilibrium trajectories indicates an ATP-dependent
stabilization of one of the protein conformers. Additionally, enthalpy
calculations suggest that entropic effects might be crucial for the
stabilization of the experimentally observed structures.


  The stability and dynamics of a double-stranded DNA (dsDNA) is affected by
the preferential occupancy of small monovalent molecular ions. Small metal and
molecular ions such as sodium and alkyl ammonium have crucial biological
functions in human body, affect the thermodynamic stability of the duplex DNA
and exhibit preferential binding. Here, using atomistic molecular dynamics
simulations we investigate the preferential binding of metal ion such as Na+
and molecular ions such as tetramethyl ammonium (TMA+) and
2-hydroxy-N,N,N-trimethylethanaminium (CHO+) to double stranded DNA. The
thermodynamic driving force for a particular molecular ion- DNA interaction is
determined by decomposing the free energy of binding into its entropic and
enthalpic contributions. Our simulations show that each of these molecular ions
preferentially binds to the minor groove of the DNA and the extent of binding
is highest for CHO+. The ion binding processes are found to be entropically
favourable. In addition, the contribution of hydrophobic effects towards the
entropic stabilization (in case of TMA+) and the effect of hydrogen bonding
contributing to enthalpic stabilization (in case of CHO+) have also been
investigated.


  Motivation: Many problems in biomedicine and other areas of the life sciences
can be characterized as control problems, with the goal of finding strategies
to change a disease or otherwise undesirable state of a biological system into
another, more desirable, state through an intervention, such as a drug or other
therapeutic treatment. The identification of such strategies is typically based
on a mathematical model of the process to be altered through targeted control
inputs. This paper focuses on processes at the molecular level that determine
the state of an individual cell, involving signaling or gene regulation. The
mathematical model type considered is that of Boolean networks. The potential
control targets can be represented by a set of nodes and edges that can be
manipulated to produce a desired effect on the system. Experimentally, node
manipulation requires technology to completely repress or fully activate a
particular gene product while edge manipulations only require a drug that
inactivates the interaction between two gene products. Results: This paper
presents a method for the identification of potential intervention targets in
Boolean molecular network models using algebraic techniques. The approach
exploits an algebraic representation of Boolean networks to encode the control
candidates in the network wiring diagram as the solutions of a system of
polynomials equations, and then uses computational algebra techniques to find
such controllers. The control methods in this paper are validated through the
identification of combinatorial interventions in the signaling pathways of
previously reported control targets in two well studied systems, a p53-mdm2
network and a blood T cell lymphocyte granular leukemia survival signaling
network.


  The article presents an application of Hidden Markov Models (HMMs) for
pattern recognition on genome sequences. We apply HMM for identifying genes
encoding the Variant Surface Glycoprotein (VSG) in the genomes of Trypanosoma
brucei (T. brucei) and other African trypanosomes. These are parasitic protozoa
causative agents of sleeping sickness and several diseases in domestic and wild
animals. These parasites have a peculiar strategy to evade the host's immune
system that consists in periodically changing their predominant cellular
surface protein (VSG). The motivation for using patterns recognition methods to
identify these genes, instead of traditional homology based ones, is that the
levels of sequence identity (amino acid and DNA sequence) amongst these genes
is often below of what is considered reliable in these methods. Among pattern
recognition approaches, HMM are particularly suitable to tackle this problem
because they can handle more naturally the determination of gene edges. We
evaluate the performance of the model using different number of states in the
Markov model, as well as several performance metrics. The model is applied
using public genomic data. Our empirical results show that the VSG genes on T.
brucei can be safely identified (high sensitivity and low rate of false
positives) using HMM.


  Possibilities for using geometry and topology to analyze statistical problems
in biology raise a host of novel questions in geometry, probability, algebra,
and combinatorics that demonstrate the power of biology to influence the future
of pure mathematics. This expository article is a tour through some biological
explorations and their mathematical ramifications. The article starts with
evolution of novel topological features in wing veins of fruit flies, which are
quantified using the algebraic structure of multiparameter persistent homology.
The statistical issues involved highlight mathematical implications of sampling
from moduli spaces. These lead to geometric probability on stratified spaces,
including the sticky phenomenon for Frechet means and the origin of this
mathematical area in the reconstruction of phylogenetic trees.


  Advances in mathematical epidemiology have led to a better understanding of
the risks posed by epidemic spreading and informed strategies to contain
disease spread. However, a challenge that has been overlooked is that, as a
disease becomes more prevalent, it can limit the availability of the capital
needed to effectively treat those who have fallen ill. Here we use a simple
mathematical model to gain insight into the dynamics of an epidemic when the
recovery of sick individuals depends on the availability of healing resources
that are generated by the healthy population. We find that epidemics spiral out
of control into "explosive" spread if the cost of recovery is above a critical
cost. This can occur even when the disease would die out without the resource
constraint. The onset of explosive epidemics is very sudden, exhibiting a
discontinuous transition under very general assumptions. We find analytical
expressions for the critical cost and the size of the explosive jump in
infection levels in terms of the parameters that characterize the spreading
process. Our model and results apply beyond epidemics to contagion dynamics
that self-induce constraints on recovery, thereby amplifying the spreading
process.


  Conformational entropy for atomic-level, three dimensional biomolecules is
known experimentally to play an important role in protein-ligand
discrimination, yet reliable computation of entropy remains a difficult
problem. Here we describe the first two accurate and efficient algorithms to
compute the conformational entropy for RNA secondary structures, with respect
to the Turner energy model, where free energy parameters are determined from UV
aborption experiments. An algorithm to compute the derivational entropy for RNA
secondary structures had previously been introduced, using stochastic context
free grammars (SCFGs). However, the numerical value of derivational entropy
depends heavily on the chosen context free grammar and on the training set used
to estimate rule probabilities. Using data from the Rfam database, we determine
that both of our thermodynamic methods, which agree in numerical value, are
substantially faster than the SCFG method. Thermodynamic structural entropy is
much smaller than derivational entropy, and the correlation between
length-normalized thermodynamic entropy and derivational entropy is moderately
weak to poor. In applications, we plot the structural entropy as a function of
temperature for known thermoswitches, determine that the correlation between
hammerhead ribozyme cleavage activity and total free energy is improved by
including an additional free energy term arising from conformational entropy,
and plot the structural entropy of windows of the HIV-1 genome. Our software
RNAentropy can compute structural entropy for any user-specified temperature,
and supports both the Turner'99 and Turner'04 energy parameters. It follows
that RNAentropy is state-of-the-art software to compute RNA secondary structure
conformational entropy. The software is available at
http://bioinformatics.bc.edu/clotelab/RNAentropy.


  We discuss a novel analysis method for reaction network systems with
polynomial or rational rate functions. This method is based on computing
tropical equilibrations defined by the equality of at least two dominant
monomials of opposite signs in the differential equations of each dynamic
variable. In algebraic geometry, the tropical equilibration problem is
tantamount to finding tropical prevarieties, that are finite intersections of
tropical hypersurfaces. Tropical equilibrations with the same set of dominant
monomials define a branch or equivalence class. Minimal branches are
particularly interesting as they describe the simplest states of the reaction
network. We provide a method to compute the number of minimal branches and to
find representative tropical equilibrations for each branch.


  Background: The length of a protein sequence is largely determined by its
function, i.e. each functional group is associated with an optimal size.
However, comparative genomics revealed that proteins length may be affected by
additional factors. In 2002 it was shown that in bacterium Escherichia coli and
the archaeon Archaeoglobus fulgidus, protein sequences with no homologs are, on
average, shorter than those with homologs. Most experts now agree that the
length distributions are distinctly different between protein sequences with
and without homologs in bacterial and archaeal genomes. In this study, we
examine this postulate by a comprehensive analysis of all annotated prokaryotic
genomes and focusing on certain exceptions.
  Results: We compared lengths distributions of having homologs proteins (HHPs)
and non-having homologs proteins (orphans or ORFans) in all currently annotated
completely sequenced prokaryotic genomes. As expected, the HHPs and ORFans have
strikingly different length distributions in almost all genomes. As previously
established, the HHPs, indeed, are, on average, longer than the ORFans, and the
length distributions for the ORFans have a relatively narrow peak, in contrast
to the HHPs, whose lengths spread over a wider range of values. However, about
thirty genomes do not obey these rules. Practically all genomes of Mycoplasma
and Ureaplasma have atypical ORFans distributions, with the mean lengths of
ORFan larger than the mean lengths of HHPs. These genera constitute over 80% of
atypical genomes.
  Conclusions: We confirmed on a ubiquitous set of genomes the previous
observation that HHPs and ORFans have different gene length distributions. We
also showed that Mycoplasmataceae genomes have distinctive distributions of
ORFans lengths. We offer several possible biological explanations of this
phenomenon.


  Theoretical and experimental evidences support the hypothesis that extremely
low-frequency electromagnetic fields can affect voltage-gated channels. Little
is known, however, about their effect on potassium channels. Kv1.3, a member of
the voltage-gated potassium channels family originally discovered in the brain,
is a key player in important biological processes including antigen-dependent
activation of T-cells during the immune response. We report that Kv1.3
expressed in CHO-K1 cells can be modulated in cell subpopulations by extremely
low frequency and relatively low intensity electromagnetic fields. In
particular, we observed that field exposure can cause an enhancement of Kv1.3
potassium conductance and that the effect lasts for several minutes after field
removal. The results contribute to put immune and nervous system responses to
extremely low-frequency electromagnetic fields into a new perspective, with
Kv1.3 playing a pivotal molecular role.
  Keywords: immunotherapy, immunomodulation, potassium channels, gating,
electromagnetic fields


  The establishment of protected areas or parks has become an important tool
for wildlife conservation. However, frequent occurrences of human-wildlife
conflict at the edges of these parks can undermine their conservation goals.
Many African protected areas have experienced concurrent declines of apex
predators alongside increases in both baboon abundance and the density of
humans living near the park boundary. Baboons then take excursions outside of
the park to raid crops for food, conflicting with the human population. We
model the interactions of mesopredators (baboons), apex predators and shared
prey in the park to analyze how four components affect the proportion of time
that mesopredators choose to crop-raid: 1) the presence of apex predators; 2)
nutritional quality of the crops; 3) mesopredator "shyness" about leaving the
park; and 4) human hunting of mesopredators. We predict that the presence of
apex predators in the park is the most effective method for controlling
mesopredator abundance, and hence significantly reduces their impact on crops.
Human hunting of mesopredators is less effective as it only occurs during
crop-raiding excursions. Furthermore, making crops less attractive, for
instance by planting crops further from the park boundary or farming less
nutritional crops, can reduce the amount of time mesopredators crop-raid.


  Caffeine (CAF) is one of the most widely and regularly consumed biologically
active substances. We use computer simulation approach to the study of CAF
activity by searching for its possible complexes with biopolymer fragments. The
principal CAF target at physiologically important concentrations refers to
adenosine receptors. It is a common opinion that CAF is a competitive
antagonist of adenosine. At the first step to molecular level elucidation of
CAF action, we have found a set of the minima of the interaction energy between
CAF and the fragments of human A1 adenosine receptor. Molecular mechanics is
the main method for the calculations of the energy of interactions between CAF
and the biopolymer fragments. We use the Monte Carlo simulation to follow
various mutual arrangements of CAF molecule near the receptor. It appears that
the deepest energy minima refer to hydrogen-bond formation of CAF with amino
acid residues involved in interactions with adenosine, its agonists and
antagonists. The results suggest that the formation of such CAF-receptor
complexes enforced by a close packing of CAF and the receptor fragments is the
reason of CAF actions on nervous system. CAF can block the atomic groups of the
adenosine repressors responsible for the interactions with adenosine, not
necessary by the formation of H bonds with them, but simply hide these groups
from the interactions with adenosine.


  Species trees represent the historical divergences of populations or species,
while gene trees trace the ancestry of individual gene copies sampled within
those populations. In cases involving rapid speciation, gene trees with
topologies that differ from that of the species tree can be most probable under
the standard multispecies coalescent model, making species tree inference more
difficult. Such anomalous gene trees are not well understood except for some
small cases. In this work, we establish one constraint that applies to trees of
any size: gene trees with "caterpillar" topologies cannot be anomalous. The
proof of this involves a new combinatorial object, called a population history,
which keeps track of the number of coalescent events in each ancestral
population.


  Genome-wide association studies (GWAS) are widely used to discover genetic
variants associated with diseases. To control false positives, all findings
from GWAS need to be verified with additional evidences, even for associations
discovered from a high power study. Replication study is a common verification
method by using independent samples. An association is regarded as true
positive with a high confidence when it can be identified in both primary study
and replication study. Currently, there is no systematic study on the behavior
of positives in the replication study when the positive results of primary
study are considered as the prior information.
  In this paper, two probabilistic measures named Reproducibility Rate (RR) and
False Irreproducibility Rate (FIR) are proposed to quantitatively describe the
behavior of primary positive associations (i.e. positive associations
identified in the primary study) in the replication study. RR is a conditional
probability measuring how likely a primary positive association will also be
positive in the replication study. This can be used to guide the design of
replication study, and to check the consistency between the results of primary
study and those of replication study. FIR, on the contrary, measures how likely
a primary positive association may still be a true positive even when it is
negative in the replication study. This can be used to generate a list of
potentially true associations in the irreproducible findings for further
scrutiny. The estimation methods of these two measures are given. Simulation
results and real experiments show that our estimation methods have high
accuracy and good prediction performance.


  One of the most intriguing questions in evolution is how organisms exhibit
suitable phenotypic variation to rapidly adapt in novel selective environments
which is crucial for evolvability. Recent work showed that when selective
environments vary in a systematic manner, it is possible that development can
constrain the phenotypic space in regions that are evolutionarily more
advantageous. Yet, the underlying mechanism that enables the spontaneous
emergence of such adaptive developmental constraints is poorly understood. How
can natural selection, given its myopic and conservative nature, favour
developmental organisations that facilitate adaptive evolution in future
previously unseen environments? Such capacity suggests a form of
\textit{foresight} facilitated by the ability of evolution to accumulate and
exploit information not only about the particular phenotypes selected in the
past, but regularities in the environment that are also relevant to future
environments. Here we argue that the ability of evolution to discover such
regularities is analogous to the ability of learning systems to generalise from
past experience. Conversely, the canalisation of evolved developmental
processes to past selective environments and failure of natural selection to
enhance evolvability in future selective environments is directly analogous to
the problem of over-fitting and failure to generalise in machine learning. We
show that this analogy arises from an underlying mechanistic equivalence by
showing that conditions corresponding to those that alleviate over-fitting in
machine learning enhance the evolution of generalised developmental
organisations under natural selection. This equivalence provides access to a
well-developed theoretical framework that enables us to characterise the
conditions where natural selection will find general rather than particular
solutions to environmental conditions.


  Nucleosome positioning is an important process required for proper genome
packing and its accessibility to execute the genetic program in a
cell-specific, timely manner. In the recent years hundreds of papers have been
devoted to the bioinformatics, physics and biology of nucleosome positioning.
The purpose of this review is to cover a practical aspect of this field, namely
to provide a guide to the multitude of nucleosome positioning resources
available online. These include almost 300 experimental datasets of genome-wide
nucleosome occupancy profiles determined in different cell types and more than
40 computational tools for the analysis of experimental nucleosome positioning
data and prediction of intrinsic nucleosome formation probabilities from the
DNA sequence. A manually curated, up to date list of these resources will be
maintained at http://generegulation.info.


  Superparamagnetism of tryptophan implying the presence of magnetic domain is
reported. The observation helps us to conceive assembly of proteins as a
physical lattice gas with multidimensional Ising character, each lattice points
assuming discrete spin states. When magnetic field is applied the equilibrium
is lost and the population density of one spin state increases (unidirectional
alignment), resulting in net magnetization. Spatial coherence between the
identical spin states further imparts a ferromagnetic memory. This effect is
observed using direct nanoscale video imaging. Out of the three proteins
ferritin serum albumin and fibrinogen, fibrinogen showed an attenuated
response, the protein being essentially one dimensional. Eventually, Ising
lattice is capable of showing ferromagnetic memory only when it has a higher
dimensional character. The study highlights possible presence of long range
spatial coherence at physiological condition and a plausible microscopic origin
of the same.


  Collective sensing by interacting cells is observed in a variety of
biological systems, and yet a quantitative understanding of how sensory
information is collectively encoded is lacking. Here we investigate the
ATP-induced calcium dynamics of monolayers of fibroblast cells that communicate
via gap junctions. Combining experiments and stochastic modeling, we find that
increasing the ATP stimulus increases the propensity for calcium oscillations
despite large cell-to-cell variability. The model further predicts that the
oscillation propensity increases not only with the stimulus, but also with the
cell density due to increased communication. Experiments confirm this
prediction, showing that cell density modulates the collective sensory
response. We further implicate cell-cell communication by coculturing the
fibroblasts with cancer cells, which we show act as "defects" in the
communication network, thereby reducing the oscillation propensity. These
results suggest that multicellular networks sit at a point in parameter space
where cell-cell communication has a significant effect on the sensory response,
allowing cells to simultaneously respond to a sensory input and to the presence
of neighbors.


  The anchor of most integral membrane proteins consists of one or several
helices spanning the lipid bilayer. The WALP peptide, GWW(LA)$_n$(L)WWA, is a
common model helix to study the fundamentals of protein insertion and folding,
as well as helix-helix association in the membrane. Its structural properties
have been illuminated in a large number of experimental and simulation studies.
In this combined coarse-grained and atomistic simulation study, we probe the
thermodynamics of a single WALP peptide, focusing on both the insertion across
the water-membrane interface, as well as folding in both water and a membrane.
The potential of mean force characterizing the peptide's insertion into the
membrane shows qualitatively similar behavior across peptides and three force
fields. However, the Martini force field exhibits a pronounced secondary
minimum for an adsorbed interfacial state, which may even become the global
minimum---in contrast to both atomistic simulations and the alternative PLUM
force field. Even though the two coarse-grained models reproduce the free
energy of insertion of individual amino acids side chains, they both
underestimate its corresponding value for the full peptide (as compared with
atomistic simulations), hinting at cooperative physics beyond the residue
level. Folding of WALP in the two environments indicates the helix as the most
stable structure, though with different relative stabilities and chain-length
dependence.


  Evolutionary hypotheses to explain the greater numbers of species in the
tropics than the temperate zone include greater age and area, higher
temperature and metabolic rates, and greater ecological opportunity. These
ideas make contrasting predictions about the relationship between speciation
processes and latitude, which I elaborate and evaluate. Available data suggest
that per capita speciation rates are currently highest in the temperate zone,
and that diversification rates (speciation minus extinction) are similar
between latitudes. In contrast, clades whose oldest analyzed dates precede the
Eocene thermal maximum, when the extent of the tropics was much greater than
today, tend to show highest speciation and diversification rates in the
tropics. These findings are consistent with age and area, which is alone among
hypotheses in predicting a time trend. Higher recent speciation rates in the
temperate zone than the tropics suggest an additional response to high
ecological opportunity associated with low species diversity. These broad
patterns are compelling but provide limited insights into underlying
mechanisms, arguing that studies of speciation processes along the latitudinal
gradient will be vital. Using threespine stickleback in depauperate northern
lakes as an example, I show how high ecological opportunity can lead to rapid
speciation. The results support a role for ecological opportunity in
speciation, but its importance in the evolution of the latitudinal gradient
remains uncertain. I conclude that per-capita evolutionary rates are no longer
higher in the tropics than the temperate zone. Nevertheless, the vast numbers
of species that have already accumulated in the tropics ensure that total rate
of species production remains highest there. Thus, tropical evolutionary
momentum helps to perpetuate the steep latitudinal biodiversity gradient.


  Identifying genes associated with complex human diseases is one of the main
challenges of human genetics and computational medicine. To answer this
question, millions of genetic variants get screened to identify a few of
importance. To increase the power of identifying genes associated with diseases
and to account for other potential sources of protein function aberrations, we
propose a novel factor-graph based model, where much of the biological
knowledge is incorporated through factors and priors. Our extensive simulations
show that our method has superior sensitivity and precision compared to
variant-aggregating and differential expression methods. Our integrative
approach was able to identify important genes in breast cancer, identifying
genes that had coding aberrations in some patients and regulatory abnormalities
in others, emphasizing the importance of data integration to explain the
disease in a larger number of patients.


  Gene regulatory circuits must contend with intrinsic noise that arises due to
finite numbers of proteins. While some circuits act to reduce this noise,
others appear to exploit it. A striking example is the competence circuit in
Bacillus subtilis, which exhibits much larger noise in the duration of its
competence events than a synthetically constructed analog that performs the
same function. Here, using stochastic modeling and fluorescence microscopy, we
show that this larger noise allows cells to exit terminal phenotypic states,
which expands the range of stress levels to which cells are responsive and
leads to phenotypic heterogeneity at the population level. This is an important
example of how noise confers a functional benefit in a genetic decision-making
circuit.


  We construct a seven-component model of the in-host dynamics of the Human
Immunodeficiency Virus Type-1 (i.e, HIV) that accounts for latent infection and
the propensity of viral mutation. A dynamical analysis is conducted and a
theorem is presented which characterizes the long time behavior of the model.
Finally, we study the effects of an antiretroviral drug and treatment
implications.


  The MAP kinase cascade is an important signal transduction system in
molecular biology for which a lot of mathematical modelling has been done. This
paper surveys what has been proved mathematically about the qualitative
properties of solutions of the ordinary differential equations arising as
models for this biological system. It focusses, in particular, on the issues of
multistability and the existence of sustained oscillations. It also gives a
concise introduction to the mathematical techniques used in this context,
bifurcation theory and geometric singular perturbation theory, as they relate
to these specific examples. In addition further directions are presented in
which the applications of these techniques could be extended in the future.


  Evolution on neutral networks of genotypes has been found in models to
concentrate on genotypes with high mutational robustness, to a degree
determined by the topology of the network. Here analysis is generalized beyond
neutral networks to arbitrary selection and parent-offspring transmission. In
this larger realm, geometric features determine mutational robustness: the
alignment of fitness with the orthogonalized eigenvectors of the mutation
matrix weighted by their eigenvalues. "House of cards" mutation is found to
preclude the evolution of mutational robustness. Genetic load is shown to
increase with increasing mutation in arbitrary single and multiple locus
fitness landscapes. The rate of decrease in population fitness can never grow
as mutation rates get higher, showing that "error catastrophes" for genotype
frequencies never cause precipitous losses of population fitness. The
"inclusive inheritance" approach taken here naturally extends these results to
a new concept of dispersal robustness.


  The prediction of phenotypic traits using high-density genomic data has many
applications such as the selection of plants and animals of commercial
interest; and it is expected to play an increasing role in medical diagnostics.
Statistical models used for this task are usually tested using
cross-validation, which implicitly assumes that new individuals (whose
phenotypes we would like to predict) originate from the same population the
genomic prediction model is trained on. In this paper we propose an approach
based on clustering and resampling to investigate the effect of increasing
genetic distance between training and target populations when predicting
quantitative traits. This is important for plant and animal genetics, where
genomic selection programs rely on the precision of predictions in future
rounds of breeding. Therefore, estimating how quickly predictive accuracy
decays is important in deciding which training population to use and how often
the model has to be recalibrated. We find that the correlation between true and
predicted values decays approximately linearly with respect to either $\F$ or
mean kinship between the training and the target populations. We illustrate
this relationship using simulations and a collection of data sets from mice,
wheat and human genetics.


  We used statistical thermodynamics of conformational fluctuations and the
elements of algebraic graph theory together with data from 2000 protein crystal
structures, and showed that folded native proteins with harmonic interactions
exhibit distribution functions each of which appear to be universal across all
proteins. The three universal distributions are: (i) the eigenvalue spectrum of
the protein graph Laplacian, (ii) the B-factor distribution of residues, and
(iii) the vibrational entropy difference per residue between the unfolded and
the folded states. The three distributions, which look independent of each
other at first sight, are strongly associated with the Rouse chain model of a
polymer as the unfolded protein. We treat the folded protein as the strongly
perturbed state of the Rouse chain. We explain the underlying factors
controlling the three distributions and discuss the differences from those of
randomly folded structures.


  Endocytosis is an essential process by which cells internalize a piece of
plasma membrane and material from the outside. In cells with turgor, pressure
opposes membrane defor- mations, and increases the amount of force that has to
be generated by the endocytic machinery. To determine this force, and calculate
the shape of the membrane, we used physical theory to model an elastic surface
under pressure. Accurate fits of experimental profiles are obtained assuming
that the coated membrane is highly rigid and preferentially curved at the
endocytic site. The forces required from the actin machinery peaks at the onset
of deformation, indicating that once invagination has been initiated,
endocytosis is unlikely to stall before completion. Coat proteins do not lower
the initiation force but may affect the process by the curvature they induce.
In the presence of isotropic curvature inducers, pulling the tip of the
invagination can trigger the formation of a neck at the base of the
invagination. Hence direct neck constriction by actin may not be required,
while its pulling role is essential. Finally, the theory shows that anisotropic
curvature effectors stabilize membrane invaginations, and the loss of
crescent-shaped BAR domain proteins such as Rvs167 could therefore trigger
membrane scission.


  Mass campaigns with antimalarial drugs are potentially a powerful tool for
local elimination of malaria, yet current diagnostic technologies are
insufficiently sensitive to identify all individuals who harbor infections. At
the same time, overtreatment of uninfected individuals increases the risk of
accelerating emergence of drug resistance and losing community acceptance.
Local heterogeneity in transmission intensity may allow campaign strategies
that respond to index cases to successfully target subpatent infections while
simultaneously limiting overtreatment. While selective targeting of hotspots of
transmission has been proposed as a strategy for malaria control, such
targeting has not been tested in the context of malaria elimination. Using
household locations, demographics, and prevalence data from a survey of four
health facility catchment areas in southern Zambia and an agent-based model of
malaria transmission and immunity acquisition, a transmission intensity was fit
to each household based on neighborhood age-dependent malaria prevalence. A set
of individual infection trajectories was constructed for every household in
each catchment area, accounting for heterogeneous exposure and immunity.
Various campaign strategies (mass drug administration, mass screen and treat,
focal mass drug administration, snowball reactive case detection, pooled
sampling, and a hypothetical serological diagnostic) were simulated and
evaluated for performance at finding infections, minimizing overtreatment,
reducing clinical case counts, and interrupting transmission. For malaria
control, presumptive treatment leads to substantial overtreatment without
additional morbidity reduction under all but the highest transmission
conditions. Selective targeting of hotspots with drug campaigns is an
ineffective tool for elimination due to limited sensitivity of available field
diagnostics.


  Infectious diseases often spread faster near their peak than would be
predicted given early data on transmission. Despite the commonality of this
phenomena, there are no known general mechanisms able to cause an exponentially
spreading dis- ease to begin spreading faster. Indeed most features of real
world social networks, e.g. clustering1,2 and community structure3, and of
human behaviour, e.g. social distancing4 and increased hygiene5, will slow
disease spread. Here, we consider a model where individuals with essential
societal roles-e.g. teachers, first responders, health-care workers, etc.- who
fall ill are replaced with healthy individuals. We refer to this process as
relational exchange. Relational exchange is also a behavioural process, but one
whose effect on disease transmission is less obvious. By incorporating this
behaviour into a dynamic network model, we demonstrate that replacing
individuals can accelerate disease transmission. Furthermore, we find that the
effects of this process are trivial when considering a standard mass-action
model, but dramatic when considering network structure. This result highlights
another critical shortcoming in mass-action models, namely their inability to
account for behavioural processes. Lastly, using empirical data, we find that
this mechanism parsimoniously explains observed patterns across more than
seventeen years of influenza and dengue virus data. We anticipate that our
findings will advance the emerging field of disease forecasting and will better
inform public health decision making during outbreaks.


  Chloroplasts regulate their growth to optimize photosynthesis. Quantitative
data shows that the ratio of total chloroplast area to mesophyll cell area is
constant across different cells within a single species, and also across
species. Wild-type chloroplasts exhibit little scatter around this trend;
highly irregularly-shaped mutant chloroplasts exhibit more scatter. Here we
propose a model motivated by a bacterial quorum-sensing model consisting of a
switch-like signalling network that turns off chloroplast growth. We calculated
the dependence of the location of the relevant saddle-node bifurcation on the
geometry of the chloroplasts. Our model exhibits a linear trend, with linearly
growing scatter dependent on chloroplast shape, consistent with the data. When
modelled chloroplasts are of a shape that grows with a constant area to volume
ratio (disks, cylinders) we find a linear trend with minimal scatter.
Chloroplasts with area and volume that do not grow proportionally (spheres)
exhibit a linear trend with additional scatter.


  Rule-based modeling is a powerful way to model kinetic interactions in
biochemical systems. Rules enable a precise encoding of biochemical
interactions at the resolution of sites within molecules, but obtaining an
integrated global view from sets of rules remains challenging. Current
automated approaches to rule visualization fail to address the complexity of
interactions between rules, limiting either the types of rules that are allowed
or the set of interactions that can be visualized simultaneously. There is a
need for scalable visualization approaches that present the information encoded
in rules in an intuitive and useful manner at different levels of detail. We
have developed new automated approaches for visualizing both individual rules
and complete rule-based models. We find that a more compact representation of
an individual rule promotes promotes understanding the model assumptions
underlying each rule. For global visualization of rule interactions, we have
developed a method to synthesize a network of interactions between sites and
processes from a rule-based model and then use a combination of user-defined
and automated approaches to compress this network into a readable form. The
resulting diagrams enable modelers to identify signaling motifs such as
cascades, feedback loops, and feed-forward loops in complex models, as we
demonstrate using several large-scale models. These capabilities are
implemented within the BioNetGen framework but the approach is equally
applicable to rule-based models specified in other formats.


  We study the folding process in the shallowly knotted protein MJ0366 within
two variants of a structure-based model. We observe that the resulting
topological pathways are much richer than identified in previous studies. In
addition to the single knot-loop events, we find novel, and dominant, two-loop
mechanisms. We demonstrate that folding takes place in a range of temperatures
and the conditions of most successful folding are at temperatures which are
higher than those required for the fastest folding. We also demonstrate that
nascent conditions are more favorable to knotting than off-ribosome folding.


  Proper folding of deeply knotted proteins has a very low success rate even in
structure-based models which favor formation of the native contacts but have no
topological bias. By employing a structure-based model, we demonstrate that
cotranslational folding on a model ribosome may enhance the odds to form
trefoil knots for protein YibK without any need to introduce any non-native
contacts. The ribosome is represented by a repulsive wall that keeps elongating
the protein. On-ribosome folding proceeds through a a slipknot conformation. We
elucidate the mechanics and energetics of its formation. We show that the
knotting probability in on-ribosome folding is a function of temperature and
that there is an optimal temperature for the process. Our model often leads to
the establishment of the native contacts without formation of the knot.


  Cellular processes do not follow deterministic rules; even in identical
environments genetically identical cells can make random choices leading to
different phenotypes. This randomness originates from fluctuations present in
the biomolecular interaction networks. Most previous work has been focused on
the intrinsic noise (IN) of these networks. Yet, especially for
high-copy-number biomolecules, extrinsic or environmental noise (EN) has been
experimentally shown to dominate the variation. Here, we develop an analytical
formalism that allows for calculation of the effect of EN on gene-expression
motifs. We introduce a method for modeling bounded EN as an auxiliary species
in the master equation. The method is fully generic and is not limited to
systems with small EN magnitudes. We focus our study on motifs that can be
viewed as the building blocks of genetic switches: a nonregulated gene, a
self-inhibiting gene, and a self-promoting gene. The role of the EN properties
(magnitude, correlation time, and distribution) on the statistics of interest
are systematically investigated, and the effect of fluctuations in different
reaction rates is compared. Due to its analytical nature, our formalism can be
used to quantify the effect of EN on the dynamics of biochemical networks and
can also be used to improve the interpretation of data from single-cell
gene-expression experiments.


  The question as to why most higher organisms reproduce sexually has remained
open despite extensive research, and has been called "the queen of problems in
evolutionary biology". Theories dating back to Weismann have suggested that the
key must lie in the creation of increased variability in offspring, causing
enhanced response to selection. Rigorously quantifying the effects of assorted
mechanisms which might lead to such increased variability, and establishing
that these beneficial effects outweigh the immediate costs of sexual
reproduction has, however, proved problematic. Here we introduce an approach
which does not focus on particular mechanisms influencing factors such as the
fixation of beneficial mutants or the ability of populations to deal with
deleterious mutations, but rather tracks the entire distribution of a
population of genotypes as it moves across vast fitness landscapes. In this
setting simulations now show sex robustly outperforming asex across a broad
spectrum of finite or infinite population models. Concentrating on the additive
infinite populations model, we are able to give a rigorous mathematical proof
establishing that sexual reproduction acts as a more efficient optimiser of
mean fitness, thereby solving the problem for this model. Some of the key
features of this analysis carry through to the finite populations case.


  Deciphering the control principles of metabolism and its interaction with
other cellular functions is central to biomedicine and biotechnology. Yet,
understanding the efficient control of metabolic fluxes remains elusive for
large-scale metabolic networks. Existing methods either require specifying a
cellular objective or are limited to small networks due to computational
complexity. Here we develop an efficient computational framework for flux
control by introducing a complete set of flux coupling relations. We analyze 23
metabolic networks from all kingdoms of life, and identify the driver reactions
facilitating their control on a large scale. We find that most unicellular
organisms require less extensive control than multicellular organisms. The
identified driver reactions are under strong transcriptional regulation in
Escherichia coli. In human cancer cells driver reactions play pivotal roles in
tumor development, representing potential therapeutic targets. The proposed
framework helps us unravel the regulatory principles of complex diseases and
design novel engineering strategies at the interface of gene regulation,
signaling, and metabolism.


  We will study a population of individuals playing the infinitely repeated
Prisoner's Dilemma under replicator dynamics. The population consists of three
kinds of individuals using the following reactive strategies: ALLD (individuals
which always defect), ATFT (almost tit-for-tat: individuals which almost always
repeat the opponent's last move) and G (generous individuals, which always
cooperate when the opponent cooperated in the last move and have a positive
probability $q$ of cooperating when they are defected). Our aim is studying in
a mathematically rigorous fashion the dynamics of a simplified version for the
computer experiment in [Nowak, Sigmund, Nature, 355, pp. 250--53, 1992]
involving 100 reactive strategies. We will see that as the generosity degree of
the G individuals varies, equilibria (rest points) of the dynamics appear or
disappear, and the dynamics changes accordingly. Not only we will prove that
the results of the experiment are true in our simplified version, but we will
have complete control on the existence or non-existence of the equilbria for
the dynamics for all possible values of the parameters, given that ATFT
individuals are close enough to TFT. For most values of the parameters the
dynamics will be completely determined.


  1) Micro-evolutionary predictions are complicated by ecological feedbacks
like density dependence, while ecological predictions can be complicated by
evolutionary change. A widely used approach in micro-evolution, quantitative
genetics, struggles to incorporate ecological processes into predictive models,
while structured population modelling, a tool widely used in ecology, rarely
incorporates evolution explicitly.
  2) In this paper we develop a flexible, general framework that links
quantitative genetics and structured population models. We use the quantitative
genetic approach to write down the phenotype as an additive map. We then
construct integral projection models for each component of the phenotype. The
dynamics of the distribution of the phenotype are generated by combining
distributions of each of its components. Population projection models can be
formulated on per generation or on shorter time steps.
  3) We introduce the framework before developing example models with
parameters chosen to exhibit specific dynamics. These models reveal (i) how
evolution of a phenotype can cause populations to move from one dynamical
regime to another (e.g. from stationarity to cycles), (ii) how additive genetic
variances and covariances (the G matrix) are expected to evolve over multiple
generations, (iii) how changing heritability with age can maintain additive
genetic variation in the face of selection and (iii) life history, population
dynamics, phenotypic characters and parameters in ecological models will change
as adaptation occurs.
  4) Our approach unifies population ecology and evolutionary biology providing
a framework allowing a very wide range of questions to be addressed. The next
step is to apply the approach to a variety of laboratory and field systems.
Once this is done we will have a much deeper understanding of eco-evolutionary
dynamics and feedbacks.


  Understanding the conditions ensuring the persistence of a population is an
issue of primary importance in population biology. The first theoretical
approach to the problem dates back to the 50's with the KiSS (after Kierstead,
Slobodkin and Skellam) model, namely a continuous reaction-diffusion equation
for a population growing on a patch of finite size $L$ surrounded by a deadly
environment with infinite mortality -- i.e. an oasis in a desert. The main
outcome of the model is that only patches above a critical size allow for
population persistence. Here, we introduce an individual-based analogue of the
KiSS model to investigate the effects of discreteness and demographic
stochasticity. In particular, we study the average time to extinction both
above and below the critical patch size of the continuous model and investigate
the quasi-stationary distribution of the number of individuals for patch sizes
above the critical threshold.


  Understanding what controls apparently random fluctuations in fish
recruitment is a major challenge in fisheries science. Our current inability to
anticipate recruitment failures has led to costly management actions and in
some cases complete fishery collapse. Time series observations of fish
recruitment reflect an interplay of underlying population processes,
environmental forcing, measurement error, and process error. Given that the
error component is often very strong, an important unresolved question is
whether any non-random signal can be uncovered in the annual fluctuations of
recruitment time series. Here, we address this fundamental question in an
analysis of 569 fish populations from a global database of recruitment. Using a
nonparametric time series analysis method, we find overwhelming evidence for
non-random dynamics operating in the recruitment process. Unlike previous
explorations of this topic, our approach does not require the specification of
a stock-recruitment model, and it can be used on recruitment estimates derived
from a wide-range of abundance estimation methods. The evidence for
non-randomness is robust across a wide range of fish families and abundance
estimation methods. We also find that the statistical support for
non-randomness increases consistently with the number of observations of a fish
stock. This result provides the encouraging news that with continued
observations, and the appropriate covariates, we should ultimately be able to
uncover the mechanistic drivers of fish recruitment.


  Malaria is a mosquito-borne, lethal disease that affects millions and kills
hundreds of thousands of people each year. In this paper, we develop a model
for allocating malaria interventions across geographic regions and time,
subject to budget constraints, with the aim of minimizing the number of
person-days of malaria infection. The model considers a range of several
conditions: climatic characteristics, treatment efficacy, distribution costs,
and treatment coverage. We couple an expanded susceptible-infected-recovered
(SIR) compartment model for the disease dynamics with an integer linear
programming (ILP) model for selecting the disease interventions. Our model
produces an intervention plan for all regions, identifying which combination of
interventions, with which level of coverage, to use in each region and year in
a five-year planning horizon. Simulations using the model yield high-level,
qualitative insights on optimal intervention policies: The optimal policy is
different when considering a five-year time horizon than when considering only
a single year, due to the effects that interventions have on the disease
transmission dynamics. The vaccine intervention is rarely selected, except if
its assumed cost is significantly lower than that predicted in the literature.
Increasing the available budget causes the number of person-days of malaria
infection to decrease linearly up to a point, after which the benefit of
increased budget starts to taper. The optimal policy is highly dependent on
assumptions about mosquito density, selecting different interventions for wet
climates with high density than for dry climates with low density, and the
interventions are found to be less effective at controlling malaria in the wet
climates when attainable intervention coverage is 60% or lower. However, when
intervention coverage of 80% is attainable, then malaria prevalence drops
quickly.


  A large class of phylogenetic networks can be obtained from trees by the
addition of horizontal edges between the tree edges. These networks are called
tree based networks. Reticulation-visible networks and child-sibling networks
are all tree based. In this work, we present a simply necessary and sufficient
condition for tree-based networks and prove that there is a universal tree
based network for each set of species such that every phylogenetic tree on the
same species is a base of this network. The existence of universal tree based
network implies that for any given set of phylogenetic trees (resp. clusters)
on the same species there exists a tree base network that display all of them.


  In many models of the Calvin cycle of photosynthesis it is observed that
there are solutions where concentrations of key substances belonging to the
cycle tend to zero at late times, a phenomenon known as overload breakdown. In
this paper we prove theorems about the existence and non-existence of solutions
of this type and obtain information on which concentrations tend to zero when
overload breakdown occurs. As a starting point we take a model of Pettersson
and Ryde-Pettersson which seems to be prone to overload breakdown and a
modification of it due to Poolman which was intended to avoid this effect.


  Ecologists and evolutionary biologists are relying on an increasingly
sophisticated set of statistical tools to describe complex natural systems. One
such tool that has gained increasing traction in the life sciences is
structural equation modeling (SEM), a variant of path analysis that resolves
complex multivariate relationships among a suite of interrelated variables. SEM
has historically relied on covariances among variables, rather than the values
of the data points themselves. While this approach permits a wide variety of
model forms, it limits the incorporation of detailed specifications. Here, I
present a fully-documented, open-source R package piecewiseSEM that builds on
the base R syntax for all current generalized linear, least-square, and mixed
effects models. I also provide two worked examples: one involving a
hierarchical dataset with non-normally distributed variables, and a second
involving phylogenetically-independent contrasts. My goal is to provide a
user-friendly and tractable implementation of SEM that also reflects the
ecological and methodological processes generating data.


  By means of the asynchronous cellular automata algorithm we study stationary
states and spatial patterning in an $SIS$ model, in which the individuals' are
attached to the vertices of a graph and their mobility is mimicked by varying
the neighbourhood size $q$. The versions with fixed $q$ and those taken at
random at each step and for each individual are studied. Numerical data on the
local behaviour of the model are mapped onto the solution of its zero
dimensional version, corresponding to the limit $q\to +\infty$ and equivalent
to the logistic growth model. This allows for deducing an explicit form of the
dependence of the fraction of infected individuals on the curing rate $\gamma$.
A detailed analysis of the appearance of spatial patterns of infected
individuals in the stationary state is performed.


  This paper analyzes, in the context of a prokaryotic cell, the stochastic
variability of the number of proteins when there is a control of gene
expression by an autoregulation scheme. The goal of this work is to estimate
the efficiency of the regulation to limit the fluctuations of the number of
copies of a given protein. The autoregulation considered in this paper relies
mainly on a negative feedback: the proteins are repressors of their own gene
expression. The efficiency of a production process without feedback control is
compared to a production process with an autoregulation of the gene expression
assuming that both of them produce the same average number of proteins. The
main characteristic used for the comparison is the standard deviation of the
number of proteins at equilibrium. With a Markovian representation and a simple
model of repression, we prove that, under a scaling regime, the repression
mechanism follows a Hill repression scheme with an hyperbolic control. An
explicit asymptotic expression of the variance of the number of proteins under
this regulation mechanism is obtained. Simulations are used to study other
aspects of autoregulation such as the rate of convergence to equilibrium of the
production process and the case where the control of the production process of
proteins is achieved via the inhibition of mRNAs.


  We present a framework to transform the problem of finding a Lyapunov
function of a Chemical Reaction Network (CRN) in concentration coordinates with
arbitrary monotone kinetics into finding a common Lyapunov function for a
linear parameter varying system in reaction coordinates. Alternative
formulations of the proposed Lyapunov function is presented also. This is
applied to reinterpret previous results by the authors on Piecewise Linear in
Rates Lyapunov functions, and to establish a link with contraction analysis.
Persistence and uniqueness of equilibria are discussed also.


  Heterozygote disadvantage is potentially a potent driver of population
genetic divergence. Also referred to as underdominance, this phenomena
describes a situation where a genetic heterozygote has a lower overall fitness
than either homozygote. Attention so far has mostly been given to
underdominance within a single population and the maintenance of genetic
differences between two populations exchanging migrants. Here we explore the
dynamics of an underdominant system in a network of multiple discrete, yet
interconnected, populations. Stability of genetic differences in response to
increases in migration in various topological networks is assessed. The network
topology can have a dominant and occasionally non-intuitive influence on the
genetic stability of the system.


  Alzheimer's disease is the most common form of dementia in the western world,
however there is no cure available for this devastating neurodegenerative
disorder. Despite clinical and experimental evidence implicating the intestinal
microbiota in a number of brain disorders, its impact on Alzheimer's disease is
not known. We generated a germ-free mouse model of Alzheimer's disease and
discovered a drastic reduction of cerebral Ab amyloid pathology when compared
to control Alzheimer's disease animals with intestinal microbiota. Sequencing
bacterial 16S rRNA from fecal samples revealed a remarkable shift in the gut
microbiota of conventionally-raised Alzheimer's disease mice as compared to
healthy, wild-type mice. Colonization of germ-free Alzheimer mice with
harvested microbiota from conventionally-raised Alzheimer mice dramatically
increased cerebral Ab pathology. In contrast, colonization with microbiota from
control wild-type mice was ineffective in increasing cerebral Ab levels. Our
results indicate a microbial involvement in the development of Alzheimer's
disease pathology, and suggest that microbiota may contribute to the
development of neurodegenerative diseases.


  Recovering a stochastic process from noisy ensembles of single particle
trajectories (SPTs) is resolved here using the Langevin equation as a model.
The massive redundancy contained in SPTs data allows recovering local
parameters of the underlying physical model. We use several parametric and
non-parametric estimators to compute the first and second moment of the process
and to recover the local drift, its derivative and the diffusion tensor. Using
a local asymptotic expansion of the estimators and computing the empirical
transition probability function, we develop here a method to deconvolve the
instrumental from the physical noise. We use numerical simulations to explore
the range of validity for the estimators. The present analysis allows
characterizing what can exactly be recovered from the statistics of
super-resolution microscopy trajectories used in molecular trafficking and
underlying cellular function.


  Motor proteins display widely different stepping patterns as they move on
microtubule tracks, from the deterministic linear or helical motion performed
by the protein kinesin to the uncoordinated random steps made by dynein. How
these different strategies produce an efficient navigation system needed to
ensure correct cellular functioning is still unclear. Here, we show by
numerical simulations that deterministic and random motor steps yield different
outcomes when random obstacles decorate the microtubule tracks: kinesin moves
faster on clean tracks but its motion is strongly hindered on decorated tracks,
while dynein is slower on clean tracks but more efficient in avoiding
obstacles. Further simulations indicate that dynein's advantage on decorated
tracks is due to its ability to step backwards. Our results explain how
different navigation strategies are employed by the cell to optimize motor
driven cargo transport.


  The unconstrained genomic DNA of bacteria forms a coil, which volume exceeds
1000 times the volume of the cell. Since prokaryotes lack a membrane-bound
nucleus, in sharp contrast with eukaryotes, the DNA may consequently be
expected to occupy the whole available volume when constrained to fit in the
cell. Still, it has been known for more than half a century that the DNA is
localized in a well defined region of the cell, called the nucleoid, which
occupies only 15% to 25% of the total volume. Although this problem has focused
the attention of many scientists for the past decades, there is still no
certainty concerning the mechanism that enables such a dramatic compaction. The
goal of this Topical Review is to take stock of our knowledge on this question
by listing all possible compaction mechanisms with the proclaimed desire to
clarify the physical principles they are based upon and discuss them in the
light of experimental results and the results of simulations based on
coarse-grained models. In particular, the fundamental differences between
psi-condensation and segregative phase separation and between the condensation
by small and long polycations are highlighted. This review suggests that the
importance of certain mechanisms, like supercoiling and the architectural
properties of DNA-bridging and DNA-bending nucleoid proteins, may have been
overestimated, whereas other mechanisms, like segregative phase separation and
the self-association of nucleoid proteins, as well as the possible role of the
synergy of two or more mechanisms, may conversely deserve more attention.


  Variation and selection are the core principles of Darwinian evolution, yet
quantitatively relating the diversity of a population to its capacity to
respond to selection is challenging. Here, we examine this problem at a
molecular level in the context of populations of partially randomized proteins
selected for binding to well-defined targets. We built several minimal protein
libraries, screened them in vitro by phage display and analyzed their response
to selection by high-throughput sequencing. A statistical analysis of the
results reveals two main findings: first, libraries with same sequence
diversity but built around different "frameworks" typically have vastly
different responses, second, the distribution of responses within a library
follows a simple scaling law. We show how an elementary probabilistic model
based on extreme value theory rationalizes these findings. Our results have
implications for designing synthetic protein libraries, for estimating the
density of functional biomolecules in sequence space, for characterizing
diversity in natural populations and for experimentally investigating the
concept of evolvability, or potential for future evolution.


  Many microbial populations rapidly adapt to changing environments with
multiple variants competing for survival. To quantify such complex evolutionary
dynamics in vivo, time resolved and genome wide data including rare variants
are essential. We performed whole-genome deep sequencing of HIV-1 populations
in 9 untreated patients, with 6-12 longitudinal samples per patient spanning
5-8 years of infection. We show that patterns of minor diversity are
reproducible between patients and mirror global HIV-1 diversity, suggesting a
universal landscape of fitness costs that control diversity. Reversions towards
the ancestral HIV-1 sequence are observed throughout infection and account for
almost one third of all sequence changes. Reversion rates depend strongly on
conservation. Frequent recombination limits linkage disequilibrium to about
100bp in most of the genome, but strong hitch-hiking due to short range linkage
limits diversity.


  Thousands of candidate human-specific regulatory sequences (HSRS) have been
identified, supporting the hypothesis that unique to human phenotypes result
from human-specific alterations of genomic regulatory networks. Here,
conservation patterns analysis of 18,364 candidate HSRS was carried out based
on definition of the sequence conservation threshold as the minimum ratio of
bases that must remap of 1.00. A total of 5,535 candidate HSRS were identified
that are: i) highly conserved in Great Apes; ii) evolved by the exaptation of
highly conserved ancestral DNA; iii) defined by either the acceleration of
mutation rates on the human lineage or the functional divergence from nonhuman
primates. The exaptation of highly conserved ancestral DNA pathway seems
mechanistically distinct from the evolution of regulatory DNA segments driven
by the species-specific expansion of transposable elements. Present analysis
supports the idea that phenotypic divergence of Homo sapiens is driven by the
evolution of human-specific genomic regulatory networks via two mechanistically
distinct pathways of creation of divergent sequences of regulatory DNA: i)
exaptation of the highly conserved ancestral regulatory DNA segments; ii)
human-specific insertions of transposable elements.


  I present a simple numerical model based on iteratively updating subgroups of
a population, individually modeled by nonnegative real numbers, by a constant
decay factor; however, at each iteration, one group is selected to instead be
updated by a constant growth factor. I discover a relationship between these
variables and their respective probabilities for a given subgroup, summarized
as the variable $c$. When $c>1$, the subgroup is found to tend towards
behaviors reminiscent of antifragility; when at least one subgroup of the
population has $c\ge1$, the population as a whole tends towards significantly
higher probabilities of "living forever," although it may first suffer a drop
in population size as less robust, fragile subgroups "die off." In concluding,
I discuss the limitations and ethics of such a model, notably the implications
of when an upper limit is placed on the growth constant, requiring a population
to facilitate an increase in the decay factor to lessen the impact of periods
of failure.


  The MAP kinase cascade is a network of enzymatic reactions arranged in
layers. In each layer occurs a multiple futile cycle of phosphorylations. The
fully phosphorylated substrate then serves as an enzyme for the layer below.
This papers focusses on the existence of parameters for which Hopf bifurcations
occur and generate periodic orbits. Furthermore it is explained how geometric
singular perturbation theory allows to generalize results from simple models to
more complex ones.


  Inclusive Fitness Theory (IFT) was proposed half a century ago by W.D.
Hamilton to explain the emergence and maintenance of cooperation between
individuals that allows the existence of society. Contemporary evolutionary
ecology identified several factors that increase inclusive fitness, in addition
to kin-selection, such as assortation or homophily, and social synergies
triggered by cooperation. Here we propose an Extend Inclusive Fitness Theory
(EIFT) that includes in the fitness calculation all direct and indirect
benefits an agent obtains by its own actions, and through interactions with kin
and with genetically unrelated individuals. This formulation focuses on the
sustainable cost/benefit threshold ratio of cooperation and on the probability
of agents sharing mutually compatible memes or genes. This broader description
of the nature of social dynamics allows to compare the evolution of cooperation
among kin and non-kin, intra- and inter-specific cooperation, co-evolution, the
emergence of symbioses, of social synergies, and the emergence of division of
labor. EIFT promotes interdisciplinary cross fertilization of ideas by allowing
to describe the role for division of labor in the emergence of social
synergies, providing an integrated framework for the study of both, biological
evolution of social behavior and economic market dynamics.


  Active stabilization in systems with zero or negative stiffness is an
essential element of a wide variety of biological processes. We study a
prototypical example of this phenomenon at a micro-scale and show how active
rigidity, interpreted as a formation of a pseudo-well in the effective energy
landscape, can be generated in an overdamped ratchet-type stochastic system. We
link the transition from negative to positive rigidity with correlations in the
noise and show that subtle differences in out-of-equilibrium driving may
compromise the emergence of a pseudo-well.


  Genome Wide Association Studies (GWAS) and eQTL analyses have produced a
large and growing number of genetic associations linked to a wide range of
human phenotypes. As of 2013, there were more than 11,000 SNPs associated with
a trait as reported in the NHGRI GWAS Catalog. However, interpreting the
functional roles played by these SNPs remains a challenge. Here we describe an
approach that uses the inherent bipartite structure of eQTL networks to place
SNPs into a functional context.
  Using genotyping and gene expression data from 163 lung tissue samples in a
study of Chronic Obstructive Pulmonary Disease (COPD) we calculated eQTL
associations between SNPs and genes and cast significant associations (FDR $<
0.1$) as links in a bipartite network. To our surprise, we discovered that the
highly-connected "hub" SNPs within the network were devoid of
disease-associations. However, within the network we identified 35 highly
modular communities, which comprise groups of SNPs associated with groups of
genes; 13 of these communities were significantly enriched for distinct
biological functions (P $ < 5 \times 10^{-4}$) including COPD-related
functions. Further, we found that GWAS-significant SNPs were enriched at the
cores of these communities, including previously identified GWAS associations
for COPD, asthma, and pulmonary function, among others. These results speak to
our intuition: rather than single SNPs influencing single genes, we see groups
of SNPs associated with the expression of families of functionally related
genes and that disease SNPs are associated with the perturbation of those
functions. These methods are not limited in their application to COPD and can
be used in the analysis of a wide variety of disease processes and other
phenotypic traits.


  Scientists have long sought to understand how vascular networks supply blood
and oxygen to cells throughout the body. Recent work focuses on principles that
constrain how vessel size changes through branching generations from the aorta
to capillaries and uses scaling exponents to quantify these changes. Prominent
scaling theories predict that combinations of these exponents explain how
metabolic, growth, and other biological rates vary with body size.
Nevertheless, direct measurements of individual vessel segments have been
limited because existing techniques for measuring vasculature are invasive,
time consuming, and technically difficult. We developed software that extracts
the length, radius, and connectivity of in vivo vessels from contrast-enhanced
3D Magnetic Resonance Angiography. Using data from 20 human subjects, we
calculated scaling exponents by four methods--two derived from local properties
of branching junctions and two from whole-network properties. Although these
methods are often used interchangeably in the literature, we do not find
general agreement between these methods, particularly for vessel lengths.
Measurements for length of vessels also diverge from theoretical values, but
those for radius show stronger agreement. Our results demonstrate that vascular
network models cannot ignore certain complexities of real vascular systems and
indicate the need to discover new principles regarding vessel lengths.


  In a complex system, the individual components are neither so tightly coupled
or correlated that they can all be treated as a single unit, nor so
uncorrelated that they can be approximated as independent entities. Instead,
patterns of interdependency lead to structure at multiple scales of
organization. Evolution excels at producing such complex structures. In turn,
the existence of these complex interrelationships within a biological system
affects the evolutionary dynamics of that system. I present a mathematical
formalism for multiscale structure, grounded in information theory, which makes
these intuitions quantitative, and I show how dynamics defined in terms of
population genetics or evolutionary game theory can lead to multiscale
organization. For complex systems, "more is different," and I address this from
several perspectives. Spatial host--consumer models demonstrate the importance
of the structures which can arise due to dynamical pattern formation.
Evolutionary game theory reveals the novel effects which can result from
multiplayer games, nonlinear payoffs and ecological stochasticity. Replicator
dynamics in an environment with mesoscale structure relates to generalized
conditionalization rules in probability theory. The idea of natural selection
"acting at multiple levels" has been mathematized in a variety of ways, not all
of which are equivalent. We will face down the confusion, using the experience
developed over the course of this thesis to clarify the situation.


  We present a modelling framework, and basic model parameterization, for the
study of DNA origami folding at the level of DNA domains. Our approach is
explicitly kinetic and does not assume a specific folding pathway. The binding
of each staple is associated with a free-energy change that depends on staple
sequence, the possibility of coaxial stacking with neighbouring domains, and
the entropic cost of constraining the scaffold by inserting staple crossovers.
A rigorous thermodynamic model is difficult to implement as a result of the
complex, multiply connected geometry of the scaffold: we present a solution to
this problem for planar origami. Coaxial stacking and entropic terms,
particularly when loop closure exponents are taken to be larger than those for
ideal chains, introduce interactions between staples. These cooperative
interactions lead to the prediction of sharp assembly transitions with notable
hysteresis that are consistent with experimental observations. We show that the
model reproduces the experimentally observed consequences of reducing staple
concentration, accelerated cooling and absent staples. We also present a
simpler methodology that gives consistent results and can be used to study a
wider range of systems including non-planar origami.


  A new general route for siRNA delivery is presented combining porous
core-shell silica nanocarriers with a modularly designed multifunctional block
copolymer. Specifically, the internal storage and release of siRNA from
mesoporous silica nanoparticles (MSN) with orthogonal core-shell surface
chemistry was investigated as a function of pore-size, pore morphology, surface
properties and pH. Very high siRNA loading capacities of up to 380 microg/mg
MSN were obtained with charge-matched amino-functionalized mesoporous cores,
and release profiles show up to 80% siRNA elution after 24 h. We demonstrate
that adsorption and desorption of siRNA is mainly driven by electrostatic
interactions, which allow for high loading capacities even in medium-sized
mesopores with pore diameters down to 4 nm in a stellate pore morphology. The
negatively charged MSN shell enabled the association with a block copolymer
containing positively charged artificial amino acids and oleic acid blocks,
which acts simultaneously as capping function and endosomal release agent. The
potential of this multifunctional delivery platform is demonstrated by highly
effective cell transfection and siRNA delivery into KB-cells. A luciferase
reporter gene knock-down of up to 90% was possible using extremely low cell
exposures with only 2.5 microg MSN containing 32 pM siRNA per 100 microL well.


  The quasi-steady state approximation and time-scale separation are commonly
applied methods to simplify models of biochemical reaction networks based on
ordinary differential equations (ODEs). The concentrations of the "fast"
species are assumed effectively to be at steady state with respect to the
"slow" species. Under this assumption the steady state equations can be used to
eliminate the "fast" variables and a new ODE system with only the slow species
can be obtained.
  We interpret a reduced system obtained by time-scale separation as the ODE
system arising from a unique reaction network, by identification of a set of
reactions and the corresponding rate functions. The procedure is graphically
based and can easily be worked out by hand for small networks. For larger
networks, we provide a pseudo-algorithm. We study properties of the reduced
network, its kinetics and conservation laws, and show that the kinetics of the
reduced network fulfil realistic assumptions, provided the original network
does. We illustrate our results using biological examples such as substrate
mechanisms, post-translational modification systems and networks with
intermediates (transient) steps.


  We employ partial integro-differential equations to model trophic interaction
in a spatially extended heterogeneous environment. Compared to classical
reaction-diffusion models, this framework allows us to more realistically
describe the situation where movement of individuals occurs on a faster time
scale than the demographic (population) time scale, and we cannot determine
population growth based on local density. However, most of the results reported
so far for such systems have only been verified numerically and for a
particular choice of model functions, which obviously casts doubts about these
findings. In this paper, we analyse a class of integro-differential
predator-prey models with a highly mobile predator in a heterogeneous
environment, and we reveal the main factors stabilizing such systems. In
particular, we explore an ecologically relevant case of interactions in a
highly eutrophic environment, where the prey carrying capacity can be formally
set to 'infinity'. We investigate two main scenarios: (i) the spatial gradient
of the growth rate is due to abiotic factors only, and (ii) the local growth
rate depends on the global density distribution across the environment (e.g.
due to non-local self-shading). For an arbitrary spatial gradient of the prey
growth rate, we analytically investigate the possibility of the predator-prey
equilibrium in such systems and we explore the conditions of stability of this
equilibrium. In particular, we demonstrate that for a Holling type I (linear)
functional response, the predator can stabilize the system at low prey density
even for an 'unlimited' carrying capacity. We conclude that the interplay
between spatial heterogeneity in the prey growth and fast displacement of the
predator across the habitat works as an efficient stabilizing mechanism.


  Noise caused by fluctuations at the molecular level is a fundamental part of
intracellular processes. While the response of biological systems to noise has
been studied extensively, there has been limited understanding of how to
exploit it to induce a desired cell state. Here we present a scalable,
quantitative method based on the Freidlin-Wentzell action to predict and
control noise-induced switching between different states in genetic networks
that, conveniently, can also control transitions between stable states in the
absence of noise. We apply this methodology to models of cell differentiation
and show how predicted manipulations of tunable factors can induce lineage
changes, and further utilize it to identify new candidate strategies for cancer
therapy in a cell death pathway model. This framework offers a systems approach
to identifying the key factors for rationally manipulating biophysical
dynamics, and should also find use in controlling other classes of noisy
complex networks.


  Protein structure prediction remains to be an open problem in bioinformatics.
There are two main categories of methods for protein structure prediction: Free
Modeling (FM) and Template Based Modeling (TBM). Protein threading, belonging
to the category of template based modeling, identifies the most likely fold
with the target by making a sequence-structure alignment between target protein
and template protein. Though protein threading has been shown to more be
successful for protein structure prediction, it performs poorly for remote
homology detection.


  We present preliminary results that quantify network robustness and fragility
of Ewing sarcoma (ES), a rare pediatric bone cancer that often exhibits de novo
or acquired drug resistance. By identifying novel proteins or pathways
susceptible to drug targeting, this formalized approach promises to improve
preclinical drug development and may lead to better treatment outcomes. Toward
that end, our network modeling focused upon the IGF-1R-PI3K-Akt-mTOR pathway,
which is of proven importance in ES. The clinical response and proteomic
networks of drug-sensitive parental cell lines and their drug-resistant
counterparts were assessed using two small molecule inhibitors for IGF-1R
(OSI-906 and NVP-ADW-742) and an mTOR inhibitor (mTORi), MK8669, such that
protein-to-protein expression networks could be generated for each group. For
the first time, mathematical modeling proves that drug resistant ES samples
exhibit higher degrees of overall network robustness (e.g., the ability of a
system to withstand random perturbations to its network configuration) to that
of their untreated or short-term (72-hour) treated samples. This was done by
leveraging previous work, which suggests that Ricci curvature, a key geometric
feature of a given network, is positively correlated to increased network
robustness. More importantly, given that Ricci curvature is a local property of
the system, it is capable of resolving pathway fragility. In this note, we
offer some encouraging yet limited insights in terms of system-level robustness
of ES and lay the foundation for scope of future work in which a complete study
will be conducted.


  We show how the standard field theoretical language based on creation and
annihilation operators may be used for a straightforward derivation of closed
master equations describing the population dynamics of multivariate stochastic
epidemic models. In order to do that, we introduce an SIR-inspired stochastic
model for hepatitis C virus epidemic, from which we obtain the time evolution
of the mean number of susceptible, infected, recovered and chronically infected
individuals in a population whose total size is allowed to change.


  We study the oscillatory behaviour of a gene regulatory network with
interlinked positive and negative feedback loop. Frequency and amplitude are
two important properties of oscillation. Studied network produces two different
modes of oscillation. In one mode (mode 1) frequency remains constant over a
wide range amplitude and in other mode (mode 2) the amplitude of oscillation
remains constant over a wide range of frequency. Our study reproduces both
features of oscillations in a single gene regulatory network and show that the
negative plus positive feedback loops in gene regulatory network offer
additional advantage. We identified the key parameters/variables responsible
for different modes of oscillation. The network is flexible in switching
between different modes by choosing appropriately the required
parameters/variables.


  The decline in extent of wild pollinators in recent years has been partly
associated with changing farm practices and in particular with increasing
pesticide use. In this paper we combine ecological modelling with economic
analysis of a single farm output under the as- sumption that both pollination
and pest control are essential inputs. We show that the drive to increase farm
output can lead to a local decline in the wild bee population. Commercial bees
are often considered an alternative to wild pollinators, but we show that their
intro- duction can lead to further decline and finally local extinction of wild
bees. The transitions between different outcomes are characterised by threshold
behaviour and are potentially difficult to predict and detect in advance. Small
changes in economic parameters (input prices) and ecological parameters (wild
bees carrying capacity and effect of pesticides on bees) can move the
economic-ecological system beyond the extinction threshold. We also show that
increasing the pesticide price or decreasing the commercial bee price might
lead to re-establishment of wild bees following their local extinction. Thus,
we demonstrate the importance of combining ecological modelling with economics
to study the provision of ecosystem services and to inform sustainable
management of ecosystem service providers.


  In evolutionary game theory, an important measure of a mutant trait
(strategy) is its ability to invade and take over an otherwise-monomorphic
population. Typically, one quantifies the success of a mutant strategy via the
probability that a randomly occurring mutant will fixate in the population.
However, in a structured population, this fixation probability may depend on
where the mutant arises. Moreover, the fixation probability is just one
quantity by which one can measure the success of a mutant; fixation time, for
instance, is another. We define a notion of homogeneity for evolutionary games
that captures what it means for two single-mutant states, i.e. two
configurations of a single mutant in an otherwise-monomorphic population, to be
"evolutionarily equivalent" in the sense that all measures of evolutionary
success are the same for both configurations. Using asymmetric games, we argue
that the term "homogeneous" should apply to the evolutionary process as a whole
rather than to just the population structure. For evolutionary matrix games in
graph-structured populations, we give precise conditions under which the
resulting process is homogeneous. Finally, we show that asymmetric matrix games
can be reduced to symmetric games if the population structure possesses a
sufficient degree of symmetry.