  Stronger selection implies faster evolution---that is, the greater the force,
the faster the change. This apparently self-evident proposition, however, is
derived under the assumption that genetic variation within a population is
primarily supplied by mutation (i.e.\ mutation-driven evolution). Here, we show
that this proposition does not actually hold for recombination-driven
evolution, i.e.\ evolution in which genetic variation is primarily created by
recombination rather than mutation. By numerically investigating population
genetics models of recombination, migration and selection, we demonstrate that
stronger selection can slow down evolution on a perfectly smooth fitness
landscape. Through simple analytical calculation, this apparently
counter-intuitive result is shown to stem from two opposing effects of natural
selection on the rate of evolution. On the one hand, natural selection tends to
increase the rate of evolution by increasing the fixation probability of fitter
genotypes. On the other hand, natural selection tends to decrease the rate of
evolution by decreasing the chance of recombination between immigrants and
resident individuals. As a consequence of these opposing effects, there is a
finite selection pressure maximizing the rate of evolution. Hence, stronger
selection can imply slower evolution if genetic variation is primarily supplied
by recombination.


  The propagation of a beneficial mutation in a spatially extended population
is usually studied using the phenomenological stochastic Fisher-Kolmogorov
(SFKPP) equation. We derive here an individual based, stochastic model founded
on the spatial Moran process where fluctuations are treated exactly. At high
selection pressure, the results of this model are different from the classical
FKPP. At small selection pressure, the front behavior can be mapped into a
Brownian motion with drift, the properties of which can be derived from
microscopic parameters of the Moran model. Finally, we show that the diffusion
coefficient and the noise amplitude of SFKPP are not independent parameters but
are both determined by the dispersal kernel of individuals.


  Background: It is necessary and essential to discovery protein function from
the novel primary sequences. Wet lab experimental procedures are not only
time-consuming, but also costly, so predicting protein structure and function
reliably based only on amino acid sequence has significant value. TATA-binding
protein (TBP) is a kind of DNA binding protein, which plays a key role in the
transcription regulation. Our study proposed an automatic approach for
identifying TATA-binding proteins efficiently, accurately, and conveniently.
This method would guide for the special protein identification with
computational intelligence strategies. Results: Firstly, we proposed novel
fingerprint features for TBP based on pseudo amino acid composition,
physicochemical properties, and secondary structure. Secondly, hierarchical
features dimensionality reduction strategies were employed to improve the
performance furthermore. Currently, Pretata achieves 92.92% TATA- binding
protein prediction accuracy, which is better than all other existing methods.
Conclusions: The experiments demonstrate that our method could greatly improve
the prediction accuracy and speed, thus allowing large-scale NGS data
prediction to be practical. A web server is developed to facilitate the other
researchers, which can be accessed at http://server.malab.cn/preTata/.


  We first recall some basic facts from the theory of discrete-time Markov
chains arising from two types neutral and non-neutral evolution models of
population genetics with constant size. We then define and analyse a version of
such models whose fluctuating total population size is conserved on average
only. In our model, the population of interest is seen as being embedded in a
frame process which is a critical Galton-Watson process. In this context, we
address problems such as extinction, fixation, size of the population at
fixation and survival probability to a bottleneck effect of the environment.
Running title: constant population size on average


  In a recent paper published in Ecosphere, their authors suggest that
extending the logistic growth model in its usual r - K parametrization to a
multi-patch environment results in undesirable properties, that were referred
to as the "perfect mixing paradox". This led the authors to recommend using the
Verhulst r - {\alpha} parametrization of the logistic model instead and
abandoning the term "carrying capacity" in the context of population dynamics.
In this study we show that the use of the logistic equation in its traditional
r - K parametrization is appropriate for representing the dynamics of
populations in a simple spatial context with passive migration among patches.
Furthermore, we show that the conclusions of the mentioned paper depend on the
specific functional form for migration rates used for their analyses. In
addition, we suggest that their specific migration model was not the best
choice since biologically realistic options exist. Moreover, these alternatives
require the same number of parameters. The model we present here is free of the
paradoxical behaviors presented previously. This allows us to conclude that the
logistic growth in its usual r - K parametrization is useful in a spatial
context if extended appropriately to a multi-patch scenario and consequently
there are no reasons to abandon the concept of carrying capacity. Ecologists
should pay great attention when using models in scenarios that are more complex
or just different from the ones for which the models were created.


  Individual computations and social interactions underlying collective
behavior in groups of animals are of great ethological, behavioral, and
theoretical interest. While complex individual behaviors have successfully been
parsed into small dictionaries of stereotyped behavioral modes, studies of
collective behavior largely ignored these findings; instead, their focus was on
inferring single, mode-independent social interaction rules that reproduced
macroscopic and often qualitative features of group behavior. Here we bring
these two approaches together to predict individual swimming patterns of adult
zebrafish in a group. We show that fish alternate between an active mode in
which they are sensitive to the swimming patterns of conspecifics, and a
passive mode where they ignore them. Using a model that accounts for these two
modes explicitly, we predict behaviors of individual fish with high accuracy,
outperforming previous approaches that assumed a single continuous computation
by individuals and simple metric or topological weighing of neighbors behavior.
At the group level, switching between active and passive modes is uncorrelated
among fish, yet correlated directional swimming behavior still emerges. Our
quantitative approach for studying complex, multi-modal individual behavior
jointly with emergent group behavior is readily extensible to additional
behavioral modes and their neural correlates, as well as to other species.


  Changepoint detection is a central problem in time series and genomic data.
For some applications, it is natural to impose constraints on the directions of
changes. One example is ChIP-seq data, for which adding an up-down constraint
improves peak detection accuracy, but makes the optimization problem more
complicated. We show how a recently proposed functional pruning technique can
be adapted to solve such constrained changepoint detection problems. This leads
to a new algorithm which can solve problems with arbitrary affine constraints
on adjacent segment means, and which has empirical time complexity that is
log-linear in the amount of data. This algorithm achieves state-of-the-art
accuracy in a benchmark of several genomic data sets, and is orders of
magnitude faster than existing algorithms that have similar accuracy. Our
implementation is available as the PeakSegPDPA function in the coseg R package,
https://github.com/tdhock/coseg


  In the immune system, T cells can quickly discriminate between foreign and
self ligands with high accuracy. There is evidence that T-cells achieve this
remarkable performance utilizing a network architecture based on a
generalization of kinetic proofreading (KPR). KPR-based mechanisms actively
consume energy to increase the specificity beyond what is possible in
equilibrium.An important theoretical question that arises is to understand the
trade-offs and fundamental limits on accuracy, speed, and dissipation (energy
consumption) in KPR and its generalization. Here, we revisit this question
through numerical simulations where we simultaneously measure the speed,
accuracy, and energy consumption of the KPR and adaptive sorting networks for
different parameter choices. Our simulations highlight the existence of a
'feasible operating regime' in the speed-energy-accuracy plane where T-cells
can quickly differentiate between foreign and self ligands at reasonable energy
expenditure. We give general arguments for why we expect this feasible
operating regime to be a generic property of all KPR-based biochemical networks
and discuss implications for our understanding of the T cell receptor circuit.


  We investigate the role of greed on the lifetime of a random-walking forager
on an initially resource-rich lattice. Whenever the forager lands on a
food-containing site, all the food there is eaten and the forager can hop
$\mathcal{S}$ more steps without food before starving. Upon reaching an empty
site, the forager comes one time unit closer to starvation. The forager is also
greedy---given a choice to move to an empty or to a food-containing site in its
local neighborhood, the forager moves preferentially towards food.
Surprisingly, the forager lifetime varies non-monotonically with greed, with
different senses of the non-monotonicity in one and two dimensions. Also
unexpectedly, the forager lifetime in one dimension has a huge peak for very
negative greed.


  The day we understand the time evolution of subcellular elements at a level
of detail comparable to physical systems governed by Newton's laws of motion
seems far away. Even so, quantitative approaches to cellular dynamics add to
our understanding of cell biology, providing data-guided frameworks that allow
us to develop better predictions about and methods for control over specific
biological processes and system-wide cell behavior. In this paper we describe
an approach to optimizing the use of transcription factors in the context of
cellular reprogramming. We construct an approximate model for the natural
evolution of a synchronized population of fibroblasts, based on data obtained
by sampling the expression of some 22,083 genes at several times along the cell
cycle. (These data are based on a colony of cells that have been cell cycle
synchronized) In order to arrive at a model of moderate complexity, we cluster
gene expression based on the division of the genome into topologically
associating domains (TADs) and then model the dynamics of the expression levels
of the TADs. Based on this dynamical model and known bioinformatics, we develop
a methodology for identifying the transcription factors that are the most
likely to be effective toward a specific cellular reprogramming task. The
approach used is based on a device commonly used in optimal control. From this
data-guided methodology, we identify a number of validated transcription
factors used in reprogramming and/or natural differentiation. Our findings
highlight the immense potential of dynamical models models, mathematics, and
data guided methodologies for improving methods for control over biological
processes.


  Native horse mucus is characterized with micro- and macrorheology and
compared to hydroxyethylcellulose (HEC) gel as a model. Both systems show
comparable viscoelastic properties on the microscale and for the HEC the
macrorheology is in good agreement with the microrheology. For the mucus, the
viscoelastic moduli on the macroscale are several orders of magnitude larger
than on the microscale. Large amplitude oscillatory shear experiments show that
the mucus responds nonlinearly at much smaller deformations than HEC. This
behavior fosters the assumption that the mucus has a foam like structure on the
microscale compared to the typical mesh like structure of the HEC, a model that
is supported by cryogenic-scanning-electron-microscopy (CSEM) images. These
images allow also to determine the relative amount of volume that is occupied
by the pores and the scaffold. Consequently, we can estimate the elastic
modulus of the scaffold. We conclude that this particular foam like
microstructure should be considered as a key factor for the transport of
particulate matter which plays a central role in mucus function with respect to
particle penetration. The mesh properties composed of very different components
are responsible for macroscopic and microscopic behavior being part of
particles fate after landing.


  Natural protein sequences contain a record of their history. A common
