  A Restricted Boltzmann Machine (RBM) is an unsupervised machine-learning
bipartite graphical model that jointly learns a probability distribution over
data and extracts their relevant statistical features. As such, RBM were
recently proposed for characterizing the patterns of coevolution between amino
acids in protein sequences and for designing new sequences. Here, we study how
the nature of the features learned by RBM changes with its defining parameters,
such as the dimensionality of the representations (size of the hidden layer)
and the sparsity of the features. We show that for adequate values of these
parameters, RBM operate in a so-called compositional phase in which visible
configurations sampled from the RBM are obtained by recombining these features.
We then compare the performance of RBM with other standard representation
learning algorithms, including Principal or Independent Component Analysis,
autoencoders (AE), variational auto-encoders (VAE), and their sparse variants.
We show that RBM, due to the stochastic mapping between data configurations and
representations, better capture the underlying interactions in the system and
are significantly more robust with respect to sample size than deterministic
methods such as PCA or ICA. In addition, this stochastic mapping is not
prescribed a priori as in VAE, but learned from data, which allows RBM to show
good performance even with shallow architectures. All numerical results are
illustrated on synthetic lattice-protein data, that share similar statistical
features with real protein sequences, and for which ground-truth interactions
are known.


  Population structure can have a significant effect on evolution. For some
systems with sufficient symmetry, analytic results can be derived within the
mathematical framework of evolutionary graph theory which relate to the outcome
of the evolutionary process. However, for more complicated heterogeneous
structures, computationally intensive methods are required such as
individual-based stochastic simulations. By adapting methods from statistical
physics, including moment closure techniques, we first show how to derive
existing homogenised pair approximation models and the exact neutral drift
model. We then develop node-level approximations to stochastic evolutionary
processes on arbitrarily complex structured populations represented by finite
graphs, which can capture the different dynamics for individual nodes in the
population. Using these approximations, we evaluate the fixation probability of
invading mutants for given initial conditions, where the dynamics follow
standard evolutionary processes such as the invasion process. Comparisons with
the output of stochastic simulations reveal the effectiveness of our
approximations in describing the stochastic processes and in predicting the
probability of fixation of mutants on a wide range of graphs. Construction of
these models facilitates a systematic analysis and is valuable for a greater
understanding of the influence of population structure on evolutionary
processes.


  Genetic inactivation of essential genes creates an evolutionary scenario
distinct from escape from drug inhibition, but the mechanisms of microbe
adaptations in such cases remain unknown. Here we inactivate E. coli
dihydrofolate reductase (DHFR) by introducing D27G,N,F chromosomal mutations in
a key catalytic residue with subsequent adaptation by serial dilutions. The
partial reversal G27->C occurred in three evolutionary trajectories.
Conversely, in one trajectory for D27G and in all trajectories for D27F,N
strains adapted to grow at very low supplement folAmix concentrations but did
not escape entirely from supplement auxotrophy. Major global shifts in
metabolome and proteome occurred upon DHFR inactivation, which were partially
reversed in adapted strains. Loss of function mutations in two genes, thyA and
deoB, ensured adaptation to low folAmix by rerouting the
2-Deoxy-D-ribose-phosphate metabolism from glycolysis towards synthesis of
dTMP. Multiple evolutionary pathways of adaptation to low folAmix converge to
highly accessible yet suboptimal fitness peak.


  The ever-increasing number of genome sequencing and resequencing projects is
a central source of insights into the ecology and evolution of non-model
organisms. An important aspect of genomics is the elucidation of sex
determination systems and identifying genes on sex chromosomes. This not only
helps reveal mechanisms behind sex determination in the species under study,
but their characteristics make sex chromosomes a unique tool for studying the
mechanisms and effects of recombination and genomic rearrangements and how they
affect adaption and selection. Despite this, many sequencing projects omit such
investigations. Here, we apply a simple method using sequencing coverage
statistics to identify scaffolds belonging to the sex chromosomes of minke
whale, and show how the sex chromosome system can be determined using coverage
statistics alone. Using publicly available data, we identify the previously
unknown sex of an Antarctic minke whale as female. We further investigate
public sequence data from the different species and sub-species of minke whale,
and classify genomic scaffolds from a published minke whale assembly as X or Y
chromosomal sequences. Our findings are consistent with previous results that
identified a handful of scaffolds as sex chromosomal, but we are able to
identify a much larger set of scaffolds, likely to represent close to the
complete sex chromosomal sequences for the minke whale. Sequence coverage
statistics provides a readily available tool for investigating the sex
determination system and locate genes on sex chromosomes. This analysis is
straightforward and can often be performed with existing resources.


  Invasive alien species are a growing threat for environment and health. They
also have a major economic impact, as they can damage many infrastructures. The
Japanese knotweed (Fallopia japonica), present in North America, Northern and
Central Europe as well as in Australia and New Zealand, is listed by the World
Conservation Union as one of the world's worst invasive species. So far, most
models have dealt with how the invasion spreads without management. This paper
aims at providing a model able to study and predict the dynamics of a stand of
Japanese knotweed taking into account mowing as a management technique. The
model we propose is stochastic and individual-based, which allows us taking
into account the behaviour of individuals depending on their size and location,
as well as individual stochasticity. We set plant dynamics parameters thanks to
a calibration with field data, and study the influence of the initial
population size, the mean number of mowing events a year and the management
project duration on mean area and mean number of crowns of stands. In
particular, our results provide the sets of parameters for which it is possible
to obtain the stand eradication, and the minimal duration of the management
project necessary to achieve this latter.


  One of the most difficult problems difficult problem in systems biology is to
discover protein-protein interactions as well as their associated functions.
The analysis and alignment of protein-protein interaction networks (PPIN),
which are the standard model to describe protein-protein interactions, has
become a key ingredient to obtain functional orthologs as well as evolutionary
conserved pathways and protein complexes. Several methods have been proposed to
solve the PPIN alignment problem, aimed to match conserved subnetworks or
functionally related proteins. However, the right balance between considering
network topology and biological information is one of the most difficult and
key points in any PPIN alignment algorithm which, unfortunately, remains
unsolved. Therefore, in this work, we propose AligNet, a new method and
software tool for the pairwise global alignment of PPIN that produces
biologically meaningful alignments and more efficient computations than
state-of-the-art methods and tools, by achieving a good balance between
structural matching and protein function conservation as well as reasonable
running times.


  Quantitative studies of cell metabolism are often based on large chemical
reaction network models. A steady state approach is suited to analyze phenomena
on the timescale of cell growth and circumvents the problem of incomplete
experimental knowledge on kinetic laws and parameters, but it shall be
supported by a correct implementation of thermodynamic constraints. In this
article we review the latter aspect highlighting its computational challenges
and physical insights. The simple introduction of Gibbs inequalities avoids the
presence of unfeasible loops allowing for correct timescale analysis but leads
to possibly non-convex feasible flux spaces, whose exploration needs efficient
algorithms. We shorty review on the implementation of thermodynamics through
variational principles in constraints based models of metabolic networks.


  To reduce experimental effort associated with directed protein evolution and
to explore the sequence space encoded by mutating multiple positions
simultaneously, we incorporate machine learning in the directed evolution
workflow. Combinatorial sequence space can be quite expensive to sample
experimentally, but machine learning models trained on tested variants provide
a fast method for testing sequence space computationally. We validate this
approach on a large published empirical fitness landscape for human GB1 binding
protein, demonstrating that machine learning-guided directed evolution finds
variants with higher fitness than those found by other directed evolution
approaches. We then provide an example application in evolving an enzyme to
produce each of the two possible product enantiomers (stereodivergence) of a
new-to-nature carbene Si-H insertion reaction. The approach predicted libraries
enriched in functional enzymes and fixed seven mutations in two rounds of
evolution to identify variants for selective catalysis with 93% and 79% ee. By
greatly increasing throughput with in silico modeling, machine learning
enhances the quality and diversity of sequence solutions for a protein
engineering problem.


  Mutational robustness quantifies the effect of random mutations on fitness.
When mutational robustness is high, most mutations do not change fitness or
have only a minor effect on it. From the point of view of fitness landscapes,
robust genotypes form neutral networks of almost equal fitness. Using
deterministic population models it has been shown that selection favors
genotypes inside such networks, which results in increased mutational
robustness. Here we demonstrate that this effect is massively enhanced by
recombination. Our results are based on a detailed analysis of mesa-shaped
fitness landscapes, where we derive precise expressions for the dependence of
the robustness on the landscape parameters for recombining and non-recombining
populations. In addition, we carry out numerical simulations on different types
of random holey landscapes as well as on an empirical fitness landscape. We
show that the mutational robustness of a genotype generally correlates with its
recombination weight, a new measure that quantifies the likelihood for the
genotype to arise from recombination. We argue that the favorable effect of
recombination on mutational robustness is a highly universal feature that may
have played an important role in the emergence and maintenance of mechanisms of
genetic exchange.


  The envelopes covering bacterial cytoplasm possess remarkable elastic
properties. They are rigid enough to resist large pressures while being
flexible enough to adapt to growth under environmental constraints. Similarly,
the virus shells play an important role in their functions. However, the
effects of mechanical properties of the outer shell in controlling and
maintaining the sizes of bacteria or viruses are unknown. Here, we present a
hydrodynamic "bubbles with shell" model, motivated by the study of bubble
stability in fluids, to demonstrate that shell rigidity and turgor pressure
control the sizes of bacteria and viruses. A dimensionless compliance
parameter, expressed in terms of the elastic modulus of the shell, its
thickness and the turgor pressure, determines membrane response to deformation
and the size of the organisms. By analyzing experiment data, we show that
bacterial and viral sizes correlate with shell elasticity, which plays a
critical role in regulating size.


  Biomolecular oscillators can function robustly in the presence of
