  We present a model for continuous cell culture coupling intra-cellular
metabolism to extracellular variables describing the state of the bioreactor,
taking into account the growth capacity of the cell and the impact of toxic
byproduct accumulation. We provide a method to determine the steady states of
this system that is tractable for metabolic networks of arbitrary complexity.
We demonstrate our approach in a toy model first, and then in a genome-scale
metabolic network of the Chinese hamster ovary cell line, obtaining results
that are in qualitative agreement with experimental observations. More
importantly, we derive a number of consequences from the model that are
independent of parameter values. First, that the ratio between cell density and
dilution rate is an ideal control parameter to fix a steady state with desired
metabolic properties invariant across perfusion systems. This conclusion is
robust even in the presence of multi-stability, which is explained in our model
by the negative feedback loop on cell growth due to toxic byproduct
accumulation. Moreover, a complex landscape of steady states in continuous cell
culture emerges from our simulations, including multiple metabolic switches,
which also explain why cell-line and media benchmarks carried out in batch
culture cannot be extrapolated to perfusion. On the other hand, we predict
invariance laws between continuous cell cultures with different parameters. A
practical consequence is that the chemostat is an ideal experimental model for
large-scale high-density perfusion cultures, where the complex landscape of
metabolic transitions is faithfully reproduced. Thus, in order to actually
reflect the expected behavior in perfusion, performance benchmarks of
cell-lines and culture media should be carried out in a chemostat.


  The ultrastructure fluctuations and complex dynamics of the multi-layered
membrane structure of myelin are fundamental for understanding and control its
formation process and its degeneration and repair in neurological diseases such
as multiple sclerosis (MS). Myelin is considered a liquid-crystal but
information are confined to its average structure due to limitations of the
available standard techniques. To overcome this limitation in this work we have
used Scanning micro X-ray Diffraction (S{\mu}XRD) which is a unique
non-invasive probe of both k-space and real space allowing to visualize
disorder in myelin with high spatial resolution in real space. We have used
this method to examine the myelin sheath in the sciatic nerve of Xenopus
laevis. Our results open could open new venues for understanding formation and
degradation of myelin.


  A review of the mechanisms of speciation is performed. The mechanisms of the
evolution of species, taking into account the feedback of the state of the
environment and mechanisms of the emergence of complexity, are considered. It
is shown that these mechanisms, at the molecular level, cannot work steadily in
terms of classical mechanics. Quantum mechanisms of changes in the genome,
based on the long-range interaction potential between biologically important
molecules, are proposed as one of possible explanation. Different variants of
interactions of the organism and environment based on molecular recognition and
leading to new species origins are considered. Experiments to verify the model
are proposed. This bio-physical study is completed by the general operational
model of based on quantum information theory. The latter is applied to model
epigenetic evolution.


  We classify research fields in biology into those on flows of materials,
energy, and information. As a representative energy transducing machinery in
biology, our research target, V1-ATPase from a bacterium Enterococcus hirae, a
typical molecular rotary motor is introduced. Structures of several
intermediates of the rotary motor are described and the molecular mechanism of
the motor converting chemical energy into mechanical energy is discussed.
Comments and considerations on the information flow in biology, especially on
the thermodynamic entropy in quantum physical and biological systems, are added
in a separate section containing the biologist friendly presentation of this
complicated question.


  The natural mortality (M) and purse-seine catchability and selectivity were
estimated for Trachurus novaezelandiae, Richardson, 1843 (yellowtail scad)-a
small inshore pelagic species harvested off south-eastern Australia. Hazard
functions were applied to two decades of data describing catches (mostly stable
at a mean +- SE of 315 +- 14 t p.a.) and effort (declining from a maximum of
2289 to 642 boat days between 1999/00 and 2015/16) and inter-dispersed (over
nine years) annual estimates of size-at-age (0+ to 18 years) to enable survival
analysis. The data were best described by a model with eight parameters,
including catchability (estimated at < 0.1 x 10-7 boat day-1), M (0.22 year-1)
and variable age-specific selection up to 6 years with a 50% retention among
5-year olds (larger than the estimated age at maturation). The low catchability
implied minimal fishing mortality by the purse-seine fleet. Ongoing monitoring
and applied gear-based studies are required to validate purse-seine
catchability and selectivity, but the data nevertheless imply T. novaezelandiae
could incur substantial additional fishing effort and, in doing, so alleviate
pressure on other regional small pelagics.


  Atomic radii and charges are two major parameters used in implicit solvent
electrostatics and energy calculations. The optimization problem for charges
and radii is under-determined, leading to uncertainty in the values of these
parameters and in the results of solvation energy calculations using these
parameters. This paper presents a new method for quantifying this uncertainty
in implicit solvation calculations of small molecules using surrogate models
based on generalized polynomial chaos (gPC) expansions. There are relatively
few atom types used to specify radii parameters in implicit solvation
calculations; therefore, surrogate models for these low-dimensional spaces
could be constructed using least-squares fitting. However, there are many more
types of atomic charges; therefore, construction of surrogate models for the
charge parameter space requires compressed sensing combined with an iterative
rotation method to enhance problem sparsity. We demonstrate the application of
the method by presenting results for the uncertainties in small molecule
solvation energies based on these approaches. The method presented in this
paper is a promising approach for efficiently quantifying uncertainty in a wide
range of force field parameterization problems, including those beyond
continuum solvation calculations.The intent of this study is to provide a way
for developers of implicit solvent model parameter sets to understand the
sensitivity of their target properties (solvation energy) on underlying choices
for solute radius and charge parameters.


  At equilibrium, thermodynamic and kinetic information can be extracted from
biomolecular energy landscapes by many techniques. However, while static,
ensemble techniques yield thermodynamic data, often only dynamic,
single-molecule techniques can yield the kinetic data that describes
transition-state energy barriers. Here we present a generalized framework based
upon dwell-time distributions that can be used to connect such static, ensemble
techniques with dynamic, single-molecule techniques, and thus characterize
energy landscapes to greater resolutions. We demonstrate the utility of this
framework by applying it to cryogenic electron microscopy (cryo-EM) and
single-molecule fluorescence resonance energy transfer (smFRET) studies of the
bacterial ribosomal pre-translocation complex. Among other benefits,
application of this framework to these data explains why two transient,
intermediate conformations of the pre-translocation complex, which are observed
in a cryo-EM study, may not be observed in several smFRET studies.


  We estimate the number of spacers in a CRISPR array of a bacterium which
maximizes its protection against a viral attack. The optimality follows from a
competition between two trends: too few distinct spacers make the bacteria
vulnerable to an attack by a virus with mutated corresponding protospacers,
while an excessive variety of spacers dilutes the number of the CRISPR
complexes armed with the most recent and thus most effective spacers. We first
evaluate the optimal number of spacers in a simple scenario of an infection by
a single viral species and later consider a more general case of multiple viral
species. We find that depending on such parameters as the concentration of
CRISPR-CAS interference complexes and its preference to arm with more recently
acquired spacers, the rate of viral mutation, and the number of viral species,
the predicted optimal array length lies within a range quite reasonable from
the viewpoint of recent experiments.


  A grand canonical Monte Carlo (MC) algorithm is presented for studying the
lattice gas model (LGM) of multiple protein sequence alignment, which
coherently combines long-range interactions and variable-length insertions. MC
simulations are used for both parameter optimization of the model and
production runs to explore the sequence subspace around a given protein family.
In this Note, I describe the details of the MC algorithm as well as some
preliminary results of MC simulations with various temperatures and chemical
potentials, and compare them with the mean-field approximation. The existence
of a two-state transition in the sequence space is suggested for the SH3 domain
family, and inappropriateness of the mean-field approximation for the LGM is
demonstrated.


  Background: Norovirus is a common cause of outbreaks of acute gastroenteritis
in health- and child-care settings, with serial outbreaks also frequently
observed aboard cruise ships. The relative contributions of environmental and
direct person-to-person transmission of norovirus has hitherto not been
quantified. Objective: We employ a novel mathematical model of norovirus
transmission, and fit the model to daily incidence data from a major norovirus
outbreak on a cruise ship, and examine the relative efficacy of potential
control strategies aimed at reducing environmental and/or direct transmission.
Results: The reproduction number for environmental and direct transmission
combined is Rtot = 11.0 [9.4,15.6], and of environmental transmission alone is
Renviron = 0.85 [0.18,2.04]. Direct transmission is overwhelmingly due to 0
passenger-to-passenger contacts, but crew can act as a reservoir of infection
from cruise-to-cruise. Implications: This is the first quantification of the
relative roles of environmental and direct transmission of norovirus. While
environmental transmission has the potential to maintain a sustained series of
outbreaks aboard a cruise ship in the absence of strict sanitation practices,
direct transmission dominates. Quarantine of ill passengers and cleaning are
likely to have little impact on final outbreak size, but intensive promotion of
good hand washing practices can prevent outbreaks.


  We introduce an Interaction and Trade-off based Eco-Evolutionary Model
(ITEEM), in which species are competing for common resources in a well-mixed
system, and their evolution in interaction trait space is subject to a
life-history trade-off between replication rate and competitive ability. We
demonstrate that the strength of the trade-off has a fundamental impact on
eco-evolutionary dynamics, as it imposes four phases of diversity, including a
sharp phase transition. Despite its minimalism, ITEEM produces without further
ad hoc features a remarkable range of observed patterns of eco-evolutionary
dynamics. Most notably we find self-organization towards structured communities
with high and sustainable diversity, in which competing species form
interaction cycles similar to rock-paper-scissors games.


  Whole genome amplification (WGA) plays an important role in sample
