  Oxygen is released to living tissues via conformational changes of hemoglobin
from R-state (oxyhemoglobin) to T-state (desoxyhemoglobin). The detailed
mechanism of this process is not yet fully understood. We have carried out
micromechanical experiments on oxyhemoglobin crystals to determine the behavior
of the Young's modulus and the internal friction for temperatures between 20 C
and 70 C. We have found that around 49 C oxyhemoglobin crystal samples undergo
a sudden and strong increase of their Young's modulus, accompanied by a sudden
decrease of the internal friction. This sudden mechanical change (proteinquake)
takes place in a partially unfolded state and precedes the full denaturation
transition at higher temperatures. The hemoglobin crystals after the
proteinquake has the same mechanical properies as the initial state at room
temperatures. We conjecture that it can be relevant for explaining the
oxygen-releasing function of native oxyhemoglobin when the temperature is
increased, e.g. due to active sport. The effect is specific for the quaternary
structure of hemoglobin, and is absent for myoglobin with only one peptide
sequence.


  We propose an extension of a standard stochastic individual-based model in
population dynamics which broadens the range of biological applications. Our
primary motivation is modelling of immunotherapy of malignant tumours. In this
context the different actors, T-cells, cytokines or cancer cells, are modelled
as single particles (individuals) in the stochastic system. The main expansions
of the model are distinguishing cancer cells by phenotype and genotype,
including environment-dependent phenotypic plasticity that does not affect the
genotype, taking into account the effects of therapy and introducing a
competition term which lowers the reproduction rate of an individual in
addition to the usual term that increases its death rate. We illustrate the new
setup by using it to model various phenomena arising in immunotherapy. Our aim
is twofold: on the one hand, we show that the interplay of genetic mutations
and phenotypic switches on different timescales as well as the occurrence of
metastability phenomena raise new mathematical challenges. On the other hand,
we argue why understanding purely stochastic events (which cannot be obtained
with deterministic models) may help to understand the resistance of tumours to
therapeutic approaches and may have non-trivial consequences on tumour
treatment protocols. This is supported through numerical simulations.


  Discoveries of new biomarkers for frequently occurring diseases are of
special importance in today's medicine. While fully developed type II diabetes
(T2D) can be detected easily, the early identification of high risk individuals
is an area of interest in T2D, too. Metagenomic analysis of the human bacterial
flora has shown subtle changes in diabetic patients, but no specific microbes
are known to cause or promote the disease. Moderate changes were also detected
in the microbial gene composition of the metagenomes of diabetic patients, but
again, no specific gene was found that is present in disease-related and
missing in healthy metagenome. However, these fine differences in microbial
taxon- and gene composition are difficult to apply as quantitative biomarkers
for diagnosing or predicting type II diabetes. In the present work we report
some nucleotide 9-mers with significantly differing frequencies in diabetic and
healthy intestinal flora. To our knowledge, it is the first time such short DNA
fragments have been associated with T2D. The automated, quantitative analysis
of the frequencies of short nucleotide sequences seems to be more feasible than
accurate phylogenetic and functional analysis, and thus it might be a promising
direction of diagnostic research.


  We present a method that can simultaneously locate positions of overlapped
multi-emitters at the theoretical-limit precision. We derive a set of simple
equations whose solution gives the maximum likelihood estimator of
multi-emitter positions. We compare the performance of our simultaneous
localization analysis with the conventional single-molecule analysis for
simulated images and show that our method can improve the time-resolution of
superresolution microscopy an order of magnitude. In particular, we derive the
information-theoretical bound on time resolution of localization-based
superresolution microscopy and demonstrate that the bound can be closely
attained by our analysis.


  Epidemic dynamics in a stochastic network of interacting epidemic centers is
considered. The epidemic and migration processes are modelled by Markov's
chains. Explicit formulas for probability distribution of the migration process
are derived. Dependence of outbreak parameters on initial parameters,
population, coupling parameters is examined analytically and numerically. The
mean field approximation for a general migration process is derived. An
approximate method allowing computation of statistical moments for networks
with highly populated centres is proposed and tested numerically.


  High-throughput RNA sequencing (RNA-seq) is now the standard method to
determine differential gene expression. Identifying differentially expressed
genes crucially depends on estimates of read count variability. These estimates
are typically based on statistical models such as the negative binomial
distribution, which is employed by the tools edgeR, DESeq and cuffdiff. Until
now, the validity of these models has usually been tested on either
low-replicate RNA-seq data or simulations. A 48-replicate RNA-seq experiment in
yeast was performed and data tested against theoretical models. The observed
gene read counts were consistent with both log-normal and negative binomial
distributions, while the mean-variance relation followed the line of constant
dispersion parameter of ~0.01. The high-replicate data also allowed for strict
quality control and screening of bad replicates, which can drastically affect
the gene read-count distribution. RNA-seq data have been submitted to ENA
archive with project ID PRJEB5348.


  Many events in the vertebrate immune system are influenced by some element of
chance. The objective of the present work is to describe affinity maturation of
B lymphocytes (in which random events are perhaps the most characteristic), and
to study a possible network model of immune memory. In our model stochastic
processes govern all events. A major novelty of this approach is that it
permits studying random variations in the immune process. Four basic components
are simulated in the model: non-immune self cells, nonself cells (pathogens), B
lymphocytes, and bone marrow cells that produce naive B lymphocytes. A point in
a generalized shape space plus the size of the corresponding population
represents nonself and non-immune self cells. On the other hand, each
individual B cell is represented by a disc that models its recognition region
in the shape space. Infection is simulated by an "injection" of nonself cells
into the system. Division of pathogens may instigate an attack of naive B
cells, which in turn may induce clonal proliferation and hypermutation in the
attacking B cells, and which eventually may slow down and stop the exponential
growth of pathogens. Affinity maturation of newly produced B cells becomes
expressed as a result of selection when the number of pathogens decreases.
Under favorable conditions, the expanded primary B cell clones may stimulate
the expansion of secondary B cell clones carrying complementary receptors to
the stimulating B cells. Like in a hall of mirrors, the image of pathogens in
the primary B cell clones then will be reflected in secondary B cell clones.
This "ping-pong" game may survive for a long time even in the absence of the
pathogen, creating a local network memory. This memory ensures that repeated
infection by the same pathogen will be eliminated more efficiently.


  Contact inhibition is the process by which cells switch from a motile growing
state to a passive and stabilized state upon touching their neighbors. When two
cells touch, an adhesion link is created between them by means of transmembrane
E-cadherin proteins. Simultaneously, their actin filaments stop polymerizing in
the direction perpendicular to the membrane and reorganize to create an apical
belt that colocalizes with the adhesion links. Here, we propose a detailed
quantitative model of the role of the cytoplasmic $\beta$-catenin and
$\alpha$-catenin proteins in this process, treated as a reaction-diffusion
system. Upon cell-cell contact, the concentration in $\alpha$-catenin dimers
increases, inhibiting actin branching and thereby reducing cellular motility
and expansion pressure. This model provides a mechanism for contact inhibition
that could explain previously unrelated experimental findings on the role
played by E-cadherin, $\beta$-catenin and $\alpha$-catenin in the cellular
phenotype and in tumorigenesis. In particular, we address the effect of a
knockout of the adenomatous polyposis coli tumor suppressor gene. Potential
direct tests of our model are discussed.


  With the advent of high-throughput wet lab technologies the amount of protein
interaction data available publicly has increased substantially, in turn
spurring a plethora of computational methods for in silico knowledge discovery
from this data. In this paper, we focus on parameterized methods for modeling
and solving complex computational problems encountered in such knowledge
discovery from protein data. Specifically, we concentrate on three relevant
problems today in proteomics, namely detection of lethal proteins, functional
modules and alignments from protein interaction networks. We propose novel
graph theoretic models for these problems and devise practical parameterized
algorithms. At a broader level, we demonstrate how these methods can be viable
alternatives for the several heurestic, randomized, approximation and
sub-optimal methods by arriving at parameterized yet optimal solutions for
these problems. We substantiate these theoretical results by experimenting on
real protein interaction data of S. cerevisiae (budding yeast) and verifying
the results using gene ontology.


  Although foot-and-mouth disease virus (FMDV) incidence has decreased in South
America over the last years, the pathogen still circulates in the region and
the risk of re-emergence in previously FMDV-free areas is a veterinary public
health concern. In this paper we merge environmental, epidemiological and
genetic data to reconstruct spatiotemporal patterns and determinants of FMDV
serotypes A and O dispersal in South America. Our dating analysis suggests that
serotype A emerged in South America around 1930, while serotype O emerged
around 1990. The rate of evolution for serotype A was significantly higher
compared to serotype O. Phylogeographic inference identified two well-connected
sub networks of viral flow, one including Venezuela, Colombia and Ecuador;
another including Brazil, Uruguay and Argentina. The spread of serotype A was
best described by geographic distances, while trade of live cattle was the
predictor that best explained serotype O spread. Our findings show that the two
serotypes have different underlying evolutionary and spatial dynamics and may
pose different threats to control programmes. Key-words: Phylogeography,
foot-and-mouth disease virus, South America, animal trade.


  Protein-peptide interactions play essential functional roles in living
organisms and their structural characterization is a hot subject of current
experimental and theoretical research. Computational modeling of the structure
of protein-peptide interactions is usually divided into two stages: prediction
of the binding site at a protein receptor surface, and then docking (and
modeling) the peptide structure into the known binding site. This paper
presents a comprehensive CABS-dock method for the simultaneous search of
binding sites and flexible protein-peptide docking, available as a users
friendly web server. We present example CABS-dock results obtained in the
default CABS-dock mode and using its advanced options that enable the user to
increase the range of flexibility for chosen receptor fragments or to exclude
user-selected binding modes from docking search. Furthermore, we demonstrate a
strategy to improve CABS-dock performance by assessing the quality of models
with classical molecular dynamics. Finally, we discuss the promising extensions
and applications of the CABS-dock method and provide a tutorial appendix for
the convenient analysis and visualization of CABS-dock results. The CABS-dock
web server is freely available at http://biocomp.chem.uw.edu.pl/CABSdock/


  We have developed an analytical, ligand-specific and scalable algorithm that