environmental perturbations, which can either be static or dynamic. While the
effect of different circuit parameters and mechanisms on the robustness to
steady perturbations has been investigated, the scenario for dynamic
perturbations is relatively unclear. To address this we use a benchmark three
protein oscillator design - the repressilator - and investigate its robustness
to pulse perturbations, computationally as well as using analytical tools of
Floquet theory. We find that the metric provided by direct computations of the
time it takes for the oscillator to settle after a pulse perturbation is
applied, correlates well with the metric provided by Floquet theory. We
investigate the parametric dependence of the Floquet metric, finding that the
parameters that increase the effective delay enhance robustness to pulse
perturbation. We find that the structural changes such as increasing the number
of proteins in a ring oscillator as well as adding positive feedback, both of
which increase effective delay, facilitates such robustness. These results
highlight such design principles, especially the role of delay, for designing
an oscillator that is robust to pulse perturbation.


  A theoretical model is presented for early evolutionary cell sorting within
cellular aggregates. The model involves an energy-saving mechanism and
principles of collective self-organization analogous to those observed in
bicycle pelotons (groups of cyclists). The theoretical framework is applied to
slime-mold slugs (Dictyostelium discoideum) and incorporated into a computer
simulation which demonstrates principally the sorting of cells between the
anterior and posterior slug regions. The simulation relies on an existing
simulation of bicycle peloton dynamics which is modified to incorporate a
limited range of cell metabolic capacities among heterogeneous cells, along
with a tunable energy-expenditure parameter, referred to as an 'output-level'
or 'starvation-level' to reflect diminishing energetic supply, proto-cellular
dynamics are modelled for three output phases: 'active', 'suffering', and
'dying or dead.' Adjusting the starvation parameter causes cell differentiation
and sorting into sub-groups within the cellular aggregate. Tuning of the
starvation parameter demonstrates how weak or expired cells shuffle backward
within the cellular aggregate.


  There are known expressions to calculate the moments of the first passage
time in Markov chains. Nevertheless, it is commonly forgotten that in most
applications the parameters of the Markov chain are constructed using estimates
based upon empirical data and in those cases the data sample size should play
an important role in estimating the variance. Here we provide a Monte Carlo
approach to estimate the first two moments of the passage time in this
situation. We illustrate this method with an example using data from the
biological field.


  The importance of semi-arid ecosystems in the global carbon cycle as sinks
for CO2 emissions has recently been highlighted. Africa is a carbon sink and
nearly half its area comprises arid and semi-arid ecosystems. However, there
are uncertainties regarding CO2 fluxes for semi-arid ecosystems in Africa,
particularly savannas and dry tropical woodlands. In order to improve on
existing remote-sensing based methods for estimating carbon uptake across
semi-arid Africa we applied and tested the recently developed plant phenology
index (PPI). We developed a PPI-based model estimating gross primary
productivity (GPP) that accounts for canopy water stress, and compared it
against three other Earth observation-based GPP models: the temperature and
greenness model, the greenness and radiation model and a light use efficiency
model. The models were evaluated against in situ data from four semi-arid sites
in Africa with varying tree canopy cover (3 to 65 percent). Evaluation results
from the four GPP models showed reasonable agreement with in situ GPP measured
from eddy covariance flux towers (EC GPP) based on coefficient of variation,
root-mean-square error, and Bayesian information criterion. The PPI-based GPP
model was able to capture the magnitude of EC GPP better than the other tested
models. The results of this study show that a PPI-based GPP model is a
promising tool for the estimation of GPP in the semi-arid ecosystems of Africa.


  We present a Bayesian hierarchical multi-view mixture model termed Symphony
that simultaneously learns clusters of cells representing cell types and their
underlying gene regulatory networks by integrating data from two views:
single-cell gene expression data and paired epigenetic data, which is
informative of gene-gene interactions. This model improves interpretation of
clusters as cell types with similar expression patterns as well as regulatory
networks driving expression, by explaining gene-gene covariances with the
biological machinery regulating gene expression. We show the theoretical
advantages of the multi-view learning approach and present a Variational EM
inference procedure. We demonstrate superior performance on both synthetic data
and real genomic data with subtypes of peripheral blood cells compared to other
methods.


  Thermoregulation in honey bee colonies during winter is thought to be
self-organised. We added mortality of individual honey bees to an existing
model of thermoregulation to account for elevated losses of bees that are
reported worldwide. The aim of analysis is to obtain a better fundamental
understanding of the consequences of individual mortality during winter. This
model resembles the well-known Keller-Segel model. In contrast to the often
studied Keller-Segel models, our model includes a chemotactic coefficient of
which the sign can change as honey bees have a preferred temperature: when the
local temperature is too low, they move towards higher temperatures, whereas
the opposite is true for too high temperatures. Our study shows that we can
distinguish two states of the colony: one in which the colony size is above a
certain critical number of bees in which the bees can keep the core temperature
of the colony above the threshold temperature, and one in which the core
temperature drops below the critical threshold and the mortality of the bees
increases dramatically, leading to a sudden death of the colony. This model
behaviour may explain the globally observed honey bee colony losses during
winter.


  The intricate shapes of biological membranes such as tubules and membrane
stacks are induced by proteins. In this article, we systematically investigate
the membrane shapes induced by arc-shaped scaffolds such as proteins and
protein complexes with coarse-grained modeling and simulations. We find that
arc-shaped scaffolds induce membrane tubules at membrane coverages larger than
a threshold of about 40%, irrespective of their arc angle. The membrane
morphologies at intermediate coverages below this tubulation threshold, in
contrast, strongly depend on the arc angle. Scaffolds with arc angles of about
60 degree akin to N-BAR domains do not change the membrane shape at coverages
below the tubulation threshold, while scaffolds with arc angles larger than
about 120 degree induce double-membrane stacks at intermediate coverages. The
scaffolds stabilize the curved membrane edges that connect the membrane stacks,
as suggested for complexes of reticulon proteins. Our results provide general
insights on the determinants of membrane shaping by arc-shaped scaffolds.


  Deciding whether a Markov matrix is embeddable (i.e. can be written as the
exponential of a rate matrix) is an open problem even for $4\times 4$ matrices.
We study the embedding problem and rate identifiability for the K80 model of
nucleotide substitution. For these $4\times 4$ matrices, we fully characterize
the set of embeddable K80 Markov matrices and the set of embeddable matrices
for which rates are identifiable. In particular, we describe an open subset of
embeddable matrices with non-identifiable rates. This set contains matrices
with positive eigenvalues and also diagonal largest in column matrices, which
might lead to consequences in parameter estimation in phylogenetics. Finally,
we compute the relative volumes of embeddable K80 matrices and of embeddable
matrices with identifiable rates. This study concludes the embedding problem
for the more general model K81 and its submodels, which had been initiated by
the last two authors in a separate work.


  Inferring the structural properties of a protein from its amino acid sequence
is a challenging yet important problem in biology. Structures are not known for
the vast majority of protein sequences, but structure is critical for
understanding function. Existing approaches for detecting structural similarity
between proteins from sequence are unable to recognize and exploit structural
patterns when sequences have diverged too far, limiting our ability to transfer
knowledge between structurally related proteins. We newly approach this problem
through the lens of representation learning. We introduce a framework that maps
any protein sequence to a sequence of vector embeddings --- one per amino acid
position --- that encode structural information. We train bidirectional long
short-term memory (LSTM) models on protein sequences with a two-part feedback
mechanism that incorporates information from (i) global structural similarity
between proteins and (ii) pairwise residue contact maps for individual
proteins. To enable learning from structural similarity information, we define
a novel similarity measure between arbitrary-length sequences of vector
embeddings based on a soft symmetric alignment (SSA) between them. Our method
is able to learn useful position-specific embeddings despite lacking direct
observations of position-level correspondence between sequences. We show
empirically that our multi-task framework outperforms other sequence-based
methods and even a top-performing structure-based alignment method when
predicting structural similarity, our goal. Finally, we demonstrate that our
learned embeddings can be transferred to other protein sequence problems,
improving the state-of-the-art in transmembrane domain prediction.


  Background: Mechanotransduction in bone cells plays a pivotal role in
osteoblast differentiation and bone remodelling. Mechanotransduction provides
the link between modulation of the extracellular matrix and intracellular
actions. By controlling the balance between the intracellular and extracellular
domains, the mechanotransduction process determines optimal functionality of
the skeletal dynamics, and it is one of the possible causes of
osteophatological diseases.
  Results: Mechanotransduction in a single osteoblast under external mechanical
perturbations has been modelled in the agent based framework to reproduce the
dynamics of the stochastic reaction diffusion process among molecules in the
cytoplasm, nuclear and extracellular domains. The amount of molecules and
fluctuations of each molecular class has been analysed in terms of recurrences
of critical events. A numerical approach has been developed to invert
subordination processes and to extract the directing processes from the
molecular signals in order to derive the distribution of recurrence of events.
  Conclusions: We observed large fluctuations enclose information hidden in the
noise which is beyond the dynamic variations of molecular baselines. Studying
the system under different parametric conditions and stimuli, the results have
shown that the waiting time distributions of each molecule are a signature of
the system's dynamics. The behaviours of the molecular waiting times change
with the changing of parameters presenting the same variation of patterns for
similar interacting molecules and identifying specific alterations for key
molecules in the mechanotransduction pathway.


  'Evolutionary rescue' is the potential for evolution to enable population
