  Ancient mitochondrial DNA has been used in a wide variety of palaeontological
and archaeological studies, ranging from population dynamics of extinct species
to patterns of domestication. Most of these studies have traditionally been
based on the analysis of short fragments from the mitochondrial control region,
analysed using PCR coupled with Sanger sequencing. With the introduction of
high-throughput sequencing, as well as new enrichment technologies, the
recovery of full mitochondrial genomes (mitogenomes) from ancient specimens has
become significantly less complicated. Here we present a protocol to build
ancient extracts into Illumina high-throughput sequencing libraries, and
subsequent Agilent array-based capture to enrich for the desired mitogenome.
Both are based on previously published protocols, with the introduction of
several improvements aimed to increase the recovery of short DNA fragments,
while keeping the cost and effort requirements low. This protocol was designed
for enrichment of mitochondrial DNA in ancient or degraded samples. However,
the protocols can be easily adapted for using for building libraries for
shotgun-sequencing of whole genomes, or enrichment of other genomic regions.


  We introduce a new representation and feature extraction method for
biological sequences. Named bio-vectors (BioVec) to refer to biological
sequences in general with protein-vectors (ProtVec) for proteins (amino-acid
sequences) and gene-vectors (GeneVec) for gene sequences, this representation
can be widely used in applications of deep learning in proteomics and genomics.
In the present paper, we focus on protein-vectors that can be utilized in a
wide array of bioinformatics investigations such as family classification,
protein visualization, structure prediction, disordered protein identification,
and protein-protein interaction prediction. In this method, we adopt artificial
neural network approaches and represent a protein sequence with a single dense
n-dimensional vector. To evaluate this method, we apply it in classification of
324,018 protein sequences obtained from Swiss-Prot belonging to 7,027 protein
families, where an average family classification accuracy of 93%+-0.06% is
obtained, outperforming existing family classification methods. In addition, we
use ProtVec representation to predict disordered proteins from structured
proteins. Two databases of disordered sequences are used: the DisProt database
as well as a database featuring the disordered regions of nucleoporins rich
with phenylalanine-glycine repeats (FG-Nups). Using support vector machine
classifiers, FG-Nup sequences are distinguished from structured protein
sequences found in Protein Data Bank (PDB) with a 99.8% accuracy, and
unstructured DisProt sequences are differentiated from structured DisProt
sequences with 100.0% accuracy. These results indicate that by only providing
sequence data for various proteins into this model, accurate information about
protein structure can be determined.


  Genome rearrangements can be modeled as $k$-breaks, which break a genome at k
positions and glue the resulting fragments in a new order. In particular,
reversals, translocations, fusions, and fissions are modeled as $2$-breaks, and
transpositions are modeled as $3$-breaks. While $k$-break rearrangements for
$k>3$ have not been observed in evolution, they are used in cancer genomics to
model chromothripsis, a catastrophic event of multiple breakages happening
simultaneously in a genome. It is known that the $k$-break distance between two
genomes (i.e., the minimum number of $k$-breaks required to transform one
genome into the other) can be computed in terms of cycle lengths in the
breakpoint graph of these genomes.
  In the current work, we address the combinatorial problem of enumerating
genomes at a given $k$-break distance from a fixed unichromosomal genome. More
generally, we enumerate genome pairs, whose breakpoint graph has a given
distribution of cycle lengths. We further show how our enumeration can be used
for uniform sampling of random genomes at a given $k$-break distance, and
describe its connection to various combinatorial objects such as Bell
polynomials.


  The availability of a large number of assembled genomes opens the way to
