  Motivation: LTR retrotransposons are mobile elements that are able, like
retroviruses, to copy and move inside eukaryotic genomes. In the present work,
we propose a branching model for studying the propagation of LTR
retrotransposons in these genomes. This model allows to take into account both
positions and degradations of LTR retrotransposons copies. In our model, the
duplication rate is also allowed to vary with the degradation level.
  Results: Various functions have been implemented in order to simulate their
spread and visualization tools are proposed. Based on these simulation tools,
we show that an accurate estimation of the parameters of this propagation model
can be performed. We applied this method to the study of the spread of the
transposable elements ROO, GYPSY, and DM412 on a chromosome of
\textit{Drosophila melanogaster}.
  Availability: Our proposal has been implemented using Python software. Source
code is freely available on the web at
https://github.com/SergeMOULIN/retrotransposons-spread.


  Exploiting pathogen genomes to reconstruct transmission represents a powerful
tool in the fight against infectious disease. However, their interpretation
rests on a number of simplifying assumptions that regularly ignore important
complexities of real data, in particular within-host evolution and non-sampled
patients.
  Here we propose a new approach to transmission inference called SCOTTI
(Structured COalescent Transmission Tree Inference). This method is based on a
statistical framework that models each host as a distinct population, and
transmissions between hosts as migration events. Our computationally efficient
implementation of this model enables the inference of host-to-host transmission
while accommodating within-host evolution and non-sampled hosts. SCOTTI is
distributed as an open source package for the phylogenetic software BEAST2.
  We show that SCOTTI can generally infer transmission events even in the
presence of considerable within-host variation, can account for the uncertainty
associated with the possible presence of non-sampled hosts, and can efficiently
use data from multiple samples of the same host, although there is some
reduction in accuracy when samples are collected very close to the infection
time.
  We illustrate the features of our approach by investigating transmission from
genetic and epidemiological data in a Foot and Mouth Disease Virus (FMDV)
veterinary outbreak in England and a Klebsiella pneumoniae outbreak in a Nepali
neonatal unit. Transmission histories inferred with SCOTTI will be important in
devising effective measures to prevent and halt transmission.


  Systematic identification of protein function is a key problem in current
biology. Most traditional methods fail to identify functionally equivalent
proteins if they lack similar sequences, structural data or extensive manual
annotations. In this thesis, I focused on feature engineering and machine
learning methods for identifying diverse classes of proteins that share
functional relatedness but little sequence or structural similarity, notably,
Neuropeptide Precursors (NPPs).
  I aim to identify functional protein classes solely using unannotated protein
primary sequences from any organism. This thesis focuses on feature
representations of whole protein sequences, sequence derived engineered
features, their extraction, frameworks for their usage by machine learning (ML)
models, and the application of ML models to biological tasks, focusing on high
level protein functions. I implemented the ideas of feature engineering to
develop a platform (called NeuroPID) that extracts meaningful features for
classification of overlooked NPPs. The platform allows mass discovery of new
NPs and NPPs. It was expanded as a webserver.
  I expanded our approach towards other challenging protein classes. This is
implemented as a novel bioinformatics toolkit called ProFET (Protein Feature
Engineering Toolkit). ProFET extracts hundreds of biophysical and sequence
derived attributes, allowing the application of machine learning methods to
proteins. ProFET was applied on many protein benchmark datasets with state of
the art performance. The success of ProFET applies to a wide range of
high-level functions such as metagenomic analysis, subcellular localization,
structure and unique functional properties (e.g. thermophiles, nucleic acid
binding).
  These methods and frameworks represent a valuable resource for using ML and
data science methods on proteins.


  Empirical studies show that preference for prevention versus treatment
remains a subject of debate. We build a paradigm model combining a utility game
for the individual-level dilemma of prevention versus treatment, and a
compartmental model for the epidemic dynamic. We assume that individuals arrive
to maximize the utility of voluntary prevention, as the epidemic reaches an
endemic level alleviated by prevention and treatment. We thus obtain an
expression for the asymptotic prevention coverage. Notably, we obtain that, if
the relative cost of prevention versus treatment is sufficiently low, epidemics
may be averted through the use of prevention alone.


  We define data analyses to monitor a change in R, the average number of
secondary cases caused by a typical infected individual. The input dataset
consists of incident cases partitioned into outbreaks, each initiated from a
single index case. We split of the input dataset into two successive subsets,
to evaluate two successive R values, according to the Bayesian paradigm. We
used the Bayes factor between the model with two different R values and that
with a single R value to justify that the change in R is statistically
significant. We validated our approach using simulated data, generated using
known R. In particular, we found that claiming two distinct R values may depend
significantly on the number of outbreaks. We then reanalyzed data previously
studied by Jansen et al. [Jansen et al. Science 301 (5634), 804], concerning
the effective reproduction number for measles in the UK, during 1995-2002. Our
analyses showed that the 1995-2002 dataset should be divided into two separate
subsets for the periods 1995-1998 and 1999-2002. In contrast, Jansen et al.
take this splitting point as input of their analysis. Our estimated effective
reproduction numbers R are in good agreement with those found by Jansen et al.
In conclusion, our methodology for detecting temporal changes in R using
outbreak-size data worked satisfactorily with both simulated and real-world
data. The methodology may be used for updating R in real time, as surveillance
outbreak data become available.


  We derive from first principles a three-dimensional theory of self-propelled
particle swarming in a viscous fluid environment. Our model predicts emergent
collective behavior that depends critically on fluid opacity, mechanism of
self-propulsion, and type of particle-particle interaction. In "clear fluids"
swimmers have full knowledge of their surroundings and can adjust their
velocities with respect to the lab frame, while in "opaque fluids," they
control their velocities only in relation to the local fluid flow. We also show
that "social" interactions that affect only a particle's propensity to swim
towards or away from neighbors induces a flow field that is qualitatively
different from the long-ranged flow fields generated by direct "physical"
interactions. The latter can be short-ranged but lead to much longer-ranged
fluid-mediated hydrodynamic forces, effectively amplifying the range over which
particles interact. These different fluid flows conspire to profoundly affect
swarm morphology, kinetically stabilizing or destabilizing swarm configurations
that would arise in the absence of fluid. Depending upon the overall
interaction potential, the mechanism of swimming (e.g., pushers or pullers),
and the degree of fluid opaqueness, we discover a number of new collective
three-dimensional patterns including flocks with prolate or oblate shapes,
recirculating peloton-like structures, and jet-like fluid flows that entrain
particles mediating their escape from the center of mill-like structures. Our
results reveal how the interplay among general physical elements influence
fluid-mediated interactions and the self-organization, mobility, and stability
of new three-dimensional swarms and suggest how they might be used to
kinetically control their collective behavior.


  We propose a series of methods and models in order to explore the Global
Burden of Disease Study and the provided healthy life expectancy HALE estimates
from the World Health Organization WHO based on the mortality mx of a
population provided in a classical life table and a mortality diagram. Our
estimates are compared with the HALE estimates for the World territories and
the WHO regions along with providing comparative results with to findings of
Chang, Molla, Truman et al. (2015) on the Differences in healthy life
expectancy for the US population by sex, race or ethnicity and geographic
region in 2008 and from Yong and Saito (2009) regarding Trends in healthy life
expectancy in Japan. From the mortality point of view we have developed a
simple model for the estimation of a characteristic parameter b related to the
healthy life years lost to disability and providing full application details
along with characteristic parameter selection and stability of the
coefficients. We also provide a direct estimation method of the parameter b
from the life tables. We straighten the importance of our methodology by
proposing and applying estimates of the parameter b by using the Gompertz and
the Weibull models. From the Health State point of view we summarize the main
points of the first exit time theory to life table data and present the basic
models starting from the first related model published by Janssen and Skiadas
(1995). Even more we develop the simpler 2-parameter health state model and an
extension of a model expressing the infant mortality to a 4-parameter model
which is the simpler model providing very good fitting on the logarithm of the
force of mortality. More important is the use of the Health State Function and
the relative impact on mortality to find an estimate for the healthy life years
lost to disability.


  The fidelity of DNA replication by DNA polymerase (DNAP) has long been an
important issue in biology. While numerous experiments have revealed details of
the molecular structure and working mechanism of DNAP which consists of both a
polymerase site and an exonuclease (proofreading) site, there were quite few
theoretical studies on the fidelity issue. The first model which explicitly
considered both sites was proposed in 1970s' and the basic idea was widely
accepted by later models. However, all these models did not systematically and
rigorously investigate the dominant factor on DNAP fidelity, i.e, the
higher-order terminal effects through which the polymerization pathway and the
proofreading pathway coordinate to achieve high fidelity. In this paper, we
propose a new and comprehensive kinetic model of DNAP based on some recent
experimental observations, which includes previous models as special cases. We
present a rigorous and unified treatment of the corresponding steady-state
kinetic equations of any-order terminal effects, and derive analytical
expressions for fidelity in terms of kinetic parameters under bio-relevant
conditions. These expressions offer new insights on how the the higher-order
terminal effects contribute substantially to the fidelity in an order-by-order
way, and also show that the polymerization-and-proofreading mechanism is
dominated only by very few key parameters. We then apply these results to
calculate the fidelity of some real DNAPs, which are in good agreements with
previous intuitive estimates given by experimentalists.


  Chemical reaction networks can be automatically generated from graph grammar