persistence in a changing environment. Even with eventual rescue, evolutionary
time lags can cause the population size to temporarily fall below a threshold
susceptible to extinction. To reduce extinction risk given human-driven global
change, conservation management can enhance populations through actions such as
captive breeding. To quantify the optimal timing of, and indicators for
engaging in, investment in temporary enhancement to enable evolutionary rescue,
we construct a model of coupled demographic-genetic dynamics given a moving
optimum. We assume 'decelerating change', as might be relevant to climate
change, where the rate of environmental change initially exceeds a rate where
evolutionary rescue is possible, but eventually slows. We analyze the optimal
control path of an intervention to avoid the population size falling below a
threshold susceptible to extinction, minimizing costs. We find that the optimal
path of intervention initially increases as the population declines, then
declines and ceases when the population growth rate becomes positive, which
lags the stabilization in environmental change. In other words, the optimal
strategy involves increasing investment even in the face of a declining
population, and positive population growth could serve as a signal to end the
intervention. In addition, a greater carrying capacity relative to the initial
population size decreases the optimal intervention. Therefore, a one-time
action to increase carrying capacity, such as habitat restoration, can reduce
the amount and duration of longer-term investment in population enhancement,
even if the population is initially lower than and declining away from the new
carrying capacity.


  Phylodynamics is a set of population genetics tools that aim at
reconstructing demographic history of a population based on molecular sequences
of individuals sampled from the population of interest. One important task in
phylodynamics is to estimate changes in (effective) population size. When
applied to infectious disease sequences such estimation of population size
trajectories can provide information about changes in the number of infections.
To model changes in the number of infected individuals, current phylodynamic
methods use non-parametric approaches, parametric approaches, and stochastic
modeling in conjunction with likelihood-free Bayesian methods. The first class
of methods yields results that are hard-to-interpret epidemiologically. The
second class of methods provides estimates of important epidemiological
parameters, such as infection and removal/recovery rates, but ignores variation
in the dynamics of infectious disease spread. The third class of methods is the
most advantageous statistically, but relies on computationally intensive
particle filtering techniques that limits its applications. We propose a
Bayesian model that combines phylodynamic inference and stochastic epidemic
models, and achieves computational tractability by using a linear noise
approximation (LNA) --- a technique that allows us to approximate probability
densities of stochastic epidemic model trajectories. LNA opens the door for
using modern Markov chain Monte Carlo tools to approximate the joint posterior
distribution of the disease transmission parameters and of high dimensional
vectors describing unobserved changes in the stochastic epidemic model
compartment sizes (e.g., numbers of infectious and susceptible individuals). We
apply our estimation technique to Ebola genealogies estimated using viral
genetic data from the 2014 epidemic in Sierra Leone and Liberia.


  We study the early stages of viral infection, and the distribution of times
to obtain a persistent infection. The virus population proliferates by entering
and reproducing inside a target cell until a sufficient number of new virus
particles are released via a burst, with a given burst size distribution, which
results in the death of the infected cell. Starting with a 2D model describing
the joint dynamics of the virus and infected cell populations, we analyze the
corresponding master equation using the probability generating function
formalism. Exploiting time-scale separation between the virus and infected cell
dynamics, the 2D model can be cast into an effective 1D model. To this end, we
solve the 1D model analytically for a particular choice of burst size
distribution. In the general case, we solve the model numerically by performing
extensive Monte-Carlo simulations, and demonstrate the equivalence between the
2D and 1D models by measuring the Kullback-Leibler divergence between the
corresponding distributions. Importantly, we find that the distribution of
infection times is highly skewed with a "fat" exponential right tail. This
indicates that there is non-negligible portion of individuals with an infection
time, significantly longer than the mean, which may have implications on when
HIV tests should be performed.


  Phylogenetic networks generalize phylogenetic trees by allowing the
modelization of events of reticulate evolution. Among the different kinds of
phylogenetic networks that have been proposed in the literature, the subclass
of binary tree-child networks is one of the most studied ones. However, very
little is known about the combinatorial structure of these networks. In this
paper we address the problem of generating all possible binary tree-child
networks with a given number of leaves in an efficient way via
reduction/augmentation operations that extend and generalize analogous
operations for phylogenetic trees and are biologically relevant. Since our
solution is recursive, this also provides us with a recurrence relation giving
an upper bound on the number of such networks.


  Microbial communities feature an immense diversity of species and this
diversity is linked with outcomes ranging from ecosystem stability to medical
prognoses. Yet the mechanisms underlying microbial diversity are under debate.
While simple resource-competition models don't allow for coexistence of a large
number of species, it was recently shown that metabolic trade-offs can allow
unlimited diversity. Does this diversity persist with more realistic,
intermittent nutrient supply? Here, we demonstrate theoretically that in serial
dilution culture, metabolic trade-offs allow for high diversity. When a small
amount of nutrient is supplied to each batch, the serial dilution dynamics
mimic a chemostat-like steady state. If more nutrient is supplied, diversity
depends on the amount of nutrient supplied due to an "early-bird" effect. The
interplay of this effect with different environmental factors and
diversity-supporting mechanisms leads to a variety of relationships between
nutrient supply and diversity, suggesting that real ecosystems may not obey a
universal nutrient-diversity relationship.


  The light-based minimum-time circadian entrainment problem for mammals,
Neurospora, and Drosophila is studied based on the mathematical models of their
circadian gene regulation. These models contain high order nonlinear
differential equations. Two model simplification methods are applied to these
high-order models: the phase response curves (PRC) and the Principal Orthogonal
Decomposition (POD). The variational calculus and a gradient descent algorithm
are applied for solving the optimal light input in the high-order models. As
the results of the gradient descent algorithm rely heavily on the initial
guesses, we use the optimal control of the PRC and the simplified model to
initialize the gradient descent algorithm. In this paper, we present: (1) the
application of PRC and direct shooting algorithm on high-order nonlinear
models; (2) a general process for solving the minimum-time optimal control
problem on high-order models; (3) the impacts of minimum-time optimal light on
circadian gene transcription and protein synthesis.


  The Third-Generation in DNA sequencing has emerged in the last few years
using new technologies that allow the production of long-read sequences.
Applications of the Third-Generation sequencing enable real-time and on-site
data production, changing the research paradigms in environmental and medical
sampling in virology. To take full advantage of large-scale data generated from
long-read sequencing, an innovation in the downstream data analysis is
necessary. Here, we discuss futuristic methods using machine learning
approaches to analyze big genetic data. Machine learning combines pattern
recognition and computational learning to perform predictive and exploratory
data analysis. In particular, deep learning is a field of machine learning that
is used to solve complex problems through artificial neural networks. Unlike
other methods, features can be learned using neural networks entirely from data
without manual specifications. We discuss the future of 21st-century virology
by presenting futuristic approaches for virus studies using real-time data
production and on-site data analysis with the Third-Generation Sequencing and
machine learning methods. We first introduce the basic concepts in conventional
statistical models and methods in virology, building gradually into the
necessity of innovating the downstream data analysis to meet the advances in
sequencing technologies. We argue that artificial neural networks can innovate
the downstream data analysis, as they can learn from big datasets without model
assumptions or feature specifications, as opposed to the current data analysis
in bioinformatics. Furthermore, we discuss how futuristic methods using
artificial neural networks combined with long-read sequences can revolutionize
virus studies, using specific examples in supervised and unsupervised settings.


  Diversity is a central concept in many fields. Despite its importance, there
is no unified methodological framework to measure diversity and its three
components of variety, balance and disparity. Current approaches take into
account disparity of the types by considering their pairwise similarities.
Pairwise similarities between types do not adequately capture total disparity,
since they fail to take into account in which way pairs are similar. Hence,
pairwise similarities do not discriminate between similarity of types in terms
of the same feature and similarity of types in terms of different features.
This paper presents an alternative approach which is based similarities of
features between types over the whole set. The proposed measure of diversity
properly takes into account the aspects of variety, balance and disparity, and
without having to set an arbitrary weight for each aspect of diversity. Based
on this measure, the 'ABC decomposition' is introduced, which provides separate
measures for the variety, balance and disparity, allowing them to enter
analysis separately. The method is illustrated by analyzing the industrial
diversity from 1850 to present while taking into account the overlap in
occupations they employ. Finally, the framework is extended to take into
account disparity considering multiple features, providing a helpful tool in
analysis of high-dimensional data.


  Living cells often need to measure chemical concentrations that vary in time.
To this end, they deploy many resources, e.g. receptors, downstream signaling
molecules, time and energy. Here, we present a theory for the optimal design of
a large class of sensing systems that need to detect time-varying signals, a
receptor driving a push-pull network. The theory is based on the concept of the
dynamic input-output relation, which describes the mapping between the current
ligand concentration and the average receptor occupancy over the past
integration time. This concept is used to develop the idea that the cell
employs its push-pull network to estimate the receptor occupancy and then uses
this estimate to infer the current ligand concentration by inverting the
dynamic input-output relation. The theory reveals that the sensing error can be
decomposed into two terms: the sampling error in the estimate of the receptor
occupancy and the dynamical error that arises because the average ligand
concentration over the past integration time may not reflect the current ligand
concentration. The theory generalizes the design principle of optimal resource
allocation previously identified for static signals, which states that in an
optimally designed sensing system the three fundamental resource classes of
receptors and their integration time, readout molecules, and energy are equally
limiting so that no resource is wasted. However, in contrast to static signals,
receptors and power cannot be traded freely against time to reach a desired
sensing precision: there exists an optimal integration time that maximizes the
sensing precision, which depends on the number of receptors, the receptor
correlation time, and the correlation time and variance of the input signal.
Applying our theory to the chemotaxis system of Escherichia coli indicates that
this bacterium has evolved to optimally sense shallow gradients.


  Partial clonality is widespread across the tree of life, but most population