preparation of low-input templates for high-throughput sequencing. Multiple
displacement amplification (MDA), a popular isothermal WGA method, suffers a
major hurdle of highly uneven amplification. Optimizations have been made in
the past by separating the reagents into numbers of tiny chambers or droplets
in microfluidic devices, which significantly improves the amplification
uniformity of MDA. However, skill barrier still exists for biological
researchers to handle chip fabrication and droplet manipulation. Here, we
present a novel MDA protocol, in-capillary MDA (icMDA), which significantly
simplifies the manipulation and improves the uniformity of amplification by
dispersing reagents in a long quasi-1D capillary tubing. We demonstrated that
icMDA is able to accurately detect SNVs with higher efficiency and sensitivity.
Moreover, this straightforward method employs neither customized instruments
nor complicated operations, making it a ready-to-use approach for most
laboratories.


  Cytoskeletal motor proteins move on filamentous tracks by converting input
chemical energy that they derive by catalyzing the hydrolysis of ATP. The
ATPase site is the analog of an engine and hydrolysis of ATP is the analog of
burning of chemical fuel. Moreover, the functional role of a segment of the
motor is analogous to that of the transmission system of an automobile that
consists of shaft, gear, clutch, etc. The operation of the engine is
intrinsically "noisy" and the motor faces a molecular "hailstorm" in the
aqueous medium. In this commemorative article, we celebrate the centenary of
Michaelis and Menten's landmark paper of 1913 and the golden jubilee of Monod
et al.'s classic paper of 1963 by highlighting their relevance in explaining
the operational mechanisms of the engine and the transmission system,
respectively, of cytoskeletal motors.


  Single molecule experiments on bacteriophages show an exponential scaling for
the dependence of mobility on the length of DNA within the capsid. It has been
suggested that this could be due to the "capstan mechanism" -- the exponential
amplification of friction forces that result when a rope is wound around a
cylinder as in a ship's capstan. Here we describe a desktop experiment that
illustrates the effect. Though our model phage is a million times larger, it
exhibits the same scaling observed in single molecule experiments.


  In this paper we consider a delayed nonlinear model of the dynamics of the
immune system against a viral infection that contains wild-type virus and one
mutant. A finite response time of the immune system was considered in order
which leads to sustained oscillatory behavior as well as chaotic behavior,
triggered by the presence of delays. We present a thoroughly numeric analysis
and some analytical results.


  We investigate the dynamics of the heterodimer autorepression loop (HAL), a
small genetic module in which a protein A acts as an auto-repressor and binds
to a second protein B to form a AB dimer. For suitable values of the rate
constants the HAL produces pulses of A alternating with pulses of B. By means
of analytical and numerical calculations, we show that the duration of A-pulses
is extremely robust against variation of the rate constants while the duration
of the B-pulses can be flexibly adjusted. The HAL is thus a minimal genetic
module generating robust pulses with tunable duration an interesting property
for cellular signalling.


  Genomes contain rare guanine-rich sequences capable of assembling into
four-stranded helical structures, termed G-quadruplexes, with potential roles
in gene regulation and chromosome stability. Their mechanical unfolding has
only been reported to date by all-atom simulations, which cannot dissect the
major physical interactions responsible for their cohesion. Here, we propose a
mesoscopic model to describe both the mechanical and thermal stability of DNA
G-quadruplexes, where each nucleotide of the structure, as well as each central
cation located at the inner channel, is mapped onto a single bead. In this
framework we are able to simulate loading rates similar to the experimental
ones, which are not reachable in simulations with atomistic resolution. In this
regard, we present single-molecule force-induced unfolding experiments by a
high-resolution optical tweezers on a DNA telomeric sequence capable of forming
a G-quadruplex conformation. Fitting the parameters of the model to the
experiments we find a correct prediction of the rupture-force kinetics and a
good agreement with previous near equilibrium measurements. Since G-quadruplex
unfolding dynamics is halfway in complexity between secondary nucleic acids and
tertiary protein structures, our model entails a nanoscale paradigm for
non-equilibrium processes in the cell.


  It is useful to have complete lists of nonisomorphic chemical reaction