detects a "signature" of the 3D binding site of a given ligand in a protein 3D
structure. The said signature is a 3D motif in the form of an irregular
tetrahedron whose vertices represent the backbone or side-chain centroids of
the amino acid residues at the binding site that physically interact with the
bound ligand atoms. The motif is determined from a set of solved training
structures, all of which bind the ligand. Just as alignment of linear amino
acid sequences enables one to determine consensus sequences in proteins, the
present method allows the determination of three-dimensional consensus
structures or "motifs" in folded proteins. Although such is accomplished by the
present method not by alignment of 3D protein structures or parts thereof
(e.g., alignment of ligand atoms from different structures) but by
nearest-neighbor analysis of ligand atoms in protein-bound forms, the same
effect, and thus the same goal, is achieved. We have applied our method to the
prediction of GTP- and ATP-binding protein families, namely, the small Ras-type
G-protein and ser/thr protein kinase families. Validation tests reveal that the
specificity of our method is nearly 100% for both protein families, and a
sensitivity of greater than 60% for the ser/thr protein kinase family and
approx. 93% for the small, Ras-type G-protein family. Further tests reveal that
our algorithm can distinguish effectively between GTP and GTP-like ligands, and
between ATP- and ATP-like ligands. The method was applied to a set of predicted
(by 123D threading) protein structures from the slime mold (D. dictyostelium)
proteome, with promising results.


  We describe two complementary methods to quantify the degree of burial of
ligand and/or ligand binding site (LBS) in a protein-ligand complex, namely,
the "cutting plane" (CP) and the "tangent sphere" (TS) methods. To construct
the CP and TS, two centroids are required: the protein molecular centroid
(global centroid, GC), and the LBS centroid (local centroid, LC). The CP is
defined as the plane passing through the LBS centroid (LC) and normal to the
line passing through the LC and the protein molecular centroid (GC). The
"anterior side" of the CP is the side not containing the GC (which the
"posterior" side does). The TS is defined as the sphere with center at GC and
tangent to the CP at LC. The percentage of protein atoms (a.) inside the TS,
and (b.) on the anterior side of the CP, are two complementary measures of
ligand or LBS burial depth since the latter is directly proportional to (b.)
and inversely proportional to (a.). We tested the CP and TS methods using a
test set of 67 well characterized protein-ligand structures (Laskowski et al.,
1996), as well as the theoretical case of an artificial protein in the form of
a cubic lattice grid of points in the overall shape of a sphere and in which
LBS of any depth can be specified. Results from both the CP and TS methods
agree very well with data reported by Laskowski et al., and results from the
theoretical case further confirm that that both methods are suitable measures
of ligand or LBS burial. Prior to this study, there were no such numerical
measures of LBS burial available, and hence no way to directly and objectively
compare LBS depths in different proteins. LBS burial depth is an important
parameter as it is usually directly related to the amount of conformational
change a protein undergoes upon ligand binding, and ability to quantify it
could allow meaningful comparison of protein dynamics and flexibility.


  Due to increased activity in high-throughput structural genomics efforts
around the globe, there has been an accumulation of experimental protein 3D
structures lacking functional annotation, thus creating a need for
structure-based protein function assignment methods. Computational prediction
of ligand binding sites (LBS) is a well-established protein function assignment
method. Here we apply the specific LBS detection algorithm we recently
described (Reyes, V.M. & Sheth, V.N., 2011; Reyes, V.M., 2015a) to some 801
functionally unannotated experimental structures in the Protein Data Bank by
screening for the binding sites (BS) of 6 biologically important ligands: GTP
in small Ras-type G-proteins, ATP in ser/thr protein kinases, sialic acid
(SIA), retinoic acid (REA), and heme-bound and unbound (free) nitric oxide
(hNO, fNO). Validation of the algorithm for the GTP- and ATP-binding sites has
been previously described in detail (ibid.); here, validation for the BSs of
the 4 other ligands shows both good specificity and sensitivity. Of the 801
structures screened, 8 tested positive for GTP binding, 61 for ATP binding, 35
for SIA binding, 132 for REA binding, 33 for hNO binding, and 10 for fNO
binding. Using the cutting plane and tangent sphere methods we described
previously, (Reyes, V.M., 2015b), we also determined the depth of burial of the
LBSs detected above and compared the values with those from the respective
training structures, and the degree of similarity between the two values taken
as a further validation of the predicted LBSs. Applying this criterion, we were
able to narrow down the predicted GTP-binding proteins to 2, the ATP-binding
proteins to 13, the SIA-binding proteins to 2, the REA-binding proteins to 14,
the hNO-binding proteins to 4, and the fNO-binding proteins to 1. We believe
this further criterion increases the confidence level of our LBS predictions.


  We report a 3D structure-based method of predicting protein-protein
interaction partners. It involves screening for pairs of tetrahedra
representing interacting amino acids at the interface of the protein-protein
complex, with one tetrahedron on each protomer. H-bonds and VDW interactions at
their interface are first determined and then interacting tetrahedral motifs
(one from each protomer) representing backbone or side chain centroids of the
interacting amino acids, are then built. The method requires that the protein
protomers be transformed first into double-centroid reduced representation
(Reyes, V.M. & Sheth, V.N., 2011; Reyes, V.M., 2015a). The method is applied to
a set of 801 protein structures in the PDB with unknown functions, which were
screened for pairs of tetrahedral motifs characteristic of nine binary
complexes, namely: (1.) RAP-Gmppnp-cRAF1 Ras-binding domain; (2.) RHOA-protein
kinase PKN/PRK1 effector domain; (3.) RAC-HOGD1; (4.) RAC-P67PHOX; (5.)
kinase-associated phosphatase (KAP)-phosphoCDK2; (6.) Ig Fc-protein A fragment
B; (7.) Ig light chain dimers; (8.) beta catenin-HTCF-4; and (9.) IL-2
homodimers. Our search method found 33, 297, 62, 63, 120, 0, 108, 16 and 504
putative complexes, respectively. After considering the degree of interface
overlap between the protomers, these numbers were significantly trimmed down to
4, 2, 1, 8, 3, 0, 1, 1 and 1, respectively. Negative and positive control
experiments indicate that the screening process has acceptable specificity and
sensitivity. The results were further validated by applying the CP and TS
methods (Reyes, V.M., 2015b) for the quantitative determination of interface
burial and inter-protomer overlap in the complex. Our method is simple, fast
and scalable, and once the partner interface 3D SMs are identified, they can be
used to computationally dock the two protomers together to form the complex.


  Background: While the importance of gene-gene interactions in human diseases
has been well recognized, identifying them has been a great challenge,
especially through association studies with millions of genetic markers and
thousands of individuals. Computationally efficient and powerful tools are in
great need for the identification of new gene-gene interactions in
high-dimensional association studies. Result: We develop C++ software for
genome-wide gene-gene interaction analyses (GWGGI). GWGGI utilizes tree-based
algorithms to search a large number of genetic markers for a disease-associated
joint association with the consideration of high-order interactions, and then
uses non-parametric statistics to test the joint association. The package
includes two functions, likelihood ratio Mann-whitney (LRMW) and Tree
Assembling Mann-whitney (TAMW).We optimize the data storage and computational
efficiency of the software, making it feasible to run the genome-wide analysis
on a personal computer. The use of GWGGI was demonstrated by using two real
data-sets with nearly 500 k genetic markers. Conclusion: Through the empirical
study, we demonstrated that the genome-wide gene-gene interaction analysis
using GWGGI could be accomplished within a reasonable time on a personal
computer (i.e., ~3.5 hours for LRMW and ~10 hours for TAMW). We also showed
that LRMW was suitable to detect interaction among a small number of genetic
variants with moderate-to-strong marginal effect, while TAMW was useful to
detect interaction among a larger number of low-marginal-effect genetic
variants.


  We examine the binding of transcription factors to DNA in terms of an
information transfer problem. The input of the noisy channel is the biophysical
signal of a factor bound to a DNA site, and the output is a distribution of
probable DNA sequences at this site. This task involves an inherent tradeoff
between the information gain and the energetics of the binding interaction -
high binding energies provide higher information gain but hinder the dynamics
of the system as factors are bound too tightly. We show that adaptation of the
binding interaction towards increasing information transfer under a general
energy constraint implies that the information gain per specific binding energy
at each base-pair is maximized. We analyze hundreds of prokaryote and eukaryote
transcription factors from various organisms to evaluate the discrimination
energies. We find that, in accordance with our theoretical argument, binding
energies nearly maximize the information gain per energy. This work suggests
the adaptation of information gain as a generic design principle of molecular
recognition systems.


  Here we study a driven lattice gas model for microtubule depolymerizing
molecular motors, where traffic jams of motors induce stochastic switching
between microtubule growth and shrinkage. We term this phenomenon
\enquote{traffic dynamic instability} because it is reminiscent of microtubule
dynamic instability [T. Mitchison and M. Kirschner, Nature 312, 237 (1984)].
The intermittent dynamics of growth and shrinking emerges from the interplay
between the arrival of motors at the microtubule tip, motor induced
depolymerization, and motor detachment from the tip. The switching dynamics
correlates with low and high motor density on the lattice. This leads to an
effectively bistable particle density in the system. A refined domain wall
theory predicts this transient appearance of different phases in the system.
The theoretical results are supported by stochastic simulations.


  There has been renewed interest in the use of sporozoite-based approaches for