study the evolution of syntenic character within a phylogenetic context. The
DeCo algorithm, recently introduced by B{\'e}rard et al. allows the computation
of parsimonious evolutionary scenarios for gene adjacencies, from pairs of
reconciled gene trees. Following the approach pioneered by Sturmfels and
Pachter, we describe how to modify the DeCo dynamic programming algorithm to
identify classes of cost schemes that generates similar parsimonious
evolutionary scenarios for gene adjacencies, as well as the robustness to
changes to the cost scheme of evolutionary events of the presence or absence of
specific ancestral gene adjacencies. We apply our method to six thousands
mammalian gene families, and show that computing the robustness to changes to
cost schemes provides new and interesting insights on the evolution of gene
adjacencies and the DeCo model.


  Structuring of DNA counterions around the double helix has been studied by
the molecular dynamics method. A DNA dodecamer d(CGCGAATTCGCG) in water
solution with the alkali metal counterions Na$^{+}$, K$^{+}$, and Cs$^{+}$ has
been simulated. The systems have been considered in the regimes without excess
salt and with different salts (0.5 M of NaCl, KCl or CsCl) added. The results
have showed that the Na$^{+}$ counterions interact with the phosphate groups
directly from outside of the double helix and via water molecules at the top
edge of DNA minor groove. The potassium ions are mostly localized in the
grooves of the double helix, and the cesium ions penetrate deeply inside the
minor groove being bonded directly to the atoms of nucleic bases. Due to the
electrostatic repulsion the chlorine ions tend to be localized at large
distances from the DNA polyanion, but some Cl$^{-}$ anions have been detected
near atomic groups of the double helix forming electrically neutral pairs with
counterions already condensed on DNA. The DNA sites, where counterions are
incorporated, are characterized by local changes of double helix structure. The
lifetime of Na$^{+}$ and K$^{+}$ in complex with DNA atomic groups is less than
0.5 ns, while in the case of the cesium ions it may reach several nanoseconds.
In this time scale, the Cs$^{+}$ counterions form a structured system of
charges in the DNA minor groove that can be considered as ionic lattice.


  Having a precise knowledge of the dispersal ability of a population in a
heterogeneous environment is of critical importance in agroecology and
conservation biology as it can provide management tools to limit the effects of
pests or to increase the survival of endangered species. In this paper, we
propose a mechanistic-statistical method to estimate space-dependent diffusion
parameters of spatially-explicit models based on stochastic differential
equations, using genetic data. Dividing the total population into
subpopulations corresponding to different habitat patches with known allele
frequencies, the expected proportions of individuals from each subpopulation at
each position is computed by solving a system of reaction-diffusion equations.
Modelling the capture and genotyping of the individuals with a statistical
approach, we derive a numerically tractable formula for the likelihood function
associated with the diffusion parameters.
  In a simulated environment made of three types of regions, each associated
with a different diffusion coefficient, we successfully estimate the diffusion
parameters with a maximum-likelihood approach. Although higher genetic
differentiation among subpopulations leads to more accurate estimations, once a
certain level of differentiation has been reached, the finite size of the
genotyped population becomes the limiting factor for accurate estimation.


  The first giant virus was identified in 2003 from a biofilm of an industrial
water-cooling tower in England. Later, numerous new giant viruses were found in
oceans and freshwater habitats, some of them having even 2,500 genes. We have
demonstrated their very likely presence in four soil samples taken from the
Kutch Desert (Gujarat, India). Here we describe a bioinformatics work-flow,
called the "Giant Virus Finder" that is capable to discover the very likely
presence of the genomes of giant viruses in metagenomic shotgun-sequenced
datasets. The new tool is applied to numerous hot and cold desert soil samples
as well as some tundra- and forest soils. We show that most of these samples
contain giant viruses, and especially many were found in the Antarctic dry
valleys. The results imply that giant viruses could be frequent not only in
aqueous habitats, but in a wide spectrum of soils on our planet.


  Bayesian approaches to variable selection have been widely used for
quantitative trait locus (QTL) mapping. The Markov chain Monte Carlo (MCMC)
algorithms for that aim are often difficult to be implemented for
high-dimensional variable selection problems, such as the ones arising in
epistatic analysis. Variational approximation is an alternative to MCMC, and
variational linear regression (VLR) is an effective solution for the variable
selection problems, but lacks accuracy in some QTL mapping problems where there
are many more variables than samples. In this paper, we propose an effective
method with aim to improve the accuracy of VLR in the case of above by
dynamically reducing components (variable or markers) with known effects (zero
or fixed). We show that the proposed method can greatly improve the accuracy of
VLR with little increase in computational cost. The method is compared with
several other variational methods used for QTL mapping, and simulation results
show that its performance is higher than those methods when applied in
high-dimensional cases.


  Many datasets exhibit a well-defined structure that can be exploited to
design faster search tools, but it is not always clear when such acceleration
is possible. Here, we introduce a framework for similarity search based on
characterizing a dataset's entropy and fractal dimension. We prove that
searching scales in time with metric entropy (number of covering hyperspheres),
if the fractal dimension of the dataset is low, and scales in space with the
sum of metric entropy and information-theoretic entropy (randomness of the
data). Using these ideas, we present accelerated versions of standard tools,
with no loss in specificity and little loss in sensitivity, for use in three
domains---high-throughput drug screening (Ammolite, 150x speedup), metagenomics
(MICA, 3.5x speedup of DIAMOND [3,700x BLASTX]), and protein structure search
(esFragBag, 10x speedup of FragBag). Our framework can be used to achieve
"compressive omics," and the general theory can be readily applied to data
science problems outside of biology.


  Coalescent histories provide lists of species tree branches on which gene
tree coalescences can take place, and their enumerative properties assist in
understanding the computational complexity of calculations central in the study
of gene trees and species trees. Here, we solve an enumerative problem left
open by Rosenberg (IEEE/ACM Transactions on Computational Biology and
Bioinformatics 10: 1253-1262, 2013) concerning the number of coalescent
histories for gene trees and species trees with a matching labeled topology
that belongs to a generic caterpillar-like family. By bringing a generating
function approach to the study of coalescent histories, we prove that for any
caterpillar-like family with seed tree $t$, the sequence $(h_n)_{n\geq 0}$
describing the number of matching coalescent histories of the $n$th tree of the
family grows asymptotically as a constant multiple of the Catalan numbers.
Thus, $h_n \sim \beta_t c_n$, where the asymptotic constant $\beta_t > 0$
depends on the shape of the seed tree $t$. The result extends a claim
demonstrated only for seed trees with at most 8 taxa to arbitrary seed trees,
expanding the set of cases for which detailed enumerative properties of
coalescent histories can be determined. We introduce a procedure that computes
from $t$ the constant $\beta_t$ as well as the algebraic expression for the
generating function of the sequence $(h_n)_{n\geq 0}$.


  We analyse interactions of Quantitative Trait Loci (QTL) in heat selected
yeast by comparing them to an unselected pool of random individuals. Here we
re-examine data on individual F12 progeny selected for heat tolerance, which
have been genotyped at 25 locations identified by sequencing a selected pool
[Parts, L., Cubillos, F. A., Warringer, J., Jain, K., Salinas, F., Bumpstead,
S. J., Molin, M., Zia, A., Simpson, J. T., Quail, M. A., Moses, A., Louis, E.
J., Durbin, R., and Liti, G. (2011). Genome research, 21(7), 1131-1138]. 960
individuals were genotyped at these locations and multi-locus genotype
frequencies were compared to 172 sequenced individuals from the original
unselected pool (a control group). Various non-random associations were found
across the genome, both within chromosomes and between chromosomes. Some of the
non-random associations are likely due to retention of linkage disequilibrium
in the F12 population, however many, including the inter-chromosomal
interactions, must be due to genetic interactions in heat tolerance. One region
of particular interest involves 3 linked loci on chromosome IV where the
central variant responsible for heat tolerance is antagonistic, coming from the
heat sensitive parent and the flanking ones are from the more heat tolerant
parent. The 3-locus haplotypes in the selected individuals represent a highly
biased sample of the population haplotypes with rare double recombinants in
high frequency. These were missed in the original analysis and would never be
seen without the multigenerational approach. We show that a statistical
analysis of entropy and information gain in genotypes of a selected population
can reveal further interactions than previously seen. Importantly this must be
done in comparison to the unselected population's genotypes to account for
inherent biases in the original population.


  Gaussian network model (GNM) is one of the most accurate and efficient
methods for biomolecular flexibility analysis. However, the systematic
generalization of the GNM has been elusive. We show that the GNM Kirchhoff
matrix can be built from the ideal low-pass filter, which is a special case of
a wide class of correlation functions underpinning the linear scaling
flexibility-rigidity index (FRI) method. Based on the mathematical structure of
correlation functions, we propose a unified framework to construct generalized
Kirchhoff matrices whose matrix inverse leads to correlation function based
GNMs, whereas, the direct inverse of the diagonal elements gives rise to FRI
method. We illustrate that correlation function based GNMs outperform the
original GNM in the B-factor prediction of a set of 364 proteins. We
demonstrate that for any given correlation function, FRI and GNM methods
provide essentially identical B-factor predictions when the scale value in the
correlation function is sufficiently large.


  Evolution is a dynamic process. The two classical forces of evolution are
mutation and selection. Assuming small mutation rates, evolution can be
predicted based solely on the fitness differences between phenotypes.
Predicting an evolutionary process under varying mutation rates as well as
varying fitness is still an open question. Experimental procedures, however, do
include these complexities along with fluctuating population sizes and
stochastic events such as extinctions. We investigate the mutational path
probabilities of systems having epistatic effects on both fitness and mutation
rates using a theoretical and computational framework. In contrast to previous
models, we do not limit ourselves to the typical strong selection, weak
mutation (SSWM)-regime or to fixed population sizes. Rather we allow epistatic
interactions to also affect mutation rates. This can lead to qualitatively
non-trivial dynamics. Pathways, that are negligible in the SSWM-regime, can
overcome fitness valleys and become accessible. This finding has the potential
to extend the traditional predictions based on the SSWM foundation and bring us
closer to what is observed in experimental systems.


  In a world in which many pressing global issues require large scale
cooperation, understanding the group size effect on cooperative behavior is a
topic of central importance. Yet, the nature of this effect remains largely
unknown, with lab experiments insisting that it is either positive or negative
or null, and field experiments suggesting that it is instead curvilinear. Here
we shed light on this apparent contradiction by considering a novel class of
public goods games inspired to the realistic scenario in which the natural
output limits of the public good imply that the benefit of cooperation
increases fast for early contributions and then decelerates. We report on a
large lab experiment providing evidence that, in this case, group size has a
curvilinear effect on cooperation, according to which intermediate-size groups
cooperate more than smaller groups and more than larger groups. In doing so,
our findings help fill the gap between lab experiments and field experiments
and suggest concrete ways to promote large scale cooperation among people.


  We develop a new thermodynamic approach to stochastic graph-rewriting. The
ingredients are a finite set of reversible graph-rewriting rules called
generating rules, a finite set of connected graphs P called energy patterns and
an energy cost function. The idea is that the generators define the qualitative
dynamics, by showing which transformations are possible, while the energy
patterns and cost function specify the long-term probability $\pi$ of any
reachable graph. Given the generators and energy patterns, we construct a
finite set of rules which (i) has the same qualitative transition system as the
generators; and (ii) when equipped with suitable rates, defines a
continuous-time Markov chain of which $\pi$ is the unique fixed point. The
construction relies on the use of site graphs and a technique of `growth
policy' for quantitative rule refinement which is of independent interest. This
division of labour between the qualitative and long-term quantitative aspects
of the dynamics leads to intuitive and concise descriptions for realistic
models (see the examples in S4 and S5). It also guarantees thermodynamical
consistency (AKA detailed balance), otherwise known to be undecidable, which is
important for some applications. Finally, it leads to parsimonious
parameterizations of models, again an important point in some applications.


  We propose two modeling approaches to describe the dynamics of ant battles,
starting from laboratory experiments on the behavior of two ant species, the
invasive Lasius neglectus and the authocthonus Lasius paralienus. This work is
mainly motivated by the need to have realistic models to predict the
interaction dynamics of invasive species. The two considered species exhibit
different fighting strategies. In order to describe the observed battle
dynamics, we start by building a chemical model considering the ants and the
fighting groups (for instance two ants of a species and one of the other one)
as a chemical species. From the chemical equations we deduce a system of
differential equations, whose parameters are estimated by minimizing the
difference between the experimental data and the model output. We model the
fluctuations observed in the experiments by means of a standard Gillespie
algorithm. In order to better reproduce the observed behavior, we adopt a
spatial agent-based model, in which ants not engaged in fighting groups move
randomly (diffusion) among compartments, and the Gillespie algorithm is used to
model the reactions inside a compartment.


  We model evolution according to an asymmetric game as occurring in multiple
finite populations, one for each role in the game, and study the effect of
subjecting individuals to stochastic strategy mutations. We show that, when
these mutations occur sufficiently infrequently, the dynamics over all
population states simplify to an ergodic Markov chain over just the pure
population states (where each population is monomorphic). This makes
calculation of the stationary distribution computationally feasible. The
transition probabilities of this embedded Markov chain involve fixation
probabilities of mutants in single populations. The asymmetry of the underlying
game leads to fixation probabilities that are derived from
frequency-independent selection, in contrast to the analogous single-population
symmetric-game case (Fudenberg and Imhof 2006). This frequency independence is
useful in that it allows us to employ results from the population genetics
literature to calculate the stationary distribution of the evolutionary
process, giving sharper, and sometimes even analytic, results. We demonstrate
the utility of this approach by applying it to a battle-of-the-sexes game, a
Crawford-Sobel signalling game, and the beer-quiche game of Cho and Kreps
(1987).


  Evolutionary game theory is used to model the evolution of competing
strategies in a population of players. Evolutionary stability of a strategy is
a dynamic equilibrium, in which any competing mutated strategy would be wiped
out from a population. If a strategy is weak evolutionarily stable, the
competing strategy may manage to survive within the network. Understanding the
network-related factors that affect the evolutionary stability of a strategy
would be critical in making accurate predictions about the behaviour of a
strategy in a real-world strategic decision making environment. In this work,
we evaluate the effect of network topology on the evolutionary stability of a
strategy. We focus on two well-known strategies known as the Zero-determinant
strategy and the Pavlov strategy. Zero-determinant strategies have been shown
to be evolutionarily unstable in a well-mixed population of players. We
identify that the Zero-determinant strategy may survive, and may even dominate
in a population of players connected through a non-homogeneous network. We
introduce the concept of `topological stability' to denote this phenomenon. We
argue that not only the network topology, but also the evolutionary process
applied and the initial distribution of strategies are critical in determining
the evolutionary stability of strategies. Further, we observe that topological
stability could affect other well-known strategies as well, such as the general
cooperator strategy and the cooperator strategy. Our observations suggest that
the variation of evolutionary stability due to topological stability of
strategies may be more prevalent in the social context of strategic evolution,
in comparison to the biological context.


  A similarity transformation is obtained between general population matrices
models of the Usher or Lefkovitch types and a simpler model, the pseudo-Leslie
model. The pseudo Leslie model is a matrix that can be decomposed in a row
matrix, which is not necessarily non-negative and a subdiagonal positive
matrix. This technique has computational advantages, since the solutions of the
iterative problem using Leslie matrices are readily obtained . In the case of
two age structured population models, one Lefkovitch and another Leslie, the
Kolmogorov-Sinai entropies are different, despite the same growth ratio of both
models. We prove that Markov matrices associated to similar population matrices
are similar.


  It has been hypothesized that the evolution of modern human cognition was
catalyzed by the development of jointly intentional modes of behaviour. From an
early age (1-2 years), human infants outperform apes at tasks that involve
collaborative activity. Specifically, human infants excel at joint action
motivated by reasoning of the form "we will do X" (shared intentions), as
opposed to reasoning of the form "I will do X [because he is doing X]"
(individual intentions). The mechanism behind the evolution of shared
intentionality is unknown. Here we formally model the evolution of jointly
intentional action and show under what conditions it is likely to have emerged
in humans. Modelling the interaction of hunter-gatherers as a coordination
game, we find that when the benefits from adopting new technologies or norms
are low but positive, the sharing of intentions does not evolve, despite being
a mutualistic behaviour that directly benefits all participants. When the
benefits from adopting new technologies or norms are high, such as may be the
case during a period of rapid environmental change, shared intentionality
evolves and rapidly becomes dominant in the population. Our results shed new
light on the evolution of collaborative behaviours.


  Molecular motors and cytoskeletal filaments work collectively most of the
time under opposing forces. This opposing force may be due to cargo carried by
motors or resistance coming from the cell membrane pressing against the
cytoskeletal filaments. Some recent studies have shown that the collective
maximum force (stall force) generated by multiple cytoskeletal filaments or
molecular motors may not always be just a simple sum of the stall forces of the
individual filaments or motors. To understand this excess or deficit in the
collective force, we study a broad class of models of both cytoskeletal
filaments and molecular motors. We argue that the stall force generated by a
group of filaments or motors is additive, that is, the stall force of $N$
number of filaments (motors) is $N$ times the stall force of one filament
(motor), when the system is in equilibrium at stall. Conversely, we show that
this additive property typically does not hold true when the system is not at
equilibrium at stall. We thus present a novel and unified understanding of the
existing models exhibiting such non-addivity, and generalise our arguments by
developing new models that demonstrate this phenomena. We also propose a
quantity similar to thermodynamic efficiency to easily predict this deviation
from stall-force additivity for filament and motor collectives.


  To understand the contact patterns of a population -- who is in contact with
whom, and when the contacts happen -- is crucial for modeling outbreaks of
infectious disease. Traditional theoretical epidemiology assumes that any
individual can meet any with equal probability. A more modern approach, network
epidemiology, assumes people are connected into a static network over which the
disease spreads. Newer yet, temporal network epidemiology, includes the time in
the contact representations. In this paper, we investigate the effect of these
successive inclusions of more information. Using empirical proximity data, we
study both outbreak sizes from unknown sources, and from known states of
ongoing outbreaks. In the first case, there are large differences going from a
fully mixed simulation to a network, and from a network to a temporal network.
In the second case, differences are smaller. We interpret these observations in
terms of the temporal network structure of the data sets. For example, a fast
overturn of nodes and links seem to make the temporal information more
important.


  Territoriality is a phenomenon exhibited throughout nature. On the individual
level, it is the processes by which organisms exclude others of the same
species from certain parts of space. On the population level, it is the
segregation of space into separate areas, each used by subsections of the
population. Proving mathematically that such individual-level processes can
cause observed population-level patterns to form is necessary for linking these
two levels of description in a non-speculative way. Previous mathematical
analysis has relied upon assuming animals are attracted to a central area. This
can either be a fixed geographical point, such as a den- or nest-site, or a
region where they have previously visited. However, recent simulation-based
studies suggest that this attractive potential is not necessary for territorial
pattern formation. Here, we construct a partial differential equation (PDE)
model of territorial interactions based on the individual-based model (IBM)
from those simulation studies. The resulting PDE does not rely on attraction to
spatial locations, but purely on conspecific avoidance, mediated via
scent-marking. We show analytically that steady-state patterns can form, as
long as (i) the scent does not decay faster than it takes the animal to
traverse the terrain, and (ii) the spatial scale over which animals detect
scent is incorporated into the PDE. As part of the analysis, we develop a
general method for taking the PDE limit of an IBM that avoids destroying any
intrinsic spatial scale in the underlying behavioral decisions.


  Recent studies reported complex post-transcriptional interplay among targets
of a common pool of microRNAs, a class of small non-coding downregulators of
gene expression. Behaving as microRNA-sponges, distinct RNA species may compete
for binding to microRNAs and coregulate each other in a dose-dependent manner.
Although previous studies in cell populations showed competition in vitro, the
detailed dynamical aspects of this process, most importantly in physiological
conditions, remains unclear. We address this point by monitoring protein
expression of two targets of a common miRNA with quantitative single-cell
measurements. In agreement with a detailed stochastic model of molecular
titration, we observed that: (i) crosstalk between targets is possible only in
particular stoichiometric conditions, (ii) a trade-off on the number of
microRNA regulatory elements may induce the coexistence of two distinct cell
populations, (iii) strong inter-targets correlations can be observed. This
phenomenology is compatible with a small amount of mRNA target molecules per
cell of the order of 10-100.


  We analyze a model statistical description of the polypeptide chain
helix-coil transition, where we take into account the specificity of its
primary sequence, as quantified by the phase space volume ratio of the number
of all accessible states to the number corresponding to a helical conformation.
The resulting transition phase diagram is then juxtaposed with the unusual
behavior of the secondary structures in Intrinsically Disordered Proteins
(IDPs) and a number of similarities are observed, even if the protein folding
is a more complex transition than the helix-coil transition. In fact, the
deficit in bulky and hydrophobic amino acids observed in IDPs, translated into
larger values of phase space volume, allows us to locate the region in
parameter space of the helix-coil transition that would correspond to the
secondary structure transformations that are intrinsic to conformational
transitions in IDPs and that is characterized by a modified phase diagram when
compared to globular proteins. Here, we argue how the nature of this modified
phase diagram, obtained from a model of the helix-coil transition in a solvent,
would illuminate the turned-out response of IDPs to the changes in the
environment conditions that follow straightforwardly from the re-entrant (cold
denaturation) branch in their folding phase diagram.


  Based on the reported data until 18 March 2015 and numerical fitting via a
simple formula of cumulative case number, we provide real-time estimation on
basic reproduction number, inflection point, peak time and final outbreak size
of ongoing Ebola outbreak in West Africa. From our simulation, we conclude that
the first wave has passed its inflection point and predict that a second
epidemic wave may appear in the near future.


  RNA molecules in living cells form what look like liquid droplets formed by
liquid/liquid phase separation. But unlike the molecules in conventional phase
separating mixtures, RNA molecules are transported by molecular motors that
consume energy and so are out of equilibrium. Motivated by this we consider
what sort of simple rules for the dynamics of model mRNA molecules lead to
liquid/liquid phase separation. We find that dynamics that slow as the local
density of molecules increases, drive the formation of liquids. We also look at
the analogous separation of the two blocks of a block copolymer, in which the
monomers of one block have dynamics that depend on the local density of
monomers of that block. We find that this block condenses and separates from
the monomers of the other block. This is a simple model of the
out-of-equilibrium domain formation found in the chromatin in the nucleus of
cells.


  Although metazoan development is conceived as resulting from gene regulatory
networks (GRNs) controlled by Hox genes, a better analogy is computer
architecture: i.e., a task accomplished in sequential steps linked to an
external referent that "counts" each step. A developmental "step" equals the
expression of genes in specific cells at specific times and telomeres represent
external "counters" wherein "counting" is a function of telomere shortening at
each cell division that permits the sequential expression of Hox genes and,
ultimately, complex form. Metazoan development thus best resembles a Turing
machine, which could be used to model the development of any metazoan.


  Subdiffusion has been proposed as an explanation of various kinetic phenomena
inside living cells. In order to fascilitate large-scale computational studies
of subdiffusive chemical processes, we extend a recently suggested mesoscopic
model of subdiffusion into an accurate and consistent reaction-subdiffusion
computational framework. Two different possible models of chemical reaction are
revealed and some basic dynamic properties are derived. In certain cases those
mesoscopic models have a direct interpretation at the macroscopic level as
fractional partial differential equations in a bounded time interval. Through
analysis and numerical experiments we estimate the macroscopic effects of
reactions under subdiffusive mixing. The models display properties observed
also in experiments: for a short time interval the behavior of the diffusion
and the reaction is ordinary, in an intermediate interval the behavior is
anomalous, and at long times the behavior is ordinary again.


  Most cancers in humans are large, measuring centimeters in diameter, composed
of many billions of cells. An equivalent mass of normal cells would be highly
heterogeneous as a result of the mutations that occur during each cell
division. What is remarkable about cancers is their homogeneity - virtually
every neoplastic cell within a large cancer contains the same core set of
genetic alterations, with heterogeneity confined to mutations that have emerged
after the last clonal expansions. How such clones expand within the
spatially-constrained three dimensional architecture of a tumor, and come to
dominate a large, pre-existing lesion, has never been explained. We here
describe a model for tumor evolution that shows how short-range migration and
cell turnover can account for rapid cell mixing inside the tumor. With it, we
show that even a small selective advantage of a single cell within a large
tumor allows the descendants of that cell to replace the precursor mass in a
clinically relevant time frame. We also demonstrate that the same mechanisms
can be responsible for the rapid onset of resistance to chemotherapy. Our model
not only provides novel insights into spatial and temporal aspects of tumor
growth but also suggests that targeting short range cellular migratory activity
could have dramatic effects on tumor growth rates.


  The main goal of Multiland is to generate neutral landscapes made of several
types of regions, with an exact control of the proportions occupied by each
type of region. An important feature of the software is that it allows a
control of the landscape fragmentation. It is intended to theoretical studies
on the effect of landscape structure in applied sciences. It has been developed
in the framework of the PEERLESS ANR project "Predictive Ecological Engineering
for Landscape Ecosystem Services and Sustainability".


  The Ebola virus is currently one of the most virulent pathogens for humans.
The latest major outbreak occurred in Guinea, Sierra Leonne and Liberia in
2014. With the aim of understanding the spread of infection in the affected
countries, it is crucial to modelize the virus and simulate it. In this paper,
we begin by studying a simple mathematical model that describes the 2014 Ebola
outbreak in Liberia. Then, we use numerical simulations and available data
provided by the World Health Organization to validate the obtained mathematical
model. Moreover, we develop a new mathematical model including vaccination of
individuals. We discuss different cases of vaccination in order to predict the
effect of vaccination on the infected individuals over time. Finally, we apply
optimal control to study the impact of vaccination on the spread of the Ebola
virus. The optimal control problem is solved numerically by using a direct
multi-shooting method.


  G protein-coupled receptors (GPCRs) constitute a large family of receptor
proteins that sense molecular signals on the exterior of a cell and activate
signal transduction pathways within the cell. Modeling how an agonist activates
such a receptor is fundamental for an understanding of a wide variety of
physiological processes and it is of tremendous value for pharmacology and drug
design. Inelastic electron tunneling spectroscopy (IETS) has been proposed as a
model for the mechanism by which olfactory GPCRs are activated by a bound
agonist. We apply this hypothesis to GPCRs within the mammalian nervous system
using quantum chemical modeling. We found that non-endogenous agonists of the
serotonin receptor share a particular IET spectral aspect both amongst each
other and with the serotonin molecule: a peak whose intensity scales with the
known agonist potencies. We propose an experiential validation of this model by
utilizing lysergic acid dimethylamide (DAM-57), an ergot derivative, and its
deuterated isotopologues; we also provide theoretical predictions for
comparison to experiment. If validated our theory may provide new avenues for
guided drug design and elevate methods of in silico potency/activity
prediction.


  The assumption that vector mortality remains constant with age is used widely
to assess malaria transmission risk and predict the public health consequences
of vector control strategies. However, laboratory studies commonly demonstrate
clear evidence of senescence, or a decrease in physiological function and
increase in vector mortality rate with age. We developed methods to integrate
available field data to understand mortality in wild Anopheles gambiae, the
most import vector of malaria in sub-Saharan Africa. We found evidence for an
increase in rates of mortality with age, a component of senescence. As
expected, we also found that overall mortality is far greater in wild cohorts
than commonly observed under protected laboratory conditions. The magnitude of
senescence increases with An. gambiae lifespan, implying that wild mosquitoes
die long before cohorts can exhibit strong senescence. We reviewed available
published mortality studies of Anopheles spp. to confirm this fundamental
prediction of aging in wild populations. Senescence becomes most apparent in
long-living mosquito cohorts, and cohorts with low extrinsic mortality, such as
those raised under protected laboratory conditions, suffer a relatively high
proportion of senescent deaths. Imprecision in estimates of vector mortality
and changes in mortality with age will severely bias models of vector borne
disease transmission risk, such as malaria, and the sensitivity of transmission
to bias increases as the extrinsic incubation period of the parasite decreases.
While we focus here on malaria, we caution that future models for
anti-vectorial interventions must therefore incorporate both realistic
mortality rates and age-dependent changes in vector mortality.


  During cell division, sister chromatids are segregated by the mitotic
spindle, a bipolar assembly of interdigitating antiparallel polar filaments
called microtubules. Establishing a stable overlap region is essential for
maintenance of bipolarity, but the underlying mechanisms are poorly understood.
Using a particle-based stochastic model, we find that the interplay of motors
and passive cross linkers can robustly generate partial overlaps between
antiparallel filaments. Our analysis shows that motors reduce the overlap in a
length-dependent manner, whereas passive cross linkers increase it
independently of the length. In addition to maintaining structural integrity,
passive cross linkers can thus also have a dynamic role for size regulation.


  Speciation is fundamental to the huge diversity of life on Earth. Evidence
suggests reproductive isolation arises most commonly in allopatry with a higher
speciation rate in small populations. Current theory does not address this
dependence in the important weak mutation regime. Here, we examine a
biophysical model of speciation based on the binding of a protein transcription
factor to a DNA binding site, and how their independent co-evolution, in a
stabilizing landscape, of two allopatric lineages leads to incompatibilities.
Our results give a new prediction for the monomorphic regime of evolution,
consistent with data, that smaller populations should develop incompatibilities
more quickly. This arises as: 1) smaller populations having a greater initial
drift load, as there are more sequences that bind poorly than well, so fewer
substitutions are needed to reach incompatible regions of phenotype space; 2)
slower divergence when the population size is larger than the inverse of
discrete differences in fitness. Further, we find longer sequences develop
incompatibilities more quickly at small population sizes, but more slowly at
large population sizes. The biophysical model thus represents a robust
mechanism of rapid reproductive isolation for small populations and large
sequences, that does not require peak-shifts or positive selection.


  Fitness landscapes are genotype to fitness mappings commonly used in