genetics models are designed for exclusively clonal or sexual organisms. This
gap hampers our understanding of the influence of clonality on evolutionary
trajectories and the interpretation of population genetics data. We performed
forward simulations of diploid populations at increasing rates of clonality
(c), analysed their relationships with genotypic (clonal richness, R, and
distribution of clonal sizes, Pareto \beta) and genetic (FIS and linkage
disequilibrium) indices, and tested predictions of c from population genetics
data through supervised machine learning. Two complementary behaviours emerged
from the probability distributions of genotypic and genetic indices with
increasing c. While the impact of c on R and Pareto \beta was easily described
by simple mathematical equations, its effects on genetic indices were
noticeable only at the highest levels (c>0.95). Consequently, genotypic indices
allowed reliable estimates of c, while genetic descriptors led to poorer
performances when c<0.95. These results provide clear baseline expectations for
genotypic and genetic diversity and dynamics under partial clonality.
Worryingly, however, the use of realistic sample sizes to acquire empirical
data systematically led to gross underestimates (often of one to two orders of
magnitude) of c, suggesting that many interpretations hitherto proposed in the
literature, mostly based on genotypic richness, should be reappraised. We
propose future avenues to derive realistic confidence intervals for c and show
that, although still approximate, a supervised learning method would greatly
improve the estimation of c from population genetics data.


  Computation of biological processes creates great promise for everyday life
and great challenges for physical scientists. Simulations of molecular dynamics
appeal to biologists as a natural extension of structural biology. Once
biologists see a structure, they want to see it move. Molecular biology shows
that a few atoms, often messenger ions like Ca$^{2+}$, control biological
function on the scale of cells, sometimes organisms. Enormously concentrated
ions (~20 M) in protein channels and enzymes can control important
characteristics of living systems, just as highly concentrated ions near
electrodes control important characteristics of electrochemical systems. The
scale differences needed to simulate all the atoms of biological cells are
$10^7$ in linear dimension, $10^{21}$ in three dimensions, $10^{9}$ in
resolution, $10^{11}$ in time, and $10^{13}$ in particle number (to deal with
concentrations of Ca$^{2+}$). $\mathbf{These}$ $\mathbf{scales}$
$\mathbf{must}$ $\mathbf{be}$ $\mathbf{dealt}$ $\mathbf{with}$
$\mathbf{simultaneously}$ if the simulation is to deal with most biological
functions. We suggest a computational approach using explicit multiscale
analysis instead of implicit simulation of all scales. Successful computation
of ions concentrated in special places may be a significant step to
understanding the defining characteristics of biological and electrochemical
systems.


  Population dynamics reflects an underlying birth-death process, where the
rates associated with different events may depend on external environmental
conditions and on the population density. A whole family of simple and popular
deterministic models (like logistic growth) support a transcritical bifurcation
point between an extinction phase and an active phase. Here we provide a
comprehensive analysis of the phases of that system, taking into account both
the endogenous demographic noise (random birth and death events) and the effect
of environmental stochasticity that causes variations in birth and death rates.
Three phases are identified: in the inactive phase the mean time to extinction
$T$ is independent of the carrying capacity $N$, and scales logarithmically
with the initial population size. In the power-law phase $T \sim N^q$ and the
exponential phase $T \sim exp(\alpha N)$. All three phases and the transitions
between them are studied in detail. The breakdown of the continuum
approximation is identified inside the power-law phase, and the accompanied
changes in decline modes are analyzed. The applicability of the emerging
picture to the analysis of ecological timeseries and to the management of
conservation efforts is briefly discussed.


  Biochemical reactions with oscillatory behavior play an essential role in
synthetic biology at the microscopic scale. Although a robust stability theory
for deterministic chemical oscillators in the macroscopic limit exists, the
dynamical stability of stochastic oscillators is an object of ongoing research.
The Chemical Master Equation along with kinetic Monte Carlo simulations
constitute the most accurate approach to modeling microscopic systems. However,
because of the challenges of solving the fully probabilistic model, most
studies in stability analysis have focused on the description of externally
disturbed oscillators. Here we apply the Maximum Entropy Principle as closure
criterion for moment equations of oscillatory networks and perform the
stability analysis of the internally disturbed Brusselator network.
Particularly, we discuss the effects of kinetic and size parameters on the
dynamics of this stochastic oscillatory system with intrinsic noise. Our
numerical experiments reveal that changes in kinetic parameters lead to
phenomenological and dynamical Hopf bifurcations, while reduced system sizes in
the oscillatory region can reverse the stochastic Hopf dynamical bifurcations
at the ensemble level. This is a unique feature of the stochastic dynamics of
oscillatory systems, with an unknown parallel in the macroscopic limit.


  Single kinesin molecular motors can processively move along a microtubule
(MT) a few micrometers on average before dissociating. However, cellular length
scales over which transport occurs are several hundred microns and more. Why
seemingly unreliable motors are used to transport cellular cargo remains poorly
understood. We propose a new theory for how low processivity, the average
length of a single bout of directed motion, can enhance cellular transport when
motors and cargoes must first diffusively self assemble into complexes. We
employ stochastic modeling to determine the effect of processivity on overall
cargo transport flux. We show that, under a wide range of physiologically
relevant conditions, possessing "infinite" processivity does not maximize flux
along MTs. Rather, we find that low processivity i.e., weak binding of motors
to MTs, is optimal. These results shed light on the relationship between
processivity and transport efficiency and offer a new theory for the
physiological benefits of low motor processivity.


  Felsenstein's classical model for Gaussian distributions on a phylogenetic
tree is shown to be a toric variety in the space of concentration matrices. We
present an exact semialgebraic characterization of this model, and we
demonstrate how the toric structure leads to exact methods for maximum
likelihood estimation. Our results also give new insights into the geometry of
ultrametric matrices.


  Enzymes speed up biochemical reactions at the core of life by as much as 15
orders of magnitude. Yet, despite considerable advances, the fine dynamical
determinants at the microscopic level of their catalytic proficiency are still
elusive. In this work, we use a powerful mathematical approach to show that
rate-promoting vibrations in the picosecond range, specifically encoded in the
3D protein structure, are localized vibrations optimally coupled to the
chemical reaction coordinates at the active site. The universality of these
features is demonstrated on a pool of more than 900 enzyme structures,
comprising a total of more than 10,000 experimentally annotated catalytic
sites. Our theory provides a natural microscopic rationale for the known subtle
structural compactness of active sites in enzymes.


  We predict disease-genes relations on the Human Interactome network using a
methodology that jointly learns functional and connectivity patterns
surrounding proteins. Contrary to other data structures, the Interactome is
characterized by high incompleteness and absence of explicit negative
knowledge, which makes predictive tasks particularly challenging. To exploit at
best latent information in the network, we propose an extended version of
random walks, named Random Watcher-Walker ($RW^2$), which is able to learn rich
representations of disease genes (or gene products) features. Our method
successfully compares with the best known system for disease gene prediction,
and other state-of-the-art graph-based methods. We perform sensitivity analysis
and apply perturbations to ensure robustness. In contrast with previous
studies, our results demonstrate that connectivity alone is not sufficient to
classify disease-related genes.


  Estimating the abundances of all $k$-mers in a set of biological sequences is
a fundamental and challenging problem with many applications in biological
analysis. While several methods have been designed for the exact or approximate
solution of this problem, they all require to process the entire dataset, that
can be extremely expensive for high-throughput sequencing datasets. While in
some applications it is crucial to estimate all $k$-mers and their abundances,
in other situations reporting only frequent $k$-mers, that appear with
relatively high frequency in a dataset, may suffice. This is the case, for
example, in the computation of $k$-mers' abundance-based distances among
datasets of reads, commonly used in metagenomic analyses. In this work, we
develop, analyze, and test, a sampling-based approach, called SAKEIMA, to
approximate the frequent $k$-mers and their frequencies in a high-throughput
sequencing dataset while providing rigorous guarantees on the quality of the
approximation. SAKEIMA employs an advanced sampling scheme and we show how the
characterization of the VC dimension, a core concept from statistical learning
theory, of a properly defined set of functions leads to practical bounds on the
sample size required for a rigorous approximation. Our experimental evaluation
shows that SAKEIMA allows to rigorously approximate frequent $k$-mers by
processing only a fraction of a dataset and that the frequencies estimated by
SAKEIMA lead to accurate estimates of $k$-mer based distances between
high-throughput sequencing datasets. Overall, SAKEIMA is an efficient and
rigorous tool to estimate $k$-mers abundances providing significant speed-ups
in the analysis of large sequencing datasets.


  Seagrass meadows are among the most important coastal ecosystems, in terms of
both spatial extent and ecosystem services, but they are also declining
worldwide. Understanding the drivers of seagrass meadow dynamics is essential
for designing sound management, conservation, and restoration strategies.
However, the poor knowledge of the effect of clonality on the population
genetics of natural populations severely limits our understanding of the
dynamics and connectivity of meadows. Recent modelling approaches have
described the expected distributions of genotypic and genetic descriptors under
increasing clonal rates. Here, in light of these recent theoretical
developments, we revisited population genetic data for 165 meadows of four
seagrass species. Contrasting shoot life span and rhizome turnover led to the
prediction that the influence of asexual reproduction will increase along a
gradient from Zostera noltii to Zostera marina, Cymodocea nodosa and Posidonia
oceanica, which should be reflected by an increasing departure from
Hardy-Weinberg equilibrium (Fis) and decreasing genotypic diversity (R). This
meta-analysis provides a nested validation of this hypothesis at both the
species and meadow scales through a significant relationship between Fis and R
within each species. By empirically demonstrating the theoretical expectations
derived from recent modelling approaches, this work calls for the use of
Hardy-Weinberg equilibrium (Fis) rather than the strongly sampling-sensitive
genotypic index R to assess the importance of clonal reproduction, at least
when the impact of selfing on Fis can be neglected. The results also emphasize
the need to revise our appraisal of the extent of clonality and its influence
on the dynamics, connectivity and evolutionary trajectory of partial asexuals
in general, including in seagrass meadows, to develop the most accurate
management strategies.


  Heterogeneity is ubiquitous in stem cells (SC), cancer cells (CS), and cancer
