  We consider interacting population systems of predator-prey type, presenting
four models of control strategies for epidemics among the prey. In particular
to contain the transmissible disease, safety niches are considered, assuming
they lessen the disease spread, but do not protect prey from predators. This
represents a novelty with respect to standard ecosystems where the refuge
prevents predators' attacks. The niche is assumed either to protect the healthy
individuals, or to hinder the infected ones to get in contact with the
susceptibles, or finally to reduce altogether contacts that might lead to new
cases of the infection. In addition a standard culling procedure is also
analysed. The effectiveness of the different strategies are compared. Probably
the environments providing a place where disease carriers cannot come in
contact with the healthy individuals, or where their contact rates are lowered,
seem to preferable for disease containment.


  The investigations of the bat White Nose Syndrome (WNS) have yet to provide
answers as to how the causative fungus Pseudogymnoascus (Geomyces) destructans
(Pd) first appeared in the Northeast and how a single clone has spread rapidly
in the US and Canada. We aimed to catalogue Pd and all other fungi (mycobiome)
by the culture-dependent (CD) and culture-independent (CI) methods in four
Mines and two Caves from the epicenter of WNS zoonotic. Six hundred sixty-five
fungal isolates were obtained by CD method including the live recovery of Pd.
Seven hundred three nucleotide sequences that met the definition of operational
taxonomic units (OTUs) were recovered by CI methods. Most OTUs belonged to
unidentified clones deposited in the databases as environmental nucleic acid
sequences (ENAS). The core mycobiome of WNS affected sites comprised of 46
species of fungi from 31 genera recovered in culture, and 17 fungal genera and
31 ENAS identified from clone libraries. Fungi such as Arthroderma spp.,
Geomyces spp., Kernia spp., Mortierella spp., Penicillium spp., and
Verticillium spp. were predominant in culture while Ganoderma spp., Geomyces
spp., Mortierella spp., Penicillium spp. and Trichosporon spp. were abundant is
clone libraries. Alpha diversity analyses from CI data revealed that fungal
community structure was highly diverse. However, the true species diversity
remains undetermined due to under sampling. The frequent recovery of Pd
indicated that the pathogen has adapted to WNS-afflicted habitats. Further,
this study supports the hypothesis that Pd is an introduced species. These
findings underscore the need for integrated WNS control measures that target
both bats and the fungal pathogen.


  In many areas of computational biology, hidden Markov models (HMMs) have been
used to model local genomic features. In particular, coalescent HMMs have been
used to infer ancient population sizes, migration rates, divergence times, and
other parameters such as mutation and recombination rates. As more loci,
sequences, and hidden states are added to the model, however, the runtime of
coalescent HMMs can quickly become prohibitive. Here we present a new algorithm
for reducing the runtime of coalescent HMMs from quadratic in the number of
hidden time states to linear, without making any additional approximations. Our
algorithm can be incorporated into various coalescent HMMs, including the
popular method PSMC for inferring variable effective population sizes. Here we
implement this algorithm to speed up our demographic inference method diCal,
which is equivalent to PSMC when applied to a sample of two haplotypes. We
demonstrate that the linear-time method can reconstruct a population size
change history more accurately than the quadratic-time method, given similar
computation resources. We also apply the method to data from the 1000 Genomes
project, inferring a high-resolution history of size changes in the European
population.


  In this paper we consider a three level food web subject to a disease
affecting the bottom prey. The resulting dynamics is much richer with respect
to the purely demographic model, in that it contains more transcritical
bifurcations, gluing together the various equilibria, as well as persistent
limit cycles, which are shown to be absent in the classical case. Finally,
bistability is discovered among some equilibria, leading to situations in which
the computation of their basins of attraction is relevant for the system
outcome in terms of its biological implications.


  Network epidemiology has become a core framework for investigating the role
of human contact patterns in the spreading of infectious diseases. In network
epidemiology represents the contact structure as a network of nodes
(individuals) connected by links (sometimes as a temporal network where the
links are not continuously active) and the disease as a compartmental model
(where individuals are assigned states with respect to the disease and follow
certain transition rules between the states). In this paper, we discuss fast
algorithms for such simulations and also compare two commonly used versions -
one where there is a constant recovery rate (the number of individuals that
stop being infectious per time is proportional to the number of such people),
the other where the duration of the disease is constant. We find that, for most
practical purposes, these versions are qualitatively the same.


  Influenza B viruses make a considerable contribution to morbidity attributed
to seasonal influenza. Currently circulating influenza B isolates are known to
belong to two antigenically distinct lineages referred to as B/Victoria and
B/Yamagata. Frequent exchange of genomic segments of these two lineages has
been noted in the past, but the observed patterns of reassortment have not been
formalized in detail. We investigate inter-lineage reassortments by comparing
phylogenetic trees across genomic segments. Our analyses indicate that of the 8
segments of influenza B viruses only PB1, PB2 and HA segments maintained
separate Victoria and Yamagata lineages and that currently circulating strains
possess PB1, PB2 and HA segments derived entirely from one or the other
lineage; other segments have repeatedly reassorted between lineages thereby
reducing genetic diversity. We argue that this difference between segments is
due to selection against reassortant viruses with mixed lineage PB1, PB2 and HA
segments. Given sufficient time and continued recruitment to the
reassortment-isolated PB1-PB2-HA gene complex, we expect influenza B viruses to
eventually undergo sympatric speciation.


  In this paper we present and analyse a simple two populations model with
migrations among two different environments. The populations interact by
competing for resources. Equilibria are investigated. A proof for the
boundedness of the populations is provided. A kind of competitive exclusion
principle for metapopulation systems is obtained. At the same time we show that
the competitive exclusion principle at the local patch level may be prevented
to hold by the migration phenomenon, i.e. two competing populations may
coexist, provided that only one of them is allowed to freely move or that
migrations for both occur just in one direction.


  Birds in a flock move in a correlated way, resulting in large polarization of
velocities. A good understanding of this collective behavior exists for linear
motion of the flock. Yet observing actual birds, the center of mass of the
group often turns giving rise to more complicated dynamics, still keeping
strong polarization of the flock. Here we propose novel dynamical equations for
the collective motion of polarized animal groups that account for correlated
turning including solely social forces. We exploit rotational symmetries and
conservation laws of the problem to formulate a theory in terms of generalized
coordinates of motion for the velocity directions akin to a Hamiltonian
formulation for rotations. We explicitly derive the correspondence between this
formulation and the dynamics of the individual velocities, thus obtaining a new
model of collective motion. In the appropriate overdamped limit we recover the
well-known Vicsek model, which dissipates rotational information and does not
allow for polarized turns. Although the new model has its most vivid success in
describing turning groups, its dynamics is intrinsically different from
previous ones in a wide dynamical regime, while reducing to the hydrodynamic
description of Toner and Tu at very large length-scales. The derived framework
is therefore general and it may describe the collective motion of any strongly
polarized active matter system.


  Network measures that reflect the most salient properties of complex
large-scale networks are in high demand in the network research community. In
this paper we adapt a combinatorial measure of negative curvature (also called
hyperbolicity) to parameterized finite networks, and show that a variety of
biological and social networks are hyperbolic. This hyperbolicity property has
strong implications on the higher-order connectivity and other topological
properties of these networks. Specifically, we derive and prove bounds on the
distance among shortest or approximately shortest paths in hyperbolic networks.
We describe two implications of these bounds to cross-talk in biological
networks, and to the existence of central, influential neighborhoods in both
biological and social networks.


  The classic model of eukaryotic gene expression requires direct spatial
contact between a distal enhancer and a proximal promoter. Recent Chromosome
Conformation Capture (3C) studies show that enhancers and promoters are
embedded in a complex network of looping interactions. Here we use a polymer
model of chromatin fiber to investigate whether, and to what extent, looping
interactions between elements in the vicinity of an enhancer-promoter pair can
influence their contact frequency. Our equilibrium polymer simulations show
that a chromatin loop, formed by elements flanking either an enhancer or a
promoter, suppresses enhancer-promoter interactions, working as an insulator. A
loop formed by elements located in the region between an enhancer and a
promoter, on the contrary, facilitates their interactions. We find that
different mechanisms underlie insulation and facilitation; insulation occurs
due to steric exclusion by the loop, and is a global effect, while facilitation
occurs due to an effective shortening of the enhancer-promoter genomic
distance, and is a local effect. Consistently, we find that these effects
manifest quite differently for in silico 3C and microscopy. Our results show
that looping interactions that do not directly involve an enhancer-promoter
pair can nevertheless significantly modulate their interactions. This
phenomenon is analogous to allosteric regulation in proteins, where a
conformational change triggered by binding of a regulatory molecule to one site
affects the state of another site.


  Plant domestication involved a process of selection through human agency of a
series of traits collectively termed the domestication syndrome. Current debate
concerns the pace at which domesticated plants emerged from cultivated wild
populations and how many genes were involved. Here we present simulations that
test how many genes could have been involved by considering the cost of
selection. We demonstrate the selection load that can be endured by populations
increases with decreasing selection coefficients and greater numbers of loci
down to values of about s = 0.005, causing a driving force that increases the
number of loci under selection. As the number of loci under selection
increases, an effect of co-selection increases resulting in individual unlinked
loci being fixed more rapidly in out-crossing populations, representing a
second driving force to increase the number of loci under selection. In
inbreeding systems co-selection results in interference and reduced rates of
fixation but does not reduce the size of the selection load that can be
endured. These driving forces result in an optimum pace of genome evolution in
which 50-100 loci are the most that could be under selection in a cultivation
regime. Furthermore, the simulations do not preclude the existence of selective
sweeps but demonstrate that they come at a cost of the selection load that can
be endured and consequently a reduction of the capacity of plants to adapt to
new environments, which may contribute to the explanation of why selective
sweeps have been so rarely detected in genome studies.


  A long genomic segment inherited by a pair of individuals from a single,