constraint in a given protein family is the ability to fold to specific
structures, and it has been shown possible to infer the main native ensemble by
analyzing covariations in extant sequences. Still, many natural proteins that
fold into the same structural topology show different stabilization energies,
and these are often related to their physiological behavior. We propose a
description for the energetic variation given by sequence modifications in
repeat proteins, systems for which the overall problem is simplified by their
inherent symmetry. We explicitly account for single amino acid and pair-wise
interactions and treat higher order correlations with a single term. We show
that the resulting force field can be interpreted with structural detail. We
trace the variations in the energetic scores of natural proteins and relate
them to their experimental characterization. The resulting energetic force
field allows the prediction of the folding free energy change for several
mutants, and can be used to generate synthetic sequences that are statistically
indistinguishable from the natural counterparts.


  Whether there are gender differences in lying has been largely debated in the
past decade. Previous studies found mixed results. To shed light on this topic,
here I report a meta-analysis of 8,728 distinct observations, collected in 65
Sender-Receiver game treatments, by 14 research groups. Following previous work
and theoretical considerations, I distinguish three types of lies: black lies,
that benefit the liar at a cost for another person; altruistic white lies, that
benefit another person at a cost for the liar; Pareto white lies, that benefit
both the liar and another person. The results show that gender differences in
lying significantly depend on the consequences of lying. Specifically: (i)
males are significantly more likely than females to tell black lies (N=4,161);
(ii) males are significantly more likely than females to tell altruistic white
(N=2,940); (iii) results are inconclusive in the case of Pareto white lies
(N=1,627).


  The Zika virus has been found in individual cases but has not been confirmed
as the cause of in the large number of cases of microcephaly in Brazil in
2015-6. Indeed, disparities between the incidence of Zika and microcephaly
across geographic locations has led to questions about the virus's role. Here
we consider whether the insecticide pyriproxyfen used in Brazilian drinking
water might be the primary cause or a cofactor. Pyriproxifen is a juvenile
hormone analog which has been shown to correspond in mammals to a number of fat
soluble regulatory molecules including retinoic acid, a metabolite of vitamin
A, with which it has cross-reactivity and whose application during development
has been shown to cause microcephaly. Methoprene, another juvenile hormone
analog approved as an insecticide in the 1970s has been shown to cause
developmental disorders in mammals. Isotretinoin is another retinoid causing
microcephaly via activation of the retinoid X receptor in developing fetuses.
We review tests of pyriproxyfen by the manufacturer Sumitomo, which actually
found some evidence for this effect, including low brain mass and
arhinencephaly in exposed rat pups. Pyriproxyfen use in Brazil is
unprecedented, never having been applied to a water supply on a large scale.
Claims that its geographical pattern of use rule it out as a cause have not
been documented or confirmed. On the other hand, the very few microcephaly
cases reported in Colombia and the wide discrepancies of incidence in different
states across Brazil despite large numbers of Zika cases undermine the claim
that Zika is the cause. Given this combination of potential molecular
mechanism, toxicological and epidemiological evidence we strongly recommend
that the use of pyriproxyfen in Brazil be suspended until the potential causal
link to microcephaly is investigated further.


  Although various norms for reciprocity-based cooperation have been suggested
that are evolutionarily stable against invasion from free riders, the process
of alternation of norms and the role of diversified norms remain unclear in the
evolution of cooperation. We clarify the co-evolutionary dynamics of norms and
cooperation in indirect reciprocity and also identify the indispensable norms
for the evolution of cooperation. Inspired by the gene knockout method, a
genetic engineering technique, we developed the norm knockout method and
clarified the norms necessary for the establishment of cooperation. The results
of numerical investigations revealed that the majority of norms gradually
transitioned to tolerant norms after defectors are eliminated by strict norms.
Furthermore, no cooperation emerges when specific norms that are intolerant to
defectors are knocked out.


  An elementary biostatistical theory based on a selectivity-variability
principle is proposed to address a question raised by Charles Darwin, namely,
how one sex of a sexually dimorphic species might tend to evolve with greater
variability than the other sex. Briefly, the theory says that if one sex is
relatively selective then from one generation to the next, more variable
subpopulations of the opposite sex will generally tend to prevail over those
with lesser variability. Moreover, the perhaps less intuitive converse also
holds: if a sex is relatively non-selective, then less variable subpopulations
of the opposite sex will prevail over those with greater variability. This
theory requires certain regularity conditions on the distributions, but makes
no assumptions about differences in means between the sexes, nor does it
presume that one sex is selective and the other non-selective. Two mathematical
models of the selectivity-variability principle are presented: a discrete-time
one-step probabilistic model of short-term behavior with an example using
normally distributed perceived fitness values; and a continuous-time
deterministic model for the long-term asymptotic behavior of the expected sizes
of the subpopulations with an example using exponentially distributed fitness
levels.


  We map a class of well-mixed stochastic models of biochemical feedback in
steady state to the mean-field Ising model near the critical point. The mapping
provides an effective temperature, magnetic field, order parameter, and heat
capacity that can be extracted from biological data without fitting or
knowledge of the underlying molecular details. We demonstrate this procedure on
fluorescence data from mouse T cells, which reveals distinctions between how
the cells respond to different drugs. We also show that the heat capacity
allows inference of absolute molecule number from fluorescence intensity. We
explain this result in terms of the underlying fluctuations and demonstrate the
generality of our work.


  The 1918 influenza pandemic was characterized by multiple epidemic waves. We
investigated into reactive social distancing, a form of behavioral responses,
and its effect on the multiple influenza waves in the United Kingdom. Two forms
of reactive social distancing have been used in previous studies: Power
function, which is a function of the proportion of recent influenza mortality
in a population, and Hill function, which is a function of the actual number of
recent influenza mortality. Using a simple epidemic model with a Power function
and one common set of parameters, we provided a good model fit for the observed
multiple epidemic waves in London boroughs, Birmingham and Liverpool. Our
approach is different from previous studies where separate models are fitted to
each city. We then applied these model parameters obtained from fitting three
cities to all 334 administrative units in England and Wales and including the
population sizes of individual administrative units. We computed the Pearson's
correlation between the observed and simulated data for each administrative
unit. We achieved a median correlation of 0.636, indicating our model
predictions perform reasonably well. Our modelling approach which requires
reduced number of parameters resulted in computational efficiency gain without
over-fitting the model. Our works have both scientific and public health
significance.


  Background: Guillain-Barr\'e Syndrome (GBS) is a common type of severe acute
paralytic neuropathy and associated with other virus infections such as dengue
fever and Zika. This study investigate the relationship between GBS, dengue,
local meteorological factors in Hong Kong and global climatic factors from
January 2000 to June 2016.
  Methods: The correlations between GBS, dengue, Multivariate El Nino Southern
Oscillation Index (MEI) and local meteorological data were explored by the
Spearman Rank correlations and cross-correlations between these time series.
Poisson regression models were fitted to identify nonlinear associations
between MEI and dengue. Cross wavelet analysis was applied to infer potential
non-stationary oscillating associations among MEI, dengue and GBS.
  Findings : An increasing trend was found for both GBS cases and imported
dengue cases in Hong Kong. We found a weak but statistically significant
negative correlation between GBS and local meteorological factors. MEI
explained over 12\% of dengue's variations from Poisson regression models.
Wavelet analyses showed that there is possible non-stationary oscillating
association between dengue and GBS from 2005 to 2015 in Hong Kong. Our study
has led to an improved understanding of the timing and relationship between
GBS, dengue and MEI.


  We discuss the gauge field theory approach to protein structure study, which
allows a natural way to introduce collective degrees of freedom and nonlinear
topological structures. Local symmetry of proteins and its breaking in the
medium is considered, what allows to derive Abelian Higgs model of protein
backbone, correct folding of which is defined by gauge symmetry breaking due
hydrophobic forces. Within this model structure of protein backbone is defined
by superposition of one-dimensional topological solitons (kinks), what allows
to reproduce the three-dimensional structure of the protein backbone with
precision up to 1A and to predict its dynamics.


  Background: Globally, influenza is a major cause of morbidity,
hospitalization and mortality. Influenza vaccination has shown substantial
protective effectiveness in the United States. We investigated state-level
patterns of coverage rates of seasonal and pandemic influenza vaccination,
among the overall population in the U.S. and specifically among children and
the elderly, from 2009/10 to 2014/15, and associations with ecological factors.
Methods and Findings: We obtained state-level influenza vaccination coverage
rates from national surveys, and state-level socio-demographic and health data
from a variety of sources. We employed a retrospective ecological study design,
and used mixed-model regression to determine the levels of ecological
association of the state-level vaccinations rates with these factors, both with
and without region as a factor for the three populations. We found that
health-care access is positively and significantly associated with mean
influenza vaccination coverage rates across all populations and models. We also
found that prevalence of asthma in adults are negatively and significantly
associated with mean influenza vaccination coverage rates in the elderly
populations. Conclusions: Health-care access has a robust, positive association
with state-level vaccination rates across different populations. This
highlights a potential population-level advantage of expanding health-care
access.


  We develop a multi-patch and multi-group model that captures the dynamics of