evolutionary biology and computer science which are closely related to spin
glass models. In this paper, we study the NK model for fitness landscapes where
the interaction scheme between genes can be explicitly defined. The focus is on
how this scheme influences the overall shape of the landscape. Our main tool
for the analysis are adaptive walks, an idealized dynamics by which the
population moves uphill in fitness and terminates at a local fitness maximum.
We use three different types of walks and investigate how their length (the
number of steps required to reach a local peak) and height (the fitness at the
endpoint of the walk) depend on the dimensionality and structure of the
landscape. We find that the distribution of local maxima over the landscape is
particularly sensitive to the choice of interaction pattern. Most quantities
that we measure are simply correlated to the rank of the scheme, which is equal
to the number of nonzero coefficients in the expansion of the fitness landscape
in terms of Walsh functions.


  Most of our understanding of bacterial chemotaxis comes from studies of
Escherichia coli. However, recent evidence suggests significant departures from
the E. coli paradigm in other bacterial species. This variation may stem from
different species inhabiting distinct environments and thus adapting to
specific environmental pressures. In particular, these complex and dynamic
environments may be poorly represented by standard experimental and theoretical
models. In this work, we study the performance of various chemotactic
strategies under a range of stochastic time- and space-varying attractant
distributions in silico. We describe a novel type of response in which the
bacterium tumbles more when attractant concentration is increasing, in contrast
to the response of E. coli, and demonstrate how this response explains the
behavior of aerobically-grown Rhodobacter sphaeroides. In this "speculator"
response, bacteria compare the current attractant concentration to the
long-term average. By tumbling persistently when the current concentration is
higher than the average, bacteria maintain their position in regions of high
attractant concentration. If the current concentration is lower than the
average, or is declining, bacteria swim away in search of more favorable
conditions. When the attractant distribution is spatially complex but
slowly-changing, this response is as effective as that of E. coli.


  The precise sequence of aminoacids plays a central role in the tertiary