descriptions, where rewrite rules model reaction patterns. Because a molecule
graph is connected and reactions in general involve multiple molecules, the
rewriting must be performed on multisets of graphs. We present a general
software package for this type of graph rewriting system, which can be used for
modelling chemical systems. The package contains a C++ library with algorithms
for working with transformation rules in the Double Pushout formalism, e.g.,
composition of rules and a domain specific language for programming graph
language generation. A Python interface makes these features easily accessible.
The package also has extensive procedures for automatically visualising not
only graphs and rewrite rules, but also Double Pushout diagrams and graph
languages in form of directed hypergraphs. The software is available as an open
source package, and interactive examples can be found on the accompanying
webpage.


  Genes and proteins regulate cellular functions through complex circuits of
biochemical reactions. Fluctuations in the components of these regulatory
networks result in noise that invariably corrupts the signal, possibly
compromising function. Here, we create a practical formalism based on ideas
introduced by Wiener and Kolmogorov (WK) for filtering noise in engineered
communications systems to quantitatively assess the extent to which noise can
be controlled in biological processes involving negative feedback. Application
of the theory, which reproduces the previously proven scaling of the lower
bound for noise suppression in terms of the number of signaling events, shows
that a tetracycline repressor-based negative-regulatory gene circuit behaves as
a WK filter. For the class of Hill-like nonlinear regulatory functions, this
type of filter provides the optimal reduction in noise. Our theoretical
approach can be readily combined with experimental measurements of response
functions in a wide variety of genetic circuits, to elucidate the general
principles by which biological networks minimize noise.


  We present a sequential Monte Carlo algorithm for Markov chain trajectories
with proposals constructed in reverse time, which is advantageous when paths
are conditioned to end in a rare set. The reverse time proposal distribution is
constructed by approximating the ratio of Green's functions in Nagasawa's
formula. Conditioning arguments can be used to interpret these ratios as
low-dimensional conditional sampling distributions of some coordinates of the
process given the others. Hence the difficulty in designing SMC proposals in
high dimension is greatly reduced. We illustrate our method on estimating an
overflow probability in a queueing model, the probability that a diffusion
follows a narrowing corridor, and the initial location of an infection in an
epidemic model on a network.


  Numerical plant models can predict the outcome of plant traits modifications
resulting from genetic variations, on plant performance, by simulating
physiological processes and their interaction with the environment.
Optimization methods complement those models to design ideotypes, i.e. ideal
values of a set of plant traits resulting in optimal adaptation for given
combinations of environment and management, mainly through the maximization of
a performance criteria (e.g. yield, light interception). As use of simulation
models gains momentum in plant breeding, numerical experiments must be
carefully engineered to provide accurate and attainable results, rooting them
in biological reality. Here, we propose a multi-objective optimization
formulation that includes a metric of performance, returned by the numerical
model, and a metric of feasibility, accounting for correlations between traits
based on field observations. We applied this approach to two contrasting
models: a process-based crop model of sunflower and a functional-structural
plant model of apple trees. In both cases, the method successfully
characterized key plant traits and identified a continuum of optimal solutions,
ranging from the most feasible to the most efficient. The present study thus
provides successful proof of concept for this enhanced modeling approach, which
identified paths for desirable trait modification, including direction and
intensity.


  A major Ebola outbreak occurs in West Africa since March 2014, being the
deadliest epidemic in history. As an infectious disease epidemiology, Ebola is
the most lethal and is moving faster than in previous outbreaks. On 8 August
2014, the World Health Organization (WHO) declared the outbreak a public health
emergency of international concern. Last update on 7 July 2015 by WHO reports
27 609 cases of Ebola with a total of 11 261 deaths. In this work, we present a
mathematical description of the spread of Ebola virus based on the SEIR
(Susceptible-Exposed-Infective-Recovered) model and optimal strategies for
Ebola control. In order to control the propagation of the virus and to predict
the impact of vaccine programmes, we investigate several strategies of optimal
control of the spread of Ebola: control infection by vaccination of
susceptible; minimize exposed and infected; reduce Ebola infection by
vaccination and education.


  We explore a definition of complexity based on logic functions, which are
widely used as compact descriptions of rules in diverse fields of contemporary
science. Detailed numerical analysis shows that (i) logic complexity is
effective in discriminating between classes of functions commonly employed in
modelling contexts; (ii) it extends the notion of canalisation, used in the
study of genetic regulation, to a more general and detailed measure; (iii) it
is tightly linked to the resilience of a function's output to noise affecting
its inputs. We demonstrate its utility by measuring it in empirical data on
gene regulation, digital circuitry, and propositional calculus. Logic
complexity is exceptionally low in these systems. The asymmetry between "on"
and "off" states in the data correlates with the complexity in a non-null way;
a model of random Boolean networks clarifies this trend and indicates a common
hierarchical architecture in the three systems.


  Living cells display a remarkable capacity to compartmentalize their
functional biochemistry. A particularly fascinating example is the cell
nucleus. Exchange of macromolecules between the nucleus and the surrounding
cytoplasm does not involve traversing a lipid bilayer membrane. Instead, large
protein channels known as nuclear pores cross the nuclear envelope and regulate
the passage of other proteins and RNA molecules. Beyond simply gating
diffusion, the system of nuclear pores and associated transport receptors is
able to generate substantial concentration gradients, at the energetic expense
of guanosine triphosphate (GTP) hydrolysis. In contrast to conventional
approaches to demixing such as reverse osmosis or dialysis, the biological
system operates continuously, without application of cyclic changes in pressure
or solvent exchange. Abstracting the biological paradigm, we examine this
transport system as a thermodynamic machine of solution demixing. Building on
the construct of free energy transduction and biochemical kinetics, we find
conditions for stable operation and optimization of the concentration gradients
as a function of dissipation in the form of entropy production.


  Over the last several decades it has been increasingly recognized that
stochastic processes play a central role in transcription. Though many
stochastic effects have been explained, the source of transcriptional bursting
(one of the most well-known sources of stochasticity) has continued to evade
understanding. Recent results have pointed to mechanical feedback as the source
of transcriptional bursting but a reconciliation of this perspective with
preexisting views of transcriptional regulation is lacking. In this letter we
present a simple phenomenological model which is able to incorporate the
traditional view of gene expression within a framework with mechanical limits
to transcription. Our model explains the emergence of universal properties of
gene expression, wherein the lower limit of intrinsic noise necessarily rises
with mean expression level.


  We consider the Combinatorial RNA Design problem, a minimal instance of RNA