recent common ancestor is said to be identical-by-descent (IBD). Shared IBD
segments have numerous applications in genetics, from demographic inference to
phasing, imputation, pedigree reconstruction, and disease mapping. Here, we
provide a theoretical analysis of IBD sharing under Markovian approximations of
the coalescent with recombination. We describe a general framework for the IBD
process along the chromosome under the Markovian models (SMC/SMC'), as well as
introduce and justify a new model, which we term the renewal approximation,
under which lengths of successive segments are independent. Then, considering
the infinite-chromosome limit of the IBD process, we recover previous results
(for SMC) and derive new results (for SMC') for the mean number of shared
segments longer than a cutoff and the fraction of the chromosome found in such
segments. We then use renewal theory to derive an expression (in Laplace space)
for the distribution of the number of shared segments and demonstrate
implications for demographic inference. We also compute (again, in Laplace
space) the distribution of the fraction of the chromosome in shared segments,
from which we obtain explicit expressions for the first two moments. Finally,
we generalize all results to populations with a variable effective size.


  Dynamic models of biochemical networks typically consist of sets of
non-linear ordinary differential equations involving states (concentrations or
amounts of the components of the network) and parameters describing the
reaction kinetics. Unfortunately, in most cases the parameters are completely
unknown or only rough estimates of their values are available. Therefore, their
values must be estimated from time-series experimental data.
  In recent years, it has been suggested that dynamic systems biology models
are universally sloppy so their parameters cannot be uniquely estimated. In
this work, we re-examine this concept, establishing links with the notions of
identifiability and experimental design. Further, considering a set of
examples, we address the following fundamental questions: i) is sloppiness
inherent to model structure?; ii) is sloppiness influenced by experimental data
or noise?; iii) does sloppiness mean that parameters cannot be identified?, and
iv) can sloppiness be modified by experimental design?
  Our results indicate that sloppiness is not equivalent to lack of structural
or practical identifiability (although they can be related), so sloppy models
can be identifiable. Therefore, drawing conclusions about the possibility of
estimating unique parameter values by sloppiness analysis can be misleading.
Checking structural and practical identifiability analyses is a better approach
to asses the uniqueness and confidence in parameter estimation.


  Let $\mbox{interval} + k v$, $\mbox{interval} + k e$, and $\mbox{interval} -
k e$ denote the classes of graphs that can be obtained from some interval graph
by adding $k$ vertices, adding $k$ edges, and deleting $k$ edges, respectively.
When $k$ is small, these graph classes are called almost interval graphs. They
are well motivated from computational biology, where the data ought to be
represented by an interval graph while we can only expect an almost interval
graph for the best. For any fixed $k$, we give linear-time algorithms for
recognizing all these classes, and in the case of membership, our algorithms
provide also a specific interval graph as evidence. When $k$ is part of the
input, these problems are also known as graph modification problems, all
NP-complete. Our results imply that they are fixed-parameter tractable
parameterized by $k$, thereby resolving the long-standing open problem on the
parameterized complexity of recognizing $\mbox{interval}+ k e$, first asked by
Bodlaender et al. [Bioinformatics, 11:49--57, 1995]. Moreover, our algorithms
for recognizing $\mbox{interval}+ k v$ and $\mbox{interval}- k e$ run in times
$O(6^k \cdot (n + m))$ and $O(8^k \cdot (n + m))$, (where $n$ and $m$ stand for
the numbers of vertices and edges respectively in the input graph,)
significantly improving the $O(k^{2k}\cdot n^3m)$-time algorithm of Heggernes
et al. [STOC 2007] and the $O(10^k \cdot n^9)$-time algorithm of Cao and Marx
[SODA 2014] respectively.


  Secondary contact between divergent populations or incipient species may
result in the exchange and introgression of genomic material. We develop a
simple DNA sequence measure, called Gmin, which is designed to identify genomic
regions experiencing introgression in a secondary contact model. Gmin is
defined as the ratio of the minimum between-population number of nucleotide
differences to the average number of between-population differences. One
advantage of Gmin is that it is computationally inexpensive relative to
model-based methods for detecting gene flow and it scales easily to the level
of whole-genome analysis. We compare the sensitivity and specificity of Gmin to
those of the widely used index of population differentiation, Fst, and suggest
a simple statistical test for identifying genomic outliers. Extensive computer
simulations demonstrate that Gmin has both greater sensitivity and specificity
for detecting recent introgression than does Fst. Furthermore, we find that the
sensitivity of Gmin is robust with respect to both the population mutation and
recombination rates. Finally, a scan of Gmin across the X chromosome of
Drosophila melanogaster identifies candidate regions of introgression between
sub-Saharan African and cosmopolitan populations that were previously missed by
other methods. These results show that Gmin is a biologically straightforward,
yet powerful, alternative to Fst, as well as to more computationally intensive
model-based methods for detecting gene flow.


  The majority of current phylogenetic comparative methods assume that the
stochastic evolutionary process is homogeneous over the phylogeny or offer
relaxations of this in rather limited and usually parameter expensive ways.
Here we make a preliminary investigation, by means of a numerical experiment,
whether the Laplace motion process can offer an alternative approach.


  Natural selection favors the more successful individuals. This is the
elementary premise that pervades common models of evolution. Under extreme
conditions, however, the process may no longer be probabilistic. Those that
meet certain conditions survive and may reproduce while others perish. By
introducing the corresponding binary birth-death dynamics to spatial
evolutionary games, we observe solutions that are fundamentally different from
those reported previously based on imitation dynamics. Social dilemmas
transform to collective enterprises, where the availability of free expansion
ranges and limited exploitation possibilities dictates self-organized growth.
Strategies that dominate are those that are collectively most apt in meeting
the survival threshold, rather than those who succeed in exploiting others for
unfair benefits. Revisiting Darwinian principles with the focus on survival
rather than imitation thus reveals the most counterintuitive ways of
reconciling cooperation with competition.


  The relationship between components of biochemical network and the resulting
dynamics of the overall system is a key focus of computational biology.
However, as these networks and resulting mathematical models are inherently
complex and non-linear, the understanding of this relationship becomes
challenging. Among many approaches, model reduction methods provide an avenue
to extract components responsible for the key dynamical features of the system.
Unfortunately, these approaches often require intuition to apply. In this
manuscript we propose a practical algorithm for the reduction of biochemical
reaction systems using fast-slow asymptotics. This method allows the ranking of
system variables according to how quickly they approach their momentary steady
state, thus selecting the fastest for a steady state approximation. We applied
this method to derive models of the Nuclear Factor kappa B network, a key
regulator of the immune response that exhibits oscillatory dynamics. Analyses
with respect to two specific solutions, which corresponded to different
experimental conditions identified different components of the system that were
responsible for the respective dynamics. This is an important demonstration of
how reduction methods that provide approximations around a specific steady
state, could be utilised in order to gain a better understanding of network
topology in a broader context.


  In mammals, most cells in the brain and peripheral tissues generate circadian
(~24hr) rhythms autonomously. These self-sustained rhythms are coordinated and
entrained by a master circadian clock in the suprachiasmatic nucleus (SCN).
Within the SCN, the individual rhythms of each neuron are synchronized through
intercellular signaling. One important feature of SCN is that the synchronized
period is close to the cell population mean of intrinsic periods. In this way,
the synchronized period of the SCN stays close to the periods of cells in
peripheral tissues. This is important for SCN to entrain cells throughout the
body. However, the mechanism that drives the period of the coupled SCN cells to
the population mean is not known. We use mathematical modeling and analysis to
show that the mechanism of transcription repression plays a pivotal role in
regulating the coupled period. Specifically, we use phase response curve
analysis to show that the coupled period within the SCN stays near the
population mean if transcriptional repression occurs via protein sequestration.
In contrast, the coupled period is far from the mean if repression occurs
through highly nonlinear Hill-type regulation (e.g. oligomer- or
phosphorylation-based repression). Furthermore, we find that the timescale of
intercellular coupling needs to be fast compared to that of intracellular
feedback to maintain the mean period. These findings reveal the important
relationship between the intracellular transcriptional feedback loop and
intercellular coupling. This relationship explains why transcriptional
repression appears to occur via protein sequestration in multicellular
organisms, mammals and Drosophila, in contrast with the phosphorylation-based
repression in unicellular organisms. That is, transition to protein
sequestration is essential for synchronizing multiple cells with a period close
to the population mean (~24hr).


  Molecular phenotypes link genomic information with organismic functions,
fitness, and evolution. Quantitative traits are complex phenotypes that depend
on multiple genomic loci. In this paper, we study the adaptive evolution of a
quantitative trait under time-dependent selection, which arises from
environmental changes or through fitness interactions with other co-evolving
phenotypes. We analyze a model of trait evolution under mutations and genetic
drift in a single-peak fitness seascape. The fitness peak performs a
constrained random walk in the trait amplitude, which determines the
time-dependent trait optimum in a given population. We derive analytical
expressions for the distribution of the time-dependent trait divergence between
populations and of the trait diversity within populations. Based on this
solution, we develop a method to infer adaptive evolution of quantitative
traits. Specifically, we show that the ratio of the average trait divergence
and the diversity is a universal function of evolutionary time, which predicts
the stabilizing strength and the driving rate of the fitness seascape. From an
information-theoretic point of view, this function measures the
macro-evolutionary entropy in a population ensemble, which determines the
predictability of the evolutionary process. Our solution also quantifies two
key characteristics of adapting populations: the cumulative fitness flux, which
measures the total amount of adaptation, and the adaptive load, which is the
fitness cost due to a population's lag behind the fitness peak.


  We consider the spread of infectious disease through contact networks of