stem cells (CSC). SC and CSC heterogeneity is manifested as diverse
sub-populations with self-renewing and unique regeneration capacity. Moreover,
the CSC progeny possesses multiple plasticity and cancerous characteristics.
Many studies have demonstrated that cancer heterogeneity is one of the greatest
obstacle for therapy. This leads to the incomplete anti-cancer therapies and
transitory efficacy. Furthermore, numerous micro-metastasis leads to the wide
spread of the tumor cells across the body which is the beginning of metastasis.
The epigenetic processes (DNA methylation or histone remodification etc.) can
provide a source for certain heterogeneity. In this study, we develop a
mathematical model to quantify the heterogeneity of SC, CSC and cancer taking
both genetic and epigenetic effects into consideration. We uncovered the roles
and physical mechanisms of heterogeneity from the three aspects (SC, CSC and
cancer). In the adiabatic regime (relatively fast regulatory binding and
effective coupling among genes), seven native states (SC, CSC, Cancer,
Premalignant, Normal, Lesion and Hyperplasia) emerge. In non-adiabatic regime
(relatively slow regulatory binding and effective weak coupling among genes),
multiple meta-stable SC, CS, CSC and differentiated states emerged which can
explain the origin of heterogeneity. In other words, the slow regulatory
binding mimicking the epigenetics can give rise to heterogeneity. Elucidating
the origin of heterogeneity and dynamical interrelationship between
intra-tumoral cells has clear clinical significance in helping to understand
the cellular basis of treatment response, therapeutic resistance, and tumor
relapse.


  Growing mixtures of annual arable crop species or genotypes is a promising
way to improve crop production without increasing agricultural inputs. To
design optimal crop mixtures, choices of species, genotypes, sowing proportion,
plant arrangement, and sowing date need to be made but field experiments alone
are not sufficient to explore such a large range of factors. Crop modeling
allows to study, understand and ultimately design cropping systems and is an
established method for sole crops. Recently, modeling started to be applied to
annual crop mixtures as well. Here, we review to what extent crop simulation
models and individual-based models are suitable to capture and predict the
specificities of annual crop mixtures. We argued that: 1) The crop mixture
spatio-temporal heterogeneity (influencing the occurrence of ecological
processes) determines the choice of the modeling approach (plant or crop
centered). 2) Only few crop models (adapted from sole crop models) and
individual-based models currently exist to simulate annual crop mixtures. 3)
Crop models are mainly used to address issues related to crop mixtures
management and to the integration of crop mixtures into larger scales such as
the rotation, whereas individual-based models are mainly used to identify plant
traits involved in crop mixture performance and to quantify the relative
contribution of the different ecological processes (niche complementarity,
facilitation, competition, plasticity) to crop mixture functioning. This review
highlights that modeling of annual crop mixtures is in its infancy and gives to
model users some important keys to choose the model based on the questions they
want to answer, with awareness of the strengths and weaknesses of each of the
modeling approaches.


  Policies to mitigate climate change and biodiversity loss often assume that
protecting carbon-rich forests provides co-benefits in terms of biodiversity,
due to the spatial congruence of carbon stocks and biodiversity at
biogeographic scales. However, it remains unclear whether this holds at the
scales relevant for management, with particularly large knowledge gaps for
temperate forests and for taxa other than trees. We built a comprehensive
dataset of Central European temperate forest structure and multi-taxonomic
diversity (beetles, birds, bryophytes, fungi, lichens, and plants) across 352
plots. We used Boosted Regression Trees to assess the relationship between
above-ground live carbon stocks and (a) taxon-specific richness, (b) a unified
multidiversity index. We used Threshold Indicator Taxa ANalysis to explore
individual species' responses to changing above-ground carbon stocks and to
detect change-points in species composition along the carbon-stock gradient.
Our results reveal an overall weak and highly variable relationship between
richness and carbon stock at the stand scale, both for individual taxonomic
groups and for multidiversity. Similarly, the proportion of win-win and
trade-off species (i.e. species favored or disadvantaged by increasing carbon
stock, respectively) varied substantially across taxa. Win-win species
gradually replaced trade-off species with increasing carbon, without clear
thresholds along the above-ground carbon gradient, suggesting that
community-level surrogates (e.g. richness) might fail to detect critical
changes in biodiversity. Collectively, our analyses highlight that leveraging
co-benefits between carbon and biodiversity in temperate forest may require
stand-scale management that prioritizes either biodiversity or carbon-in order
to maximize co-benefits at broader scales. Importantly, this contrasts with
tropical forests, where climate [...]


  Ediacaran fossils are now largely known in different parts of the world.
However, some countries are poorly documented on these remains of a still
enigmatic life. Thus, rare fossils from the Neoproterozoic Histria Formation of
central Dobrogea (Romania) have been reported. Two specimens with discoid
imprints are described here in detail and assigned to the typical Ediacaran
species Beltanelliformis brunsae Menner in Keller et al., 1974. This
paleontological development confirms both the large geographical distribution
of this species and the Ediacaran age of the Histria Formation. KEY WORDS Flat
discoid imprints, Ediacaran biota, Precambrian, Dobrogea, Romania.
R{\'E}SUM{\'E}


  Antifragility is a property from which systems are able to resist stress and
furthermore benefit from it. Even though antifragile dynamics is found in
various real-world complex systems where multiple subsystems interact with each
other, the attribute has not been quantitatively explored yet in those complex
systems which can be regarded as multilayer networks. Here we study how the
multilayer structure affects the antifragility of the whole system. By
comparing single-layer and multilayer Boolean networks based on our recently
proposed antifragility measure, we found that the multilayer structure
facilitated the production of antifragile systems. Our measure and findings
will be useful for various applications such as exploring properties of
biological systems with multilayer structures and creating more antifragile
engineered systems.


  We study the multiscale structure of the Jain-Krishna adaptive network model.
This model describes the co-evolution of a set of continuous-time autocatalytic
ordinary differential equations and its underlying discrete-time graph
structure. The graph dynamics is governed by deletion of vertices with
asymptotically weak concentrations of prevalence and then re-insertion of
vertices with new random connections. In this work we prove several results
about convergence of the continuous-time dynamics to equilibrium points.
Furthermore, we motivate via formal asymptotic calculations several conjectures
regarding the discrete-time graph updates. In summary, our results clearly show
that there are several time scales in the problem depending upon system
parameters, and that analysis can be carried out in certain singular limits.
This shows that for the Jain-Krishna model, and potentially many other adaptive
network models, a mixture of deterministic and/or stochastic multiscale methods
is a good approach to work towards a rigorous mathematical analysis.


  Some bacteria and archaea possess an immune system, based on the CRISPR-Cas
mechanism, that confers adaptive immunity against phage. In such species,
individual bacteria maintain a "cassette" of viral DNA elements called spacers
as a memory of past infections. The typical cassette contains a few dozen
spacers. Given that bacteria can have very large genomes, and since having more
spacers should confer a better memory, it is puzzling that so little genetic
space would be devoted by bacteria to their adaptive immune system. Here, we
identify a fundamental trade-off between the size of the bacterial immune
repertoire and effectiveness of response to a given threat, and show how this
tradeoff imposes a limit on the optimal size of the CRISPR cassette.


  Public health surveillance systems often fail to detect emerging infectious
diseases, particularly in resource limited settings. By integrating relevant
clinical and internet-source data, we can close critical gaps in coverage and
accelerate outbreak detection. Here, we present a multivariate algorithm that
uses freely available online data to provide early warning of emerging
influenza epidemics in the US. We evaluated 240 candidate predictors and found
that the most predictive combination does \textit{not} include surveillance or
electronic health records data, but instead consists of eight Google search and
Wikipedia pageview time series reflecting changing levels of interest in
influenza-related topics. In cross validation on 2010-2016 data, this algorithm
sounds alarms an average of 16.4 weeks prior to influenza activity reaching the
Center for Disease Control and Prevention (CDC) threshold for declaring the
start of the season. In an out-of-sample test on data from the rapidly-emerging
fall wave of the 2009 H1N1 pandemic, it recognized the threat five weeks in
advance of this surveillance threshold. Simpler algorithms, including fixed
week-of-the-year triggers, lag the optimized alarms by only a few weeks when
detecting seasonal influenza, but fail to provide early warning in the 2009
pandemic scenario. This demonstrates a robust method for designing next
generation outbreak detection algorithms. By combining scan statistics with
machine learning, it identifies tractable combinations of data sources (from
among thousands of candidates) that can provide early warning of emerging
infectious disease threats worldwide.


  Genome-wide association studies (GWAS) have been widely used to examine the
association between single nucleotide polymorphisms (SNPs) and complex traits,
where both the sample size n and the number of SNPs p can be very large.
Recently, cross-trait polygenic risk score (PRS) method has gained extremely
popular for assessing genetic correlation of complex traits based on GWAS
summary statistics (e.g., SNP effect size). However, empirical evidence has
shown a common bias phenomenon that even highly significant cross-trait PRS can
only account for a very small amount of genetic variance (R^2 often <1%). The
aim of this paper is to develop a novel and powerful method to address the bias
phenomenon of cross-trait PRS. We theoretically show that the estimated genetic
correlation is asymptotically biased towards zero when complex traits are
highly polygenic/omnigenic. When all p SNPs are used to construct PRS, we show
that the asymptotic bias of PRS estimator is independent of the unknown number
of causal SNPs m. We propose a consistent PRS estimator to correct such
asymptotic bias. We also develop a novel estimator of genetic correlation which
is solely based on two sets of GWAS summary statistics. In addition, we
investigate whether or not SNP screening by GWAS p-values can lead to improved
estimation and show the effect of overlapping samples among GWAS. Our results
may help demystify and tackle the puzzling "missing genetic overlap" phenomenon
of cross-trait PRS for dissecting the genetic similarity of closely related
heritable traits. We illustrate the finite sample performance of our
bias-corrected PRS estimator by using both numerical experiments and the UK
Biobank data, in which we assess the genetic correlation between brain white
matter tracts and neuropsychiatric disorders.


  Multiple lines of evidence point to one or more moderately nearby supernovae,
with the strongest signal ~2.6 Ma. We build on previous work to argue for the
likelihood of cosmic ray ionization of the atmosphere and electron cascades
leading to more frequent lightning, and therefore an increase in nitrate
deposition and in wildfires. The potential exists for a large increase in the
pre-human nitrate flux onto the surface, which has previously been argued to
lead to CO2 drawdown and cooling of the climate. Evidence for increased
wildfires exists in an increase in soot and carbon deposits over the relevant
period. The wildfires would have contributed to the transition from forest to
savanna in northeast Africa, long argued to have been a factor in the evolution
of hominin bipedalism.


  In recent studies, non-coding protein RNAs have been identified as microRNA
that can be used as biomarkers for early diagnosis and treatment of cancer,
that decrease mortality in cancer. A microRNA may target hundreds or thousands
of genes and a gene may regulate several microRNAs, so determining which
microRNA is associated with which cancer is a big challenge. Many computational
methods have been performed to detect micoRNAs association with cancer, but
more effort is needed with higher accuracy. Increasing research has shown that
relationship between microRNAs and TFs play a significant role in the diagnosis
of cancer. Therefore, we developed a new computational framework (CAMIRADA) to
identify cancer-related microRNAs based on the relationship between microRNAs
and disease genes (DG) in the protein network, the functional relationships
between microRNAs and Transcription Factors (TF) on the co-expression network,
and the relationship between microRNAs and the Differential Expression Gene
(DEG) on co-expression network. The CAMIRADA was applied to assess breast
cancer data from two HMDD and miR2Disease databases. In this study, the AUC for
the 65 microRNAs of the top of the list was 0.95, which was more accurate than
the similar methods used to detect microRNAs associated with the cancer artery.


  This review focuses on microbiomes associated with residues within the
context of other microbial habitats in cereal-producing agroecosystems such as
phyllosphere or rhizosphere. We connect residue microbiome with the survival of
residue-borne fungal plant pathogens, thus combining knowledge in microbial
ecology and epidemiology, two disciplines still not sufficiently connected. We
provide an overview of the impact of residues on cereal disease epidemics and
how dynamic interactions between microbial communities of non-buried residues
during their degradation, along with soil and multitude of abiotic factors, can
contribute to innovative disease management strategies, including
next-generation microbiome-based biocontrol strategies. Starting from the
classical but still relevant view of crop residues as a source of pathogen
inoculum, we first consider possibilities for limiting the amount of residues
on the soil surface to reduce the pathogen pressure. We then describe residues
as a transient half-plant/half-soil compartment constituting a key fully
fledged microbial ecosystem: in other words, an ecotone which deserves special
attention. We focus on microbial communities, the changes in these communities
over time and the factors influencing them. Finally, we discuss how the
interactions between the microbial communities and the pathogens present on
residues could be used: identification of keystone taxa and beneficial
assemblages, then preservation of these taxa by adapted agronomic practices or
development of synthetic communities, rather than the introduction of a single
exogenous biocontrol species designed as a 'treatment product'.


  We propose a fractional order model for HIV/AIDS transmission. Local and
uniform stability of the fractional order model is studied. The theoretical
results are illustrated through numerical simulations.


  This work concerns the question of how two important dynamical properties,
oscillations and bistability, emerge in an important biological signaling
network. Specifically, we consider a model for dual-site phosphorylation and
dephosphorylation of extracellular signal-regulated kinase (ERK). We prove that
oscillations persist even as the model is greatly simplified (reactions are
made irreversible and intermediates are removed). Bistability, however, is much
less robust -- this property is lost when intermediates are removed or even
when all reactions are made irreversible. Moreover, bistability is
characterized by the presence of two reversible, catalytic reactions: as other
reactions are made irreversible, bistability persists as long as one or both of
the specified reactions is preserved. Finally, we investigate the maximum
number of steady states, aided by a network's "mixed volume" (a concept from
convex geometry). Taken together, our results shed light on the question of how
oscillations and bistability emerge from a limiting network of the ERK network
-- namely, the fully processive dual-site network -- which is known to be
globally stable and therefore lack both oscillations and bistability. Our
proofs are enabled by a Hopf bifurcation criterion due to Yang, analyses of
Newton polytopes arising from Hurwitz determinants, and recent
characterizations of multistationarity for networks having a steady-state
parametrization.


  Resource sharing outside the kinship bonds is rare. Besides humans, it occurs