design where one must produce an RNA sequence that adopts a given secondary
structure as its minimal free-energy structure. We consider two free-energy
models where the contributions of base pairs are additive and independent: the
purely combinatorial Watson-Crick model, which only allows equally-contributing
A -- U and C -- G base pairs, and the real-valued Nussinov-Jacobson model,
which associates arbitrary energies to A -- U, C -- G and G -- U base pairs. We
first provide a complete characterization of designable structures using
restricted alphabets and, in the four-letter alphabet, provide a complete
characterization for designable structures without unpaired bases. When
unpaired bases are allowed, we characterize extensive classes of
(non-)designable structures, and prove the closure of the set of designable
structures under the stutter operation. Membership of a given structure to any
of the classes can be tested in $\Theta$(n) time, including the generation of a
solution sequence for positive instances. Finally, we consider a
structure-approximating relaxation of the design, and provide a $\Theta$(n)
algorithm which, given a structure S that avoids two trivially non-designable
motifs, transforms S into a designable structure constructively by adding at
most one base-pair to each of its stems.


  Using graphical methods based on a `lookdown' and pruned version of the {\em
ancestral selection graph}, we obtain a representation of the type distribution
of the ancestor in a two-type Wright-Fisher population with mutation and
selection, conditional on the overall type frequency in the old population.
This extends results from Lenz, Kluth, Baake, and Wakolbinger (Theor. Pop.
Biol., 103 (2015), 27-37) to the case of heavy-tailed offspring, directed by a
reproduction measure $\Lambda$. The representation is in terms of the
equilibrium tail probabilities of the line-counting process $L$ of the graph.
We identify a strong pathwise Siegmund dual of $L$, and characterise the
equilibrium tail probabilities of $L$ in terms of hitting probabilities of the
dual process.


  The dynamics of an association of interactive biological species is studied
theoretically. We explore a mean field approximation to describe the temporal
evolution of an ecological system with the basic prey-predator interspecies
relation, as well as an approximation to introduce time correlations in the
dynamics. We start by discussing the solution of the Volterra-Lotka model in a
mean field approximation based in an analogy with the Weiss solution to the
Ising model for ferromagnetic materials. In order to explore the effects of
long-range time correlations, we describe the time evolution of the system
within a kind of Bragg-Williams approximation. This approach allows us to
evaluate a characteristic life-time of the ecosystem. This quantity could be
very useful to discuss the time evolution of the system under a wide diversity
of environmental conditions of the ecosystem which is not usually considered.
We discuss the general trends of the temporal evolution of the association with
some data from real ecosystems.


  Given a random RNA secondary structure, $S$, we study RNA sequences having
fixed ratios of nuclotides that are compatible with $S$. We perform this
analysis for RNA secondary structures subject to various base pairing rules and
minimum arc- and stack-length restrictions. Our main result reads as follows:
in the simplex of the nucleotide ratios there exists a convex region in which,
in the limit of long sequences, a random structure a.a.s.~has compatible
sequence with these ratios and outside of which a.a.s.~a random structure has
no such compatible sequence. We localize this region for RNA secondary
structures subject to various base pairing rules and minimum arc- and
stack-length restrictions. In particular, for {\bf GC}-sequences having a ratio
of {\bf G} nucleotides smaller than $1/3$, a random RNA secondary structure
without any minimum arc- and stack-length restrictions has a.a.s.~no such
compatible sequence. For sequences having a ratio of {\bf G} nucleotides larger
than $1/3$, a random RNA secondary structure has a.a.s. such compatible
sequences. We discuss our results in the context of various families of RNA
structures.


  The goal of these lectures is to review some mathematical aspects of random
tree models used in evolutionary biology to model gene trees or species trees.
We start with stochastic models of tree shapes (finite trees without edge
lengths), culminating in the $\beta$-family of Aldous' branching models. We
next introduce real trees (trees as metric spaces) and show how to study them
through their contour, provided they are properly measured and ordered. We then
focus on the reduced tree, or coalescent tree, which is the tree spanned by
individuals/species alive at the same fixed time. We show how reduced trees,
like any compact ultrametric space, can be represented in a simple way via the
so-called comb metric. Beautiful examples of random combs include the Kingman
coalescent and coalescent point processes. We end up displaying some recent
biological applications of coalescent point processes to the inference of
species diversification, to conservation biology and to epidemiology.


  Tick-Borne diseases can be transmitted via non-systemic (NS) transmission.
This occurs when tick gets the infection by co-feeding with infected ticks on
the same host resulting in a direct pathogen transmission between the vectors,
without infecting the host. This transmission is peculiar, as it does not
require any systemic infection of the host. The NS transmission is the main
efficient transmission for the persistence of the Tick-Borne Encephalitis virus
in nature. By describing the heterogeneous ticks aggregation on hosts through a
\hyphenation{dynamical} bipartite graphs representation, we are able to
mathematically define the NS transmission and to depict the epidemiological
conditions for the pathogen persistence. Despite the fact that the underlying
network is largely fragmented, analytical and computational results show that
the larger is the variability of the aggregation, and the easier is for the
pathogen to persist in the population.


  Food webs -- networks of predators and prey -- have long been known to
exhibit "intervality": species can generally be ordered along a single axis in
such a way that the prey of any given predator tend to lie on unbroken compact
intervals. Although the meaning of this axis -- identified with a "niche"
dimension -- has remained a mystery, it is assumed to lie at the basis of the
highly non-trivial structure of food webs. With this in mind, most trophic
network modelling has for decades been based on assigning species a niche value
by hand. However, we argue here that intervality should not be considered the
cause but rather a consequence of food-web structure. First, analysing a set of
$46$ empirical food webs, we find that they also exhibit {\it predator}
intervality: the predators of any given species are as likely to be contiguous
as the prey are, but in a different ordering. Furthermore, this property is not
exclusive of trophic networks: several networks of genes, neurons, metabolites,
cellular machines, airports, and words are found to be approximately as
interval as food webs. We go on to show that a simple model of food-web
assembly which does not make use of a niche axis can nevertheless generate
significant intervality. Therefore, the niche dimension (in the sense used for
food-web modelling) could in fact be the consequence of other, more fundamental
structural traits, such as trophic coherence. We conclude that a new approach
to food-web modelling is required for a deeper understanding of ecosystem
assembly, structure and function, and propose that certain topological features
thought to be specific of food webs are in fact common to many complex
networks.


  This paper is not (or at least not only) about human infant mortality. In
line with reliability theory, "infant" will refer here to the time interval
following birth during which the mortality (or failure) rate decreases. This
definition provides a systems science perspective in which birth constitutes a
sudden transition which falls within the field of application of the "Transient
Shock" (TS) conjecture put forward in Richmond et al. (2016c). This conjecture
provides predictions about the timing and shape of the death rate peak. (i) It
says that there will be a death rate spike whenever external conditions change
abruptly and drastically. (ii) It predicts that after a steep rising there will
be a much longer hyperbolic relaxation process. These predictions can be tested
by considering living organisms for which birth is a multi-step process. Thus,
for fish there are three states: egg, yolk-sac phase, young adult. The TS
conjecture predicts a mortality spike at the end of the yolk-sac phase, and
this timing is indeed confirmed by observation. Secondly, the hyperbolic nature
of the relaxation process can be tested using high accuracy Swiss statistics
which give postnatal death rates from one hour after birth up to the age of 10
years. It turns out that since the 19th century despite a great overall
reduction in infant mortality, the shape of the age-specific death rate has
remained basically unchanged. This hyperbolic pattern is not specific to
humans. It can also be found in small primates as recorded in the archives of
zoological gardens. Our ultimate objective is to set up a chain of cases which
starts from simple systems and then moves up step by step to more complex
organisms. The cases discussed here can be seen as initial landmarks.


  Poisson-Boltzmann (PB) model is one of the most popular implicit solvent
models in biophysical modeling and computation. The ability of providing
accurate and reliable PB estimation of electrostatic solvation free energy,
$\Delta G_{\text{el}}$, and binding free energy, $\Delta\Delta G_{\text{el}}$,
is of tremendous significance to computational biophysics and biochemistry.
Recently, it has been warned in the literature (Journal of Chemical Theory and
Computation 2013, 9, 3677-3685) that the widely used grid spacing of $0.5$ \AA
$ $ produces unacceptable errors in $\Delta\Delta G_{\text{el}}$ estimation
with the solvent exclude surface (SES). In this work, we investigate the grid
dependence of our PB solver (MIBPB) with SESs for estimating both electrostatic
solvation free energies and electrostatic binding free energies. It is found
that the relative absolute error of $\Delta G_{\text{el}}$ obtained at the grid
spacing of $1.0$ \AA $ $ compared to $\Delta G_{\text{el}}$ at $0.2$ \AA $ $
averaged over 153 molecules is less than 0.2\%. Our results indicate that the
use of grid spacing $0.6$ \AA $ $ ensures accuracy and reliability in
$\Delta\Delta G_{\text{el}}$ calculation. In fact, the grid spacing of $1.1$
\AA $ $ appears to deliver adequate accuracy for high throughput screening.


  In many species, sucrose en route out of the leaf migrates from
photosynthetically active mesophyll cells into the phloem down its
concentration gradient via plasmodesmata, i.e., symplastically. In some of
these plants the process is entirely passive, but in others phloem sucrose is
actively converted into larger sugars, raffinose and stachyose, and segregated
(trapped), thus raising total phloem sugar concentration to a level higher than
in the mesophyll. Questions remain regarding the mechanisms and selective
advantages conferred by both of these symplastic loading processes. Here we
present an integrated model - including local and global transport and the
kinetics of oligomerization - for passive and active symplastic loading. We
also propose a physical model of transport through the plasmodesmata. With
these models, we predict that: 1) relative to passive loading, oligomerization
of sucrose in the phloem, even in the absence of segregation, lowers the sugar
content in the leaf required to achieve a given export rate and accelerates
export for a given concentration of sucrose in the mesophyll; and 2)
segregation of oligomers and the inverted gradient of total sugar content can
be achieved for physiologically reasonable parameter values, but even higher
export rates can be accessed in scenarios in which polymers are allowed to
diffuse back into the mesophyll. We discuss these predictions in relation to
further studies aimed at the clarification of loading mechanisms, fitness of
active and passive symplastic loading, and potential targets for engineering
improved rates of export.


  The precision of concentration sensing is improved when cells communicate.
Here we derive the physical limits to concentration sensing for cells that
communicate over short distances by directly exchanging small molecules
(juxtacrine signaling), or over longer distances by secreting and sensing a
diffusive messenger molecule (autocrine signaling). In the latter case, we find
that the optimal cell spacing can be large, due to a tradeoff between
maintaining communication strength and reducing signal cross-correlations. This
leads to the surprising result that autocrine signaling allows more precise
sensing than juxtacrine signaling for sufficiently large populations. We
compare our results to data from a wide variety of communicating cell types.


  As an example of topic where biology and physics meet, we present the issue
of protein folding and stability, and the development of thermodynamics-based
bioinformatics tools that predict the stability and thermal resistance of
proteins and the change of these quantities upon amino acid substitutions.
These methods are based on knowledge-driven statistical potentials, derived
from experimental protein structures using the inverse Boltzmann law. We also
describe an application of these predictors, which contributed to the
understanding of the mechanisms of aggregation of a particular protein known to
cause a neuronal disease.


  The songs and calls of many bird species, like human speech, form distinct
regional dialects. We suggest that the process of dialect formation is
analogous to the physical process of magnetic domain formation. We take the
coastal breeding grounds of the Puget Sound white crowned sparrow as an
example. Previous field studies suggest that birds of this species learn
multiple songs early in life, and when establishing a territory for the first
time, retain one of these dialects in order to match the majority of their
neighbours. We introduce a simple lattice model of the process, showing that
this matching behaviour can produce single dialect domains provided the death
rate of adult birds is sufficiently low. We relate death rate to thermodynamic
temperature in magnetic materials, and calculate the critical death rate by
analogy with the Ising model. Using parameters consistent with the known
behavior of these birds we show that coastal dialect domain shapes may be
explained by viewing them as low temperature "stripe states".


  We discuss the transition paths in a coupled bistable system consisting of
interacting multiple identical bistable motifs. We propose a simple model of
coupled bistable gene circuits as an example, and show that its transition
paths are bifurcating. We then derive a criterion to predict the bifurcation of
transition paths in a generalized coupled bistable system. We confirm the
validity of the theory for the example system by numerical simulation. We also
demonstrate in the example system that, if the steady states of individual gene
circuits are not changed by the coupling, the bifurcation pattern is not
dependent on the number of gene circuits. We further show that the transition
rate exponentially decreases with the number of gene circuits when the
transition path does not bifurcate, while a bifurcation facilitates the
transition by lowering the quasi-potential energy barrier.


  This paper introduces a theoretical framework for the analysis and control of
the stochastic susceptible-infected-removed (SIR) spreading process over a
network of heterogeneous agents. In our analysis, we analyze the exact
networked Markov process describing the SIR model, without resorting to
mean-field approximations, and introduce a convex optimization framework to
find an efficient allocation of resources to contain the expected number of
accumulated infections over time. Numerical simulations are presented to
illustrate the effectiveness of the obtained results.


  Umbrella sampling efficiently yields equilibrium averages that depend on
exploring rare states of a model by biasing simulations to windows of
coordinate values and then combining the resulting data with physical
weighting. Here, we introduce a mathematical framework that casts the step of
combining the data as an eigenproblem. The advantage to this approach is that
it facilitates error analysis. We discuss how the error scales with the number
of windows. Then, we derive a central limit theorem for averages that are
obtained from umbrella sampling. The central limit theorem suggests an
estimator of the error contributions from individual windows, and we develop a
simple and computationally inexpensive procedure for implementing it. We
demonstrate this estimator for simulations of the alanine dipeptide and show
that it emphasizes low free energy pathways between stable states in comparison
to existing approaches for assessing error contributions. We discuss the
possibility of using the estimator and, more generally, the eigenvector method
for umbrella sampling to guide adaptation of the simulation parameters to
accelerate convergence.


  Despite major advances in quantitative approaches to natural resource
management, there has been resistance to using these tools in the actual
practice of managing ecological populations. Given a managed system and a set
of assumptions, translated into a model, optimization methods can be used to
solve for the most cost effective management actions. However, when the
underlying assumptions are not met, such methods can potentially lead to poor
decisions. Managers who develop decisions based on past experience and
judgment, without the aid of mathematical models, can potentially learn about
the system and develop flexible management strategies. However, these
strategies are often based on subjective criteria and equally invalid and often
unstated assumptions. Given the drawbacks of both methods, it is unclear
whether simple quantitative models improve environmental decision making over
expert opinion. In this paper, we explore how well students, using their
experience and judgment, manage simulated fishery populations in an online
computer game and compare their management outcomes to the performance of
model-based decisions. We consider harvest decisions generated using four
different quantitative models: 1. the model used to produce the simulated
population dynamics observed in the game, with all underlying parameter values
known [a control], 2. the same model, but with unknown parameter values that
must be estimated during the game from observed data, 3. models that are
structurally different from those used to simulate the population dynamics and
4. a model that ignores age structure. Humans on average performed much worse
than the models in cases 1 - 3. When the models ignored age structure, they
generated poorly performing management decisions, but still outperformed
students using experience and judgment 66 percent of the time.


  Abiotic factors are fundamental drivers of the dynamics of wild marine fish
populations. Identifying and quantifying their influence on species targeted by
the fishing industry is difficult and very important for managing fisheries in
a changing climate. Using multiple regression, we investigated the influence of
both temperature and rainfall on the variability of recruitment of a tropical
species, the brown tiger prawn (Penaeus esculentus), in Moreton Bay which is
located near the southern limit of its distribution on the east coast of
Australia. A step-wise selection between environmental variables identified
that variations in recruitment from 1990 to 2014 were best explained by a
combination of temperature and spawning stock biomass. Temperature explains 35%
of recruitment variability and spawning stock biomass 33%. This analysis
suggests that increasing temperatures have increased recruitment of brown tiger
prawn in Moreton Bay.


  In order to increase their robustness against environmental fluctuations,
many biological populations have developed bet-hedging mechanisms in which the
population `bets' against the presence of prolonged favorable environmental
conditions by having a few individual behaving as if they sensed a threatening
or stressful environment. As a result, the population (as a whole) increases
its chances of surviving environmental fluctuations in the long term, while
sacrificing short-term performance. In this paper, we propose a theoretical
framework, based on Markov jump linear systems, to model and evaluate the
performance of bet-hedging strategies in the presence of stochastic
fluctuations. We illustrate our results using numerical simulations.


  Background As more regions approach malaria elimination, understanding how