Configuration Model type. We assume that the disease spreads through contacts
and infected individuals recover into an immune state. We discuss a number of
existing mathematical models used to investigate this system, and show
relations between the underlying assumptions of the models. In the process we
offer simplifications of some of the existing models. The distinctions between
the underlying assumptions are subtle, and in many if not most cases this
subtlety is irrelevant. Indeed, under appropriate conditions the models are
equivalent. We compare the benefits and disadvantages of the different models,
and discuss their application to other populations (\emph{e.g.,} clustered
networks). Finally we discuss ongoing challenges for network-based epidemic
modeling.


  We consider a class of evolution equations describing population dynamics in
the presence of a carrying capacity depending on the population with delay. In
an earlier work, we presented an exhaustive classification of the logistic
equation where the carrying capacity is linearly dependent on the population
with a time delay, which we refer to as the "linear delayed carrying capacity"
model. Here, we generalize it to the case of a nonlinear delayed carrying
capacity. The nonlinear functional form of the carrying capacity characterizes
the delayed feedback of the evolving population on the capacity of their
surrounding by either creating additional means for survival or destroying the
available resources. The previously studied linear approximation for the
capacity assumed weak feedback, while the nonlinear form is applicable to
arbitrarily strong feedback. The nonlinearity essentially changes the behavior
of solutions to the evolution equation, as compared to the linear case. All
admissible dynamical regimes are analyzed, which can be of the following types:
punctuated unbounded growth, punctuated increase or punctuated degradation to a
stationary state, convergence to a stationary state with sharp reversals of
plateaus, oscillatory attenuation, everlasting fluctuations, everlasting
up-down plateau reversals, and divergence in finite time. The theorem is proved
that, for the case characterizing the evolution under gain and competition,
solutions are always bounded, if the feedback is destructive. We find that even
a small noise level profoundly affects the position of the finite-time
singularities. Finally, we demonstrate the feasibility of predicting the
critical time of solutions having finite-time singularities from the knowledge
of a simple quadratic approximation of the early time dynamics.


  Programs based on hash tables and Burrows-Wheeler are very fast for mapping
short reads to genomes but have low accuracy in the presence of mismatches and
gaps. Such reads can be aligned accurately with the Smith-Waterman algorithm
but it can take hours and days to map millions of reads even for bacteria
genomes. We introduce a GPU program called MaxSSmap with the aim of achieving
comparable accuracy to Smith-Waterman but with faster runtimes. Similar to most
programs MaxSSmap identifies a local region of the genome followed by exact
alignment. Instead of using hash tables or Burrows-Wheeler in the first part,
MaxSSmap calculates maximum scoring subsequence score between the read and
disjoint fragments of the genome in parallel on a GPU and selects the highest
scoring fragment for exact alignment. We evaluate MaxSSmap's accuracy and
runtime when mapping simulated Illumina E.coli and human chromosome one reads
of different lengths and 10\% to 30\% mismatches with gaps to the E.coli genome
and human chromosome one. We also demonstrate applications on real data by
mapping ancient horse DNA reads to modern genomes and unmapped paired reads
from NA12878 in 1000 genomes. We show that MaxSSmap attains comparable high
accuracy and low error to fast Smith-Waterman programs yet has much lower
runtimes. We show that MaxSSmap can map reads rejected by BWA and NextGenMap
with high accuracy and low error much faster than if Smith-Waterman were used.
On short read lengths of 36 and 51 both MaxSSmap and Smith-Waterman have lower
accuracy compared to at higher lengths. On real data MaxSSmap produces many
alignments with high score and mapping quality that are not given by NextGenMap
and BWA. The MaxSSmap source code is freely available from
http://www.cs.njit.edu/usman/MaxSSmap.


  In this paper, the exact analytical solution of the
Susceptible-Infected-Recovered (SIR) epidemic model is obtained in a parametric
form. By using the exact solution we investigate some explicit models
corresponding to fixed values of the parameters, and show that the numerical
solution reproduces exactly the analytical solution. We also show that the
generalization of the SIR model, including births and deaths, described by a
nonlinear system of differential equations, can be reduced to an Abel type
equation. The reduction of the complex SIR model with vital dynamics to an Abel
type equation can greatly simplify the analysis of its properties. The general
solution of the Abel equation is obtained by using a perturbative approach, in
a power series form, and it is shown that the general solution of the SIR model
with vital dynamics can be represented in an exact parametric form.


  For many viruses assembly and budding occur simultaneously during virion
formation. Understanding the mechanisms underlying this process could promote
biomedical efforts to block viral propagation and enable use of capsids in
nanomaterials applications. To this end, we have performed molecular dynamics
simulations on a coarse-grained model that describes virus assembly on a
fluctuating lipid membrane. Our simulations show that the membrane can promote
association of adsorbed subunits through dimensional reduction, but also can
introduce barriers that inhibit complete assembly. We find several mechanisms,
including one not anticipated by equilibrium theories, by which membrane
microdomains, such as lipid rafts, can enhance assembly by reducing these
barriers. We show how this predicted mechanism can be experimentally tested.
Furthermore, the simulations demonstrate that assembly and budding depend
crucially on the system dynamics via multiple timescales related to membrane
deformation, protein diffusion, association, and adsorption onto the membrane.


  Social ecological systems are often difficult to investigate and manage
because of their inherent complexity1. Small variations in external drivers can
lead to abrupt changes associated with instabilities and bifurcations in the
underlying dynamics2-4. Anticipating critical transitions and divergence from
the present state of the system is particularly crucial to the prevention or
mitigation of the effects of unwanted and irreversible changes5-10. Recent
research in ecology has focused on leading indicators of regime shift in
ecosystems characterized by one state variable5,7,11,12. The case of systems
with several mutually interacting components, however, has remained poorly
investigated13, while the connection between network stability and research on
indicators for loss of resilience has been elusive14. Here we develop a
theoretical framework to analyze early warning signs of instability and regime
shift in social ecological networks. We provide analytical expressions for a
set of precursors of instability in social ecological systems with additive
noise for a variety of network structures. In particular, we show that the
covariance matrix of the dynamics can effectively anticipate the emergence of
instability. We also compare signals of early warning based on the dynamics of
suitably selected nodes, to indicators based on the integrated behavior of the
whole network. We find that the performances of these indicators are affected
by the network structure and the type of interaction among nodes. These results
provide new advances in multidimensional early warning analysis and offer a
framework to evaluate the resilience of social ecological networks.


  Within food webs, species can be partitioned into groups according to various
criteria. Two notions have received particular attention: trophic groups, which
have been used for decades in the ecological literature, and more recently,
modules. The relationship between these two group definitions remains unknown
in empirical food webs because they have so far been studied separately. While
recent developments in network theory have led to efficient methods for
detecting modules in food webs, the determination of trophic groups (sets of
species that are functionally similar) is based on subjective expert knowledge.
Here, we develop a novel algorithm for trophic group detection. We apply this
method to several well-resolved empirical food webs, and show that aggregation
into trophic groups allows the simplification of food webs while preserving
their information content. Furthermore, we reveal a 2-level hierarchical
structure where modules partition food webs into large bottom-top trophic
pathways whereas trophic groups further partition these pathways into sets of
species with similar trophic connections. Bringing together trophic groups and
modules provides new perspectives to the study of dynamical and functional
consequences of food-web structure, bridging topological analysis and dynamical
systems. Trophic groups have a clear ecological meaning in terms of trophic
similarity, and are found to provide a trade-off between network complexity and
information loss.


  We discuss an autocatalytic reaction system: the cyclic competition model A1
+ A2 --> 2 A2, A2 + A3 --> 2 A3, A3 + A4 --> 2 A4, A4 + A1 --> 2 A1), as well
as its neutral counterpart. Migrations are introduced into the model. When
stochastic phenomena are taken into account, a phase transition between a
``fixation'' and a ``neutral'' regime is observed. In the ``fixation'' regime,
species A1 and A3 form an alliance against species A2 and A4, and the final
state is one in which one of the symbiotic pairs has won. The odd--even
``coarse--grained'' systems is mapped onto the two--species neutral (Kimura)
model. In the ``neutral'' regime, all four species survive for long
(evolutionary) times. The analytical results are checked against computer
simulations of the model. The model is generalized for n species. Also, a
generalized version of the Volterra model is analysed.


  The metazoan genome is replicated in precise cell lineage specific temporal
order. However, the mechanism controlling this orchestrated process is poorly
understood as no molecular mechanisms have been identified that actively
regulate the firing sequence of genome replication. Here we develop a
mechanistic model of genome replication capable of predicting, with accuracy
rivaling experimental repeats, observed empirical replication timing program in
humans. In our model, replication is initiated in an uncoordinated
(time-stochastic) manner at well-defined sites. The model contains, in addition
to the choice of the genomic landmark that localizes initiation, only a single
adjustable parameter of direct biological relevance: the number of replication
forks. We find that DNase hypersensitive sites are optimal and independent
determinants of DNA replication initiation. We demonstrate that the DNA
replication timing program in human cells is a robust emergent phenomenon that,
by its very nature, does not require a regulatory mechanism determining a
proper replication initiation firing sequence.


  The Web of Life (www.web-of-life.es) provides a graphical user interface,
based on Google Maps, for easily visualizing and downloading data on ecological
networks of species interactions. It is designed and implemented in a
relational database management system, allowing sophisticated and user-friendly
searching across networks. Users can access the database by any web browser
using a variety of operating systems. Data can be downloaded in several common
formats, and a data transmission webservice in JavaScript Object Notation is
also provided.


  The compartmental models used to study epidemic spreading often assume the