an infectious disease when the host is structured into an arbitrary number of
groups and interacts into an arbitrary number of patches where the infection
takes place. In this framework, we model host mobility that depends on its
epidemiological status, by a Lagrangian approach. This framework is applied to
a general SEIRS model and the basic reproduction number $\mathcal{R_0}$ is
derived. The effects of heterogeneity in groups, patches and mobility patterns
on $\mathcal{R_0}$ and disease prevalence are explored. Our results show that
for a fixed number of groups, the basic reproduction number increases with
respect to the number of patches and the host mobility patterns. Moreover, when
the mobility matrix of susceptible individuals is of rank one, the basic
reproduction number is explicitly determined and was found to be independent of
the latter. The cases where mobility matrices are of rank one capture important
modeling scenarios. Additionally, we study the global analysis of equilibria
for some special cases. Numerical simulations are carried out to showcase the
ramifications of mobility pattern matrices on disease prevalence and basic
reproduction number.


  Variation in behavior among group members often impacts collective outcomes.
Individuals may vary both in the task that they perform and in the persistence
with which they perform each task. Although both the distribution of
individuals among tasks and differences among individuals in behavioral
persistence can each impact collective behavior, we do not know if and how they
jointly affect collective outcomes. Here we use a detailed computational model
to examine the joint impact of colony-level distribution among tasks and
behavioral persistence of individuals, specifically their fidelity to
particular resource sites, on the collective tradeoff between exploring for new
resources and exploiting familiar ones. We developed an agent-based model of
foraging honey bees, parameterized by data from 5 colonies, in which we
simulated scouts, who search the environment for new resources, and individuals
who are recruited by the scouts to the newly found resources, i.e., recruits.
We found that for each value of persistence there is a different optimal ratio
of scouts to recruits that maximizes resource collection by the colony.
Furthermore, changes to the persistence of scouts induced opposite effects from
changes to the persistence of recruits on the collective foraging of the
colony. The proportion of scouts that resulted in the most resources collected
by the colony decreased as the persistence of recruits increased. However, this
optimal proportion of scouts increased as the persistence of scouts increased.
Thus, behavioral persistence and task participation can interact to impact a
colony's collective behavior in orthogonal directions. Our work provides new
insights and generates new hypotheses into how variation in behavior at both
the individual and colony levels jointly impact the trade-off between exploring
for new resources and exploiting familiar ones.


  Phylogenetic inference aims to reconstruct the evolutionary relationships of
different species based on genetic (or other) data. Discrete characters are a
particular type of data, which contain information on how the species should be
grouped together. However, it has long been known that some characters contain
more information than others. For instance, a character that assigns the same
state to each species groups all of them together and so provides no insight
into the relationships of the species considered. At the other extreme, a
character that assigns a different state to each species also conveys no
phylogenetic signal. In this manuscript, we study a natural combinatorial
measure of the information content of an individual character and analyse
properties of characters that provide the maximum phylogenetic information,
particularly, the number of states such a character uses and how the different
states have to be distributed among the species or taxa of the phylogenetic
tree.


  Uncovering the mechanisms that control size, growth, and division rates of
systems reproducing through binary division means understanding basic
principles of their life cycle. Recent work has focused on how division rates
are regulated in bacteria and yeast, but this question has not yet been
addressed in more complex, multicellular organisms. We have acquired a unique
large-scale data set on the growth and asexual reproduction of two freshwater
planarian species, Dugesia japonica and Dugesia tigrina, which reproduce by
transverse fission and succeeding regeneration of head and tail pieces into new
planarians. We show that generation-dependent memory effects in planarian
reproduction need to be taken into account to accurately capture the
experimental data. To achieve this, we developed a new additive model that
mixes multiple size control strategies based on planarian size, growth, and
time between divisions. Our model quantifies the proportions of each strategy
in the mixed dynamics, revealing the ability of the two planarian species to
utilize different strategies in a coordinated manner for size control.
Additionally, we found that head and tail offspring of both species employ
different mechanisms to monitor and trigger their reproduction cycles. Thus, we
find a diversity of strategies not only between species but also between heads
and tails within species. Our additive model provides two advantages over
existing 2D models that fit a multivariable splitting rate function to the data
for size control: Firstly, it can be fit to relatively small data sets and can
thus be applied to systems where available data is limited. Secondly, it
enables new biological insights because it explicitly shows the contributions
of different size control strategies for each offspring type.


  We study the persistence and extinction of species in a simple food chain
that is modelled by a Lotka-Volterra system with environmental stochasticity.
There exist sharp results for deterministic Lotka-Volterra systems in the
literature but few for their stochastic counterparts. The food chain we analyze
consists of one prey and $n-1$ predators. The $j$th predator eats the $j-1$th
species and is eaten by the $j+1$th predator; this way each species only
interacts with at most two other species - the ones that are immediately above
or below it in the trophic chain. We show that one can classify, based on an
explicit quantity depending on the interaction coefficients of the system,
which species go extinct and which converge to their unique invariant
probability measure. Our work can be seen as a natural extension of the
deterministic results of Gard and Hallam '79 to a stochastic setting.
  As one consequence we show that environmental stochasticity makes species
more likely to go extinct. However, if the environmental fluctuations are
small, persistence in the deterministic setting is preserved in the stochastic
system. Our analysis also shows that the addition of a new apex predator makes,
as expected, the different species more prone to extinction.
  Another novelty of our analysis is the fact that we can describe the behavior
the system when the noise is degenerate. This is relevant because of the
possibility of strong correlations between the effects of the environment on
the different species.


  Ribbons are topological objects of biological and technological importance.
Here, we study the folding of thick ribbons with hydrophobic surfaces in a bad
solvent in regimes in which either the ribbon's thickness or the solvent
molecule size is not vanishingly small compared to the ribbon's width.
Extensive Monte Carlo simulations show that ribbons of various lengths and with
a small stiffness adopt several distinct configurations as the ground state
that include rolled (Archimedean spiral), curled, twisted and globule
conformations. Analytic and numerical calculations based on the consideration
of putative ground states lead to phase diagrams that qualitatively agree with
the simulation results. A symmetry breaking of the planar rolled configuration
in favor of the elongated twisted and the globular ribbons is observed on
increasing the solvent size. Interestingly, the twisted ribbon is found as the
ground state in the absence of any energetic preference for twisting. We show
that the twist of the DNA double helix structure can be stabilized when modeled
as a hydrophobic thick ribbon even in the limit of vanishing solvent size.


  We present an analysis of an avian flu model that yields insight into the
role of different transmission routes in the recurrence of avian influenza
epidemics. Recent modelling work on avian influenza in wild bird populations
takes into account demographic stochasticity and highlights the importance of
environmental transmission in determining the outbreak periodicity, but only
for a weak between-host transmission rate. We determine the relative
contribution of environmental and direct transmission routes to the intensity
of outbreaks. We use an approximation method to simulate noise sustained
oscillations in a stochastic avian flu model with environmental and direct
transmission routes. We see that the oscillations are governed by the product
of a rotation and a slowly varying standard Ornstein-Uhlenbeck process (i.e.,
mean-reverting process). The intrinsic frequency of the damped deterministic
version of the system predicts the dominant period of outbreaks. We show, using
analytic computation of the intrinsic frequency and theoretical power spectral
density, that the outbreak periodicity can be explained in terms of either or
both types of transmission.The amplitude of outbreaks tends to be high when
both types of transmission are strong.


  Avalanches of electrochemical activity in brain networks have been
empirically reported to obey scale-invariant behavior --characterized by
power-law distributions up to some upper cut-off-- both in vitro and in vivo.
Elucidating whether such scaling laws stem from the underlying neural dynamics
operating at the edge of a phase transition is a fascinating possibility, as
systems poised at criticality have been argued to exhibit a number of important
functional advantages. Here we employ a well-known model for neural dynamics
with synaptic plasticity, to elucidate an alternative scenario in which
neuronal avalanches can coexist, overlapping in time, but still remaining
scale-free. Remarkably their scale-invariance does not stem from underlying
criticality nor self-organization at the edge of a continuous phase transition.
Instead, it emerges from the fact that perturbations to the system exhibit a
neutral drift --guided by demographic fluctuations-- with respect to endogenous
spontaneous activity. Such a neutral dynamics --similar to the one in neutral
theories of population genetics-- implies marginal propagation of activity,
characterized by power-law distributed causal avalanches. Importantly, our
results underline the importance of considering causal information --on which
neuron triggers the firing of which-- to properly estimate the statistics of
avalanches of neural activity. We discuss the implications of these findings
both in modeling and to elucidate experimental observations, as well as its
possible consequences for actual neural dynamics and information processing in
actual neural networks.


  In evolutionary biology, phylogenetic networks are constructed to represent
the evolution of species in which reticulate events are thought to have
occurred, such as recombination and hybridization. It is therefore useful to
have efficiently computable metrics with which to systematically compare such
networks. Through developing an optimal algorithm to enumerate all trinets
displayed by a level-1 network (a type of network that is slightly more general
than an evolutionary tree), here we propose a cubic-time algorithm to compute
the trinet distance between two level-1 networks. Employing simulations, we
also present a comparison between the trinet metric and the so-called
Robinson-Foulds phylogenetic network metric restricted to level-1 networks. The
algorithms described in this paper have been implemented in JAVA and are freely
available at https://www.uea.ac.uk/computing/TriLoNet.


  Variations in mitochondrial genes are usually considered to infer
phylogenies. However some of these genes are lesser constraint than other ones,
and thus may blur the phylogenetic signals shared by the majority of the
mitochondrial DNA sequences. To investigate such effects, in this research
work, the molecular phylogeny of the genus Taenia is studied using 14 coding
sequences extracted from mitochondrial genomes of 17 species. We constructed
16,384 trees, using a combination of 1 up to 14 genes. We obtained 131
topologies, and we showed that only four particular instances were relevant.
Using further statistical investigations, we then extracted a particular
topology, which displays more robustness properties.


  Robustness of spatial pattern against perturbations is an indispensable