different interventions interact to reduce transmission becomes critical. The
Lake Kariba area of Southern Province, Zambia, is part of a multi-country
elimination effort and presents a particular challenge as it is an
interconnected region of variable transmission intensities.
  Methods In 2012-13, six rounds of mass-screen-and-treat drug campaigns were
carried out in the Lake Kariba region. A spatial dynamical model of malaria
transmission in the Lake Kariba area, with transmission and climate modeled at
the village scale, was calibrated to the 2012-13 prevalence survey data, with
case management rates, insecticide-treated net usage, and drug campaign
coverage informed by surveillance. The model was used to simulate the effect of
various interventions implemented in 2014-22 on reducing regional transmission,
achieving elimination by 2022, and maintaining elimination through 2028.
  Findings The model captured the spatio-temporal trends of decline and rebound
in malaria prevalence in 2012-13 at the village scale. Simulations predicted
that elimination required repeated mass drug administrations coupled with
simultaneous increase in net usage. Drug campaigns targeted only at high-burden
areas were as successful as campaigns covering the entire region.
  Interpretation Elimination in the Lake Kariba region is possible through
coordinating mass drug campaigns with high-coverage vector control. Targeting
regional hotspots is a viable alternative to global campaigns when human
migration within an interconnected area is responsible for maintaining
transmission in low-burden areas.


  Inferring gene regulatory networks is an important problem in systems
biology. However, these networks can be hard to infer from experimental data
because of the inherent variability in biological data as well as the large
number of genes involved. We propose a fast, simple method for inferring
regulatory relationships between genes from knockdown experiments in the NIH
LINCS dataset by calculating posterior probabilities, incorporating prior
information. We show that the method is able to find previously identified
edges from TRANSFAC and JASPAR and discuss the merits and limitations of this
approach.


  A modified, deterministic SEIR model is developed for the 2014 Ebola epidemic
occurring in the West African nations of Guinea, Liberia, and Sierra Leone. The
model describes the dynamical interaction of susceptible and infected
populations, while accounting for the effects of hospitalization and the spread
of disease through interactions with deceased, but infectious, individuals.
Using data from the World Health Organization (WHO), parameters within the
model are fit to recent estimates of infected and deceased cases from each
nation. The model is then analyzed using these parameter values. Finally,
several metrics are proposed to determine which of these nations is in greatest
need of additional resources to combat the spread of infection. These include
local and global sensitivity metrics of both the infected population and the
basic reproduction number with respect to rates of hospitalization and proper
burial.


  We consider ecological difference equations of the form $X_{t+1}^i =X_t^i
A_i(X_t)$ where $X_t^i$ is a vector of densities corresponding to the
subpopulations of species $i$ (e.g. subpopulations of different ages or living
in different patches), $X_t=(X_t^1,X_t^2,\dots,X_t^m)$ is state of the entire
community, and $A_i(X_t)$ are matrices determining the update rule for species
$i$. These equations are permanent if they are dissipative and the extinction
set $\{X: \prod_i \|X^i\|=0\}$ is repelling. If permanence persists under
perturbations of the matrices $A_i(X)$, the equations are robustly permanent.
We provide sufficient and necessary conditions for robust permanence in terms
of Lyapunov exponents for invariant measures supported by the extinction set.
Applications to ecological and epidemiological models are given.


  We study a modified record process where the $k$'th record in a series of
independent and identically distributed random variables is defined recursively
through the condition $Y_k > Y_{k-1} - \delta_{k-1}$ with a deterministic
sequence $\delta_k > 0$ called the handicap. For constant $\delta_k \equiv
\delta$ and exponentially distributed random variables it has been shown in
previous work that the process displays a phase transition as a function of
$\delta$ between a normal phase where the mean record value increases
indefinitely and a stationary phase where the mean record value remains bounded
and a finite fraction of all entries are records (Park \textit{et al} 2015 {\it
Phys. Rev.} E \textbf{91} 042707). Here we explore the behavior for general
probability distributions and decreasing and increasing sequences $\delta_k$,
focusing in particular on the case when $\delta_k$ matches the typical spacing
between subsequent records in the underlying simple record process without
handicap. We find that a continuous phase transition occurs only in the
exponential case, but a novel kind of first order transition emerges when
$\delta_k$ is increasing. The problem is partly motivated by the dynamics of
evolutionary adaptation in biological fitness landscapes, where $\delta_k$
corresponds to the change of the deterministic fitness component after $k$
mutational steps. The results for the record process are used to compute the
mean number of steps that a population performs in such a landscape before
being trapped at a local fitness maximum.


  Workers in insect societies are sometimes observed to kill male eggs of other
workers, a phenomenon known as worker policing. We perform a mathematical
analysis of the evolutionary dynamics of policing. We investigate the selective
forces behind policing for both dominant and recessive mutations for different
numbers of matings of the queen. The traditional, relatedness-based argument
suggests that policing evolves if the queen mates with more than two males, but
does not evolve if the queen mates with a single male. We derive precise
conditions for the invasion and stability of policing alleles. We find that the
relatedness-based argument is not robust with respect to small changes in
colony efficiency caused by policing. We also calculate evolutionarily singular
strategies and determine when they are evolutionarily stable. We use a
population genetics approach that applies to dominant or recessive mutations of
any effect size.


  Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), linked
with CRISPR associated (CAS) genes, play a profound role in the interactions
between phage and their bacterial hosts. It is now well understood that
CRISPR-CAS systems can confer adaptive immunity against bacteriophage
infections. However, the possibility of failure of CRISPR immunity may lead to
a productive infection by the phage (cell lysis) or lysogeny. Recently,
CRISPR-CAS genes have been implicated in changes to group behaviour, including
biofilm formation, of the bacterium Pseudomonas aeruginosa when lysogenized.
For lysogens with a CRISPR system, another recent experimental study suggests
that bacteriophage re-infection of previously lysogenized bacteria may lead to
cell death. Thus CRISPR immunity can have complex effects on phage-host-lysogen
interactions, particularly in a biofilm. In this contribution, we develop and
analyse a series of models to elucidate and disentangle these interactions.
From a therapeutic standpoint, CRISPR immunity increases biofilm resistance to
phage therapy. Our models predict that lysogens may be able to displace
CRISPR-immune bacteria in a biofilm, and thus suggest strategies to eliminate
phage resistant biofilms.


  Cells may control fluctuations in protein levels by means of negative