same susceptibility for all individuals, and are therefore, agnostic about the
effects that differences in susceptibility can have on epidemic spreading. Here
we show that--for the SIS model--differential susceptibility can make networks
more vulnerable to the spread of diseases when the correlation between a node's
degree and susceptibility are positive, and less vulnerable when this
correlation is negative. Moreover, we show that networks become more likely to
contain a pocket of infection when individuals are more likely to connect with
others that have similar susceptibility (the network is segregated). These
results show that the failure to include differential susceptibility to
epidemic models can lead to a systematic over/under estimation of fundamental
epidemic parameters when the structure of the networks is not independent from
the susceptibility of the nodes or when there are correlations between the
susceptibility of connected individuals.


  Biological diversity is known to play an important role in generating and
maintaining ecosystem productivity and other functions, and has consequently
become a central focus of many efforts to preserve ecosystem services.
Theoretical parallels suggest the diversity of fishing fleets may have a
similarly important role in determining the productivity and ecological impacts
of fisheries, but this possibility has rarely been explored. Here I present
theoretical analyses showing that the diversity of m\'etiers - combinations of
technology, target species, and fishing grounds - and technical efficiencies in
a fishing fleet have important impacts on the productivity, profitability, and
ecological impacts of fisheries, particularly mixed-stock or multispecies
fisheries. Diversification of m\'etiers can increase yields and reduce threats
to weak stocks in both managed and unmanaged multispecies fisheries.
Diversification of technical efficiencies creates opportunities for larger
profits in managed fisheries, but often decreases yields and worsens impacts on
weak stocks in unmanaged fisheries. These results suggest that the potential
impact of management may be highest in fisheries with diverse fleets.


  A crop can be represented as a biotechnical system in which components are
either chosen (cultivar, management) or given (soil, climate) and whose
combination generates highly variable stress patterns and yield responses.
Here, we used modeling and simulation to predict the crop phenotypic plasticity
resulting from the interaction of plant traits (G), climatic variability (E)
and management actions (M). We designed two in silico experiments that compared
existing and virtual sunflower cultivars (Helianthus annuus L.) in a target
population of cropping environments by simulating a range of indicators of crop
performance. Optimization methods were then used to search for GEM combinations
that matched desired crop specifications. Computational experiments showed that
the fit of particular cultivars in specific environments is gradually
increasing with the knowledge of pedo-climatic conditions. At the regional
scale, tuning the choice of cultivar impacted crop performance the same
magnitude as the effect of yearly genetic progress made by breeding. When
considering virtual genetic material, designed by recombining plant traits,
cultivar choice had a greater positive impact on crop performance and
stability. Results suggested that breeding for key traits conferring plant
plasticity improved cultivar global adaptation capacity whereas increasing
genetic diversity allowed to choose cultivars with distinctive traits that were
more adapted to specific conditions. Consequently, breeding genetic material
that is both plastic and diverse may improve yield stability of agricultural
systems exposed to climatic variability. We argue that process-based modeling
could help enhancing spatial management of cultivated genetic diversity and
could be integrated in functional breeding approaches.


  A feature often observed in epidemiological networks is significant
heterogeneity in degree. A popular modelling approach to this has been to
consider large populations with highly heterogeneous discrete contact rates.
This paper defines an individual-level non-Markovian stochastic process that
converges on standard ODE models of such populations in the appropriate
asymptotic limit. A generalised Sellke construction is derived for this model,
and this is then used to consider final outcomes in the case where
heterogeneity follows a truncated Zipf distribution.


  A topological RNA structure is derived from a diagram and its shape is
obtained by collapsing the stacks of the structure into single arcs and by
removing any arcs of length one. Shapes contain key topological, information
and for fixed topological genus there exist only finitely many such shapes. We
shall express topological RNA structures as unicellular maps, i.e. graphs
together with a cyclic ordering of their half-edges. In this paper we prove a
bijection of shapes of topological RNA structures. We furthermore derive a
linear time algorithm generating shapes of fixed topological genus. We derive
explicit expressions for the coefficients of the generating polynomial of these
shapes and the generating function of RNA structures of genus $g$. Furthermore
we outline how shapes can be used in order to extract essential information of
RNA structure databases.


  We investigate the behaviour of the genealogy of a Wright-Fisher population