property of developmental processes for organisms, which need to adapt to
changing environments. Although specific mechanisms for this robustness have
been extensively investigated, little is known about a general mechanism for
achieving robustness in reaction-diffusion systems. Here, we propose a buffered
reaction-diffusion system, in which active states of chemicals mediated by
buffer molecules contribute to reactions, and demonstrate that robustness of
the pattern wavelength is achieved by the dynamics of the buffer molecule. This
robustness is analytically explained as a result of the scaling properties of
the buffered system, which also lead to a reciprocal relationship between the
wavelength's robustness and the plasticity of the spatial phase upon external
perturbations. Finally, we explore the relevance of this reciprocity to
biological systems.


  We consider a recent coinfection model for Tuberculosis (TB), Human
Immunodeficiency Virus (HIV) infection and Acquired Immunodeficiency Syndrome
(AIDS) proposed in [Discrete Contin. Dyn. Syst. 35 (2015), no. 9, 4639--4663].
We introduce and analyze a multiobjective formulation of an optimal control
problem, where the two conflicting objectives are: minimization of the number
of HIV infected individuals with AIDS clinical symptoms and coinfected with
AIDS and active TB; and costs related to prevention and treatment of HIV and/or
TB measures. The proposed approach eliminates some limitations of previous
works. The results of the numerical study provide comprehensive insights about
the optimal treatment policies and the population dynamics resulting from their
implementation. Some nonintuitive conclusions are drawn. Overall, the
simulation results demonstrate the usefulness and validity of the proposed
approach.


  We investigate the formation of stable ecological networks where many species
share the same resource. We show that such stable ecosystem naturally occurs as
a result of extinctions. We obtain an analytical relation for the number of
coexisting species and find a relation describing how many species that may go
extinct as a result of a sharp environmental change. We introduce a special
parameter that is a combination of species traits and resource characteristics
used in the model formulation. This parameter describes the pressure on system
to converge, by extinctions. When that stress parameter is large we obtain that
the species traits concentrate at some values. This stress parameter is thereby
a parameter that determines the level of final biodiversity of the system.
Moreover, we show that dynamics of this limit system can be described by simple
differential equations.


  Finding influential spreaders of information and disease in networks is an
important theoretical problem, and one of considerable recent interest. It has
been almost exclusively formulated as a node-ranking problem -- methods for
identifying influential spreaders rank nodes according to how influential they
are. In this work, we show that the ranking approach does not necessarily work:
the set of most influential nodes depends on the number of nodes in the set.
Therefore, the set of $n$ most important nodes to vaccinate does not need to
have any node in common with the set of $n+1$ most important nodes. We propose
a method for quantifying the extent and impact of this phenomenon, and show
that it is common in both empirical and model networks.


  The complicated, evolving landscape of cancer mutations poses a formidable
challenge to identify cancer genes among the large lists of mutations typically
generated in NGS experiments. The ability to prioritize these variants is
therefore of paramount importance. To address this issue we developed
OncoScore, a text-mining tool that ranks genes according to their association
with cancer, based on available biomedical literature. Receiver operating
characteristic curve and the area under the curve (AUC) metrics on manually
curated datasets confirmed the excellent discriminating capability of OncoScore
(OncoScore cut-off threshold = 21.09; AUC = 90.3%, 95% CI: 88.1-92.5%),
indicating that OncoScore provides useful results in cases where an efficient
prioritization of cancer-associated genes is needed.


  We study a model of seed dispersal that considers the inclusion of an animal
disperser moving diffusively, feeding on fruits and transporting the seeds,
which are later deposited and capable of germination. The dynamics depends on
several population parameters of growth, decay, harvesting, transport,
digestion and germination. In particular, the deposition of transported seeds
at places away from their collection sites produces a delay in the dynamics,
whose effects are the focus of this work. Analytical and numerical solutions of
different simplified scenarios show the existence of travelling waves. The
effect of zoochory is apparent in the increase of the velocity of these waves.
The results support the hypothesis of the relevance of animal mediated seed
dispersion when trying to understand the origin of the high rates of vegetable
invasion observed in real systems.


  To the known fact that Parsimony method sometimes fails on the problem of
inferring species trees from gene trees, here we proved that no mater of what
topology the true 9-taxon and greater species tree is the only thing one needs
to break down Parsimony is to have in this species tree three consecutive inner
edges not going through the root but perhaps ending on it with lengths T1, T2,
T3 of some proportions. Obviously, the probability to meet these lengths is
growing in general with the size of species tree. Therefore, Parsimony may be
applied only when the described lengths of edges cannot be met in the tree.


  Herein (the first part of my work), I debunk the long-standing hypotheses
that explain mitochondrial oxidative phosphorylation. Simple calculations point
out that mitochondria are highly proton-deficient microcosms and therefore,
elaborate proton pump machinery are not tenable. Further, other elements like
the elaborate electron transport chain, chemiosmosis, rotary ATP synthesis,
etc. are also critically evaluated to point out that such complicated systems
are non-viable. The communication necessitates a new explanatory paradigm for
cellular respiration. In the second part of my work, I have put forward a
viable alternative explanatory paradigm for mitochondrial oxidative
phosphorylation.


  Via a concomitant communication (the first part of my work), I have
conclusively debunked the prevailing explanations for mitochondrial oxidative
phosphorylation and established the need for a novel rationale to account for
the reaction paradigm. Towards the same, murburn concept is hereby floated as a
viable explanation (in the second part of my work). It is proposed that the
inner mitochondrial membrane (harboring the various metal and flavin enzyme
complexes) serves as means to confine and stabilize radical reactions, which
effectively couple and bring about ATP synthesis in the proton-deficient
microcosm. The proposed scheme is un-ordered and favored by Ockham's razor and
evolutionary perspectives. Murburn concept is a paradigm-shift in biochemistry
because it advocates that diffusible reactive (oxygen) species are mainstay of
routine cellular metabolic process within the mitochondria.


  A growing body of empirical evidence indicates that social and cooperative
behavior can be affected by cognitive and neurological factors, suggesting the
existence of state-based decision-making mechanisms that may have emerged by
evolution. Motivated by these observations, we propose a simple mechanism of
anonymous network interactions identified as a form of generalized reciprocity
- a concept organized around the premise "help anyone if helped by someone",
and study its dynamics on random graphs. In the presence of such mechanism, the
evolution of cooperation is related to the dynamics of the levels of
investments (i.e. probabilities of cooperation) of the individual nodes
engaging in interactions. We demonstrate that the propensity for cooperation is
determined by a network centrality measure here referred to as neighborhood
importance index and discuss relevant implications to natural and artificial
systems. To address the robustness of the state-based strategies to an invasion
of defectors, we additionally provide an analysis which redefines the results
for the case when a fraction of the nodes behave as unconditional defectors.


  Chemotherapy is a class of cancer treatment that uses drugs to kill cancer
cells. A typical chemotherapeutic protocol consists of several drugs delivered
in cycles of three weeks. We present mathematical analyses demonstrating the
existence of a maximum time between cycles of chemotherapy for a protocol to be
effective. A mathematical equation is derived, which relates such a maximum
time with the variables that govern the kinetics of the tumor and those
characterizing the chemotherapeutic treatment. Our results suggest that there
are compelling arguments supporting the use of dose-dense protocols. Finally,
we discuss the limitations of these protocols and suggest an alternative.


  Relations among species in ecosystems can be represented by complex networks
where both negative (competition) and positive (mutualism) interactions are
concurrently present. Recently, it has been shown that many ecosystems can be
cast into mutualistic networks, and that nestedness reduces effective
inter-species competition, thus facilitating mutually beneficial interactions
and increasing the number of coexisting species or the biodiversity. However,
current approaches neglect the structure of inter-species competition by
adopting a mean-field perspective that does not deal with competitive
interactions properly. Here, we introduce a framework based on the concept of
multilayer networks, which naturally accounts for both mutualism and
competition. Hence, we abandon the mean field hypothesis and show, through a
dynamical population model and numerical simulations, that there is an
intricate relation between competition and mutualism. Specifically, we show
that when all interactions are taken into account, mutualism does not have the
same consequences on the evolution of specialist and generalist species. This
leads to a non-trivial profile of biodiversity in the parameter space of
competition and mutualism. Our findings emphasize how the simultaneous
consideration of positive and negative interactions can contribute to our
understanding of the delicate trade-offs between topology and biodiversity in
ecosystems and call for a reconsideration of previous findings in theoretical
ecology, as they may affect the structural and dynamical stability of
mutualistic systems.


  Pre-exposure prophylaxis (PrEP) consists in the use of an antiretroviral
medication to prevent the acquisition of HIV infection by uninfected
individuals and has recently demonstrated to be highly efficacious for HIV
prevention. We propose a new epidemiological model for HIV/AIDS transmission
including PrEP. Existence, uniqueness and global stability of the disease free
and endemic equilibriums are proved. The model with no PrEP is calibrated with
the cumulative cases of infection by HIV and AIDS reported in Cape Verde from
1987 to 2014, showing that it predicts well such reality. An optimal control
problem with a mixed state control constraint is then proposed and analyzed,
where the control function represents the PrEP strategy and the mixed
constraint models the fact that, due to PrEP costs, epidemic context and
program coverage, the number of individuals under PrEP is limited at each
instant of time. The objective is to determine the PrEP strategy that satisfies
the mixed state control constraint and minimizes the number of individuals with
pre-AIDS HIV-infection as well as the costs associated with PrEP. The optimal
control problem is studied analytically. Through numerical simulations, we
demonstrate that PrEP reduces HIV transmission significantly.


  An Elementary Flux Mode (EFM) is a pathway with minimum set of reactions that