autoregulation, where transcription factors bind DNA sites to repress their own
production. Theoretical studies have assumed a single binding site for the
repressor, while in most species it is found that multiple binding sites are
arranged in clusters. We study a stochastic description of negative
autoregulation with multiple binding sites for the repressor. We find that
increasing the number of binding sites induces regular bursting of gene
products. By tuning the threshold for repression, we show that multiple binding
sites can also suppress fluctuations. Our results highlight possible roles for
the presence of multiple binding sites of negative autoregulators.


  Distributions of abundances or frequencies play an important role in many
fields of science, from biology to sociology, as does the R\'enyi entropy,
which measures the diversity of a statistical ensemble. We derive a
mathematical relation between the abundance distribution and the R\'enyi
entropy, by analogy with the equivalence of ensembles in thermodynamics. The
abundance distribution is mapped onto the density of states, and the R\'enyi
entropy to the free energy. The two quantities are related in the thermodynamic
limit by a Legendre transform, by virtue of the equivalence between the
micro-canonical and canonical ensembles. In this limit, we show how the R\'enyi
entropy can be constructed geometrically from rank-frequency plots. This
mapping predicts that non-concave regions of the rank-frequency curve should
result in kinks in the R\'enyi entropy as a function of its order. We
illustrate our results on simple examples, and emphasize the limitations of the
equivalence of ensembles when a thermodynamic limit is not well defined. Our
results help choose reliable diversity measures based on the experimental
accuracy of the abundance distributions in particular frequency ranges.


  During the lifecycle of a virus, viral proteins and other components
self-assemble to form a symmetric protein shell called a capsid. This assembly
process is subject to multiple competing constraints, including the need to
form a thermostable shell while avoiding kinetic traps. It has been proposed
that viral assembly satisfies these constraints through allosteric regulation,
including the interconversion of capsid proteins among conformations with
different propensities for assembly. In this article we use computational and
theoretical modeling to explore how such allostery affects the assembly of
icosahedral shells. We simulate assembly under a wide range of protein
concentrations, protein binding affinities, and two different mechanisms of
allosteric control. We find that, above a threshold strength of allosteric
control, assembly becomes robust over a broad range of subunit binding
affinities and concentrations, allowing the formation of highly thermostable
capsids. Our results suggest that allostery can significantly shift the range
of protein binding affinities that lead to successful assembly, and thus should
be accounted for in models that are used to estimate interaction parameters
from experimental data.


  We propose a nonlocal Poisson-Fermi model for ionic solvent that includes ion
size effects and polarization correlations among water molecules in the
calculation of electrostatic potential. It includes the previous Poisson-Fermi
models as special cases, and its solution is the convolution of a solution of
the corresponding nonlocal Poisson dielectric model with a Yukawa-type kernel
function. Moreover, the Fermi distribution is shown to be a set of optimal
ionic concentration functions in the sense of minimizing an electrostatic
potential free energy. Finally, numerical results are reported to show the
difference between a Poisson-Fermi solution and a corresponding Poisson
solution.


  We present a multiscale approach to model diffusion in a crowded environment
and its effect on the reaction rates. Diffusion in biological systems is often
modeled by a discrete space jump process in order to capture the inherent noise
of biological systems, which becomes important in the low copy number regime.
To model diffusion in the crowded cell environment efficiently, we compute the
jump rates in this mesoscopic model from local first exit times, which account
for the microscopic positions of the crowding molecules, while the diffusing
molecules jump on a coarser Cartesian grid. We then extract a macroscopic
description from the resulting jump rates, where the excluded volume effect is
modeled by a diffusion equation with space dependent diffusion coefficient. The
crowding molecules can be of arbitrary shape and size and numerical experiments
demonstrate that those factors together with the size of the diffusing molecule
play a crucial role on the magnitude of the decrease in diffusive motion. When
correcting the reaction rates for the altered diffusion we can show that
molecular crowding either enhances or inhibits chemical reactions depending on
local fluctuations of the obstacle density.


  This paper provides an overview and a tutorial of molecular clock dating
using MrBayes, which is a software for Bayesian inference of phylogeny. Two
modern approaches, total-evidence dating and node dating, are demonstrated
using a dataset of Hymenoptera with molecular sequences and morphological
characters. The similarity and difference of the two methods are compared and
discussed. Besides, a non-clock analysis is performed on the same dataset to
compare with the molecular clock dating analyses.


  We present a mathematical analysis of the early detection of Ebola virus. The
propagation of the virus is analysed by using a Susceptible, Infected,
Recovered (SIR) model. In order to provide useful predictions about the
potential transmission of the virus, we analyse and simulate the SIR model with
vital dynamics, by adding demographic effects and an induced death rate. Then,
we compute the equilibria of the model. The numerical simulations confirm the
theoretical analysis. Our study describes the 2015 detection of Ebola virus in
Guinea, the parameters of the model being identified from the World Health
Organization data. Finally, we consider an optimal control problem of the
propagation of the Ebola virus, minimizing the number of infected individuals
while taking into account the cost of vaccination.


  Next-generation RNA sequencing (RNA-seq) technology has been widely used to
assess full-length RNA isoform abundance in a high-throughput manner. RNA-seq
data offer insight into gene expression levels and transcriptome structures,
enabling us to better understand the regulation of gene expression and
fundamental biological processes. Accurate isoform quantification from RNA-seq
data is challenging due to the information loss in sequencing experiments. A
recent accumulation of multiple RNA-seq data sets from the same tissue or cell
type provides new opportunities to improve the accuracy of isoform
quantification. However, existing statistical or computational methods for
multiple RNA-seq samples either pool the samples into one sample or assign
equal weights to the samples when estimating isoform abundance. These methods
ignore the possible heterogeneity in the quality of different samples and could
result in biased and unrobust estimates. In this article, we develop a method,
which we call "joint modeling of multiple RNA-seq samples for accurate isoform
quantification" (MSIQ), for more accurate and robust isoform quantification by
integrating multiple RNA-seq samples under a Bayesian framework. Our method
aims to (1) identify a consistent group of samples with homogeneous quality and
(2) improve isoform quantification accuracy by jointly modeling multiple
RNA-seq samples by allowing for higher weights on the consistent group. We show
that MSIQ provides a consistent estimator of isoform abundance, and we
demonstrate the accuracy and effectiveness of MSIQ compared with alternative
methods through simulation studies on D. melanogaster genes. We justify MSIQ's
advantages over existing approaches via application studies on real RNA-seq
data from human embryonic stem cells, brain tissues, and the HepG2 immortalized
cell line.


  The now classical replicator equation describes a wide variety of biological
phenomena, including those in theoretical genetics, evolutionary game theory,
or in the theories of the origin of life. Among other questions, the permanence
of the replicator equation is well studied in the local, well-mixed case.
Inasmuch as the spatial heterogeneities are key to understanding the species
coexistence at least in some cases, it is important to supplement the classical
theory of the non-distributed replicator equation with a spatially explicit
framework. One possible approach, motivated by the porous medium equation, is
introduced. It is shown that the solutions to the spatially heterogeneous
replicator equation may evolve to equilibrium states that have a bounded
support, and, moreover, that these solutions are of paramount importance for
the overall system permanence, which is shown to be a more commonplace
phenomenon for the spatially explicit equation if compared with the local
model.


  Summary: The development of automated servers to predict the
three-dimensional structure of proteins has seen much progress over the years.
These servers make modeling simpler, but largely exclude users from the
process. We present an Interactive Modeling Pipeline (IMPi) for homology
modeling. The pipeline simplifies the modeling process and reduces the workload
required by the user, while still allowing engagement from the user during
every step. Default parameters are given for each step, which can either be
modified or supplemented with additional external input. As such, it has been
designed for users of varying levels of experience with homology modeling.
  Availability and implementation: The IMPi site is free for non-commercial use
and can be accessed at https://impi.rubi.ru.ac.za.
  Supplementary information: Documentation available at
https://impi.rubi.ru.ac.za/#documentation.


  Recently, we have proposed a new illness-death model that comprises a state
of undiagnosed chronic disease preceding the diagnosed disease. Based on this
model, the question arises how case-finding can be assessed in the presence of
mortality from all these states. We simulate two scenarios of different
performance of case-finding and apply several indices to assess case-finding in
both scenarios. One of the prevalence based indices leads to wrong conclusions.
Some indices are partly insensitive to distinguish the quality of case-finding.
The incidence based indices perform well. If possible, incidence based indices
should be preferred.


  Self-assembly of proteins is a biological phenomenon which gives rise to
spontaneous formation of amyloid fibrils or polymers. The starting point of
this phase, called nucleation exhibits an important variability among
replicated experiments.To analyse the stochastic nature of this phenomenon, one
of the simplest models considers two populations of chemical components:
monomers and polymerised monomers. Initially there are only monomers. There are
two reactions for the polymerization of a monomer: either two monomers collide
to combine into two polymerised monomers or a monomer is polymerised after the
encounter of a polymerised monomer. It turns out that this simple model does
not explain completely the variability observed in the experiments. This paper
investigates extensions of this model to take into account other mechanisms of
the polymerization process that may have impact an impact on fluctuations.The
first variant consists in introducing a preliminary conformation step to take
into account the biological fact that, before being polymerised, a monomer has
two states, regular or misfolded. Only misfolded monomers can be polymerised so
that the fluctuations of the number of misfolded monomers can be also a source
of variability of the number of polymerised monomers. The second variant
represents the reaction rate $\alpha$ of spontaneous formation of a polymer as
of the order of $N^{-\nu}$, with $\nu$ some positive constant. First and second
order results for the starting instant of nucleation are derived from these
limit theorems. The proofs of the results rely on a study of a stochastic
averaging principle for a model related to an Ehrenfest urn model, and also on
a scaling analysis of a population model.


  Coevolutionary games cast players that may change their strategies as well as