model under the influence of a strong seed-bank effect. More precisely, we
consider a simple seed-bank age distribution with two atoms, leading to either
classical or long genealogical jumps (the latter modeling the effect of
seed-dormancy). We assume that the length of these long jumps scales like a
power $N^\beta$ of the original population size $N$, thus giving rise to a
`strong' seed-bank effect. For a certain range of $\beta$, we prove that the
ancestral process of a sample of $n$ individuals converges under a
non-classical time-scaling to Kingman's $n-$coalescent. Further, for a wider
range of parameters, we analyze the time to the most recent common ancestor of
two individuals analytically and by simulation.


  The antibody repertoire of each individual is continuously updated by the
evolutionary process of B cell receptor mutation and selection. It has recently
become possible to gain detailed information concerning this process through
high-throughput sequencing. Here, we develop modern statistical molecular
evolution methods for the analysis of B cell sequence data, and then apply them
to a very deep short-read data set of B cell receptors. We find that the
substitution process is conserved across individuals but varies significantly
across gene segments. We investigate selection on B cell receptors using a
novel method that side-steps the difficulties encountered by previous work in
differentiating between selection and motif-driven mutation; this is done
through stochastic mapping and empirical Bayes estimators that compare the
evolution of in-frame and out-of-frame rearrangements. We use this new method
to derive a per-residue map of selection, which provides a more nuanced view of
the constraints on framework and variable regions.


  The emergence of multicellularity and developmental programs are among the
major problems of evolutionary biology. Traditionally, research in this area
has been based on the combination of data analysis and experimental work on one
hand and theoretical approximations on the other. A third possibility is
provided by computer simulation models, which allow to both simulate reality
and explore alternative possibilities. These in silico models offer a powerful
window to the possible and the actual by means of modeling how virtual cells
and groups of cells can evolve complex interactions beyond a set of isolated
entities. Here we present several examples of such models, each one
illustrating the potential for artificial modeling of the transition to
multicellularity.


  Expression of cellular genes is regulated by binding of transcription factors
to their promoter, either activating or inhibiting transcription of a gene.
Particularly interesting is the case when the expressed protein regulates its
own transcription. In this paper the features of this self-regulating process
are investigated. In the here presented model the gene can be in two states.
Either a protein is bound to its promoter or not. The steady state
distributions of protein during and at the end of both states are analyzed.
Moreover a powerful numerical method based on the corresponding master equation
to compute the protein distribution in the steady state is presented and
compared to an already existing method. Additionally the special case of
self-regulation, in which protein can only be produced, if one of these
proteins is bound to the promoter region, is analyzed. Furthermore a
self-regulating gene is compared to a similar gene, which also has two states
and produces the same amount of proteins but is not regulated by its
protein-product.


  The twenty protein coding amino acids are found in proteomes with different
relative abundances. The most abundant amino acid, leucine, is nearly an order
of magnitude more prevalent than the least abundant amino acid, cysteine. Amino
acid metabolic costs differ similarly, constraining their incorporation into
proteins. On the other hand, sequence diversity is necessary for protein
folding, function and evolution. Here we present a simple model for a
cost-diversity trade-off postulating that natural proteomes minimize amino acid
metabolic flux while maximizing sequence entropy. The model explains the
relative abundances of amino acids across a diverse set of proteomes. We found
that the data is remarkably well explained when the cost function accounts for
amino acid chemical decay. More than one hundred proteomes reach comparable
solutions to the trade-off by different combinations of cost and diversity.
Quantifying the interplay between proteome size and entropy shows that
proteomes can get optimally large and diverse.


  Cell systems consist of a huge number of various molecules that display
specific patterns of interactions, which have a determining influence on the
cell's functioning. In general, such complexity is seen to increase with the
complexity of the organism, with a concomitant increase of the accuracy and
specificity of the cellular processes. The question thus arises how the
complexification of systems - modeled here by simple interacting birth-death
type processes - can lead to a reduction of the noise - described by the
variance of the number of molecules. To gain understanding of this issue, we
investigated the difference between a single system containing molecules that
are produced and degraded, and the same system - with the same average number
of molecules - connected to a buffer. We modeled these systems using Ito
stochastic differential equations in discrete time, as they allow
straightforward analytical developments. In general, when the molecules in the
system and the buffer are positively correlated, the variance on the number of
molecules in the system is found to decrease compared to the equivalent system
without a buffer. Only buffers that are too noisy by themselves tend to
increase the noise in the main system. We tested this result on two model
cases, in which the system and the buffer contain proteins in their active and
inactive state, or protein monomers and homodimers. We found that in the second
test case, where the interconversion terms are non-linear in the number of
molecules, the noise reduction is much more pronounced; it reaches up to 20%
reduction of the Fano factor with the parameter values tested in numerical
simulations on an unperturbed birth-death model. We extended our analysis to
two arbitrary interconnected systems.


  Successful sequencing experiments require judicious sample selection.
However, this selection must often be performed on the basis of limited
preliminary data. Predicting the statistical properties of the final sample
based on preliminary data can be challenging, because numerous uncertain model
assumptions may be involved. Here, we ask whether we can predict ``omics"
variation across many samples by sequencing only a fraction of them. In the
infinite-genome limit, we find that a pilot study sequencing $5\%$ of a
population is sufficient to predict the number of genetic variants in the
entire population within $6\%$ of the correct value, using an estimator
agnostic to demography, selection, or population structure. To reach similar
accuracy in a finite genome with millions of polymorphisms, the pilot study
would require about $15\%$ of the population. We present computationally
efficient jackknife and linear programming methods that exhibit substantially
less bias than the state of the art when applied to simulated data and
sub-sampled 1000 Genomes Project data. Extrapolating based on the NHLBI Exome
Sequencing Project data, we predict that $7.2\%$ of sites in the capture region
would be variable in a sample of $50,000$ African-Americans, and $8.8\%$ in a
European sample of equal size. Finally, we show how the linear programming
method can also predict discovery rates of various genomic features, such as
the number of transcription factor binding sites across different cell types.


  Heritable differences in gene expression between individuals are an important
source of phenotypic variation. The question of how closely the effects of
genetic variation on protein levels mirror those on mRNA levels remains open.
Here, we addressed this question by using ribosome profiling to examine how
genetic differences between two strains of the yeast S. cerevisiae affect
translation. Strain differences in translation were observed for hundreds of
genes. Allele specific measurements in the diploid hybrid between the two
strains revealed roughly half as many cis-acting effects on translation as were
observed for mRNA levels. In both the parents and the hybrid, most effects on
translation were of small magnitude, such that the direction of an mRNA
difference was typically reflected in a concordant footprint difference. The
relative importance of cis and trans acting variation on footprint levels was
similar to that for mRNA levels. There was a tendency for translation to cause
larger footprint differences than expected given the respective mRNA
differences. This is in contrast to translational differences between yeast
species that have been reported to more often oppose than reinforce mRNA
differences. Finally, we catalogued instances of premature translation
termination in the two yeast strains and also found several instances where
erroneous reference gene annotations lead to apparent nonsense mutations that
in fact reside outside of the translated gene body. Overall, genetic influences
on translation subtly modulate gene expression differences, and translation
does not create strong discrepancies between genetic influences on mRNA and
protein levels.


  Traditional biomedical approaches treat diseases in isolation, but the
importance of synergistic disease interactions is now recognized. As a first
step we present and analyze a simple coinfection model for two diseases
affecting simultaneously a population. The host population is affected by the
\emph{primary disease}, a long-term infection whose dynamics is described by a
SIS model with demography, which facilitates individuals acquiring a second
disease, \emph{secondary (or \emph{opportunistic}) disease}. The secondary
disease is instead a short-term infection affecting only the primary-infected
individuals. Its dynamics is also represented by a SIS model with no
demography. To distinguish between short and long-term infection the complete
model is written as a two time scales system. The primary disease acts at the
slow time scale while the secondary disease does at the fast one, allowing a
dimension reduction of the system and making its analysis tractable. We show
that an opportunistic disease outbreak might change drastically the outcome of
the primary epidemic process, although it does among the outcomes allowed by
the primary disease. We have found situations in which either acting on the
opportunistic disease transmission or recovery rates or controlling the
susceptible and infected population size allow to eradicate/promote disease
endemicity.


  Cooperative behaviors are defined as the production of common goods
benefitting all members of the community at the producer's cost. They could
seem to be in contradiction with natural selection, as non-cooperators have an
increased fitness compared to cooperators. Understanding the emergence of
cooperation has necessitated the development of concepts and models (inclusive
fitness, multilevel selection, ...) attributing deterministic advantages to
this behavior. In contrast to these models, we show here that cooperative
behaviors can emerge by taking into account the stochastic nature of
evolutionary dynamics : when cooperative behaviors increase the carrying
capacity of the habitat, they also increase the genetic drift against
non-cooperators. Using the Wright-Fisher models of population genetics, we
compute exactly this increased genetic drift and its consequences on the
fixation probability of both types of individuals. This computation leads to a
simple criterion: cooperative behavior dominates when the relative increase in
carrying capacity of the habitat caused by cooperators is higher than the
selection pressure against them. This is a purely stochastic effect with no
deterministic interpretation.


  Cells generally convert external nutrient resources to support metabolismand
growth. Understanding the thermodynamic efficiency of this conversion is
essential to determine the general characteristics of cellular growth. Using a
simple protocell model with catalytic reaction dynamics to synthesize the
necessary enzyme and membrane components from nutrients, the entropy production
per unit cell-volume growth is calculated analytically and numerically based on
the rate equation for chemical kinetics and linear non-equilibrium
thermodynamics. The minimal entropy production per unit cell growth is found to
be achieved at a non-zero nutrient uptake rate, rather than at a quasi-static
limit as in the standard Carnot engine. This difference appears because the
equilibration mediated by the enzyme exists only within cells that grow through
enzyme and membrane synthesis. Optimal nutrient uptake is also confirmed by
protocell models with many chemical components synthesized through a catalytic
reaction network. The possible relevance of the identified optimal uptake to
optimal yield for cellular growth is also discussed.


  The Rock-Paper-Scissors (RPS) game and its generalizations with ${\cal S}>3$
species are well studied models for cyclically interacting populations. Four
is, however, the minimum number of species that, by allowing other interactions
beyond the single, cyclic loop, breaks both the full intransitivity of the food
graph and the one predator, one prey symmetry. L\"utz {\it et al} (J. Theor.
Biol. {\bf 317} (2013) 286) have shown the existence, on a square lattice, of
two distinct phases, with either four or three coexisting species. In both
phases, each agent is eventually replaced by one of its predators but these
strategy oscillations remain localized as long as the interactions are short
ranged. Distant regions may be either out of phase or cycling through different
food web subloops (if any). Here we show that upon replacing a minimum fraction
$Q_c$ of the short range interactions by long range ones, there is a Hopf
bifurcation and global oscillations become stable. Surprisingly, to build such
long distance, global synchronization, the four species coexistence phase
requires less long range interactions than the three species phase, while one
would naively expect the contrary. Moreover, deviations from highly homogeneous
conditions ($\chi=0$ or 1) increase $Q_c$ and the more heterogeneous is the
food web, the harder the synchronization is. By further increasing $Q$, while
the three species phase remains stable, the four species one has a transition
to an absorbing, single species state. The existence of a phase with global
oscillations for ${\cal S}>3$, when the interaction graph has multiple subloops
and several possible local cycles, lead to the conjecture that global
oscillations are a general characteristic, even for large, realistic food webs.


  Nanoparticles (NPs) engineered for biomedical applications are meant to be in
contact with protein-rich physiological fluids. These proteins are usually
adsorbed onto the NP surface, forming a swaddling layer called protein corona
that influences cell internalization. We present a study on protein adsorption
onto different magnetic NPs (MNPs) when immersed in cell culture medium, and
how these changes affect the cellular uptake. Two colloids with magnetite cores
of 25 nm, same hydrodynamic size and opposite surface charge were in situ
coated with (a) positive polyethyleneimine (PEI-MNPs) and (b) negative
poly(acrylic acid) (PAA-MNPs). After few minutes of incubation in cell culture
medium the wrapping of the MNPs by protein adsorption resulted in a 5-fold size
increase. After 24 h of incubation large MNP-protein aggregates with
hydrodynamic sizes 1500 to 3000 nm (PAA-MNPs and PEI-MNPs respectively) were
observed. Each cluster contained an estimated number of magnetic cores between
450 and 1000, indicating the formation of large aggregates with a "plum
pudding" structure of MNPs embedded into a protein network of negative surface
charge irrespective of the MNP_core charge. We demonstrated that PEI-MNPs are
incorporated in much larger amounts than the PAA-MNPs units. Quantitative
analysis showed that SH-SY5Y cells can incorporate 100 per cent of the added
PEI-MNPs up to about 100 pg per cell, whereas for PAA-MNPs the uptake was less
than 50 percent. The final cellular distribution showed also notable
differences regarding partial attachment to the cell membrane. These results
highlight the need to characterize the final properties of MNPs after protein
adsorption in biological media, and demonstrate the impact of these properties
on the internalization mechanisms in neural cells.


  Here we report on the use of double-nanohole (DNH) optical tweezers as a
label-free and free-solution single-molecule probe for protein-DNA
interactions. Using this approach, we demonstrate the unzipping of individual
10 base pair DNA-hairpins, and quantify how tumor suppressor p53 protein delays
the unzipping. From the Arrhenius behavior, we find the energy barrier to
unzipping introduced by p53 to be $2\times 10^{-20}$ J, whereas cys135ser
mutant p53 does not show suppression of unzipping, which gives clues to its
functional inability to suppress tumor growth. This transformative approach to
single molecule analysis allows for ultra-sensitive detection and
quantification of protein-DNA interactions to revolutionize the fight against
genetic diseases.


  In this paper, we propose a mathematical model for HIV-1 infection with
intracellular delay. The model examines a viral-therapy for controlling
infections through recombining HIV-1 virus with a genetically modified virus.
For this model, the basic reproduction number $\mathcal{R}_0$ are identified
and its threshold properties are discussed. When $\mathcal{R}_0 < 1$, the
infection-free equilibrium $E_0$ is globally asymptotically stable. When
$\mathcal{R}_0 > 1$, $E_0$ becomes unstable and there occurs the
single-infection equilibrium $E_s$, and $E_0$ and $E_s$ exchange their
stability at the transcritical point $\mathcal{R}_0 =1$. If $1< \mathcal{R}_0 <
R_1$, where $R_1$ is a positive constant explicitly depending on the model
parameters, $E_s$ is globally asymptotically stable, while when $\mathcal{R}_0
> R_1$, $E_s$ loses its stability to the double-infection equilibrium $E_d$.
There exist a constant $R_2$ such that $E_d$ is asymptotically stable if
$R_1<\mathcal R_0 < R_2$, and $E_s$ and $E_d$ exchange their stability at the
transcritical point $\mathcal{R}_0 =R_1$. We use one numerical example to
determine the largest range of $\mathcal R_0$ for the local stability of $E_d$
and existence of Hopf bifurcation. Some simulations are performed to support
the theoretical results. These results show that the delay plays an important
role in determining the dynamic behaviour of the system. In the normal range of
values, the delay may change the dynamic behaviour quantitatively, such as
greatly reducing the amplitudes of oscillations, or even qualitatively changes
the dynamical behaviour such as revoking oscillating solutions to equilibrium
solutions. This suggests that the delay is a very important fact which should
not be missed in HIV-1 modelling.


  Motivation: Recent advances in high-throughput sequencing (HTS) have made it
possible to monitor genomes in great detail. New experiments not only use HTS
to measure genomic features at one time point but to monitor them changing over
time with the aim of identifying significant changes in their abundance. In
population genetics, for example, allele frequencies are monitored over time to
detect significant frequency changes that indicate selection pressures.
Previous attempts at analysing data from HTS experiments have been limited as
they could not simultaneously include data at intermediate time points,
replicate experiments and sources of uncertainty specific to HTS such as
sequencing depth.
  Results: We present the beta-binomial Gaussian process (BBGP) model for
ranking features with significant non-random variation in abundance over time.
The features are assumed to represent proportions, such as proportion of an
alternative allele in a population. We use the beta-binomial model to capture
the uncertainty arising from finite sequencing depth and combine it with a
Gaussian process model over the time series. In simulations that mimic the
features of experimental evolution data, the proposed method clearly
outperforms classical testing in average precision of finding selected alleles.
We also present simulations exploring different experimental design choices and
results on real data from Drosophila experimental evolution experiment in
temperature adaptation.
  Availability: R software implementing the test is available at
https://github.com/handetopa/BBGP


  Mammalian brains span about 4 orders of magnitude in cortical volume and have
to operate in different environments that require diverse behavioral skills.
Despite these geometric and behavioral diversities, the examination of cerebral
cortex across species reveals that it contains a substantial number of
conserved characteristics that are associated with neuroanatomy and metabolism,
i.e. with neuronal connectivity and function. Some of these cortical constants
or invariants have been known for a long time but not sufficiently appreciated,
and others were only recently discovered. The focus of this review is to
present the cortical invariants and discuss their role in the efficient
information processing. Global conservation in neuroanatomy and metabolism, as
well as their correlated regional and developmental variability suggest that
these two parallel systems are mutually coupled. It is argued that energetic
constraint on cortical organization can be strong if cerebral blood supplied is
either below or above a certain level, and it is rather soft otherwise.
Moreover, because maximization or minimization of parameters associated with
cortical connectivity, function and cost often leads to conflicts in design, it
is argued that the architecture of the cerebral cortex is a result of
structural and functional compromises.


  We present a new, nucleotide-level model for RNA, oxRNA, based on the
coarse-graining methodology recently developed for the oxDNA model of DNA. The
model is designed to reproduce structural, mechanical and thermodynamic
properties of RNA, and the coarse-graining level aims to retain the relevant
physics for RNA hybridization and the structure of single- and double-stranded
RNA. In order to explore its strengths and weaknesses, we test the model in a
range of nanotechnological and biological settings. Applications explored
include the folding thermodynamics of a pseudoknot, the formation of a kissing
loop complex, the structure of a hexagonal RNA nanoring, and the unzipping of a
hairpin motif. We argue that the model can be used for efficient simulations of
the structure of systems with thousands of base pairs, and for the assembly of
systems of up to hundreds of base pairs. The source code implementing the model
is released for public use.


  Massive research efforts are now underway to develop a cure for HIV
infection, allowing patients to discontinue lifelong combination antiretroviral
therapy (ART). New latency-reversing agents (LRAs) may be able to purge the
persistent reservoir of latent virus in resting memory CD4+ T cells, but the
degree of reservoir reduction needed for cure remains unknown. Here we use a
stochastic model of infection dynamics to estimate the efficacy of LRA needed
to prevent viral rebound after ART interruption. We incorporate clinical data
to estimate population-level parameter distributions and outcomes. Our findings
suggest that approximately 2,000-fold reductions are required to permit a
majority of patients to interrupt ART for one year without rebound and that
rebound may occur suddenly after multiple years. Greater than 10,000-fold
reductions may be required to prevent rebound altogether. Our results predict
large variation in rebound times following LRA therapy, which will complicate
clinical management. This model provides benchmarks for moving LRAs from the
lab to the clinic and can aid in the design and interpretation of clinical
trials. These results also apply to other interventions to reduce the latent
reservoir and can explain the observed return of viremia after months of
apparent cure in recent bone marrow transplant recipients and an
immediately-treated neonate.


  Efforts to reconstruct phylogenetic trees and understand evolutionary
processes depend fundamentally on stochastic models of speciation and mutation.
The simplest continuous-time model for speciation in phylogenetic trees is the
Yule process, in which new species are "born" from existing lineages at a
constant rate. Recent work has illuminated some of the structural properties of
Yule trees, but it remains mostly unknown how these properties affect sequence
and trait patterns observed at the tips of the phylogenetic tree. Understanding
the interplay between speciation and mutation under simple models of evolution
is essential for deriving valid phylogenetic inference methods and gives
insight into the optimal design of phylogenetic studies. In this work, we
derive the probability distribution of interspecies covariance under Brownian
motion and Ornstein-Uhlenbeck models of phenotypic change on a Yule tree. We
compute the probability distribution of the number of mutations shared between
two randomly chosen taxa in a Yule tree under discrete Markov mutation models.
Our results suggest summary measures of phylogenetic information content,
illuminate the correlation between site patterns in sequences or traits of
related organisms, and provide heuristics for experimental design and
reconstruction of phylogenetic trees.


  Motivation: Recombination rates vary considerably at the fine scale within
mammalian genomes, with the majority of recombination occurring within hotspots
of ~2 kb in width. We present a method for inferring the location of
recombination hotspots from patterns of linkage disequilibrium within samples
of population genetic data. Results: Using simulations, we show that our method
has hotspot detection power of approximately 50-60%, but depending on the
magnitude of the hotspot. The false positive rate is between 0.24 and 0.56
false positives per Mb for data typical of humans. Availability:
http://github.com/auton1/LDhot


  Advances in experimental techniques and in theoretical models have improved
our understanding of protein crystallization. But they have also left open
questions regarding the protein phase behavior and self-assembly kinetics, such
as why (nearly) identical crystallization conditions sometimes result in the
formation of different crystal forms. Here, we develop a patchy particle model
with competing sets of patches that provides a microscopic explanation of this
phenomenon. We identify different regimes in which one or two crystal forms can
coexist with a low-density fluid. Using analytical approximations, we extend
our findings to different crystal phases, providing a general framework for
treating protein crystallization when multiple crystal forms compete. Our
results also suggest different experimental routes for targeting a specific
crystal form, and for reducing the dynamical competition between the two forms,
thus facilitating protein crystal assembly.


  In order to model real ecological systems one has to consider many species
that interact in complex ways. However, most of the recent theoretical studies
have been restricted to few species systems with rather trivial interactions.
The few studies dealing with larger number of species and/or more complex
interaction schemes are mostly restricted to numerical explorations. In this
paper we determine, starting from the deterministic mean-field rate equations,
for large classes of systems the space of coexistence fixed points at which
biodiversity is maximal. For systems with a single coexistence fixed point we
derive complex Ginzburg-Landau equations that allow to describe space-time
pattern realized in two space dimensions. For selected cases we compare the
theoretical predictions with the pattern observed in numerical simulations.


  A new model to investigate environmental effects of genetically
distinguishable predators is presented. The Holling type II response function,
modelling feeding satiation, leads to persistent system's oscillations, as in
classical population models. An almost complete classification of the cases
arising in the Routh-Hurwitz stability conditions mathematically characterizes
the paper. It is instrumental as a guideline in the numerical experiments
leading to the findings on the limit cycles. This result extends what found in
an earlier parallel investigation containing a standard bilinear response
function.


  In this paper we consider a competition system in which two diseases spread
by contact. We characterize the system behavior, establishing that only some
configurations are possible. In particular we discover that coexistence of the
two strains is not possible, under the assumptions of the model. A number of
transcritical bifurcations relate the more relevant system's equilibria.
Bistability is shown between a situation in which only the disease-unaffected
population thrives and another one containing only the second population with
endemic disease. An accurate computation of the separating surface of the
basins of attraction of these two mutually exclusive equilibria is obtained via
novel results in approximation theory.


  Models of coordinated behavior of populations living in the same environment
are introduced for the cases when they either compete with each other, or they
both gain by mutual interactions, or finally when one hunts the other one. The
equilibria of the systems are analysed, showing that in some cases the
populations may both disappear. Coexistence leads to global asymptotic
stability for symbiotic populations, or to Hopf bifurcations for predator-prey
systems. Finally, a new very interesting phenomenon is discovered in one of
these rather simple models. Indeed tristability may be achieved in the
competition case, for which competitive exclusion is allowed to occur together
with populations coexistence.


  In this paper we consider a simple two-patch model in which a population
affected by a disease can freely move. We assume that the capacity of the
interconnected paths is limited, and thereby influencing the migration rates.
Possible habitat disruptions due to human activities or natural events are
accounted for. The demographic assumptions prevent the ecosystem to be wiped
out, and the disease remains endemic in both populated patches at a stable
equilibrium, but possibly also with an oscillatory behavior in the case of
unidirectional migrations. Interestingly, if infected cannot migrate, it is
possible that one patch becomes disease-free. This fact could be exploited to
keep disease-free at least part of the population.


  The possible dynamics of an ecosystem with three interconnected patches among
which one population can migrate are investigated, paying attention to the
consequences that possible disruption in the communicating paths can cause.


  We continue the investigations of an ecosystem where a epidemic-affected
population can move between two connected patches, by considering what happens
to the system when the migration paths are interrupted in one direction, or
when the infected are not able to exert the effort for migrating into the other
patch.


  Ribonucleotide reductases are encoded by many viruses, but without other
enzymes of nucleotide metabolism of no obvious use. A look at the enzymes'
molecular properties and their possible mutator action may give clues.


  Social insects provide an excellent platform to investigate flow of
information in regulatory systems since their successful social organization is
essentially achieved by effective information transfer through complex
connectivity patterns among the colony members. Network representation of such
behavioural interactions offers a powerful tool for structural as well as
dynamical analysis of the underlying regulatory systems. In this paper, we
focus on the dominance interaction networks in the tropical social wasp
\textit{Ropalidia marginata} - a species where behavioural observations
indicate that such interactions are principally responsible for the transfer of
information between individuals about their colony needs, resulting in a
regulation of their own activities. Our research reveals that the dominance
networks of \textit{R. marginata} are structurally similar to a class of
naturally evolved information processing networks, a fact confirmed also by the
predominance of a specific substructure - the `feed-forward loop' - a key
functional component in many other information transfer networks. The dynamical
analysis through Boolean modeling confirms that the networks are sufficiently
stable under small fluctuations and yet capable of more efficient information
transfer compared to their randomized counterparts. Our results suggest the
involvement of a common structural design principle in different biological
regulatory systems and a possible similarity with respect to the effect of
selection on the organization levels of such systems. The findings are also
consistent with the hypothesis that dominance behaviour has been shaped by
natural selection to co-opt the information transfer process in such social
insect species, in addition to its primal function of mediation of reproductive
competition in the colony.


  Limited dispersal distance, whether associated with vegetative growth or