structure of proteins and their functional properties. The Hydrophobic-Polar
lattice models have provided valuable insights regarding the energy landscape.
We demonstrate here the isomorphism between the protein sequences and
designable structures for two and three dimensional lattice proteins of very
long aminoacid chains using exact enumerations and intuitive considerations.We
emphasize that the topological arrangement of the aminoacid residues alone is
adequate to deduce the designable and non-designable sequences without explicit
recourse to energetics and degeneracies. The results indicate the computational
feasibility of realistic lattice models for proteins in two and three
dimensions and imply that the fundamental principle underlying the designing of
structures is the connectivity of the hydrophobic and polar residues.


  The thermodynamical stability of DNA minicircles is investigated by means of
path integral techniques. Hydrogen bonds between base pairs on complementary
strands can be broken by thermal fluctuations and temporary fluctuational
openings along the double helix are essential to biological functions such as
transcription and replication of the genetic information. Helix unwinding and
bubble formation patterns are computed in circular sequences with variable
radius in order to analyze the interplay between molecule size and appearance
of helical disruptions. The latter are found in minicircles with $< 100$ base
pairs and appear as a strategy to soften the stress due to the bending and
torsion of the helix.


  The estimation of covariance matrices of gene expressions has many
applications in cancer systems biology. Many gene expression studies, however,
are hampered by low sample size and it has therefore become popular to increase
sample size by collecting gene expression data across studies. Motivated by the
traditional meta-analysis using random effects models, we present a
hierarchical random covariance model and use it for the meta-analysis of gene
correlation networks across 11 large-scale gene expression studies of diffuse
large B-cell lymphoma (DLBCL). We suggest to use a maximum likelihood estimator
for the underlying common covariance matrix and introduce an EM algorithm for
estimation. By simulation experiments comparing the estimated covariance
matrices by cophenetic correlation and Kullback-Leibler divergence the
suggested estimator showed to perform better or not worse than a simple pooled
estimator. In a posthoc analysis of the estimated common covariance matrix for
the DLBCL data we were able to identify novel biologically meaningful gene
correlation networks with eigengenes of prognostic value. In conclusion, the
method seems to provide a generally applicable framework for meta-analysis,
when multiple features are measured and believed to share a common covariance
matrix obscured by study dependent noise.


  Phylogeographic methods aim to infer migration trends and the history of
sampled lineages from genetic data. Applications of phylogeography are broad,
and in the context of pathogens include the reconstruction of transmission
histories and the origin and emergence of outbreaks. Phylogeographic inference
based on bottom-up population genetics models is computationally expensive, and
as a result faster alternatives based on the evolution of discrete traits have
become popular. In this paper, we show that inference of migration rates and
root locations based on discrete trait models is extremely unreliable and
sensitive to biased sampling. To address this problem, we introduce BASTA
(BAyesian STructured coalescent Approximation), a new approach implemented in
BEAST2 that combines the accuracy of methods based on the structured coalescent
with the computational efficiency required to handle more than just few
populations. We illustrate the potentially severe implications of poor model
choice for phylogeographic analyses by investigating the zoonotic transmission
of Ebola virus. Whereas the structured coalescent analysis correctly infers
that successive human Ebola outbreaks have been seeded by a large unsampled
non-human reservoir population, the discrete trait analysis implausibly
concludes that undetected human-to-human transmission has allowed the virus to
persist over the past four decades. As genomics takes on an increasingly
prominent role informing the control and prevention of infectious diseases, it
will be vital that phylogeographic inference provides robust insights into
transmission history.


  We propose a general model to study the interplay between spatial dispersal