in chimpanzee, wild dogs and hyenas as well as in vampire bats. Resource
sharing is an instance of animal cooperation, where an animal gives away part
of the resources that it owns for the benefit of a recipient. Taking
inspiration from blood-sharing in vampire bats, here show the emergence of
generosity in a Markov game, which couples the resource sharing between two
players with the gathering task of that resource. At variance with the
classical evolutionary models for cooperation, the optimal strategies of this
game can be potentially learned by animals during their life-time. The players
act greedily, that is, they try to individually maximize only their personal
income. Nonetheless, the analytical solution of the model shows that three non
trivial optimal behaviours emerge depending on conditions. Besides the obvious
case when players are selfish in their choice of resource division, there are
conditions under which both players are generous. Moreover, we also found a
range of situations in which one selfish player exploits another generous
individual, for the satisfaction of both players. Our results show that
resource sharing is favoured by three factors: a long time horizon over which
the players try to optimize their own game, the similarity among players in
their ability of performing the resource-gathering task, as well as by the
availability of resources in the environment. These concurrent requirements
lead to identify necessary conditions for the emergence of generosity.


  Stock management should be guided by assessment models that, among others,
need to be fed by reliable data of catch and effort. However, precise data are
difficult to obtain in heterogeneous fisheries. Specifically, small scale,
multi gear, multispecies fisheries are dynamic systems where fishers may lively
change fishing strategy conditioned by multiple drivers. Provided that some
stocks can be shared by several metiers, a precise categorization of metiers
should be the first step toward metier specific estimates of catch and effort,
which in turn would allow a better understanding of the system dynamics. Here
we propose an approach for predicting the metier of any given fishing trip from
its landing records. This approach combines the knowledge of expert fishers
with the existing sales register of landings in Mallorca. It successfully
predicts metiers for all the 162815 small scale fishery fishing trips from
Mallorca between 2004 and 2015. The largest effort is invested in the metiers
Cuttlefish Fish and Spiny lobster, landings peak for Cuttlefish Fish and
Dolphinfish and revenues for Spiny lobster and Dolphinfish. Metier predictions
also allowed us to describe the temporal trends experienced by each metier and
to characterize the species that are specific to each metier. Seasonal
variability is by far more relevant than between year variability, which
confirms that at least some fishers are adopting a rotation cycle of metiers
along the year. Effort, landings and gross revenues decreased in the last 12
years. The approach proposed is also applicable to any other fishery for which
the metier for a fishing trip sample is known, but relying on fishers expertise
points more directly to fishers intention. Thus, metier predictions produced
with the proposed approach are closer to the actual uses of fishers, providing
better grounds for an improved management.


  In this paper, we systematically review weighted persistent homology (WPH)
models and their applications in biomolecular data analysis. Essentially, the
weight value, which reflects physical, chemical and biological properties, can
be assigned to vertices (atom centers), edges (bonds), or higher order
simplexes (cluster of atoms), depending on the biomolecular structure,
function, and dynamics properties. Further, we propose the first localized
weighted persistent homology (LWPH). Inspired by the great success of element
specific persistent homology (ESPH), we do not treat biomolecules as an
inseparable system like all previous weighted models, instead we decompose them
into a series of local domains, which may be overlapped with each other. The
general persistent homology or weighted persistent homology analysis is then
applied on each of these local domains. In this way, functional properties,
that are embedded in local structures, can be revealed. Our model has been
applied to systematically studying DNA structures. It has been found that our
LWPH based features can be used to successfully discriminate the A-, B-, and
Z-types of DNA. More importantly, our LWPH based PCA model can identify two
configurational states of DNA structure in ion liquid environment, which can be
revealed only by the complicated helical coordinate system. The great
consistence with the helical-coordinate model demonstrates that our model
captures local structure variations so well that it is comparable with
geometric models. Moreover, geometric measurements are usually defined in very
local regions. For instance, the helical-coordinate system is limited to one or
two basepairs. However, our LWPH can quantitatively characterize structure
information in local regions or domains with arbitrary sizes and shapes, where
traditional geometrical measurements fail.


  We devise a method to detect and estimate forces in a heterogeneous
environment based on experimentally recorded stochastic trajectories. In
particular, we focus on systems modeled by the heterogeneous overdamped
Langevin equation. Here, the observed drift includes a "spurious" force term
when the diffusivity varies in space. We show how Bayesian inference can be
leveraged to reliably infer forces by taking into account such spurious forces
of unknown amplitude as well as experimental sources of error. The method is
based on marginalizing the force posterior over all possible spurious force
contributions. The approach is combined with a Bayes factor statistical test
for the presence of forces. The performance of our method is investigated
analytically, numerically and tested on experimental data sets. The main
results are obtained in a closed form allowing for direct exploration of their
properties and fast computation. The method is incorporated into TRamWAy, an
open-source software platform for automated analysis of biomolecule
trajectories.


  Gene regulatory networks must relay information from extracellular signals to
downstream genes in an efficient, timely and coherent manner. Many complex
functional tasks such as the immune response require system-wide broadcasting
of information not to one but to many genes carrying out distinct functions
whose dynamical binding and unbinding characteristics are widely distributed.
In such broadcasting networks, the intended target sites are also often dwarfed
in number by the even more numerous non-functional binding sites. Taking the
genetic regulatory network of NF\k{appa}B as an exemplary system we explore the
impact of having numerous distributed sites on the stochastic dynamics of
oscillatory broadcasting genetic networks pointing out how resonances in
binding cycles control the network's specificity and performance. We also show
that active kinetic regulation of binding and unbinding through molecular
stripping of DNA bound transcription factors can lead to a higher coherence of
gene co-expression and synchronous clearance.


  Modern biological tools have made it possible to unequivocally demonstrate