malaria vaccination and controlled human infections, and several sets of human
challenge studies have recently completed. A study undertaken in Tanzania and
published in 2014 found dose-dependence between 10,000 and 25,000 sporozoite
doses, as well as divergent times-to-parasitemia relative to earlier studies in
European volunteers. However, this analysis shows that these conclusions are
based upon suboptimal analytical methods; with more optimal analysis, there is
no evidence for dose-dependence within this dose range; and more importantly,
no evidence for differences in event times between Dutch and Tanzanian study
sites. While these finding do not impact the reported safety and tolerability
of PfSPZ, they highlight critical issues that should be comprehensively
considered in future challenge studies.


  In recent years, many variants of percolation have been used to study network
structure and the behavior of processes spreading on networks. These include
bond percolation, site percolation, $k$-core percolation, bootstrap
percolation, the generalized epidemic process, and the Watts Threshold Model
(WTM). We show that --- except for bond percolation --- each of these processes
arises as a special case of the WTM and bond percolation arises from a small
modification. In fact "heterogeneous $k$-core percolation", a corresponding
"heterogeneous bootstrap percolation" model, and the generalized epidemic
process are completely equivalent to one another and the WTM. We further show
that a natural generalization of the WTM in which individuals "transmit" or
"send a message" to their neighbors with some probability less than $1$ can be
reformulated in terms of the WTM, and so this apparent generalization is in
fact not more general. Finally, we show that in bond percolation, finding the
set of nodes in the component containing a given node is equivalent to finding
the set of nodes activated if that node is initially activated and the node
thresholds are chosen from the appropriate distribution. A consequence of these
results is that mathematical techniques developed for the WTM apply to these
other models as well, and techniques that were developed for some particular
case may in fact apply much more generally.


  Hardy and Higgins first proposed the amyloid cascade hypothesis in 1992,
stating that the decrease in neuronal function observed in Alzheimer's Disease
(AD) is due to a process initiated by the oligomerization of amyloid-beta
peptides. One hypothesis states that toxicity arises from the aggregation of
amyloid-beta into a pore structure, which can then puncture the brain cell
membrane; this allow toxic calcium ions to flood through the opening, causing
eventual cell death. In 2007, neurobiologist Ruth Nussinov calculated the three
pore sizes most likely to occur within the brain. Based on her findings, we
constructed a method to determine the time it takes for a cell to die after the
cell is punctured by the pore. Our findings have shown that cell death occurs
within one second after the oligomer makes contact with the cell. We believe
this is important because instant cell death has been one criticism of
Nussinov's model, and we have calculated a concrete time value for that
criticism. We identify two potential deficiencies with our model that could be
improved: first, we treat the calcium as an ideal gas, which it is not; second,
we assume that the pores are static (i.e. constantly open), while recent
developments suggest that they may open and close dynamically.


  We investigate the conformations of DNA-like stiff chains, characterized by
contour length ($L$) and persistence length ($l_p$), in a variety of crowded
environments containing mono disperse soft spherical (SS) and spherocylindrical
(SC) particles, mixture of SS and SC, and a milieu mimicking the composition of
proteins in $E. coli.$ cytoplasm. The stiff chain, whose size modestly
increases in SS crowders up to $\phi\approx 0.1$, is considerably more compact
at low volume fractions ($\phi \leq 0.2$) in monodisperse SC particles than in
a medium containing SS particles. A 1:1 mixture of SS and SC crowders induces
greater chain compaction than the pure SS or SC crowders at the same $\phi$
with the effect being highly non-additive. We also discover a counter-intuitive
result that polydisperse crowding environment, mimicking the composition of a
cell lysate, swells the DNA-like polymer, which is in stark contrast to the
size reduction of flexible polymer in the same milieu. Trapping of the stiff
chain in a fluctuating tube-like environment created by large-sized crowders
explains the dramatic increase in size and persistence length of the stiff
chain. In the polydisperse medium, mimicking the cellular environment, the size
of the DNA (or related RNA) is determined by $L/l_p$. At low $L/l_p$ the size
of the polymer is unaffected whereas there is a dramatic swelling at
intermediate value of $L/l_p$. We use these results to provide insights into
recent experiments on crowding effects on RNA, and also make testable
predictions.


  I consider the many ways in which evolved information-flows are restricted
and metabolic resources protected and hidden -- the thesis of living phenomena
as evolutionary cryptosystems. I present the information theory of secrecy
systems and discuss mechanisms acquired by evolved lineages that encrypt
sensitive heritable information with random keys. I explore the idea that
complexity science is a cryptographic discipline as "frozen accidents", or
various forms of regularized randomness, historically encrypt adaptive
dynamics.


  The Lenski experiment investigates the long-term evolution of bacterial
populations. Its design allows the direct comparison of the reproductive
fitness of an evolved strain with its founder ancestor. It was observed by
Wiser et al. (2013) that the relative fitness over time increases sublinearly,
a behaviour which is commonly attributed to effects like clonal interference or
epistasis. In this paper we present an individual-based probabilistic model
that captures essential features of the design of the Lenski experiment. We
assume that each beneficial mutation increases the individual reproduction rate
by a fixed amount, which corresponds to the absence of epistasis in the
continuous-time (intraday) part of the model, but leads to an epistatic effect
in the discrete-time (interday) part of the model. Using an approximation by
near-critical Galton-Watson processes, we prove that under some assumptions on
the model parameters which exclude clonal interference, the relative fitness
process converges, after suitable rescaling, in the large population limit to a
power law function.


  We propose the entropy of random Markov trajectories originating and
terminating at a state as a measure of the stability of a state of a Markov
process. These entropies can be computed in terms of the entropy rates and
stationary distributions of Markov processes. We apply this definition of
stability to local maxima and minima of the stationary distribution of the
Moran process with mutation and show that variations in population size,
mutation rate, and strength of selection all affect the stability of the
stationary extrema.


  When double-stranded DNA molecules are heated, or exposed to denaturing
agents, the two strands get separated. The statistical physics of this process
has a long history, and is commonly described in term of the Poland-Scheraga
(PS) model. Crucial to this model is the configurational entropy for a melted
region (compared to the entropy of an intact region of the same size),
quantified by the loop factor. In this study we investigate how confinement
affects the DNA melting transition, by using the loop factor for an ideal
Gaussian chain. By subsequent numerical solutions of the PS model, we
demonstrate that the melting temperature depends on the persistence lengths of
single-stranded and double-stranded DNA. For realistic values of the
persistence lengths the melting temperature is predicted to decrease with
decreasing channel diameter. We also demonstrate that confinement broadens the
melting transition. These general findings hold for the three scenarios
investigated: namely 1. homo-DNA, i.e. identical basepairs along the DNA
molecule; 2. random sequence DNA, and 3. "real" DNA, here T4 phage DNA. We show
that cases 2 and 3 in general give rise to broader transitions than case 1.
Case 3 exhibits a similar phase transition as case 2 provided the random
sequence DNA has the same ratio of AT to GC basepairs. A simple analytical
estimate for the shift in melting temperature is provided as a function of
nanochannel diameter. For homo-DNA, we also present an analytical prediction of
the melting probability as a function of temperature.


  The editing of transcribed RNA by other molecules such that the form of the
final product differs from that specified in the corresponding DNA sequence is
ubiquitous. This paper uses an abstract, tunable Boolean genetic regulatory
network model to explore aspects of RNA editing. In particular, it is shown how
dynamically altering expressed sequences via a guide RNA-inspired mechanism can
be selected for by simulated evolution under various single and multicellular
scenarios.


  The Calculus of Wrapped Compartments is a framework based on stochastic
multiset rewriting in a compartmentalised setting originally developed for the
modelling and analysis of biological interactions. In this paper, we propose to
use this calculus for the description of ecological systems and we provide the
modelling guidelines to encode within the calculus some of the main
interactions leading ecosystems evolution. As a case study, we model the
distribution of height of Croton wagneri, a shrub constituting the endemic
predominant species of the dry ecosystem in southern Ecuador. In particular, we
consider the plant at different altitude gradients (i.e. at different
temperature conditions), to study how it adapts under the effects of global
climate change.


  An RNA-seq experiment with 48 biological replicates in each of 2 conditions
was performed to determine the number of biological replicates ($n_r$)
required, and to identify the most effective statistical analysis tools for
identifying differential gene expression (DGE). When $n_r=3$, seven of the nine
tools evaluated give true positive rates (TPR) of only 20 to 40 percent. For
high fold-change genes ($|log_{2}(FC)|\gt2$) the TPR is $\gt85$ percent. Two
tools performed poorly; over- or under-predicting the number of differentially
expressed genes. Increasing replication gives a large increase in TPR when
considering all DE genes but only a small increase for high fold-change genes.
Achieving a TPR $\gt85$% across all fold-changes requires $n_r\gt20$. For
future RNA-seq experiments these results suggest $n_r\gt6$, rising to
$n_r\gt12$ when identifying DGE irrespective of fold-change is important. For
$6 \lt n_r \lt 12$, superior TPR makes edgeR the leading tool tested. For $n_r
\ge12$, minimizing false positives is more important and DESeq outperforms the
other tools.


  Malaria is a major contributor to health burdens throughout the regions where
it is endemic. Historically, it was believed that there was limited morbidity
and essentially no mortality associated with Plasmodium vivax; however,
evidence from diverse settings now suggests that infections with P. vivax can
be both severe and fatal. This awareness has highlighted a critical gap: the
vast majority of research has been directed towards P. falciparum, leading to a
decades-long neglect of epidemiological and clinical studies of P. vivax. There
exists a large body of historical data on human experimental infections with P.
vivax; these studies in controlled settings provided a wealth of wide-ranging
statements based on expert opinion, which form the basis for much of what is
currently known about P. vivax. In this thesis, portions of this evidence-base
have been re-examined using modern epidemiological analyses with two aims: to
critically examine this accumulated knowledge base, and to inform current
research agendas towards global malaria elimination for all species of
Plasmodium. Chapter 2 examines geographic variation in the epidemiology of P.
vivax, especially the timing of incubation periods and of relapses, by origin
of the parasites. Chapter 3 re-assesses the impact of sporozoite dosage upon
incubation and pre-patent periods; Chapter 4 provides well-defined mathematical
distributions for incubation and relapses periods in experimental infections,
and explores their epidemiological impacts using simple transmission models.
Chapter 5 examines the epidemiology of mixed-strain P. vivax infections and
compares these results with studies in murine models and general ecological
theory; and Chapter 6 clarifies the origin of the Madagascar strain of P.
vivax.


  The Poisson Nernst-Planck equations for charge concentration and electric