localized reproduction, induces spatial clustering and, in turn, focuses
ecological interactions at the neighborhood scale. In particular, most invasive
plants are clonal, cluster through vegetative propagation, and compete locally.
Dispersal limitation implies that invasive spread occurs as advance of an
ecological interface between invader and resident species. Interspecific
competition along the interface produces random variation in the extent of
invasive growth. Development of these random fluctuations, termed stochastic
roughening, will often structure the interface as a self-affine fractal; a
series of power-law scaling relationships follows as a result. For a diverse
array of local growth processes exhibiting both forward and lateral
propagation, the extent of invader advance becomes spatially correlated along
the interface, and the width of the interface (the area where invader and
resident compete directly) increases as a power function of time. Once
roughening equilibrates statistically, interface width and the location of the
most advanced invader (the "front-runner") beyond the mean incursion should
both increase as a power function of interface length. To test these
predictions, we let white clover (Trifolium repens) invade ryegrass (Lolium
perenne) experimentally. Spatial correlation developed as anticipated, and both
interface width and the front-runner's lead scaled as a power law of length.
However, the scaling exponents differed, likely a consequence of clover's
growth morphology. The theory of kinetic roughening offers a new framework for
understanding causes and consequences of spatial pattern in between-species
interaction, and indicates when interface measures at a local scale predict
properties of an invasive front at extended spatial scales.


  We present a palette of brightly fluorescent genetically encoded voltage