the deep relationship among species in terms of genes and basic molecular
mechanisms. In addition, results from genetic, physical and physiological
approaches applied to individual model cell systems have led to the mapping of
large networks of macromolecular interactions with similar general properties.
Although gene, mechanism and network structure similarities among species tend
to suggest the existence of important constraints applied to organisms, it is
surprising that the elucidation of the precise general mechanisms by which
natural selection could have operated appears to have taken a back seat in
most, if not all, studies concerning the evolutionary development of cell
macromolecular networks. Herein, a possible explanation is presented for how
cells could have evolved into the sophisticated and integrated organisms we
know today, using natural selection. Based on the new concept of gene toxicity,
essentially applied to unicellular organisms with their rapid multiplication
rate, it is proposed that natural selection likely exerts its effects in two
opposing and complementary directions to integrate a new gene function into the
host system; a simple and basic qualitative model for the buildup of cell
macromolecular networks is suggested. Finally, as a radical consequence, it is
proposed that complex life evolution has proceeded through a long alternation
of unicellular and multicellular states.


  Maximum likelihood estimators are used extensively to estimate unknown
parameters of stochastic trait evolution models on phylogenetic trees. Although
the MLE has been proven to converge to the true value in the independent-sample
case, we cannot appeal to this result because trait values of different species
are correlated due to shared evolutionary history. In this paper, we consider a
$2$-state symmetric model for a single binary trait and investigate the
theoretical properties of the MLE for the transition rate in the large-tree
limit. Here, the large-tree limit is a theoretical scenario where the number of
taxa increases to infinity and we can observe the trait values for all species.
Specifically, we prove that the MLE converges to the true value under some
regularity conditions. These conditions ensure that the tree shape is not too
irregular, and holds for many practical scenarios such as trees with bounded
edges, trees generated from the Yule (pure birth) process, and trees generated
from the coalescent point process. Our result also provides an upper bound for
the distance between the MLE and the true value.


  Protein molecules can be approximated by discrete polygonal chains of amino
acids. Standard topological tools can be applied to the smoothening of the
polygons to introduce a topological classification of proteins, for example,
using the self-linking number of the corresponding framed curves. In this paper
we add new details to the standard classification. Known definitions of the
self-linking number apply to non-singular framings: for example, the Frenet
framing cannot be used if the curve has inflection points. Meanwhile in the
discrete proteins the special points are naturally resolved. Consequently, a
separate integer topological characteristics can be introduced, which takes
into account the intrinsic features of the special points. For large number of
proteins we compute integer topological indices associated with the
singularities of the Frenet framing. We show how a version of the Calugareanu's
theorem is satisfied for the associated self-linking number of a discrete
curve. Since the singularities of the Frenet framing correspond to the
structural motifs of proteins, we propose topological indices as a technical
tool for the description of the folding dynamics of proteins.


  Several technological limitations of traditional silicon based computing are
leading towards the paradigm shift, from silicon to carbon, in computational
world. Among the unconventional modes of computing evolved in past several
decades, DNA computing has been considered to be quite promising in solving
computational and reasoning problems by using DNA strands. Along with the
sequential operations, the huge parallelism of DNA computing methodologies
engaging numerous numbers of DNA strands induce the consideration of concurrent
high-level formalisms. In this paper we have reviewed the algebraic explanation
of concurrent DNA processes using DNA strand algebra, process calculus and DNA
strand graph. We have demonstrated the application of syntax and semantics of
the illustrated methodologies in the domains of reasoning and theorem proving
with resolution refutation. Finally, we have presented DNA cryptography as one
of the prominent areas for the future scope of research work where DNA strand
algebra can be used as formal modelling tool to authenticate the security,
logic and reasoning of the existing protocols.


  The evolution of dispersal rate is studied with a model of several local
populations linked by dispersal. Three dispersal strategies are considered
where all, half, or none of the offspring disperse. The spatial scale (number
of patches) and the temporal scale (probability of local extinction) of the
environment are critical in determining the selective advantage of the
different dispersal strategies. The results from the simulations suggest that
an interaction between group selection and individual selection results in a
different outcome in relation to the spatial and temporal scales of the
environment. Such an interaction is able to maintain a polymorphism in
dispersal strategies. The maintenance of this polymorphism is also
scale-dependent. This study suggests a mechanism for the short-term evolution
of dispersal, and provides a testable prediction of this hypothesis, namely
that loss of dispersal abilities should be more frequent in spatially more
continuous environments, or in temporally more stable environments.


  We introduce an agent-based model describing a
susceptible-infectious-susceptible (SIS) system of humans and mosquitoes to
predict malaria epidemiological scenarios in realistic biological conditions.
Emphasis is given to the transition from endemic behavior to eradication of
malaria transmission induced by combined drug therapies acting on both the
gametocytemia reduction and on the selective mosquito mortality during parasite
development in the mosquito. Our mathematical framework enables to uncover the
critical values of the parameters characterizing the effect of each drug
therapy. Moreover, our results provide quantitative evidence of what is
empirically known: interventions combining gametocytemia reduction through the
use of gametocidal drugs, with the selective action of ivermectin during
parasite development in the mosquito, may actively promote disease eradication
in the long run. In the agent model, the main properties of human-mosquito
interactions are implemented as parameters and the model is validated by
comparing simulations with real data of malaria incidence collected in the
endemic malaria region of Chimoio in Mozambique. Finally, we discuss our
findings in light of current drug administration strategies for malaria
prevention, that may interfere with human-to-mosquito transmission process.


  In solution, UV-vis spectroscopy is often used to investigate structural
changes in biomolecules (i.e., nucleic acids), owing to changes in the
environment of their chromophores (i.e., the nucleobases). Here we address
whether action spectroscopy could achieve the same for gas-phase ions, while
taking the advantage of additional spectrometric separation of complex
mixtures. We therefore systematically studied the action spectroscopy of
homo-base 6-mer DNA strands (dG6, dA6, dC6, dT6) and discuss the results in
light of gas-phase structures validated by ion mobility spectrometry and
infrared ion spectroscopy, of electron binding energies measured by
photoelectron spectroscopy, and of calculated electronic photo-absorption
spectra. When UV photons interact with oligonucleotide polyanions, two main
actions may take place: (1) fragmentation and (2) electron detachment. The
action spectra reconstructed from fragmentation follow the absorption spectra
well, and result from multiple cycles of absorption and internal conversion.
The action spectra reconstructed from the electron photodetachment (ePD)


  In an iterated game between two players, there is much interest in
characterizing the set of feasible payoffs for both players when one player
uses a fixed strategy and the other player is free to switch. Such
characterizations have led to extortionists, equalizers, partners, and rivals.
Most of those studies use memory-one strategies, which specify the
probabilities to take actions depending on the outcome of the previous round.
Here, we consider "reactive learning strategies," which gradually modify their
propensity to take certain actions based on past actions of the opponent. Every
linear reactive learning strategy, $\mathbf{p}^{\ast}$, corresponds to a memory
one-strategy, $\mathbf{p}$, and vice versa. We prove that for evaluating the
region of feasible payoffs against a memory-one strategy,
$\mathcal{C}\left(\mathbf{p}\right)$, we need to check its performance against
at most $11$ other strategies. Thus, $\mathcal{C}\left(\mathbf{p}\right)$ is
the convex hull in $\mathbb{R}^{2}$ of at most $11$ points. Furthermore, if
$\mathbf{p}$ is a memory-one strategy, with feasible payoff region
$\mathcal{C}\left(\mathbf{p}\right)$, and $\mathbf{p}^{\ast}$ is the
corresponding reactive learning strategy, with feasible payoff region
$\mathcal{C}\left(\mathbf{p}^{\ast}\right)$, then
$\mathcal{C}\left(\mathbf{p}^{\ast}\right)$ is a subset of
$\mathcal{C}\left(\mathbf{p}\right)$. Reactive learning strategies are
therefore powerful tools in restricting the outcomes of iterated games.


  Classic concepts of genetic (gene) diversity (heterozygosity) such as Nei
(1973: PNAS) and Nei and Li (1979: PNAS) nucleotide diversity were defined
within the context of populations. Although variations are often measured in
population context, the basic carriers of variation are individuals. Hence,
measuring variations such as SNP of individual against a reference genome,
which has been ignored currently, is certainly of its own right. Indeed,
similar practice has been a tradition in ecology, where the basic framework of
diversity measure is individual community sample. We propose to use
Renyi-entropy-derived Hill numbers to define SNP (single nucleotide
polymorphism) diversity (including alpha-, beta-, and gamma-diversities) and
similarity profiles. Hill numbers are derived from Renyi entropy, of which
Shannon entropy is a special case and which have found widely applications
including measuring the quantum information entanglement, wealth distribution
in economics and ecological diversity. The newly proposed SNP diversity not
only complements the existing genetic diversity concepts by offering
individual-level metrics, but also offers building blocks for comparative
genetic analysis at higher levels. The profile concept also helps to resolve a
dilemma in measuring diversity: the choice from various diversity indexes,
because diversity profile unifies some of the most commonly used indexes (as
special cases) with different diversity orders (along the rareness-commonness
spectrum of gene mutations). Finally, the profiles can be estimated with
rarefaction approach, which may help to relieve some effect of insufficient
sequencing coverage.


  Orthology is a central concept in evolutionary and comparative genomics, used