and environment spatiotemporal fluctuations in metapopulation dynamics. An
ecological landscape of favorable patches is generated like a L\'{e}vy dust,
which allows to build a range of patterns, from dispersed to clustered ones.
Locally, the dynamics is driven by a canonical model for the evolution of the
population density, consisting of a logistic expression plus multiplicative
noises. Spatial coupling is introduced by means of two spreading mechanisms:
diffusive dispersion and selective migration driven by patch suitability. We
focus on the long-time population size as a function of habitat configurations,
environment fluctuations and coupling schemes. We obtain the conditions, that
the spatial distribution of favorable patches and the coupling mechanisms must
fulfill, to grant population survival. The fundamental phenomenon that we
observe is the positive feedback between environment fluctuations and spatial
spread preventing extinction.


  Dynamics of p53 is known to play important roles in the regulation of cell
fate decisions in response to various stresses, and PDCD5 functions as a
co-activator of p53 to modulate the p53 dynamics. In the present paper, we
investigate how p53 dynamics are modulated by PDCD5 during the DNA damage
response using methods of bifurcation analysis and potential landscape. Our
results reveal that p53 activities can display rich dynamics under different
PDCD5 levels, including monostability, bistability with two stable steady
states, oscillations, and co-existence of a stable steady state and an
oscillatory state. Physical properties of the p53 oscillations are further
shown by the potential landscape, in which the potential force attracts the
system state to the limit cycle attractor, and the curl flux force drives the
coherent oscillation along the cyclic. We also investigate the effect of PDCD5
efficiency on inducing the p53 oscillations. We show that Hopf bifurcation is
induced by increasing the PDCD5 efficiency, and the system dynamics show clear
transition features in both barrier height and energy dissipation when the
efficiency is close to the bifurcation point. This study provides a global
picture of how PDCD5 regulates p53 dynamics via the interaction with the
p53-Mdm2 oscillator and can be helpful in understanding the complicate p53
dynamics in a more complete p53 pathway.


  The tumor suppressor p53 plays a central role in cell fate decisions after
DNA damage. Programmed Cell Death 5 (PDCD5) is known to interact with the p53
pathway to promote cell apoptosis. Recombinant human PDCD5 can significantly
sensitize different cancers to chemotherapies. In the present paper, we
construct a computational model that includes PDCD5 interactions in the p53
signaling network and study the effects of PDCD5 on p53-mediated cell fate
decisions during the DNA damage response. Our results revealed that PDCD5
functions as a co-activator of p53 that regulates p53-dependent cell fate
decisions via the mediation of p53 dynamics. The effects of PDCD5 are
dose-dependent such that p53 can display either sustained or pulsed dynamics at
different PDCD5 levels. Moreover, PDCD5 regulates caspase-3 activation via two
mechanisms during the two phases of sustained and pulsed p53 dynamics. This
study provides insights regarding how PDCD5 functions as a regulator of the p53
pathway and might be helpful for increasing our understanding of the molecular
mechanisms by which PDCD5 can be used to treat cancers.


  The objective of this research was to identify factors affecting tree damage
in the historical Minot flood of 2011. We hypothesized that tree height,
identity, origin, and maximum water height affect in the severity of damage
sustained by a tree in a flood event. All these factors were significant but
highly interactive. The results from this research can influence planting
practices in valleys and other flood prone areas to mitigate future damage.


  This paper investigates the dynamics of infectious diseases with a
non-exponentially distributed infectious period. This is achieved by
considering a multi-stage infection model on networks. Using pairwise
approximation with a standard closure, a number of important characteristics of
disease dynamics are derived analytically, including the final size of an
epidemic and a threshold for epidemic outbreaks, and it is shown how these
quantities depend on disease characteristics, as well as the number of disease
stages. Stochastic simulations of dynamics on networks are performed and
compared to the results of pairwise models for several realistic examples of
infectious diseases to illustrate the role played by the number of stages in
the disease dynamics. These results show that a higher number of disease stages
results in faster epidemic outbreaks with a higher peak prevalence and a larger
final size of the epidemic. The agreement between the pairwise and simulation
methods is excellent in the cases we consider.


  In this paper, we develop a multi-group epidemic framework via virtual
dispersal where the risk of infection is a function of the residence time and
local environmental risk. This novel approach eliminates the need to define and
measure contact rates that are used in the traditional multi-group epidemic
models with heterogeneous mixing. We apply this approach to a general $n$-patch
SIS model whose basic reproduction number $\mathcal R_0 $ is computed as a
function of a patch residence-times matrix $\mathbb P$. Our analysis implies
that the resulting $n$-patch SIS model has robust dynamics when patches are
strongly connected: there is a globally stable endemic equilibrium when
$\mathcal R_0>1 $ while the disease free equilibrium is globally stable when
$\mathcal R_0\leq1 $. Our further analysis indicates that the dispersal
behavior described by the residence-times matrix $\mathbb P$ has profound
effects on the disease dynamics at the single patch level with consequences
that proper dispersal behavior along with the local environmental risk can
either promote or eliminate the endemic in particular patches. Our work
highlights the impact of residence times matrix if the patches are not strongly
connected. Our framework can be generalized in other endemic and disease
outbreak models. As an illustration, we apply our framework to a two-patch SIR
single outbreak epidemic model where the process of disease invasion is
connected to the final epidemic size relationship. We also explore the impact
of disease prevalence driven decision using a phenomenological modeling
approach in order to contrast the role of constant versus state dependent
$\mathbb P$ on disease dynamics.


  HIV-1 can disseminate between susceptible cells by two mechanisms: cell-free
infection following fluid-phase diffusion of virions and by highly-efficient
direct cell-to-cell transmission at immune cell contacts. The contribution of
this hybrid spreading mechanism, which is also a characteristic of some
important computer worm outbreaks, to HIV-1 progression in vivo remains
unknown. Here we present a new mathematical model that explicitly incorporates
the ability of HIV-1 to use hybrid spreading mechanisms and evaluate the
consequences for HIV-1 pathogenenesis. The model captures the major phases of
the HIV-1 infection course of a cohort of treatment naive patients and also
accurately predicts the results of the Short Pulse Anti-Retroviral Therapy at
Seroconversion (SPARTAC) trial. Using this model we find that hybrid spreading
is critical to seed and establish infection, and that cell-to-cell spread and
increased CD4+ T cell activation are important for HIV-1 progression. Notably,
the model predicts that cell-to-cell spread becomes increasingly effective as
infection progresses and thus may present a considerable treatment barrier.
Deriving predictions of various treatments' influence on HIV-1 progression
highlights the importance of earlier intervention and suggests that treatments
effectively targeting cell-to-cell HIV-1 spread can delay progression to AIDS.
This study suggests that hybrid spreading is a fundamental feature of HIV
infection, and provides the mathematical framework incorporating this feature
with which to evaluate future therapeutic strategies.


  Background Elimination of malaria can only be achieved through removal of all
vectors or complete depletion of the infectious reservoir in humans.
Mechanistic models can be built to synthesize diverse observations from the
field collected under a variety of conditions and subsequently used to query
the infectious reservoir in great detail. Methods The EMOD model of malaria
transmission was calibrated to prevalence, incidence, asexual parasite density,
gametocyte density, infection duration, and infectiousness data from 9 study
sites. The infectious reservoir was characterized by diagnostic detection limit
and age group over a range of transmission intensities with and without case
management and vector control. Mass screen-and-treat drug campaigns were tested
for likelihood of achieving elimination. Results The composition of the
infectious reservoir by diagnostic threshold is similar over a range of
transmission intensities, and higher intensity settings are biased toward
infections in children. Recent ramp-ups in case management and use of
insecticide-treated bednets reduce the infectious reservoir and shift the
composition toward submicroscopic infections. Mass campaigns with antimalarial
drugs are highly effective at interrupting transmission if deployed shortly
after ITN campaigns. Conclusions Low density infections comprise a substantial
portion of the infectious reservoir. Proper timing of vector control, seasonal
variation in transmission intensity, and mass drug campaigns allows lingering
population immunity to help drive a region toward elimination.


  Antimalarial drugs are a powerful tool for malaria control and elimination.
Artemisinin-based combination therapies (ACTs) can reduce transmission when
widely distributed in a campaign setting. Modelling mass antimalarial campaigns
can elucidate how to most effectively deploy drug-based interventions and
quantitatively compare the effects of cure, prophylaxis, and
transmission-blocking in suppressing parasite prevalence. A previously
established agent-based model that includes innate and adaptive immunity was
used to simulate malaria infections and transmission. Pharmacokinetics of
artemether, lumefantrine, dihydroartemisinin, piperaquine, and primaquine were
modelled with a double-exponential distribution-elimination model including
weight-dependent parameters and age-dependent dosing. Drug killing of asexual
parasites and gametocytes was calibrated to clinical data. Mass distribution of
ACTs and primaquine was simulated with seasonal mosquito dynamics at a range of
transmission intensities. A single mass campaign with antimalarial drugs is
insufficient to permanently reduce malaria prevalence when transmission is
high. Current diagnostics are insufficiently sensitive to accurately identify
asymptomatic infections, and mass-screen-and-treat campaigns are much less
efficacious than mass drug administrations. Improving campaign coverage leads
to decreased prevalence one month after the end of the campaign, while
increasing compliance lengthens the duration of protection against reinfection.
Use of a long-lasting prophylactic as part of a mass drug administration
regimen confers the most benefit under conditions of high transmission and
moderately high coverage. Addition of primaquine can reduce prevalence but
exerts its largest effect when coupled with a long-lasting prophylactic.


  For many membrane proteins the determination of their topology remains a
