  SARS-CoV is believed to originate from civets and was thought to have been
eliminated as a threat after the 2003 outbreak. Here, we show that human
SARS-CoV (huSARS-CoV) originated directly from bats, rather than civets, by a
cross-species jump in 1991, and formed a human-adapted strain in 1998. Since
then huSARS-CoV has evolved further into highly virulent strains with genotype
T and a 29-nt deletion mutation, and weakly virulent strains with genotype C
but without the 29-nt deletion. The former can cause pneumonia in humans and
could be the major causative pathogen of the SARS outbreak, whereas the latter
might not cause pneumonia in humans, but evolved the ability to co-utilize
civet ACE2 as an entry receptor, leading to interspecies transmission between
humans and civets. Three crucial time points - 1991, for the cross-species jump
from bats to humans; 1998, for the formation of the human-adapted SARS-CoV; and
2003, when there was an outbreak of SARS in humans - were found to associate
with anomalously low annual precipitation and high temperatures in Guangdong.
Anti-SARS-CoV sero-positivity was detected in 20% of all the samples tested
from Guangzhou children who were born after 2005, suggesting that weakly
virulent huSARS-CoVs might still exist in humans. These existing but undetected
SARS-CoVs have a large potential to evolve into highly virulent strains when
favorable climate conditions occur, highlighting a potential risk for the
reemergence of SARS.


  Following recent discoveries of colocalization of downstream-regulating genes
in living cells, the impact of the spatial distance between such genes on the
kinetics of gene product formation is increasingly recognized. We here show
from analytical and numerical analysis that the distance between a
transcription factor (TF) gene and its target gene drastically affects the
speed and reliability of transcriptional regulation in bacterial cells. For an
explicit model system we develop a general theory for the interactions between
a TF and a transcription unit. The observed variations in regulation efficiency
are linked to the magnitude of the variation of the TF concentration peaks as a
function of the binding site distance from the signal source. Our results
support the role of rapid binding site search for gene colocalization and
emphasize the role of local concentration differences.


  Dispersal is a key ecological process, that enables local populations to form
spatially extended systems called metapopulations. In the present study, we
investigate how dispersal affects the linear stability of a general
single-species metapopulation model. We discuss both the influence of local
within-patch dynamics and the effects of various dispersal behaviors on
stability. We find that positive density-dependent dispersal and positive
density-dependent settlement are destabilizing dispersal behaviors while
negative density-dependent dispersal and negative density-dependent settlement
are stabilizing. It is also shown that dispersal has a stabilizing impact on
heterogeneous metapopulations that correlates positively with the number of
patches and the connectance of metapopulation networks.


  Nowadays, chaos theory and nonlinear dynamics lack research focuses. Here we
mention a few major open problems: 1. an effective description of chaos and
turbulence, 2. rough dependence on initial data, 3. arrow of time, 4. the
paradox of enrichment, 5. the paradox of pesticides, 6. the paradox of
plankton.


  The synthetic construction of intracellular circuits is frequently hindered
by a poor knowledge of appropriate kinetics and precise rate parameters. Here,
we use generalized modeling (GM) to study the dynamical behavior of topological
models of a family of hybrid metabolic-genetic circuits known as
"metabolators." Under mild assumptions on the kinetics, we use GM to
analytically prove that all explicit kinetic models which are topologically
analogous to one such circuit, the "core metabolator," cannot undergo Hopf
bifurcations. Then, we examine more detailed models of the metabolator.
Inspired by the experimental observation of a Hopf bifurcation in a
synthetically constructed circuit related to the core metabolator, we apply GM
to identify the critical components of the synthetically constructed
metabolator which must be reintroduced in order to recover the Hopf
bifurcation. Next, we study the dynamics of a re-wired version of the core
metabolator, dubbed the "reverse" metabolator, and show that it exhibits a
substantially richer set of dynamical behaviors, including both local and
global oscillations. Prompted by the observation of relaxation oscillations in
the reverse metabolator, we study the role that a separation of genetic and
metabolic time scales may play in its dynamics, and find that widely separated
time scales promote stability in the circuit. Our results illustrate a generic
pipeline for vetting the potential success of a potential circuit design,
simply by studying the dynamics of the corresponding generalized model.


  In this work we propose a model for gene expression based on the theory of
random dynamical systems (RDS) and show that it has a "modularity property" in
the following sense: given any collection of genes that are linked in a
transcriptional network, if each of them is individually described by a certain
class of RDS then there is a natural, and essentially unique, prescription for
coupling them together, respecting the network topology, in such a way that the
collective system formed by all genes is a RDS as well. Moreover, the class of
RDS used to describe the individual genes is flexible enough to account for a
wide range of stochastic behaviors within the realm of stationary processes.


  We model a general, hierarchically organized tissue by a multi compartment
approach, allowing any number of mutations within a cell. We derive closed
solutions for the deterministic clonal dynamics and the reproductive capacity
of single clones. Our results hold for the average dynamics in a hierarchical
tissue characterized by an arbitrary combination of proliferation parameters.


  I compute exact partition function zeros of the Wako-Saito-Mu\~noz-Eaton
model for various secondary structural elements and for two proteins, 1BBL and
1I6C, using both analytic and numerical methods. Two-state and barrierless
downhill folding transitions can be distinguished by a gap in the distribution
of zeros at the positive real axis.


  The estimation of mutation probabilities and relative fitnesses in
fluctuation analysis is based on the unrealistic hypothesis that the
single-cell times to division are exponentially distributed. Using the
classical Luria-Delbr\"{u}ck distribution outside its modelling hypotheses
induces an important bias on the estimation of the relative fitness. The model
is extended here to any division time distribution. Mutant counts follow a
generalization of the Luria-Delbr\"{u}ck distribution, which depends on the
mean number of mutations, the relative fitness of normal cells compared to
mutants, and the division time distribution of mutant cells. Empirical
probability generating function techniques yield precise estimates both of the
mean number of mutations and the relative fitness of normal cells compared to
mutants. In the case where no information is available on the division time
distribution, it is shown that the estimation procedure using constant division
times yields more reliable results. Numerical results both on observed and
simulated data are reported.


  Next-generation sequencing technologies generate millions of short sequence
