  Following the recent publication of our book on Exploring the Health State of
a Population by Dynamic Modeling Methods in The Springer Series on Demographic
Methods and Population Analysis (DOI 10.1007/978-3-319-65142-2) we provide this
brief presentation of the main findings and improvements regarding the Health
State of a Population. (See at: http://www.springer.com/gp/book/9783319651415).
Here the brief history of the Health State or Health Status curves for
individuals and populations is presented including the main references and
important figures along with an illustrated Poster (see Figure 13 and
http://www.smtda.net/demographics2018.html). Although the Survival Curve is
known as long as the life tables have introduced, the Health State Curve was
calculated after the introduction of the advanced stochastic theory of the
first exit time. The health state curve is illustrated in several graphs either
as a fit curve to data or produced after a large number of stochastic
realizations. The Health State, the Life Expectancy and the age at mean zero
health state are also estimated. Keywords: Health State and Survival Curves,
Health status of a population, First exit time stochastic theory, stochastic
simulations of health state, Age at Maximum Curvature, Healthy Life Expectancy
and HALE, Standard Deviation, Health State Curves, Maximum human lifespan and
other.


  Thanks to rapidly evolving sequencing techniques, the amount of genomic data
at our disposal is growing increasingly large. Determining the gene structure
is a fundamental requirement to effectively interpret gene function and
regulation. An important part in that determination process is the
identification of translation initiation sites. In this paper, we propose a
novel approach for automatic prediction of translation initiation sites,
leveraging convolutional neural networks that allow for automatic feature
extraction. Our experimental results demonstrate that we are able to improve
the state-of-the-art approaches with a decrease of 75.2% in false positive rate
and with a decrease of 24.5% in error rate on chosen datasets. Furthermore, an
in-depth analysis of the decision-making process used by our predictive model
shows that our neural network implicitly learns biologically relevant features
from scratch, without any prior knowledge about the problem at hand, such as
the Kozak consensus sequence, the influence of stop and start codons in the
sequence and the presence of donor splice site patterns. In summary, our
findings yield a better understanding of the internal reasoning of a
convolutional neural network when applying such a neural network to genomic
data.


  The extinction of species is a core process that affects the diversity of
life on Earth. One way of investigating the causes and consequences of
extinctions is to build conceptual ecological models, and to use the dynamical
outcomes of such models to provide quantitative formalization of changes to
Earth's biosphere. In this paper we propose and study a conceptual resource
model that describes a simple and easily understandable mechanism for resource
competition, generalizes the well-known Huisman and Weissing model, and takes
into account species self-regulation, extinctions, and time dependence of
resources. We use analytical investigations and numerical simulations to study
the dynamics of our model under chaotic and periodic climate oscillations, and
show that the stochastic dynamics of our model exhibit strong dependence on
initial parameters. We also demonstrate that extinctions in our model are
inevitable if an ecosystem has the maximal possible biodiversity and uses the
maximal amount of resources. Our conceptual modeling provides theoretical
support for suggestions that non-linear processes were important during major
extinction events in Earth history.


  Mathematicians have always been attracted to the field of genetics. I am
especially interested in the mathematical aspects of research on homosexuality.
Certain studies show that male homosexuality may have a genetic component that
is correlated with female fertility. Other studies show the existence of the
fraternal birth order effect, that is the correlation of homosexuality with the
number of older brothers.
  This paper is devoted to the mathematical aspects of how these two phenomena
are interconnected. In particular, I show that the fraternal birth order effect
produces a correlation between homosexuality and maternal fecundity. Vice
versa, I show that the correlation between homosexuality and female fecundity
implies the increase of the probability of the younger brothers being
homosexual.


  Pure and homogeneous biological macromolecules (i.e. proteins, nucleic acids,
protein-protein or protein-nucleic acid complexes, and functional assemblies
such as ribosomes and viruses) are the key for consistent and reliable
biochemical and biophysical measurements, as well as for reproducible
crystallizations, best crystal diffraction properties, and exploitable electron
microscopy images. Highlights: Pure and homogeneous macromolecules are the key
for the best experimental results; They warrant the consistency and the
reliability of biochemical and biophysical data; They give more reproducible
crystallography and electron microscopy results as well.


  The theory of life history evolution provides a powerful framework to
understand the evolutionary dynamics of pathogens in both epidemic and endemic
situations. This framework, however, relies on the assumption that pathogen
populations are very large and that one can neglect the effects of demographic
stochasticity. Here we expand the theory of life history evolution to account
for the effects of finite population size on the evolution of pathogen
virulence. We show that demographic stochasticity introduces additional
evolutionary forces that can qualitatively affect the dynamics and the
evolutionary outcome. We discuss the importance of the shape of pathogen
fitness landscape and host heterogeneity on the balance between mutation,
selection and genetic drift. In particular, we discuss scenarios where finite
population size can dramatically affect classical predictions of deterministic
models. This analysis reconciles Adaptive Dynamics with population genetics in
finite populations and thus provides a new theoretical toolbox to study
life-history evolution in realistic ecological scenarios.


  Network growth as described by the Duplication-Divergence model proposes a
simple general idea for the evolution dynamics of natural networks. In
particular it is an alternative to the well known Barab\'asi-Albert model when
applied to protein-protein interaction networks. In this work we derive a
master equation for the node degree distribution of networks growing via
Duplication and Divergence and we obtain an expression for the total number of
links and for the degree distribution as a function of the number of nodes.
Using algebra tools we investigate the degree distribution asymptotic behavior.
Analytic results show that the network nodes average degree converges if the
total mutation rate is greater than 0.5 and diverges otherwise. Treating
original and duplicated node mutation rates as independent parameters has no
effect on this result. However, difference in these parameters results in a
slower rate of convergence and in different degree distributions. The more
different these parameters are, the denser the tail of the distribution. We
compare the solutions obtained with simulated networks. These results are in
good agreement with the expected values from the derived expressions. The
method developed is a robust tool to investigate other models for network
growing dynamics.


  We show that many topological features of level-1 species networks are
identifiable from the distribution of the gene tree quartets under the network
multi-species coalescent model. In particular, every cycle of size at least 4
and every hybrid node in a cycle of size at least 5 is identifiable. This is a
step toward justifying the inference of such networks which was recently
implemented by Sol\'is-Lemus and An\'e. We show additionally how to compute
quartet concordance factors for a network in terms of simpler networks, and
explore some circumstances in which cycles of size 3 and hybrid nodes in
4-cycles can be detected.


  DNA shape readout is an important mechanism of target site recognition by
transcription factors, in addition to the sequence readout. Several models of
transcription factor-DNA binding which consider DNA shape have been developed
in recent years. We present a new biophysical model of protein-DNA interaction
by considering the DNA shape features, which is based on a neighbour
dinucleotide dependency model BayesPI2. The parameters of the new model are
restricted to a subspace spanned by the 2-mer DNA shape features, which
allowing a biophysical interpretation of the new parameters as
position-dependent preferences towards certain values of the features. Using
the new model, we explore the variation of DNA shape preferences in several
transcription factors across cancer cell lines and cellular conditions. We find
evidence of DNA shape variations at FOXA1 binding sites in MCF7 cells after
treatment with steroids. The new model is useful for elucidating finer details
of transcription factor-DNA interaction. It may be used to improve the
prediction of cancer mutation effects in the future.


  Recently, a framework considering RNA sequences and their RNA secondary
structures as pairs, led to some information-theoretic perspectives on how the
semantics encoded in RNA sequences can be inferred. In this context, the
pairing arises naturally from the energy model of RNA secondary structures.
Fixing the sequence in the pairing produces the RNA energy landscape, whose
partition function was discovered by McCaskill. Dually, fixing the structure
induces the energy landscape of sequences. The latter has been considered for
designing more efficient inverse folding algorithms. We present here the
Hamming distance filtered, dual partition function, together with a Boltzmann
sampler using novel dynamic programming routines for the loop-based energy
model. The time complexity of the algorithm is $O(h^2n)$, where $h,n$ are
Hamming distance and sequence length, respectively, reducing the time
complexity of samplers, reported in the literature by $O(n^2)$. We then present
two applications, the first being in the context of the evolution of natural
sequence-structure pairs of microRNAs and the second constructing neutral
paths. The former studies the inverse fold rate (IFR) of sequence-structure
pairs, filtered by Hamming distance, observing that such pairs evolve towards
higher levels of robustness, i.e.,~increasing IFR. The latter is an algorithm
that constructs neutral paths: given two sequences in a neutral network, we
employ the sampler in order to construct short paths connecting them,
consisting of sequences all contained in the neutral network.


  It is acknowledged that co-evolutionary nucleotide-nucleotide interactions
are essential for RNA structures and functions. Currently, direct coupling
analysis (DCA) infers nucleotide contacts in a sequence from its homologous
sequence alignment across different species. DCA and similar approaches that
use sequence information alone usually yield a low accuracy, especially when
the available homologous sequences are limited. Here we present a new method
that incorporates a Restricted Boltzmann Machine (RBM) to augment the
information on sequence co-variations with structural patterns in contact
inference. We thus name our method DIRECT that stands for Direct Information
REweighted by Contact Templates. Benchmark tests demonstrate that DIRECT
produces a substantial enhancement of 13% in accuracy on average for contact
prediction in comparison to the traditional DCA. These results suggest that
DIRECT could be used for improving predictions of RNA tertiary structures and
functions. The source codes and dataset of DIRECT are available at http://
http://zhao.phy.ccnu.edu.cn:8122/DIRECT/index.html.


  The rate of recombination affects the mode of molecular evolution. In
high-recombining sequence, the targets of selection are individual genetic
loci; under low recombination, selection collectively acts on large,
genetically linked genomic segments. Selection under linkage can induce clonal
interference, a specific mode of evolution by competition of genetic clades
within a population. This mode is well known in asexually evolving microbes,
but has not been traced systematically in an obligate sexual organism. Here we
show that the Drosophila genome is partitioned into two modes of evolution: a
local interference regime with limited effects of genetic linkage, and an
interference condensate with clonal competition. We map these modes by
differences in mutation frequency spectra, and we show that the transition
between them occurs at a threshold recombination rate that is predictable from
genomic summary statistics. We find the interference condensate in segments of
low-recombining sequence that are located primarily in chromosomal regions
flanking the centromeres and cover about 20% of the Drosophila genome.
Condensate regions have characteristics of asexual evolution that impact gene
function: the efficacy of selection and the speed of evolution are lower and
the genetic load is higher than in regions of local interference. Our results
suggest that multicellular eukaryotes can harbor heterogeneous modes and tempi
of evolution within one genome. We argue that this variation generates
selection on genome architecture.


  Machine Learning techniques have been used to teach computer programs how to
play games as complicated as Chess and Go. These were achieved using powerful
tools such as Neural Networks and Parallel Computing on Supercomputers. In this
paper, we define a model of populational growth and evolution based on the idea
of Reinforcement Learning, but using only the 3 sources stated in the title
processed on a low-tier laptop. The model correctly predicts the development of
a population around food sources and their migration in search of a new one
when the known ones become saturated. Additionally, we compared our model to a
pure random one and the population number was fitted to a logistic function for
two interesting evolutions of the system.


  Multiple mechanisms in the HIV lifecycle play a role in its ability to evade
therapy and become a chronic, difficult-to-treat infection. Within its major
cellular target, the activated T cell, many steps occur between viral entry and
viral burst, including reverse transcription of viral RNA, integration of the
viral DNA in the host genome, viral transcription, splicing, translation, host
and viral regulation, and viral packaging. These steps exploit complex networks
of macromolecular interactions that exhibit various forms of stochastic
behavior. While each of the steps of HIV infection have been individually
studied extensively, the combinatorial contribution of rare events in each of
the steps, and how series of these rare events lead to different infection
phenotypes, are not well understood. The complexity of these processes render
experimental study challenging. Therefore, we have built a comprehensive
computational model of this large system, by collating the community's
knowledge of the infection process. It is a stochastic model where rates of
different events in the system are represented as probabilities of the event
occurring in a timestep of the simulation. This model enables an understanding
of the noise and variation in the system. The model also facilitates a
dissected understanding of each small part of the large complex system, and its
impact on the overall system dynamics.


  Competition is the main driver of population dynamics, which shapes the
genetic composition of populations and the assembly of ecological communities.
Neutral models assume that all the individuals are equivalent and that the
dynamics is governed by demographic (shot) noise, with a steady state species
abundance distribution (SAD) that reflects a mutation-extinction equilibrium.
Recently, many empirical and theoretical studies emphasized the importance of
environmental variations that affect coherently the relative fitness of entire
populations. Here we consider two generic time-averaged neutral models, in both
the relative fitness of each species fluctuates independently in time but its
mean is zero. The first (model A) describes a system with local competition and
linear fitness-dependence of the birth-death rates, while in the second (model
B) the competition is global and the fitness dependence is nonlinear. Due to
this nonlinearity, model B admits a noise-induced stabilization mechanism that
facilitates the invasion of new mutants. A self-consistent mean-field approach
is used to reduce the multi-species problem to two-species dynamics, and the
large-$N$ asymptotics of the emerging set of Fokker-Planck equations is
presented and solved. Our analytic expressions are shown to fit the SADs
obtained from extensive Monte-Carlo simulations and from numerical solutions of
the corresponding master equations.


  The Luruaco Lake located in the Department of Atl\'antico, Colombia, is
damaged by the discharge of untreated sewage, bringing risks to the health of
all who use its waters. The present study aims to perform the numerical
simulation of the concentration dynamics of fecal coliforms in the lake. The
simulation of the hydrodynamic flow is carried out by means of a
two-dimensional horizontal (2DH) model, given by a Navier-Stokes system. The
simulation of fecal coliform transport is described by a
convective-dispersive-reactive equation. These equations are solved numerically
by the Finite Difference Method (FDM) and the Mark and Cell (MAC) method, in
generalized coordinates. Regarding the construction of the computational mesh
of the Luruaco Lake, the cubic spline and multiblock methods were used. The
results obtained in the simulations allow a better understanding of the
dynamics of fecal coliforms in the Luruaco Lake, showing the more polluted
regions. They can also advise public agencies on identifying the emitters of
pollutants in the lake and on developing an optimal treatment for the recovery
of the polluted environment.


  Optical molecular sensing techniques are often limited by the refractive
index change associated with the probed interactions. In this work, we present
a closed form analytical model to estimate the magnitude of optical refractive
index change arising from protein-protein interactions. The model, based on the
Maxwell Garnett effective medium theory and first order chemical kinetics
serves as a general framework for estimating the detection limits of optical
sensing of molecular interactions. The model is applicable to situations where
one interacting species is immobilized to a surface, as commonly done, or to
emerging techniques such as Back-Scattering Interferometry (BSI) where both
interacting species are un-tethered. Our findings from this model point to the
strong role of as yet unidentified factors in the origin of the BSI signal
resulting in significant deviation from linear optical response.


  We consider a population evolving under mutation and selection. The genotype
of an individual is a word of length $\ell$ over a finite alphabet. Mutations
occur during reproduction, independently on each locus; the fitness depends on
the Hamming class (the distance to a reference sequence $w^*$). Evolution is
driven according to the classical Wright--Fisher process. We focus on the
proportion of the different classes under the invariant measure of the process.
We consider the regime where the length of the genotypes $\ell$ goes to
infinity, and both the population size and the inverse of the mutation rate are
of order $\ell$. We prove the existence of a critical curve, which depends both
on the population size and the mutation rate. Below the critical curve, the
proportion of any fixed class converges to $0$, whereas above the curve, it
converges to a positive quantity, for which we give an explicit formula.


  The adaptive immune system responds to foreign invaders by activating T
cells. However, this response is perilous if T cells are activated by the wrong
signal or if they remain activated unduly long after the threat has been
eliminated. It is therefore important that T cells get activated only by the
right kind of signals and for the right duration. The dominant theory in
immunology over recent decades has been that a T cell must receive at least two
signals before it can become activated. It is, however, unclear whether and how
this two-signal requirement ensures that T cell activation is provoked only by
the right signals and that the response is just long enough. Here, we propose
that the two signal requirement induces a novel coherent feedforward motif
whose properties align with those preferred for an ideal immune response.
Further consideration of the interaction between helper and regulatory T cells
induces a composite feedforward-feedback motif, analysis of which generates a
condition for a healthy concentration of activated T helper cells.


  Major histocompatibility complex class two (MHC-II) molecules are
trans-membrane proteins and key components of the cellular immune system. Upon
recognition of foreign peptides expressed on the MHC-II binding groove, helper
T cells mount an immune response against invading pathogens. Therefore,
mechanistic identification and knowledge of physico-chemical features that
govern interactions between peptides and MHC-II molecules is useful for the
design of effective epitope-based vaccines, as well as for understanding of
immune responses. In this paper, we present a comprehensive trans-allelic
prediction model, a generalized version of our previous biophysical model, that
can predict peptide interactions for all three human MHC-II loci (HLA-DR,
HLA-DP and HLA-DQ), using both peptide sequence data and structural information
of MHC-II molecules. The advantage of this approach over other machine learning
models is that it offers a simple and plausible physical explanation for
peptide-MHC-II interactions. We train the model using a benchmark experimental
dataset, and measure its predictive performance using novel data. Despite its
relative simplicity, we find that the model has comparable performance to the
state-of-the-art method. Focusing on the physical bases of peptide-MHC binding,
we find support for previous theoretical predictions about the contributions of
certain binding pockets to the binding energy. Additionally, we find that
binding pockets P 4 and P 5 of HLA-DP, which were not previously considered as
primary anchors, do make strong contributions to the binding energy. Together,
the results indicate that our model can serve as a useful complement to
alternative approaches to predicting peptide-MHC interactions.


  This chapter gives a graceful introduction to problem of protein three-
dimensional structure prediction, and focuses on how to make structural sense
out of a single input sequence with unknown structure, the 'query' or 'target'
sequence. We give an overview of the different classes of modelling techniques,
notably template-based and template free. We also discuss the way in which
structural predictions are validated within the global com- munity, and
elaborate on the extent to which predicted structures may be trusted and used
in practice. Finally we discuss whether the concept of a sin- gle fold
pertaining to a protein structure is sustainable given recent insights. In
short, we conclude that the general protein three-dimensional structure
prediction problem remains unsolved, especially if we desire quantitative
predictions. However, if a homologous structural template is available in the
PDB model or reasonable to high accuracy may be generated.


  This chapter deals with approaches for protein three-dimensional structure
prediction, starting out from a single input sequence with unknown struc- ture,
the 'query' or 'target' sequence. Both template based and template free
modelling techniques are treated, and how resulting structural models may be
selected and refined. We give a concrete flowchart for how to de- cide which
modelling strategy is best suited in particular circumstances, and which steps
need to be taken in each strategy. Notably, the ability to locate a suitable
structural template by homology or fold recognition is crucial; without this
models will be of low quality at best. With a template avail- able, the quality
of the query-template alignment crucially determines the model quality. We also
discuss how other, courser, experimental data may be incorporated in the
modelling process to alleviate the problem of missing template structures.
Finally, we discuss measures to predict the quality of models generated.


  Urban outputs often scale superlinearly with city population. A difficulty in
understanding the mechanism of this phenomenon is that different outputs differ
considerably in their scaling behaviors. Here, we formulate a physics-based
model for the origin of superlinear scaling in urban outputs by treating human
interaction as a random process. Our model suggests that the increased
likelihood of finding required collaborations in a larger population can
explain this superlinear scaling, which our model predicts to be non-power-law.
Moreover, the extent of superlinearity should be greater for activities that
require more collaborators. We test this model using a novel dataset for seven
crime types and find strong support.


  Non-Gaussian diffusion is commonly considered as a result of fluctuating
diffusivity, which is correlated in time or in space or both. In this work, we
investigate the non-Gaussian diffusion in static disordered media via a
quenched trap model, where the diffusivity is spatially correlated. Several
unique effects due to quenched disorder are reported. We analytically estimate
the diffusion coefficient $D_{\text{dis}}$ and its fluctuation over samples of
finite size. We show a mechanism of population splitting in the non-Gaussian
diffusion. It results in a sharp peak in the distribution of displacement
$P(x,t)$ around $x=0$, that has frequently been observed in experiments. We
examine the fidelity of the coarse-grained diffusion map, which is
reconstructed from particle trajectories. Finally, we propose a procedure to
estimate the correlation length in static disordered environments, where the
information stored in the sample-to-sample fluctuation has been utilized.


  In long-term potentiation (LTP), one of the most studied types of neural
plasticity, synaptic strength is persistently increased in response to
stimulation. Although a number of different proteins have been implicated in
the sub-cellular molecular processes underlying induction and maintenance of
LTP, the precise mechanisms remain unknown. A particular challenge is to
demonstrate that a proposed molecular mechanism can provide the level of
stability needed to maintain memories for months or longer, in spite of the
fact that many of the participating molecules have much shorter life spans.
Here we present a computational model that combines simulations of several
biochemical reactions that have been suggested in the LTP literature and show
that the resulting system does exhibit the required stability. At the core of
the model are two interlinked feedback loops of molecular reactions, one
involving the atypical protein kinase PKM{\zeta} and its messenger RNA, the
other involving PKM{\zeta} and GluA2-containing AMPA receptors. We demonstrate
that robust bistability - stable equilibria both in the synapse's potentiated
and unpotentiated states - can arise from a set of simple molecular reactions.
The model is able to account for a wide range of empirical results, including
induction and maintenance of late-phase LTP, cellular memory reconsolidation
and the effects of different pharmaceutical interventions.


  De novo prediction of protein folding is an open scientific challenge. Many
folding models and force fields have been developed, yet all face difficulties
converging to native conformations. Hydrophobicity scales (HSs) play a crucial
role in such simulations as they define the energetic interactions between
protein residues, thus determining the energetically favorable conformation.
While many HSs have been developed over the years using various methods, it is
surprising that the scales show very weak consensus in their assignment of
hydrophobicity indexes to the various residues. In this work, several HSs are
systematically assessed via atomistic Monte Carlo simulation of folding of
small proteins, by converting the HSs of interest into residue-residue contact
energy matrices. HSs that poorly preserve native structures of proteins were
tuned by applying a linear transformation. Subsequently, folding simulations
were used to examine the ability of the HSs to correctly fold the proteins from
a random initial conformation. Root mean square deviation (RMSD) and energy of
the proteins during folding were sampled and used to define an ER-score, as the
correlation between the 2-dimensional energy-RMSD (ER) histogram with 50%
lowest energy conformations and the ER histogram with 50% lowest RMSD
conformations. Thus, we were able to compare the ability of the different HSs
to predict de novo protein folding quantitatively.


  Amyloid fibrils are stable aggregates of misfolded proteins and polypeptides
that are insoluble and resistant to protease activity. Abnormal formation of
amyloid fibrils in vivo may lead to neurodegenerative disorders and other
systemic amyloidosis such as Alzheimer's, Parkinson's, and atherosclerosis.
Because of their clinical importance amyloids are found under intense
scientific research. Amyloidogenic sequences of short polypeptide segments
within proteins are responsible for the transformation of correctly folded
proteins into parts of larger amyloid fibrils. The {\alpha}-helical secondary
structure is believed to host many amyloidogenic sequences and be a key player
in different stages of the amyloidogenesis process. Most of the studies on
amyloids focus on the role of amyloidogenic sequences. The focus of this study
is the relation between amyloidogenicity and the structure of the amyloidogenic
{\alpha}-helical sequence. We have previously shown that the {\alpha}-helical
conformation may be expressed by two parameters ({\theta} and \{rho}) that form
orthogonal coordinates based on the Ramachandran dihedrals ({\phi} and {\psi})
and provide an illuminating interpretation of the {\alpha}-helical
conformation. By performing statistical analysis on {\alpha}-helical
conformations found in the protein data bank, an apparent relation between
{\alpha}-helical conformation, as expressed by {\theta} and \{rho}, and
amyloidogenicity is revealed. Remarkably, random amino acid sequences, whose
helical structure was obtained from the most probably dihedral angles as
obtained from PDB data, revealed the same dependency of amyloidogenicity,
suggesting the importance of {\alpha}-helical structure as opposed to sequence.


  Discovering disease pathways, which can be defined as sets of proteins
associated with a given disease, is an important problem that has the potential
to provide clinically actionable insights for disease diagnosis, prognosis, and
treatment. Computational methods aid the discovery by relying on
protein-protein interaction (PPI) networks. They start with a few known
disease-associated proteins and aim to find the rest of the pathway by
exploring the PPI network around the known disease proteins. However, the
success of such methods has been limited, and failure cases have not been well
understood. Here we study the PPI network structure of 519 disease pathways. We
find that 90% of pathways do not correspond to single well-connected components
in the PPI network. Instead, proteins associated with a single disease tend to
form many separate connected components/regions in the network. We then
evaluate state-of-the-art disease pathway discovery methods and show that their
performance is especially poor on diseases with disconnected pathways. Thus, we
conclude that network connectivity structure alone may not be sufficient for
disease pathway discovery. However, we show that higher-order network
structures, such as small subgraphs of the pathway, provide a promising
direction for the development of new methods.


  We investigate the symmetry of circular genome rearrangement models, discuss
the implementation of a new representation-theoretic method of calculating
evolutionary distances between circular genomes, and give the results of some
initial calculations for genomes with up to 11 regions.


  Co-infecting parasites and pathogens remain a leading challenge for global
public health due to their consequences for individual-level infection risk and
disease progression. However, a clear understanding of the population-level
consequences of co-infection is lacking. Here, we constructed a model that
includes three individual-level effects of co-infection: mortality, fecundity,
and transmission. We used the model to investigate how these individual-level
consequences of co-infection scale up to produce population-level infection
patterns. To parameterize this model, we conducted a four-year cohort study in
African buffalo to estimate the individual-level effects of co-infection with
two bacterial pathogens, bovine tuberculosis (BTB) and brucellosis, across a
range of demographic and environmental contexts. At the individual-level, our
empirical results identified BTB as a risk factor for acquiring brucellosis,
but we found no association between brucellosis and the risk of acquiring BTB.
Both infections were associated with reductions in survival and neither
infection was associated with reductions in fecundity. Results of the model
reproduce co-infection patterns in the data and predict opposite impacts of
co-infection at individual and population scales: whereas BTB facilitated
brucellosis infection at the individual-level, our model predicts the presence
of brucellosis to have a strong negative impact on BTB at the population-level.
In modeled populations where brucellosis is present, the endemic prevalence and
basic reproduction number (Ro) of BTB were lower than in populations without
brucellosis. Therefore, these results provide a data-driven example of
competition between co-infecting pathogens that occurs when one pathogen
facilitates secondary infections at the individual level.


  Essential to each other, growth and exploration are jointly observed in
populations, be it alive such as animals and cells or inanimate such as goods
and money. But their ability to move, crucial to cope with uncertainty and
optimize returns, is tempered by the space/time properties of the environment.
We investigate how the environment shape optimal growth and population
distribution in such conditions. We uncover a trade-off between risks and
returns by revisiting a common growth model over general graphs. Our results
reveal a rich and nuanced picture: fruitful strategies commonly lead to risky
positions, but this tension may nonetheless be alleviated by the geometry of
the explored space. The applicability of our conclusions is subsequently
illustrated over an empirical study of financial data.


  Motivation: Identifying genomic variants is an essential step for connecting
genotype and phenotype. The usual approach consists of statistical inference of
variants from alignments of sequencing reads. State-of-the-art variant callers
can resolve a wide range of different variant types with high accuracy.
However, they require that all read alignments be available from the beginning
of variant calling and be sorted by coordinates. Sorting is computationally
expensive, both memory- and speed-wise, and the resulting pipelines suffer from
storing and retrieving large alignments files from external memory. Therefore,
there is interest in developing methods for resource-efficient variant calling.
  Results: We present Ococo, the first program capable of inferring variants in
a real-time, as read alignments are fed in. Ococo inputs unsorted alignments
from a stream and infers single-nucleotide variants, together with a genomic
consensus, using statistics stored in compact several-bit counters. Ococo
provides a fast and memory-efficient alternative to the usual variant calling.
It is particularly advantageous when reads are sequenced or mapped
progressively, or when available computational resources are at a premium.


  Collective motion is an intriguing phenomenon, especially considering that it
arises from a set of simple rules governing local interactions between
individuals. In theoretical models, these rules are normally \emph{assumed} to
take a particular form, possibly constrained by heuristic arguments. We propose
a new class of models, which describe the individuals as \emph{agents}, capable
of deciding for themselves how to act and learning from their experiences. The
local interaction rules do not need to be postulated in this model, since they
\emph{emerge} from the learning process. We apply this ansatz to a concrete
scenario involving marching locusts, in order to model the phenomenon of
density-dependent alignment. We show that our learning agent-based model can
account for a Fokker-Planck equation that describes the collective motion and,
most notably, that the agents can learn the appropriate local interactions,
requiring no strong previous assumptions on their form. These results suggest
that learning agent-based models are a powerful tool for studying a broader
class of problems involving collective motion and animal agency in general.


  Network analysis has driven key developments in research on animal behaviour
by providing quantitative methods to study the social structures of animal
groups and populations. A recent formalism, known as \emph{multilayer network
analysis}, has advanced the study of multifaceted networked systems in many
disciplines. It offers novel ways to study and quantify animal behaviour as
connected 'layers' of interactions. In this article, we review common questions
in animal behaviour that can be studied using a multilayer approach, and we
link these questions to specific analyses. We outline the types of behavioural
data and questions that may be suitable to study using multilayer network
analysis. We detail several multilayer methods, which can provide new insights
into questions about animal sociality at individual, group, population, and
evolutionary levels of organisation. We give examples for how to implement
multilayer methods to demonstrate how taking a multilayer approach can alter
inferences about social structure and the positions of individuals within such
a structure. Finally, we discuss caveats to undertaking multilayer network
analysis in the study of animal social networks, and we call attention to
methodological challenges for the application of these approaches. Our aim is
to instigate the study of new questions about animal sociality using the new
toolbox of multilayer network analysis.


  Conjugation is a process in which bacteria exchange DNA through a physical
connection (conjugative junction) between mating cells. Despite its
significance for processes such as the spread of antibiotic resistance, the
role of physical forces in conjugation is poorly understood. Here we use
computer models to show that the conjugative junction not only serves as a link
to transfer the DNA but it also mechanically stabilises the mating pair which
significantly increases the conjugation rate. We discuss the importance of our
findings for biological evolution and suggest experiments to validate them.


  RNA-binding proteins (RBPs) play crucial roles in many biological processes,
e.g. gene regulation. Computational identification of RBP binding sites on RNAs
are urgently needed. In particular, RBPs bind to RNAs by recognizing sequence
motifs. Thus, fast locating those motifs on RNA sequences is crucial and
time-efficient for determining whether the RNAs interact with the RBPs or not.
In this study, we present an attention based convolutional neural network,
iDeepA, to predict RNA-protein binding sites from raw RNA sequences. We first
encode RNA sequences into one-hot encoding. Next, we design a deep learning
model with a convolutional neural network (CNN) and an attention mechanism,
which automatically search for important positions, e.g. binding motifs, to
learn discriminant high-level features for predicting RBP binding sites. We
evaluate iDeepA on publicly gold-standard RBP binding sites derived from
CLIP-seq data. The results demonstrate iDeepA achieves comparable performance
with other state-of-the-art methods.


  Transcription factors are able to associate to their binding sites on DNA
faster than the physical limit posed by diffusion. Such high association rates
can be achieved by alternating between three-dimensional diffusion and
one-dimensional sliding along the DNA chain, a mechanism dubbed Facilitated
Diffusion. By studying a collection of transcription factor binding sites of
Escherichia coli from the RegulonDB database and of Bacillus subtilis from
DBTBS, we reveal a funnel in the binding energy landscape around the target
sequences. We show that such a funnel is linked to the presence of gradients of
AT in the base composition in the DNA region around the binding sites. An
extensive computational study of the stochastic sliding process along the
energetic landscapes obtained from the database shows that the funnel can
significantly enhance the probability of transcription factors to find their
target sequences when sliding in their proximity. We demonstrate that this
enhancement leads to a speed-up of the association process.


  Cooperation is a difficult proposition in the face of Darwinian selection.
Those that defect have an evolutionary advantage over cooperators who should
therefore die out. However, spatial structure enables cooperators to survive
through the formation of homogeneous clusters, which is the hallmark of network
reciprocity. Here we go beyond this traditional setup and study the
spatiotemporal dynamics of cooperation in a population of populations. We use
the prisoner's dilemma game as the mathematical model and show that considering
several populations simultaneously give rise to fascinating spatiotemporal
dynamics and pattern formation. Even the simplest assumption that strategies
between different populations are payoff-neutral with one another results in
the spontaneous emergence of cyclic dominance, where defectors of one
population become prey of cooperators in the other population, and vice versa.
Moreover, if social interactions within different populations are characterized
by significantly different temptations to defect, we observe that defectors in
the population with the largest temptation counterintuitively vanish the
fastest, while cooperators that hang on eventually take over the whole
available space. Our results reveal that considering the simultaneous presence
of different populations significantly expands the complexity of evolutionary
dynamics in structured populations, and it allow us to understand the stability
of cooperation under adverse conditions that could never be bridged by network
reciprocity alone.


  Living cells dynamically modulate the local morphologies of their actin
cytoskeletons to perform biological functions, including force transduction,
intracellular transport, and cell division. A major challenge is to understand
how diverse structures of the actin cytoskeleton are assembled from a limited
set of molecular building blocks. Here we study the spontaneous self-assembly
of a minimal model of cytoskeletal materials, consisting of semiflexible actin
filaments, crosslinkers, and molecular motors. Using coarse-grained
simulations, we demonstrate that by changing concentrations and kinetics of
crosslinkers and motors we can generate three distinct structural phases of
actomyosin assemblies: bundled, polarity-sorted, and contracted. We introduce
new metrics to distinguish these structural phases and demonstrate their
functional roles. We find that the binding kinetics of motors and crosslinkers
can be tuned to optimize contractile force generation, motor transport, and
mechanical response. By quantitatively characterizing the relationships between
modes of cytoskeletal self-assembly, the resulting structures, and their
functional consequences, our work suggests new principles for the design of
active materials.


  We present an elaborate framework for formally modelling pathways in chemical
reaction networks on a mechanistic level. Networks are modelled mathematically
as directed multi-hypergraphs, with vertices corresponding to molecules and
hyperedges to reactions. Pathways are modelled as integer hyperflows and we
expand the network model by detailed routing constraints. In contrast to the
more traditional approaches like Flux Balance Analysis or Elementary Mode
analysis we insist on integer-valued flows. While this choice makes it
necessary to solve possibly hard integer linear programs, it has the advantage
that more detailed mechanistic questions can be formulated. It is thus possible
to query networks for general transformation motifs, and to automatically
enumerate optimal and near-optimal pathways. Similarities and differences
between our work and traditional approaches in metabolic network analysis are
discussed in detail. To demonstrate the applicability of the mathematical
framework to real-life problems we first explore the design space of possible
non-oxidative glycolysis pathways and show that recent manually designed
pathways can be further optimised. We then use a model of sugar chemistry to
investigate pathways in the autocatalytic formose process. A graph
transformation-based approach is used to automatically generate the reaction
networks of interest.


  Recent data indicate one or more moderately nearby supernovae in the early
Pleistocene, with additional events likely in the Miocene. This has motivated
more detailed computations, using new information about the nature of
supernovae and the distances of these events to describe in more detail the
sorts of effects that are indicated at the Earth. This short
communication/review is designed to describe some of these effects so that they
may possibly be related to changes in the biota around these times.


  Distinguishing between the capacity of ecosystems to generate ecosystem
services (ES) and the actual use of these service (ES flow) in ES assessment
and mapping is important to develop an understanding of the sustainability of
ES use. This study assesses the spatial variation in ES capacity and flow in
the Mediterranean small island state of Malta. The services included in this
study were crop provisioning, beekeeping and honey production, fodder and
livestock production, crop pollination, air quality regulation, and aesthetic
ES. This assessment develops different spatial models, which make use of
available datasets, causal relationships between datasets, including a
generated land use land cover (LULC) map, and statistical models and indicators
based on direct measurements. Individual ES indicators were mapped to visualise
and compare their spatial patterns across the case study area. Subsequently, an
analysis of ES associations and bundles was carried out using Pearson
parametric correlation test, for both ES capacity and flow indicators generated
from this study, and through Principal Component Analysis. Results demonstrate
several significant synergistic interactions between ES capacity and flow in
rural landscapes characterised with agricultural and semi-natural LULC
categories, indicating high landscape multifunctionality. In contrast,
predominantly urban areas tend to be characterised with a low ecosystem
capacity and ES flow, suggesting that ES delivery in the landscapes of the
study area is determined by land use intensity. These findings support the
notion that multifunctional rural landscapes provide multiple ES, making an
important contribution to human well-being, and that land use planning that
develops green infrastructure in urban areas can significantly contribute to
support biodiversity and ES delivery.


  Quantum calculations on the voltage sensing domain (VSD) of the Kv1.2
potassium channel (pdb: 3Lut)have been carried out on a 904 atoms subset of the
VSD, plus 24 water molecules. Side chains pointing away from the center of the
VSD were truncated; S1,S2,S3 end atoms were were fixed (all calculations); S4
end atoms could be fixed or free. Open conformations (membrane potentials >= 0)
closely match the known X-ray structure of the open state with salt bridges in
the in the VSD not ionized (H+ on the acid) whether S4 end atoms were fixed or
free (slightly closer fixed than free).The S4 segment backbone, free or not,
moves less than 2.5 A for positive to negative membrane potential switches, not
entirely in the expected direction, leaving H+ motion as the principal
component of the gating current. Groups of 3 - 5 side chains are important for
proton transport, based on the calculations. A proton transfers from tyrosine
(Y266), through arginine (R300), to glutamate (E183), accounting for
approximately 20 - 25% of the gating charge. Clusters of amino acids that can
transfer protons (acids, bases, tyrosine, histidine) are the main paths for
proton transport. A group of five amino acids, bounded by the conserved
aromatic F233, appears to exchange a proton. Dipole rotations may also
contribute. A proton path (calculations still in progress) is proposed for the
remainder of the VSD, suggesting a hypothesis for a complete gating mechanism.


  In the present work, we study the properties of biological networks by
applying analogous notions of fundamental concepts in Riemannian geometry and
optimal mass transport to discrete networks described by weighted graphs.
Specifically, we employ possible generalizations of the notion of Ricci
curvature on Riemannian manifold to discrete spaces in order to infer certain
robustness properties of the networks of interest. We compare three possible
discrete notions of Ricci curvature (Olivier Ricci curvature, Bakry-\'Emery
Ricci curvature, and Forman Ricci curvature) on some model and biological
networks. While the exact relationship of each of the three definitions of
curvature to one another is still not known, they do yield similar results on
our biological networks of interest. These notions are initially defined on
positively weighted graphs; however, Forman-Ricci curvature can also be defined
on directed positively weighted networks. We will generalize this notion of
directed Forman Ricci curvature on the network to a form that also considers
the signs of the directions (e. g., activator and repressor in transcription
networks). We call this notion the \emph{signed-control Ricci curvature}. Given
real biological networks are almost always directed and possess positive and
negative controls, the notion of signed-control curvature can elucidate the
network-based study of these complex networks. Finally, we compare the results
of these notions of Ricci curvature on networks to some known network measures,
namely, betweenness centrality and clustering coefficients on graphs.


  Here we show that deciding whether two rooted binary phylogenetic trees on
the same set of taxa permit a cherry-picking sequence, a special type of
elimination order on the taxa, is NP-complete. This improves on an earlier
result which proved hardness for eight or more trees. Via a known equivalence
between cherry-picking sequences and temporal phylogenetic networks, our result
proves that it is NP-complete to determine the existence of a temporal
phylogenetic network that contains topological embeddings of both trees. The
hardness result also greatly strengthens previous inapproximability results for
the minimum temporal-hybridization number problem. This is the optimization
version of the problem where we wish to construct a temporal phylogenetic
network that topologically embeds two given rooted binary phylogenetic trees
and that has a minimum number of indegree-2 nodes, which represent events such
as hybridization and horizontal gene transfer. We end on a positive note,
pointing out that fixed parameter tractability results in this area are likely
to ensure the continued relevance of the temporal phylogenetic network model.


  This preprint presents a web app (essentially a web page-based program) with
which two or more users (peers) can view and handle 3D molecular structures in
a concurrent, interactive way through their web browsers. This means they can
share orientation and zoom level, commands and other operations in almost real
time over the Internet through standard web pages. This web app, open source
and built with the open source components JSmol for molecular visualization and
Peer.js for WebRTC connection, provides a practical tool for online
collaboration and teaching at a distance. More broadly, it illustrates the
strong integrability of technologies for client-side web programming, and paves
the way for similar apps for concurrent work in other disciplines. Web app is
available at:
http://lucianoabriata.altervista.org/jsinscience/concurrent-jsmol/concurrent-jsmol-visualization.html


  Oscillations play a major role in a number of biological systems, from
predator-prey models of ecology to circadian clocks. In this paper we focus on
the question of whether oscillations exist within dual-site phosphorylation
systems. Previously, Wang and Sontag showed, using monotone systems theory,
that the Michaelis-Menten (MM) approximation of the distributive and sequential
dual-site phosphorylation system lacks oscillations. However, biological
systems are generally not purely distributive; there is generally some
processive behavior as well. Accordingly, this paper focuses on the MM
approximation of a general sequential dual-site phosphorylation system that
contains both processive and distributive components, termed the composite
system. Expanding on the methods of Bozeman and Morales, we preclude
oscillations in the MM approximation of the composite system. This implies the
lack of oscillations in the MM approximations of the processive and
distributive systems, shown previously, as well as in the MM approximation of
the partially processive and partially distributive mixed-mechanism system.


  Noise plays a fundamental role in a wide variety of physical and biological
dynamical systems. It can arise from an external forcing or due to random
dynamics internal to the system. It is well established that even weak noise
can result in large behavioral changes such as transitions between or escapes
from quasi-stable states. These transitions can correspond to critical events
such as failures or extinctions that make them essential phenomena to
understand and quantify, despite the fact that their occurrence is rare. This
article will provide an overview of the theory underlying the dynamics of rare
events for stochastic models along with some example applications.


  We review the major progress in the rigorous analysis of the classical
quasispecies model that usually comes in two related but different forms: the
Eigen model and the Crow--Kimura model. The model itself was formulated almost
50 years ago, and in its stationary form represents an easy to formulate
eigenvalue problem. Notwithstanding the simplicity of the problem statement, we
still lack full understanding of the behavior of the mean population fitness
and the quasispecies distribution for an arbitrary fitness landscape. Our main
goal in this review is two-fold: First, to highlight a number of impressive
mathematical results, including some of the recent ones, which pertain to the
mathematical development of the quasispecies theory. Second, to emphasize that,
despite these 50 years of vigorous research, there are still very natural both
biological and mathematical questions that remain to be addressed within the
quasispecies framework. Our hope is that at least some of the approaches we
review in this text can be of help for anyone embarking on further analysis of
the quasispecies model.


  We investigated the impact of hydrodynamic interactions (HI) on protein
folding using a coarse-grained model. The extent of the impact of hydrodynamic
interactions, whether it accelerates, retards, or has no effect on protein
folding, has been controversial. Together with a theoretical framework of the
energy landscape theory (ELT) for protein folding that describes the dynamics
of the collective motion with a single reaction coordinate across a folding
barrier, we compared the kinetic effects of HI on the folding rates of two
protein models that use a chain of single beads with distinctive topologies: a
64-residue alpha/beta chymotrypsin inhibitor 2 (CI2) protein, and a 57-residue
beta-barrel alpha-spectrin src-Homology 3 domain (SH3) protein. When comparing
the protein folding kinetics simulated with Brownian dynamics in the presence
of HI to that in the absence of HI, we find that the effect of HI on protein
folding appears to have a crossover behavior about the folding temperature.
Meaning that at a temperature greater than the folding temperature, the
enhanced friction from the hydrodynamic solvents between the beads in an
unfolded configuration results in lowered folding rate; conversely, at a
temperature lower than the folding temperature, HI accelerates folding by the
backflow of solvent toward the native folded state. Additionally, the extent of
acceleration depends on the topology of a protein: for a protein like CI2,
where its folding nucleus is rather diffuse in a transition state, HI channels
the formation of contacts by favoring a major folding pathway in a complex free
energy landscape, thus accelerating folding. For a protein like SH3, where its
folding nucleus is already specific and less diffuse, HI matters less at a
temperature lower than the folding temperature. Our findings provide further
theoretical insight to protein folding kinetic experiments and simulations.


  Affinity maturation produces antibodies that bind antigens with high
specificity by accumulating mutations in the antibody sequence. Mapping out the
antibody-antigen affinity landscape can give us insight into the accessible
paths during this rapid evolutionary process. By developing a carefully
controlled null model for noninteracting mutations, we characterized epistasis
in affinity measurements of a large library of antibody variants obtained by
Tite-Seq, a recently introduced Deep Mutational Scan method yielding physical
values of the binding constant. We show that representing affinity as the
binding free energy minimizes epistasis. Yet, we find that epistatically
interacting sites contribute substantially to binding. In addition to negative
epistasis, we report a large amount of beneficial epistasis, enlarging the
space of high-affinity antibodies as well as their mutational accessibility.
These properties suggest that the degeneracy of antibody sequences that can
bind a given antigen is enhanced by epistasis - an important property for
vaccine design.


  Mitochondria of brown adipocyte (BA) are the main intracellular sites for
thermogenesis, which have been targeted for therapy to reduce obesity. However,
there are long-standing critique and debates about the ability of raising
cellular temperature by endogenous thermogenesis. Currently, the wrong
theoretical model gave about the five orders of magnitude less than facts.
Here, based on the first law of thermodynamics and thermal diffusion equation,
the deduced theoretical model of mitochondrial thermogenesis satisfies Laplace
equation, and is a special case for thermal diffusion equation. The model
settles the long-standing questioning about the ability of raising cellular
temperature by endogenous thermogenesis, and explains the thermogenic
characteristics of brown adipocyte. The model and calculations also suggest
that the number of free available proton is the major limiting factor for
endogenous thermogenesis and its speed.


  Recently, the minimum number of reticulation events that is required to
simultaneously embed a collection P of rooted binary phylogenetic trees into a
so-called temporal network has been characterized in terms of cherry-picking
sequences. Such a sequence is a particular ordering on the leaves of the trees
in P. However, it is well-known that not all collections of phylogenetic trees
have a cherry-picking sequence. In this paper, we show that the problem of
deciding whether or not P has a cherry-picking sequence is NP-complete for when
P contains at least eight rooted binary phylogenetic trees. Moreover, we use
automata theory to show that the problem can be solved in polynomial time if
the number of trees in P and the number of cherries in each such tree are
bounded by a constant.


  Throughout the last decade, we have seen much progress towards characterising
and computing the minimum hybridisation number for a set P of rooted
phylogenetic trees. Roughly speaking, this minimum quantifies the number of
hybridisation events needed to explain a set of phylogenetic trees by
simultaneously embedding them into a phylogenetic network. From a mathematical
viewpoint, the notion of agreement forests is the underpinning concept for
almost all results that are related to calculating the minimum hybridisation
number for when |P|=2. However, despite various attempts, characterising this
number in terms of agreement forests for |P|>2 remains elusive. In this paper,
we characterise the minimum hybridisation number for when P is of arbitrary
size and consists of not necessarily binary trees. Building on our previous
work on cherry-picking sequences, we first establish a new characterisation to
compute the minimum hybridisation number in the space of tree-child networks.
Subsequently, we show how this characterisation extends to the space of all
rooted phylogenetic networks. Moreover, we establish a particular hardness
result that gives new insight into some of the limitations of agreement
forests.


  Phylogenetic networks are an extension of phylogenetic trees which are used
to represent evolutionary histories in which reticulation events (such as
recombination and hybridization) have occurred. A central question for such
networks is that of identifiability, which essentially asks under what
circumstances can we reliably identify the phylogenetic network that gave rise
to the observed data? Recently, identifiability results have appeared for
networks relative to a model of sequence evolution that generalizes the
standard Markov models used for phylogenetic trees. However, these results are
quite limited in terms of the complexity of the networks that are considered.
In this paper, by introducing an alternative probabilistic model for evolution
along a network that is based on some ground-breaking work by Thatte for
pedigrees, we are able to obtain an identifiability result for a much larger
class of phylogenetic networks (essentially the class of so-called tree-child
networks). To prove our main theorem, we derive some new results for
identifying tree-child networks combinatorially, and then adapt some techniques
developed by Thatte for pedigrees to show that our combinatorial results imply
identifiability in the probabilistic setting. We hope that the introduction of
our new model for networks could lead to new approaches to reliably construct
phylogenetic networks.


  Little is known about clouded leopards (Neofelis nebulosa), who have a
vulnerable population that extends across southern Asia. We reviewed the
literature and synthesized what is known about their ecology and behavior. Much
of the published literature either note detections within and on the edges of
their range, or are anecdotal observations, many of which are decades if not
over a century old. Clouded leopards are a medium-sized felid, with distinctive
cloud-shape markings, and notably long canines relative to skull size.
Estimates for population densities range from 0.58 to 6.53 individuals per 100
km2. Only 7 clouded leopards have been tracked via radio-collars, and home
range estimates range from 33.6-39.7 km2 for females and 35.5-43.5 km2 for
males. Most accounts describe clouded leopards as nocturnal, but radio
telemetry studies showed that clouded leopards have arrhythmic activity
patterns, with highest activity in the morning followed by evening crepuscular
hours. There has never been a targeted study of clouded leopard diet, but
observations show that they consume a variety of animals, including ungulates,
primates, and rodents. We encourage future study of their population density
and range to inform conservation efforts, and ecological studies in order to
understand the species and its ecological niche.


  The spread of invasive species to new areas threatens the stability of
ecosystems and causes major economic losses in agriculture and forestry. We
propose a novel approach to minimizing the spread of an invasive species given
a limited intervention budget. We first model invasive species propagation
using Hawkes processes, and then derive closed-form expressions for
characterizing the effect of an intervention action on the invasion process. We
use this to obtain an optimal intervention plan based on an integer programming
formulation, and compare the optimal plan against several
ecologically-motivated heuristic strategies used in practice. We present an
empirical study of two variants of the invasive control problem: minimizing the
final rate of invasions, and minimizing the number of invasions at the end of a
given time horizon. Our results show that the optimized intervention achieves
nearly the same level of control that would be attained by completely
eradicating the species, with a 20% cost saving. Additionally, we design a
heuristic intervention strategy based on a combination of the density and life
stage of the invasive individuals, and find that it comes surprisingly close to
the optimized strategy, suggesting that this could serve as a good rule of
thumb in invasive species management.


  The largely intrinsically disordered phenylalanine-glycine-rich nucleoporins
(FG Nups) underline a selectivity mechanism, which enables the rapid
translocation of transport factors (TFs) through the nuclear pore complexes
(NPCs). Conflicting models of NPC transport have assumed that FG Nups undergo
different conformational transitions upon interacting with TFs. To selectively
characterize conformational changes in FG Nups induced by TFs we performed
small-angle neutron scattering (SANS) with contrast matching. Conformational
ensembles derived SANS data indicate an increase in the overall size of FG Nups
is associated with TF interaction. Moreover, the organization of the FG motif
in the interacting state is consistent with prior experimental analyses
defining that FG motifs undergo conformational restriction upon interacting
with TFs. These results provide structural insights into a highly dynamic
interaction and illustrate how functional disorder imparts rapid and selective
FG Nup / TF interactions.


  The spread of an infectious disease is known to change people's behavior,
which in turn affects the spread of disease. Adaptive network models that
account for both epidemic and behavior change have found oscillations, but in
an extremely narrow region of the parameter space, which contrasts with
intuition and available data. In this paper we propose a simple SIS epidemic
model on an adaptive network with time-delayed rewiring, and show that
oscillatory solutions are now present in a wide region of the parameter space.
Altering the transmission or rewiring rates reveals the presence of an endemic
bubble - an enclosed region of the parameter space where oscillations are
observed.


  We study quasi-species and closely related evolutionary dynamics like the
replicator-mutator equation in high dimensions. In particular, we show that
under certain conditions the fitness of almost all quasi-species becomes
independent of mutational probabilities and the initial frequency distributions
of the sequences in high dimensional sequence spaces. This result is the
consequence of the concentration of measure on a high dimensional hypersphere
and its extension to Lipschitz functions known %knows as the Levy's Lemma.
Therefore, evolutionary dynamics almost always yields the same value for
fitness of the quasi-species, independent of the mutational process and initial
conditions, and is quite robust to mutational changes and fluctuations in
initial conditions. Our results naturally extend to any Lipschitz function
whose input parameters are the frequencies of individual constituents of the
quasi-species. This suggests that the functional capabilities of high
dimensional quasi-species are robust to fluctuations in the mutational
probabilities and initial conditions. We discuss the consequences of our study
for the replicator-mutator equation.


  Human beings are talkative. What advantage did their ancestors find in
communicating so much? Numerous authors consider this advantage to be "obvious"
and "enormous". If so, the problem of the evolutionary emergence of language
amounts to explaining why none of the other primate species evolved anything
even remotely similar to language. What I propose here is to reverse the
picture. On closer examination, language resembles a losing strategy. Competing
for providing other individuals with information, sometimes striving to be
heard, makes apparently no sense within a Darwinian framework. At face value,
language as we can observe it should never have existed or should have been
counter-selected. In other words, the selection pressure that led to language
is still missing. The solution I propose consists in regarding language as a
social signaling device that developed in a context of generalized insecurity
that is unique to our species. By talking, individuals advertise their
alertness and their ability to get informed. This hypothesis is shown to be
compatible with many characteristics of language that otherwise are left
unexplained.


  The sample frequency spectrum (SFS), which describes the distribution of
mutant alleles in a sample of DNA sequences, is a widely used summary statistic
in population genetics. The expected SFS has a strong dependence on the
historical population demography and this property is exploited by popular
statistical methods to infer complex demographic histories from DNA sequence
data. Most, if not all, of these inference methods exhibit pathological
behavior, however. Specifically, they often display runaway behavior in
optimization, where the inferred population sizes and epoch durations can
degenerate to 0 or diverge to infinity, and show undesirable sensitivity of the
inferred demography to perturbations in the data. The goal of this paper is to
provide theoretical insights into why such problems arise. To this end, we
characterize the geometry of the expected SFS for piecewise-constant
demographic histories and use our results to show that the aforementioned
pathological behavior of popular inference methods is intrinsic to the geometry
of the expected SFS. We provide explicit descriptions and visualizations for a
toy model with sample size 4, and generalize our intuition to arbitrary sample
sizes n using tools from convex and algebraic geometry. We also develop a
universal characterization result which shows that the expected SFS of a sample
of size n under an arbitrary population history can be recapitulated by a
piecewise-constant demography with only k(n) epochs, where k(n) is between n/2
and 2n-1. The set of expected SFS for piecewise-constant demographies with
fewer than k(n) epochs is open and non-convex, which causes the above phenomena
for inference from data.


  We seek models for the genotype evolution of agricultural animals, animals
involved in primary production processes. Classical models for genotype
evolution have tended to be very simple in order that analytic methods may be
employed in their study. Unfortunately these models fail to describe processes
in artificially controlled populations including agricultural livestock. It is
particularly important {to describe such processes} in order to make better use
of the massive genotyping data becoming available. We describe an approach to
stochastically modeling the dynamics of a biallelic polymorphism herds under
conditions of controlled mating {and restriction of herds size from above}. The
system of stochastic differential equations that we propose is based on jump
diffusion processes to provide an effective platform for Monte Carlo
simulation. A feature of the model is the division of the {population} into a
{\it main herd} comprised of animals involved in the production process and a
{\it replacement herd} of {animals not currently in the production process,
typically juvenile animals}. This feature allows for exploring different
strategies for adding replacement animals to the main herd without altering the
part of the model concerned with the dynamics of the main herd. A discrete-time
version of the model has been developed which reflects the typical practice of
New Zealand dairy herd management. Our Monte Carlo simulation has demonstrated
that an isolated deme whose size is bounded above demonstrate size
stabilization at a level less than the control limit, it is looks like partial
extinction, the effect being well known in classic models.


  We analyze a multi-type age dependent model for cell populations subject to
unidirectional motion, in both a stochastic and deterministic framework. Cells
are distributed into successive layers; they may divide and move irreversibly
from one layer to the next. We adapt results on the large-time convergence of
PDE systems and branching processes to our context, where the Perron-Frobenius
or Krein-Rutman theorem can not be applied. We derive explicit analytical
formulas for the asymptotic cell number moments, and the stable age
distribution. We illustrate these results numerically and we apply them to the
study of the morphodynamics of ovarian follicles. We prove the structural
parameter identifiability of our model in the case of age independent division
rates. Using a set of experimental biological data, we estimate the model
parameters to fit the changes in the cell numbers in each layer during the
early stages of follicle development.


  The period and amplitude of biomolecular oscillators are functionally
important properties in multiple contexts. For a biomolecular oscillator, the
overall constraints in how tuning of amplitude affects period, and vice versa,
are generally unclear. Here we investigate this co-variation of the period and
amplitude in mathematical models of biomolecular oscillators using both
simulations and analytical approximations. We computed the amplitude-period
co-variation of eleven benchmark biomolecular oscillators as their parameters
were individually varied around a nominal value, classifying the various
co-variation patterns such as a simultaneous increase/ decrease in period and
amplitude. Next, we repeated the classification using a power norm-based
amplitude metric, to account for the amplitudes of the many biomolecular
species that may be part of the oscillations, finding largely similar trends.
Finally, we calculate "scaling laws" of period-amplitude co-variation for a
subset of these benchmark oscillators finding that as the approximated period
increases, the upper bound of the amplitude increases, or reaches a constant
value. Based on these results, we discuss the effect of different parameters on
the type of period-amplitude co-variation as well as the difficulty in
achieving an oscillation with large amplitude and small period.


  Dengue is the most prevalent arthropod-borne viral disease. Clinical symptoms
of dengue virus (DENV) infection range from classical mild dengue fever to
severe, life-threatening dengue shock syndrome. However, most DENV infections
cause few or no symptoms. Asymptomatic DENV-infected patients provide a unique
opportunity to decipher the host immune responses leading to virus elimination
without negative impact on t v 'health. We used an integrated approach of
transcriptional profiling and immunological analysis comparing a Cambodian
population of strictly asymptomatic viremic individuals with clinical dengue
patients. Whereas inflammatory pathways and innate immune responses were
similar between asymptomatic individuals and clinical dengue patients,
expression of proteins related to antigen presentation and subsequent T and B
cell activation pathways were differentially regulated, independent of viral
load or previous DENV infection. Feedback mechanisms controlled the immune
response in asymptomatic viremic individuals as demonstrated by increased
activation of T cell apoptosis-related pathways and Fc$\gamma$RIIB signaling
associated with decreased anti-DENV specific antibody concentrations. Taken
together, our data illustrate that symptom-free DENV infection in children is
determined by increased activation of the adaptive immune compartment and
proper control mechanisms leading to elimination of viral infection without
excessive immune activation, having implications for novel vaccine development
strategies.


  Cooperators benefit others with paying costs. Evolution of cooperation
crucially depends on the cost-benefit ratio of cooperation, denoted as $c$. In
this work, we investigate the infinitely repeated prisoner's dilemma for
various values of $c$ with four of the representative memory-one strategies,
i.e., unconditional cooperation, unconditional defection, tit-for-tat, and
win-stay-lose-shift. We consider replicator dynamics which deterministically
describes how the fraction of each strategy evolves over time in an
infinite-sized well-mixed population in the presence of implementation error
and mutation among the four strategies. Our finding is that this
three-dimensional continuous-time dynamics exhibits chaos through a bifurcation
sequence similar to that of a logistic map as $c$ varies. If mutation occurs
with rate $\mu \ll 1$, the position of the bifurcation sequence on the $c$ axis
is numerically found to scale as $\mu^{0.1}$, and such sensitivity to $\mu$
suggests that mutation may have non-perturbative effects on evolutionary paths.
It demonstrates how the microscopic randomness of the mutation process can be
amplified to macroscopic unpredictability by evolutionary dynamics.


  Individual heterogeneity is a key characteristic of many real-world systems,
from organisms to humans. However its role in determining the system's
collective dynamics is typically not well understood. Here we study how
individual heterogeneity impacts the system network dynamics by comparing
linking mechanisms that favor similar or dissimilar individuals. We find that
this heterogeneity-based evolution can drive explosive network behavior and
dictates how a polarized population moves toward consensus. Our model shows
good agreement with data from both biological and social science domains. We
conclude that individual heterogeneity likely plays a key role in the
collective development of real-world networks and communities, and cannot be
ignored.


  For a recently derived pairwise model of network epidemics with non-Markovian
recovery, we prove that under some mild technical conditions on the
distribution of the infectious periods, smaller variance in the recovery time
leads to higher reproduction number, and consequently to a larger epidemic
outbreak, when the mean infectious period is fixed. We discuss how this result
is related to various stochastic orderings of the distributions of infectious
periods. The results are illustrated by a number of explicit stochastic
simulations, suggesting that their validity goes beyond regular networks.


  Experimental studies of microbial communities routinely reveal that they have
multiple stable states. While each of these states is generally resilient,
certain perturbations such as antibiotics, probiotics and diet shifts, result
in transitions to other states. Can we reliably both predict such stable states
as well as direct and control transitions between them? Here we present a new
conceptual model inspired by the stable marriage problem in game theory and
economics in which microbial communities naturally exhibit multiple stable
states, each state with a different species' abundance profile. Our model's
core ingredient is that microbes utilize nutrients one at a time while
competing with each other. Using only two ranked tables, one with microbes'
nutrient preferences and one with their competitive abilities, we can determine
all possible stable states as well as predict inter-state transitions,
triggered by the removal or addition of a specific nutrient or microbe.
Further, using an example of 7 Bacteroides species common to the human gut
utilizing 9 polysaccharides, we predict that mutual complementarity in nutrient
preferences enables these species to coexist at high abundances.


  We propose generative neural network methods to generate DNA sequences and
tune them to have desired properties. We present three approaches: creating
synthetic DNA sequences using a generative adversarial network; a DNA-based
variant of the activation maximization ("deep dream") design method; and a
joint procedure which combines these two approaches together. We show that
these tools capture important structures of the data and, when applied to
designing probes for protein binding microarrays, allow us to generate new
sequences whose properties are estimated to be superior to those found in the
training data. We believe that these results open the door for applying deep
generative models to advance genomics research.


  Punishment and partner switching are two well-studied mechanisms that support
the evolution of cooperation. Observation of human behaviour suggests that the
extent to which punishment is adopted depends on the usage of alternative
mechanisms, including partner switching. In this study, we investigate the
combined effect of punishment and partner switching in evolutionary prisoner's
dilemma games conducted on a network. In the model, agents are located on the
network and participate in the prisoner's dilemma games with punishment. In
addition, they can opportunistically switch interaction partners to improve
their payoff. Our Monte Carlo simulation showed that a large frequency of
punishers is required to suppress defectors when the frequency of partner
switching is low. In contrast, cooperation is the most abundant strategy when
the frequency of partner switching is high regardless of the strength of
punishment. Interestingly, cooperators become abundant not because they avoid
the cost of inflicting punishment and earn a larger average payoff per game but
rather because they have more numerous opportunities to be referred as a role
agent by defectors. Our results imply that the fluidity of social relationships
has a profound effect on the adopted strategy in maintaining cooperation.


  Many cellular processes are tightly connected to the dynamics of microtubules
(MTs). While in neuronal axons MTs mainly regulate intracellular trafficking,
they participate in cytoskeleton reorganization in many other eukaryotic cells,
enabling the cell to efficiently adapt to changes in the environment. We show
that the functional differences of MTs in different cell types and regions is
reflected in the dynamic properties of MT tips. Using plus-end tracking
proteins EB1 to monitor growing MT plus-ends, we show that MT dynamics and life
cycle in axons of human neurons significantly differ from that of fibroblast
cells. The density of plus-ends, as well as the rescue and catastrophe
frequencies increase while the growth rate decreases toward the fibroblast cell
margin. This results in a rather stable filamentous network structure and
maintains the connection between nucleus and membrane. In contrast, plus-ends
are uniformly distributed along the axons and exhibit diverse polymerization
run times and spatially homogeneous rescue and catastrophe frequencies, leading
to MT segments of various lengths. The probability distributions of the
excursion length of polymerization and the MT length both follow nearly
exponential tails, in agreement with the analytical predictions of a two-state
model of MT dynamics.


  Transmembrane receptor clustering is a ubiquitous phenomenon in pro- and
eukaryotic cells to physically sense receptor/ligand interactions and
subsequently translate an exogenous signal into a cellular response. Despite
that receptor cluster formation has been described for a wide variety of
receptors, ranging from chemotactic receptors in bacteria to growth factor and
neurotransmitter receptors in mammalian cells, a mechanistic understanding of
the underlying molecular processes is still puzzling. In an attempt to fill
this gap we followed a combined experimental and theoretical approach by
dissecting and modulating cargo binding, internalization and cellular response
mediated by KDEL receptors (KDELRs) at the mammalian cell surface after
interaction with a model cargo/ligand. Using a fluorescent variant of ricin
toxin A chain as KDELR-ligand, eGFP-RTA (H/KDEL), we demonstrate that cargo
binding induces dose-dependent receptor cluster formation at and subsequent
internalization from the membrane which is associated and counteracted by
anterograde and microtubule-assisted receptor transport to preferred docking
sites at the plasma membrane. By means of analytical arguments and extensive
numerical simulations we show that cargo-synchronized receptor transport from
and to the membrane is causative for KDELR/cargo cluster formation at the
mammalian cell surface.


  Technological advances have led to a strong increase in the number of data
collection efforts aimed at measuring co-presence of individuals at different
spatial resolutions. It is however unclear how much co-presence data can inform
us on actual face-to-face contacts, of particular interest to study the
structure of a population in social groups or for use in data-driven models of
information or epidemic spreading processes. Here, we address this issue by
leveraging data sets containing high resolution face-to-face contacts as well
as a coarser spatial localisation of individuals, both temporally resolved, in
various contexts. The co-presence and the face-to-face contact temporal
networks share a number of structural and statistical features, but the former
is (by definition) much denser than the latter. We thus consider several
down-sampling methods that generate surrogate contact networks from the
co-presence signal and compare them with the real face-to-face data. We show
that these surrogate networks reproduce some features of the real data but are
only partially able to identify the most central nodes of the face-to-face
network. We then address the issue of using such down-sampled co-presence data
in data-driven simulations of epidemic processes, and in identifying efficient
containment strategies. We show that the performance of the various sampling
methods strongly varies depending on context. We discuss the consequences of
our results with respect to data collection strategies and methodologies.


  Phylogenomics heavily relies on well-curated sequence data sets that consist,
for each gene, exclusively of 1:1-orthologous. Paralogs are treated as a
dangerous nuisance that has to be detected and removed. We show here that this
severe restriction of the data sets is not necessary. Building upon recent
advances in mathematical phylogenetics we demonstrate that gene duplications
convey meaningful phylogenetic information and allow the inference of plausible
phylogenetic trees, provided orthologs and paralogs can be distinguished with a
degree of certainty. Starting from tree-free estimates of orthology, cograph
editing can sufficiently reduce the noise in order to find correct
event-annotated gene trees. The information of gene trees can then directly be
translated into constraints on the species trees. While the resolution is very
poor for individual gene families, we show that genome-wide data sets are
sufficient to generate fully resolved phylogenetic trees, even in the presence
of horizontal gene transfer.
  We demonstrate that the distribution of paralogs in large gene families
contains in itself sufficient phylogenetic signal to infer fully resolved
species phylogenies. This source of phylogenetic information is independent of
information contained in orthologous sequences and is resilient against
horizontal gene transfer. An important consequence is that phylogenomics data
sets need not be restricted to 1:1 orthologs.


  In large scale genetic association studies, a primary aim is to test for
association between genetic variants and a disease outcome. The variants of
interest are often rare, and appear with low frequency among subjects. In this
situation, statistical tests based on standard asymptotic results do not
adequately control the Type I error rate, especially if the case:control ratio
is unbalanced. In this paper, we propose the use of permutation and approximate
unconditional tests for testing association with rare variants. We use novel
analytical calculations to efficiently approximate the true Type I error rate
under common study designs, and in numerical studies show that the proposed
classes of tests significantly improve upon standard testing methods. We also
illustrate our methods in data from a recent case-control study, for genetic
causes of a severe side-effect of a common drug treatment.


  The biological processes of cellular decision making and differentiation
involve a plethora of signalling pathways and gene regulatory circuits. These
networks, in their turn, exhibit a multitude of motifs playing crucial parts in
regulating network activity. Here, we compare the topological placement of
motifs in gene regulatory and signalling networks and find that it suggests
different evolutionary strategies in motif distribution for distinct cellular
subnetworks.


  Allosteric transcription factors undergo binding events both at their inducer
binding sites as well as at distinct DNA binding domains, and it is often
difficult to disentangle the structural and functional consequences of these
two classes of interactions. In this work, we compare the ability of two
statistical mechanical models - the Monod-Wyman-Changeux (MWC) and the
Koshland-N\'emethy-Filmer (KNF) models of protein conformational change - to
characterize the multi-step activation mechanism of the broadly acting
cyclic-AMP receptor protein (CRP). We first consider the allosteric transition
resulting from cyclic-AMP binding to CRP, then analyze how CRP binds to its
operator, and finally investigate the ability of CRP to activate gene
expression. In light of these models, we examine data from a beautiful recent
experiment that created a single-chain version of the CRP homodimer, thereby
enabling each subunit to be mutated separately. Using this construct, six
mutants were created using all possible combinations of the wild type subunit,
a D53H mutant subunit, and an S62F mutant subunit. We demonstrate that both the
MWC and KNF models can explain the behavior of all six mutants using a small,
self-consistent set of parameters. In comparing the results, we find that the
MWC model slightly outperforms the KNF model in the quality of its fits, but
more importantly the parameters inferred by the MWC model are more in line with
structural knowledge of CRP. In addition, we discuss how the conceptual
framework developed here for CRP enables us to not merely analyze data
retrospectively, but has the predictive power to determine how combinations of
mutations will interact, how double mutants will behave, and how each construct
would regulate gene expression.


  Receptor-induced apoptosis is a complex signal transduction pathway involving
numerous protein/protein interactions and post-transcriptional modifications.
The response to death receptor stimulation varies significantly from one cell
line to another and even from one cell to another within a given cell line. In
this context, it is often difficult to assess whether the molecular mechanisms
identified so far are sufficient to explain the rich quantitative observations
now available, and to detect possible gaps in our understanding. This is
precisely where computational systems biology approaches may contribute. In
this chapter, we review studies done in this direction, focusing on those that
provided significant insight on the functioning of this complex pathway by
tightly integrating experimental and computational approaches.


  In plant-pollinator communities many pollinators are potential generalists
and their preferences for certain plants can change quickly in response to
changes in plant and pollinator densities. These changes in preferences affect
coexistence within pollinator guilds as well as within plant guilds. Using a
mathematical model, we study how adaptations of pollinator preferences
influence population dynamics of a two-plant-two-pollinator community
interaction module. Adaptation leads to coexistence between generalist and
specialist pollinators, and produces complex plant population dynamics,
involving alternative stable states and discrete transitions in the plant
community. Pollinator adaptation also leads to plant-plant apparent
facilitation that is mediated by changes in pollinator preferences. We show
that adaptive pollinator behavior reduces niche overlap and leads to
coexistence by specialization on different plants. Thus, this article documents
how adaptive pollinator preferences for plants change the structure and
coexistence of plant-pollinator communities.


  Structure based ligand discovery is one of the most successful approaches for
augmenting the drug discovery process. Currently, there is a notable shift
towards machine learning (ML) methodologies to aid such procedures. Deep
learning has recently gained considerable attention as it allows the model to
"learn" to extract features that are relevant for the task at hand. We have
developed a novel deep neural network estimating the binding affinity of
ligand-receptor complexes. The complex is represented with a 3D grid, and the
model utilizes a 3D convolution to produce a feature map of this
representation, treating the atoms of both proteins and ligands in the same
manner. Our network was tested on the CASF "scoring power" benchmark and Astex
Diverse Set and outperformed classical scoring functions. The model, together
with usage instructions and examples, is available as a git repository at
http://gitlab.com/cheminfIBB/pafnucy


  Background. Protein dihedral angles provide a detailed description of protein
local conformation. Predicted dihedral angles can be used to narrow down the
conformational space of the whole polypeptide chain significantly, thus aiding
protein tertiary structure prediction. However, direct angle prediction from
sequence alone is challenging.
  Method. In this study, we present a novel method to predict real-valued
angles by combining clustering and deep learning. That is, we first generate
certain clusters of angles (each assigned a label) and then apply a deep
residual neural network to predict the label posterior probability. Finally, we
output real-valued prediction by a mixture of the clusters with their predicted
probabilities. At the same time, we also estimate the bound of the prediction
errors at each residue from the predicted label probabilities.
  Result. In this article, we present a novel method (named RaptorX-Angle) to
predict real-valued angles by combining clustering and deep learning. Tested on
a subset of PDB25 and the targets in the latest two Critical Assessment of
protein Structure Prediction (CASP), our method outperforms the existing
state-of-art method SPIDER2 in terms of Pearson Correlation Coefficient (PCC)
and Mean Absolute Error (MAE). Our result also shows approximately linear
relationship between the real prediction errors and our estimated bounds. That
is, the real prediction error can be well approximated by our estimated bounds.
  Conclusions. Our study provides an alternative and more accurate prediction
of dihedral angles, which may facilitate protein structure prediction and
functional study.


  Hyperpolarized 13C-MRI allows real time observation of metabolism in vivo.
Imaging sequences have been developed to follow the metabolism of [1-13C]
pyruvate and extract reaction kinetics, which can show tumour treatment
response. We applied the fitting model and algorithm for the imaging data of
mice tumour models and determined error estimates for the parameters of
interest. Data was least-squares fitted onto a two-site exchange model in
MATLAB, followed by statistic computation to assess model performance.
Inference through the application of MCMC was also performed. The modelling and
inference process extracted quantitative information satisfactorily and
reproducibly, demonstrating metabolic activity and intratumour heterogeneity.
Finally, novel fitting methods were evaluated and further recommendations were
made.


  Rhino populations are at a critical level due to the demand for rhino horn
and the subsequent poaching. Wildlife managers attempt to secure rhinos with
approaches to devalue the horn, the most common of which is dehorning. Game
theory has been used to examine the interaction of poachers and wildlife
managers where a manager can either `dehorn' their rhinos or leave the horn
attached and poachers may behave `selectively' or `indiscriminately'. The
approach described in this paper builds on this previous work and investigates
the interactions between the poachers. We build an evolutionary game theoretic
model and determine which strategy is preferred by a poacher in various
different populations of poachers. The purpose of this work is to discover
whether conditions which encourage the poachers to behave selectively exist,
that is, they only kill those rhinos with full horns.
  The analytical results show that full devaluation of all rhinos will likely
lead to indiscriminate poaching. In turn it shows that devaluing of rhinos can
only be effective when implemented along with a strong disincentive framework.
This paper aims to contribute to the necessary research required for informed
discussion about the lively debate on legalising rhino horn trade.