challenge for methods like X-ray crystallography and nuclear magnetic resonance
(NMR) spectroscopy. Electron paramagnetic resonance (EPR) spectroscopy has
evolved as an alternative technique to study structure and dynamics of membrane
proteins. The present study demonstrates the feasibility of membrane protein
topology determination using limited EPR distance and accessibility
measurements. The BCL::MP-Fold (BioChemical Library membrane protein fold)
algorithm assembles secondary structure elements (SSEs) in the membrane using a
Monte Carlo Metropolis (MCM) approach. Sampled models are evaluated using
knowledge-based potential functions and agreement with the EPR data and a
knowledge-based energy function. Twenty-nine membrane proteins of up to 696
residues are used to test the algorithm. The RMSD100 value of the most accurate
model is better than 8 {\AA} for twenty-seven, better than 6 {\AA} for
twenty-two and better than 4 {\AA} for fifteen out of twenty-nine proteins,
demonstrating the algorithms ability to sample the native topology. The average
enrichment could be improved from 1.3 to 2.5, showing the improved
discrimination power by using EPR data.


  Persistent homology has been devised as a promising tool for the topological
simplification of complex data. However, it is computationally intractable for
large data sets. In this work, we introduce multiresolution persistent homology
for tackling large data sets. Our basic idea is to match the resolution with
the scale of interest so as to create a topological microscopy for the
underlying data. We utilize flexibility-rigidity index (FRI) to access the
topological connectivity of the data set and define a rigidity density for the
filtration analysis. By appropriately tuning the resolution, we are able to
focus the topological lens on a desirable scale. The proposed multiresolution
topological analysis is validated by a hexagonal fractal image which has three
distinct scales. We further demonstrate the proposed method for extracting
topological fingerprints from DNA and RNA molecules. In particular, the
topological persistence of a virus capsid with 240 protein monomers is
successfully analyzed which would otherwise be inaccessible to the normal point
cloud method and unreliable by using coarse-grained multiscale persistent
homology. The proposed method has also been successfully applied to the protein
domain classification, which is the first time that persistent homology is used
for practical protein domain analysis, to our knowledge. The proposed
multiresolution topological method has potential applications in arbitrary data
sets, such as social networks, biological networks and graphs.


  One of the key characteristics of cancer cells is an increased phenotypic
plasticity, driven by underlying genetic and epigenetic perturbations. However,
at a systems-level it is unclear how these perturbations give rise to the
observed increased plasticity. Elucidating such systems-level principles is key
for an improved understanding of cancer. Recently, it has been shown that
signaling entropy, an overall measure of signaling pathway promiscuity, and
computable from integrating a sample's gene expression profile with a protein
interaction network, correlates with phenotypic plasticity and is increased in
cancer compared to normal tissue. Here we develop a computational framework for
studying the effects of network perturbations on signaling entropy. We
demonstrate that the increased signaling entropy of cancer is driven by two
factors: (i) the scale-free (or near scale-free) topology of the interaction
network, and (ii) a subtle positive correlation between differential gene
expression and node connectivity. Indeed, we show that if protein interaction
networks were random graphs, described by Poisson degree distributions, that
cancer would generally not exhibit an increased signaling entropy. In summary,
this work exposes a deep connection between cancer, signaling entropy and
interaction network topology.


  Cell signaling, gene expression, and metabolism are affected by cell-cell
heterogeneity and random changes in the environment. The effects of such
fluctuations on cell signaling and gene expression have recently been studied
intensively using single-cell experiments. In metabolism heterogeneity may be
particularly important because it may affect synchronisation of metabolic
oscillations, an important example of cell-cell communication. This
synchronisation is notoriously difficult to describe theoretically as the
example of glycolytic oscillations shows: neither is the mechanism of
glycolytic synchronisation understood nor the role of cell-cell heterogeneity.
To pin down the mechanism and to assess its robustness and universality we have
experimentally investigated the entrainment of glycolytic oscillations in
individual yeast cells by periodic external perturbations. We find that
oscillatory cells synchronise through phase shifts and that the mechanism is
insensitive to cell heterogeneity (robustness) and similar for different types
of external perturbations (universality).


  Intracellular calcium is regulated in part by the release of Ca$^{2+}$ ions
from the endoplasmic reticulum via inositol-4,5-triphosphate receptor (IP$_3$R)
channels (among other possibilities such as RyR and L-type calcium channels).
The resulting dynamics are highly diverse, lead to local calcium "puffs" as
well as global waves propagating through cells, as observed in {\it Xenopus}
oocytes, neurons, and other cell types. Local fluctuations in the number of
calcium ions play a crucial role in the onset of these features. Previous
modeling studies of calcium puff dynamics stemming from IP$_3$R channels have
predominantly focused on stochastic channel models coupled to deterministic
diffusion of ions, thereby neglecting local fluctuations of the ion number.
Tracking of individual ions is computationally difficult due to the scale
separation in the Ca$^{2+}$ concentration when channels are in the open or
closed states. In this paper, a spatial multiscale model for investigating of
the dynamics of puffs is presented. It couples Brownian motion (diffusion) of
ions with a stochastic channel gating model. The model is used to analyze
calcium puff statistics. Concentration time traces as well as channel state
information are studied. We identify the regime in which puffs can be found and
develop a mean-field theory to extract the boundary of this regime. Puffs are
only possible when the time scale of channel inhibition is sufficiently large.
Implications for the understanding of puff generation and termination are
discussed.


  Keratin are among the most abundant proteins in epithelial cells. Functions
of the keratin network in cells are shaped by their dynamical organization.
Using a collection of experimentally-driven mathematical models, different
hypotheses for the turnover and transport of the keratin material in epithelial
cells are tested. The interplay between turnover and transport and their
effects on the keratin organization in cells are hence investigated by
combining mathematical modeling and experimental data. Amongst the collection
of mathematical models considered, a best model strongly supported by
experimental data is identified. Fundamental to this approach is the fact that
optimal parameter values associated with the best fit for each model are
established. The best candidate among the best fits is characterized by the
disassembly of the assembled keratin material in the perinuclear region and an
active transport of the assembled keratin. Our study shows that an active
transport of the assembled keratin is required to explain the experimentally
observed keratin organization.


  Network-based methods are playing an increasingly important role in drug
design. Our main question in this paper was whether the efficiency of drug
target proteins to spread perturbations in the human interactome is larger if
the binding drugs have side effects, as compared to those which have no
reported side effects. Our results showed that in general, drug targets were
better spreaders of perturbations than non-target proteins, and in particular,
targets of drugs with side effects were also better spreaders of perturbations
than targets of drugs having no reported side effects in human protein-protein
interaction networks. Colorectal cancer-related proteins were good spreaders
and had a high centrality, while type 2 diabetes-related proteins showed an
average spreading efficiency and had an average centrality in the human
interactome. Moreover, the interactome-distance between drug targets and
disease-related proteins was higher in diabetes than in colorectal cancer. Our
results may help a better understanding of the network position and dynamics of
drug targets and disease-related proteins, and may contribute to develop
additional, network-based tests to increase the potential safety of drug
candidates.


  The canonical view of the interactions between viruses and their microbial
hosts presumes that changes in host and virus fate require the initiation of
infection of a host by a virus. That is, first virus particles diffuse randomly
outside of host cells, then the virus genome enters the target host cell, and
only then do intracellular dynamics and regulation of virus and host cell fate
unfold. Intracellular dynamics may lead to the death of the host cell and
release of viruses, to the elimination of the virus genome through cellular
defense mechanisms, or the integration of the virus genome with the host as a
chromosomal or extra-chromosomal element. Here we revisit this canonical view,
inspired by recent experimental findings of Bautista and colleagues (mBio,
2015) in which the majority of target host cells can be induced into a dormant
state when exposed to either active or de-activated viruses, even when viruses
are present at low relative titer. We propose that both the qualitative
phenomena and the quantitative time-scales of dormancy induction can be
reconciled given the hypothesis that cellular physiology can be altered by
contact on the surface of host cells rather than strictly by infection. We
develop a biophysical model of contact-mediated dynamics involving virus
particles and target cells. We show how in this model virus particles can
catalyze - extracellularly - cellular transformations amongst many cells, even
if they ultimately infect only one (or none). We discuss implications of the
present biophysical model relevant to the study of virus-microbe interactions
more generally.


  The phase coherence relaxation times as long as $T_2\sim830-1030\pm20$ ns
were measured for the [2Fe-2S] cluster in the intrinsic protein environment.
This relaxation corresponds to a relatively long lasting coherence of the
low-spin $S=1/2$ state. For this biological cluster, the phase coherence
relaxation time was significantly affected by the nuclear hyperfine
interactions of $^{14}$N with $I=1$. After labeling the cluster environments
uniformly with the $^{15}$N isotope, $T_2$ exceeds $\sim1.1-1.4\pm0.1~\mu$s at
the canonical orientations. This is already an order of magnitude longer than
the duration of a single-spin qubit manipulation $\sim10-100$ ns. The transient
nutation experiment corresponds to the coherent manipulation of the electron
spin.


  Statistical phylogenetic inference methods use tree rearrangement operations
to perform either hill-climbing local search or Markov chain Monte Carlo across
tree topologies. The canonical class of such moves are the
subtree-prune-regraft (SPR) moves that remove a subtree and reattach it
somewhere else via the cut edge of the subtree. Phylogenetic trees and such
moves naturally form the vertices and edges of a graph, such that tree search
algorithms perform a (potentially stochastic) traversal of this SPR graph.
Despite the centrality of such graphs in phylogenetic inference, rather little
is known about their large-scale properties. In this paper we learn about the
rooted-tree version of the graph, known as the rSPR graph, by calculating the
Ricci-Ollivier curvature for pairs of vertices in the rSPR graph with respect
to two simple random walks on the rSPR graph. By proving theorems and direct
calculation with novel algorithms, we find a remarkable diversity of different
curvatures on the rSPR graph for pairs of vertices separated by the same
distance. We confirm using simulation that degree and curvature have the
expected impact on mean access time distributions, demonstrating relevance of
these curvature results to stochastic tree search. This indicates significant
structure of the rSPR graph beyond that which was previously understood in
terms of pairwise distances and vertex degrees; a greater understanding of
curvature could ultimately lead to improved strategies for tree search.


  In this commentary I utilize the general methods of the mathematical theory