indicators (GEVIs) with excitation and emission peaks spanning the visible
spectrum, sensitivities from 6 - 10% Delta F/F per 100 mV, and half-maximal
response times from 1 - 7 ms. A fluorescent protein is fused to an
Archaerhodopsin-derived voltage sensor. Voltage-induced shifts in the
absorption spectrum of the rhodopsin lead to voltage-dependent nonradiative
quenching of the appended fluorescent protein. Through a library screen, we
identified linkers and fluorescent protein combinations which reported neuronal
action potentials in cultured rat hippocampal neurons with a single-trial
signal-to-noise ratio from 6.6 to 11.6 in a 1 kHz imaging bandwidth at modest
illumination intensity. The freedom to choose a voltage indicator from an array
of colors facilitates multicolor voltage imaging, as well as combination with
other optical reporters and optogenetic actuators.


  The spreading of evolutionary novelties across populations is the central
element of adaptation. Unless population are well-mixed (like bacteria in a
shaken test tube), the spreading dynamics not only depends on fitness
differences but also on the dispersal behavior of the species. Spreading at a
constant speed is generally predicted when dispersal is sufficiently
short-ranged. However, the case of long-range dispersal is unresolved: While it
is clear that even rare long-range jumps can lead to a drastic speedup, it has
been difficult to quantify the ensuing stochastic growth process. Yet such
knowledge is indispensable to reveal general laws for the spread of modern
human epidemics, which is greatly accelerated by aviation. We present a simple
iterative scaling approximation supported by simulations and rigorous bounds
that accurately predicts evolutionary spread for broad distributions of long
distance dispersal. In contrast to the exponential laws predicted by
deterministic "mean-field" approximations, we show that the asymptotic growth
is either according to a power-law or a stretched exponential, depending on the
tails of the dispersal kernel. More importantly, we provide a full
time-dependent description of the convergence to the asymptotic behavior which
can be anomalously slow and is needed even for long times. Our results also
apply to spreading dynamics on networks with a spectrum of long-range links
under certain conditions on the probabilities of long distance travel and are
thus relevant for the spread of epidemics.


  Mass spectrometry-based metabolomic analysis depends upon the identification
of spectral peaks by their mass and retention time. Statistical analysis that
follows the identification currently relies on one main peak of each compound.
However, a compound present in the sample typically produces several spectral
peaks due to its isotopic properties and the ionization process of the mass
spectrometer device. In this work, we investigate the extent to which these
additional peaks can be used to increase the statistical strength of
differential analysis.
  We present a Bayesian approach for integrating data of multiple detected
peaks that come from one compound. We demonstrate the approach through a
simulated experiment and validate it on ultra performance liquid
chromatography-mass spectrometry (UPLC-MS) experiments for metabolomics and
lipidomics. Peaks that are likely to be associated with one compound can be
clustered by the similarity of their chromatographic shape. Changes of
concentration between sample groups can be inferred more accurately when
multiple peaks are available.
  When the sample-size is limited, the proposed multi-peak approach improves
the accuracy at inferring covariate effects. An R implementation, data and the
supplementary material are available at
http://research.ics.aalto.fi/mi/software/peakANOVA/ .


  Motivation: Identification of flexible regions of protein structures is
important for understanding of their biological functions. Recently, we have
developed a fast approach for predicting protein structure fluctuations from a
single protein model: the CABS-flex. CABS-flex was shown to be an efficient
alternative to conventional all-atom molecular dynamics (MD). In this work, we
evaluate CABS-flex and MD predictions by comparison with protein structural
variations within NMR ensembles.
  Results: Based on a benchmark set of 140 proteins, we show that the relative
fluctuations of protein residues obtained from CABS-flex are well correlated to
those of NMR ensembles. On average, this correlation is stronger than that
between MD and NMR ensembles. In conclusion, CABS-flex is useful and
complementary to MD in predicting of protein regions that undergo
conformational changes and the extent of such changes.


  Recent improvements in high-throughput genotyping and sequencing technologies
have afforded the collection of massive, genome-wide datasets of DNA
information from hundreds of thousands of individuals. These datasets, in turn,
provide unprecedented opportunities to reconstruct the history of human
populations and detect genotype-phenotype association. Recently developed
computational methods can identify long-range chromosomal segments that are
identical across samples, and have been transmitted from common ancestors that
lived tens to hundreds of generations in the past. These segments reveal
genealogical relationships that are typically unknown to the carrying
individuals. In this work, we demonstrate that such identical-by-descent (IBD)
segments are informative about a number of relevant population genetics
features: they enable the inference of details about past population size
fluctuations, migration events, and they carry the genomic signature of natural
selection. We derive a mathematical model, based on coalescent theory, that
allows for a quantitative description of IBD sharing across purportedly
unrelated individuals, and develop inference procedures for the reconstruction
of recent demographic events, where classical methodologies are statistically
underpowered. We analyze IBD sharing in several contemporary human populations,
including representative communities of the Jewish Diaspora, Kenyan Maasai
samples, and individuals from several Dutch provinces, in all cases retrieving
evidence of fine-scale demographic events from recent history. Finally, we
expand the presented model to describe distributions for those sites in IBD
shared segments that harbor mutation events, showing how these may be used for
the inference of mutation rates in humans and other species.


  Phylogenetic stochastic mapping is a method for reconstructing the history of
trait changes on a phylogenetic tree relating species/organisms carrying the
trait. State-of-the-art methods assume that the trait evolves according to a
continuous-time Markov chain (CTMC) and work well for small state spaces. The
computations slow down considerably for larger state spaces (e.g. space of
codons), because current methodology relies on exponentiating CTMC
infinitesimal rate matrices -- an operation whose computational complexity
grows as the size of the CTMC state space cubed. In this work, we introduce a
new approach, based on a CTMC technique called uniformization, that does not
use matrix exponentiation for phylogenetic stochastic mapping. Our method is
based on a new Markov chain Monte Carlo (MCMC) algorithm that targets the
distribution of trait histories conditional on the trait data observed at the
tips of the tree. The computational complexity of our MCMC method grows as the
size of the CTMC state space squared. Moreover, in contrast to competing matrix
exponentiation methods, if the rate matrix is sparse, we can leverage this
sparsity and increase the computational efficiency of our algorithm further.
Using simulated data, we illustrate advantages of our MCMC algorithm and
investigate how large the state space needs to be for our method to outperform
matrix exponentiation approaches. We show that even on the moderately large
state space of codons our MCMC method can be significantly faster than
currently used matrix exponentiation methods.


  We prove the nonlinear local stability of Dirac masses for a kinetic model of
alignment of particles on the unit sphere, each point of the unit sphere
representing a direction. A population concentrated in a Dirac mass then
corresponds to the global alignment of all individuals. The main difficulty of
this model is the lack of conserved quantities and the absence of an energy
that would decrease for any initial condition. We overcome this difficulty
thanks to a functional which is decreasing in time in a neighborhood of any
Dirac mass (in the sense of the Wasserstein distance). The results are then
extended to the case where the unit sphere is replaced by a general Riemannian
manifold.


  Several living systems have been examined for their apparent optimization of