are functional in steady-state constrained space. Due to the high computational
complexity of calculating EFMs, different approaches have been proposed to find
these flux-balanced pathways. In this paper, an approach to find a subset of
EFMs is proposed based on a graph data model. The given metabolic network is
mapped to the graph model and decisions for reaction inclusion can be made
based on metabolites and their associated reactions. This notion makes the
approach more convenient to categorize the output pathways. Implications of the
proposed method on metabolic networks are discussed.


  Understanding the relationship between spontaneous stochastic fluctuations
and the topology of the underlying gene regulatory network is of fundamental
importance for the study of single-cell stochastic gene expression. Here by
solving the analytical steady-state distribution of the protein copy number in
a general kinetic model of stochastic gene expression with nonlinear feedback
regulation, we reveal the relationship between stochastic fluctuations and
feedback topology at the single-molecule level, which provides novel insights
into how and to what extent a feedback loop can enhance or suppress molecular
fluctuations. Based on such relationship, we also develop an effective method
to extract the topological information of a gene regulatory network from
single-cell gene expression data. The theory is demonstrated by numerical
simulations and, more importantly, validated quantitatively by single-cell data
analysis of a synthetic gene circuit integrated in human kidney cells.


  Consumer demand for plant and animal products threatens many populations with
extinction. The anthropogenic Allee effect (AAE) proposes that such extinctions
can be caused by prices for wildlife products increasing with species rarity.
This price-rarity relationship creates financial incentives to extract the last
remaining individuals of a population, despite higher search and harvest costs.
The AAE has become a standard approach for conceptualizing the threat of
economic markets on endangered species. Despite its potential importance for
conservation, AAE theory is based on a simple graphical model with limited
analysis of possible population trajectories. By specifying a general class of
functions for price-rarity relationships, we show that the classic theory can
understate the risk of species extinction. AAE theory proposes that only
populations below a critical Allee threshold will go extinct due to increasing
price-rarity relationships. Our analysis shows that this threshold can be much
higher than the original theory suggests, depending on initial harvest effort.
More alarmingly, even species with population sizes above this Allee threshold,
for which AAE predicts persistence, can be destined to extinction. Introducing
even a minimum price for harvested individuals, close to zero, can cause large
populations to cross the classic anthropogenic Allee threshold on a trajectory
towards extinction. These results suggest that traditional AAE theory may give
a false sense of security when managing large harvested populations.


  Infectious disease outbreaks recapitulate biology: they emerge from the
multi-level interaction of hosts, pathogens, and their shared environment. As a
result, predicting when, where, and how far diseases will spread requires a
complex systems approach to modeling. Recent studies have demonstrated that
predicting different components of outbreaks--e.g., the expected number of
cases, pace and tempo of cases needing treatment, demand for prophylactic
equipment, importation probability etc.--is feasible. Therefore, advancing both
the science and practice of disease forecasting now requires testing for the
presence of fundamental limits to outbreak prediction. To investigate the
question of outbreak prediction, we study the information theoretic limits to
forecasting across a broad set of infectious diseases using permutation entropy
as a model independent measure of predictability. Studying the predictability
of a diverse collection of historical outbreaks--including, chlamydia, dengue,
gonorrhea, hepatitis A, influenza, measles, mumps, polio, and whooping
cough--we identify a fundamental entropy barrier for infectious disease time
series forecasting. However, we find that for most diseases this barrier to
prediction is often well beyond the time scale of single outbreaks. We also
find that the forecast horizon varies by disease and demonstrate that both
shifting model structures and social network heterogeneity are the most likely
mechanisms for the observed differences across contagions. Our results
highlight the importance of moving beyond time series forecasting, by embracing
dynamic modeling approaches, and suggest challenges for performing model
selection across long time series. We further anticipate that our findings will
contribute to the rapidly growing field of epidemiological forecasting and may
relate more broadly to the predictability of complex adaptive systems.


  We develope a two-species exclusion process with a distinct pair of entry and
exit sites for each species of rigid rods. The relatively slower forward
stepping of the rods in an extended bottleneck region, located in between the
two entry sites, controls the extent of interference of the co-directional flow
of the two species of rods. The relative positions of the sites of entry of the
two species of rods with respect to the location of the bottleneck are
motivated by a biological phenomenon. However, the primary focus of the study
here is to explore the effects of the interference of the flow of the two
species of rods on their spatio-temporal organization and the regulations of
this interference by the extended bottleneck. By a combination of mean-field
theory and computer simulation we calculate the flux of both species of rods
and their density profiles as well as the composite phase diagrams of the
system. If the bottleneck is sufficiently stringent some of the phases become
practically unrealizable although not ruled out on the basis of any fundamental
physical principle. Moreover the extent of suppression of flow of the
downstream entrants by the flow of the upstream entrants can also be regulated
by the strength of the bottleneck. We speculate on the possible implications of
the results in the context of the biological phenomenon that motivated the
formulation of the theoretical model.


  The diversity revealed by large scale genomics in microbiology is calling
into question long held beliefs about genome stability, evolutionary rate, even
the definition of a species. MacArthur and Wilson's theory of insular
biogeography provides an explanation for the diversity of macroscopic animal
and plant species as a consequence of the associated hierarchical web of
species interdependence. We report a large scale study of microbial diversity
that reveals that the cumulative number of genes discovered increases with the
number of genomes studied as a simple power law. This result is demonstrated
for three different genera comparing over 15,000 isolates. We show that this
power law is formally related to the MacArthur-Wilson exponent, suggesting the
emerging diversity of microbial genotypes arises because the scale independent
behavior first reported by MacArthur and Wilson extends down to the scale of
microbes and their genes. Assessing the depth of available whole genome
sequences implies a dynamically changing core genome, suggesting that
traditional taxonomic classifications should be replaced with a quasispecies
model that captures the diversity and dynamic exchange of genes. We report
Species population "clouds" in a defined microbiome, with scale invariance
extending down to the level of single-nucleotide polymorphisms (SNPs).


  We describe and analyse Levenberg-Marquardt methods for solving systems of
nonlinear equations. More specifically, we propose an adaptive formula for the
Levenberg-Marquardt parameter and analyse the local convergence of the method
under H\"{o}lder metric subregularity of the function defining the equation and
H\"older continuity of its gradient mapping. Further, we analyse the local
convergence of the method under the additional assumption that the
\L{}ojasiewicz gradient inequality holds. We finally report encouraging
numerical results confirming the theoretical findings for the problem of
computing moiety conserved steady states in biochemical reaction networks. This
problem can be cast as finding a solution of a system of nonlinear equations,
where the associated mapping satisfies the \L{}ojasiewicz gradient inequality
assumption.


  The plant hormones brassinosteroid (BR) and gibberellin (GA) have important
roles in a wide range of processes involved in plant growth and development. In
this paper we derive and analyse new mathematical models for the BR signalling
pathway and for the crosstalk between the BR and GA signalling pathways. To
analyse the effects of spatial heterogeneity of the signalling processes, along
with spatially-homogeneous ODE models we derive coupled PDE-ODE systems
modelling the temporal and spatial dynamics of molecules involved in the
signalling pathways. The values of the parameters in the model for the BR
signalling pathway are determined using experimental data on the gene
expression of BR biosynthetic enzymes. The stability of steady state solutions
of our mathematical model, shown for a wide range of parameters, can be related
to the BR homeostasis. Our results for the crosstalk model suggest that the
interaction between transcription factors BZR and DELLA exerts more influence
on the dynamics of the signalling pathways than BZR-mediated biosynthesis of
GA, suggesting that the interaction between transcription factors may
constitute the principal mechanism of the crosstalk between the BR and GA
signalling pathways. In general, perturbations in the GA signalling pathway
have larger effects on the dynamics of components of the BR signalling pathway
than perturbations in the BR signalling pathway on the dynamics of the GA
pathway. The perturbation in the crosstalk mechanism also has a larger effect
on the dynamics of the BR pathway than of the GA pathway. Large changes in the
dynamics of the GA signalling processes can be observed only in the case where
there are disturbances in both the BR signalling pathway and the crosstalk
mechanism.


  During asymmetric cell division, dynein generates forces, which position the
spindle to reflect polarity and ensure correct daughter cell fates. The
transient cortical localization of dynein raises the question of its targeting.
We found that it accumulates at the microtubule plus-ends like in budding
yeast, indirectly hitch-hiking on $\text{EBP-2}^{\text{EB1}}$ likely via
dynactin. Importantly, this mechanism, which modestly accounts for cortical
forces, does not transport dynein, which displays the same binding/unbinding
dynamics as $\text{EBP-2}^{\text{EB1}}$. At the cortex, dynein tracks can be
classified as having either directed or diffusive-like motion. Diffusive-like
tracks reveal force-generating dyneins. Their densities are higher on the
posterior tip of the embryos, where $\text{GPR-1/2}^{\text{LGN}}$ concentrate,
but their durations are symmetric. Since dynein flows to the cortex are
non-polarized, we suggest that this posterior enrichment increases dynein
binding, thus accounts for the force imbalance reflecting polarity, and
supplements the regulation of mitotic progression via the non-polarized
detachment rate.


  We here present SIMLR (Single-cell Interpretation via Multi-kernel LeaRning),
an open-source tool that implements a novel framework to learn a
sample-to-sample similarity measure from expression data observed for
heterogenous samples. SIMLR can be effectively used to perform tasks such as
dimension reduction, clustering, and visualization of heterogeneous populations
of samples. SIMLR was benchmarked against state-of-the-art methods for these
three tasks on several public datasets, showing it to be scalable and capable
of greatly improving clustering performance, as well as providing valuable
insights by making the data more interpretable via better a visualization.
Availability and Implementation
  SIMLR is available on GitHub in both R and MATLAB implementations.