potential in a ball is a model of electro-diffusion of ions in the head of a
neuronal dendritic spine. We study the relaxation and the steady state when an
initial charge of ions is injected into the ball. The steady state equation is
similar to the Liouville-Gelfand-Bratu-type equation with the difference that
the boundary condition is Neumann, not Dirichlet and there a minus sign in the
exponent of the exponential term. The entire boundary is impermeable to the
ions and the electric field satisfies the compatibility condition of Poisson's
equation. We construct a steady radial solution and find that the potential is
maximal in the center and decreases toward the boundary. We study the limit of
large charge in dimension 1,2 and 3. For the case of a small absorbing window
in the sphere, we find the escape rate of an ion from the steady density.


  Starting with Darwin, biologists have asked how populations evolve from a low
fitness state that is evolutionarily stable to a high fitness state that is
not. Specifically of interest is the emergence of cooperation and
multicellularity where the fitness of individuals often appears in conflict
with that of the population. Theories of social evolution and evolutionary game
theory have produced a number of fruitful results employing two-state two-body
frameworks. In this study we depart from this tradition and instead consider a
multi-player, multi-state evolutionary game, in which the fitness of an agent
is determined by its relationship to an arbitrary number of other agents. We
show that populations organize themselves in one of four distinct phases of
interdependence depending on one parameter, selection strength. Some of these
phases involve the formation of specialized large-scale structures. We then
describe how the evolution of independence can be manipulated through various
external perturbations.


  Intracellular biochemical networks fluctuate dynamically due to various
internal and external sources of fluctuation. Dissecting the fluctuation into
biologically relevant components is important for understanding how a cell
controls and harnesses noise and how information is transferred over apparently
noisy intracellular networks. While substantial theoretical and experimental
advancement on the decomposition of fluctuation was achieved for feedforward
networks without any loop, we still lack a theoretical basis that can
consistently extend such advancement to feedback networks. The main obstacle
that hampers is the circulative propagation of fluctuation by feedback loops.
In order to define the relevant quantity for the impact of feedback loops for
fluctuation, disentanglement of the causally interlocked influence between the
components is required. In addition, we also lack an approach that enables us
to infer non-perturbatively the influence of the feedback to fluctuation as the
dual reporter system does in the feedforward network. In this work, we resolve
these problems by extending the work on the fluctuation decomposition and the
dual reporter system. For a single-loop feedback network with two components,
we define feedback loop gain as the feedback efficiency that is consistent with
the fluctuation decomposition for feedforward networks. Then, we clarify the
relation of the feedback efficiency with the fluctuation propagation in an
open-looped FF network. Finally, by extending the dual reporter system, we
propose a conjugate feedback and feedforward system for estimating the feedback
efficiency only from the statistics of the system non-perturbatively.


  Phenotype-switching with and without sensing environment is a ubiquitous
strategy of organisms to survive in fluctuating environment. Fitness of a
population of organisms with phenotype-switching may be constrained and
restricted by hidden relations as the entropy production in a thermal system
with and without sensing and feedback is well-characterized via fluctuation
relations (FRs) . In this work, we derive such FRs of fitness together with an
underlying information-theoretic structure in selection. By using path-integral
formulation of a multi-phenotype population dynamics, we clarify that the
optimal switching strategy is characterized as a consistency condition for
time-forward and backward path probabilities. Within the formulation, the
selection is regarded as passive information compression, and the loss of
fitness from the optimal strategy is shown to satisfy various FRs that
constrain the average and fluctuation of the loss. These results are naturally
extended to the situation that organisms can use an environmental signal by
actively sensing the environment. FRs of fitness gain by sensing are derived in
which the multivariate mutual information among the phenotype, the environment
and the signal plays the role to quantify the relevant information in the
signal for fitness gain.


  A complexity-theoretic approach to studying biological networks is proposed.
A simple graph representation is used where molecules (DNA, RNA, proteins and
chemicals) are vertices and relations between them are directed and signed
(promotional (+) or inhibitory (-)) edges. Based on this model, the problem of
network evolution (NE) is defined formally as an optimization problem and
subsequently proven to be fundamentally hard (NP-hard) by means of reduction
from the Knapsack problem (KP). Second, for empirical validation, various
biological networks of experimentally-validated interactions are compared
against randomly generated networks with varying degree distributions. An NE
instance is created using a given real or synthetic (random) network. After
being reverse-reduced to a KP instance, each NE instance is fed to a KP solver
and the average achieved knapsack value-to-weight ratio is recorded from
multiple rounds of simulated evolutionary pressure. The results show that
biological networks (and synthetic networks of similar degree distribution)
achieve the highest ratios at maximal evolutionary pressure and minimal error
tolerance conditions. The more distant (in degree distribution) a synthetic
network is from biological networks the lower its achieved ratio. The results
shed light on how computational intractability has shaped the evolution of
biological networks into their current topology.


  We study an elasticity model for compressed protein monolayers or particle
rafts at a liquid interface. Based on the microscopic view of hard-core
particles with soft shells, a bead-spring model is formulated and analyzed in
terms of continuum elasticity theory. The theory can be applied, for example,
to hydrophobin-coated air-water interfaces or, more generally, to liquid
interfaces coated with an adsorbed monolayer of interacting hard-core
particles. We derive constitutive relations for such particle rafts and
describe the buckling of compressed planar liquid interfaces as well as their
apparent Poisson ratio. We also use the constitutive relations to obtain shape
equations for pendant or buoyant capsules attached to a capillary, and to
compute deflated shapes of such capsules. A comparison with capsules obeying
the usual Hookean elasticity (without hard cores) reveals that the hard cores
trigger capsule wrinkling. Furthermore, it is shown that a shape analysis of
deflated capsules with hard-core/soft-shell elasticity gives apparent elastic
moduli which can be much higher than the original values if Hookean elasticity
is assumed.


  We consider a model for Darwinian evolution in an asexual population with a
large but non-constant populations size characterized by a natural birth rate,
a logistic death rate modelling competition and a probability of mutation at
each birth event. In the present paper, we study the long-term behavior of the
system in the limit of large population $(K\to \infty)$ size, rare mutations
$(u\to 0)$, and small mutational effects $(\sigma\to 0)$, proving convergence
to the canonical equation of adaptive dynamics (CEAD). In contrast to earlier
works, e.g. by Champagnat and M\'el\'eard, we take the three limits
simultaneously, i.e. $u=u_K$ and $\sigma=\sigma_K$, tend to zero with $K$,
subject to conditions that ensure that the time-scale of birth and death events
remains separated from that of successful mutational events. This slows down
the dynamics of the microscopic system and leads to serious technical
difficulties that requires the use of completely different methods. In
particular, we cannot use the law of large numbers on the diverging time needed
for fixation to approximate the stochastic system with the corresponding
deterministic one. To solve this problem we develop a "stochastic Euler scheme"
based on coupling arguments that allows to control the time evolution of the
stochastic system over time-scales that diverge with $K$.


  A central goal of synthetic biology is to design sophisticated synthetic
cellular circuits that can perform complex computations and information
processing tasks in response to specific inputs. The tremendous advances in our
ability to understand and manipulate cellular information processing networks
raises several fundamental physics questions: How do the molecular components
of cellular circuits exploit energy consumption to improve information
processing? Can one utilize ideas from thermodynamics to improve the design of
synthetic cellular circuits and modules? Here, we summarize recent theoretical
work addressing these questions. Energy consumption in cellular circuits serves
five basic purposes: (1) increasing specificity, (2) manipulating dynamics, (3)
reducing variability, (4) amplifying signal, and (5) erasing memory. We
demonstrate these ideas using several simple examples and discuss the
implications of these theoretical ideas for the emerging field of synthetic
biology. We conclude by discussing how it may be possible to overcome these
limitations using "post-translational" synthetic biology that exploits
reversible protein modification.


  Over the past several decades there has been a proliferation of
epidemiological models with ordinary derivatives replaced by fractional
derivatives in an an-hoc manner. These models may be mathematically interesting
but their relevance is uncertain. Here we develop an SIR model for an epidemic,
including vital dynamics, from an underlying stochastic process. We show how
fractional differential operators arise naturally in these models whenever the
recovery time from the disease is power law distributed. This can provide a
model for a chronic disease process where individuals who are infected for a
long time are unlikely to recover. The fractional order recovery model is shown
to be consistent with the Kermack-McKendrick age-structured SIR model and it
reduces to the Hethcote-Tudor integral equation SIR model. The derivation from
a stochastic process is extended to discrete time, providing a stable numerical
method for solving the model equations. We have carried out simulations of the
fractional order recovery model showing convergence to equilibrium states. The
number of infecteds in the endemic equilibrium state increases as the
fractional order of the derivative tends to zero.


  The identification of low-energy conformers for a given molecule is a
fundamental problem in computational chemistry and cheminformatics. We assess
here a conformer search that employs a genetic algorithm for sampling the
low-energy segment of the conformation space of molecules. The algorithm is
designed to work with first-principles methods, facilitated by the
incorporation of local optimization and blacklisting conformers to prevent
repeated evaluations of very similar solutions. The aim of the search is not
only to find the global minimum, but to predict all conformers within an energy
window above the global minimum. The performance of the search strategy is: (i)
evaluated for a reference data set extracted from a database with amino acid
dipeptide conformers obtained by an extensive combined force field and
first-principles search and (ii) compared to the performance of a systematic
search and a random conformer generator for the example of a drug-like ligand
with 43 atoms, 8 rotatable bonds and 1 cis/trans bond.


  The antibiotics time machine is an optimization question posed by Mira