reads, which are usually aligned to a reference genome. In many applications,
the key information required for downstream analysis is the number of reads
mapping to each genomic feature, for example to each exon or each gene. The
process of counting reads is called read summarization. Read summarization is
required for a great variety of genomic analyses but has so far received
relatively little attention in the literature. We present featureCounts, a read
summarization program suitable for counting reads generated from either RNA or
genomic DNA sequencing experiments. featureCounts implements highly efficient
chromosome hashing and feature blocking techniques. It is considerably faster
than existing methods (by an order of magnitude for gene-level summarization)
and requires far less computer memory. It works with either single or
paired-end reads and provides a wide range of options appropriate for different
sequencing applications. featureCounts is available under GNU General Public
License as part of the Subread (http://subread.sourceforge.net) or Rsubread
(http://www.bioconductor.org) software packages.


  The evolutionary dynamics of human Influenza A virus presents a challenging
theoretical problem. An extremely high mutation rate allows the virus to
escape, at each epidemic season, the host immune protection elicited by
previous infections. At the same time, at each given epidemic season a single
quasi-species, that is a set of closely related strains, is observed. A
non-trivial relation between the genetic (i.e., at the sequence level) and the
antigenic (i.e., related to the host immune response) distances can shed light
into this puzzle. In this paper we introduce a model in which, in accordance
with experimental observations, a simple interaction rule based on spatial
correlations among point mutations dynamically defines an immunity space in the
space of sequences. We investigate the static and dynamic structure of this
space and we discuss how it affects the dynamics of the virus-host interaction.
Interestingly we observe a staggered time structure in the virus evolution as
in the real Influenza evolutionary dynamics.


  Estimating population models from uncertain observations is an important
problem in ecology. Perretti et al. observed that standard Bayesian state-space
solutions to this problem may provide biased parameter estimates when the
underlying dynamics are chaotic. Consequently, forecasts based on these
estimates showed poor predictive accuracy compared to simple "model-free"
methods, which lead Perretti et al. to conclude that "Model-free forecasting
outperforms the correct mechanistic model for simulated and experimental data".
However, a simple modification of the statistical methods also suffices to
remove the bias and reverse their results.


  In the recent paper (2011) Faas and co-workers claimed that calmodulin could
directly detect Ca2+ signals by acting as extremely fast calcium buffer.
However the results and their interpretation raise serious doubts about this
conclusion.


  Understanding how virus capsids assemble around their nucleic acid (NA)
genomes could promote efforts to block viral propagation or to reengineer
capsids for gene therapy applications. We develop a coarse-grained model of
capsid proteins and NAs with which we investigate assembly dynamics and
thermodynamics. In contrast to recent theoretical models, we find that capsids
spontaneously `overcharge'; i.e., the negative charge of the NA exceeds the
positive charge on capsid. When applied to specific viruses, the optimal NA
lengths closely correspond to the natural genome lengths. Calculations based on
linear polyelectrolytes rather than base-paired NAs underpredict the optimal
length, demonstrating the importance of NA structure to capsid assembly. These
results suggest that electrostatics, excluded volume, and NA tertiary structure
are sufficient to predict assembly thermodynamics and that the ability of
viruses to selectively encapsidate their genomic NAs can be explained, at least
in part, on a thermodynamic basis.


  The influence of hydroxyectoine on the properties of the aqueous solution in
presence of DPPC lipid bilayers is studied via semi-isotropic constant pressure
(NPT) Molecular Dynamics simulations. We investigate the solvent-co-solute
behavior in terms of Kirkwood-Buff integrals as well as hydrogen bond life
times for an increasing hydroxyectoine concentration up to 0.148 mol/L. The
observed preferential exclusion mechanism identifies hydroxyectoine as a
kosmotropic osmolyte. Our findings in regards to the DPPC lipid bilayer
indicate an increase of the surface pressure as well as the solvent accessible
surface area in presence of higher hydroxyectoine concentrations. The results
are in agreement to the outcome of recent experiments. With this study, we are
able to validate the visibility of co-solute-solute-solvent effects for low and
physiologically relevant osmolyte concentrations.


  Many bioinformatics problems, such as sequence alignment, gene prediction,
phylogenetic tree estimation and RNA secondary structure prediction, are often
affected by the "uncertainty" of a solution; that is, the probability of the
solution is extremely small. This situation arises for estimation problems on
high-dimensional discrete spaces in which the number of possible discrete
solutions is immense. In the analysis of biological data or the development of
prediction algorithms, this uncertainty should be handled carefully and
appropriately. In this review, I will explain several methods to combat this
uncertainty, presenting a number of examples in bioinformatics. The methods
include (i) avoiding point estimation, (ii) maximum expected accuracy (MEA)
estimations, and (iii) several strategies to design a pipeline involving
several prediction methods. I believe that the basic concepts and ideas
described in this review will be generally useful for estimation problems in
various areas of bioinformatics.


  We have developed a novel method to rapidly obtain homologous genomic data
for phylogenetics directly from next-generation sequencing reads without the
use of a reference genome. This software, called SISRS, avoids the time
consuming steps of de novo whole genome assembly, genome-genome alignment, and
annotation. For simulations SISRS is able to identify large numbers of loci
containing variable sites with phylogenetic signal. For genomic data from apes,
SISRS identified thousands of variable sites, from which we produced an
accurate phylogeny. Finally, we used SISRS to identify phylogenetic markers
that we used to estimate the phylogeny of placental mammals. We recovered
phylogenies from multiple datasets that were consistent with previous
conflicting estimates of the relationships among mammals. SISRS is open source
and freely available at https://github.com/rachelss/SISRS.


  In the present study, the origin of the anti-quorum sensing (QS) activities
of several members of a recently synthesized and in vitro tested class of
lactone and thiolactone based inhibitors were computationally investigated.
Docking and molecular dynamic (MD) simulations and binding free energy
calculations were carried out to reveal the exact binding and inhibitory
profiles of these compounds. The higher in vitro activity of the lactone series
relative to their thiolactone isosteres was verified based on estimating the
binding energies, the docking scores and monitoring the stability of the
complexes produced in the MD simulations. The strong electrostatic contribution
to the binding energies may be responsible for the higher inhibitory activity
of the lactone with respect to the thiolactone series. The results of this
study help to understand the anti-QS properties of lactone-based inhibitors and
provide important information that may assist in the synthesis of novel QS
inhibitors.


  Motivation: The construction of statistics for summarizing posterior samples
returned by a Bayesian phylogenetic study has so far been hindered by the poor
geometric insights available into the space of phylogenetic trees, and ad hoc
methods such as the derivation of a consensus tree makeup for the
ill-definition of the usual concepts of posterior mean, while bootstrap methods
mitigate the absence of a sound concept of variance. Yielding satisfactory
results with sufficiently concentrated posterior distributions, such methods
fall short of providing a faithful summary of posterior distributions if the
data do not offer compelling evidence for a single topology.
  Results: Building upon previous work of Billera et al., summary statistics
such as sample mean, median and variance are defined as the geometric median,
Fr\'echet mean and variance, respectively. Their computation is enabled by
recently published works, and embeds an algorithm for computing shortest paths
in the space of trees. Studying the phylogeny of a set of plants, where several
tree topologies occur in the posterior sample, the posterior mean balances
correctly the contributions from the different topologies, where a consensus
tree would be biased. Comparisons of the posterior mean, median and consensus
trees with the ground truth using simulated data also reveals the benefits of a
sound averaging method when reconstructing phylogenetic trees.


  A protocol for the identification of ancestry informative markers (AIMs) from
genome-wide single nucleotide polymorphism (SNP) data is proposed. The protocol
consists of three main steps: (a) identification of potential positive
selection regions via Fst extremity measurement, (b) SNP screening via
two-stage attribute selection and (c) classification model construction using a
naive Bayes classifier. The two-stage attribute selection is composed of a
newly developed round robin symmetrical uncertainty ranking technique and a
wrapper embedded with a naive Bayes classifier. The protocol has been applied
to the HapMap Phase II data. Two AIM panels, which consist of 10 and 16 SNPs
that lead to complete classification between CEU, CHB, JPT and YRI populations,
are identified. Moreover, the panels are at least four times smaller than those
reported in previous studies. The results suggest that the protocol could be
useful in a scenario involving a larger number of populations.


  The karyotype ontology describes the human chromosome complement as
determined cytogenetically, and is designed as an initial step toward the goal
of replacing the current system which is based on semantically meaningful
strings. This ontology uses a novel, semi-programmatic methodology based around
the tawny library to construct many classes rapidly. Here, we describe our use
case, methodology and the event-based approach that we use to represent
karyotypes.
  The ontology is available at http://www.purl.org/ontolink/karyotype/. The
clojure code is available at http://code.google.com/p/karyotype-clj/.


  We introduce a method to estimate the complexity function of symbolic
dynamical systems from a finite sequence of symbols. We test such complexity
estimator on several symbolic dynamical systems whose complexity functions are
known exactly. We use this technique to estimate the complexity function for
genomes of several organisms under the assumption that a genome is a sequence
produced by a (unknown) dynamical system. We show that the genome of several
organisms share the property that their complexity functions behaves
exponentially for words of small length $\ell$ ($0\leq \ell \leq 10$) and
linearly for word lengths in the range $11 \leq \ell \leq 50$. It is also found
that the species which are phylogenetically close each other have similar
complexity functions calculated from a sample of their corresponding coding
regions.


  Robustness and sensitivity of responses generated by cell signaling networks
has been associated with survival and evolvability of organisms. However,
existing methods analyzing robustness and sensitivity of signaling networks
ignore the experimentally observed cell-to-cell variations of protein
abundances and cell functions or contain ad hoc assumptions. We propose and
apply a data driven Maximum Entropy (MaxEnt) based method to quantify
robustness and sensitivity of Escherichia coli (E. coli) chemotaxis signaling
network. Our analysis correctly rank orders different models of E. coli
chemotaxis based on their robustness and suggests that parameters regulating
cell signaling are evolutionary selected to vary in individual cells according
to their abilities to perturb cell functions. Furthermore, predictions from our
approach regarding distribution of protein abundances and properties of
chemotactic responses in individual cells based on cell population averaged
data are in excellent agreement with their experimental counterparts. Our
approach is general and can be used to evaluate robustness as well as generate
predictions of single cell properties based on population averaged experimental
data in a wide range of cell signaling systems.


  Nucleosome core particle is a dynamic structure -- DNA may transiently peel
off the histone octamer surface due to thermal fluctuations or the action of
chromatin remodeling enzymes. Partial DNA unwrapping enables easier access of
DNA-binding factors to their target sites and thus may provide a dominant
pathway for effecting rapid and robust access to DNA packaged into chromatin.
Indeed, a recent high-resolution map of distances between neighboring
nucleosome dyads in \emph{S.cerevisiae} shows that at least 38.7\% of all
nucleosomes are partially unwrapped. The extent of unwrapping follows a
stereotypical pattern in the vicinity of genes, reminiscent of the canonical
pattern of nucleosome occupancy in which nucleosomes are depleted over
promoters and well-positioned over coding regions. To explain these
observations, we developed a biophysical model which employs a 10-11 base pair
periodic nucleosome energy profile. The profile, based on the pattern of
histone-DNA contacts in nucleosome crystal structures and the idea of linker
length discretization, accounts for both nucleosome unwrapping and higher-order
chromatin structure. Our model reproduces the observed genome-wide distribution
of inter-dyad distances, and accounts for patterns of nucleosome occupancy and
unwrapping around coding regions. At the same time, our approach explains
\emph{in vitro} measurements of accessibility of nucleosome-covered binding
sites, and of nucleosome-induced cooperativity between DNA-binding factors. We
are able to rule out several alternative scenarios of nucleosome unwrapping as
inconsistent with the genomic data.


  Interactions between individual members of the B-cell lymphoma 2 (Bcl-2)
family of proteins form a regulatory network governing mitochondrial outer
membrane permeabilization (MOMP). Bcl-2 family initiated MOMP causes release of
the inter-membrane pro-apoptotic proteins to cytosol and creates a cytosolic
environment suitable for the executionary phase of apoptosis. Using
mathematical modeling and computational simulation, we have analyzed the
response of this regulatory network caused by the up-/downregulation of each
individual member of the Bcl-2 family. As a result, a non-linear
stimulus-response emerged, the characteristics of which are associated with
bistability and switch-like behavior. Using the principal component analysis
(PCA) we have shown that the Bcl-2 family classifies the random combinations of
inputs into two distinct classes, and responds to these by one of the two
qualitatively distinct outputs. As we showed, the emergence of this behavior
requires specific organization of the interactions between particular Bcl-2
proteins.


  Genetic association studies, in particular the genome-wide association study
design, have provided a wealth of novel insights into the aetiology of a wide
range of human diseases and traits. The next challenge consists of
understanding the molecular basis of these associations. The integration of
multiple association datasets, including gene expression datasets, can
contribute to this goal. We have developed a novel statistical methodology to
assess whether two association signals are consistent with a shared causal
variant. An application is the integration of disease scans with expression
quantitative trait locus (eQTL) studies, but any pair of GWAS datasets can be
integrated in this framework. We demonstrate the value of the approach by
reanalysing a gene expression dataset in 966 liver samples with a published
meta-analysis of lipid traits including >100, 000 individuals of European
ancestry. Combining all lipid biomarkers, our reanalysis supported 29 out of 38
reported colocalisation results with eQTLs and identified 14 new colocalisation
results, highlighting the value of a formal statistical test. In two cases of
reported eQTL-lipid pairs (IFT172, TBKBP1) for which our analysis suggests that
the eQTL pattern is not consistent with the lipid association, we identify
alternative colocalisation results with GCKR and KPNB1, indicating that these
genes are more likely to be causal in these genomic intervals. A key feature of
the method is the ability to derive the output statistics from single SNP
summary statistics, hence making it possible to perform systematic
meta-analysis type comparisons across multiple GWAS datasets
(http://coloc.cs.ucl.ac.uk/coloc/). Our methodology provides information about
candidate causal genes in associated intervals and has direct implications for
the understanding of complex diseases and the design of drugs to target disease
pathways.


  The base pair fluctuations and helix untwisting are examined for a circular
molecule. A realistic mesoscopic model including twisting degrees of freedom
and bending of the molecular axis is proposed. The computational method, based
on path integral techniques, simulates a distribution of topoisomers with
various twist numbers and finds the energetically most favorable molecular
conformation as a function of temperature. The method can predict helical
repeat, openings loci and bubble sizes for specific sequences in a broad
temperature range. Some results are presented for a short DNA circle recently
identified in mammalian cells.


  Phytophthora infestans, the cause of potato late blight, is infamous for
having triggered the Irish Great Famine in the 1840s. Until the late 1970s, P.
infestans diversity outside of its Mexican center of origin was low, and one
scenario held that a single strain, US-1, had dominated the global population
for 150 years; this was later challenged based on DNA analysis of historical
herbarium specimens. We have compared the genomes of 11 herbarium and 15 modern
strains. We conclude that the nineteenth century epidemic was caused by a
unique genotype, HERB-1, that persisted for over 50 years. HERB-1 is distinct
from all examined modern strains, but it is a close relative of US-1, which
replaced it outside of Mexico in the twentieth century. We propose that HERB-1
and US-1 emerged from a metapopulation that was established in the early 1800s
outside of the species' center of diversity.


  In his highly influential view of evolution, G. G. Simpson hypothesized that
clades of species evolve in adaptive zones, defined as collections of niches
occupied by species with similar traits and patterns of habitat use. Simpson
hypothesized that species enter new adaptive zones in one of three ways:
extinction of competitor species, dispersal to a new geographic region, or the
evolution of a key trait that allows species to exploit resources in a new way.
However, direct tests of Simpson's hypotheses for the entry into new adaptive
zones remain elusive. Here we evaluate the fit of a Simpsonian model of jumps
between adaptive zones to phylogenetic comparative data. We use a novel
statistical approach to show that anoles, a well-studied adaptive radiation of
Caribbean lizards, have evolved by a series of evolutionary jumps in trait
evolution. Furthermore, as Simpson predicted, trait axes strongly tied to
habitat specialization show jumps that correspond with the evolution of key
traits and/or dispersal between islands in the Greater Antilles. We conclude
that jumps are commonly associated with major adaptive shifts in the
evolutionary radiation of anoles.


  Mathematical methods provide useful framework for the analysis and design of
complex systems. In newer contexts such as biology, however, there is a need to
both adapt existing methods as well as to develop new ones. Using a combination
of analytical and computational approaches, we adapt and develop the method of
describing functions to represent the input-output responses of biomolecular
signalling systems. We approximate representative systems exhibiting various
saturating and hysteretic dynamics in a way that is better than the standard
linearization. Further, we develop analytical upper bounds for the
computational error estimates. Finally, we use these error estimates to augment
the limit cycle analysis with a simple and quick way to bound the predicted
oscillation amplitude. These results provide system approximations that can add
more insight into the local behaviour of these systems than standard
linearization, compute responses to other periodic inputs, and to analyze limit
cycles.


  Kin selection theory is a kind of causal analysis. The initial form of kin
selection ascribed cause to costs, benefits, and genetic relatedness. The
theory then slowly developed a deeper and more sophisticated approach to
partitioning the causes of social evolution. Controversy followed because
causal analysis inevitably attracts opposing views. It is always possible to
separate total effects into different component causes. Alternative causal
schemes emphasize different aspects of a problem, reflecting the distinct
goals, interests, and biases of different perspectives. For example, group
selection is a particular causal scheme with certain advantages and significant
limitations. Ultimately, to use kin selection theory to analyze natural
patterns and to understand the history of debates over different approaches,
one must follow the underlying history of causal analysis. This article
describes the history of kin selection theory, with emphasis on how the causal
perspective improved through the study of key patterns of natural history, such
as dispersal and sex ratio, and through a unified approach to demographic and
social processes. Independent historical developments in the multivariate
analysis of quantitative traits merged with the causal analysis of social
evolution by kin selection.


  Ecological trade-offs between species are often invoked to explain species
coexistence in ecological communities. However, few mathematical models have
been proposed for which coexistence conditions can be characterized explicitly
in terms of a trade-off. Here we present a model of a plant community which
allows such a characterization. In the model plant species compete for sites
where each site has a fixed stress condition. Species differ both in stress
tolerance and competitive ability. Stress tolerance is quantified as the
fraction of sites with stress conditions low enough to allow establishment.
Competitive ability is quantified as the propensity to win the competition for
empty sites. We derive the deterministic, discrete-time dynamical system for
the species abundances. We prove the conditions under which plant species can
coexist in a stable equilibrium. We show that the coexistence conditions can be
characterized graphically, clearly illustrating the trade-off between stress
tolerance and competitive ability. We compare our model with a recently
proposed, continuous-time dynamical system for a tolerance-fecundity trade-off
in plant communities, and we show that this model is a special case of the
continuous-time version of our model.


  The spread of certain diseases can be promoted, in some cases substantially,
by prior infection with another disease. One example is that of HIV, whose
immunosuppressant effects significantly increase the chances of infection with
other pathogens. Such coinfection processes, when combined with nontrivial
structure in the contact networks over which diseases spread, can lead to
complex patterns of epidemiological behavior. Here we consider a mathematical
model of two diseases spreading through a single population, where infection
with one disease is dependent on prior infection with the other. We solve
exactly for the sizes of the outbreaks of both diseases in the limit of large
population size, along with the complete phase diagram of the system. Among
other things, we use our model to demonstrate how diseases can be controlled
not only by reducing the rate of their spread, but also by reducing the spread
of other infections upon which they depend.


  Two novel reduced dimensionality (RD) experiments -(4,3)D-hNCOcaNH and
(4,3)D-hNcoCANH- have been presented here to facilitate the backbone resonance
assignment of proteins both in terms of its accuracy and speed. The experiments
basically represent an improvisation of previously reported HN(C)N experiment
[Panchal et. al., J. Biomol. NMR. (2002), 20 (2), 135-147] and exploit the
simple reduced dimensionality NMR concept [Szyperski et. al. (2002), Proc.
Natl. Acad. Sci. U.S.A. 99(12), 8009-8014] to achieve (a) higher dispersion and
resolution along the co-evolved F1 dimension and (b) rapid identification of
sequential HSQC peaks on its F2(15N)- F3(1H) planes. The current implementation
is based on the fact that the linear combination of 15N and 13CO/13Ca chemical
shifts offers relatively better dispersion and randomness compared to the
individual chemical shifts; thus enables the assignment of crowded HSQC spectra
by resolving the ambiguities generally encountered in HNCN based assignment
protocol because of amide 15N shift degeneracy. Additionally, each of these
experiments enables assignment of backbone 13CO/13Ca resonances as well.
Overall, the reduced dimensionality tailored HNCN experiments presented here
will be of immense value for various structural and functional proteomics
studies by NMR; particularly of intrinsically or partially/unstructured
proteins and medium sized (MW ~12-15 kDa) folded proteins. The experiments,
like any other experiment that yields protein assignments, would be extremely
valuable for protein folding and drug discovery programs as well.


  We develop an analytical approach to the susceptible-infected-susceptible
(SIS) epidemic model that allows us to unravel the true origin of the absence
of an epidemic threshold in heterogeneous networks. We find that a delicate
balance between the number of high degree nodes in the network and the
topological distance between them dictates the existence or absence of such a
threshold. In particular, small-world random networks with a degree
distribution decaying slower than an exponential have a vanishing epidemic
threshold in the thermodynamic limit.


  Semantic similarity measures (SSMs) refer to a set of algorithms used to
quantify the similarity of two or more terms belonging to the same ontology.
Ontology terms may be associated to concepts, for instance in computational
biology gene and proteins are associated with terms of biological ontologies.
Thus, SSMs may be used to quantify the similarity of genes and proteins
starting from the comparison of the associated annotations. SSMs have been
recently used to compare genes and proteins even on a system level scale. More
recently some works have focused on the building and analysis of Semantic
Similarity Networks (SSNs) i.e. weighted networks in which nodes represents
genes or proteins while weighted edges represent the semantic similarity score
among them. SSNs are quasi-complete networks, thus their analysis presents
different challenges that should be addressed. For instance, the need for the
introduction of reliable thresholds for the elimination of meaningless edges
arises. Nevertheless, the use of global thresholding methods may produce the
elimination of meaningful nodes, while the use of local thresholds may
introduce biases. For these aims, we introduce a novel technique, based on
spectral graph considerations and on a mixed global-local focus. The
effectiveness of our technique is demonstrated by using markov clustering for
the extraction of biological modules. We applied clustering to simplified
networks demonstrating a considerable improvements with respect to the original
ones.


  This study reveals that, in strong coulomb coupling regime, bending a
straight and fully overcharged DNA (up to its maximal acceptance by multivalent
counterions) to a circle releases some of the adsorbed (correlated)counterions
but still remains fully overcharged. This phenomenon seems to be inherent to
the minimum energy state of a DNA. By definition, the total electrostatic
potential energy of a macroion-counterion system reaches to its lowest point at
maximal acceptance of overcharging counterions that ensures the most stable
conformation. This intermediate phenomenon of release of cations from DNA
surface due to bending can be taken into account in theoretical modeling of
some ionic concentration dependent physico-chemical aspects of DNA solutions in
strong Coulomb coupling regimes.


  Plasmodium falciparum is the parasite that causes the most severe form of
malaria. Currently, science has been established about its cellular structures,
its metabolic processes, and even the molecular structures of its intrinsic
membrane proteins responsible for transporting water, nutrient, and waste
molecules across the parasite plasma membrane (PPM). I hypothesize that
Plasmodium falciparum has an Achilles' heel that can be attacked with
erythritol, the well-known sweetener that is classified as generally safe. Most
organisms have in their cell membrane two types of water-channel proteins:
aquaporins to maintain hydro-homeostasis across the membrane and
aquaglyceroporins to uptake glycerols etc. In contrast, P. falciparum has only
one type of such proteins---the multi-functional aquaglyceroporin (PfAQP)
expressed in the PPM---to do both jobs. Moreover, the parasite also uses PfAQP
to excrete its metabolic wastes (ammonia included) produced at a very high rate
in the blood stage. This extremely high efficiency of the bug using one protein
for multiple essential tasks makes the parasite fatally vulnerable. Erythritol
in the blood stream can kill the parasite by clogging up its PfAQP channel that
needs to be open for maintaining hydro-homeostasis and for excreting toxic
wastes across the bug's PPM. In vitro tests are to measure the growth/death
rate of P. falciparum in blood with various erythritol concentrations. In vivo
experiments are to administer groups of infected mice with various doses of
erythritol and monitor the parasite growth levels from blood samples drawn from
each group. Clinic trials can be performed to observe the added effects of
administering to patients erythritol along with the known drugs because
erythritol was classified as a safe food ingredient.


  We investigate the effects of local population structure in
reaction-diffusion processes representing a contact process (CP) on
metapopulations represented as complex networks. Considering a model in which
the nodes of a large scale network represent local populations defined in terms
of a homogeneous graph, we show by means of extensive numerical simulations
that the critical properties of the reaction-diffusion system are independent
of the local population structure, even when this one is given by a ordered
linear chain. This independence is confirmed by the perfect matching between
numerical critical exponents and the results from a heterogeneous mean field
theory suited, in principle, to describe situations of local homogeneous
mixing. The analysis of several variations of the reaction-diffusion process
allow to conclude the independence from population structure of the critical
properties of CP-like models on metapopulations, and thus of the universality
of the reaction-diffusion description of this kind of models.


  Using the Helmholtz decomposition of the vector field of folding fluxes in a
two-dimensional space of collective variables, a potential of the driving force
for protein folding is introduced. The potential has two components. One
component is responsible for the source and sink of the folding flows, which
represent, respectively, the unfolded states and the native state of the
protein, and the other, which accounts for the flow vorticity inherently
generated at the periphery of the flow field, is responsible for the
canalization of the flow between the source and sink. The theoretical
consideration is illustrated by calculations for a model $\beta$-hairpin
protein.


  This paper introduces the class of "strongly endotactic networks", a subclass
of the endotactic networks introduced by G. Craciun, F. Nazarov, and C. Pantea.
The main result states that the global attractor conjecture holds for
complex-balanced systems that are strongly endotactic: every trajectory with
positive initial condition converges to the unique positive equilibrium allowed
by conservation laws. This extends a recent result by D. F. Anderson for
systems where the reaction diagram has only one linkage class (connected
component). The results here are proved using differential inclusions, a
setting that includes power-law systems. The key ideas include a perspective on
reaction kinetics in terms of combinatorial geometry of reaction diagrams, a
projection argument that enables analysis of a given system in terms of systems
with lower dimension, and an extension of Birch's theorem, a well-known result
about intersections of affine subspaces with manifolds parameterized by
monomials.


  Models of Stock Recruitment Relationships (SRRs) are often used to predict
fish population dynamics. Commonly used SRRs include the Ricker, Beverton-Holt,
and Cushing functional forms, which differ primarily by the degree of density
dependent effects (compensation). The degree of compensation determines whether
recruitment respectively decreases, saturates, or increases at high levels of
spawning stock biomass. In 1982 J.G. Shepherd united these dynamics into a
single model, where the degree of compensation is determined by a single
parameter, however the difficulty in relating this parameter to biological data
has limited its usefulness. Here we use a generalized modeling framework to
show that the degree of compensation can be related directly to the functional
elasticity of growth, which is a general quantity that measures the change in
recruitment relative to a change in biomass, irrespective of the specific SRR.
We show that the elasticity of growth can be calculated from short-term
fluctuations in fish biomass, is robust to observation error, and can be used
to determine general attributes of the SRR in both continuous time production
models, as well as discrete time age-structured models. This framework may be
particularly useful if fisheries time-series data are limited, and not
conducive to determining functional relationships using traditional methods of
statistical best-fit.


  The identification and validation of gene-gene interactions is a major
challenge in human studies. Here, we explore an approach for studying epistasis
in humans using a Drosophila melanogaster model of neonatal diabetes mellitus.
Expression of mutant preproinsulin, hINSC96Y, in the eye imaginal disc mimics
the human disease activating conserved cell stress response pathways leading to
cell death and reduction in eye area. Dominant-acting variants in wild-derived
inbred lines from the Drosophila Genetics Reference Panel produce a continuous,
highly heritable, distribution of eye degeneration phenotypes. A genome-wide
association study (GWAS) in 154 sequenced lines identified 29 candidate SNPs in
16 loci with P < 10-5 including one SNP in an intron of the gene sulfateless
(sfl) which exceeded a conservative genome-wide significance threshold of P =
0.05 level (-log10 P > 7.62). RNAi knock-downs of sfl enhanced the eye
degeneration phenotype in a mutant-hINS-dependent manner. sfl encodes a protein
required for sulfation of the glycosaminoglycan, heparan sulfate. Two
additional genes in the heparan sulfate (HS) biosynthetic pathway (tout velu,
ttv and brother of tout velu, botv) also modified the eye phenotype, suggesting
a link between HS-modified proteins and cellular responses to misfolded
proteins. Finally, intronic variants marking the QTL were associated with
decreased sfl expression, a result consistent with that predicted by RNAi
studies. The ability to create a model of human genetic disease in the fly, map
a QTL by GWAS to a specific gene (and noncoding variant), validate its
contribution to disease with available genetic resources, and experimentally
link the variant to a molecular mechanism, demonstrate the many advantages
Drosophila holds in determining the genetic underpinnings of human disease.


  One of the major developments in recent years in the search for missing
heritability of human phenotypes is the adoption of linear mixed-effects models
(LMMs) to estimate heritability due to genetic variants which are not
significantly associated with the phenotype. A variant of the LMM approach has
been adapted to case-control studies and applied to many major diseases by Lee
et al. (2011), successfully accounting for a considerable portion of the
missing heritability. For example, for Crohn's disease their estimated
heritability was 22% compared to 50-60% from family studies. In this letter we
propose to estimate heritability of disease directly by regression of phenotype
similarities on genotype correlations, corrected to account for ascertainment.
We refer to this method as genetic correlation regression (GCR). Using GCR we
estimate the heritability of Crohn's disease at 34% using the same data. We
demonstrate through extensive simulation that our method yields unbiased
heritability estimates, which are consistently higher than LMM estimates.
Moreover, we develop a heuristic correction to LMM estimates, which can be
applied to published LMM results. Applying our heuristic correction increases
the estimated heritability of multiple sclerosis from 30% to 52.6%.


  We characterize different cell states, related to cancer and ageing
phenotypes, by a measure of entropy of network ensembles, integrating gene
expression values and protein interaction networks. The entropy measure
estimates the parameter space available to the network ensemble, that can be
interpreted as the level of plasticity of the system for high entropy values
(the ability to change its internal parameters, e.g. in response to
environmental stimuli), or as a fine tuning of the parameters (that restricts
the range of possible parameter values) in the opposite case. This approach can
be applied at different scales, from whole cell to single biological functions,
by defining appropriate subnetworks based on a priori biological knowledge,
thus allowing a deeper understanding of the cell processes involved. In our
analysis we used specific network features (degree sequence, subnetwork
structure and distance between gene profiles) to obtain informations at
different biological scales, providing a novel point of view for the
integration of experimental transcriptomic data and a priori biological
knowledge, but the entropy measure can also highlight other aspects of the
biological systems studied depending on the constraints introduced in the model
(e.g. community structures).


  Mutualistic communities have an internal structure that makes them resilient
to external per- turbations. Late research has focused on their stability and
the topology of the relations between the different organisms to explain the
reasons of the system robustness. Much less attention has been invested in
analyzing the systems dynamics. The main population models in use are modifi-
cations of the logistic equation with additional terms to account for the
benefits produced by the interspecific interactions. These models have
shortcomings as the so called r - K formulation of logistic equation diverges
under some conditions. In this work, we introduce a model for population
dynamics under mutualism inspired by the logistic equation but avoiding
singularities. The model is mathematically simpler than the widely used type II
models, although it shows similar complexity in terms of fixed points and
stability of the dynamics. Furthermore, each term of our model has a more
direct ecological interpretation, which can facilitate the measurement of the
rates involved in field campaigns. We perform an analytical stability analysis
and numerical simulations to study the model behavior in more general
interaction scenarios including tests of the resilience of its dynamics under
external perturbations. Despite its simplicity, our results indicate that the
model dynamics shows an important richness that can be used to gain further
insights in the dynamics of mutualistic communities.


  We apply formal measures of emergence, self-organization, homeostasis,
autopoiesis and complexity to an aquatic ecosystem; in particular to the
physiochemical component of an Arctic lake. These measures are based on
information theory. Variables with an homogeneous distribution have higher
values of emergence, while variables with a more heterogeneous distribution
have a higher self-organization. Variables with a high complexity reflect a
balance between change (emergence) and regularity/order (self-organization). In
addition, homeostasis values coincide with the variation of the winter and
summer seasons. Autopoiesis values show a higher degree of independence of
biological components over their environment. Our approach shows how the
ecological dynamics can be described in terms of information.


  Myoglobin modulates the binding of diatomic molecules to its heme group via
hydrogen-bonding and steric interactions with neighboring residues, and is an
important benchmark for computational studies of biomolecules. We have
performed calculations on the heme binding site and a significant proportion of
the protein environment (more than 1000 atoms) using linear-scaling density
functional theory and the DFT+U method to correct for self-interaction errors
associated with localized 3d states. We confirm both the hydrogen-bonding
nature of the discrimination effect (3.6 kcal/mol) and assumptions that the
relative strain energy stored in the protein is low (less than 1 kcal/mol). Our
calculations significantly widen the scope for tackling problems in drug design
and enzymology, especially in cases where electron localization, allostery or
long-ranged polarization influence ligand binding and reaction.


  Our aim is to provide an overview of half-saturation constants reported in
literature and to explore their consistency with body size. In many ecological
models, intake of nutrients by plants and consumption of food by animals is
considered to be a hyperbolic function of the nutrient concentration and the
food density, respectively. However, data on the concentration (or density) at
which half of the maximum intake rate is reached are scarce, limiting the
applicability of the computational models. The meta-analysis was conducted on
literature published worldwide. Most studies focused on algae and
invertebrates, whereas some included fish, birds and mammals. The
half-saturation constants obtained were linked to body size using ordinary
regression analysis. The observed trends were compared to those noted in
reviews on other density parameters. Half-saturation constants for different
clades range within one or two orders of magnitude. Although these constants
are inherently variable, exploring allometric relationships across different
taxa helps to improve consistent parameterization of ecological models.


  Suppose we label the vertices of a tree by positive integers. The weight of
an edge is defined by a monotonically increasing function of the absolute value
of the difference of the labels of its endpoints. We define the total cost of
the labeling to be the sum of weight of all the edges.The problem we consider
is that of determining for a given tree G and given a labeling of the leaves of
G the minimum total cost labellings of G. In this paper we present an algorithm
that works for any cost function satisfies the condition of monotony mentioned
above. In a case of the function defined as the absolute value of the
difference of the labels the fast algorithm is presented.


  The emergence of self-replication and information transmission in life's
origin remains unexplained despite extensive research on the topic. A
hypothesis explaining the transition from a simple organic world to a complex
RNA world is offered here based on physical factors in hydrothermal vent
systems. An interdisciplinary approach is taken using techniques from
thermodynamics, fluid dynamics, oceanography, statistical mechanics, and
stochastic processes to examine nucleic acid dynamics and kinetics in a
hydrothermal vent from first principles. Analyses are carried out using both
analytic and computational methods and confirm the plausibility of a reaction
involving the PCR-like assembly of ribonucleotides. The proposal is put into
perspective with established theories on the origin of life and more generally
the onset of order and information transmission in prebiotic systems. A
biomimicry application of this hypothetical process to PCR technology is
suggested and its viability is evaluated in a rigorous logical analysis.
Optimal temperature curves begin to be established using Monte Carlo
simulation, variational calculus, and Fourier analysis. The converse argument
is also made but qualitatively, asserting that the success of such a
modification to PCR would in turn reconfirm the biological theory.


  Searching experiments conducted in different virtual environments over a
gender balanced group of people revealed a gender irrelevant scale-free spread
of searching activity on large spatiotemporal scales. We have suggested and
solved analytically a simple statistical model of the coherent-noise type
describing the exploration-exploitation trade-off in humans ("should I stay or
should I go"). The model exhibits a variety of saltatory behaviours, ranging
from Levy flights occurring under uncertainty to Brownian walks performed by a
treasure hunter confident of the eventual success.


  We explore the dynamical large-deviations of a lattice heteropolymer model of
a protein by means of path sampling of trajectories. We uncover the existence
of non-equilibrium dynamical phase-transitions in ensembles of trajectories
between active and inactive dynamical phases, whose nature depends on
properties of the interaction potential. When the full heterogeneity of
interactions due to the amino-acid sequence is preserved, as in a fully
interacting model or in a heterogeneous version of the G\={o} model where only
native interactions are considered, the transition is between the equilibrium
native state and a highly native but kinetically trapped state. In contrast,
for the homogeneous G\={o} model, where there is a single native energy and the
sequence plays no role, the dynamical transition is a direct consequence of the
static bi-stability between unfolded and native states. In the heterogeneous
case the native-active and native-inactive states, despite their static
similarity, have widely varying dynamical properties, and the transition
between them occurs even in lattice proteins whose sequences are designed to
make them optimal folders.


  What is aging? Mechanistic answers to this question remain elusive despite
decades of research. Here, we propose a mathematical model of cellular aging
based on a model gene interaction network. Our network model is made of only
non-aging components - the biological functions of gene interactions decrease
with a constant mortality rate. Death of a cell occurs in the model when an
essential gene loses all of its interactions to other genes, equivalent to the
deletion of an essential gene. Gene interactions are stochastic based on a
binomial distribution. We show that the defining characteristic of biological
aging, the exponential increase of mortality rate over time, can arise from
this gene network model during the early stage of aging. Hence, we demonstrate
that cellular aging is an emergent property of this model network. Our model
predicts that the rate of aging, defined by the Gompertz coefficient, is
approximately proportional to the average number of active interactions per
gene and that the stochastic heterogeneity of gene interactions is an important
factor in the dynamics of the aging process. This theoretic framework offers a
mechanistic foundation for the pleiotropic nature of aging and can provide
insights on cellular aging.


  We investigate performance and optimization of the "digital" bioanalytical
response. Specifically, we consider the recently introduced approach of a
partial input conversion into inactive compounds, resulting in the "branch
point effect" similar to that encountered in biological systems. This
corresponds to an "intensity filter," which can yield a binary-type
sigmoid-response output signal of interest in information and signal processing
and in biosensing applications. We define measures for optimizing the response
for information processing applications based on the kinetic modeling of the
enzymatic reactions involved, and apply the developed approach to the recently
published data for glucose detection.


  Many biochemical and industrial applications involve complicated networks of
simultaneously occurring chemical reactions. Under the assumption of mass
action kinetics, the dynamics of these chemical reaction networks are governed
by systems of polynomial ordinary differential equations. The steady states of
these mass action systems have been analysed via a variety of techniques,
including elementary flux mode analysis, algebraic techniques (e.g. Groebner
bases), and deficiency theory. In this paper, we present a novel method for
characterizing the steady states of mass action systems. Our method explicitly
links a network's capacity to permit a particular class of steady states,
called toric steady states, to topological properties of a related network
called a translated chemical reaction network. These networks share their
reaction stoichiometries with their source network but are permitted to have
different complex stoichiometries and different network topologies. We apply
the results to examples drawn from the biochemical literature.


  During embryonic development tissue morphogenesis and signaling are tightly
coupled. It is therefore important to simulate both tissue morphogenesis and
signaling simultaneously in in silico models of developmental processes. The
resolution of the processes depends on the questions of interest. As part of
this chapter we will introduce different descriptions of tissue morphogenesis.
In the most simple approximation tissue is a continuous domain and tissue
expansion is described according to a pre-defined function of time (and
possibly space). In a slightly more advanced version the expansion speed and
direction of the tissue may depend on a signaling variable that evolves on the
domain. Both versions will be referred to as 'prescribed growth'. Alternatively
tissue can be regarded as incompressible fluid and can be described with
Navier-Stokes equations. Local cell expansion, proliferation, and death are
then incorporated by a source term. In other applications the cell boundaries
may be important and cell-based models must be introduced. Finally, cells may
move within the tissue, a process best described by agent-based models.


  We develop a framework for optimizing a novel approach to extending the
linear range of bioanalytical systems and biosensors by utilizing two enzymes
with different kinetic responses to the input chemical as their substrate. Data
for the flow-injection amperometric system devised for detection of lysine
based on the function of L-Lysine-alpha-Oxidase and Lysine-2-monooxygenase are
analyzed. Lysine is a homotropic substrate for the latter enzyme. We elucidate
the mechanism for extending the linear response range and develop optimization
techniques for future applications of such systems.


  Computing based on biochemical processes is a newest rapidly developing field
of unconventional information and signal processing. In this paper we present
results of our research in the field of biochemical computing and summarize the
obtained numerical and experimental data for implementations of the standard
two-input OR and AND gates with double-sigmoid shape of the output signal. This
form of response was obtained as a function of the two inputs in each of the
realized biochemical systems. The enzymatic gate processes in the first system
were activated with two chemical inputs and resulted in optically detected
chromogen oxidation, which happens when either one or both of the inputs are
present. In this case, the biochemical system is functioning as the OR gate. We
demonstrate that the addition of a "filtering" biocatalytic process leads to a
considerable reduction of the noise transmission factor and the resulting gate
response has sigmoid shape in both inputs. The second system was developed for
functioning as an AND gate, where the output signal was activated only by a
simultaneous action of two enzymatic biomarkers. This response can be used as
an indicator of liver damage, but only if both of these of the inputs are
present at their elevated, pathophysiological values of concentrations. A
kinetic numerical model was developed and used to estimate the range of
parameters for which the experimentally realized logic gate is close to
optimal. We also analyzed the system to evaluate its noise-handling properties.


  We consider the Moran model in continuous time with two types, mutation, and
selection. We concentrate on the ancestral line and its stationary type
distribution. Building on work by Fearnhead (J. Appl. Prob. 39 (2002), 38-54)
and Taylor (Electron. J. Probab. 12 (2007), 808-847), we characterise this
distribution via the fixation probability of the offspring of all individuals
of favourable type (regardless of the offspring's types). We concentrate on a
finite population and stay with the resulting discrete setting all the way
through. This way, we extend previous results and gain new insight into the
underlying particle picture.


  The competition between two ecologically similar species that use the same
resources and differ from each other only in the type of spatial motion they
undergo is studied. The latter is assumed to be described either by Brownian
motion or L\'evy flights. Competition is taken into account by assuming that
individuals reproduce in a density-dependent fashion. It is observed that no
influence of the type of motion occurs when the two species are in a well-mixed
unstructured state. However, as soon as the species develop spatial clustering,
the one forming more concentrated clusters gets a competitive advantage and
eliminates the other. Similar competitive advantage would occur between walkers
of the same type but with different diffusivities if this leads also to
different clustering. Coexistence of both species is also possible under
certain conditions.


  Statistical properties of the site frequency spectrum associated with
Lambda-coalescents are our objects of study. In particular, we derive
recursions for the expected value, variance, and covariance of the spectrum,
extending earlier results of Fu (1995) for the classical Kingman coalescent.
Estimating coalescent parameters introduced by certain Lambda-coalescents for
datasets too large for full likelihood methods is our focus. The recursions for
the expected values we obtain can be used to find the parameter values which
give the best fit to the observed frequency spectrum. The expected values are
also used to approximate the probability a (derived) mutation arises on a
branch subtending a given number of leaves (DNA sequences), allowing us to
apply a pseudo-likelihood inference to estimate coalescence parameters
associated with certain subclasses of Lambda coalescents. The properties of the
pseudo-likelihood approach are investigated on simulated as well as real mtDNA
datasets for the high fecundity Atlantic cod (\emph{Gadus morhua}). Our results
for two subclasses of Lambda coalescents show that one can distinguish these
subclasses from the Kingman coalescent, as well as between the
Lambda-subclasses, even for moderate sample sizes.


  We present a semi-empirical (PM6-based) computational method for
systematically estimating the effect of all possible single mutants, within a
certain radius of the active site, on the barrier height of an enzymatic
reaction. The intent of this method is not a quantitative prediction of the
barrier heights, but rather to identify promising mutants for further
computational or experimental study. The method is applied to identify
promising single and double mutants of Bacillus circulans xylanase (BCX) with
increased hydrolytic activity for the artificial substrate ortho-nitrophenyl
\beta-xylobioside (ONPX$_2$). The estimated reaction barrier for wild-type (WT)
BCX is 18.5 kcal/mol, which is in good agreement with the experimental
activation free energy value of 17.0 kcal/mol extracted from the observed
k$_\text{cat}$ using transition state theory (Joshi et al., Biochemistry 2001,
40, 10115). The PM6 reaction profiles for eight single point mutations are
recomputed using FMO-MP2/PCM/6-31G(d) single points. PM6 predicts an increase
in barrier height for all eight mutants while FMO predicts an increase for six
of the eight mutants. Both methods predict that the largest change in barrier
occurs for N35F, where PM6 and FMO predict a 9.0 and 15.8 kcal/mol increase,
respectively. We thus conclude that PM6 is sufficiently accurate to identify
promising mutants for further study. We prepared a set of all theoretically
possible (342) single mutants in which every amino acid of the active site
(except for the catalytically active residues E78 and E172) was mutated to
every other amino acid. Based on results from the single mutants we construct a
set of 111 double mutants consisting of all possible pairs of single mutants
with the lowest barrier for a particular position and compute their reaction
profile. None of the mutants have, to our knowledge, been prepared
experimentally[...].


  Adaptation is a central topic in theoretical biology, of practical importance
for analyzing drug resistance mutations. Several authors have used arguments
based on extreme value theory in their work on adaptation. There are
complications with these approaches if fitness is additive (meaning that
fitness effects of mutations sum), or whenever there is more additivity than
what one would expect in an uncorrelated fitness landscape. However, the
approaches have been used in published work, even in situations with
substantial amounts of additivity. In particular, extreme value theory has been
used in discussions on evolutionary predictability. We say that evolution is
predictable if the use of a particular drug at different locations tends lead
to the same resistance mutations. Evolutionary predictability depends on the
probabilities of mutational trajectories. Arguments about probabilities based
on extreme value theory can be misleading. Additivity may cause errors in
estimates of the probabilities of some mutational trajectories by a factor 20
even for rather small examples. We show that additivity gives systematic errors
so as to exaggerate the differences between the most and the least likely
trajectory. As a result of this bias, evolution may appear more predictable
than it is. From a broader perspective, our results suggest that approaches
which depend on the Orr-Gillespie theory are likely to give misleading results
for realistic fitness landscapes whenever one considers adaptation in several
steps.


  RNA-seq allows detection and precise quantification of transcripts, provides
comprehensive understanding of exon/intron boundaries, aids discovery of
alternatively spliced isoforms and fusion transcripts along with measurement of
allele-specific expression. Researchers interested in studying and constructing
transcriptomes, especially for non-model species, often face the conundrum of
choosing from a number of available de novo and genome-guided assemblers. A
comprehensive comparative study is required to assess and evaluate their
efficiency and sensitivity for transcript assembly, reconstruction and
recovery. None of the popular assembly tools in use today achieves requisite
sensitivity, specificity or recovery of full-length transcripts on its own.
Hence, it is imperative that methods be developed in order to augment
assemblies generated from multiple tools, with minimal compounding of error.
Here, we present an approach to combinatorially augment transciptome assembly
based on a rigorous comparative study of popular de novo and genome-guided
transcriptome assembly tools.


  The main signalling proteins that control early kidney branching have been
defined. Yet the underlying mechanism is still elusive. We have previously
shown that a Schnakenberg-type Turing mechanism can recapitulate the branching
and protein expression patterns in wildtype and mutant lungs, but it is unclear
whether this mechanism would extend to other branched organs that are regulated
by other proteins. Here we show that the GDNF-RET regulatory interaction gives
rise to a Schnakenberg-type Turing model that reproduces the observed budding
of the ureteric bud from the Wolffian duct, its invasion into the mesenchyme,
and the observed branching pattern. The model also recapitulates all relevant
protein expression patterns in wild-type and mutant mice. The lung and kidney
models are both based on a particular receptor- ligand interaction and require:
(1) cooperative binding of ligand and receptor, (2) a lower diffusion
coefficient for the receptor than for the ligand, and (3) an increase in the
receptor concentration in response to receptor-ligand binding (by enhanced
transcription, more recycling or similar). These conditions are met also by
other receptor-ligand systems. We propose that ligand-receptor based Turing
patterns represent a general mechanism to control branching morphogenesis and
other developmental processes.


  Transcriptional delay can significantly impact the dynamics of gene networks.
Here we examine how such delay affects bistable systems. We investigate several
stochastic models of bistable gene networks and find that increasing delay
dramatically increases the mean residence times near stable states. To explain
this, we introduce a non-Markovian, analytically tractable reduced model. The
model shows that stabilization is the consequence of an increased number of
failed transitions between stable states. Each of the bistable systems that we
simulate behaves in this manner.


  Metabolic control analysis is a biochemical formalism defined by Kacser and
Burns in 1973, and given firm mathematical basis by Reder in 1988. The
algorithm defined by Reder for calculating the control matrices is still used
by software programs today, but is only valid for some biochemical models. We
show that, with slight modification, the algorithm may be applied to all
models.


  The identification of the H3K4 trimethylase, PRDM9, as the gene responsible
for recombination hotspot localization has provided considerable insight into
the mechanisms by which recombination is initiated in mammals. However,
uniquely amongst mammals, canids appear to lack a functional version of PRDM9
and may therefore provide a model for understanding recombination that occurs
in the absence of PRDM9, and thus how PRDM9 functions to shape the
recombination landscape. We have constructed a fine-scale genetic map from
patterns of linkage disequilibrium assessed using high-throughput sequence data
from 51 free-ranging dogs, Canis lupus familiaris. While broad-scale properties
of recombination appear similar to other mammalian species, our fine-scale
estimates indicate that canine highly elevated recombination rates are observed
in the vicinity of CpG rich regions including gene promoter regions, but show
little association with H3K4 trimethylation marks identified in spermatocytes.
By comparison to genomic data from the Andean fox, Lycalopex culpaeus, we show
that biased gene conversion is a plausible mechanism by which the high CpG
content of the dog genome could have occurred.


  Consider a population where individuals give birth at constant rate during
their lifetimes to i.i.d. copies of themselves. Individuals bear clonally
inherited types, but (neutral) mutations may happen at the birth events. The
smallest subtree containing the genealogy of all the extant individuals at a
fixed time \tau, is called the coalescent point process. We enrich this process
with the history of the mutations that appeared over time, and call it the
marked coalescent point process. With the help of limit theorems for L\'evy
processes with marked jumps established in a previous work (arXiv:1305.6245),
we prove the convergence of the marked coalescent point process with large
population size and two possible regimes for the mutations - one of them being
a classical rare mutation regime, towards a multivariate Poisson point process.
This Poisson point process can be described as the coalescent point process of
the limiting population at \tau, with mutations arising as inhomogeneous
regenerative sets along the lineages. Its intensity measure is further
characterized thanks to the excursion theory for spectrally positive L\'evy
processes. In the rare mutations asymptotic, mutations arise as the image of a
Poisson process by the ladder height process of a L\'evy process with infinite
variation, and in the particular case of the critical branching process with
exponential lifetimes, the limiting object is the Poisson point process of the
depths of excursions of the Brownian motion, with Poissonian mutations on the
lineages.


  This review presents how R, the popular statistical environment and
programming language, can be used in the frame of proteomics data analysis. A
short introduction to R is given, with special emphasis on some of the features
that make R and its add-on packages a premium software for sound and
reproducible data analysis. The reader is also advised on how to find relevant
R software for proteomics. Several use cases are then presented, illustrating
data input/output, quality control, quantitative proteomics and data analysis.
Detailed code and additional links to extensive documentation are available in
the freely available companion package RforProteomics.


  We present a novel approach allowing the study of rare events like fixation
under fluctuating environments, modeled as extrinsic noise, in evolutionary
processes characterized by the dominance of one species. Our treatment consists
of mapping the system onto an auxiliary model, exhibiting metastable species
coexistence, that can be analyzed semiclassically. This approach enables us to
study the interplay between extrinsic and demographic noise on the statistics
of interest. We illustrate our theory by considering the paradigmatic
prisoner's dilemma game whose evolution is described by the probability that
cooperators fixate the population and replace all defectors. We analytically
and numerically demonstrate that extrinsic noise may drastically enhance the
cooperation fixation probability and even change its functional dependence on
the population size. These results, which generalize earlier works in
population genetics, indicate that extrinsic noise may help sustain and promote
a much higher level of cooperation than static settings.


  The expression of genes usually follows a two-step procedure. First, a gene
(encoded in the genome) is transcribed resulting in a strand of (messenger)
RNA. Afterwards, the RNA is translated into protein. Classically, this gene
expression is modeled using a Markov jump process including activation and
deactivation of the gene, transcription and translation rates together with
degradation of RNA and protein. We extend this model by adding delays (with
arbitrary distributions) to transcription and translation. Such delays can
e.g.\ mean that RNA has to be transported to a different part of a cell before
translation can be initiated. Already in the classical model, production of RNA
and protein come in bursts by activation and deactivation of the gene,
resulting in a large variance of the number of RNA and proteins in equilibrium.
We derive precise formulas for this second-order structure with the model
including delay in equilibrium. As a general fact, the delay decreases the
variance of the number of RNA and proteins.


  Starting from the lead structure we have identified in our previous works, we
are extending our insight understanding of its potential inhibitory effect
against both EGFR and HER2 receptors. Herein and using extended molecular
dynamic simulations and different scoring techniques, we are providing
plausible explanations for the observed inhibitory effect. Also, we are
comparing the binding mechanism in addition to the dynamics of binding with two
other approved inhibitors against EGFR (Lapatinib) and HER2 (SYR). Based on
this information, we are also designing and in silico screening new potential
inhibitors sharing the same scaffold of the lead structure. We have chosen the
best scoring inhibitor for additional in silico investigation against both the
wild-type and T790M mutant strain of EGFR. It seems that certain substitution
pattern guarantees the binding to the conserved water molecule commonly
observed with kinase crystal structures. Also, the new inhibitors seem to form
a stable interaction with the mutant strain as a direct consequence of their
enhanced ability to form additional interactions with binding site residues.


  The realization that complex systems such as ecological communities can
collapse or shift regimes suddenly and without rapid external forcing poses a
serious challenge to our understanding and management of the natural world. The
potential to identify early warning signals that would allow researchers and
managers to predict such events before they happen has therefore been an
invaluable discovery that offers a way forward in spite of such seemingly
unpredictable behavior. Research into early warning signals has demonstrated
that it is possible to define and detect such early warning signals in advance
of a transition in certain contexts. Here we describe the pattern emerging as
research continues to explore just how far we can generalize these results. A
core of examples emerges that shares three properties: the phenomenon of rapid
regime shifts, a pattern of 'critical slowing down' that can be used to detect
the approaching shift, and a mechanism of bifurcation driving the sudden
change. As research has expanded beyond these core examples, it is becoming
clear that not all systems that show regime shifts exhibit critical slowing
down, or vice versa. Even when systems exhibit critical slowing down,
statistical detection is a challenge. We review the literature that explores
these edge cases and highlight the need for (a) new early warning behaviors
that can be used in cases where rapid shifts do not exhibit critical slowing
down, (b) the development of methods to identify which behavior might be an
appropriate signal when encountering a novel system; bearing in mind that a
positive indication for some systems is a negative indication in others, and
(c) statistical methods that can distinguish between signatures of early
warning behaviors and noise.


  Queues formed by social wasps to inherit the dominant position in the nest
are analyzed by using a transient quasi-birth-and-death process. We show that
the extended nest life time due to division of labor between queen and helpers
has a big impact for the nest productivity.


  Motivation: Transcriptome sequencing has long been the favored method for
quickly and inexpensively obtaining the sequences for a large number of genes
from an organism with no reference genome. With the rapidly increasing
throughputs and decreasing costs of next generation sequencing, RNA-Seq has
gained in popularity; but given the typically short reads (e.g. 2 x 90 bp
paired ends) of this technol- ogy, de novo assembly to recover complete or
full-length transcript sequences remains an algorithmic challenge. Results: We
present SOAPdenovo-Trans, a de novo transcriptome assembler designed
specifically for RNA-Seq. Its performance was evaluated on transcriptome
datasets from rice and mouse. Using the known transcripts from these
well-annotated genomes (sequenced a decade ago) as our benchmark, we assessed
how SOAPdenovo- Trans and two other popular software handle the practical
issues of alternative splicing and variable expression levels. Our conclusion
is that SOAPdenovo-Trans provides higher contiguity, lower redundancy, and
faster execution. Availability and Implementation: Source code and user manual
are at http://sourceforge.net/projects/soapdenovotrans/ Contact:
xieyl@genomics.cn or bgi-soap@googlegroups.com


  Movement behavior of an indicator species, zebrafish (Danio rerio), was
analyzed with one- and two-individual groups before and after treatment with a
toxic chemical, formaldehyde, at a low concentration (1 ppm). After the
boundary area had been determined based on experimental data, intermittency was
defined as the probability distributions of the shadowing time during which
data were above a pre-determined threshold and were obtained from experimental
time-series data on forces and the inter-distances for one and two individuals.
Overall intermittencies were similar in the boundary and central areas.
However, the intermittencies were remarkably different between the one- and the
two-individual groups: the single line was used to fit the data for the
one-individual group whereas two phases were observed with breakpoints
(approximately 10 seconds in logarithm) in the exponential fitting curves for
the two-individual group. A difference in the probability distributions of
shadowing time was observed "before" and "after" treatment for different areas.
Intermittency patterns before and after treatment were contrasted in the center
for the one-individual group whereas the difference was observed in the
boundary for two-individual group. The intermittencies for the inter-distances
of two individuals in the boundary and central areas were markedly different
before and after treatment. When the differences between the intermittencies in
the boundary and the central areas and between "before" and "after" treatment
are considered, the distribution patterns of the shadowing time (scaling
behaviors or intermittency patterns) should be a useful means of bio-monitoring
to detect contaminants in the environment.


  The accumulation of beneficial mutations on many competing genetic
backgrounds in rapidly adapting populations has a striking impact on
evolutionary dynamics. This effect, known as clonal interference, causes
erratic fluctuations in the frequencies of observed mutations, randomizes the
fixation times of successful mutations, and leaves distinct signatures on
patterns of genetic variation. Here, we show how this form of `genetic draft'
affects the forward-time dynamics of site frequencies in rapidly adapting
asexual populations. We calculate the probability that mutations at individual
sites shift in frequency over a characteristic timescale, extending Gillespie's
original model of draft to the case where many strongly selected beneficial
mutations segregate simultaneously. We then derive the sojourn time of mutant
alleles, the expected fixation time of successful mutants, and the site
frequency spectrum of beneficial and neutral mutations. We show how this form
of draft affects inferences in the McDonald-Kreitman test, and how it relates
to recent observations that some aspects of genetic diversity are described by
the Bolthausen-Sznitman coalescent in the limit of very rapid adaptation.
Finally, we describe how our method can be extended to model evolution on
fitness landscapes that include some forms of epistasis, such as landscapes
that are partitioned into two or more incompatible evolutionary trajectories.


  Phylogenomic analyses are subject to bias from compositional convergence and
noise from horizontal gene transfer (HGT). Compositional convergence is a
likely cause of controversy regarding phylogeny of the SAR11 group of
Alphaproteobacteria that have extremely streamlined, A+T-biased genomes. While
careful modeling can reduce artifacts caused by convergence, the most
consistent and robust phylogenetic signal in genomes may lie distributed among
encoded functional features that govern macromolecular interactions. Here we
develop a novel phyloclassification method based on signatures derived from
bioinformatically defined tRNA Class-Informative Features (CIFs). tRNA CIFs are
enriched for features that underlie tRNA-protein interactions. Using a simple
tRNA-CIF-based phyloclassifier, we obtained results consistent with those of
bias-corrected whole proteome phylogenomic studies, rejecting monophyly of
SAR11 and affiliating most strains with Rhizobiales with strong statistical
support. Yet SAR11 and Rickettsiales tRNA genes share distinct patterns of
A+T-richness, as expected from their elevated genomic A+T compositions. Using
conventional supermatrix methods on total tRNA sequence data, we could recover
the artifactual result of a monophyletic SAR11 grouping with Rickettsiales.
Thus tRNA CIF-based phyloclassification is more robust to base content
convergence than supermatrix phylogenomics on whole tRNA sequences. Also, given
the notoriously promiscuous HGT of aminoacyl-tRNA synthetases, tRNA CIF-based
phyloclassification may be relatively robust to HGT of network components. We
describe how unique features of tRNA-protein interaction networks facilitate
the mining of traits governing macromolecular interactions from genomic data,
and discuss why interaction-governing traits may be especially useful to solve
difficult problems in microbial classification and phylogeny.


  Evolutionary game dynamics with two 2-strategy games in a finite population
has been investigated in this study. Traditionally, frequency-dependent
evolutionary dynamics are modeled by deterministic replicator dynamics under
the assumption that the population size is infinite. However, in reality,
population sizes are finite. Recently, stochastic processes in finite
populations have been introduced into evolutionary games in order to study
finite size effects in evolutionary game dynamics. However, most of these
studies focus on populations playing only single games. In this study, we
investigate a finite population with two games and show that a finite
population playing two games tends to evolve toward a specific direction to
form particular linkages between the strategies of the two games.


  Understanding the dynamics of metapopulations close to extinction is of vital
importance for management. Levins-like models, in which local patches are
treated as either occupied or empty, have been used extensively for this
purpose, but they ignore the important role of local population dynamics. In
this paper, we consider a stochastic metapopulation model where local
populations follow a Ricker dynamics, and use this framework to investigate the
behaviour of the metapopulation at the brink of extinction. As long as
dispersal rates are not too large, the system is shown to have a time evolution
consistent with Levins' dynamics. We derive analytical expressions for the
colonisation and extinction rates (c and e) in Levins-type models in terms of
reproduction, survival, and dispersal parameters of the local populations,
providing an avenue to parameterising Levins-like models from the type of
information on local demography that is available for a number of species. To
facilitate applying our results, we provide a numerical implementation for
computing c and e.


  To identify genetic changes underlying dog domestication and reconstruct
their early evolutionary history, we analyzed novel high-quality genome
sequences of three gray wolves, one from each of three putative centers of dog
domestication, two ancient dog lineages (Basenji and Dingo) and a golden jackal
as an outgroup. We find dogs and wolves diverged through a dynamic process
involving population bottlenecks in both lineages and post-divergence gene
flow, which confounds previous inferences of dog origins. In dogs, the
domestication bottleneck was severe involving a 17 to 49-fold reduction in
population size, a much stronger bottleneck than estimated previously from less
intensive sequencing efforts. A sharp bottleneck in wolves occurred soon after
their divergence from dogs, implying that the pool of diversity from which dogs
arose was far larger than represented by modern wolf populations. Conditional
on mutation rate, we narrow the plausible range for the date of initial dog
domestication to an interval from 11 to 16 thousand years ago. This period
predates the rise of agriculture, implying that the earliest dogs arose
alongside hunter-gathers rather than agriculturists. Regarding the geographic
origin of dogs, we find that surprisingly, none of the extant wolf lineages
from putative domestication centers are more closely related to dogs, and the
sampled wolves instead form a sister monophyletic clade. This result, in
combination with our finding of dog-wolf admixture during the process of
domestication, suggests a re-evaluation of past hypotheses of dog origin is
necessary. Finally, we also detect signatures of selection, including evidence
for selection on genes implicated in morphology, metabolism, and neural
development. Uniquely, we find support for selective sweeps at regulatory sites
suggesting gene regulatory changes played a critical role in dog domestication.


  The construction of a network of cell-to-cell contacts makes it possible to
characterize the patterns and spatial organisation of tissues. Such networks
are highly dynamic, depending on the changes of the tissue architecture caused
by cell division, death and migration. Local competitive and cooperative
cell-to-cell interactions influence the choices cells make. We review the
literature on quantitative data of epithelial tissue topology and present a
dynamical network model that can be used to explore the evolutionary dynamics
of a two dimensional tissue architecture with arbitrary cell-to-cell
interactions. In particular, we show that various forms of experimentally
observed types of interactions can be modelled using game theory. We discuss a
model of cooperative and non-cooperative cell-to-cell communication that can
capture the interplay between cellular competition and tissue dynamics. We
conclude with an outlook on the possible uses of this approach in modelling
tumorigenesis and tissue homeostasis.


  Here we use Drosophila melanogaster to create a genetic model of human
permanent neonatal diabetes mellitus and present experimental results
describing dimensions of this complexity. The approach involves the transgenic
expression of a misfolded mutant of human preproinsulin, hINSC96Y, which is a
cause of the disease. When expressed in fly imaginal discs, hINSC96Y causes a
reduction of adult structures, including the eye, wing and notum. Eye imaginal
discs exhibit defects in both the structure and arrangement of ommatidia. In
the wing, expression of hINSC96Y leads to ectopic expression of veins and
mechano-sensory organs, indicating disruption of wild type signaling processes
regulating cell fates. These readily measurable disease phenotypes are
sensitive to temperature, gene dose and sex. Mutant (but not wild type)
proinsulin expression in the eye imaginal disc induces IRE1-mediated Xbp1
alternative splicing, a signal for endoplasmic reticulum stress response
activation, and produces global change in gene expression. Mutant hINS
transgene tester strains, when crossed to stocks from the Drosophila Genetic
Reference Panel produces F1 adults with a continuous range of disease
phenotypes and large broad-sense heritability. Surprisingly, the severity of
mutant hINS-induced disease in the eye is not correlated with that in the notum
in these crosses, nor with eye reduction phenotypes caused by the expression of
two dominant eye mutants acting in two different eye development pathways, Drop
(Dr) or Lobe (L) when crossed into the same genetic backgrounds. The tissue
specificity of genetic variability for mutant hINS-induced disease thus has its
own distinct signature. The genetic dominance of disease-specific phenotypic
variability makes this approach amenable to genome-wide association study
(GWAS) in a simple F1 screen of natural variation.


  Often a screening or selection experiment targets a cell or tissue, which
presents many possible molecular targets and identifies a correspondingly large
number of ligands. We describe a statistical method to extract an estimate of
the complexity or richness of the set of molecular targets from competition
experiments between distinguishable ligands, including aptamers derived from
combinatorial experiments (SELEX or phage display). In simulations, the
nonparametric statistic provides a robust estimate of complexity from a 100x100
matrix of competition experiments, which is clearly feasible in high-throughput
format. The statistic and method are potentially applicable to other ligand
binding situations.


  A recent study reported that there is evidence life may have originated prior
to the formation of the Earth. That conclusion was based on a regression
analysis of a certain data set involving evolution of functional genome size
across major phyla. Here it is shown that if measurement errors and
"confidence" intervals are taken into account, then the regression analysis of
the same data set leads to conclusions that allow for life to have appeared
after the formation of the Earth.


  Migration is a fundamental trait in humans and animals. Recent studies
investigated the effect of migration on the evolution of cooperation, showing
that contingent migration favors cooperation in spatial structures. In those
studies, only local migration to immediate neighbors was considered, while
long-range migration has not been considered yet, partly because the long-range
migration has been generally regarded as harmful for cooperation as it would
bring the population to a well-mixed state that favors defection. Here, we
studied the effects of adaptive long-range migration on the evolution of
cooperation through agent-based simulations of a spatial Prisoner's Dilemma
game where individuals can jump to a farther site if they are surrounded by
more defectors. Our results show that adaptive long-range migration strongly
promotes cooperation, especially under conditions where the temptation to
defect is considerably high. These findings demonstrate the significance of
adaptive long-range migration for the evolution of cooperation.


  A particular case of the famous quasispecies model - the Crow-Kimura model
with a permutation invariant fitness landscape - is investigated. Using the
fact that the mutation matrix in the case of a permutation invariant fitness
landscape has a special tridiagonal form, a change of the basis is suggested
such that in the new coordinates a number of analytical results can be
obtained. In particular, using the eigenvectors of the mutation matrix as the
new basis, we show that the quasispecies distribution approaches a binomial one
and give simple estimates for the speed of convergence. Another consequence of
the suggested approach is a parametric solution to the system of equations
determining the quasispecies. Using this parametric solution we show that our
approach leads to exact asymptotic results in some cases, which are not covered
by the existing methods. In particular, we are able to present not only the
limit behavior of the leading eigenvalue (mean population fitness), but also
the exact formulas for the limit quasispecies eigenvector for special cases.
For instance, this eigenvector has a geometric distribution in the case of the
classical single peaked fitness landscape. On the biological side, we propose a
mathematical definition, based on the closeness of the quasispecies to the
binomial distribution, which can be used as an operational definition of the
notorious error threshold. Using this definition, we suggest two approximate
formulas to estimate the critical mutation rate after which the quasispecies
delocalization occurs.


  How high should be the rate of immigration into a stochastic population in
order to significantly reduce the probability of observing the population
extinct? Is there any relation between the population size distributions with
and without immigration? Under what conditions can one justify the simple patch
occupancy models which ignore the population distribution and its dynamics in a
patch, and treat a patch simply as either occupied or empty? We address these
questions by exactly solving a simple stochastic model obtained by adding a
steady immigration to a variant of the Verhulst model: a prototypical model of
an isolated stochastic population.


  A mechanism of the replication of proto-RNAs in oligomer world is proposed.
The replication is carried out by a minimum cycle which is sustained by a
ligase and a helicase. We expect that such a cycle actually worked in the
primordial soup and can be constructed in vitro. By computer simulation the
products of the replication acquires diversity and complexity. Such diversity
and complexity are the bases of the evolution.


  We show that numerical approximations of Kolmogorov complexity (K) applied to
graph adjacency matrices capture some group-theoretic and topological
properties of graphs and empirical networks ranging from metabolic to social
networks. That K and the size of the group of automorphisms of a graph are
correlated opens up interesting connections to problems in computational
geometry, and thus connects several measures and concepts from complexity
science. We show that approximations of K characterise synthetic and natural
networks by their generating mechanisms, assigning lower algorithmic randomness
to complex network models (Watts-Strogatz and Barabasi-Albert networks) and
high Kolmogorov complexity to (random) Erdos-Renyi graphs. We derive these
results via two different Kolmogorov complexity approximation methods applied
to the adjacency matrices of the graphs and networks. The methods used are the
traditional lossless compression approach to Kolmogorov complexity, and a
normalised version of a Block Decomposition Method (BDM) measure, based on
algorithmic probability theory.