to relate corresponding genes in different species. In particular, orthologs
are needed to infer species trees. In this chapter, we introduce the
fundamental concepts of orthology relationships and orthologous groups,
including some non-trivial (and thus commonly misunderstood) implications.
Next, we review some of the main methods and resources used to identify
orthologs. The final part of the chapter discusses the impact of orthology
methods on species phylogeny inference, drawing lessons from several recent
comparative studies.


  Many bacterial species produce toxins that inhibit their competitors. We
model this phenomenon by extending classic two-species Lotka-Volterra
competition in one spatial dimension to incorporate toxin production by one
species. Considering solutions comprising two adjacent single-species colonies,
we show how the toxin inhibits the susceptible species near the interface
between the two colonies. Moreover, a sufficiently effective toxin inhibits the
susceptible species to such a degree that an `inhibition zone' is formed
separating the two colonies. In the special case of truly non-motile bacteria,
i.e. with zero bacterial diffusivity, we derive analytical expressions
describing the bacterial distributions and size of the inhibition zone. In the
more general case of weakly motile bacteria, i.e. small bacterial diffusivity,
these two-colony solutions become travelling waves. We employ numerical methods
to show that the wavespeed is dependent upon both interspecific competition and
toxin strength; precisely which colony expands at the expense of the other
depends upon the choice of parameter values. In particular, a sufficiently
effective toxin allows the producer to expand at the expense of the
susceptible, with a wavespeed magnitude that is bounded above as the toxin
strength increases. This asymptotic wavespeed is independent of interspecific
competition and due to the formation of the inhibition zone; when the colonies
are thus separated, there is no longer direct competition between the two
species and the producer can invade effectively unimpeded by its competitor. We
note that the minimum toxin strength required to produce an inhibition zone
increases rapidly with increasing bacterial diffusivity, suggesting that even
moderately motile bacteria must produce very strong toxins if they are to
benefit in this way.


  Parasitological examinations of White-tailed Lapwing Vanellus leucurus
(Charadriiformes: Charadriidae) collected from Dadu district (Sindh Province,
Southern Pakistan) were carried out at the Department of Zoology, University of
Sindh, Jamshoro. A total of 64 trematodes belonging to the genus
Echinoparyphium [1], were collected from the small intestine of two individual
hosts. On the basis of the number and varying size of collar spines, number and
size of ventral spines, body shape, arrangement of genital organs, distribution
of vitellaria and other diagnostic characters, these trematodes were identified
as Echinoparyphium recurvatum [2]. Our study provides the first Pakistan record
of the trematode E. recurvatum with also the first record of V. leucurus as
avian definitive host.


  The utility of DNA sequence substrings (k-mers) in alignment-free
phylogenetic classification, including that of bacteria and viruses, is
increasingly recognized. However, its biological basis eludes many twenty-first
century practitioners. A path from the nineteenth century recognition of the
informational basis of heredity to the modern era can be discerned. Crick's DNA
"unpairing postulate" predicted that recombinational pairing of homologous DNAs
during meiosis would be mediated by short k-mers in the loops of stem-loop
structures extruded from classical duplex helices. The complementary "kissing"
duplex loops - like tRNA anticodon-codon k-mer duplexes - would seed a more
extensive pairing that would then extend until limited by lack of homology or
other factors. Indeed, this became the principle behind alignment-based methods
that assessed similarity by degree of DNA-DNA reassociation in vitro. These are
now seen as less sensitive than alignment-free methods that are closely
consistent, both theoretically and mechanistically, with chromosomal
anti-recombination models for the initiation of divergence into new species.
The analytical power of k-mer differences supports the theses that evolutionary
advance sometimes serves the needs of nucleic acids (genomes) rather than
proteins (genes), and that such differences have often played a role in early
speciation events.


  Single Molecule Real-Time (SMRT) sequencing is a recent advancement of Next
Gen technology developed by Pacific Bio (PacBio). It comes with an explosion of
long and noisy reads demanding cutting edge research to get most out of it. To
deal with the high error probability of SMRT data, a novel contextual Locality
Sensitive Hashing (conLSH) based algorithm is proposed in this article, which
can effectively align the noisy SMRT reads to the reference genome. Here,
sequences are hashed together based not only on their closeness, but also on
similarity of context. The algorithm has $\mathcal{O}(n^{\rho+1})$ space
requirement, where $n$ is the number of sequences in the corpus and $\rho$ is a
constant. The indexing time and querying time are bounded by $\mathcal{O}(
\frac{n^{\rho+1} \cdot \ln n}{\ln \frac{1}{P_2}})$ and $\mathcal{O}(n^\rho)$
respectively, where $P_2 > 0$, is a probability value. This algorithm is
particularly useful for retrieving similar sequences, a widely used task in
biology. The proposed conLSH based aligner is compared with rHAT, popularly
used for aligning SMRT reads, and is found to comprehensively beat it in speed
as well as in memory requirements. In particular, it takes approximately
$24.2\%$ less processing time, while saving about $70.3\%$ in peak memory
requirement for H.sapiens PacBio dataset.


  Agricultural losses to pest represent an important challenge in a global
warming scenario. Intercropping is an alternative farming practice that
promotes pest control without the use of chemical pesticides. Here we develop a
mathematical model to study epidemic spreading and control in intercropped
agricultural fields as a sustainable pest management tool for agriculture. The
model combines the movement of aphids transmitting a virus in an agricultural
field, the spatial distribution of plants in the intercropped field, and the
presence of `trap crops' in an epidemiological Susceptible-Infected-Removed
(SIR) model. Using this model we study several intercropping arrangements
without and with trap crops and find a new intercropping arrangement that
improves significantly pest management in agricultural fields respect to the
commonly used intercrop systems.


  HIV virus replication is a time-related process that includes several stages.
Focusing on the core steps, RNA polymerase II transcripts in an early stage
pre-mRNA containing regulator proteins (i.e nef,tat,rev,vif,vpr,vpu), which are
completely spliced by the spliceosome complex (0.9kb and 1.8kb) and exported to
the ribosome for protein synthesis. These splicing and export processes are
regulated by tat protein, which binds on Trans-activation response (TAR)
element, and by rev protein, which binds to the Rev-responsive Element (RRE).
As long as these regulators are synthesized, splicing is progressively
inhibited (from 4.0kb to 9.0kb) and mRNAs are translated into structural and
enzymatic proteins (env, gag-pol). During this RNAPII scanning and splicing,
around 40 different multi-cystronic mRNA have been produced. Long-read
sequencing has been applied to the HIV-1 virus genome (type HXB2CG) with the
HIV.pro software, a fortran 90 code for simulating the virus replication cycle,
specially RNAPII transcription, exon/intron splicing and ribosome protein
synthesis, including the frameshift at gag/pol gene and the ribosome pause at
env gene. All HIV-1 virus proteins have been identified as far as other ORFs.
As observed, tat/rev protein regulators have different length depending on the
splicing cleavage site: tat protein varies from 224aa to a final state of 72aa,
whereas rev protein from 25aa to 27aa, with a maximum of 119aa. Furthermore,
several ORFs coding for small polypeptides sPEP (less than 10 amino acids) and
for other unidentified proteins have been localised with unknown functionality.
The detailed analysis of the HIV virus replication and the virus proteomics are
important for identifying which antigens are presented by macrophages to CD4
cells, for localizing reactive epitopes or for creating transfer vectors to
develop new HIV vaccines and effective therapies.


  We introduce and analyze two general dynamical models for unidirectional
movement of particles along a circular chain and an open chain of sites. The
models include a soft version of the simple exclusion principle, that is, as
the density in a site increases the effective entry rate into this site
decreases. This allows to model and study the evolution of "traffic jams" of
particles along the chain. A unique feature of these two new models is that
each site along the chain can have a different size.
  Although the models are nonlinear, they are amenable to rigorous asymptotic
analysis. In particular, we show that the dynamics always converges to a
steady-state, and that the steady-state densities along the chain and the
steady-state output flow rate from the chain can be derived from the spectral
properties of a suitable matrix, thus eliminating the need to numerically
simulate the dynamics until convergence. This spectral representation also
allows for powerful sensitivity analysis, i.e. understanding how a change in
one of the parameters in the models affects the steady-state.
  We show that the site sizes and the transition rates from site to site play
different roles in the dynamics, and that for the purpose of maximizing the
steady-state output (or production) rate the site sizes are more important than
the transition rates. We also show that the problem of finding parameter values
that maximize the production rate is tractable.
  We believe that the models introduced here can be applied to study various
natural and artificial processes including ribosome flow during mRNA
translation, the movement of molecular motors along filaments of the
cytoskeleton, pedestrian and vehicular traffic, evacuation dynamics, and more.


  Mathematical models of spatial population dynamics typically focus on the
interplay between dispersal events and birth/death processes. However, for many
animal communities, significant arrangement in space can occur on shorter
timescales, where births and deaths are negligible. This phenomenon is
particularly prevalent in populations of larger, vertebrate animals who often
reproduce only once per year or less. To understand spatial arrangements of
animal communities on such timescales, we use a class of diffusion-taxis
equations for modelling inter-population movement responses between $N \geq 2$
populations. These systems of equations incorporate the effect on animal
movement of both the current presence of other populations and the memory of
past presence encoded either in the environment or in the minds of animals. We
give general criteria for the spontaneous formation of both stationary and
oscillatory patterns, via linear pattern formation analysis. For $N=2$, we
classify completely the pattern formation properties using a combination of
linear analysis and non-linear energy functionals. In this case, the only
patterns that can occur asymptotically in time are stationary. However, for $N
\geq 3$, oscillatory patterns can occur asymptotically, giving rise to a
sequence of period-doubling bifurcations leading to patterns with no obvious
regularity, a hallmark of chaos. Our study highlights the importance of
understanding between-population animal movement for understanding spatial
species distributions, something that is typically ignored in species
distribution modelling, and so develops a new paradigm for spatial population
dynamics.