\latin{et al.} on the design of antibiotic treatment plans to minimize
antibiotic resistance. The problem is a variation of the Markov decision
process. These authors asked if the problem can be solved efficiently. In this
paper, we show that this problem is NP-hard in general.


  We present a novel approach to RNA-Seq quantification that is near optimal in
speed and accuracy. Software implementing the approach, called kallisto, can be
used to analyze 30 million unaligned paired-end RNA-Seq reads in less than 5
minutes on a standard laptop computer while providing results as accurate as
those of the best existing tools. This removes a major computational bottleneck
in RNA-Seq analysis.


  The transition to turbulence exhibits remarkable spatio-temporal behavior
that continues to defy detailed understanding. Near the onset to turbulence in
pipes, transient turbulent regions decay either directly or, at higher Reynolds
numbers through splitting, with characteristic time-scales that exhibit a
super-exponential dependence on Reynolds number. Here we report numerical
simulations of transitional pipe flow, showing that a zonal flow emerges at
large scales, activated by anisotropic turbulent fluctuations; in turn, the
zonal flow suppresses the small-scale turbulence leading to stochastic
predator-prey dynamics. We show that this "ecological" model of transitional
turbulence reproduces the super-exponential lifetime statistics and
phenomenology of pipe flow experiments. Our work demonstrates that a fluid on
the edge of turbulence is mathematically analogous to an ecosystem on the edge
of extinction, and provides an unbroken link between the equations of fluid
dynamics and the directed percolation universality class.


  Neutral dynamics, where taxa are assumed to be demographically equivalent and
their abundance is governed solely by the stochasticity of the underlying
birth-death process, has proved itself as an important minimal model that
accounts for many empirical datasets in genetics and ecology. However, the
restriction of the model to demographic [${\cal{O}} ({\sqrt N})$] noise yields
relatively slow dynamics that appears to be in conflict with both short-term
and long-term characteristics of the observed systems. Here we analyze two of
these problems - age size relationships and species extinction time - in the
framework of a neutral theory with both demographic and environmental
stochasticity. It turns out that environmentally induced variations of the
demographic rates control the long-term dynamics and modify dramatically the
predictions of the neutral theory with demographic noise only, yielding much
better agreement with empirical data. We consider two prototypes of "zero mean"
environmental noise, one which is balanced with regard to the arithmetic
abundance, another balanced in the logarithmic (fitness) space, study their
species lifetime statistics and discuss their relevance to realistic models of
community dynamics.


  There has been a long-standing and at times fractious debate whether complex