their networks of interaction. In this paper a framework is introduced for
describing coevolutionary game dynamics by landscape models. It is shown that
coevolutionary games invoke dynamic landscapes. Numerical experiments are shown
for a prisoner's dilemma (PD) and a snow drift (SD) game that both use either
birth-death (BD) or death-birth (DB) strategy updating. The resulting
landscapes are analyzed with respect to modality and ruggedness


  We study a simple model of a foraging animal that modifies the substrate on
which it moves. This substrate provides its only resource, and the forager
manage it by taking a limited portion at each visited site. The resource
recovers its value after the visit following a relaxation law. We study
different scenarios to analyze the efficiency of the managing strategy,
corresponding to control the bite size. We observe the non trivial emergence of
a home range, that is visited in a periodic way. The duration of the
corresponding cycles and the transient until it emerges is affected by the bite
size. Our results show that the most efficient use of the resource, measured as
the balance between gathering and travelled distance, corresponds to foragers
that take larger portions but without exhausting the resource. We also analyze
the use of space determining the number of attractors of the dynamics, and we
observe that it depends on the bite size and the recovery time of the resource.


  In this paper, we propose and analyze a new fractional order malaria model
with temporary immunity. We extend the existing malaria model to include
fractional derivatives. Basic reproduction number is computed using the next
generation matrix method. Model equilibria are determined and their local
asymptotic stability analysis are considered. An efficient Adams-type
predictor-corrector method is used to solve the initial value problem. Finally,
phase plane portraits are generated to describe the dynamics between
susceptible fractions of human hosts and infected fractions of mosquito
vectors.


  In the era of big data, transformation of biomedical big data into valuable
knowledge has been one of the most important challenges in bioinformatics. Deep
learning has advanced rapidly since the early 2000s and now demonstrates
state-of-the-art performance in various fields. Accordingly, application of
deep learning in bioinformatics to gain insight from data has been emphasized
in both academia and industry. Here, we review deep learning in bioinformatics,
presenting examples of current research. To provide a useful and comprehensive
perspective, we categorize research both by the bioinformatics domain (i.e.,
omics, biomedical imaging, biomedical signal processing) and deep learning
architecture (i.e., deep neural networks, convolutional neural networks,
recurrent neural networks, emergent architectures) and present brief
descriptions of each study. Additionally, we discuss theoretical and practical
issues of deep learning in bioinformatics and suggest future research
directions. We believe that this review will provide valuable insights and
serve as a starting point for researchers to apply deep learning approaches in
their bioinformatics studies.


  Site-directed mutagenesis refers to a man-made molecular biology method that
is used to make genetic alterations in the DNA sequence of a gene of interest.
But based on our recently published experimental findings, we propose that
natural site-directed mutagenesis might exist in the eukaryotic cells, which is
triggered by harmful agents and co-directed by special transcription hotspots
and mutation-contained intranuclear primers.


  Mutation rates and fitness costs of deleterious mutations are difficult to
measure in vivo but essential for a quantitative understanding of evolution.
Using whole genome deep sequencing data from longitudinal samples during
untreated HIV-1 infection, we estimated mutation rates and fitness costs in
HIV-1 from the temporal dynamics of genetic variation. At approximately neutral
sites, mutations accumulate with a rate of 1.2 x 10^-5 per site per day, in
agreement with the rate measured in cell cultures. The rate from G to A is
largest, followed by the other transitions C to T, T to C, and A to G, while
transversions are more rare. At non-neutral sites, most mutations reduce virus
replication; using a model of mutation selection balance, we estimated the
fitness cost of mutations at every site in the HIV-1 genome. About half of all
nonsynonymous mutations have large fitness costs (greater than 10\%), while
most synonymous mutations have costs below 1\%. The cost of synonymous
mutations is especially low in most of gag and pol, while much higher costs are
observed in important RNA structures and regulatory regions. The intrapatient
fitness cost estimates are consistent across multiple patients, suggesting that
the deleterious part of the fitness landscape is universal and explains a large
fraction of global HIV-1 group M diversity.


  By deploying myelofibrosis as the disease context, I wish to propose that
increased availability of Tmp21 (an NFAT gene target) induces aberrant protein
secretion from the ER contributing to pathological consequences, which has not
been elucidated before. Primary myelofibrosis is now mainly considered as an
advanced stage of BCR-ABL1 negative myeloproliferative neoplasms (MPN), which
otherwise include polycythemia vera and essential thrombocythemia.
Myelofibrosis is defined by an increased insoluble collagen fiber deposition in
the bone marrow2 and harbors chronic inflammation as an important component in
disease progression.


  Rapid development of modern sequencing platforms enabled an unprecedented
growth of protein families databases. The abundance of sets composed of
hundreds of thousands sequences is a great challenge for multiple sequence
alignment algorithms. In the article we introduce FAMSA, a new progressive
algorithm designed for fast and accurate alignment of thousands of protein
sequences. Its features include the utilisation of longest common subsequence
measure for determining pairwise similarities, a novel method of gap costs
evaluation, and a new iterative refinement scheme. Importantly, its
implementation is highly optimised and parallelised to make the most of modern
computer platforms. Thanks to the above, quality indicators, namely
sum-of-pairs and total-column scores, show FAMSA to be superior to competing
algorithms like Clustal Omega or MAFFT for datasets exceeding a few thousand of
sequences. The quality does not compromise time and memory requirements which
are an order of magnitude lower than that of existing solutions. For example, a
family of 415 519 sequences was analysed in less than two hours and required
only 8GB of RAM.
  FAMSA is freely available at http://sun.aei.polsl.pl/REFRESH/famsa.


  Bacteria communicate using external chemical signals called autoinducers (AI)
in a process known as quorum sensing (QS). QS efficiency is reduced by both
limitations of AI diffusion and potential interference from neighboring
strains. There is thus a need for predictive theories of how spatial community
structure shapes information processing in complex microbial ecosystems. As a
step in this direction, we apply a reaction-diffusion model to study
autoinducer signaling dynamics in a single-species community as a function of
the spatial distribution of colonies in the system. We predict a dynamical
transition between a local quorum sensing (LQS) regime, with the AI signaling
dynamics primarily controlled by the local population densities of individual
colonies, and a global quorum sensing (GQS) regime, with the dynamics being
dependent on collective inter-colony diffusive interactions. The crossover
between LQS to GQS is intimately connected to a tradeoff between the signaling
network's latency, or speed of activation, and its throughput, or the total
spatial range over which all the components of the system communicate.


  Proteomic matrix-assisted laser desorption/ionisation (MALDI) linear