of heterogeneous populations in order to point out an omission in the analysis
of the mathematical model in [Dwyer et al (2000), Am Nat, 156(2):105--120],
which led to the conclusion in [Elderd et al (2008) Am Nat, 172(6):829--842]
that the original model must be replaced with an alternative one because of the
new data. I show that more thorough understanding of the underlying model
allows twitching the model parameters to account for the observed data.


  In a population of size N, adaptive evolution is 2N times faster under
Mendelian inheritance than the rate implied by Victorian theories of heredity
and evolution.


  A model of an Ebola epidemic is developed with infected individuals
structured according to disease age. The transmission of the infection is
tracked by disease age through an initial incubation (exposed) phase, followed
by an infectious phase with variable transmission infectiousness. The removal
of infected individuals is dependent on disease age, with three types of
removal rates: (1) removal due to hospitalization (isolation), (2) removal due
to mortality separate from hospitalization, and (3) removal due to recovery
separate from hospitalization. The model is applied to the Ebola epidemics in
Sierra Leone and Guinea. Model simulations indicate that successive stages of
increased and earlier hospitalization of cases have resulted in mitigation of
the epidemics.


  Charge hydration asymmetry (CHA)--a characteristic dependence of hydration
free energy on the sign of the solute charge--quantifies the asymmetric
response of water to electric field at microscopic level. Accurate estimates of
CHA are critical for understanding hydration effects ubiquitous in chemistry
and biology. However, measuring hydration energies of charged species is
fraught with significant difficulties, which lead to unacceptably large (up to
300%) variation in the available estimates of the CHA effect. We circumvent
these difficulties by developing a framework which allows us to extract and
accurately estimate the intrinsic propensity of water to exhibit CHA from
accurate experimental hydration free energies of neutral polar molecules.
Specifically, from a set of 504 small molecules we identify two pairs that are
analogous, with respect to CHA, to the K+/F- pair--a classical probe for the
effect. We use these "CHA-conjugate" molecule pairs to quantify the intrinsic
charge-asymmetric response of water to the microscopic charge perturbations:
the asymmetry of the response is strong, ~50% of the average hydration free
energy of these molecules. The ability of widely used classical water models to
predict hydration energies of small molecules correlates with their ability to
predict CHA.


  We present p53-MDM2-Glucose model to study spatio-temporal properties of the
system induced by glucose. The variation in glucose concentration level
triggers the system at different states, namely, oscillation death
(stabilized), sustain and damped oscillations which correspond to various
cellular states. The transition of these states induced by glucose is phase
transition like behaviour. We also found that the intrinsic noise in stochastic
system helps the system to stabilize more effectively. Further, the amplitude
of $p53$ dynamics with the variation of glucose concentration level follows
power law behaviour, $A_s(k)\sim k^\gamma$, where, $\gamma$ is a constant.


  We study the efficiency of synchronization in ensembles of identical coupled
stochastic oscillator systems. By deriving a chemical Langevin equation, we
measure the rate at which the systems synchronize. The rate at which the
difference in the Hilbert phases of the systems evolve provides a suitable
order parameter, and a 2--dimensional recurrence plot further facilitates the
analysis of stochastic synchrony. We find that a global mean--field coupling
effects the most rapid approach to global synchrony, and that when the number
of "information carrying" molecular species increases, the rate of synchrony
increases. The Langevin analysis is complemented by numerical simulations.


  Aligning reads to a reference sequence is a fundamental step in numerous
bioinformatics pipelines. As a consequence, the sensitivity and precision of
the mapping tool, applied with certain parameters to certain data, can
critically affect the accuracy of produced results (e.g., in variant calling
applications). Therefore, there has been an increasing demand of methods for
comparing mappers and for measuring effects of their parameters.
  Read simulators combined with alignment evaluation tools provide the most
straightforward way to evaluate and compare mappers. Simulation of reads is
accompanied by information about their positions in the source genome. This
information is then used to evaluate alignments produced by the mapper.
Finally, reports containing statistics of successful read alignments are
created.
  In default of standards for encoding read origins, every evaluation tool has
to be made explicitly compatible with the simulator used to generate reads. In
order to solve this obstacle, we have created a generic format RNF (Read Naming
Format) for assigning read names with encoded information about original
positions.
  Futhermore, we have developed an associated software package RNF containing
two principal components. MIShmash applies one of popular read simulating tools
(among DwgSim, Art, Mason, CuReSim etc.) and transforms the generated reads
into RNF format. LAVEnder evaluates then a given read mapper using simulated
reads in RNF format. A special attention is payed to mapping qualities that
serve for parametrization of ROC curves, and to evaluation of the effect of
read sample contamination.


  Translation is a central biological process by which proteins are synthesized
from genetic information contained within mRNAs. Here we study the kinetics of
translation at molecular level through a stochastic simulation model. The model
explicitly include RNA sequences, ribosome dynamics, tRNA pool and biochemical
reactions in the translation elongation. The results show that the translation
efficiency is mainly limited by the available ribosome number, translation
initiation and the translation elongation time. The elongation time is
log-normal distribution with mean and variance determined by both the codon
saturation and the process of aa-tRNA selection at each codon binding.
Moreover, our simulations show that the translation accuracy exponentially
decreases with the sequence length. These results suggest that aa-tRNA
competition is crucial for both translation elongation, translation efficiency
and the accuracy, which in turn determined the effective protein production
rate of correct proteins. Our results improve the dynamical equation of protein
production with a delay differential equation which is dependent on sequence
informations through both the effective production rate and the distribution of
elongation time.


  Traditional disease surveillance systems suffer from several disadvantages,
including reporting lags and antiquated technology, that have caused a movement
towards internet-based disease surveillance systems. Internet systems are
particularly attractive for disease outbreaks because they can provide data in
near real-time and can be verified by individuals around the globe. However,
most existing systems have focused on disease monitoring and do not provide a
data repository for policy makers or researchers. In order to fill this gap, we
analyzed Wikipedia article content.
  We demonstrate how a named-entity recognizer can be trained to tag case
counts, death counts, and hospitalization counts in the article narrative that
achieves an F1 score of 0.753. We also show, using the 2014 West African Ebola
virus disease epidemic article as a case study, that there are detailed time
series data that are consistently updated that closely align with ground truth
data.
  We argue that Wikipedia can be used to create the first community-driven
open-source emerging disease detection, monitoring, and repository system.


  If deleterious mutations near a fitness maximum in a spatially distributed
population are sufficiently frequent or detrimental, the population can undergo
a fitness collapse, similarly to the Muller's ratchet effect in well-mixed
populations. Recent studies of one-dimensional habitats (e.g., the frontier of
a two-dimensional range expansion) have shown that the onset of the fitness
collapse is described by a directed percolation phase transition with its
associated critical exponents. We consider population fitness collapse in
three-dimensional range expansions with both inflating and fixed-size frontiers
(applicable to, e.g., expanding and treadmilling spherical tumors,
respectively). We find that the onset of fitness collapse in these two cases
obeys different scaling laws, and that competition between species at the
frontier leads to a deviation from directed percolation scaling. As in
two-dimensional range expansions, inflating frontiers modify the critical
behavior by causally disconnecting well-separated portions of the population.


  We introduce an extended version of oxDNA, a coarse-grained model of DNA
designed to capture the thermodynamic, structural and mechanical properties of
single- and double-stranded DNA. By including explicit major and minor grooves,
and by slightly modifying the coaxial stacking and backbone-backbone
interactions, we improve the ability of the model to treat large
(kilobase-pair) structures such as DNA origami which are sensitive to these
geometric features. Further, we extend the model, which was previously
parameterised to just one salt concentration ([Na$^+$]=0.5M), so that it can be
used for a range of salt concentrations including those corresponding to
physiological conditions. Finally, we use new experimental data to parameterise
the oxDNA potential so that consecutive adenine bases stack with a different
strength to consecutive thymine bases, a feature which allows a more accurate
treatment of systems where the flexibility of single-stranded regions is
important. We illustrate the new possibilities opened up by the updated model,
oxDNA2, by presenting results from simulations of the structure of large DNA
objects and by using the model to investigate some salt-dependent properties of
DNA.


  Despite recent technological advances in genomic sciences, our understanding
of cancer progression and its driving genetic alterations remains incomplete.
Here, we introduce TiMEx, a generative probabilistic model for detecting
patterns of various degrees of mutual exclusivity across genetic alterations,
which can indicate pathways involved in cancer progression. TiMEx explicitly
accounts for the temporal interplay between the waiting times to alterations
and the observation time. In simulation studies, we show that our model
outperforms previous methods for detecting mutual exclusivity. On large-scale
biological datasets, TiMEx identifies gene groups with strong functional
biological relevance, while also proposing many new candidates for biological
validation. TiMEx possesses several advantages over previous methods, including
a novel generative probabilistic model of tumorigenesis, direct estimation of
the probability of mutual exclusivity interaction, computational efficiency, as
well as high sensitivity in detecting gene groups involving low-frequency
alterations. R code is available at www.cbg.bsse.ethz.ch/software/TiMEx.


  Predators often consume multiple prey and by mutually subsidizing a shared
predator, the prey may reciprocally harm each other. When predation levels are
high, this apparent competition can culminate in a prey species being
displaced. Coupling quantitative genetics and Lotka-Volterra models, we study
how predator evolution alters this and other ecological outcomes. These models
account for a trade-off between the predator's attack rates on two prey
species. We provide a mathematical characterization of a strong form of
persistence--permanence--for which there is a global attractor bounded away
from extinction. When the evolutionary dynamics occur at a sufficiently slower
time scale than the ecological dynamics, we also characterize attractors and
their basins' of attraction using singular perturbation theory and a graphical
approach to the eco-evolutionary dynamics. Our results show that
eco-evolutionary feedbacks can mediate permanence at intermediate trade-offs in
the attack rates. However, at strong trade-offs, permanence is lost. Despite
this loss of permanence, there can be attractors supporting coexistence. These
attractors, however, may coincide with attractors at which the predator is
excluded. Our results highlight that evo-evolutionary feedbacks can alter
community structure by mediating coexistence or leading to trait-dependent
alternative stable states.


  Next Generation Sequencing (NGS) technologies generate large amounts of short