Furthermore, it is also available as an R package on http://bioconductor.org.


  Background: Mustard aphid is a major pest of Brassica oilseeds. No source for
aphid resistance is presently available in Brassica juncea . A wild crucifer,
Brassica fruticulosa is known to be resistant to mustard aphid. An artificially
synthesized amphiploid, AD-4 (B. fruticulosa x B. rapa var. brown sarson) was
developed for use as a bridge species to transfer fruticulosa resistance to B.
juncea. Using the selfed backcross we could select a large number of lines with
resistance to mustard aphid. This paper reports cytogenetic stability of
introgression lines, molecular evidence for alien introgression and their
reaction to mustard aphid infestation. Results: Majority of introgression lines
had expected euploid chromosome number(2n= 36), showed normal meiosis and high
pollen grain fertility. Well-distributed and transferable simple-sequence
repeats (SSR) markers for all the 18 B. juncea chromosomes helped to
characterize introgression events. Average proportions of recipient and donor
genome in the substitution lines were 49.72 and 35.06%, respectively. Minimum
alien parent genome presence (27.29%) was observed in the introgression line,
Ad3K-280 . Introgressed genotypes also varied for their resistance responses to
mustard aphid infestations under artificial release conditions for two
continuous seasons. Some of the test genotypes showed consistent resistant
reaction. Conclusions: B.juncea-fruticulosa introgression set may prove to be a
very powerful breeding tool for aphid resistance related QTL/gene discovery and
fine mapping of the desired genes/QTLs to facilitate marker assisted transfer
of identified gene(s) for mustard aphid resistance in the background of
commercial mustard genotypes.


  A six year field study was conducted from 2001 2002 to 2006 2007 at Punjab
Agricultural University, Ludhiana, India to study the losses in seed yield of
different Brassica species (B. juncea, B. napus, B. carinata, B. rapa and Eruca
sativa) by the infestation of insect pests. The experiment was conducted in two
different sets viz. protected/sprayed and unprotected, in a randomized block
design, with three replications. Data on the infestation of insect pests, and
seed yield were recorded at weekly intervals and at harvest, respectively. The
loss in seed yield, due to mustard aphid and cabbage caterpillar, varied from
6.5 to 26.4 per cent. E. sativa suffered the least loss in seed yield and
harboured the minimum population of mustard aphid (2.1 aphids/plant) and
cabbage caterpillar (2.4 larvae/plant). On the other hand, B. carinata was
highly susceptible to the cabbage caterpillar (26.2 larvae/plant) and suffered
the maximum yield loss (26.4%).


  There is a long-standing experimental observation that the melting of
topologically constrained DNA, such as circular-closed plasmids, is less abrupt
than that of linear molecules. This finding points to an intriguing role of
topology in the physics of DNA denaturation, which is however poorly
understood. Here, we shed light on this issue by combining large-scale Brownian
Dynamics simulations with an analytically solvable phenomenological Landau mean
field theory. We find that the competition between melting and supercoiling
leads to phase coexistence of denatured and intact phases at the single
molecule level. This coexistence occurs in a wide temperature range, thereby
accounting for the broadening of the transition. Finally, our simulations show
an intriguing topology-dependent scaling law governing the growth of
denaturation bubbles in supercoiled plasmids, which can be understood within
the proposed mean field theory.


  Essentially all biology is active and dynamic. Biological entities
autonomously sense, com- pute, and respond using energy-coupled ratchets that
can produce force and do work. The cytoskeleton, along with its associated
proteins and motors, is a canonical example of biological active matter, which
is responsible for cargo transport, cell motility, division, and morphol- ogy.
Prior work on cytoskeletal active matter systems showed either extensile or
contractile dynamics. Here, we demonstrate a cytoskeletal system that can
control the direction of the network dynamics to be either extensile,
contractile, or static depending on the concentration of filaments or transient
crosslinkers through systematic variation of the crosslinker or micro- tubule
concentrations. Based off these new observations and our previously published
results, we created a simple one-dimensional model of the interaction of
filaments within a bundle. Despite its simplicity, our model recapitulates the
observed activities of our experimental sys- tem, implying that the dynamics of
our finite networks of bundles are driven by the local filament-filament
interactions within the bundle. Finally, we show that contractile phases can
result in autonomously motile networks that resemble cells. Our experiments and
model allow us to gain a deeper understanding of cytoskeletal dynamics and
provide a stepping stone for designing active, autonomous systems that could
potentially dynamically switch states.


  Alignment of large genomic sequences is a fundamental task in computational
genome analysis. Most methods for genomic alignment use high-scoring local
alignments as {\em anchor points} to reduce the search space of the alignment
procedure. Speed and quality of these methods therefore depend on the
underlying anchor points. Herein, we propose to use {\em Filtered Spaced Word
Matches} to calculate anchor points for genome alignment. To evaluate this
approach, we used these anchor points in the the widely used alignment pipeline
{\em Mugsy}. For distantly related sequence sets, we could substantially
improve the quality of alignments produced by {\em Mugsy}.


  We introduce a new dominance concept consisting of three new dominance
metrics based on Lloyd's (1967) mean crowding index. The new metrics link
communities and species, whereas existing ones are applicable only to
communities. Our community-level metric is a function of Simpson's diversity
index. For species, our metric quantifies the difference between community
dominance and the dominance of a virtual community whose mean population size
(per species) equals the population size of the focal species. The new metrics
have at least two immediate applications: (i) acting as proxies for diversity
in diversity-stability modeling (ii) replacing population abundance in
reconstructing species dominance networks. The first application is
demonstrated here using data from a longitudinal study of the human vaginal
microbiome, and provides new insights relevant for microbial community
stability and disease etiology.


  Deciphering the links between amino acid sequence and amyloid fibril
formation is key for understanding protein misfolding diseases. Here we use
Monte Carlo simulations to study aggregation of short peptides in a
coarse-grained model with hydrophobic-polar (HP) amino acid sequences and
correlated side chain orientations for hydrophobic contacts. A significant
heterogeneity is observed in the aggregate structures and in the thermodynamics
of aggregation for systems of different HP sequences and different number of
peptides. Fibril-like ordered aggregates are found for several sequences that
contain the common HPH pattern while other sequences may form helix bundles or
disordered aggregates. A wide variation of the aggregation transition
temperatures among sequences, even among those of the same hydrophobic
fraction, indicates that not all sequences undergo aggregation at a presumable
physiological temperature. The transition is found to be the most cooperative
for sequences forming fibril-like structures. For a fibril-prone sequence, it
is shown that fibril formation follows the nucleation and growth mechanism.
Interestingly, a binary mixture of peptides of an aggregation-prone and a
non-aggregation-prone sequence shows association and conversion of the latter
to the fibrillar structure. Our study highlights the role of sequence in
selecting fibril-like aggregates and also the impact of structural template on
fibril formation by peptides of unrelated sequences.


  The dynamics of a mosquito population depends heavily on climatic variables
such as temperature and precipitation. Since climate change models predict that
global warming will impact on the frequency and intensity of rainfall, it is
important to understand how these variables affect the mosquito populations. We
present a model of the dynamics of a {\it Culex quinquefasciatus} mosquito
population that incorporates the effect of rainfall and use it to study the
influence of the number of rainy days and the mean monthly precipitation on the
maximum yearly abundance of mosquitoes $M_{max}$. Additionally, using a
fracturing process, we investigate the influence of the variability in daily
rainfall on $M_{max}$. We find that, given a constant value of monthly
precipitation, there is an optimum number of rainy days for which $M_{max}$ is
a maximum. On the other hand, we show that increasing daily rainfall
variability reduces the dependence of $M_{max}$ on the number of rainy days,
leading also to a higher abundance of mosquitoes for the case of low mean
monthly precipitation. Finally, we explore the effect of the rainfall in the
months preceding the wettest season, and we obtain that a regimen with high
precipitations throughout the year and a higher variability tends to advance
slightly the time at which the peak mosquito abundance occurs, but could
significantly change the total mosquito abundance in a year.


  Many discrete models of biological networks rely exclusively on Boolean
variables and many tools and theorems are available for analysis of strictly
Boolean models. However, multilevel variables are often required to account for
threshold effects, in which knowledge of the Boolean case does not generalise
straightforwardly. This motivated the development of conversion methods for
multilevel to Boolean models. In particular, Van Ham's method has been shown to
yield a one-to-one, neighbour and regulation preserving dynamics, making it the
de facto standard approach to the problem. However, Van Ham's method has
several drawbacks: most notably, it introduces vast regions of "non-admissible"
states that have no counterpart in the multilevel, original model. This raises
special difficulties for the analysis of interaction between variables and
circuit functionality, which is believed to be central to the understanding of
dynamic properties of logical models. Here, we propose a new multilevel to
Boolean conversion method, with software implementation. Contrary to Van Ham's,
our method doesn't yield a one-to-one transposition of multilevel trajectories,
however, it maps each and every Boolean state to a specific multilevel state,
thus getting rid of the non-admissible regions and, at the expense of
(apparently) more complicated, "parallel" trajectories. One of the prominent
features of our method is that it preserves dynamics and interaction of
variables in a certain manner. As a demonstration of the usability of our
method, we apply it to construct a new Boolean counter-example to the
well-known conjecture that a local negative circuit is necessary to generate
sustained oscillations. This result illustrates the general relevance of our
method for the study of multilevel logical models.


  In this work we study the limit distribution of an appropriately normalized