time-of-flight (TOF) mass spectrometry (MS) may be used to produce protein
profiles from biological samples with the aim of discovering biomarkers for
disease. However, the raw protein profiles suffer from several sources of bias
or systematic variation which need to be removed via pre-processing before
meaningful downstream analysis of the data can be undertaken. Baseline
subtraction, an early pre-processing step that removes the non-peptide signal
from the spectra, is complicated by the following: (i) each spectrum has, on
average, wider peaks for peptides with higher mass-to-charge ratios (m/z), and
(ii) the time-consuming and error-prone trial-and-error process for optimising
the baseline subtraction input arguments. With reference to the aforementioned
complications, we present an automated pipeline that includes (i) a novel
`continuous' line segment algorithm that efficiently operates over data with a
transformed m/z-axis to remove the relationship between peptide mass and peak
width, and (ii) an input-free algorithm to estimate peak widths on the
transformed m/z scale. The automated baseline subtraction method was deployed
on six publicly available proteomic MS datasets using six different m/z-axis
transformations. Optimality of the automated baseline subtraction pipeline was
assessed quantitatively using the mean absolute scaled error (MASE) when
compared to a gold-standard baseline subtracted signal. Near-optimal baseline
subtraction was achieved using the automated pipeline. The advantages of the
proposed pipeline include informed and data specific input arguments for
baseline subtraction methods, the avoidance of time-intensive and subjective
piecewise baseline subtraction, and the ability to automate baseline
subtraction completely. Moreover, individual steps can be adopted as
stand-alone routines.


  Protein solubility and conformational stability are a result of a balance of
interactions both within a protein and between protein and solvent. The
electrostatic solvation free energy of oligoglycines, models for the peptide
backbone, becomes more favorable with an increasing length, yet longer peptides
collapse due to the formation of favorable intrapeptide interactions between CO
dipoles, in some cases without hydrogen bonds. The strongly repulsive solvent
cavity formation is balanced by van der Waals attractions and electrostatic
contributions. In order to investigate the competition between solvent
exclusion and charge interactions we simulate the collapse of a long
oligoglycine comprised of 15 residues while scaling the charges on the peptide
from zero to fully charged. We examine the effect this has on the
conformational properties of the peptide. We also describe the approximate
thermodynamic changes that occur during the scaling both in terms of
intrapeptide potentials and peptide-water potentials, and estimate the
electrostatic solvation free energy of the system.


  Large complexes of classical particles play central roles in biology, in
polymer physics, and in other disciplines. However, physics currently lacks
mathematical methods for describing such complexes in terms of component
particles, interaction energies, and assembly rules. Here we describe a Fock
space structure that addresses this need, as well as diagrammatic methods that
facilitate the use of this formalism. These methods can dramatically simplify
the equations governing both equilibrium and non-equilibrium stochastic
chemical systems. A mathematical relationship between the set of all complexes
and a list of rules for complex assembly is also identified.


  Known as one of the hallmarks of cancer [30], cancer cell invasion of human
body tissue is a complicated spatio-temporal multiscale process which enables a
localised solid tumour to transform into a systemic, metastatic and fatal
disease. This process explores and takes advantage of the reciprocal relation
that solid tumours establish with the extracellular matrix (ECM) components and
other multiple distinct cell types from the surrounding microenvironment.
Through the secretion of various proteolytic enzymes such as matrix
metalloproteinases (MMP) or the urokinase plasminogen activator (uPA), the
cancer cell population alters the configuration of the surrounding ECM
composition and overcomes the physical barriers to ultimately achieve local
cancer spread into the surrounding tissue. The active interplay between the
tissue-scale tumour dynamics and the molecular mechanics of the involved
proteolytic enzymes at the cell-scale underlines the biologically multiscale
character of invasion, and raises the challenge of modelling this process with
an appropriate multiscale approach. In this paper, we present a new two-scale
moving boundary model of cancer invasion that explores the tissue scale tumour
dynamics in conjunction with the molecular dynamics of the urokinase
plasminogen activation system. Building on the multiscale moving boundary
method proposed in [58], the modelling that we propose here allows us to study
the changes in tissue scale tumour morphology caused by the cell-scale uPA
micro-dynamics occurring along the invasive edge of the tumour. Our
computational simulation results demonstrate a range of heterogeneous dynamics
which are qualitatively similar to the invasive growth patterns observed in a
number of different types of cancer, such as the tumour infiltrative growth
patterns discussed in [33].


  We study a model of a population making a binary decision based on
information spreading within the population, which is fully connected or
covering a square grid. We assume that a fraction of the population wants to
make the choice of the majority, while the rest want to make the minority
choice. This resembles opinion spreading with "contrarian" agents, but has the
game theoretic aspect that agents try to optimize their own situation in ways
that are incompatible with the common good. When this fraction is less than
1/2, the population can efficiently self-organize to a state where agents get
what they want -- the majority (i.e. the majority seekers) have one opinion,
the minority seekers have the other. If the fraction is larger than 1/2, there
is a frustration in the population that dramatically changes the dynamics. In
this region, the population converges, through some distinct phases, to a state
of approximately equal-sized opinions. Just over the threshold the state of the
population is furthest from the collectively optimal solution.


  An innovative theoretical framework for stochastic dynamics based on a
decomposition of a stochastic differential equation (SDE) has been developed
with an evident advantage in connecting deterministic and stochastic dynamics,
as well as useful applications in physics, engineering, chemistry and biology.
It introduces the A-type stochastic integration for SDE beyond traditional
Ito's or Stratonovich's interpretation. Serious question on its uniqueness was
recently raised. We provide here both mathematical and physical demonstrations
that the uniqueness is guaranteed. Such discussion leads to a better
understanding on the robustness of the novel framework. We also discuss the
limitation of a related approach of obtaining potential function from steady
state distribution.


  Intracellular transport is an essential function in eucaryotic cells,
facilitated by motor proteins - proteins converting chemical energy into
kinetic energy. It is known that motor proteins work in teams enabling
unidirectional and bidirectional transport of intracellular cargo over long
distances. Disruptions of the underlying transport mechanisms, often caused by
mutations that alter single motor characteristics, are known to cause
neurodegenerative diseases. For example, phosphorylation of kinesin motor
domain at the serine residue is implicated in Huntington's disease, with a
recent study of phosphorylated and phosphomimetic serine residues indicating
lowered single motor stalling forces. In this article we report the effects of
mutations of this nature on transport properties of cargo carried by multiple
wild-type and mutant motors. Results indicate that mutants with altered stall
forces might determine the average velocity and run-length even when they are
outnumbered by wild type motors in the ensemble. It is shown that mutants gain
a competitive advantage and lead to an increase in expected run-length when
load on the cargo is in the vicinity of the mutant's stalling force or a
multiple of its stalling force. A separate contribution of this article is the
development of a semi-analytic method to analyze transport of cargo by multiple
motors of multiple types. The technique determines transition rates between
various relative configurations of motors carrying the cargo using the
transition rates between various absolute configurations. This enables exact
computation of average velocity and run-length. It can also be used to
introduce alterations of various single motor parameters to model a mutation
and to deduce effects of such alterations on the transport of a common cargo by
multiple motors. Our method is easily implementable and we provide a software
package for general use.


  The human immune system depends on a highly diverse collection of
antibody-making B cells. B cell receptor sequence diversity is generated by a
random recombination process called "rearrangement" forming progenitor B cells,
then a Darwinian process of lineage diversification and selection called
"affinity maturation." The resulting receptors can be sequenced in high
throughput for research and diagnostics. Such a collection of sequences
contains a mixture of various lineages, each of which may be quite numerous, or
may consist of only a single member. As a step to understanding the process and
result of this diversification, one may wish to reconstruct lineage membership,
i.e. to cluster sampled sequences according to which came from the same
rearrangement events. We call this clustering problem "clonal family
inference." In this paper we describe and validate a likelihood-based framework
for clonal family inference based on a multi-hidden Markov Model (multi-HMM)
framework for B cell receptor sequences. We describe an agglomerative algorithm
to find a maximum likelihood clustering, two approximate algorithms with
various trade-offs of speed versus accuracy, and a third, fast algorithm for
finding specific lineages. We show that under simulation these algorithms
greatly improve upon existing clonal family inference methods, and that they
also give significantly different clusters than previous methods when applied
to two real data sets.


  Motivation: DNA methylation is an intensely studied epigenetic mark, yet its
functional role is incompletely understood. Attempts to quantitatively
associate average DNA methylation to gene expression yield poor correlations
outside of the well-understood methylation-switch at CpG islands.
  Results: Here we use probabilistic machine learning to extract higher order
features associated with the methylation profile across a defined region. These
features quantitate precisely notions of shape of a methylation profile,
capturing spatial correlations in DNA methylation across genomic regions. Using
these higher order features across promoter-proximal regions, we are able to
construct a powerful machine learning predictor of gene expression,
significantly improving upon the predictive power of average DNA methylation
levels. Furthermore, we can use higher order features to cluster
promoter-proximal regions, showing that five major patterns of methylation
occur at promoters across different cell lines, and we provide evidence that
methylation beyond CpG islands may be related to regulation of gene expression.
Our results support previous reports of a functional role of spatial
correlations in methylation patterns, and provide a mean to quantitate such
features for downstream analyses.
  Availability: https://github.com/andreaskapou/BPRMeth


  Growth in static and controlled environments such as a Petri dish can be used
to study the spatial population dynamics of microorganisms. However, natural
populations such as marine microbes experience fluid advection and often grow
up in heterogeneous environments. We investigate a generalized
Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation describing single species
population subject to a constant flow field and quenched random spatially
inhomogeneous growth rates with a fertile overall growth condition. We
analytically and numerically demonstrate that the non-equilibrium steady-state
population density develops a flow-driven striation pattern. The striations are
highly asymmetric with a longitudinal correlation length that diverges linearly
with the flow speed and a transverse correlation length that approaches a
finite velocity-independent value. Linear response theory is developed to study
the statistics of the steady states. Theoretical predictions show excellent
agreement with the numerical steady states of the generalized FKPP equation
obtained from Lattice Boltzmann simulations. These findings suggest that,
although the growth disorder can be spatially uncorrelated, correlated
population structures with striations emerge naturally at sufficiently strong
advection.


  In phylogenetics, phylogenetic trees are rooted binary trees, whereas
phylogenetic networks are rooted arbitrary acyclic digraphs. Edges are directed
away from the root and leaves are uniquely labeled with taxa in phylogenetic
networks. For the purpose of validating evolutionary models, biologists check
whether or not a phylogenetic tree is contained in a phylogenetic network on
the same taxa. This tree containment problem is known to be NP-complete. A
phylogenetic network is reticulation-visible if every reticulation node
separates the root of the network from some leaves. We answer an open problem
by proving that the problem is solvable in quadratic time for
reticulation-visible networks. The key tool used in our answer is a powerful
decomposition theorem. It also allows us to design a linear-time algorithm for
the cluster containment problem for networks of this type and to prove that
every galled network with n leaves has 2(n-1) reticulation nodes at most.


  Reconstructing ancestral gene orders in a given phylogeny is a classical
problem in comparative genomics. Most existing methods compare conserved
features in extant genomes in the phylogeny to define potential ancestral gene
adjacencies, and either try to reconstruct all ancestral genomes under a global
evolutionary parsimony criterion, or, focusing on a single ancestral genome,
use a scaffolding approach to select a subset of ancestral gene adjacencies,
generally aiming at reducing the fragmentation of the reconstructed ancestral
genome. In this paper, we describe an exact algorithm for the Small Parsimony
Problem that combines both approaches. We consider that gene adjacencies at
internal nodes of the species phylogeny are weighted, and we introduce an
objective function defined as a convex combination of these weights and the
evolutionary cost under the Single-Cut-or-Join (SCJ) model. The weights of
ancestral gene adjacencies can e.g. be obtained through the recent availability
of ancient DNA sequencing data, which provide a direct hint at the genome
structure of the considered ancestor, or through probabilistic analysis of gene
adjacencies evolution. We show the NP-hardness of our problem variant and
propose a Fixed-Parameter Tractable algorithm based on the Sankoff-Rousseau
dynamic programming algorithm that also allows to sample co-optimal solutions.
We apply our approach to mammalian and bacterial data providing different
degrees of complexity. We show that including adjacency weights in the
objective has a significant impact in reducing the fragmentation of the
reconstructed ancestral gene orders.


  We introduce a method for predicting RNA folding pathways, with an
application to the most important RNA tetraloops. The method is based on the
idea that ensembles of three-dimensional fragments extracted from
high-resolution crystal structures are heterogeneous enough to describe
metastable as well as intermediate states. These ensembles are first validated
by performing a quantitative comparison against available solution NMR data of
a set of RNA tetranucleotides. Notably, the agreement is better with respect to
the one obtained by comparing NMR with extensive all-atom molecular dynamics
simulations. We then propose a procedure based on diffusion maps and Markov
models that makes it possible to obtain reaction pathways and their relative
probabilities from fragment ensembles. This approach is applied to study the
helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG
families. The results give detailed insights into the folding mechanism that
are compatible with available experimental data and clarify the role of
intermediate states observed in previous simulation studies. The method is
computationally inexpensive and can be used to study arbitrary conformational
transitions.


  We consider adaptive change of diet of a predator population that switches
its feeding between two prey populations. We develop a novel 1 fast--3 slow
dynamical system to describe the dynamics of the three populations amidst
continuous but rapid evolution of the predator's diet choice. The two extremes
at which the predator's diet is composed solely of one prey correspond to two
branches of the three-branch critical manifold of the fast--slow system. By
calculating the points at which there is a fast transition between these two
feeding choices (i.e., branches of the critical manifold), we prove that the
system has a two-parameter family of periodic orbits for sufficiently large
separation of the time scales between the evolutionary and ecological dynamics.
Using numerical simulations, we show that these periodic orbits exist, and that
their phase difference and oscillation patterns persist, when ecological and
evolutionary interactions occur on comparable time scales. Our model also
exhibits periodic orbits that agree qualitatively with oscillation patterns
observed in experimental studies of the coupling between rapid evolution and
ecological interactions.


  We investigate the effects of cooperativity between contagion processes that
spread and persist in a host population. We propose and analyze a dynamical
model in which individuals that are affected by one transmissible agent $A$
exhibit a higher than baseline propensity of being affected by a second agent
$B$ and vice versa. The model is a natural extension of the traditional SIS
(Susceptible-Infected-Susceptible) model used for modeling single contagion
processes. We show that cooperativity changes the dynamics of the system
considerably when cooperativity is strong. The system exhibits discontinuous
phase transitions not observed in single agent contagion, multi-stability, a
separation of the traditional epidemic threshold into different thresholds for
inception and extinction as well as hysteresis. These properties are robust and
are corroborated by stochastic simulations on lattices and generic network
topologies. Finally, we investigate wave propagation and transients in a
spatially extended version of the model and show that especially for
intermediate values of baseline reproduction ratios the system is characterized
by various types of wave-front speeds. The system can exhibit spatially
heterogeneous stationary states for some parameters and negative front speeds
(receding wave fronts). The two agent model can be employed as a starting point
for more complex contagion processes, involving several interacting agents, a
model framework particularly suitable for modeling the spread and dynamics of
microbiological ecosystems in host populations.


  MicroRNAs (miRNAs) are short sequences of ribonucleic acids that control the
expression of target messenger RNAs (mRNAs) by binding them. Robust prediction
of miRNA-mRNA pairs is of utmost importance in deciphering gene regulations but
has been challenging because of high false positive rates, despite a deluge of
computational tools that normally require laborious manual feature extraction.
This paper presents an end-to-end machine learning framework for miRNA target
prediction. Leveraged by deep recurrent neural networks-based auto-encoding and
sequence-sequence interaction learning, our approach not only delivers an
unprecedented level of accuracy but also eliminates the need for manual feature
extraction. The performance gap between the proposed method and existing
alternatives is substantial (over 25% increase in F-measure), and deepTarget
delivers a quantum leap in the long-standing challenge of robust miRNA target
prediction.


  Motivation: The MinION device by Oxford Nanopore is the first portable
sequencing device. MinION is able to produce very long reads (reads over
100~kBp were reported), however it suffers from high sequencing error rate. In
this paper, we show that the error rate can be reduced by improving the base
calling process.
  Results: We present the first open-source DNA base caller for the MinION
sequencing platform by Oxford Nanopore. By employing carefully crafted
recurrent neural networks, our tool improves the base calling accuracy compared
to the default base caller supplied by the manufacturer. This advance may
further enhance applicability of MinION for genome sequencing and various
clinical applications.
  Availability: DeepNano can be downloaded at
http://compbio.fmph.uniba.sk/deepnano/.
  Contact: boza@fmph.uniba.sk


  Integrating population genetics into community ecology theory is a major goal
in ecology and evolution, but analyzing the resulting models is computationally
daunting. Here we describe sPEGG ($\underline{\textrm{s}}\textrm{imulating}$
$\underline{\textrm{P}}\textrm{henotypic}$
$\underline{\textrm{E}}\textrm{volution}$ on
$\underline{\textrm{G}}\textrm{eneral Purpose}$
$\underline{\textrm{G}}\textrm{raphics Processing Units}$ (GPGPUs)), an
open-source, multi-species forward-time population genetics simulator. Using a
single commodity GPGPU instead of a single central processor, we find sPEGG can
accelerate eco-evolutionary simulations by a factor of over 200, comparable to
performance on a small-to-medium sized computer cluster.


  Most antigenically novel and evolutionarily successful strains of seasonal
influenza A (H3N2) originate in East, South, and Southeast Asia. To understand
this pattern, we simulated the ecological and evolutionary dynamics of
influenza in a host metapopulation representing the temperate north, tropics,
and temperate south. Although seasonality and air traffic are frequently used
to explain global migratory patterns of influenza, we find that other factors
may have a comparable or greater impact. Notably, a region's basic reproductive
number ($R_0$) strongly affects the antigenic evolution of its viral population
and the probability that its strains will spread and fix globally: a 17-28%
higher $R_0$ in one region can explain the observed patterns. Seasonality, in
contrast, increases the probability that a tropical (less seasonal) population
will export evolutionarily successful strains but alone does not predict that
these strains will be antigenically advanced. The relative sizes of different
host populations, their birth and death rates, and the region in which H3N2
first appears affect influenza's phylogeography in different but relatively
minor ways. These results suggest general principles that dictate the spatial
dynamics of antigenically evolving pathogens and offer predictions for how
changes in human ecology might affect influenza evolution.


  Traditionally, studies of coevolving systems have considered cases where a
parasite may inhabit only a single host. The case where a parasite may infect
many hosts, widespread parasitism, has until recently gained little traction.
This is due in part to the computational complexity involved in reconstructing
the coevolutionary histories where parasites may infect only a single host,
which is NP-Hard. Allowing parasites to inhabit more than one host has been
seen to only further compound this computationally intractable problem.
Recently however, well-established algorithms for estimating the problem
instance where a parasite may infect only a single host have been extended to
handle widespread parasites. Although this has offered significant progress, it
has been noted that these algorithms poorly handle parasites that inhabit
phylogenetically distant hosts.
  In this work we extend these previous algorithms to handle cases where
parasites inhabit phylogenetically distant hosts using an additional
evolutionary event which we call spread. Our new framework is shown to infer
significantly more congruent coevolutionary histories compared to existing
methods over both synthetic and biological data sets. We then apply the newly
proposed algorithm, which we call WiSPA (WideSpread Parasitism Analyser), to
the well studied coevolutionary system of Primates and Enterobius (pinworms),
where existing methods have been unable to reconcile the widespread parasitism
present without permitting additional divergence events. Using WiSPA and the
new biological event, spread, we provide the first statistically significant
coevolutionary hypothesis for this system.


  Recent advances in nonequilibrium statistical physics have provided
unprecedented insight into the thermodynamics of dynamic processes. The author
recently used these advances to extend Landauer's semi-formal reasoning
concerning the thermodynamics of bit erasure, to derive the minimal free energy
required to implement an arbitrary computation. Here, I extend this analysis,
deriving the minimal free energy required by an organism to run a given
(stochastic) map $\pi$ from its sensor inputs to its actuator outputs. I use
this result to calculate the input-output map $\pi$ of an organism that
optimally trades off the free energy needed to run $\pi$ with the phenotypic
fitness that results from implementing $\pi$. I end with a general discussion
of the limits imposed on the rate of the terrestrial biosphere's information
processing by the flux of sunlight on the Earth.


  Graph theory and graph mining constitute rich fields of computational
techniques to study the structures, topologies and properties of graphs. These
techniques constitute a good asset in bioinformatics if there exist efficient
methods for transforming biological data into graphs. In this paper, we present
Protein Graph Repository (PGR), a novel database of protein 3D-structures
transformed into graphs allowing the use of the large repertoire of graph
theory techniques in protein mining. This repository contains graph
representations of all currently known protein 3D-structures described in the
Protein Data Bank (PDB). PGR also provides an efficient online converter of
protein 3D-structures into graphs, biological and graph-based description,
pre-computed protein graph attributes and statistics, visualization of each
protein graph, as well as graph-based protein similarity search tool. Such
repository presents an enrichment of existing online databases that will help
bridging the gap between graph mining and protein structure analysis. PGR data
and features are unique and not included in any other protein database. The
repository is available at http://wjdi.bioinfo.uqam.ca/.


  Many areas of agriculture rely on honey bees to provide pollination services
and any decline in honey bee numbers can impact on global food security. In
order to understand the dynamics of honey bee colonies we present a discrete
time marked renewal process model for the size of a colony. We demonstrate that
under mild conditions this attains a stationary distribution that depends on
the distribution of the numbers of eggs per batch, the probability an egg
hatches and the distributions of the times between batches and bee lifetime.
This allows an analytic examination of the effect of changing these quantities.
We then extend this model to cyclic annual effects where for example the
numbers of eggs per batch and {the probability an egg hatches} may vary over
the year.