and large systems can be stable. In ecology, the so-called `diversity-stability
debate' arose because mathematical analyses of ecosystem stability were either
specific to a particular model (leading to results that were not general), or
chosen for mathematical convenience, yielding results unlikely to be meaningful
for any interesting realistic system. May's work, and its subsequent
elaborations, relied upon results from random matrix theory, particularly the
circular law and its extensions, which only apply when the strengths of
interactions between entities in the system are assumed to be independent and
identically distributed (i.i.d.). Other studies have optimistically generalised
from the analysis of very specific systems, in a way that does not hold up to
closer scrutiny. We show here that this debate can be put to rest, once these
two contrasting views have been reconciled --- which is possible in the
statistical framework developed here. Here we use a range of illustrative
examples of dynamical systems to demonstrate that (i) stability probability
cannot be summarily deduced from any single property of the system (e.g. its
diversity), and (ii) our assessment of stability depends on adequately
capturing the details of the systems analysed. Failing to condition on the
structure of dynamical systems will skew our analysis and can, even for very
small systems, result in an unnecessarily pessimistic diagnosis of their
stability.


  The behaviour of a high dimensional stochastic system described by a Chemical
Master Equation (CME) depends on many parameters, rendering explicit simulation
an inefficient method for exploring the properties of such models. Capturing
their behaviour by low-dimensional models makes analysis of system behaviour
tractable. In this paper, we present low dimensional models for the
noise-induced excitable dynamics in Bacillus subtilis, whereby a key protein
ComK, which drives a complex chain of reactions leading to bacterial
competence, gets expressed rapidly in large quantities (competent state) before
subsiding to low levels of expression (vegetative state). These rapid reactions
suggest the application of an adiabatic approximation of the dynamics of the
regulatory model that, however, lead to competence durations that are incorrect
by a factor of 2. We apply a modified version of an iterative functional
procedure that faithfully approximates the time-course of the trajectories in
terms of a 2-dimensional model involving proteins ComK and ComS. Furthermore,
in order to describe the bimodal bivariate marginal probability distribution
obtained from the Gillespie simulations of the CME, we introduce a tunable
multiplicative noise term in a 2-dimensional Langevin model whose stationary
state is described by the time-independent solution of the corresponding
Fokker-Planck equation.


  Globally disruptive events include asteroid/comet impacts, large igneous
provinces and glaciations, all of which have been considered as contributors to
mass extinctions. Understanding the overall relationship between the timings of
the largest extinctions and their potential proximal causes remains one of
science's great unsolved mysteries. Cycles of about 60 million years in both
fossil diversity and environmental data suggest external drivers such as the
passage of the Solar System through the galactic plane. While cyclic phenomena
are recognised statistically, a lack of coherent mechanisms and a failure to
link key events has hampered wider acceptance of multi-million year periodicity
and its relevance to earth science and evolution. The generation of a robust
predictive model of timings, with a clear plausible primary mechanism, would
signal a paradigm shift. Here, we present a model of the timings of globally
disruptive events and a possible explanation of their ultimate cause. The
proposed model is a symmetrical pattern of 63 million-year sequences around a
central value, interpreted as the occurrence of events along, and parallel to,
the galactic midplane. The symmetry is consistent with multiple dark matter
disks, aligned parallel to the midplane. One implication of the precise pattern
of timings and the underlying physical model is the ability to predict future
events, such as a major extinction in one to two million years.


  Some cells have to take decision based on the quality of surroundings
ligands, almost irrespective of their quantity, a problem we name "absolute
discrimination". An example of absolute discrimination is recognition of
not-self by immune T Cells. We show how the problem of absolute discrimination
can be solved by a process called "adaptive sorting". We review several
implementations of adaptive sorting, as well as its generic properties such as
antagonism. We show how kinetic proofreading with negative feedback implements
an approximate version of adaptive sorting in the immune context. Finally, we
revisit the decision problem at the cell population level, showing how
phenotypic variability and feedbacks between population and single cells are
crucial for proper decision.


  Human {\alpha}-synuclein (SNCA) is a 140-amino-acid protein belonging to the
three-member synuclein family. It has been extensively studied due to its
misfolding/aggregation in and genetic linkage to neurodegenerative diseases,
especially Parkinson's disease (PD). To better understand its biology, models
of SNCA toxicity have been developed in budding yeast over the past decade,
which have yielded insights into the protein's modes of action in specific
pathways and potential therapeutic targets. Given that the synuclein gene
family is not present in yeast, an extensive homology search was undertaken to
determine if any yeast protein may possess structural homology to SNCA and
whose native biology may shed more light on SNCA's pathomechanism in
eukaryotes. We identified Zeo1, a membrane-associated protein involved in the
cell wall integrity (CWI) pathway, as a candidate structural homolog. We show
that Zeo1 overexpression is toxic in yeast and, similar to SNCA, localizes to
lipid membranes. A number of biochemical similarities between Zeo1 and SNCA
also become apparent in light of this potential structural connection.
Moreover, the yeast PKC1 gene, a kinase acting as a downstream signaling hub in
the CWI pathway, rescues both SNCA- and Zeo1-induced toxicities. Using the same
homology search methods that identified Zeo1, we show that Pkc1 has a hybrid
structural similarity to PINK1 and PARIS, two mitochondrial PD-implicated
proteins not generally linked directly to synuclein-specific pathobiology.
Overall, this proof-of-concept study shows the potential utility of hitherto
uncharacterized cross-species structural homologs, identified using comparative
proteome-wide structure prediction algorithms, in shedding light on abstruse
connections among disease-relevant proteins.


  Quantitative methods for studying biodiversity have been traditionally rooted
in the classical theory of finite frequency tables analysis. However, with the
help of modern experimental tools, like high throughput sequencing, we now
begin to unlock the outstanding diversity of genomic data in plants and animals
reflective of the long evolutionary history of our planet. This molecular data
often defies the classical frequency/contingency tables assumptions and seems
to require sparse tables with very large number of categories and highly
unbalanced cell counts, e.g., following heavy tailed distributions (for
instance, power laws). Motivated by the molecular diversity studies, we propose
here a frequency-based framework for biodiversity analysis in the asymptotic
regime where the number of categories grows with sample size (an infinite
contingency table). Our approach is rooted in information theory and based on
the Gaussian limit results for the effective number of species (the Hill
numbers) and the empirical Renyi entropy and divergence. We argue that when
applied to molecular biodiversity analysis our methods can properly account for
the complicated data frequency patterns on one hand and the practical sample
size limitations on the other. We illustrate this principle with two specific
RNA sequencing examples: a comparative study of T-cell receptor populations and
a validation of some preselected molecular hepatocellular carcinoma (HCC)
markers.


  Moment-closure techniques are commonly used to generate low-dimensional
deterministic models to approximate the average dynamics of stochastic systems
on networks. The quality of such closures is usually difficult to asses and the
relationship between model assumptions and closure accuracy are often
difficult, if not impossible, to quantify. Here we carefully examine some
commonly used moment closures, in particular a new one based on the concept of
maximum entropy, for approximating the spread of epidemics on networks by
reconstructing the probability distributions over triplets based on those over
pairs. We consider various models (SI, SIR, SEIR and Reed-Frost-type) under
Markovian and non-Markovian assumption characterising the latent and infectious
periods. We initially study two special networks, namely the open triplet and
closed triangle, for which we can obtain analytical results. We then explore
numerically the exactness of moment closures for a wide range of larger motifs,
thus gaining understanding of the factors that introduce errors in the
approximations, in particular the presence of a random duration of the
infectious period and the presence of overlapping triangles in a network. We
also derive a simpler and more intuitive proof than previously available
concerning the known result that pair-based moment closure is exact for the
Markovian SIR model on tree-like networks under pure initial conditions. We
also extend such a result to all infectious models, Markovian and
non-Markovian, in which susceptibles escape infection independently from each
infected neighbour and for which infectives cannot regain susceptible status,
provided the network is tree-like and initial conditions are pure. This works
represent a valuable step in deepening understanding of the assumptions behind
moment closure approximations and for putting them on a more rigorous
mathematical footing.


  The abundance of a species' population in an ecosystem is rarely stationary,
often exhibiting large fluctuations over time. Using historical data on marine
species, we show that the year-to-year fluctuations of population growth rate
obey a well-defined double-exponential (Laplace) distribution. This striking
regularity allows us to devise a stochastic model despite seemingly irregular
variations in population abundances. The model identifies the effect of reduced
growth at low population density as a key factor missed in current approaches
of population variability analysis and without which extinction risks are
severely underestimated. The model also allows us to separate the effect of
demographic stochasticity and show that single-species growth rates are
dominantly determined by stochasticity common to all species. This
dominance---and the implications it has for interspecies correlations,
including co-extinctions---emphasizes the need of ecosystem-level management
approaches to reduce the extinction risk of the individual species themselves.


  Cooperative behavior is widespread in nature, even though cooperating
individuals always run the risk to be exploited by free-riders. Population
structure effectively promotes cooperation given that a threshold in the level
of cooperation was already reached. However, the question how cooperation can
emerge from a single mutant, which cannot rely on a benefit provided by other
cooperators, is still puzzling. Here, we investigate this question for a
well-defined but generic situation based on typical life-cycles of microbial
populations where individuals regularly form new colonies followed by growth
phases. We analyze two evolutionary mechanisms favoring cooperative behavior
and study their strength depending on the inoculation size and the length of a
life-cycle. In particular, we find that population bottlenecks followed by
exponential growth phases strongly increase the survival and fixation
probabilities of a single cooperator in a free-riding population.


  The viscous liquid surrounding a hair bundle dissipates energy and dampens
oscillations, which poses a fundamental physical challenge to the high
sensitivity and sharp frequency selectivity of hearing. To identify the
mechanical forces at play, we constructed a detailed finite-element model of
the hair bundle. Based on data from the hair bundle of the bullfrog's sacculus,
this model treats the interaction of stereocilia both with the surrounding
liquid and with the liquid in the narrow gaps between the individual
stereocilia. The investigation revealed that grouping stereocilia in a bundle
dramatically reduces the total drag. During hair-bundle deflections, the tip
links potentially induce drag by causing small but very dissipative relative
motions between stereocilia; this effect is offset by the horizontal top
connectors that restrain such relative movements at low frequencies. For higher
frequencies the coupling liquid is sufficient to assure that the hair bundle
moves as a unit with a low total drag. This work reveals the mechanical
characteristics originating from hair-bundle morphology and shows
quantitatively how a hair bundle is adapted for sensitive mechanotransduction.


  DNA strand displacement systems have proven themselves to be fertile
substrates for the design of programmable molecular machinery and circuitry.
Domain-level reaction enumerators provide the foundations for molecular
programming languages by formalizing DNA strand displacement mechanisms and
modeling interactions at the "domain" level - one level of abstraction above
models that explicitly describe DNA strand sequences. Unfortunately, the
most-developed models currently only treat pseudo-linear DNA structures, while
many systems being experimentally and theoretically pursued exploit a much
broader range of secondary structure configurations. Here, we describe a new
domain-level reaction enumerator that can handle arbitrary non-pseudoknotted
secondary structures and reaction mechanisms including association and
dissociation, 3-way and 4-way branch migration, and direct as well as remote
toehold activation. To avoid polymerization that is inherent when considering
general structures, we employ a time-scale separation technique that holds in
the limit of low concentrations. This also allows us to "condense" the detailed
reactions by eliminating fast transients, with provable guarantees of
correctness for the set of reactions and their kinetics. We hope that the new
reaction enumerator will be used in new molecular programming languages,
compilers, and tools for analysis and verification that treat a wider variety
of mechanisms of interest to experimental and theoretical work. We have
implemented this enumerator in Python, and it is included in the DyNAMiC
Workbench Integrated Development Environment.


  We propose a honeybee-mite-virus model that incorporates (1) parasitic
interactions between honeybees and the Varroa mites; (2) five virus
transmission terms between honeybees and mites at different stages of Varroa
mites: from honeybees to honeybees, from adult honeybees to phoretic mites,
from honeybee brood to reproductive mites, from reproductive mites to honeybee
brood, and from honeybees to phoretic mites; and (3) Allee effects in the
honeybee population generated by its internal organization such as division of
labor. We provide completed local and global analysis for the full system and
its subsystems. Our analytical and numerical results allow us have a better
understanding of the synergistic effects of parasitism and virus infections on
honeybee population dynamics and its persistence. Interesting findings from our
work include: (a) Due to Allee effects experienced by the honeybee population,
initial conditions are essential for the survival of the colony. (b) Low adult
honeybee to brood ratios have destabilizing effects on the system, generate
fluctuated dynamics, and potentially lead to a \emph{catastrophic event} where
both honeybees and mites suddenly become extinct. This catastrophic event could
be potentially linked to Colony Collapse Disorder (CCD) of honeybee colonies.
(c) Virus infections may have stabilizing effects on the system, and could make
disease more persistent in the presence of parasitic mites. Our model
illustrates how the synergy between the parasitic mites and virus infections
consequently generates rich dynamics including multiple attractors where all
species can coexist or go extinct depending on initial conditions. Our findings
may provide important insights on honeybee diseases and parasites and how to
best control them.


  Foraging movements of predator play an important role in population dynamics
of prey-predator interactions, which have been considered as mechanisms that
contribute to spatial self-organization of prey and predator. In nature, there
are many examples of prey-predator interactions where prey is immobile while
predator disperses between patches non-randomly through different factors such
as stimuli following the encounter of a prey. In this work, we formulate a
Rosenzweig-MacArthur prey-predator two patch model with mobility only in
predator and the assumption that predators move towards patches with more
concentrated prey-predator interactions. We provide completed local and global
analysis of our model. Our analytical results combined with bifurcation
diagrams suggest that: (1) dispersal may stabilize or destabilize the coupled
system; (2) dispersal may generate multiple interior equilibria that lead to
rich bistable dynamics or may destroy interior equilibria that lead to the
extinction of predator in one patch or both patches; (3) Under certain
conditions, the large dispersal can promote the permanence of the system. In
addition, we compare the dynamics of our model to the classic two patch model
to obtain a better understanding how different dispersal strategies may have
different impacts on the dynamics and spatial patterns.


  Contact tracing is an important control strategy for containing Ebola
epidemics. From a theoretical perspective, explicitly incorporating contact
tracing with disease dynamics presents challenges, and population level effects
of contact tracing are difficult to determine. In this work, we formulate and
analyze a mechanistic SEIR type outbreak model which considers the key features
of contact tracing, and we characterize the impact of contact tracing on the
effective reproduction number, $\mathcal R_e$, of Ebola. In particular, we
determine how relevant epidemiological properties such as incubation period,
infectious period and case reporting, along with varying monitoring protocols,
affect the efficacy of contact tracing. In the special cases of either perfect
monitoring of traced cases or perfect reporting of all cases, we derive simple
formulae for the critical proportion of contacts that need to be traced in
order to bring the effective reproduction number $\mathcal R_e$ below one.
Also, in either case, we show that $\mathcal R_e$ can be expressed completely
in terms of observable reported case/tracing quantities, namely $\mathcal
R_e=k\dfrac{(1-q)}{q}+k_m$ where $k$ is the number of secondary traced infected
contacts per primary untraced reported case, $k_m$ is the number of secondary
traced infected contacts per primary traced reported case and $(1-q)/q$ is the
odds that a reported case is not a traced contact. These formulae quantify
contact tracing as both an intervention strategy that impacts disease spread
and a probe into the current epidemic status at the population level. Data from
the West Africa Ebola outbreak is utilized to form real-time estimates of
$\mathcal R_e$, and inform our projections of the impact of contact tracing,
and other control measures, on the epidemic trajectory.


  We show that a bacteria and bacteriophage system with either a perfectly
nested or a one-to-one infection network is permanent, a.k.a uniformly
persistent, provided that bacteria that are superior competitors for nutrient
devote the least to defence against infection and the virus that are the most
efficient at infecting host have the smallest host range. By ensuring that the
density-dependent reduction in bacterial growth rates are independent of
bacterial strain, we are able to arrive at the permanence conclusion sought by
Jover et al (J. Theor. Biol. 332:65-77, 2013). The same permanence results hold
for the one-to-one infection network considered by Thingstad (Limnol Oceanogr
45:1320-1328, 2000) but without virus efficiency ordering. Additionally we show
the global stability for the nested infection network, and the global dynamics
for the one-to-one network.


  As a simple model for criminal behavior, the traditional two-strategy
inspection game yields counterintuitive results that fail to describe empirical
data. The latter shows that crime is often recurrent, and that crime rates do
not respond linearly to mitigation attempts. A more apt model entails ordinary
people who neither commit nor sanction crime as the third strategy besides the
criminals and punishers. Since ordinary people free-ride on the sanctioning
efforts of punishers, they may introduce cyclic dominance that enables the
coexistence of all three competing strategies. In this setup ordinary
individuals become the biggest impediment to crime abatement. We therefore also
consider heterogeneous punisher strategies, which seek to reduce their
investment into fighting crime in order to attain a more competitive payoff. We
show that this diversity of punishment leads to an explosion of complexity in
the system, where the benefits and pitfalls of criminal behavior are revealed
in the most unexpected ways. Due to the raise and fall of different alliances
no less than six consecutive phase transitions occur in dependence on solely
the temptation to succumb to criminal behavior, leading the population from
ordinary people-dominated across punisher-dominated to crime-dominated phases,
yet always failing to abolish crime completely.


  Cellular reprogramming, the conversion of one cell type to another, has
fundamentally transformed our conception of cell types. Cellular reprogramming
induces global changes in gene expression involving hundreds of transcription
factors and thousands of genes and understanding how cells globally alter their
gene expression profile during reprogramming is an open problem. Here we
reanalyze time-series data on cellular reprogramming from differentiated cell
types to induced pluripotent stem cells (iPSCs) and show that gene expression
dynamics during reprogramming follow a simple one-dimensional reaction
coordinate. This reaction coordinate is independent of both the time it takes
to reach the iPSC state as well as the details of experimental protocol used.
Using Monte-Carlo simulations, we show that such a reaction coordinate emerges
naturally from epigenetic landscape models of cell identity where cellular
reprogramming is viewed as a "barrier-crossing" between the starting and ending
cell fates. The model also provides gene-level insight into reprogramming
dynamics and resolves a debate in the stem cell field about the different
phases of reprogramming dynamics. Overall, our analysis and model suggest that
gene expression dynamics during reprogramming follow a canonical trajectory
consistent with the idea of an "optimal path" in gene expression space for
reprogramming.


  In Misrepair-accumulation theory, we have proposed a Misrepair mechanism for
interpreting aging. However, Misrepair mechanism is also important in
biological adaptation. Misrepair is a strategy of repair for increasing the
surviving chance of an organism when it suffers from severe injuries. As a
surviving strategy, Misrepair mechanism plays also an important role in
individual adaptation, species survival, and species differentiation. Firstly,
Misrepair of an injury is one of the manners of individual adaptation; and
Misrepair mechanism gives an organism a great potential in adapting to
changeable and destructive environment. Secondly, Misrepair mechanism is
important in maintaining and enlarging the diversity of genome DNAs of a
species; and a large diversity of genome DNAs is essential for the adaptation
and the differentiation of a species in different environments. On one hand,
somatic Misrepairs are essential for maintaining the sufficient number of
individuals in a species, which are the vectors of different genome DNAs. On
the other hand, Misrepair of DNA is a source of DNA mutations and the DNA
Misrepairs in germ cells may contribute to the diversity of genome DNAs in a
species. In conclusion, Misrepair mechanism is a mechanism essential for
individual adaptation, species adaptation, and species evolution.


  Helicases, involved in a number of cellular functions, are motors that
translocate along singlestranded nucleic acid and couple the motion to
unwinding double-strands of a duplex nucleic acid. The junction between double
and single strands creates a barrier to the movement of the helicase, which can
be manipulated in vitro by applying mechanical forces directly on the nucleic
acid strands. Single molecule experiments have demonstrated that the unwinding
velocities of some helicases increase dramatically with increase in the
external force, while others show little response. In contrast, the unwinding
processivity always increases when the force increases. The differing responses
of the unwinding velocity and processivity to force has lacked explanation. By
generalizing a previous model of processive unwinding by helicases, we provide
a unified framework for understanding the dependence of velocity and
processivity on force and the nucleic acid sequence. We predict that the
sensitivity of unwinding processivity to external force is a universal feature
that should be observed in all helicases. Our prediction is illustrated using
T7 and NS3 helicases as case studies. Interestingly, the increase in unwinding
processivity with force depends on whether the helicase forces base pair
opening by direct interaction or if such a disruption occurs spontaneously due
to thermal uctuations. Based on the theoretical results, we propose that
proteins like single-strand binding proteins associated with helicases in the
replisome, may have co-evolved with helicases to increase the unwinding
processivity even if the velocity remains unaffected.


  Several algorithms for RNA inverse folding have been used to design synthetic
riboswitches, ribozymes and thermoswitches, whose activity has been
experimentally validated. The RNAiFold software is unique among approaches for
inverse folding in that (exhaustive) constraint programming is used instead of
heuristic methods. For that reason, RNAiFold can generate all sequences that
fold into the target structure, or determine that there is no solution.
RNAiFold 2.0 is a complete overhaul of RNAiFold 1.0, rewritten from the now
defunct COMET language to C++. The new code properly extends the capabilities
of its predecessor by providing a user-friendly pipeline to design synthetic
constructs having the functionality of given Rfam families. In addition, the
new software supports amino acid constraints, even for proteins translated in
different reading frames from overlapping coding sequences; moreover, structure
compatibility/incompatibility constraints have been expanded. With these
features, RNAiFold 2.0 allows the user to design single RNA molecules as well
as hybridization complexes of two RNA molecules.
  The web server, source code and linux binaries are publicly accessible at
http://bioinformatics.bc.edu/clotelab/RNAiFold2.0


  The sample frequency spectrum (SFS) of DNA sequences from a collection of
individuals is a summary statistic which is commonly used for parametric
inference in population genetics. Despite the popularity of SFS-based inference
methods, currently little is known about the information-theoretic limit on the
estimation accuracy as a function of sample size. Here, we show that using the
SFS to estimate the size history of a population has a minimax error of at
least $O(1/\log s)$, where $s$ is the number of independent segregating sites
used in the analysis. This rate is exponentially worse than known convergence
rates for many classical estimation problems in statistics. Another surprising
aspect of our theoretical bound is that it does not depend on the dimension of
the SFS, which is related to the number of sampled individuals. This means
that, for a fixed number $s$ of segregating sites considered, using more
individuals does not help to reduce the minimax error bound. Our result
pertains to populations that have experienced a bottleneck, and we argue that
it can be expected to apply to many populations in nature.


  We develop a cross-platform open-source Java application (BACOM2) with
graphic user interface (GUI), and users also can use a XML file to set the
parameters of algorithm model, file paths and the dataset of paired samples.
BACOM2 implements the new entire pipeline of copy number change analysis for
heterogeneous cancer tissues, including extraction of raw copy number signals
from CEL files of paired samples, attenuation correction, identification of
balanced AB-genotype loci, copy number detection and segmentation, global
baseline calculation and absolute normalization, differentiation of deletion
types, estimation of the normal tissue fraction and correction of normal tissue
contamination. BACOM2 focuses on the common tools for data preparation and
absolute normalization for copy number analysis of heterogeneous cancer
tissues. The software provides an additional choice for scientists who require
a user-friendly, high-speed processing, cross-platform computing environment
for large copy number data analysis.


  A conjecture of Bandelt and Dress states that the maximum quartet distance
between any two phylogenetic trees on $n$ leaves is at most $(\frac 23
+o(1))\binom{n}{4}$. Using the machinery of flag algebras we improve the
currently known bounds regarding this conjecture, in particular we show that
the maximum is at most $(0.69 +o(1))\binom{n}{4}$. We also give further
evidence that the conjecture is true by proving that the maximum distance
between caterpillar trees is at most $(\frac 23 +o(1))\binom{n}{4}$.


  Gradient sensing requires at least two measurements at different points in
space. These measurements must then be communicated to a common location to be
compared, which is unavoidably noisy. While much is known about the limits of
measurement precision by cells, the limits placed by the communication are not
understood. Motivated by recent experiments, we derive the fundamental limits
to the precision of gradient sensing in a multicellular system, accounting for
communication and temporal integration. The gradient is estimated by comparing
a "local" and a "global" molecular reporter of the external concentration,
where the global reporter is exchanged between neighboring cells. Using the
fluctuation-dissipation framework, we find, in contrast to the case when
communication is ignored, that precision saturates with the number of cells
independently of the measurement time duration, since communication establishes
a maximum lengthscale over which sensory information can be reliably conveyed.
Surprisingly, we also find that precision is improved if the local reporter is
exchanged between cells as well, albeit more slowly than the global reporter.
The reason is that while exchange of the local reporter weakens the comparison,
it decreases the measurement noise. We term such a model "regional
excitation--global inhibition" (REGI). Our results demonstrate that fundamental
sensing limits are necessarily sharpened when the need to communicate
information is taken into account.


  To examine the co-evolution of quantitative traits in hosts and parasites, we
present and study a co-evolutionary model of a social parasite-host system that
incorporates (1) ecological dynamics that feed back into their co-evolutionary
outcomes; (2) variation in whether the parasite is obligate or facultative; and
(3) Holling Type II functional responses between host and parasite, which are
particularly suitable for social parasites that face time costs for host
location and its social manipulation. We perform local and global analyses for
the co-evolutionary model and the corresponding ecological model. In the
absence of evolution, the facultative parasite system can have one, two, or
three interior equilibria, while the obligate parasite system can have either
one or three interior equilibria. Multiple interior equilibria result in rich
dynamics with multiple attractors. The ecological system, in particular, can
exhibit bi-stability between the facultative-parasite-only equilibrium and the
interior coexistence equilibrium when it has two interior equilibria. Our
findings suggest that: (a) The host and parasite can select different
strategies that may result in local extinction of one species. These strategies
can have convergence stability (CS), but may not be evolutionary stable
strategies (ESS); (b) The host and its facultative (or obligate) parasite can
have ESS that drive the host (or the obligate parasite) extinct locally; (c)
Trait functions play an important role in the CS of both boundary and interior
equilibria, as well as their ESS; and (d) A small variance in the trait
difference that measures parasitism efficiency can destabilize the
co-evolutionary system, and generate evolutionary arms-race dynamics with
different host-parasite fluctuating patterns.


  One essential ingredient of evolutionary theory is the concept of fitness as
a measure for a species' success in its living conditions. Here, we quantify
the effect of environmental fluctuations onto fitness by analytical
calculations on a general evolutionary model and by studying corresponding
individual-based microscopic models. We demonstrate that not only larger growth
rates and viabilities, but also reduced sensitivity to environmental
variability substantially increases the fitness. Even for neutral evolution,
variability in the growth rates plays the crucial role of strongly reducing the
expected fixation times. Thereby, environmental fluctuations constitute a
mechanism to account for the effective population sizes inferred from genetic
data that often are much smaller than the census population size.


  Emergence is a phenomenon taken for granted in science but also still not
well understood. We have developed a model of artificial genetic evolution
intended to allow for emergence on genetic, population and social levels. We
present the details of the current state of our environment, agent, and
reproductive models. In developing our models we have relied on a principle of
using non-linear systems to model as many systems as possible including
mutation and recombination, gene-environment interaction, agent metabolism,
agent survival, resource gathering and sexual reproduction. In this paper we
review the genetic dynamics that have emerged in our system including
genotype-phenotype divergence, genetic drift, pseudogenes, and gene
duplication. We conclude that emergence-focused design in complex system
simulation is necessary to reproduce the multilevel emergence seen in the
natural world.


  In this article the notion of metabolic turnover is revisited in the light of
recent results of out-of-equilibrium thermodynamics. By means of Monte Carlo
methods we perform an exact uniform sampling of the steady state fluxes in a
genome scale metabolic network of E Coli from which we infer the metabolites
turnover times. However the latter are inferred from net fluxes, and we argue
that this approximation is not valid for enzymes working nearby thermodynamic
equilibrium. We recalculate turnover times from total fluxes by performing an
energy balance analysis of the network and recurring to the fluctuation
theorem. We find in many cases values one of order of magnitude lower, implying
a faster picture of intermediate metabolism.


  Clinical studies have indicated that in malignant surveillances fluctuations
in the population of certain effector T-cell repertoire become suppressed.
Motivated by such observations and in an attempt to quantify adaptive human
response to pathogens, we define an immune response function (IMRF) in terms of
mean square fluctuations of T-cell concentrations. We employ a recently
developed kinetic model of T-cell regulation that contains the essential
immunosuppressive effects of vitamin-D. We employ Gillespie algorithm to make
the first study of fluctuations along the stochastic trajectories. This
fluctuation-based IMRF can differentiate responses of different individuals
after pathogenic incursion both under healthy and disease conditions. We find
that relative fluctuations in T-cells (and hence IMRF) are different in
strongly regulated (malignant prone) and weakly regulated (autoimmune prone)
regions. The cross-over from one steady state (weakly regulated) to the other
(strongly regulated) is accompanied by a divergence-like growth in the
fluctuation of both the effector and regulatory T-cell concentration over a
wide range of pathogenic stimulation, displaying a dynamical phase transition
like behavior. The growth in fluctuation in this desired immune response regime
is found to arise from an intermittent fluctuation between regulatory and
effector T-cells that results in a bimodal distribution of population of each,
indicating bistability. The signature of intermittent behavior is further
confirmed by calculating the power spectrum of the corresponding fluctuation of
time correlation function. The calculated time correlation functions of
fluctuations show that the slow fluctuation causes the bistabilty in healthy
state. Thus, in diseases diagnosis process, such steady state response
parameters can provide immense information which might become helpful to define
an immune status.


  In this paper, we study the distribution and behaviour of internal equilibria
in a $d$-player $n$-strategy random evolutionary game where the game payoff
matrix is generated from normal distributions. The study of this paper reveals
and exploits interesting connections between evolutionary game theory and
random polynomial theory. The main novelties of the paper are some qualitative
and quantitative results on the expected density, $f_{n,d}$, and the expected
number, $E(n,d)$, of (stable) internal equilibria. Firstly, we show that in
multi-player two-strategy games, they behave asymptotically as $\sqrt{d-1}$ as
$d$ is sufficiently large. Secondly, we prove that they are monotone functions
of $d$. We also make a conjecture for games with more than two strategies.
Thirdly, we provide numerical simulations for our analytical results and to
support the conjecture. As consequences of our analysis, some qualitative and
quantitative results on the distribution of zeros of a random Bernstein
polynomial are also obtained.


  Food webs represent the set of consumer-resource interactions among a set of
species that co-occur in a habitat, but most food web studies have omitted
parasites and their interactions. Recent studies have provided conflicting
evidence on whether including parasites changes food web structure, with some
suggesting that parasitic interactions are structurally distinct from those
among free-living species while others claim the opposite. Here, we describe a
principled method for understanding food web structure that combines an
efficient optimization algorithm from statistical physics called parallel
tempering with a probabilistic generalization of the empirically well-supported
food web niche model. This generative model approach allows us to rigorously
estimate the degree to which interactions that involve parasites are
statistically distinguishable from interactions among free-living species,
whether parasite niches behave similarly to free-living niches, and the degree
to which existing hypotheses about food web structure are naturally recovered.
We apply this method to the well-studied Flensburg Fjord food web and show that
while predation on parasites, concomitant predation of parasites, and parasitic
intraguild trophic interactions are largely indistinguishable from free-living
predation interactions, parasite-host interactions are different. These results
provide a powerful new tool for evaluating the impact of classes of species and
interactions on food web structure to shed new light on the roles of parasites
in food webs


  Michaelis-Menten equation is a basic equation of enzyme kinetics and gives an
acceptable approximation of real chemical reaction processes. Analyzing the
derivation of this equation yields the fact that its good performance of
approximating real reaction processes is due to Michaelis-Menten curve (15).
This curve is derived from Quasi-Steady-State Assumption(QSSA), which has been
proved always true and called Quasi-Steady-State Law by Banghe Li et al [19].
  Here, we found a quartic equation A(S,E)=0 (22), which gives more accurate
approximation of the reaction process in two aspects: during the quasi-steady
state of a reaction, Michaelis-Menten curve approximates the reaction well,
while our quartic equation $A(S,E)=0$ gives better approximation; near the end
of the reaction, our equation approaches the end of the reaction with a tangent
line same to that of the reaction, while Michaelis-Menten curve does not. In
addition, our quartic equation A(S,E)=0 differs to Michaelis-Menten curve less
than the order of $1/S^3$ as S approaches $+\infty$.
  By considering the above merits of A(S,E)=0, we suggest it as a replacement
of Michaelis-Menten curve. Intuitively, this new equation is more complex and
harder to understand. But, just because its complexity, it provides more
information about the rate constants than Michaelis-Menten curve does.
  Finally, we get a better replacement of the Michaelis-Menten equation by
combing A(S,E)=0 and the equation $dP/dt=k_2C(t)$.


  The increasing observation of mutual exclusivity correlations among cancer
gene mutations is a key component for identifying driver events or pathways in
cancer genome analysis. Here we report a rigorous statistical method to compute
an exact $p$-value for the beyond-pairwise mutual exclusivity or co-occurrence
relationships among cancer gene mutations by enumerating a null distribution of
overlapping mutations across more than two genes. The validity and the
advantage of our method is explicitly demonstrated in both cancer gene
mutations and simulation data through the comparison to the permutation test.


  Accurate subtyping or classification of breast cancer is important for
ensuring proper treatment of patients and also for understanding the molecular
mechanisms driving this disease. While there have been several gene signatures
proposed in the literature to classify breast tumours, these signatures show
very low overlaps, different classification performance, and not much relevance
to the underlying biology of these tumours. Here we evaluate DNA-damage
response (DDR) and cell cycle pathways, which are critical pathways implicated
in a considerable proportion of breast tumours, for their usefulness and
ability in breast tumour subtyping. We think that subtyping breast tumours
based on these two pathways could lead to vital insights into molecular
mechanisms driving these tumours. Here, we performed a systematic evaluation of
DDR and cell-cycle pathways for subtyping of breast tumours into the five known
intrinsic subtypes. Homologous Recombination (HR) pathway showed the best
performance in subtyping breast tumours, indicating that HR genes are strongly
involved in all breast tumours. Comparisons of pathway based signatures and two
standard gene signatures supported the use of known pathways for breast tumour
subtyping. Further, the evaluation of these standard gene signatures showed
that breast tumour subtyping, prognosis and survival estimation are all closely
related. Finally, we constructed an all-inclusive super-signature by combining
(union of) all genes and performing a stringent feature selection, and found it
to be reasonably accurate and robust in classification as well as prognostic
value. Adopting DDR and cell cycle pathways for breast tumour subtyping
achieved robust and accurate breast tumour subtyping, and constructing a
super-signature which contains feature selected mix of genes from these
molecular pathways as well as clinical aspects is valuable in clinical
practice.


  Background Next Generation Sequencing (NGS) has dramatically enhanced our
ability to sequence genomes, but not to assemble them. In practice, many
published genome sequences remain in the state of a large set of contigs. Each
contig describes the sequence found along some path of the assembly graph,
however, the set of contigs does not record all the sequence information
contained in that graph. Although many subsequent analyses can be performed
with the set of contigs, one may ask whether mapping reads on the contigs is as
informative as mapping them on the paths of the assembly graph. Currently, one
lacks practical tools to perform mapping on such graphs. Results Here, we
propose a formal definition of mapping on a de Bruijn graph, analyse the
problem complexity which turns out to be NP-complete, and provide a practical
solution.We propose a pipeline called GGMAP (Greedy Graph MAPping). Its novelty
is a procedure to map reads on branching paths of the graph, for which we
designed a heuristic algorithm called BGREAT (de Bruijn Graph REAd mapping
Tool). For the sake of efficiency, BGREAT rewrites a read sequence as a
succession of unitigs sequences. GGMAP can map millions of reads per CPU hour
on a de Bruijn graph built from a large set of human genomic reads.
Surprisingly, results show that up to 22% more reads can be mapped on the graph
but not on the contig set. Conclusions Although mapping reads on a de Bruijn
graph is complex task, our proposal offers a practical solution combining
efficiency with an improved mapping capacity compared to assembly-based mapping
even for complex eukaryotic data. Availability: github.com/Malfoy/BGREAT
Keywords: Read mapping; De bruijn graphs; NGS; NP-completeness