cophenetic index of the pure-birth tree conditioned on $n$ contemporary tips.
We show that this normalized phylogenetic balance index is a submartingale that
converges almost surely and in L2. We link our work with studies on trees
without branch lengths and show that in this case the limit distribution is a
contraction-type distribution, similar to the Quicksort limit distribution. In
the continuous branch case we suggest approximations to the limit distribution.
We propose heuristic methods of simulating from these distributions and it may
be observed that these algorithms result in reasonable tails. Therefore, we
propose a way based on the quantiles of the derived distributions for
hypothesis testing, whether an observed phylogenetic tree is consistent with
the pure-birth process. Simulating a sample by the proposed heuristics is
rapid, while exact simulation (simulating the tree and then calculating the
index) is a time-consuming procedure. We conduct a power study to investigate
how well the cophenetic indices detect deviations from the Yule tree and apply
the methodology to empirical phylogenies.


  We study a birth and death model for the adapatation of a sexual population
to an environment. The population is structured by a phenotypical trait, and,
possibly, an age variable. Recombination is modeled by Fisher's infinitesimal
operator. We prove the existence of principal eigenelements for the
corresponding eigenproblem. As the infinitesimal operator is 1-homogeneous but
nor linear nor monotone, the general Krein-Rutman theory cannot be applied to
this problem.


  The present paper deals with a prey-predator model with prey refuge
proportion to both species and independent harvesting of each species. Our
study shows that using refuge as control, it can break the limit circle of the
system and reach the required state of equilibrium level. It is established the
optimal harvesting policy. The boundedness, feasibility of interior equilibria,
bionomic equilibrium have been determined. The main observation is that
coefficient of refuge plays an important role in regulating the dynamics of the
present system. Moreover the variation of the coefficient of refuge change the
system from stable to unstable and vice-versa. Some numerical illustration are
given in order to support of our analytical and theoretical findings.


  Human populations have a complex history of introgression and of changing
population size. Human genetic variation has been affected by both these
processes, so that inference of past population size depends upon the pattern
of gene flow and introgression among past populations. One remarkable aspect of
human population history as inferred from genetics is a consistent "wave" of
larger effective population size, prior to the bottlenecks and expansions of
the last 100,000 years. Here I carry out a series of simulations to investigate
how introgression and gene flow from genetically divergent ancestral
populations affect the inference of ancestral effective population size. Both
introgression and gene flow from an extinct, genetically divergent population
consistently produce a wave in the history of inferred effective population
size. The time and amplitude of the wave reflect the time of origin of the
genetically divergent ancestral populations and the strength of introgression
or gene flow. These results demonstrate that even small fractions of
introgression or gene flow from ancient populations may have large effects on
the inference of effective population size.


  An SEIRS epidemic with disease fatalities is introduced in a growing
population (modelled as a super-critical linear birth and death process). The
study of the initial phase of the epidemic is stochastic, while the analysis of
the major outbreaks is deterministic. Depending on the values of the
parameters, the following scenarios are possible. i) The disease dies out
quickly, only infecting few; ii) the epidemic takes off, the \textit{number} of
infected individuals grows exponentially, but the \textit{fraction} of infected
individuals remains negligible; iii) the epidemic takes off, the
\textit{number} of infected grows initially quicker than the population, the
disease fatalities diminish the growth rate of the population, but it remains
super critical, and the \emph{fraction} of infected go to an endemic
equilibrium; iv) the epidemic takes off, the \textit{number} of infected
individuals grows initially quicker than the population, the diseases
fatalities turn the exponential growth of the population to an exponential
decay.


  Biological populations are subject to fluctuating environmental conditions.
Different adaptive strategies can allow them to cope with these fluctuations:
specialization to one particular environmental condition, adoption of a
generalist phenotype that compromise between conditions, or population-wise
diversification (bet-hedging). Which strategy provides the largest selective
advantage in the long run depends on the range of accessible phenotypes and the
statistics of the environmental fluctuations. Here, we analyze this problem in
a simple mathematical model of population growth. First, we review and extend a
graphical method to identify the nature of the optimal strategy when the
environmental fluctuations are uncorrelated. Temporal correlations in
environmental fluctuations open up new strategies that rely on memory but are
mathematically challenging to study: we present here new analytical results to
address this challenge. We illustrate our general approach by analyzing optimal
adaptive strategies in the presence of trade-offs that constrain the range of
accessible phenotypes. Our results extend several previous studies and have
applications to a variety of biological phenomena, from antibiotic resistance
in bacteria to immune responses in vertebrates.


  Among several quantitative invariants found in evolutionary genomics, one of
the most striking is the scaling of the overall abundance of proteins, or
protein domains, sharing a specific functional annotation across genomes of
given size. The size of these functional categories change, on average, as
power-laws in the total number of protein-coding genes. Here, we show that such
regularities are not restricted to the overall behavior of high-level
functional categories, but also exist systematically at the level of single
evolutionary families of protein domains. Specifically, the number of proteins
within each family follows family-specific scaling laws with genome size.
Functionally similar sets of families tend to follow similar scaling laws, but
this is not always the case. To understand this systematically, we provide a
comprehensive classification of families based on their scaling properties.
Additionally, we develop a quantitative score for the heterogeneity of the
scaling of families belonging to a given category or predefined group. Under
the common reasonable assumption that selection is driven solely or mainly by
biological function, these findings point to fine-tuned and interdependent
functional roles of specific protein domains, beyond our current functional
annotations. This analysis provides a deeper view on the links between
evolutionary expansion of protein families and the functional constraints
shaping the gene repertoire of bacterial genomes.


  Background: Over the past few decades, numerous forecasting methods have been
proposed in the field of epidemic forecasting. Such methods can be classified
into different categories such as deterministic vs. probabilistic, comparative
methods vs. generative methods, and so on. In some of the more popular
comparative methods, researchers compare observed epidemiological data from
early stages of an outbreak with the output of proposed models to forecast the
future trend and prevalence of the pandemic. A significant problem in this area
is the lack of standard well-defined evaluation measures to select the best
algorithm among different ones, as well as for selecting the best possible
configuration for a particular algorithm.
  Results: In this paper, we present an evaluation framework which allows for
combining different features, error measures, and ranking schema to evaluate
forecasts. We describe the various epidemic features (Epi-features) included to
characterize the output of forecasting methods and provide suitable error
measures that could be used to evaluate the accuracy of the methods with
respect to these Epi-features. We focus on long-term predictions rather than
short-term forecasting and demonstrate the utility of the framework by
evaluating six forecasting methods for predicting influenza in the United
States. Our results demonstrate that different error measures lead to different
rankings even for a single Epi-feature. Further, our experimental analyses show
that no single method dominates the rest in predicting all Epi-features, when
evaluated across error measures. As an alternative, we provide various
consensus ranking schema that summarizes individual rankings, thus accounting
for different error measures. We believe that a comprehensive evaluation
framework, as presented in this paper, will add value to the computational
epidemiology community.


  While many short read assemblers attempt to simplify the de Brujin graph by
identifying and resolving variant-induced bubbles to produce a haploid mosaic
result, this approach is only viable when variants are relatively rare and the
bubbles are well defined in a graph context. We observed that diploid genomes
with very high levels of heterozygosity fail to display well-resolved bubble
structures in a typical assembly graph and thus result in highly fragmented and
incomplete assemblies. Here we present an enhancement of Meraculous2 algorithm,
called Meraculous-2D, which preserves haplotypes across variant sites and
generates accurate assembly of highly heterozygous diploid genomes. Preserving
and taking advantage of the allelic variation throughout the assembly process
allows reconstructing both haplomes at once, without the need to pick arbitrary
paths through bubble structures. We also enhanced the original diploidy
resolution method of Meraculous2 to maintain and report phased haplotype
variant information.


  Electrostatic interactions play a fundamental role in the structure and
function of proteins. Due to ionizable amino acid residues present on the
solvent-exposed surfaces of proteins, the protein charge is not constant but
varies with the changes in the environment -- most notably, the pH of the
surrounding solution. We study the effects of pH on the charge of four globular
proteins by expanding their surface charge distributions in terms of
multipoles. The detailed representation of the charges on the proteins is in
this way replaced by the magnitudes and orientations of the multipole moments
of varying order. Focusing on the three lowest-order multipoles -- the total
charge, dipole, and quadrupole moment -- we show that the value of pH
influences not only their magnitudes, but more notably and importantly also the
spatial orientation of their principal axes. Our findings imply important
consequences for the study of protein-protein interactions and the assembly of
both proteinaceous shells and patchy colloids with dissociable charge groups.


  Causal ordering of key events in the cell cycle is essential for proper
functioning of an organism. Yet, it remains a mystery how a specific temporal
program of events is maintained despite ineluctable stochasticity in the
biochemical dynamics which dictate timing of cellular events. We propose that
if a change of cell fate is triggered by the {\em time-integral} of the
underlying stochastic biochemical signal, rather than the original signal, then
a dramatic improvement in temporal specificity results. Exact analytical
results for stochastic models of hourglass-timers and pendulum-clocks, two
important paradigms for biological timekeeping, elucidate how temporal
specificity is achieved through time-integration. En route, we introduce a
natural representation for time-integrals of stochastic processes, provide an
analytical prescription for evaluating corresponding first-passage-time
distributions, and uncover a mechanism by which a population of identical cells
can spontaneously bifurcate into subpopulations of early and late responders,
depending on hierarchy of timescales in the dynamics. Moreover, our approach
reveals how time-integration of stochastic signals may be realized
biochemically, through a simple chemical reaction scheme.


  Collapsin response mediator protein CRMP2 (gene: DPYSL2) is crucial for