read data for many different organisms. The fact that NGS reads are generally
short makes it challenging to assemble the reads and reconstruct the original
genome sequence. For clustering genomes using such NGS data, word-count based
alignment-free sequence comparison is a promising approach, but for this
approach, the underlying expected word counts are essential.
  A plausible model for this underlying distribution of word counts is given
through modelling the DNA sequence as a Markov chain (MC). For single long
sequences, efficient statistics are available to estimate the order of MCs and
the transition probability matrix for the sequences. As NGS data do not provide
a single long sequence, inference methods on Markovian properties of sequences
based on single long sequences cannot be directly used for NGS short read data.
  Here we derive a normal approximation for such word counts. We also show that
the traditional Chi-square statistic has an approximate gamma distribution,
using the Lander-Waterman model for physical mapping. We propose several
methods to estimate the order of the MC based on NGS reads and evaluate them
using simulations. We illustrate the applications of our results by clustering
genomic sequences of several vertebrate and tree species based on NGS reads
using alignment-free sequence dissimilarity measures. We find that the
estimated order of the MC has a considerable effect on the clustering results,
and that the clustering results that use a MC of the estimated order give a
plausible clustering of the species.


  Many biological questions, including the estimation of deep evolutionary
histories and the detection of remote homology between protein sequences, rely
upon multiple sequence alignments (MSAs) and phylogenetic trees of large
datasets. However, accurate large-scale multiple sequence alignment is very
difficult, especially when the dataset contains fragmentary sequences. We
present UPP, an MSA method that uses a new machine learning technique - the
Ensemble of Hidden Markov Models - that we propose here. UPP produces highly
accurate alignments for both nucleotide and amino acid sequences, even on
ultra-large datasets or datasets containing fragmentary sequences. UPP is
available at https://github.com/smirarab/sepp.


  Across many fields, a problem of interest is to predict the transition rates
between nodes of a network, given limited stationary state and dynamical
information. We give a solution using the principle of Maximum Caliber. We find
the transition rate matrix by maximizing the path entropy of a random walker on
the network constrained to reproducing a stationary distribution and a few
dynamical averages. A main finding here is that when constrained only by the
mean jump rate, the rate matrix is given by a square-root dependence of the
rate, $\omega_{ab} \propto \sqrt{p_b/p_a}$, on $p_a$ and $p_b$, the stationary
state populations at nodes a and b. We give two examples of our approach.
First, we show that this method correctly predicts the correlated rates in a
biochemical network of two genes, where we know the exact results from prior
simulation. Second, we show that it correctly predicts rates of peptide
conformational transitions, when compared to molecular dynamics simulations.
This method can be used to infer large numbers of rates on known networks where
smaller numbers of steady-state node populations are known.


  Metagenomics research has accelerated the studies of microbial organisms,
providing insights into the composition and potential functionality of various
microbial communities. Metatranscriptomics (studies of the transcripts from a
mixture of microbial species) and other meta-omics approaches hold even greater
promise for providing additional insights into functional and regulatory
characteristics of the microbial communities. Current metatranscriptomics
projects are often carried out without matched metagenomic datasets (of the
same microbial communities). For the projects that produce both
metatranscriptomic and metagenomic datasets, their analyses are often not
integrated. Metagenome assemblies are far from perfect, partially explaining
why metagenome assemblies are not used for the analysis of metatranscriptomic
datasets. Here we report a reads mapping algorithm for mapping of short reads
onto a de Bruijn graph of assemblies. A hash table of junction k-mers (k-mers
spanning branching structures in the de Bruijn graph) is used to facilitate
fast mapping of reads to the graph. We developed an application of this mapping
algorithm: a reference based approach to metatranscriptome assembly using
graphs of metagenome assembly as the reference. Our results show that this new
approach (called TAG) helps to assemble substantially more transcripts that
otherwise would have been missed or truncated because of the fragmented nature
of the reference metagenome. TAG was implemented in C++ and has been tested
extensively on the linux platform. It is available for download as open source
at http://omics.informatics.indiana.edu/TAG.


  Campylobacter jejuni and Campylobacter coli are the biggest causes of
bacterial gastroenteritis in the developed world, with human infections
typically arising from zoonotic transmission associated with infected meat,
especially poultry. Because this organism is not thought to survive well
outside of the gut, host associated populations are genetically isolated to
varying degrees. Therefore the likely origin of most Campylobacter strains can
be determined by host-associated variation in the genome. This is instructive
for characterizing the source of human infection at the population level.
However, some common strains appear to have broad host ranges, hindering source
attribution. Whole genome sequencing has the potential to reveal fine-scale
genetic structure associated with host specificity within each of these
strains.
  We found that rates of zoonotic transmission among animal host species in
ST-21, ST-45 and ST-828 clonal complexes were so high that the signal of host
association is all but obliterated. We attributed 89% of clinical cases to a
chicken source, 10% to cattle and 1% to pig. Our results reveal that common
strains of C. jejuni and C. coli infectious to humans are adapted to a
generalist lifestyle, permitting rapid transmission between different hosts.
Furthermore, they show that the weak signal of host association within these
complexes presents a challenge for pinpointing the source of clinical
infections, underlining the view that whole genome sequencing, powerful though
it is, cannot substitute for intensive sampling of suspected transmission
reservoirs.


  Continuous-time Markov processes over finite state-spaces are widely used to
model dynamical processes in many fields of natural and social science. Here,
we introduce an maximum likelihood estimator for constructing such models from
data observed at a finite time interval. This estimator is dramatically more
efficient than prior approaches, enables the calculation of deterministic
confidence intervals in all model parameters, and can easily enforce important
physical constraints on the models such as detailed balance. We demonstrate and
discuss the advantages of these models over existing discrete-time Markov
models for the analysis of molecular dynamics simulations.


  The detection of genomic structural variation (SV) has advanced tremendously
in recent years due to progress in high-throughput sequencing technologies.
Novel sequence insertions, insertions without similarity to a human reference
genome, have received less attention than other types of SVs due to the
computational challenges in their detection from short read sequencing data,
which inherently involves de novo assembly. De novo assembly is not only
computationally challenging, but also requires high-quality data. While the
reads from a single individual may not always meet this requirement, using
reads from multiple individuals can increase power to detect novel insertions.
We have developed the program PopIns, which can discover and characterize
non-reference insertions of 100 bp or longer on a population scale. In this
paper, we describe the approach we implemented in PopIns. It takes as input a
reads-to-reference alignment, assembles unaligned reads using a standard
assembly tool, merges the contigs of different individuals into high-confidence
sequences, anchors the merged sequences into the reference genome, and finally
genotypes all individuals for the discovered insertions. Our tests on simulated
data indicate that the merging step greatly improves the quality and
reliability of predicted insertions and that PopIns shows significantly better
recall and precision than the recent tool MindTheGap. Preliminary results on a
data set of 305 Icelanders demonstrate the practicality of the new approach.
The source code of PopIns is available from http://github.com/bkehr/popins.


  It is well known that in contrast to the Prisoner's Dilemma, the snowdrift
game can lead to a stable coexistence of cooperators and cheaters. Recent
theoretical evidence on the snowdrift game suggests that gradual evolution for
individuals choosing to contribute in continuous degrees can result in the
social diversification to a 100% contribution and 0% contribution through
so-called evolutionary branching. Until now, however, game-theoretical studies
have shed little light on the evolutionary dynamics and consequences of the
loss of diversity in strategy. Here we analyze continuous snowdrift games with
quadratic payoff functions in dimorphic populations. Subsequently, conditions
are clarified under which gradual evolution can lead a population consisting of
those with 100% contribution and those with 0% contribution to merge into one
species with an intermediate contribution level. The key finding is that the
continuous snowdrift game is more likely to lead to assimilation of different
cooperation levels rather than maintenance of diversity. Importantly, this
implies that allowing the gradual evolution of cooperative behavior can
facilitate social inequity aversion in joint ventures that otherwise could
cause conflicts that are based on commonly accepted notions of fairness.


  We study the properties of the potential overlap between two networks $A,B$
sharing the same set of $N$ nodes (a two-layer network) whose respective degree
distributions $p_A(k), p_B(k)$ are given. Defining the overlap coefficient
$\alpha$ as the Jaccard index, we derive upper bounds for the minimum and
maximum overlap coefficient in terms of $p_A(k)$, $p_B(k)$ and $N$. We also
present an algorithm based on cross-rewiring of links to obtain a two-layer
network with any prescribed $\alpha$ inside the permitted range. Finally, to
illustrate the importance of the overlap for the dynamics of interacting
contagious processes, we derive a mean-field model for the spread of an SIS
epidemic with awareness against infection over a two-layer network, containing
$\alpha$ as a parameter. A simple analytical relationship between $\alpha$ and
the basic reproduction number follows. Stochastic simulations are presented to
assess the accuracy of the upper bounds of $\alpha$ and the predictions of the
mean-field epidemic model.


  Protein distributions measured under a broad set of conditions in bacteria
and yeast were shown to exhibit a common skewed shape, with variances depending
quadratically on means. For bacteria these properties were reproduced by
temporal measurements of protein content, showing accumulation and division
across generations. Here we present a stochastic growth-and-division model with
feedback which captures these observed properties. The limiting copy number
distribution is calculated exactly, and a single parameter is found to
determine the distribution shape and the variance-to-mean relation. Estimating
this parameter from bacterial temporal data reproduces the measured
distribution shape with high accuracy, and leads to predictions for future
experiments.