structure and function for quantum behavior at biological length scales.
Orthodox type II endonucleases, the largest class of restriction enzymes,
recognize four-to-eight base pair sequences of palindromic DNA, cut both
strands symmetrically, and act without an external metabolite such as ATP.
While it is known that these enzymes induce strand breaks by attacking
phosphodiester bonds, what remains unclear is the mechanism by which cutting
occurs in concert at the catalytic centers. Previous studies indicate the
primacy of intimate DNA contacts made by the specifically bound enzyme in
coordinating the two synchronized cuts. We propose that collective electronic
behavior in the DNA helix generates coherent oscillations, quantized through
boundary conditions imposed by the endonuclease, that provide the energy
required to break two phosphodiester bonds. Such quanta may be preserved in the
presence of thermal noise and electromagnetic interference through decoherence
shielding, the specific complex's exclusion of water and ions surrounding the
helix. Clamping energy imparted by the enzyme decoherence shield is comparable
with zero-point modes of the dipole-dipole oscillations in the DNA recognition
sequence. The palindromic mirror symmetry of this sequence should conserve
parity during the process. Experimental data corroborate that symmetric
bond-breaking ceases when the symmetry of the endonuclease complex is violated,
or when environmental parameters are perturbed far from biological optima.
Persistent correlation between states in DNA sequence across spatial
separations of any length--a characteristic signature of quantum
entanglement--may be explained by such a physical mechanism.


  We study the parameter sensitivity of hetero-polymeric DNA within the purview
of DNA breathing dynamics. The degree of correlation between the mean bubble
size and the model parameters are estimated for this purpose for three
different DNA sequences. The analysis leads us to a better understanding of the
sequence dependent nature of the breathing dynamics of hetero-polymeric DNA.
Out of the fourteen model parameters for DNA stability in the statistical
Poland-Scheraga approach, the hydrogen bond interaction
$\epsilon_{hb}(\mathtt{AT})$ for an $\mathtt{AT}$ base pair and the ring factor
$\xi$ turn out to be the most sensitive parameters. In addition, the stacking
interaction $\epsilon_{st}(\mathtt{TA}-\mathtt{TA})$ for an
$\mathtt{TA}-\mathtt{TA}$ nearest neighbor pair of base-pairs is found to be
the most sensitive one among all stacking interactions. Moreover, we also
establish that the nature of stacking interaction has a deciding effect on the
DNA breathing dynamics, not the number of times a particular stacking
interaction appears in a sequence. We show that the sensitivity analysis can be
used as an effective measure to guide a stochastic optimization technique to
find the kinetic rate constants related to the dynamics as opposed to the case
where the rate constants are measured using the conventional unbiased way of
optimization.


  Proteomics will celebrate its 20th year in 2014. In this relatively short
period of time, it has invaded most areas of biology and its use will probably
continue to spread in the future. These two decades have seen a considerable
increase in the speed and sensitivity of protein identification and
characterization, even from complex samples. Indeed, what was a challenge
twenty years ago is now little more than a daily routine. Although not
completely over, the technological challenge now makes room to another
challenge, which is the best possible appraisal and exploitation of proteomic
data to draw the best possible conclusions from a biological point of view. The
point developed in this paper is that proteomic data are almost always
fragmentary. This means in turn that although better than an mRNA level, a
protein level is often insufficient to draw a valid conclusion from a
biological point of view, especially in a world where PTMs play such an
important role. This means in turn that transformation of proteomic data into
biological data requires an important intermediate layer of functional
validation, i.e. not merely the confirmation of protein abundance changes by
other methods, but a functional appraisal of the biological consequences of the
protein level changes highlighted by the proteomic screens.


  RNA secondary structure prediction and classification are two important
problems in the field of RNA biology. Here, we propose a new permutation based
approach to create logical non-disjoint clusters of different secondary
structures of a single class or type. Many different types of techniques exist
to classify RNA secondary structure data but none of them have ever used
permutation based approach which is very simple and yet powerful. We have
written a small JAVA program to generate permutation, apply our algorithm on
those permutations and analyze the data and create different logical clusters.
We believe that these clusters can be utilized to untangle the mystery of RNA
secondary structure and analyze the development patterns of unknown RNA.


  Mutualism is a biological interaction mutually beneficial for both species
involved, such as the interaction between plants and their pollinators. Real
mutualistic communities can be understood as weighted bipartite networks and
they present a nested structure and truncated power law degree and strength
distributions. We present a novel link aggregation model that works on a
strength-preferential attachment rule based on the Individual Neutrality
hypothesis. The model generates mutualistic networks with emergent nestedness
and truncated distributions. We provide some analytical results and compare the
simulated and empirical network topology. Upon further improving the shape of
the distributions, we have also studied the role of forbidden interactions on
the model and found that the inclusion of forbidden links does not prevent for
the appearance of super-generalist species. A Python script with the model
algorithms is available.


  We report on the observation of a change in the bend angle of an RNA kissing
complex upon Rop binding using single-molecular-pair FRET. The angular
relationship between the dyes, rather than the distance between them, is shown
to be responsible for the observed change in energy transfer. It has long been
thought that Rop increases the bend angle of the R1inv-R2inv complex upon
binding, but this has never been directly observed. In contrast, we find an
increase in FRET upon the addition of Rop that is shown via modeling to be
consistent with a decrease in the bend angle of the complex of
$-15^{\circ}\pm7^{\circ}$. The model predicts FRET from dye trajectories
generated using molecular dynamics simulations of Cy3 and Cy5 attached to $5'$
terminal cytosine or guanosine on RNA. While FRET is commonly used to observe
global changes in molecular structure attributed to changes in the distance
between dyes, it is rarely, if ever, used to elucidate angular changes. Subtle
global changes in molecular structure upon binding are generally difficult to
discern using NMR or crystallography, but FRET is well suited to observe these
changes because of its sensitivity at interdye distance around the F\"orster
radius, $R_F\approx5$ nm. While FRET is often referred to and used as a
"molecular ruler" for distances around $R_F$, for dye pairs that have minimal
rotational freedom, FRET can also be used to observe changes in structure for
which there is no significant change in distance between the dyes. This is the
case for the R1inv-R2inv kissing complex studied here. This complex is derived
from the RNA I - RNA II complex in E. coli. RNA II is a primer for replication
of the ColE1 plasmid, its function is modulated by interaction with RNA II. Rop
is known to stabilize the complex, and it is also known to bind kissing
complexes in a structure, but not sequence, dependent fashion.


  Recent studies have shown that human populations have experienced a complex
demographic history, including a recent epoch of rapid population growth that
led to an excess in the proportion of rare genetic variants in humans today.
This excess can impact the burden of private mutations for each individual,
defined here as the proportion of heterozygous variants in each newly sequenced
individual that are novel compared to another large sample of sequenced
individuals. We calculated the burden of private mutations predicted by
different demographic models, and compared with empirical estimates based on
data from the NHLBI Exome Sequencing Project and data from the Neutral Regions
(NR) dataset. We observed a significant excess in the proportion of private
mutations in the empirical data compared with models of demographic history
without a recent epoch of population growth. Incorporating recent growth into
the model provides a much improved fit to empirical observations. This
phenomenon becomes more marked for larger sample sizes. The proportion of
private mutations is additionally increased by purifying selection, which
differentially affect mutations of different functional annotations. These
results have important implications to the design and analysis of
sequencing-based association studies of complex human disease as they pertain
to private and very rare variants.


  Recent emergence of next-generation DNA sequencing technology has enabled
acquisition of genetic information at unprecedented scales. In order to
determine the genetic blueprint of an organism, sequencing platforms typically
employ so-called shotgun sequencing strategy to oversample the target genome
with a library of relatively short overlapping reads. The order of nucleotides
in the reads is determined by processing the acquired noisy signals generated
by the sequencing instrument. Assembly of a genome from potentially erroneous
short reads is a computationally daunting task even in the scenario where a
reference genome exists. Errors and gaps in the reference, and perfect repeat
regions in the target, further render the assembly challenging and cause
inaccuracies. In this paper, we formulate the reference-guided sequence
assembly problem as the inference of the genome sequence on a bipartite graph
and solve it using a message-passing algorithm. The proposed algorithm can be
interpreted as the well-known classical belief propagation scheme under a
certain prior. Unlike existing state-of-the-art methods, the proposed algorithm
combines the information provided by the reads without needing to know
reliability of the short reads (so-called quality scores). Relation of the
message-passing algorithm to a provably convergent power iteration scheme is
discussed. To evaluate and benchmark the performance of the proposed technique,
we find an analytical expression for the probability of error of a genie-aided
maximum a posteriori (MAP) decision scheme. Results on both simulated and
experimental data demonstrate that the proposed message-passing algorithm
outperforms commonly used state-of-the-art tools, and it nearly achieves the
performance of the aforementioned MAP decision scheme.


  In this paper we investigate a susceptible-infected-susceptible (SIS)
epidemic model describing data dissemination in opportunistic networks with
heterogeneous setting of transmission parameters. We obtained the estimation of
the final epidemic size assuming that amount of data transferred between
network nodes possesses a Pareto distribution, implying scale-free properties.
In this context, more heterogeneity in susceptibility means the less severe
epidemic progression, and, on the contrary, more heterogeneity in infectivity
leads to more severe epidemics -- assuming that the other parameter (either
heterogeneity or susceptibility) stays fixed. The results are general enough
and can be useful in general epidemic theory for estimating the epidemic
progression for diseases with no significant acquired immunity -- in the cases
where Pareto distribution holds.