neuronal development. The homotetrameric CRMP2 complex is regulated via two
mechanisms, first by phosphorylation at, and second by reduction and oxidation
of the Cys504 residues of two adjacent subunits. Here, we analyzed the effects
of this redox switch on the protein in vitro combined with force field
molecular dynamics (MD). Earlier X-ray data contain the structure of the rigid
body of the molecule but lack the flexible C-terminus with the important sites
for phosphorylation and redox regulation. An in silico model for this part was
established by replica exchange simulations and homology modelling, which is
consistent with results gained from CD spectroscopy with recombinant protein.
Thermofluor data indicated that the protein aggregates at bivalent ion
concentrations below 200 mM. In simulations the protein surface was covered at
these conditions by large amounts of ions, which most likely prevent
aggregation. A tryptophan residue (Trp295) in close proximity to the forming
disulfide allowed the measurement of the structural relaxation of the rigid
body upon reduction by fluorescent quenching. We were also able to determine
the second order rate constant of CRMP2 oxidation by H2O2. The simulated
solvent accessible surface of the hydroxyl group of Ser518 significantly
increased upon reduction of the disulfide bond. Our results give first detailed
insight in the profound structural changes of the tetrameric CRMP2 due to
oxidation and indicate a tightly connected regulation by phosphorylation and
redox modification.


  Long-read sequencing has enabled the de novo assembly of several mammalian
genomes, but with high cost in computing. Here, we demonstrated de novo
assembly of mammalian genome using long reads in an efficient and inexpensive
workstation.


  In living cells, biochemical reactions are catalyzed by specific enzymes and
connect to one another by sharing substrates and products, forming complex
networks. In our previous studies, we established a framework determining the
responses to enzyme perturbations only from network topology, and then proved a
theorem, called the law of localization, explaining response patterns in terms
of network topology. In this paper, we generalize these results to reaction
networks with conserved concentrations, which allows us to study any reaction
systems. We also propose novel network characteristics quantifying robustness.
We compare E. coli metabolic network with randomly rewired networks, and find
that the robustness of the E. coli network is significantly higher than that of
the random networks.


  The Poland-Scheraga model for DNA denaturation, besides playing a central
role in applications, has been widely studied in the physical and mathematical
literature over the past decades. More recently a natural generalization has
been introduced in the biophysics literature to overcome the limits of the
original model, namely to allow an excess of bases -- i.e. a different length
of the two single stranded DNA chains -- and to allow slippages in the chain
pairing. The increased complexity of the model is reflected in the appearance
of configurational transitions when the DNA is in double stranded form. In a
previous work of two of the authors the generalized Poland-Scheraga model has
been analyzed thanks to a representation in terms of a bivariate renewal
process. In this work we exploit this representation farther and fully
characterize the path properties of the system, making therefore explicit the
geometric structures -- and the configurational transitions -- that are
observed when the polymer is in the double stranded form. What we prove is
that, when the excess of bases is not absorbed in a homogeneous fashion along
the double stranded chain, then it either condensates in a single macroscopic
loop or it accumulates into an unbound single strand free end.


  Thermodynamic bulk measurements of binding reactions critically rely on the
validity of the law of mass action and the assumption of a dilute solution. Yet
important biological systems such as allosteric ligand-receptor binding,
macromolecular crowding, or misfolded molecules may not follow this fundamental
law and require a particular reaction model. Here we introduce a fluctuation
theorem for ligand binding and an experimental approach using single-molecule
force-spectroscopy to determine binding energies, selectivity and allostery of
nucleic acids, proteins and peptides in a model-independent fashion. This work
extends the use of fluctuation theorems beyond unimolecular folding reactions,
bridging the thermodynamics of small systems and the basic laws of chemical
equilibrium.


  Transcription of the genetic message encoded chemically in the sequence of
the DNA template is carried out by a molecular machine called RNA polymerase
(RNAP). Backward or forward slippage of the nascent RNA with respect to the DNA
template strand give rise to a transcript that is, respectively, longer or
shorter than the corresponding template. We model a RNAP as a "Tape-copying
Turing machine" (TCTM) where the DNA template is the input tape while the
nascent RNA strand is the output tape. Although the TCTM always steps forward
the process is assumed to be stochastic that has a probability of occurrence
per unit time. The time taken by a TCTM for each single successful forward
stepping on the input tape, during which the output tape suffers lengthening or
shortening by $n$ units because of backward or forward slippage, is a random
variable; we report some of the statistical characteristics of this time by
using the formalism for calculation of the distributions of {\it first-passage
time}. The results are likely to find applications in the analysis of
experimental data on "programmed" transcriptional error caused by
transcriptional slippage which is a mode of "recoding" of genetic information.


  Because of the limitations of classical silicon based computational
technology, several alternatives to traditional method in form of
unconventional computing have been proposed. In this paper we will focus on DNA
computing which is showing the possibility of excellence for its massive
parallelism, potential for information storage, speed and energy efficiency. In
this paper we will describe how syllogistic reasoning by DNA tweezers can be
presented by the semantics of process calculus and DNA strand graph. Syllogism
is an essential ingredient for commonsense reasoning of an individual. This
paper enlightens the procedure to deduce a precise conclusion from a set of
propositions by using formal language theory in form of process calculus and
the expressive power of DNA strand graph.


  The colonic mucus layer is a dynamic and complex structure formed by secreted
and transmembrane mucins, which are high-molecular-weight and heavily
glycosylated proteins. Colonic mucus consists of a loose outer layer and a
dense epithelium-attached layer. The outer layer is inhabited by various
representatives of the human gut microbiota (HGM). Glycans of the colonic mucus
can be used by the HGM as a source of carbon and energy when dietary fibers are
not sufficiently available. Here, we analyzed 397 individual HGM genomes to
identify pathways for the cleavage of host-synthetized mucin glycans to
monosaccharides as well as for the catabolism of the derived monosaccharides.
Our key results are as follows: (i) Genes for the cleavage of mucin glycans
were found in 86% of the analyzed genomes, whereas genes for the catabolism of
derived monosaccharides were found in 89% of the analyzed genomes. (ii)
Comparative genomic analysis identified four alternative forms of the
monosaccharide-catabolizing enzymes and four alternative forms of
monosaccharide transporters. (iii) Eighty-five percent of the analyzed genomes
may be involved in exchange pathways for the monosaccharides derived from
cleaved mucin glycans. (iv) The analyzed genomes demonstrated different
abilities to degrade known mucin glycans. Generally, the ability to degrade at
least one type of mucin glycan was predicted for 81% of the analyzed genomes.
(v) Eighty-two percent of the analyzed genomes can form mutualistic pairs that
are able to degrade mucin glycans and are not degradable by any of the paired
organisms alone. Taken together, these findings provide further insight into
the inter-microbial communications of the HGM as well as into host-HGM
interactions.


  Artificial neural networks (ANNs) have gained a well-deserved popularity
among machine learning tools upon their recent successful applications in
image- and sound processing and classification problems. ANNs have also been
applied for predicting the family or function of a protein, knowing its residue
sequence. Here we present two new ANNs with multi-label classification ability,
showing impressive accuracy when classifying protein sequences into 698 UniProt
families (AUC=99.99%) and 983 Gene Ontology classes (AUC=99.45%).


  Carnivores are important components of ecosystems with wide-ranging effects
on ecological communities.We studied the carnivore community in the Apostle
Islands National Lakeshore (APIS), where the presence, distribution, and
populations of carnivores was largely unknown. We developed a systematic method
to deploy camera traps across a grid while targeting fine-scale features to
maximize carnivore detection (Appendix 1), including systematic methods for
organizing and tagging the photo data (Appendix 2). We deployed 88 cameras on
13 islands from 2014-2016. We collected 92,694 photographs across 18,721 trap
nights, including 3,591 wildlife events and 1,070 carnivore events. We had a
mean of 6.6 cameras per island (range 2-30), and our camera density averaged
1.23 (range 0.74-3.08) cameras/ km2. We detected 27 species and 10 terrestrial
carnivores, including surprising detections of American martens (Martes
americana) and gray wolves (Canis lupus). The mean richness of carnivores on an
island was 3.23 (range 0-10). The best single variable to explain carnivore
richness on the Apostle Islands was island size, while the best model was
island size (positive correlation) and distance from mainland (negative
correlation) (R2 = 0.92). Relative abundances for carnivores ranged from a low
of 0.01 for weasels (Mustela spp.) to a high of 2.64 for black bears (Ursus
americanus), and the relative abundance of a species was significantly
correlated with the number of islands on which they were found. Carnivore
occupancy ranged from lows of 0.09 for gray wolves and 0.11 for weasels to a
high of 0.82 for black bears. Fuller understanding of APIS ecology will require
on-going monitoring of carnivores to evaluate temporal dynamics as well as
related ecological evaluations (e.g. small mammal dynamics, plant community
dynamics) to understand trophic effects.


  In this article we study the treewidth of the \emph{display graph}, an
auxiliary graph structure obtained from the fusion of phylogenetic (i.e.,
evolutionary) trees at their leaves. Earlier work has shown that the treewidth
of the display graph is bounded if the trees are in some formal sense
topologically similar. Here we further expand upon this relationship. We
analyse a number of reduction rules which are commonly used in the
phylogenetics literature to obtain fixed parameter tractable algorithms. In
some cases (the \emph{subtree} reduction) the reduction rules behave similarly
with respect to treewidth, while others (the \emph{cluster} reduction) behave
very differently, and the behaviour of the \emph{chain reduction} is
particularly intriguing because of its link with graph separators and forbidden
minors. We also show that the gap between treewidth and Tree Bisection and
Reconnect (TBR) distance can be infinitely large, and that unlike, for example,
planar graphs the treewidth of the display graph can be as much as linear in
its number of vertices. On a slightly different note we show that if a display
graph is formed from the fusion of a phylogenetic network and a tree, rather
than from two trees, the treewidth of the display graph is bounded whenever the
tree can be topologically embedded ("displayed") within the network. This opens
the door to the formulation of the display problem in Monadic Second Order
Logic (MSOL). A number of other auxiliary results are given. We conclude with a
discussion and list a number of open problems.