networks (CRNs) of a given size, with or without various restrictions. One may,
for example, be interested in exploring how often certain dynamical behaviours
occur in small CRNs, or wish to find examples to illustrate some aspect of the
theory. In such cases, it is natural to examine one representative from each
isomorphism class of CRNs. Inspired by the related project of Deckard et al
(Enumeration and online library of mass-action reaction networks.
http://arXiv.org/abs/0901.3067, 2009), this document outlines the methodology
involved in listing all CRNs in various classes of interest including, for
example, general CRNs, dynamically nontrivial CRNs, weakly reversible CRNs,
fully open CRNs, etc. The accompanying data (i.e., lists of nonisomorphic CRNs
in the various classes) is at https://reaction-networks.net/networks/. Note
that both document and data are work in progress.


  The incubation period of a disease is the time between an initiating
pathologic event and the onset of symptoms. For typhoid fever, polio, measles,
leukemia and many other diseases, the incubation period is highly variable.
Some affected people take much longer than average to show symptoms, leading to
a distribution of incubation periods that is right skewed and often
approximately lognormal. Although this statistical pattern was discovered more
than sixty years ago, it remains an open question to explain its ubiquity. Here
we propose an explanation based on evolutionary dynamics on graphs. For simple
models of a mutant or pathogen invading a network-structured population of
healthy cells, we show that skewed distributions of incubation periods emerge
for a wide range of assumptions about invader fitness, competition dynamics,
and network structure. The skewness stems from stochastic mechanisms associated
with two classic problems in probability theory: the coupon collector and the
random walk. Unlike previous explanations that rely crucially on heterogeneity,
our results hold even for homogeneous populations. Thus, we predict that two
equally healthy individuals subjected to equal doses of equally pathogenic
agents may, by chance alone, show remarkably different time courses of disease.


  Post-translational modification (PTM) of proteins plays a key role in signal
transduction, and hence significant effort has gone toward understanding how
PTM networks process information. This involves, on the theory side, analyzing
the dynamical systems arising from such networks. Which networks are, for
instance, bistable? Which networks admit sustained oscillations? Which
parameter values enable such behaviors? In this Perspective, we highlight
recent progress in this area and point out some important future directions.
Along the way, we summarize several techniques for analyzing general networks,
such as eliminating variables to obtain steady-state parametrizations, and
harnessing results on how incorporating intermediates affects dynamics.


  The advent of high--throughput transcription profiling technologies has
enabled identification of genes and pathways associated with disease, providing
new avenues for precision medicine. A key challenge is to analyze this data in
the context of the regulatory networks and pathways that control cellular
processes, while still obtaining insights that can be used to design new
diagnostic and therapeutic interventions. While classical differential
expression analysis provides specific and hence targetable gene-level insights,
it does not include any systems-level information. On the other hand, pathway
analyses integrate systems-level information with expression data, but are
often limited in their ability to indicate specific molecular targets. We
introduce GeneSurrounder, an analysis method that takes into account the
complex structure of interaction networks to identify specific genes that
disrupt pathway activity in a disease-specific manner. GeneSurrounder
integrates transcriptomic data and pathway network information in a novel
two-step procedure to detect genes that (i) appear to influence the expression
of other genes local to it in the network and (ii) are part of a subnetwork of
differentially expressed genes. Combined, this evidence can be used to pinpoint
specific genes that have a mechanistic role in the phenotype of interest.
Applying GeneSurrounder to three distinct ovarian cancer studies using a global
KEGG network, we show that our method is able to identify biologically relevant
genes and genes missed by single-gene association tests, integrate pathway and
expression data, and yield more consistent results across multiple studies of
the same phenotype than competing methods.


  The goal of this note is to present a general approach to define the net
reproduction function for a large class of nonlinear physiologically structured
population models. In particular, we are going to show that this can be
achieved in a natural way by reformulating a nonlinear problem as a family of
linear ones; each of the linear problems describing the evolution of the
population in a different, but constant environment. The reformulation of a
nonlinear population model as a family of linear ones is a new approach, and
provides an elegant way to study qualitative questions, for example the
existence of positive steady states. To define the net reproduction number for
any fixed (constant) environment, i.e. for the linear models, we use a fairly
recent spectral theoretic result, which characterizes the connection between
the spectral bound of an unbounded operator and the spectral radius of a
corresponding bounded operator. For nonlinear models, varying the environment
naturally leads to a net reproduction function.


  We present a simple and efficient method for prediction of transcription
factor binding sites from DNA sequence. Our method computes a random
approximation of a convolutional kernel feature map from DNA sequence and then
learns a linear model from the approximated feature map. Our method outperforms
state-of-the-art deep learning methods on five out of six test datasets from
the ENCODE consortium, while training in less than one eighth the time.


  Molecular variants of vitamin B12, siderophores and glycans occur. To take up
variant forms, bacteria may express an array of receptors. The gut microbe
Bacteroides thetaiotaomicron has three different receptors to take up variants
of vitamin B12 and 88 receptors to take up various glycans. The design of
receptor arrays reflects key processes that shape cellular evolution.
Competition may focus each species on a subset of the available nutrient
diversity. Some gut bacteria can take up only a narrow range of carbohydrates,
whereas species such as B.~thetaiotaomicron can digest many different complex
glycans. Comparison of different nutrients, habitats, and genomes provide
opportunity to test hypotheses about the breadth of receptor arrays. Another
important process concerns fluctuations in nutrient availability. Such
fluctuations enhance the value of cellular sensors, which gain information
about environmental availability and adjust receptor deployment. Bacteria often
adjust receptor expression in response to fluctuations of particular
carbohydrate food sources. Some species may adjust expression of uptake
receptors for specific siderophores. How do cells use sensor information to
control the response to fluctuations? That question about regulatory wiring
relates to problems that arise in control theory and artificial intelligence.
Control theory clarifies how to analyze environmental fluctuations in relation
to the design of sensors and response systems. Recent advances in deep learning
studies of artificial intelligence focus on the architecture of regulatory
wiring and the ways in which complex control networks represent and classify
environmental states. I emphasize the similar design problems that arise in
cellular evolution, control theory, and artificial intelligence. I connect
those broad concepts to testable hypotheses for bacterial uptake of B12,
siderophores and glycans.


  Numerical studies of the May-Leonard model for cyclically competing species
exhibit spontaneous spatial structures in the form of spirals. It is desirable
to obtain a simple coarse-grained evolution equation describing spatio-temporal
pattern formation in such spatially extended stochastic population dynamics
models. Extending earlier work on the corresponding deterministic system, we
derive the complex Ginzburg-Landau equation as the effective representation of
the fully stochastic dynamics of this paradigmatic model for cyclic dominance
near its Hopf bifurcation, and for small fluctuations in the three-species
coexistence regime. The internal stochastic reaction noise is accounted for
through the Doi-Peliti coherent-state path integral formalism, and subsequent
mapping to three coupled non-linear Langevin equations. This analysis provides
constraints on the model parameters that allow time scale separation and in
consequence a further reduction to just two coarse-grained slow degrees of
freedom.


  Phylogenetics has seen an steady increase in substitution model complexity,
which requires increasing amounts of computational power to compute
likelihoods. This model complexity motivates strategies to approximate the
likelihood functions for branch length optimization and Bayesian sampling. In
this paper, we develop an approximation to the one-dimensional likelihood
function as parametrized by a single branch length. This new method uses a
four-parameter surrogate function abstracted from the simplest phylogenetic
likelihood function, the binary symmetric model. We show that it offers a
surrogate that can be fit over a variety of branch lengths, that it is
applicable to a wide variety of models and trees, and that it can be used
effectively as a proposal mechanism for Bayesian sampling. The method is
implemented as a stand-alone open-source C library for calling from
phylogenetics algorithms; it has proven essential for good performance of our
online phylogenetic algorithm sts.


  Background: National or local laws, norms or regulations (sometimes and in
some countries) require medical providers to report notifiable diseases to
public health authorities. Reporting, however, is almost always incomplete.
This is due to a variety of reasons, ranging from not recognizing the diseased
to failures in the technical or administrative steps leading to the final
official register in the disease notification system. The reported fraction
varies from 9% to 99% and is strongly associated with the disease being
reported.
  Methods: In this paper we propose a method to approximately estimate the full
prevalence (and any other variable or parameter related to transmission
intensity) of infectious diseases. The model assumes incomplete notification of
incidence and allows the estimation of the non-notified number of infections
and it is illustrated by the case of hepatitis C in Brazil. The method has the
advantage that it can be corrected iteratively by comparing its findings with
empirical results.
  Results: The application of the model for the case of hepatitis C in Brazil
resulted in a prevalence of notified cases that varied between 163,902 and
169,382 cases; a prevalence of non-notified cases that varied between 1,433,638
and 1,446,771; and a total prevalence of infections that varied between
1,597,540 and 1,616,153 cases.
  Conclusions: We conclude that that the model proposed can be useful for
estimation of the actual magnitude of endemic states of infectious diseases,
particularly for those where the number of notified cases is only the tip of
the iceberg. In addition, the method can be applied to other situations, such
as the well known underreported incidence of criminality (for example rape),
among others.


  We discuss two threshold phenomena in blood coagulation dynamics using a
simplified model. This perspective of the underlying complex phenomena is
expected to aid the understanding and characterization of many blood
coagulation pathologies with altered protein dynamics.


  Motivation: Inferring the structure of gene regulatory networks from
high--throughput datasets remains an important and unsolved problem. Current
methods are hampered by problems such as noise, low sample size, and incomplete
characterizations of regulatory dynamics, leading to networks with missing and
anomalous links. Integration of prior network information (e.g., from pathway
databases) has the potential to improve reconstructions.
  Results: We developed a semi--supervised network reconstruction algorithm
that enables the synthesis of information from partially known networks with
time course gene expression data. We adapted PLS-VIP for time course data and
used reference networks to simulate expression data from which null
distributions of VIP scores are generated and used to estimate edge
probabilities for input expression data. By using simulated dynamics to
generate reference distributions, this approach incorporates previously known
regulatory relationships and links the network to the dynamics to form a
semi-supervised approach that discovers novel and anomalous connections. We
applied this approach to data from a sleep deprivation study with KEGG pathways
treated as prior networks, as well as to synthetic data from several DREAM
challenges, and find that it is able to recover many of the true edges and
identify errors in these networks, suggesting its ability to derive posterior
networks that accurately reflect gene expression dynamics.


  Motivation
  Protein fold recognition is an important problem in structural
bioinformatics. Almost all traditional fold recognition methods use sequence
(homology) comparison to indirectly predict the fold of a tar get protein based
on the fold of a template protein with known structure, which cannot explain
the relationship between sequence and fold. Only a few methods had been
developed to classify protein sequences into a small number of folds due to
methodological limitations, which are not generally useful in practice.
  Results
  We develop a deep 1D-convolution neural network (DeepSF) to directly classify
any protein se quence into one of 1195 known folds, which is useful for both
fold recognition and the study of se quence-structure relationship. Different
from traditional sequence alignment (comparison) based methods, our method
automatically extracts fold-related features from a protein sequence of any
length and map it to the fold space. We train and test our method on the
datasets curated from SCOP1.75, yielding a classification accuracy of 80.4%. On
the independent testing dataset curated from SCOP2.06, the classification
accuracy is 77.0%. We compare our method with a top profile profile alignment
method - HHSearch on hard template-based and template-free modeling targets of
CASP9-12 in terms of fold recognition accuracy. The accuracy of our method is
14.5%-29.1% higher than HHSearch on template-free modeling targets and
4.5%-16.7% higher on hard template-based modeling targets for top 1, 5, and 10
predicted folds. The hidden features extracted from sequence by our method is
robust against sequence mutation, insertion, deletion and truncation, and can
be used for other protein pattern recognition problems such as protein
clustering, comparison and ranking.


  We introduce a new random graph model motivated by biological questions
relating to speciation. This random graph is defined as the stationary
distribution of a Markov chain on the space of graphs on $\{1, \ldots, n\}$.
The dynamics of this Markov chain is governed by two types of events: vertex
duplication, where at constant rate a pair of vertices is sampled uniformly and
one of these vertices loses its incident edges and is rewired to the other
vertex and its neighbors; and edge removal, where each edge disappears at
constant rate. Besides the number of vertices $n$, the model has a single
parameter $r_n$.
  Using a coalescent approach, we obtain explicit formulas for the first
moments of several graph invariants such as the number of edges or the number
of complete subgraphs of order $k$. These are then used to identify five
non-trivial regimes depending on the asymptotics of the parameter $r_n$. We
derive an explicit expression for the degree distribution, and show that under
appropriate rescaling it converges to classical distributions when the number
of vertices goes to infinity. Finally, we give asymptotic bounds for the number
of connected components, and show that in the sparse regime the number of edges
is Poissonian.


  Until recently, much of the microbial world was hidden from view. A global
research effort has changed this, unveiling and quantifying microbial diversity
across enormous range of critically-important contexts, from the human
microbiome, to plant-soil interactions, to marine life. Yet what has remained
largely hidden is the interplay of ecological and evolutionary processes that
led to the diversity we observe in the present day. We introduce a theoretical
framework to quantify the effect of ecological innovations in microbial
evolutionary history, using a new, coarse-grained approach that is robust to
the incompleteness and ambiguities in microbial community data. Applying this
methodology, we identify a balance of gradual, ongoing diversification and
rapid bursts across a vast range of microbial habitats. Moreover, we find
universal quantitative similarities in the tempo of diversification,
independent of habitat type.


  In this study, we construct an evolutionary model of a population of DNA
sequences interacting with the surrounding environment on the topological
monoid A* of strings on the alphabet A = { a, c, g, t }. A partial differential
equation governing the evolution of the DNA population is derived as a kind of
diffusion equation on A*. Analyzing the constructed model in a theoretical
manner, we present conditions for sympatric speciation, the possibility of
which continues to be discussed. It is shown that under other same conditions
one condition determines whether sympatric speciation occurs or the DNA
population continues to move around randomly in a subset of A*. We next
demonstrate that the population maintains a kind of equlibrium state under
certain conditions. In this situation, the population remains nearly unchanged
and does not differentiate even if it can differentiate into others.
Furthermore, we calculate the probability of sympatric speciation and the time
expected to elapse before it.


  The use of differential equations (ODE) is one of the most promising
approaches to network inference. The success of ODE-based approaches has,
however, been limited, due to the difficulty in estimating parameters and by
their lack of scalability. Here we introduce a novel method and pipeline to
reverse engineer gene regulatory networks from gene expression of time series
and perturbation data based upon an improvement on the calculation scheme of
the derivatives and a pre-filtration step to reduce the number of possible
links. The method introduces a linear differential equation model with adaptive
numerical differentiation that is scalable to extremely large regulatory
networks. We demonstrate the ability of this method to outperform current
state-of-the-art methods applied to experimental and synthetic data using test
data from the DREAM4 and DREAM5 challenges. Our method displays greater
accuracy and scalability. We benchmark the performance of the pipeline with
respect to data set size and levels of noise. We show that the computation time
is linear over various network sizes.


  In this paper we work out a parameterization of the environment noise within
the Haken-Strobl-Reinenker (HSR) model for the PE545 light-harvesting complex
based on atomic-level quantum mechanics/molecular mechanics (QM/MM)
simulations. We use this approach to investigate the role of different auto-
and cross-correlations in the HSR noise tensor, confirming that site-energy
autocorrelations (pure dephasing) terms dominate the noise-induced exciton
mobility enhancement, followed by site energy-coupling cross-correlations for
specific triplets of pigments. Interestingly, several cross-correlations of the
latter kind, together with coupling-coupling cross-correlations, display clear
low-frequency signatures in their spectral densities in the region 30-70
inverse centimeters. These slow components lie at the limits of validity of the
HSR approach, requiring that environmental fluctuations be faster than typical
exciton transfer time scales. We show that a simple coarse-grained
elastic-network-model (ENM) analysis of the PE545 protein naturally spotlights
collective normal modes in this frequency range, that represent specific
concerted motions of the subnetwork of cysteines that are covalenty linked to
the pigments. This analysis strongly suggests that protein scaffolds in
light-harvesting complexes are able to express specific collective,
low-frequency normal modes providing a fold-rooted blueprint of exciton
transport pathways. We speculate that ENM-based mixed quantum classical
methods, such as Ehrenfest dynamics, might be promising tools to disentangle
the fundamental designing principles of these dynamical processes in natural
and artificial light-harvesting structures.


  A protein is traditionally visualised as a piecewise linear discrete curve,
and its geometry is conventionally characterised by the extrinsically
determined Ramachandran angles. However, a protein backbone has also two
independent intrinsic geometric structures, due to the peptide planes and the
side chains. Here we adapt and develop modern 3D virtual reality techniques to
scrutinize the atomic geometry along a protein backbone, in the vicinity of a
peptide plane. For this we compare backbone geometry-based (extrinsic) and
structure-based (intrinsic) coordinate systems, and as an example we inspect
the trans and cis peptide planes. We reveal systematics in the way how a cis
peptide plane deforms the neighbouring atomic geometry, and we develop a
virtual reality based visual methodology that can identify the presence of a
cis peptide plane from the arrangement of atoms in its vicinity. Our approach
can easily detect exceptionally placed atoms in crystallographic structures.
Thus it can be employed as a powerful visual refinement tool which is
applicable also in the case when resolution of the protein structure is limited
and whenever refinement is needed. As concrete examples we identify a number of
crystallographic protein structures in Protein Data Bank (PDB) that display
exceptional atomic positions around their cis peptide planes.


  The Indo-Pacific region contains a unique mix of opportunities for the
development and use of genetic-pest-management, gene-drive, and gene-drive-like
technologies. Here I collectively refer to these technologies as Evolutionary
Genetic Engineering (EGE). Indo-Pacific Islands have some of the world's
highest rates of endemism and extinction---species and entire ecosystems are at
risk. This threat to the natural world is coupled with the burden of human
diseases, many of which are new and emerging or neglected tropical diseases.
The same factors which have led to high rates of endemism also, in some ways,
make this region an ideal testing ground for some types of EGE's. There is
great potential for positive humanitarian, economic, and conservation
applications of EGE's. However, these types of new technologies will be
initially viewed from the perspective of the recent history of a loss of self
determination, issues of social justice, and the testing of new technologies
(e.g., biocontrol, agricultural, nuclear) in the Indo-Pacific---a region of the
world that is still extensively colonized and controlled by Western Nations.
Experience with successes and failures in related technologies suggests a path
to move forward---a set of eight recommendations---to maximize the potential
payoffs and minimize unintended negative effects of EGE's.


  A system-level framework of complex microbe-microbe and host-microbe chemical
cross-talk would help elucidate the role of our gut microbiota in health and
disease. Here we report a literature-curated interspecies network of the human
gut microbiota, called NJS16. This is an extensive data resource composed of
~570 microbial species and 3 human cell types metabolically interacting through
>4,400 small-molecule transport and macromolecule degradation events. Based on
the contents of our network, we develop a mathematical approach to elucidate
representative microbial and metabolic features of the gut microbial community
in a given population, such as a disease cohort. Applying this strategy to
microbiome data from type 2 diabetes patients reveals a context-specific
infrastructure of the gut microbial ecosystem, core microbial entities with
large metabolic influence, and frequently-produced metabolic compounds that
might indicate relevant community metabolic processes. Our network presents a
foundation towards integrative investigations of community-scale microbial
activities within the human gut.


  In this paper, the stability analysis of an anaerobic digestion process is
presented. The analysis is performed using a simplified mathematical model,
which includes explicit temperature and pH dependence on kinetic growth rates.
A detailed nonlinear and bifurcation analyses are performed in order to study
the effects of temperature and pH parameters on the process behaviour. In
addition, both safety and optimal operation regions of the bioreactor are
established. Based on bifurcation diagrams, it is observed that the washout
condition occurs by combining a fold bifurcation and a transcritical
bifurcation. Consequently, to prevent the washout and guarantee optimal
operation conditions of the bioreactor, a risk criterion oriented to monitor
the bioprocess on-line is proposed, allowing to detect the system
destabilisation.


  Circadian clocks are biochemical oscillators that allow organisms to estimate
the time of the day. These oscillators are inherently noisy due to the discrete
nature of the reactants and the stochastic character of their interactions. To
keep these oscillators in sync with the daily day-night rhythm in the presence
of noise, circadian clocks must be coupled to the dark-light cycle. In this
manuscript, we study the entrainment of phase oscillators as a function of the
intrinsic noise in the system. Using stochastic simulations, we compute the
optimal coupling strength, intrinsic frequency and shape of the phase-response
curve, that maximize the mutual information between the phase of the clock and
time. We show that the optimal coupling strength and intrinsic frequency
increase with the noise, but that the shape of the phase-response curve varies
non-monotonically with the noise: in the low-noise regime, it features a
deadzone that increases in width as the noise increases, while in the
high-noise regime, the width decreases with the noise. These results arise from
a trade-off between maximizing stability -noise suppression - and maximizing
linearity of the input-output, i.e. time- phase, relation. We also show that
three analytic approximations -the linear-noise approximation, the
phase-averaging method, and linear-response theory- accurately describe
different regimes of the coupling strength and the noise.


  Epidemiological models for the spread of pathogens in a population are
usually only able to describe a single pathogen. This makes their application
unrealistic in cases where multiple pathogens with similar symptoms are
spreading concurrently within the same population. Here we describe a method
which makes possible the application of multiple single-strain models under
minimal conditions. As such, our method provides a bridge between theoretical
models of epidemiology and data-driven approaches for modeling of influenza and
other similar viruses.
  Our model extends the Susceptible-Infected-Recovered model to higher
dimensions, allowing the modeling of a population infected by multiple viruses.
We further provide a method, based on an overcomplete dictionary of feasible
realizations of SIR solutions, to blindly partition the time series
representing the number of infected people in a population into individual
components, each representing the effect of a single pathogen.
  We demonstrate the applicability of our proposed method on five years of
seasonal influenza-like illness (ILI) rates, estimated from Twitter data. We
demonstrate that our method describes, on average, 44\% of the variance in the
ILI time series. The individual infectious components derived from our model
are matched to known viral profiles in the populations, which we demonstrate
matches that of independently collected epidemiological data. We further show
that the basic reproductive numbers ($R0$) of the matched components are in
range known for these pathogens.
  Our results suggest that the proposed method can be applied to other
pathogens and geographies, providing a simple method for estimating the
parameters of epidemics in a population.


  Determining the principal energy pathways for allosteric communication in
biomolecules, that occur as a result of thermal motion, remains challenging due
to the intrinsic complexity of the systems involved. Graph theory provides an
approach for making sense of such complexity, where allosteric proteins can be
represented as networks of amino acids. In this work, we establish the
eigenvector centrality metric in terms of the mutual information, as a mean of
elucidating the allosteric mechanism that regulates the enzymatic activity of
proteins. Moreover, we propose a strategy to characterize the range of the
physical interactions that underlie the allosteric process. In particular, the
well known enzyme, imidazol glycerol phosphate synthase (IGPS), is utilized to
test the proposed methodology. The eigenvector centrality measurement
successfully describes the allosteric pathways of IGPS, and allows to pinpoint
key amino acids in terms of their relevance in the momentum transfer process.
The resulting insight can be utilized for refining the control of IGPS
activity, widening the scope for its engineering. Furthermore, we propose a new
centrality metric quantifying the relevance of the surroundings of each
residue. In addition, the proposed technique is validated against experimental
solution NMR measurements yielding fully consistent results. Overall, the
methodologies proposed in the present work constitute a powerful and cost
effective strategy to gain insight on the allosteric mechanism of proteins.


  We analyse the unfolding pathway of biomolecules comprising several
independent modules in pulling experiments. In a recently proposed model, a
critical velocity $v_{c}$ has been predicted, such that for pulling speeds
$v>v_{c}$ it is the module at the pulled end that opens first, whereas for
$v<v_{c}$ it is the weakest. Here, we introduce a variant of the model that is
closer to the experimental setup, and discuss the robustness of the emergence
of the critical velocity and of its dependence on the model parameters. We also
propose a possible experiment to test the theoretical predictions of the model,
which seems feasible with state-of-art molecular engineering techniques.


  Influenza remains a significant burden on health systems. Effective responses
rely on the timely understanding of the magnitude and the evolution of an
outbreak. For monitoring purposes, data on severe cases of influenza in England
are reported weekly to Public Health England. These data are both readily
available and have the potential to provide valuable information to estimate
and predict the key transmission features of seasonal and pandemic influenza.
We propose an epidemic model that links the underlying unobserved influenza
transmission process to data on severe influenza cases. Within a Bayesian
framework, we infer retrospectively the parameters of the epidemic model for
each seasonal outbreak from 2012 to 2015, including: the effective reproduction
number; the initial susceptibility; the probability of admission to intensive
care given infection; and the effect of school closure on transmission. The
model is also implemented in real time to assess whether early forecasting of
the number of admission to intensive care is possible. Our model of admissions
data allows reconstruction of the underlying transmission dynamics revealing:
increased transmission during the season 2013/14 and a noticeable effect of
Christmas school holiday on disease spread during season 2012/13 and 2014/15.
When information on the initial immunity of the population is available,
forecasts of the number of admissions to intensive care can be substantially
improved. Readily available severe case data can be effectively used to
estimate epidemiological characteristics and to predict the evolution of an
epidemic, crucially allowing real-time monitoring of the transmission and
severity of the outbreak.


  We investigate spatially inhomogeneous versions of the stochastic
Lotka-Volterra model for predator-prey competition and coexistence by means of
Monte Carlo simulations on a two-dimensional lattice with periodic boundary
conditions. To study boundary effects for this paradigmatic population dynamics
system, we employ a simulation domain split into two patches: Upon setting the
predation rates at two distinct values, one half of the system resides in an
absorbing state where only the prey survives, while the other half attains a
stable coexistence state wherein both species remain active. At the domain
boundary, we observe a marked enhancement of the predator population density.
The predator correlation length displays a minimum at the boundary, before
reaching its asymptotic constant value deep in the active region. The frequency
of the population oscillations appears only very weakly affected by the
existence of two distinct domains, in contrast to their attenuation rate, which
assumes its largest value there. We also observe that boundary effects become
less prominent as the system is successively divided into subdomains in a
checkerboard pattern, with two different reaction rates assigned to neighboring
patches. When the domain size becomes reduced to the scale of the correlation
length, the mean population densities attain values that are very similar to
those in a disordered system with randomly assigned reaction rates drawn from a
bimodal distribution.


  In recent years the role of epidemic models in informing public health
policies has progressively grown. Models have become increasingly realistic and
more complex, requiring the use of multiple data sources to estimate all
quantities of interest. This review summarises the different types of
stochastic epidemic models that use evidence synthesis and highlights current
challenges.


  The "Kill the Winner" hypothesis is an attempt to address the problem of
diversity in biology. It argues that host-specific predators control the
population of each prey, preventing a winner from emerging and thus maintaining
the coexistence of all species in the system. We develop a stochastic model for
the "Kill the Winner" paradigm and show that the stable coexistence state of
the deterministic "Kill the Winner" model is destroyed by demographic
stochasticity, through a cascade of extinction events. We formulate an
individual-level stochastic model in which predator-prey coevolution promotes
high diversity of the ecosystem by generating a persistent population flux of
species.


  Proliferating cells properly divide into their daughter cells through a
process that is mediated by kinetochores, protein-complexes that assemble at
the centromere of each sister chromatid. Each kinetochore has to establish a
tight bipolar attachment to the spindle apparatus before sister-chromatid
separation is initiated. The Spindle Assembly Checkpoint (SAC) links the
biophysical attachment status of the kinetochores to mitotic progression, and
ensures that even a single misaligned kinetochore keeps the checkpoint active.
The mechanism by which this is achieved is still elusive. Current computational
models of the human SAC disregard important biochemical properties by omitting
any kind of feedback loop, proper kinetochore signals, and other spatial
properties such as the stability of the system and diffusion effects. To allow
for more realistic in silico study of the dynamics of the SAC model, a minimal
mathematical framework for SAC activation and silencing is introduced. A
nonlinear ordinary differential equation model successfully reproduces
bifurcation signaling switches with attachment of all 92 kinetochores and
activation of APC/C by kinetochore-driven feedback. A partial differential
equation model and mathematical linear stability analyses indicate the
influence of diffusion and system stability. The conclusion is that
quantitative models of the human SAC should account for the positive feedback
on APC/C activation driven by the kinetochores which is essential for SAC
silencing. Experimental diffusion coefficients for MCC sub-complexes are found
to be insufficient for rapid APC/C inhibition. The presented analysis allows
for systems-level understanding of mitotic control and the minimal new model
can function as a basis for developing further quantitative-integrative models
of the cell division cycle


  The managed-metabolism hypothesis suggests that a cooperation barrier must be
overcome if self-producing chemical organizations are to transition from
non-life to life. This barrier prevents un-managed, self-organizing,
autocatalytic networks of molecular species from individuating into complex,
cooperative organizations. The barrier arises because molecular species that
could otherwise make significant cooperative contributions to the success of an
organization will often not be supported within the organization, and because
side reactions and other free-riding processes will undermine cooperation. As a
result, the barrier seriously limits the possibility space that can be explored
by un-managed organizations, impeding individuation, complex functionality and
the transition to life. The barrier can be overcome comprehensively by
appropriate management which implements a system of evolvable constraints. The
constraints support beneficial co-operators and suppress free riders. In this
way management can manipulate the chemical processes of an autocatalytic
organization, producing novel processes that serve the interests of the
organization as a whole and that could not arise and persist spontaneously in
an un-managed chemical organization. Management self-organizes because it is
able to capture some of the benefits that are produced when its management of
an autocatalytic organization promotes beneficial cooperation. Selection
therefore favours the emergence of managers that take over and manage chemical
organizations so as to overcome the cooperation barrier. The managed-metabolism
hypothesis shows that if management is to overcome the cooperation barrier
comprehensively, its interventions must be digitally coded. In this way, the
hypothesis accounts for the two-tiered structure of all living cells in which a
digitally-coded genetic apparatus manages an analogically-informed metabolism.


  Cooperation is ubiquitous in every level of living organisms. It is known
that spatial (network) structure is a viable mechanism for cooperation to
evolve. Until recently, it has been difficult to predict whether cooperation
can evolve at a network (population) level. To address this problem, Pinheiro
et al. proposed a numerical metric, called Average Gradient of Selection (AGoS)
in 2012. AGoS can characterize and forecast the evolutionary fate of
cooperation at a population level. However, stochastic mutation of strategies
was not considered in the analysis of AGoS. Here we analyzed the evolution of
cooperation using AGoS where mutation may occur to strategies of individuals in
networks. Our analyses revealed that mutation always has a negative effect on
the evolution of cooperation regardless of the fraction of cooperators and
network structures. Moreover, we found that mutation affects the fitness of
cooperation differently on different social network structures.


  Development of new medications is a very lengthy and costly process. Finding
novel indications for existing drugs, or drug repositioning, can serve as a
useful strategy to shorten the development cycle. In this study, we present an
approach to drug discovery or repositioning by predicting indication for a
particular disease based on expression profiles of drugs, with a focus on
applications in psychiatry. Drugs that are not originally indicated for the
disease but with high predicted probabilities serve as good candidates for
repurposing. This framework is widely applicable to any chemicals or drugs with
expression profiles measured, even if the drug targets are unknown. It is also
highly flexible as virtually any supervised learning algorithms can be used. We
applied this approach to identify repositioning opportunities for schizophrenia
as well as depression and anxiety disorders. We applied various
state-of-the-art machine learning (ML) approaches for prediction, including
deep neural networks, support vector machines (SVM), elastic net, random forest
and gradient boosted machines. The performance of the five approaches did not
differ substantially, with SVM slightly outperformed the others. However,
methods with lower predictive accuracy can still reveal literature-supported
candidates that are of different mechanisms of actions. As a further
validation, we showed that the repositioning hits are enriched for psychiatric
medications considered in clinical trials. Notably, many top repositioning hits
are supported by previous preclinical or clinical studies. Finally, we propose
that ML approaches may provide a new avenue to explore drug mechanisms via
examining the variable importance of gene features.


  How cooperation can evolve between players is an unsolved problem of biology.
Here we use Hamiltonian dynamics of models of the Ising type to describe
populations of cooperating and defecting players to show that the equilibrium
fraction of cooperators is given by the expectation value of a thermal
observable akin to a magnetization. We apply the formalism to the Public Goods
game with three players, and show that a phase transition between cooperation
and defection occurs that is equivalent to a transition in one-dimensional
Ising crystals with long-range interactions. We then investigate the effect of
punishment on cooperation and find that punishment plays the role of a magnetic
field that leads to an "alignment" between players, thus encouraging
cooperation. We suggest that a thermal Hamiltonian picture of the evolution of
cooperation can generate other insights about the dynamics of evolving groups
by mining the rich literature of critical dynamics in low-dimensional spin
systems.


  Phylogenetic networks are becoming increasingly popular in phylogenetics
since they have the ability to describe a wider range of evolutionary events
than their tree counterparts. In this paper, we study Markov models on
phylogenetic networks and their associated geometry. We restrict our attention
to large-cycle networks, networks with a single undirected cycle of length at
least four. Using tools from computational algebraic geometry, we show that the
semi-directed network topology is generically identifiable for Jukes-Cantor
large-cycle network models.


  Aging in Caenorhabditis elegans is controlled, in part, by the insulin-like
signaling and heat shock response pathways. Following thermal stress,
expression levels of small heat shock protein 16.2 show a spatial patterning
across the 20 intestinal cells that reside along the length of the worm. Here,
we present a hypothesized mechanism that could lead to this patterned response
and develop a mathematical model of this system to test our hypothesis. We
propose that the patterned expression of heat shock protein is caused by a
diffusion-driven instability within the pseudocoelom, or fluid-filled cavity,
that borders the intestinal cells in C. elegans. This instability is due to the
interactions between two classes of insulin like peptides that serve
antagonistic roles. We examine output from the developed model and compare it
to experimental data on heat shock protein expression. Furthermore, we use the
model to gain insight on possible biological parameters in the system. The
model presented is capable of producing patterns similar to what is observed
experimentally and provides a first step in mathematically modeling
aging-related mechanisms in C. elegans.


  Background: There is limited information about the role of different age
groups, particularly subgroups of school-age children and younger adults in
propagating influenza epidemics.
  Methods: For a communicable disease outbreak, some subpopulations may play a
disproportionate role during the ascent of the outbreak due to increased
susceptibility and/or contact rates. Such subpopulations can be identified by
considering the proportion that cases in a subpopulation represent among all
cases in the population occurring before the epidemic peak (Bp), the
corresponding proportion after the epidemic peak (Ap), to calculate the
relative risk for a subpopulation, RR=Bp/Ap. We estimated RR for several age
groups using data on reported influenza A cases in Germany between 2002-2017.
  Results: Children aged 14-17y had the highest RR estimates for 7 out of 15
influenza A epidemics in the data, including the 2009 pandemic, and the large
2016/17, 2008/09, and 2006/07 seasons. Children aged 10-13y had the highest RR
estimates during 3 epidemics, including the large 2014/15 and 2004/05 seasons.
Children aged 6-9y had the highest RR estimates during two epidemics, including
the large 2012/13 season. Children aged 2-5y had the highest RR estimate during
the moderate 2015/16 season; adults aged 18-24y had the highest RR estimate
during the small 2005/06 season; adults aged 25-34y had the highest RR estimate
during the large, 2002/03 season.
  Conclusions: Our results support the prominent role of all school-age
children, particularly the oldest ones, in propagating influenza epidemics in
the community. We note that national vaccination coverage levels among older
school-age children were lower than among younger school-age children during
the recent influenza seasons in the US, and influenza vaccination program in
England has not been phased in yet for secondary school students.


  We present a computational parameter estimation tool aimed at assisting
molecule characterization applying nonlinear least-squares optimization to a
specific closed-form solution of a capillary electrophoretic transport problem.
Numerical estimates for the respective physicochemical parameters are
generated. Among other features, the tool facilitates automated and
user-controllable hypothesis testing and model selection. It can also be used
to perform fully automated batch-processing in true concurrent mode.


  Plants depend on the signaling of the phytohormone auxin for their
development and for responding to environmental perturbations. The associated
biomolecular signaling network involves a negative feedback at the level of the
Aux/IAA proteins which mediate the influence of auxin (the signal) on the ARF
transcription factors (the drivers of the response). To probe the role of this
feedback, we consider alternative in silico signaling networks implementing
different operating principles. By a comparative analysis, we find that the
presence of a negative regulatory feedback loop allows the system to have a far
larger sensitivity in its dynamical response to auxin. At the same time, this
sensitivity does not prevent the system from being highly resilient. Given this
insight, we reconsider previously published models and build a new quantitative
and calibrated biomolecular model of auxin signaling.


  In species reproducing both sexually and asexually clones are often more
common in recently established populations. Earlier studies have suggested that
this pattern arises from natural selection favouring asexual recruitment in
young populations. Alternatively, as we show here, this pattern may result from
stochastic processes during species-range expansions. We model a dioecious
species expanding into a new area in which all individuals are capable of both
sexual and asexual reproduction, and all individuals have equal survival rates
and dispersal distances. Even under conditions that eventually favour sexual
recruitment, colonisation starts with an asexual wave. Long after colonisation
is completed, a sexual wave erodes clonal dominance. If individuals reproduce
more than one season, and with only local dispersal, a few large clones
typically dominate for thousands of reproductive seasons. Adding occasional
long-distance dispersal, more dominant clones emerge, but they persist for a
shorter period of time. The general mechanism involved is simple: edge effects
at the expansion front favour asexual (uniparental) recruitment where potential
mates are rare. Specifically, our stochastic model makes detailed predictions
different from a selection model, and comparing these with empirical data from
a postglacially established seaweed species (Fucus radicans) shows that in this
case a stochastic mechanism is strongly supported.


  With the tremendous increase in the amount of biological literature,
developing automated methods for extracting big data from papers, building
models and explaining big mechanisms becomes a necessity. We describe here our
approach to translating machine reading outputs, obtained by reading bio-
logical signaling literature, to discrete models of cellular networks. We use
out- puts from three different reading engines, and describe our approach to
translating their different features, using examples from reading cancer
literature. We also outline several issues that still arise when assembling
cellular network models from state-of-the-art reading engines. Finally, we
illustrate the details of our approach with a case study in pancreatic cancer.


  The quasispecies model introduced by Eigen in 1971 has close connections with
the isometry group of the space of binary sequences relative to the Hamming
distance metric. Generalizing this observation we introduce an abstract
quasispecies model on a finite metric space $X$ together with a group of
isometries $\Gamma$ acting transitively on $X$. We show that if the domain of
the fitness function has a natural decomposition into the union of $t$
$G$-orbits, $G$ being a subgroup of $\Gamma$, then the dominant eigenvalue of
the evolutionary matrix satisfies an algebraic equation of degree at most
$t\cdot {\rm rk}_{\mathbf Z} R$, where $R$ is what we call the orbital ring.
The general theory is illustrated by two examples, in both of which $X$ is
taken to be the metric space of vertices of a regular polytope with the "edge"
metric; namely, the case of a regular $m$-gon and of a hyperoctahedron are
considered.


  RNA viruses comprise vast populations of closely related, but highly
genetically diverse, entities known as quasispecies. Understanding the
mechanisms by which this extreme diversity is generated and maintained is
fundamental when approaching viral persistence and pathobiology in infected
hosts. In this paper, we access quasispecies theory through a mathematical
model based on the theory of multitype branching processes, to better
understand the roles of mechanisms resulting in viral diversity, persistence
and extinction. We accomplish this understanding by a combination of
computational simulations and the theoretical analysis of the model. In order
to perform the simulations, we have implemented the mathematical model into a
computational platform capable of running simulations and presenting the
results in a graphical format in real time. Among other things, we show that
the establishment of virus populations may display four distinct regimes from
its introduction into new hosts until achieving equilibrium or undergoing
extinction. Also, we were able to simulate different fitness distributions
representing distinct environments within a host which could either be
favorable or hostile to the viral success. We addressed the most used
mechanisms for explaining the extinction of RNA virus populations called lethal
mutagenesis and mutational meltdown. We were able to demonstrate a
correspondence between these two mechanisms implying the existence of a
unifying principle leading to the extinction of RNA viruses.


  The volume occupied by the unconstrained genomic DNA of prokaryotes in saline
solutions is thousand times larger than the cell. Moreover, it is not separated
from the rest of the cell by a membrane. Nevertheless, it occupies only a small
fraction of the cell called the nucleoid. The mechanisms leading to such
compaction are the matter of ongoing debates. The present work aims at
exploring a newly proposed mechanism, according to which the formation of the
nucleoid would result from the demixing of the DNA and non-binding globular
macromolecules of the cytoplasm, like ribosomes. To this end, a coarse-grained
model of prokaryotic cells was developed and demixing was analyzed as a
function of the size and number of crowders. The model suggests that compaction
of the DNA is actually governed by the volume occupancy ratio of the crowders
and remains weak almost up to the jamming critical density. Strong compaction
is however observed just before jamming, suggesting that crowding and
electrostatic repulsion work synergetically in this limit. Finally, simulations
performed with crowders with different sizes indicate that the DNA and the
largest crowders demix preferentially. Together with the recent observation of
the gradual compaction of long DNA molecules upon increase of the concentration
of BSA proteins and silica nanoparticles, this work supports the demixing
mechanism as a key player for the formation of the nucleoid.


  During the alphavirus lifecycle, a nucleocapsid core buds through the cell
membrane to acquire an outer envelope of lipid membrane and viral
glycoproteins. However, the presence of a nucleocapsid core is not required for
assembly of infectious particles. To determine the role of the nucleocapsid
core, we develop a coarse-grained computational model with which we investigate
budding dynamics as a function of glycoprotein and nucleocapsid interactions,
as well as budding in the absence of a nucleocapsid. We find that there is a
transition between glycoprotein-directed budding and nucleocapsid-directed
budding which occurs above a threshold strength of nucleocapsid interactions.
The simulations predict that glycoprotein-directed budding leads to
significantly increased size polydispersity and particle polymorphism. This
polydispersity can be qualitatively explained by a theoretical model accounting
for the competition between bending energy of the membrane and the glycoprotein
shell. The simulations also show that the geometry of a budding particle leads
to a barrier to subunit diffusion, which can result in a stalled, partially
budded state. We present a phase diagram for this and other morphologies of
budded particles. Comparison of these structures against experiments could
establish bounds on whether budding is directed by glycoprotein or nucleocapsid
interactions. Although our model is motivated by alphaviruses, we discuss
implications of our results for other enveloped viruses.


  Biomedical research results are being published at a high rate, and with
existing search engines, the vast amount of published work is usually easily
accessible. However, reproducing published results, either experimental data or
observations is often not viable. In this work, we propose a framework to
overcome some of the issues of reproducing previous research, and to ensure
re-usability of published information. We present here a framework that
utilizes the results from state-of-the-art biomedical literature mining,
biological system modeling and analysis techniques, and provides means to
scientists to assemble and reason about information from voluminous, fragmented
and sometimes inconsistent literature. The overall process of automated
reading, assembly and reasoning can speed up discoveries from the order of
decades to the order of hours or days. Our framework described here allows for
rapidly conducting thousands of in silico experiments that are designed as part
of this process.


  Most of zoonoses are multi-host parasites with multiple transmission routes
that are usually investigated separately despite their potential interplay. As
a unifying framework for modelling parasite spread through different paths of
infection, we suggest "ecomultiplex" networks, i.e. multiplex networks
representing interacting animal communities with (i) spatial structure and (ii)
metabolic scaling. We exploit this ecological framework for testing potential
control strategies for $T. cruzii$ spread in two real-world ecosystems. Our
investigation highlights two interesting results. Firstly, the ecomultiplex
topology can be as efficient as more data-demanding epidemiological measures in
identifying which species facilitate parasite spread. Secondly, the interplay
between predator-prey and host-parasite interactions leads to a phenomenon of
parasite amplification in which top predators facilitate $T. cruzii$ spread,
offering theoretical interpretation of previous empirical findings. Our
approach is broadly applicable and could provide novel insights in designing
immunisation strategies for pathogens with multiple transmission routes in
real-world ecosystems.


  In this paper we study a nonlinear age structured model of locust population
dynamics with variable time of egg incubation that describes the phase
polyphenism and behaviour of desert locust, Schistocerca gregaria. The analysis
of asymptotical stability of trivial and nontrivial equilibriums of autonomous
system allows us to derive the conditions and understand the particularities of
bidirectional phase transitions between solitarious and gregarious. Simulation
of locust population dynamics with different conditions of phase transitions
for autonomous and non-autonomous dynamical systems of our two phase
age-structured competitive model with time delay exhibits the features of
behaviour of documented outbreak dynamics of Schistocerca gregaria.


  In biodiversity conservation it is often necessary to prioritize the species
to conserve. Existing approaches to prioritization, e.g. the Fair Proportion
Index and the Shapley Value, are based on phylogenetic trees and rank species
according to their contribution to overall phylogenetic diversity. However, in
many cases evolution is not treelike and thus, phylogenetic networks have come
to the fore as a generalization of phylogenetic trees, allowing for the
representation of non-treelike evolutionary events, such as horizontal gene
transfer or hybridization. Here, we extend the concepts of phylogenetic
diversity and phylogenetic diversity indices from phylogenetic trees to
phylogenetic networks. On the one hand, we consider the treelike content of a
phylogenetic network, e.g. the (multi)set of phylogenetic trees displayed by a
network and the LSA tree associated with it. On the other hand, we derive the
phylogenetic diversity of subsets of taxa and biodiversity indices directly
from the internal structure of the network. Furthermore, we introduce our
software package NetDiversity, which was implemented in Perl and allows for the
calculation of all generalized measures of phylogenetic diversity and
generalized phylogenetic diversity indices established in this note that are
independent of inheritance probabilities. We apply our methods to a phylogentic
network representing the evolutionary relationships among swordtails and
platyfishes (Xiphophorus: Poeciliidae), a group of species characterized by
widespread hybridization.


  As computer scientists working in bioinformatics/computational biology, we
often face the challenge of coming up with an algorithm to answer a biological
question. This occurs in many areas, such as variant calling, alignment, and
assembly. In this tutorial, we use the example of the genome assembly problem
to demonstrate how to go from a question in the biological realm to a solution
in the computer science realm. We show the modeling process step-by-step,
including all the intermediate failed attempts.


  The use of beta-solenoid proteins as functionalizable, nanoscale,
self-assembling molecular building blocks may have many applications, including
templating the growth of wires or higher-dimensional structures. By
understanding their mechanical strengths, we can efficiently design the
proteins for specific functions. We present a study of the mechanical
properties of seven beta-solenoid proteins using GROMACS molecular dynamics
software to produce force/torque-displacement data, implement umbrella sampling
of bending/twisting trajectories, produce Potentials of Mean Force (PMFs),
extract effective spring constants, and calculate rigidities for two bending
and two twisting directions for each protein. We examine the differences
between computing the strength values from force/torque-displacement data alone
and PMF data, and show how higher precision estimates can be obtained from the
former. In addition to the analysis of the methods, we report estimates for the
bend/twist persistence lengths for each protein, which range from 0.5-3.4
$\mu$m. We note that beta-solenoid proteins with internal disulfide bridges do
not enjoy enhanced bending or twisting strength, and that the strongest
correlate with bend/twist rigidity is the number of hydrogen bonds per turn. In
addition, we compute estimates of the Young's modulus ($Y$) for each protein,
which range from $Y$ = 3.5 to 7.2 GPa.


  One of the main aims of phylogenetics is to reconstruct the \enquote{Tree of
Life}. In this respect, different methods and criteria are used to analyze DNA
sequences of different species and to compare them in order to derive the
evolutionary relationships of these species. Maximum Parsimony is one such
criterion for tree reconstruction and, it is the one which we will use in this
paper. However, it is well-known that tree reconstruction methods can lead to
wrong relationship estimates. One typical problem of Maximum Parsimony is long
branch attraction, which can lead to statistical inconsistency. In this work,
we will consider a blockwise approach to alignment analysis, namely so-called
$k$-tuple analyses. For four taxa it has already been shown that
$k$-tuple-based analyses are statistically inconsistent if and only if the
standard character-based (site-based) analyses are statistically inconsistent.
So, in the four-taxon case, going from individual sites to $k$-tuples does not
lead to any improvement. However, real biological analyses often consider more
than only four taxa. Therefore, we analyze the case of five taxa for $2$- and
$3$-tuple-site data and consider alphabets with two and four elements. We show
that the equivalence of single-site data and $k$-tuple-site data then no longer
holds. Even so, we can show that Maximum Parsimony is statistically
inconsistent for $k$-tuple site data and five taxa.


  Stability of dynamical systems is a central topic with applications in
widespread areas such as economy, biology, physics and mechanical engineering.
The dynamics of nonlinear systems may completely change due to perturbations
forcing the solution to jump from a safe state into another, possibly
dangerous, attractor. Such phenomena can not be traced by the widespread local
stability and resilience measures, based on linearizations, accounting only for
arbitrary small perturbations. Using numerical estimates of the size and shape
of the basin of attraction, as well as the systems returntime to the attractor
after given a perturbation, we construct simple nonlocal stability and
resilience measures that record a systems ability to tackle both large and
small perturbations. We demonstrate our approach on the Solow-Swan model of
economic growth, an electro-mechanical system as well as on a stage-structured
population model, and conclude that the suggested measures detect dynamic
behaviour, crucial for a systems stability and resilience, which can be
completely missed by local measures. The presented measures are also easy to
implement on a standard laptop computer. We believe that our approach will
constitute an important step towards filling a current gap in the literature by
putting forward and explaining simple ideas and methods, and by delivering
explicit constructions of several promising nonlocal stability and resilience
measures.


  We propose an alternative to the prevailing two origin of life narratives,
one based on a replicator first hypothesis, and one based on a metabolism first
hypothesis. Both hypotheses have known difficulties: All known evolvable
molecular replicators such as RNA require complex chemical (enzymatic)
machinery for the replication process. Likewise, contemporary cellular
metabolisms require several enzymatically catalyzed steps, and it is difficult
to identify a non-enzymatic path to their realization. We propose that there
must have been precursors to both replication and metabolism that enable a form
of selection to take place through action of simple chemical and physical
processes. We model a concrete example of such a process, repeated
sequestration of binary molecular combinations after exposure to an environment
with a broad distribution of chemical components, as might be realized
experimentally in in a repeated wet-dry cycle. We show that the repeated
sequestration dynamics results in a selective amplification of a very small
subset of molecular species present in the environment, thus providing a
candidate primordial selection process.


  Current understanding of the critical outbreak condition on temporal networks
relies on approximations (time scale separation, discretization) that may bias
the results. We propose a theoretical framework to compute the epidemic
threshold in continuous time through the infection propagator approach. We
introduce the {\em weak commutation} condition allowing the interpretation of
annealed networks, activity-driven networks, and time scale separation into one
formalism. Our work provides a coherent connection between discrete and
continuous time representations applicable to realistic scenarios.


  We examine a mathematical question concerning the reconstruction accuracy of
the Fitch algorithm for reconstructing the ancestral sequence of the most
recent common ancestor given a phylogenetic tree and sequence data for all taxa
under consideration. In particular, for the symmetric 4-state substitution
model which is also known as Jukes-Cantor model, we answer affirmatively a
conjecture of Li, Steel and Zhang which states that for any ultrametric
phylogenetic tree and a symmetric model, the Fitch parsimony method using all
terminal taxa is more accurate, or at least as accurate, for ancestral state
reconstruction than using any particular terminal taxon or any particular pair
of taxa. This conjecture had so far only been answered for two-state data by
Fischer and Thatte. Here, we focus on answering the biologically more relevant
case with four states, which corresponds to ancestral sequence reconstruction
from DNA or RNA data.


  Several recent studies have tackled the issue of optimal network immunization
by providing efficient criteria to identify key nodes to be removed in order to
break apart a network, thus preventing the occurrence of extensive epidemic
outbreaks. Yet, although the efficiency of those criteria has been demonstrated
also in empirical networks, preventive immunization is rarely applied to
real-world scenarios, where the usual approach is the a posteriori attempt to
contain epidemic outbreaks using quarantine measures. Here we compare the
efficiency of prevention with that of quarantine in terms of the tradeoff
between the number of removed and saved nodes on both synthetic and empirical
topologies. We show how, consistent with common sense, but contrary to common
practice, in many cases preventing is better than curing: depending on network
structure, rescuing an infected network by quarantine could become inefficient
soon after the first infection.


  To explore the mechanistic relationships between ageing, frailty and
mortality, we developed a computational model in which possible health
attributes are represented by the nodes of a complex network. Each node can be
either damaged (i.e. a deficit) or undamaged. Damage of connected nodes
facilitates further local damage. Our model recovers the known patterns of
frailty and mortality without any programmed ageing. It helps us to understand
how the observed maximum of the frailty index (FI) might arise, and allows us
to start to understand how health deficits accumulate. Large model populations
allow us to exploit new analytic tools, including information theory. This will
allow us to systematically characterize the effects of sudden changes in the
health trajectories of individuals and serve as a way to evaluate large
clinical and population databases.


  We study an optimal control problem for a non-autonomous SEIRS model with
incidence given by a general function of the infective, the susceptible and the
total population, and with vaccination and treatment as control variables. We
prove existence and uniqueness results for our problem and, for the case of
mass-action incidence, we present some simulation results designed to compare
an autonomous and corresponding periodic model, as well as the controlled
versus uncontrolled models.


  This paper formulates two 3D stochastic differential equations (SDEs) of two
microbial populations in a chemostat competing over a single substrate. The two
models have two distinct noise sources. One is general noise whereas the other
is dilution rate induced noise. Nonlinear Monod growth rates are assumed and
the paper is mainly focused on the parameter values where coexistence is
present deterministically. Nondimensionalising the equations around the point
of intersection of the two growth rates leads to a large parameter which is the
nondimensional substrate feed. This in turn is used to perform an asymptotic
analysis leading to a reduced 2D system of equations describing the dynamics of
the populations on and close to a line of steady states retrieved from the
deterministic stability analysis. That reduced system allows the formulation of
a spatially 2D Fokker-Planck equation which when solved numerically admits
results similar to those from simulation of the SDEs. Contrary to previous
suggestions, one particular population becomes dominant at large times.
Finally, we brie y explore the case where death rates are added.


  Different subsystems of organisms adapt over many time scales, such as rapid
changes in the nervous system (learning), slower morphological and neurological
change over the lifetime of the organism (postnatal development), and change
over many generations (evolution). Much work has focused on instantiating
learning or evolution in robots, but relatively little on development. Although
many theories have been forwarded as to how development can aid evolution, it
is difficult to isolate each such proposed mechanism. Thus, here we introduce a
minimal yet embodied model of development: the body of the robot changes over
its lifetime, yet growth is not influenced by the environment. We show that
even this simple developmental model confers evolvability because it allows
evolution to sweep over a larger range of body plans than an equivalent
non-developmental system, and subsequent heterochronic mutations 'lock in' this
body plan in more morphologically-static descendants. Future work will involve
gradually complexifying the developmental model to determine when and how such
added complexity increases evolvability.


  Living cells respond to spatial signals. Signal transmission to the cell
interior often involves the release of second messengers like $Ca^{2+}$ . They
will eventually trigger a physiological response by activating kinases that in
turn activate target proteins through phosphorylation. Here, we investigate
theoretically how positional information can be accurately read out by protein
phosphorylation in spite of rapid second messenger diffusion. We find that
accuracy is increased by binding of the kinases to the cell membrane prior to
phosphorylation and by increasing the rate of $Ca^{2+}$ loss from the cell
interior. These findings could explain some salient features of conventional
protein kinases C.


  In biology phenotypic switching is a common bet-hedging strategy in the face
of uncertain environmental conditions. Existing mathematical models often focus
on periodically changing environments to determine the optimal phenotypic
response. We focus on the case in which the environment switches randomly
between discrete states. Starting from an individual-based model we derive
stochastic differential equations to describe the dynamics, and obtain
analytical expressions for the mean instantaneous growth rates based on the
theory of piecewise deterministic Markov processes. We show that optimal
phenotypic responses are non-trivial for slow and intermediate environmental
processes, and systematically compare the cases of periodic and random
environments. The best response to random switching is more likely to be
heterogeneity than in the case of deterministic periodic environments, net
growth rates tend to be higher under stochastic environmental dynamics. The
combined system of environment and population of cells can be interpreted as
host-pathogen interaction, in which the host tries to choose environmental
switching so as to minimise growth of the pathogen, and in which the pathogen
employs a phenotypic switching optimised to increase its growth rate. We
discuss the existence of Nash-like mutual best-response scenarios for such
host-pathogen games.


  The Flow Decomposition problem, which asks for the smallest set of weighted
paths that "covers" a flow on a DAG, has recently been used as an important
computational step in transcript assembly. We prove the problem is in FPT when
parameterized by the number of paths by giving a practical linear fpt
algorithm. Further, we implement and engineer a Flow Decomposition solver based
on this algorithm, and evaluate its performance on RNA-sequence data.
Crucially, our solver finds exact solutions while achieving runtimes
competitive with a state-of-the-art heuristic. Finally, we contextualize our
design choices with two hardness results related to preprocessing and weight
recovery. Specifically, $k$-Flow Decomposition does not admit polynomial
kernels under standard complexity assumptions, and the related problem of
assigning (known) weights to a given set of paths is NP-hard.


  In an iterated two-person game, for instance prisoner's dilemma or the
snowdrift game, there exist strategies that force the payoffs of the opponents
to be equal. These equalizer strategies form a subset of the more general
zero-determinant strategies that unilaterally set the payoff of an opponent. A
challenge in the attempts to understand the role of these strategies in the
evolution of animal behavior is the lack of iterations in the fights for mating
opportunities or territory control. We show that an arbitrary two-parameter
strategy may possess a corresponding equalizer strategy which produces the same
result: statistics of the fight outcomes in the contests with mutants are the
same for each of these two strategies. Therefore, analyzing only the equalizer
strategy space may be sufficient to predict animal behavior if nature, indeed,
reduces (marginalizes) complex strategies to equalizer strategy space. The
work's main finding is that there is a unique equalizer strategy that predicts
fight outcomes without mutual cooperation. The lack of mutual cooperation is a
common trait in conflict escalation contests that generally require a clear
winner. In addition this unique strategy does not assess information of the
opponent's state. The method bypasses the standard analysis of evolutionary
stability. The results fit well the observations of combat between male bowl
and doily spiders and support an empirical assumption of the war of attrition
model that the species use only information regarding their own state during
conflict escalation.


  During evolution of microorganisms genomes underwork have different changes
in their lengths, gene orders, and gene contents. Investigating these
structural rearrangements helps to understand how genomes have been modified
over time. Some elements that play an important role in genome rearrangements
are called insertion sequences (ISs), they are the simplest types of
transposable elements (TEs) that widely spread within prokaryotic genomes. ISs
can be defined as DNA segments that have the ability to move (cut and paste)
themselves to another location within the same chromosome or not. Due to their
ability to move around, they are often presented as responsible of some of
these genomic recombination. Authors of this research work have regarded this
claim, by checking if a relation between insertion sequences (ISs) and genome
rearrangements can be found. To achieve this goal, a new pipeline that combines
various tools have firstly been designed, for detecting the distribution of
ORFs that belongs to each IS category. Secondly, links between these predicted
ISs and observed rearrangements of two close genomes have been investigated, by
seeing them with the naked eye, and by using computational approaches. The
proposal has been tested on 18 complete bacterial genomes of Pseudomonas
aeruginosa, leading to the conclusion that IS3 family of insertion sequences
are related to genomic inversions.


  Variability in drug efficacy and adverse effects are observed in clinical
practice. While the extent of genetic variability in classical pharmacokinetic
genes is rather well understood, the role of genetic variation in drug targets
is typically less studied. Based on 60,706 human exomes from the ExAC dataset,
we performed an in-depth computational analysis of the prevalence of
functional-variants in in 806 drug-related genes, including 628 known drug
targets. We find that most genetic variants in these genes are very rare (f <
0.1%) and thus likely not observed in clinical trials. Overall, however, four
in five patients are likely to carry a functional-variant in a target for
commonly prescribed drugs and many of these might alter drug efficacy. We
further computed the likelihood of 1,236 FDA approved drugs to be affected by
functional-variants in their targets and show that the patient-risk varies for
many drugs with respect to geographic ancestry. A focused analysis of
oncological drug targets indicates that the probability of a patient carrying
germline variants in oncological drug targets is with 44% high enough to
suggest that not only somatic alterations, but also germline variants carried
over into the tumor genome should be included in therapeutic decision-making.


  Insertion Sequences (ISs) are small DNA segments that have the ability of
moving themselves into genomes. These types of mobile genetic elements (MGEs)
seem to play an essential role in genomes rearrangements and evolution of
prokaryotic genomes, but the tools that deal with discovering ISs in an
efficient and accurate way are still too few and not totally precise. Two main
factors have big effects on IS discovery, namely: genes annotation and
functionality prediction. Indeed, some specific genes called "transposases" are
enzymes that are responsible of the production and catalysis for such
transposition, but there is currently no fully accurate method that could
decide whether a given predicted gene is either a real transposase or not. This
is why authors of this article aim at designing a novel pipeline for ISs
detection and classification, which embeds the most recently available tools
developed in this field of research, namely OASIS (Optimized Annotation System
for Insertion Sequence) and ISFinder database (an up-to-date and accurate
repository of known insertion sequences). As this latter depend on predicted
coding sequences, the proposed pipeline will encompass too various kinds of
bacterial genes annotation tools (that is, Prokka, BASys, and Prodigal). A
complete IS detection and classification pipeline is then proposed and tested
on a set of 23 complete genomes of Pseudomonas aeruginosa. This pipeline can
also be used as an investigator of annotation tools performance, which has led
us to conclude that Prodigal is the best software for IS prediction. A deepen
study regarding IS elements in P.aeruginosa has then been conducted, leading to
the conclusion that close genomes inside this species have also a close numbers
of IS families and groups.