  We establish a link between Wakeley et al's (2012) cyclical pedigree model
from population genetics and a randomized directed configuration model (DCM)
considered by Cooper and Frieze (2004). We then exploit this link in
combination with asymptotic results for the in-degree distribution of the
corresponding DCM to compute the asymptotic size of the largest strongly
connected component $S^N$ (where $N$ is the population size) of the DCM resp.
the pedigree. The size of the giant component can be characterized explicitly
(amounting to approximately $80 \%$ of the total populations size) and thus
contributes to a reduced `pedigree effective population size'. In addition, the
second largest strongly connected component is only of size $O(\log N)$.
Moreover, we describe the size and structure of the `domain of attraction' of
$S^N$. In particular, we show that with high probability for any individual the
shortest ancestral line reaches $S^N$ after $O(\log \log N)$ generations, while
almost all other ancestral lines take at most $O(\log N)$ generations.


  RNA sequencing (RNA-seq) enables characterization and quantification of
individual transcriptomes as well as detection of patterns of allelic
expression and alternative splicing. Current RNA-seq protocols depend on
high-throughput short-read sequencing of cDNA. However, as ongoing advances are
rapidly yielding increasing read lengths, a technical hurdle remains in
identifying the degree to which differences in read length influence various
transcriptome analyses. In this study, we generated two paired-end RNA-seq
datasets of differing read lengths (2x75 bp and 2x262 bp) for lymphoblastoid
cell line GM12878 and compared the effect of read length on transcriptome
analyses, including read-mapping performance, gene and transcript
quantification, and detection of allele-specific expression (ASE) and
allele-specific alternative splicing (ASAS) patterns. Our results indicate
that, while the current long-read protocol is considerably more expensive than
short-read sequencing, there are important benefits that can only be achieved
with longer read length, including lower mapping bias and reduced ambiguity in
assigning reads to genomic elements, such as mRNA transcript. We show that
these benefits ultimately lead to improved detection of cis-acting regulatory
and splicing variation effects within individuals.


  Some recent results have introduced external-memory algorithms to compute
self-indexes of a set of strings, mainly via computing the Burrows-Wheeler
Transform (BWT) of the input strings. The motivations for those results stem
from Bioinformatics, where a large number of short strings (called reads) are
routinely produced and analyzed. In that field, a fundamental problem is to
assemble a genome from a large set of much shorter samples extracted from the
unknown genome. The approaches that are currently used to tackle this problem
are memory-intensive. This fact does not bode well with the ongoing increase in
the availability of genomic data. A data structure that is used in genome
assembly is the string graph, where vertices correspond to samples and arcs
represent two overlapping samples. In this paper we address an open problem: to
design an external-memory algorithm to compute the string graph.


  The binding of proteins onto DNA contributes to the shaping and packaging of
genome as well as to the expression of specific genetic messages. With a view
to understanding the interplay between the presence of proteins and the
deformation of DNA involved in such processes, we developed a new method to
minimize the elastic energy of DNA fragments at the mesoscale level. Our method
makes it possible to obtain the optimal pathways of protein-decorated DNA
molecules for which the terminal base pairs are spatially constrained. We focus
in this work on the deformations induced by selected architectural proteins on
circular DNA. We report the energy landscapes of DNA minicircles subjected to
different levels of torsional stress and containing one or two proteins as
functions of the chain length and spacing between the proteins. Our results
reveal cooperation between the elasticity of the double helix and the
structural distortions of DNA induced by bound proteins. We find that the
imposed mechanical stress influences the placement of proteins on DNA and that,
the proteins, in turn, modulate the mechanical stress and thereby broadcast
their presence along DNA.


  Since its introduction, landscape genomics has developed quickly with the
increasing availability of both molecular and topo-climatic data. The current
challenges of the field mainly involve processing large numbers of models and
disentangling selection from demography. Several methods address the latter,
either by estimating a neutral model from population structure or by inferring
simultaneously environmental and demographic effects. Here we present
Sam$\beta$ada, an integrated approach to study signatures of local adaptation,
providing rapid processing of whole genome data and enabling assessment of
spatial association using molecular markers. Specifically, candidate loci to
adaptation are identified by automatically assessing genome-environment
associations. In complement, measuring the Local Indicators of Spatial
Association (LISA) for these candidate loci allows to detect whether similar
genotypes tend to gather in space, which constitutes a useful indication of the
possible kinship relationship between individuals. In this paper, we also
analyze SNP data from Ugandan cattle to detect signatures of local adaptation
with Sam$\beta$ada, BayEnv, LFMM and an outlier method (FDIST approach in
Arlequin) and compare their results. Sam$\beta$ada is an open source software
for Windows, Linux and MacOS X available at \url{http://lasig.epfl.ch/sambada}


  Transport networks play a key role across four realms of eukaryotic life:
slime molds, fungi, plants, and animals. In addition to the developmental
algorithms that build them, many also employ adaptive strategies to respond to
stimuli, damage, and other environmental changes. We model these adapting
network architectures using a generic dynamical system on weighted graphs and
find in simulation that these networks ultimately develop a hierarchical
organization of the final weighted architecture accompanied by the formation of
a system-spanning backbone. In addition, we find that the long term
equilibration dynamics exhibit glassy behavior characterized by long periods of
slow changes punctuated by bursts of reorganization events.


  The growing usage of nanoparticles of zinc sulfide as quantum dots and
biosensors calls for a theoretical assessment of interactions of ZnS with
biomolecules. We employ the molecular-dynamics-based umbrella sampling method
to determine potentials of mean force for 20 single amino acids near the ZnS
(110) surface in aqueous solutions. We find that five amino acids do not bind
at all and the binding energy of the remaining amino acids does not exceed 4.3
kJ/mol. Such energies are comparable to those found for ZnO (and to hydrogen
bonds in proteins) but the nature of the specificity is different. Cysteine can
bind with ZnS in a covalent way, e.g. by forming the disulfide bond with S in
the solid. If this effect is included within a model incorporating the Morse
potential, then the potential well becomes much deeper - the binding energy is
close to 98 kJ/mol. We then consider tryptophan cage, a protein of 20 residues,
and characterize its events of adsorption to ZnS. We demonstrate the relevance
of interactions between the amino acids in the selection of optimal adsorbed
conformations and recognize the key role of cysteine in generation of lasting
adsorption. We show that ZnS is more hydrophobic than ZnO and that the density
profile of water is quite different than that forming near ZnO - it has only a
minor articulation into layers. Furthermore, the first layer of water is
disordered and mobile.


  We determine potentials of the mean force for interactions of amino acids
with four common surfaces of ZnO in aqueous solutions. The method involves
all-atom molecular dynamics simulations combined with the umbrella sampling
technique. The profiled nature of the density of water with the strongly
adsorbed first layer affects the approach of amino acids to the surface and
generates either repulsion or weak binding. The largest binding energy is found
for tyrosine interacting with the surface in which the Zn ions are at the top.
It is equal to 7 kJ/mol which is comparable to that of the hydrogen bonds in a
protein. This makes the adsorption of amino acids onto the ZnO surface much
weaker than onto the well studied surface of gold. Under vacuum, binding
energies are more than 40 times stronger (for one of the surfaces). The precise
manner in which water molecules interact with a given surface influences the
binding energies in a way that depends on the surface. Among the four
considered surfaces the one with Zn at the top is recognized as binding almost
all amino acids with an average binding energy of 2.60 kJ/mol. Another (O at
the top) is non-binding for most amino acids. For binding situations the
average energy is 0.66 kJ/mol. The remaining two surfaces bind nearly as many
amino acids as they do not and the average binding energies are 1.46 and 1.22
kJ/mol. For all of the surfaces the binding energies vary between amino acids
significantly: the dispersion in the range of 68-154% of the mean. A small
protein is shown to adsorb onto ZnO only intermittently and with only a small
deformation. Various adsorption events lead to different patterns in mobilities
of amino acids within the protein.


  We calculate potentials of the mean force for twenty amino acids in the
vicinity of the (111) surface of gold, for several dipeptides, and for some
analogs of the side chains, using molecular dynamics simulations and the
umbrella sampling method. We compare results obtained within three different
force fields: one hydrophobic (for a contaminated surface) and two hydrophilic.
All of these fields lead to good binding with very different specificities and
different patterns in the density and polarization of water. The covalent bond
with the sulfur atom on cysteine is modeled by the Morse potential. We
demonstrate that binding energies of dipeptides are different than the combined
binding energies of their amino-acidic components. For the hydrophobic gold,
adsorption events of a small protein are driven by attraction to the strongest
binding amino acids. This is not so in the hydrophilic cases - a result of
smaller specificities combined with the difficulty for proteins, but not for
single amino acids, to penetrate the first layer of water. The properties of
water near the surface sensitively depend on the force field.


  Population dynamics and in particular microbial population dynamics, though
they are complex but also intrinsically discrete and random, are conventionally
represented as deterministic differential equations systems. We propose to
revisit this approach by complementing these classic formalisms by stochastic
formalisms and to explain the links between these representations in terms of
mathematical analysis but also in terms of modeling and numerical simulations.
We illustrate this approach on the model of chemostat.


  As an equilibrium refinement of the Nash equilibrium, evolutionarily stable
strategy (ESS) is a key concept in evolutionary game theory and has attracted
growing interest. An ESS can be either a pure strategy or a mixed strategy.
Even though the randomness is allowed in mixed strategy, the selection
probability of pure strategy in a mixed strategy may fluctuate due to the
impact of many factors. The fluctuation can lead to more uncertainty. In this
paper, such uncertainty involved in mixed strategy has been further taken into
consideration: a belief strategy is proposed in terms of Dempster-Shafer
evidence theory. Furthermore, based on the proposed belief strategy, a
belief-based ESS has been developed. The belief strategy and belief-based ESS
can reduce to the mixed strategy and mixed ESS, which provide more realistic
and powerful tools to describe interactions among agents.


  Exchange type chromosome aberrations (ETCAs) are rearrangements of the genome
that occur when chromosomes break and the resulting fragments rejoin with other
fragments from other chromosomes. ETCAs are commonly observed in cancer cells
and in cells exposed to radiation. The frequency of these chromosome
rearrangements is correlated with their spatial proximity, therefore it can be
used to infer the three dimensional organization of the genome. Extracting
statistical significance of spatial proximity from cancer and radiation data
has remained somewhat elusive because of the sparsity of the data. We here
propose a new approach to study the three dimensional organization of the
genome using algebraic statistics. We test our method on a published data set
of irradiated human blood lymphocyte cells. We provide a rigorous method for
testing the overall organization of the genome, and in agreement with previous
results we find a random relative positioning of chromosomes with the exception
of the chromosome pairs \{1,22\} and \{13,14\} that have a significantly larger
number of ETCAs than the rest of the chromosome pairs suggesting their spatial
proximity. We conclude that algebraic methods can successfully be used to
analyze genetic data and have potential applications to larger and more complex
data sets.


  We study the total least squares (TLS) problem that generalizes least squares
regression by allowing measurement errors in both dependent and independent
variables. TLS is widely used in applied fields including computer vision,
system identification and econometrics. The special case when all dependent and
independent variables have the same level of uncorrelated Gaussian noise, known
as ordinary TLS, can be solved by singular value decomposition (SVD). However,
SVD cannot solve many important practical TLS problems with realistic noise
structure, such as having varying measurement noise, known structure on the
errors, or large outliers requiring robust error-norms. To solve such problems,
we develop convex relaxation approaches for a general class of structured TLS
(STLS). We show both theoretically and experimentally, that while the plain
nuclear norm relaxation incurs large approximation errors for STLS, the
re-weighted nuclear norm approach is very effective, and achieves better
accuracy on challenging STLS problems than popular non-convex solvers. We
describe a fast solution based on augmented Lagrangian formulation, and apply
our approach to an important class of biological problems that use population
average measurements to infer cell-type and physiological-state specific
expression levels that are very hard to measure directly.


  Networks are a unifying framework for modeling complex systems and network
inference problems are frequently encountered in many fields. Here, I develop
and apply a generative approach to network inference (RCweb) for the case when
the network is sparse and the latent (not observed) variables affect the
observed ones. From all possible factor analysis (FA) decompositions explaining
the variance in the data, RCweb selects the FA decomposition that is consistent
with a sparse underlying network. The sparsity constraint is imposed by a novel
method that significantly outperforms (in terms of accuracy, robustness to
noise, complexity scaling, and computational efficiency) Bayesian methods and
MLE methods using l1 norm relaxation such as K-SVD and l1--based sparse
principle component analysis (PCA). Results from simulated models demonstrate
that RCweb recovers exactly the model structures for sparsity as low (as
non-sparse) as 50% and with ratio of unobserved to observed variables as high
as 2. RCweb is robust to noise, with gradual decrease in the parameter ranges
as the noise level increases.


  By means of computer simulations of a coarse-grained DNA model we show that
the DNA hairpin zippering dynamics is anomalous, i.e. the characteristic time T
scales non-linearly with N, the hairpin length: T ~ N^a with a>1. This is in
sharp contrast with the prediction of the zipper model for which T ~ N. We show
that the anomalous dynamics originates from an increase in the friction during
zippering due to the tension built in the closing strands. From a simple
polymer model we get a = 1+ nu = 1.59 with nu the Flory exponent, a result
which is in agreement with the simulations. We discuss transition path times
data where such effects should be detected.


  Applying a method to reconstruct a phylogenetic tree from random data
provides a way to detect whether that method has an inherent bias towards
certain tree `shapes'. For maximum parsimony, applied to a sequence of random
2-state data, each possible binary phylogenetic tree has exactly the same
distribution for its parsimony score. Despite this pleasing and slightly
surprising symmetry, some binary phylogenetic trees are more likely than others
to be a most parsimonious (MP) tree for a sequence of $k$ such characters, as
we show. For $k=2$, and unrooted binary trees on six taxa, any tree with a
caterpillar shape has a higher chance of being an MP tree than any tree with a
symmetric shape. On the other hand, if we take any two binary trees, on any
number of taxa, we prove that this bias between the two trees vanishes as the
number of characters grows. However, again there is a twist: MP trees on six
taxa are more likely to have certain shapes than a uniform distribution on
binary phylogenetic trees predicts, and this difference does not appear to
dissipate as $k$ grows.


  Ecological networks exhibit non-random structural patterns, such as
modularity and nestedness, which indicate ecosystem stability, species
diversity, and connectance. Such structure-stability relationships are well
known. However, another important perspective is less well understood: the
relationship between the environment and structure. Inspired by theoretical
studies that suggest that network structure can change due to environmental
variability, we collected data on a number of empirical food webs and
mutualistic networks and evaluated the effect of climatic seasonality on
ecological network structure. As expected, we found that climatic seasonality
affects ecological network structure. In particular, an increase in modularity
due to climatic seasonality was observed in food webs; however, it is debatable
whether this occurs in mutualistic networks. Interestingly, the type of
climatic seasonality that affects network structure differs with ecosystem
type. Rainfall and temperature seasonality influence freshwater food webs and
mutualistic networks, respectively; food webs are smaller, and more modular,
with increasing rainfall seasonality. Mutualistic networks exhibit a higher
diversity (particularly of animals) with increasing temperature seasonality.
These results confirm the theoretical prediction that stability increases with
greater perturbation. Although these results are still debatable because of
several limitations in the data analysis, they may enhance our understanding of
environment-structure relationships.


  Understanding the regulation and structure of ribosomes is essential to
understanding protein synthesis and its deregulation in disease. While
ribosomes are believed to have a fixed stoichiometry among their core ribosomal
proteins (RPs), some experiments suggest a more variable composition. Testing
such variability requires direct and precise quantification of RPs. We used
mass-spectrometry to directly quantify RPs across monosomes and polysomes of
mouse embryonic stem cells (ESC) and budding yeast. Our data show that the
stoichiometry among core RPs in wild-type yeast cells and ESC depends both on
the growth conditions and on the number of ribosomes bound per mRNA.
Furthermore, we find that the fitness of cells with a deleted RP-gene is
inversely proportional to the enrichment of the corresponding RP in polysomes.
Together, our findings support the existence of ribosomes with distinct protein
composition and physiological function.


  Gaining a genomic perspective on phylogeny requires the collection of data
from many putatively independent loci collected across the genome. Among
insects, an increasingly common approach to collecting this class of data
involves transcriptome sequencing, because few insects have high-quality genome
sequences available; assembling new genomes remains a limiting factor; the
transcribed portion of the genome is a reasonable, reduced subset of the genome
to target; and the data collected from transcribed portions of the genome are
similar in composition to the types of data with which biologists have
traditionally worked (e.g., exons). However, molecular techniques requiring RNA
as a template are limited to using very high quality source materials, which
are often unavailable from a large proportion of biologically important insect
samples. Recent research suggests that DNA-based target enrichment of conserved
genomic elements offers another path to collecting phylogenomic data across
insect taxa, provided that conserved elements are present in and can be
collected from insect genomes. Here, we identify a large set (n$=$1510) of
ultraconserved elements (UCE) shared among the insect order Hymenoptera. We use
in silico analyses to show that these loci accurately reconstruct relationships
among genome-enabled Hymenoptera, and we design a set of baits for enriching
these loci that researchers can use with DNA templates extracted from a variety
of sources. We use our UCE bait set to enrich an average of 721 UCE loci from
30 hymenopteran taxa, and we use these UCE loci to reconstruct phylogenetic
relationships spanning very old ($\geq$220 MYA) to very young ($\leq$1 MYA)
divergences among hymenopteran lineages. In contrast to a recent study
addressing hymenopteran phylogeny using transcriptome data, we found ants to be
sister to all remaining aculeate lineages with complete support.


  Within the context of evolution, an altruistic act that benefits the
receiving individual at the expense of the acting individual is a puzzling
phenomenon. An extreme form of altruism can be found in colicinogenic E. coli.
These suicidal altruists explode, releasing colicins that kill unrelated
individuals, which are not colicin resistant. By committing suicide, the
altruist makes it more likely that its kin will have less competition. The
benefits of this strategy rely on the number of competitors and kin nearby. If
the organism explodes at an inopportune time, the suicidal act may not harm any
competitors. Communication could enable organisms to act altruistically when
environmental conditions suggest that that strategy would be most beneficial.
Quorum sensing is a form of communication in which bacteria produce a protein
and gauge the amount of that protein around them. Quorum sensing is one means
by which bacteria sense the biotic factors around them and determine when to
produce products, such as antibiotics, that influence competition. Suicidal
altruists could use quorum sensing to determine when exploding is most
beneficial, but it is challenging to study the selective forces at work in
microbes. To address these challenges, we use digital evolution (a form of
experimental evolution that uses self-replicating computer programs as
organisms) to investigate the effects of enabling altruistic organisms to
communicate via quorum sensing. We found that quorum-sensing altruists killed a
greater number of competitors per explosion, winning competitions against
non-communicative altruists. These findings indicate that quorum sensing could
increase the beneficial effect of altruism and the suite of conditions under
which it will evolve.


  Summary: We present a new method to incrementally construct the FM-index for
both short and long sequence reads, up to the size of a genome. It is the first
algorithm that can build the index while implicitly sorting the sequences in
the reverse (complement) lexicographical order without a separate sorting step.
The implementation is among the fastest for indexing short reads and the only
one that practically works for reads of averaged kilobases in length.
  Availability and implementation: https://github.com/lh3/ropebwt2
  Contact: hengli@broadinstitute.org


  The reconstruction of transmission trees for epidemics from genetic data has
been the subject of some recent interest. It has been demonstrated that the
transmission tree structure can be investigated by augmenting internal nodes of
a phylogenetic tree constructed using pathogen sequences from the epidemic with
information about the host that held the corresponding lineage. In this paper,
we note that this augmentation is equivalent to a correspondence between
transmission trees and partitions of the phylogenetic tree into connected
subtrees each containing one tip, and provide a framework for Markov Chain
Monte Carlo inference of phylogenies that are partitioned in this way, giving a
new method to co-estimate both trees. The procedure is integrated in the
existing phylogenetic inference package BEAST.


  Given a sample of genome sequences from an asexual population, can one
predict its evolutionary future? Here we demonstrate that the branching
patterns of reconstructed genealogical trees contains information about the
relative fitness of the sampled sequences and that this information can be used
to predict successful strains. Our approach is based on the assumption that
evolution proceeds by accumulation of small effect mutations, does not require
species specific input and can be applied to any asexual population under
persistent selection pressure. We demonstrate its performance using historical
data on seasonal influenza A/H3N2 virus. We predict the progenitor lineage of
the upcoming influenza season with near optimal performance in 30% of cases and
make informative predictions in 16 out of 19 years. Beyond providing a tool for
prediction, our ability to make informative predictions implies persistent
fitness variation among circulating influenza A/H3N2 viruses.


  Tuberculosis (TB) and human immunodeficiency virus (HIV) can be considered a
deadly human syndemic. In this article, we formulate a model for TB and HIV
transmission dynamics. The model considers both TB and acquired immune
deficiency syndrome (AIDS) treatment for individuals with only one of the
infectious diseases or both. The basic reproduction number and equilibrium
points are determined and stability is analyzed. Through simulations, we show
that TB treatment for individuals with only TB infection reduces the number of
individuals that become co-infected with TB and HIV/AIDS, and reduces the
diseases (TB and AIDS) induced deaths. Analogously, the treatment of
individuals with only AIDS also reduces the number of co-infected individuals.
Further, TB-treatment for co-infected individuals in the active and latent
stage of TB disease, implies a decrease of the number of individuals that
passes from HIV-positive to AIDS.


  The contact map of a protein fold is a graph that represents the patterns of
contacts in the fold. It is known that the contact map can be decomposed into
stacks and queues. RNA secondary structures are special stacks in which the
degree of each vertex is at most one and each arc has length at least two.
Waterman and Smith derived a formula for the number of RNA secondary structures
of length $n$ with exactly $k$ arcs. H\"{o}ner zu Siederdissen et al. developed
a folding algorithm for extended RNA secondary structures in which each vertex
has maximum degree two. An equation for the generating function of extended RNA
secondary structures was obtained by M\"{u}ller and Nebel by using a
context-free grammar approach, which leads to an asymptotic formula. In this
paper, we consider $m$-regular linear stacks, where each arc has length at
least $m$ and the degree of each vertex is bounded by two. Extended RNA
secondary structures are exactly $2$-regular linear stacks. For any $m\geq 2$,
we obtain an equation for the generating function of the $m$-regular linear
stacks. For given $m$, we can deduce a recurrence relation and an asymptotic
formula for the number of $m$-regular linear stacks on $n$ vertices. To
establish the equation, we use the reduction operation of Chen, Deng and Du to
transform an $m$-regular linear stack to an $m$-reduced zigzag (or alternating)
stack. Then we find an equation for $m$-reduced zigzag stacks leading to an
equation for $m$-regular linear stacks.


  Long range dependence or long memory is a feature of many processes in the
natural world, which provides important insights on the underlying mechanisms
that generate the observed data. The usual tools available to characterize the
phenomenon are mostly based on second order correlations. However, the long
memory effects may not be evident at the level of second order correlations and
may require a deeper analysis of the nature of the stochastic processes.
  After a short review of the notions and tools used to characterize long range
dependence, we analyse data related to the abundance of exploited fish
populations which provides an example of higher order long range dependence.


  Pathogenic bacteria present a large disease burden on human health. Control
of these pathogens is hampered by rampant lateral gene transfer, whereby
pathogenic strains may acquire genes conferring resistance to common
antibiotics. Here we introduce tools from topological data analysis to
characterize the frequency and scale of lateral gene transfer in bacteria,
focusing on a set of pathogens of significant public health relevance. As a
case study, we examine the spread of antibiotic resistance in Staphylococcus
aureus. Finally, we consider the possible role of the human microbiome as a
reservoir for antibiotic resistance genes.


  We introduce a model for microtubule mechanics containing lateral bonds
between dimers in neighboring protofilaments, bending rigidity of dimers, and
repulsive interactions between protofilaments modeling steric constraints to
investigate the influence of mechanical forces on hydrolysis and catastrophes.
We use the allosteric dimer model, where tubulin dimers are characterized by an
equilibrium bending angle, which changes from $0^\circ$ to $22^\circ$ by
hydrolysis of a dimer. This also affects the lateral interaction and bending
energies and, thus, the mechanical equilibrium state of the microtubule. As
hydrolysis gives rise to conformational changes in dimers, mechanical forces
also influence the hydrolysis rates by mechanical energy changes modulating the
hydrolysis rate. The interaction via the microtubule mechanics then gives rise
to correlation effects in the hydrolysis dynamics, which have not been taken
into account before. Assuming a dominant influence of mechanical energies on
hydrolysis rates, we investigate the most probable hydrolysis pathways both for
vectorial and random hydrolysis. Investigating the stability with respect to
lateral bond rupture, we identify initiation configurations for catastrophes
along the hydrolysis pathways and values for a lateral bond rupture force. If
we allow for rupturing of lateral bonds between dimers in neighboring
protofilaments above this threshold force, our model exhibits avalanche-like
catastrophe events.


  The Min proteins from Escherichia coli can self-organize into traveling waves
on supported lipid bilayers. In Proc. Natl. Acad. Sci. USA 109, 15283 (2012) we
showed that these waves are guided along the boundaries of membrane patches. We
introduced an effective two-dimensional model reproducing the observed
patterns. In arXiv:1403.5934v1, Jacob Halatek and Erwin Frey contest the
ability of our effective two-dimensional model to describe the dynamics of Min
proteins on patterned supported lipid bilayers. We thank Halatek and Frey for
their interest in our work and for again highlighting the importance of
dimensionality and geometry for pattern formation by the Min proteins. Here we
reply in detail to the objections by Halatek and Frey and show that (1) our
effective two-dimensional model reproduces the observed patterns on isolated
patches and that (2) a three-dimensional version of our model produces similar
patterns on square patches.


  Bacteria produce a range of proteolytic enzymes, for which a number human
equivalent or structurally similar examples exist and the primary focus of this
study was to analyse the published literature to find proteolytic enzymes,
specifically endoproteses and to examine the similarity in the substrates that
they act on so as to predict a suitable structural motif which can be used as
the basis for preparation of useful prodrug carriers against diseases caused by
specific bacteria like Salmonella. Also, the similarities between the bacterial
proteases and the action of human matrix metalloproteinases (MMPs), together
with the MMP-like activity of bacterial endoproteases to activate human MMPs,
were also analysed. This information was used to try to identify substrates on
which the MMPs and bacterial proteases act, to aid the design of oligopeptide
prodrug carriers to treat cancer and its metastatic spread. MMPs are greatly
involved in cancer growth and progression, a few MMPs and certain proteases
share a similar type of activity in degrading the extra cellular matrix (ECM)
and substrates including gelatin. Our primary targets of study were to identify
the proteases and MMPs that facilitate the migration of bacteria and growth of
tumour cells respectively. The study was thus a two-way approach to study the
substrate specificity of both bacterial proteases and MMPs, thereby to help in
characterisation of their substrates. Various bioinformatics tools were used in
the characterisation of the proteases and substrates as well as in the
identification of possible binding sites and conserved regions in a range of
candidate proteins. Central to this project was the salmonella derived PgtE
surface protease that has been shown recently to act upon the pro-forms of
human MMP-9.


  Many real world networks contain a statistically surprising number of certain
subgraphs, called network motifs. In the prevalent approach to motif analysis,
network motifs are detected by comparing subgraph frequencies in the original
network with a statistical null model. In this paper we propose an alternative
approach to motif analysis where network motifs are defined to be connectivity
patterns that occur in a subgraph cover that represents the network using
minimal total information. A subgraph cover is defined to be a set of subgraphs
such that every edge of the graph is contained in at least one of the subgraphs
in the cover. Some recently introduced random graph models that can incorporate
significant densities of motifs have natural formulations in terms of subgraph
covers and the presented approach can be used to match networks with such
models. To prove the practical value of our approach we also present a
heuristic for the resulting NP-hard optimization problem and give results for
several real world networks.


  This paper continues the analysis, started in [2, 3], of a class of
degenerate elliptic operators defined on manifolds with corners, which arise in
Population Biology. Using techniques pioneered by J. Moser, and extended and
refined by L. Saloff-Coste, Grigoryan, and Sturm, we show that weak solutions
to the parabolic problem defined by a sub-class of these operators, which
consists of those that can be defined by Dirichlet forms and have non-vanishing
transverse vector field, satisfy a Harnack inequality. This allows us to
conclude that the solutions to these equations belong, for positive times, to
the natural anisotropic Holder spaces, and also leads to upper and, in some
cases, lower bounds for the heat kernels of these operators. These results
imply that these operators have a compact resolvent when acting on $C^0$ or
$L^2.$ The proof relies upon a scale invariant Poincare inequality that we
establish for a large class of weighted Dirichlet forms, as well as estimates
to handle certain mildly singular perturbation terms. The weights that we
consider are neither Ahlfors regular, nor do they generally belong to the
Muckenhaupt class $A_2.$


  Background/Purpose: In the context of the European Biodiversity policy, the
Green Infrastructure Strategy is one supporting tool to mitigate fragmentation,
inter-alia to increase the spatial and functional connectivity between
protected and unprotected areas. The Joint Research Centre has developed an
integrated model to provide a macro-scale set of indices to evaluate the
connectivity of the Natura 2000 network, which forms the backbone of a Green
Infrastructure for Europe. The model allows a wide assessment and comparison to
be performed across countries in terms of structural (spatially connected or
isolated sites) and functional connectivity (least-cost distances between sites
influenced by distribution, distance and land cover).
  Main conclusion: The Natura 2000 network in Europe shows differences among
countries in terms of the sizes and numbers of sites, their distribution as
well as distances between sites. Connectivity has been assessed on the basis of
a 500 m average inter-site distance, roads and intensive land use as barrier
effects as well as the presence of "green" corridors. In all countries the
Natura 2000 network is mostly made of sites which are not physically connected.
Highest functional connectivity values are found for Spain, Slovakia, Romania
and Bulgaria. The more natural landscape in Sweden and Finland does not result
in high inter-site network connectivity due to large inter-site distances. The
distribution of subnets with respect to roads explains the higher share of
isolated subnets in Portugal than in Belgium.


  The development of reliable methods for restoring susceptibility after
antibiotic resistance arises has proven elusive. A greater understanding of the
relationship between antibiotic administration and the evolution of resistance
is key to overcoming this challenge. Here we present a data-driven mathematical
approach for developing antibiotic treatment plans that can reverse the
evolution of antibiotic resistance determinants. We have generated adaptive
landscapes for 16 genotypes of the TEM beta-lactamase that vary from the wild
type genotype TEM-1 through all combinations of four amino acid substitutions.
We determined the growth rate of each genotype when treated with each of 15
beta-lactam antibiotics. By using growth rates as a measure of fitness, we
computed the probability of each amino acid substitution in each beta-lactam
treatment using two different models named the Correlated Probability Model
(CPM) and the Equal Probability Model (EPM). We then performed an exhaustive
search through the 15 treatments for substitution paths leading from each of
the 16 genotypes back to the wild type TEM-1. We identified those treatment
paths that returned the highest probabilities of selecting for reversions of
amino acid substitutions and returning TEM to the wild type state. For the CPM
model, the optimized probabilities ranged between 0.6 and 1.0. For the EPM
model, the optimized probabilities ranged between 0.38 and 1.0. For cyclical
CPM treatment plans in which the starting and ending genotype was the wild
type, the probabilities were between 0.62 and 0.7. Overall this study shows
that there is promise for reversing the evolution of resistance through
antibiotic treatment plans.


  Diffusion processes on trees are commonly used in evolutionary biology to
model the joint distribution of continuous traits, such as body mass, across
species. Estimating the parameters of such processes from tip values presents
challenges because of the intrinsic correlation between the observations
produced by the shared evolutionary history, thus violating the standard
independence assumption of large-sample theory. For instance Ho and An\'e
\cite{HoAne13} recently proved that the mean (also known in this context as
selection optimum) of an Ornstein-Uhlenbeck process on a tree cannot be
estimated consistently from an increasing number of tip observations if the
tree height is bounded. Here, using a fruitful connection to the so-called
reconstruction problem in probability theory, we study the convergence rate of
parameter estimation in the unbounded height case. For the mean of the process,
we provide a necessary and sufficient condition for the consistency of the
maximum likelihood estimator (MLE) and establish a phase transition on its
convergence rate in terms of the growth of the tree. In particular we show that
a loss of $\sqrt{n}$-consistency (i.e., the variance of the MLE becomes
$\Omega(n^{-1})$, where $n$ is the number of tips) occurs when the tree growth
is larger than a threshold related to the phase transition of the
reconstruction problem. For the covariance parameters, we give a novel,
efficient estimation method which achieves $\sqrt{n}$-consistency under natural
assumptions on the tree.


  Plants create sugar in the mesophyll cells of their leaves by photosynthesis.
This sugar, mostly sucrose, has to be loaded via the bundle sheath into the
phloem vascular system (the sieve elements), where it is distributed to growing
parts of the plant. We analyze the feasibility of a particular loading
mechanism, active symplasmic loading, also called the polymer trap mechanism,
where sucrose is transformed into heavier sugars, such as raffinose and
stachyose, in the intermediary-type companion cells bordering the sieve
elements in the minor veins of the phloem. Keeping the heavier sugars from
diffusing back requires that the plasmodesmata connecting the bundle sheath
with the intermediary cell act as extremely precise filters, which are able to
distinguish between molecules that differ by less than 20% in size. In our
modeling, we take into account the coupled water and sugar movement across the
relevant interfaces, without explicitly considering the chemical reactions
transforming the sucrose into the heavier sugars. Based on the available data
for plasmodesmata geometry, sugar concentrations and flux rates, we conclude
that this mechanism can in principle function. We find that the water flow
through the plasmodesmata, which has not been quantified before, contributes
only 10-20% to the sucrose flux into the intermediary cells, while the main
part is transported by diffusion. On the other hand, the subsequent sugar
translocation into the sieve elements would very likely be carried
predominantly by bulk water flow through the plasmodesmata. Thus, in contrast
to apoplasmic loaders, all the necessary water for phloem translocation would
be supplied in this way with no need for additional water uptake across the
plasma membranes of the phloem.


  We present a simple model to study L\'{e}vy-flight foraging in a finite
landscape with countable targets. In our approach, foraging is a step-based
exploratory random search process with a power-law step-size distribution $P(l)
\propto l^{-\mu}$. We find that, when the termination is regulated by a finite
number of steps $N$, the optimum value of $\mu$ that maximises the foraging
efficiency can vary substantially in the interval $\mu \in (1,3)$, depending on
the landscape features (landscape size and number of targets). We further
demonstrate that subjective returning can be another significant factor that
affects the foraging efficiency in such context. Our results suggest that
L\'{e}vy-flight foraging may arise through an interaction between the
environmental context and the termination of exploitation, and particularly
that the number of steps can play an important role in this scenario which is
overlooked by most previous work. Our study not only provides a new perspective
on L\'{e}vy-flight foraging, but also opens new avenues for investigating the
interaction between foraging dynamics and environment as well as offers a
realistic framework for analysing animal movement patterns from empirical data.


  Many species of plants are found in regions to which they are alien and their
global distribution has been found to exhibit several remarkable
patterns,characterised by exponential functions of the kind that could arise
through versions of MacArthur's broken stick. We show here that these various
patterns are all quantitatively reproduced by a simple algorithm, in terms of a
single parameter- a single stick to be broken. This algorithm admits a
biological interpretation in terms of niche structures fluctuating with time
and productivity; with sites and species highly idiosyncratic. Technically,
this is an application of statistical mechanics to ecology quite different from
the familiar application to species abundance distributions.


  We have performed single molecule stretching experiments and dynamic light
scattering (DLS) in order to characterize the interaction between the DNA
molecule and the fluorescent stain GelRed. The results from single molecule
stretching show that the persistence length of the DNA-GelRed complexes
increases as the ligand concentration increases up to some critical
concentration, then decreasing for higher concentrations. The contour length of
the complexes, on the other hand, increases monotonically as a function of
GelRed concentration, suggesting that intercalation is the main binding
mechanism. In order to characterize the phys- ical chemistry of the
interaction, we use the McGhee-von Hippel binding isotherm to extract the
physicochemical parameters of the interaction from the contour length data. The
DLS experiments were performed to study the changes of the effective size of
the DNA-GelRed complexes, measured by the hydrodynamic radius, as a function of
ligand concentration. We found a qualitative agreement between the results
obtained from the two techniques by com- paring the behaviors of the
hydrodynamics radius and the radius of gyration, since this last quantity can
be expressed as a function of mechanical parameters determined from the
stretching experiments.


  Data-based mathematical modeling of biochemical reaction networks, e.g. by
nonlinear ordinary differential equation (ODE) models, has been successfully
applied. In this context, parameter estimation and uncertainty analysis is a
major task in order to assess the quality of the description of the system by
the model. Recently, a broadened eigenvalue spectrum of the Hessian matrix of
the objective function covering orders of magnitudes was observed and has been
termed as sloppiness. In this work, we investigate the origin of sloppiness
from structures in the sensitivity matrix arising from the properties of the
model topology and the experimental design. Furthermore, we present strategies
using optimal experimental design methods in order to circumvent the sloppiness
issue and present non-sloppy designs for a benchmark model.


  We study a general class of birth-and-death processes with state space
$\mathbb{N}$ that describes the size of a population going to extinction with
probability one. This class contains the logistic case. The scale of the
population is measured in terms of a `carrying capacity' $K$. When $K$ is
large, the process is expected to stay close to its deterministic equilibrium
during a long time but ultimately goes extinct. Our aim is to quantify the
behavior of the process and the mean time to extinction in the quasi-stationary
distribution as a function of $K$, for large $K$. We also give a quantitative
description of this quasi-stationary distribution. It turns out to be close to
a Gaussian distribution centered about the deterministic long-time equilibrium,
when $K$ is large. Our analysis relies on precise estimates of the maximal
eigenvalue, of the corresponding eigenvector and of the spectral gap of a
self-adjoint operator associated with the semigroup of the process.


  Traditionally, the evolution of cooperation has been studied on single,
isolated networks. Yet a player, especially in human societies, will typically
be a member of many different networks, and those networks will play a
different role in the evolutionary process. Multilayer networks are therefore
rapidly gaining on popularity as the more apt description of a networked
society. With this motivation, we here consider 2-layer scale-free networks
with all possible combinations of degree mixing, wherein one network layer is
used for the accumulation of payoffs and the other is used for strategy
updating. We find that breaking the symmetry through assortative mixing in one
layer and/or disassortative mixing in the other layer, as well as preserving
the symmetry by means of assortative mixing in both layers, impedes the
evolution of cooperation. We use degree-dependent distributions of strategies
and cluster-size analysis to explain these results, which highlight the
importance of hubs and the preservation of symmetry between multilayer networks
for the successful resolution of social dilemmas.


  The contemporary forestry preoccupation with non-declining even-flow during
yield simulations detracts from more important questions about the constraints
that should bind the end of a simulation. Whilst long simulations help to
convey a sense of sustainability, they are inferior to stronger indicators such
as the optimal state and binding conditions at the end of a simulation.
Rigorous definitions of sustainability that constrain the terminal state should
allow flexibility in the planning horizon and relaxation of non-declining
even-flow, allowing both greater economic efficiency and better environmental
outcomes. Suitable definitions cannot be divorced from forest type and
management objectives, but should embrace concepts that ensure the anticipated
value of the next harvest, the continuity of growing stock, and in the case of
uneven-aged management, the adequacy of regeneration.


  The Disrupted in Schizophrenia Gene 1 (DISC1) plays a role in both neural
signalling and development and is associated with schizophrenia, although its
links to altered brain structure and function in this disorder are not fully
established. Here we have used structural and functional MRI to investigate
links with six DISC1 single nucleotide polymorphisms (SNPs). We employed a
brain-wide association analysis (BWAS) together with a Jacknife internal
validation approach in 46 schizophrenia patients and 24 matched healthy control
subjects. Results from structural MRI showed significant associations between
all six DISC1 variants and gray matter volume in the precuneus, post-central
gyrus and middle cingulate gyrus. Associations with specific SNPs were found
for rs2738880 in the left precuneus and right post-central gyrus, and rs1535530
in the right precuneus and middle cingulate gyrus. Using regions showing
structural associations as seeds a resting-state functional connectivity
analysis revealed significant associations between all 6 SNPS and connectivity
between the right precuneus and inferior frontal gyrus. The connection between
the right precuneus and inferior frontal gyrus was also specifically associated
with rs821617. Importantly schizophrenia patients showed positive correlations
between the six DISC-1 SNPs associated gray matter volume in the left precuneus
and right post-central gyrus and negative symptom severity. No correlations
with illness duration were found. Our results provide the first evidence
suggesting a key role for structural and functional connectivity associations
between DISC1 polymorphisms and the precuneus in schizophrenia.


  Adenosine 5'-triphosphate (ATP) is the nearly ubiquitous "energy currency" of
living organisms, and thus is a crucial participant in the majority of
enzymatic reactions. The standard models in enzyme kinetics generally ignore
the temporal dynamics of ATP because it is assumed to be present in large
excess. However, this assumption may not hold in many situations of cellular
stress where ATP concentrations may be comparable to substrate levels. Here, we
demonstrate the importance of ATP concentration on the dynamics of multi-enzyme
reactions by explicit consideration of ATP as a secondary substrate for an
enzyme. We apply our model to the mitogen-activated protein (MAP) kinase
cascade, which is involved in the regulation of a vast range of cellular
activities. We show that three fundamental features of this signaling network
--- (i) duration of response, (ii) signal amplification, and (iii)
ultrasensitivity to stimulus concentration --- are all dependent on ATP
concentration. Our results indicate that the concentration of ATP regulates the
response of the MAP kinase activation network, and potentially suggests another
possible mechanism for disruption of the cascade in pathogenic states.


  Calmodulin, the ubiquitous calcium-activated second messenger in eukaryotes,
is an extremely versatile molecule involved in many biological processes:
muscular contraction, synaptic plasticity, circadian rhythm, and cell cycle,
among others. The protein is structurally organised into two globular lobes,
joined by a flexible linker. Calcium modulates calmodulin activity by favoring
a conformational transition of each lobe from a closed conformation to an open
conformation. Most targets have a strong preference for one conformation over
the other, and depending on the free calcium concentration in a cell,
particular sets of targets will preferentially interact with calmodulin. In
turn, targets can increase or decrease the calcium affinity of the calmodulin
molecules to which they bind. Interestingly, experiments with the tryptic
fragments showed that most targets have a much lower affinity for the N-lobe
than for the C-lobe. Hence, the latter predominates in the formation of most
calmodulin-target complexes. We showed that a relatively simple allosteric
mechanism, based the classic MWC model, can capture the observed modulation of
both the isolated C-lobe, and intact calmodulin, by individual targets.
Moreover, our model can be naturally extended to study how the calcium affinity
of a single pool of calmodulin is modulated by a mixture of competing targets
in vivo.


  In biochemical networks, reactions often occur on disparate timescales and
can be characterized as either "fast" or "slow." The quasi-steady state
approximation (QSSA) utilizes timescale separation to project models of
biochemical networks onto lower-dimensional slow manifolds. As a result, fast
elementary reactions are not modeled explicitly, and their effect is captured
by non-elementary reaction rate functions (e.g. Hill functions). The accuracy
of the QSSA applied to deterministic systems depends on how well timescales are
separated. Recently, it has been proposed to use the non-elementary rate
functions obtained via the deterministic QSSA to define propensity functions in
stochastic simulations of biochemical networks. In this approach, termed the
stochastic QSSA, fast reactions that are part of non-elementary reactions are
not simulated, greatly reducing computation time. However, it is unclear when
the stochastic QSSA provides an accurate approximation of the original
stochastic simulation. We show that, unlike the deterministic QSSA, the
validity of the stochastic QSSA does not follow from timescale separation
alone, but also depends on the sensitivity of the non-elementary reaction rate
functions to changes in the slow species. The stochastic QSSA becomes more
accurate when this sensitivity is small. Different types of QSSAs result in
non-elementary functions with different sensitivities, and the total QSSA
results in less sensitive functions than the standard or the pre-factor QSSA.
We prove that, as a result, the stochastic QSSA becomes more accurate when
non-elementary reaction functions are obtained using the total QSSA. Our work
provides a novel condition for the validity of the QSSA in stochastic
simulations of biochemical reaction networks with disparate timescales.


  Population dynamics are often subject to random independent changes in the
environment. For the two strategy stochastic replicator dynamic, we assume that
stochastic changes in the environment replace the payoffs and variance. This is
modeled by a continuous time Markov chain in a finite atom space. We establish
conditions for this dynamic to have an analogous characterization of the
long-run behavior to that of the deterministic dynamic. To create intuition, we
first consider the case when the Markov chain has two states. A very natural
extension to the general finite state space of the Markov chain will be given.


  In this paper, we study an adaptive spatial network. We consider an SIS
(susceptible-infectedsusceptible) epidemic on the network, with a link/contact
rewiring process constrained by spatial proximity. In particular, we assume
that susceptible nodes break links with infected nodes independently of
distance, and reconnect at random to susceptible nodes available within a given
radius. By systematically manipulating this radius we investigate the impact of
rewiring on the structure of the network and characteristics of the epidemic.
We adopt a step-by-step approach whereby we first study the impact of rewiring
on the network structure in the absence of an epidemic, then with nodes
assigned a disease status but without disease dynamics, and finally running
network and epidemic dynamics simultaneously. In the case of no labelling and
no epidemic dynamics, we provide both analytic and semi-analytic formulas for
the value of clustering achieved in the network. Our results also show that the
rewiring radius and the network's initial structure have a pronounced effect on
the endemic equilibrium, with increasingly large rewiring radiuses yielding
smaller disease prevalence.


  We consider the effects of social learning on the individual learning and
genetic evolution of a colony of artificial agents capable of genetic,
individual and social modes of adaptation. We confirm that there is strong
selection pressure to acquire traits of individual learning and social learning
when these are adaptive traits. We show that selection pressure for learning of
either kind can supress selection pressure for reproduction or greater fitness.
We show that social learning differs from individual learning in that it can
support a second evolutionary system that is decoupled from the biological
evolutionary system. This decoupling leads to an emergent interaction where
immature agents are more likely to engage in learning activities than mature
agents.


  We present a theoretical model for the description of the adsorption kinetics
of globular proteins onto charged core-shell microgel particles based on
Dynamic Density Functional Theory (DDFT). This model builds on a previous
description of protein adsorption thermodynamics [Yigit \textit{et al},
Langmuir 28 (2012)], shown to well interpret the available calorimetric
experimental data of binding isotherms. In practice, a spatially-dependent
free-energy functional including the same physical interactions is built, and
used to study the kinetics via a generalised diffusion equation. To test this
model, we apply it to the case study of Lysozyme adsorption on PNIPAM coated
nanoparticles, and show that the dynamics obtained within DDFT is consistent
with that extrapolated from experiments. We also perform a systematic study of
the effect of various parameters in our model, and investigate the loading
dynamics as a function of proteins' valence and hydrophobic adsorption energy,
as well as their concentration and that of the nanoparticles. Although we
concentrated here on the case of adsorption for a single protein type, the
model's generality allows to study multi-component system, providing a reliable
instrument for future studies of competitive and cooperative adsorption effects
often encountered in protein adsorption experiments.


  DNA looping has been observed to enhance and suppress transcriptional noise
but it is uncertain which of these two opposite effects is to be expected for
given conditions. Here, we derive analytical expressions for the main
quantifiers of transcriptional noise in terms of the molecular parameters and
elucidate the role of DNA looping. Our results rationalize paradoxical
experimental observations and provide the first quantitative explanation of
landmark individual-cell measurements at the single molecule level on the
classical lac operon genetic system [Choi et al., Science 322, 442-446 (2008)].


  Understanding the forces shaping ecological communities is crucially
important to basic science and conservation. In recent years, considerable
progress was made in explaining communities using simple and general models,
with neutral theory as a prominent example. However, while successful in
explaining static patterns such as species abundance distributions, the neutral
theory was criticized for making unrealistic predictions of fundamental dynamic
patterns. Here we incorporate environmental stochasticity into the neutral
framework, and show that the resulting generalized neutral theory is capable of
predicting realistic patterns of both population and community dynamics.
Applying the theory to real data (the tropical forest of Barro-Colorado
Island), we find that it better fits the observed distribution of short-term
fluctuations, the temporal scaling of such fluctuations, and the decay of
compositional similarity with time, than the original theory, while retaining
its power to explain static patterns of species abundance. Importantly,
although the proposed theory is neutral (all species are functionally
equivalent) and stochastic, it is a niche-based theory in the sense that
species differ in their demographic responses to environmental variation. Our
results show that this integration of niche forces and stochasticity within a
minimalistic neutral framework is highly successful in explaining fundamental
static and dynamic characteristics of ecological communities.


  Many aspects of macroevolutionary theory and our understanding of biotic
responses to global environmental change derive from literature-based
compilations of palaeontological data. Existing manually assembled databases
are, however, incomplete and difficult to assess and enhance. Here, we develop
and validate the quality of a machine reading system, PaleoDeepDive, that
automatically locates and extracts data from heterogeneous text, tables, and
figures in publications. PaleoDeepDive performs comparably to humans in complex
data extraction and inference tasks and generates congruent synthetic
macroevolutionary results. Unlike traditional databases, PaleoDeepDive produces
a probabilistic database that systematically improves as information is added.
We also show that the system can readily accommodate sophisticated data types,
such as morphological data in biological illustrations and associated textual
descriptions. Our machine reading approach to scientific data integration and
synthesis brings within reach many questions that are currently underdetermined
and does so in ways that may stimulate entirely new modes of inquiry.


  A central challenge in the study of protein evolution is the identification
of historic amino acid sequence changes responsible for creating novel
functions observed in present-day proteins. To address this problem, we
developed a new method to identify and rank amino acid mutations in ancestral
protein sequences according to their function-shifting potential. Our approach
scans the changes between two reconstructed ancestral sequences in order to
find (1) sites with sequence changes that significantly deviate from our
model-based probabilistic expectations, (2) sites that demonstrate extreme
changes in mutual information, and (3) sites with extreme gains or losses of
information content. By taking the overlaps of these statistical signals, the
method accurately identifies cryptic evolutionary patterns that are often not
obvious when examining only the conservation of modern-day protein sequences.
We validated this method with a training set of previously-discovered
function-shifting mutations in three essential protein families in animals and
fungi, whose evolutionary histories were the prior subject of systematic
molecular biological investigation. Our method identified the known
function-shifting mutations in the training set with a very low rate of false
positive discovery. Further, our approach significantly outperformed other
methods that use variability in evolutionary rates to detect functional loci.
The accuracy of our approach indicates it could be a useful tool for generating
specific testable hypotheses regarding the acquisition of new functions across
a wide range of protein families.


  The existence of a die-out threshold (different from the classic
disease-invasion one) defining a region of slow extinction of an epidemic has
been proved elsewhere for susceptible-aware-infectious-susceptible models
without awareness decay, through bifurcation analysis. By means of an
equivalent mean-field model defined on regular random networks, we interpret
the dynamics of the system in this region and prove that the existence of
bifurcation for this second epidemic threshold crucially depends on the absence
of awareness decay. We show that the continuum of equilibria that characterizes
the slow die-out dynamics collapses into a unique equilibrium when a constant
rate of awareness decay is assumed, no matter how small, and that the resulting
bifurcation from the disease-free equilibrium is equivalent to that of standard
epidemic models. We illustrate these findings with continuous-time stochastic
simulations on regular random networks with different degrees. Finally, the
behaviour of solutions with and without decay in awareness is compared around
the second epidemic threshold for a small rate of awareness decay.


  In this manuscript we propose a mathematical framework to couple
transcription and translation in which mRNA production is described by a set of
master equations while the dynamics of protein density is governed by a random
differential equation. The coupling between the two processes is given by a
stochastic perturbation whose statistics satisfies the master equations. In
this approach, from the knowledge of the analytical time dependent distribution
of mRNA number, we are able to calculate the dynamics of the probability
density of the protein population.


  Wolumes is a fast and stand-alone computer program written in standard C that
allows the measure of atom volumes in proteins. Its algorithm is a simple
discretization of the space by means of a grid of points at 0.75 Angstroms from
each other and it uses a set of van der Waals radii optimized for protein
atoms. By comparing the computed values with distributions derived from a
non-redundant subset of the Protein Data Bank, the new methods allows to
identify atoms and residues abnormally large/small. The source code is freely
available, together with some examples.


  The detection of local genomic signals using high-throughput DNA sequencing
data can be cast as a problem of scanning a Poisson random field for local
changes in the rate of the process. We propose a likelihood-based framework for
for such scans, and derive formulas for false positive rate control and power
calculations. The framework can also accommodate mixtures of Poisson processes
to deal with over-dispersion. As a specific, detailed example, we consider the
detection of insertions and deletions by paired-end DNA-sequencing. We propose
several statistics for this problem, compare their power under current
experimental designs, and illustrate their application on an Illumina Platinum
Genomes data set.


  Cellular signaling involves the transmission of environmental information
through cascades of stochastic biochemical reactions, inevitably introducing
noise that compromises signal fidelity. Each stage of the cascade often takes
the form of a kinase-phosphatase push-pull network, a basic unit of signaling
pathways whose malfunction is linked with a host of cancers. We show this
ubiquitous enzymatic network motif effectively behaves as a Wiener-Kolmogorov
(WK) optimal noise filter. Using concepts from umbral calculus, we generalize
the linear WK theory, originally introduced in the context of communication and
control engineering, to take nonlinear signal transduction and discrete
molecule populations into account. This allows us to derive rigorous
constraints for efficient noise reduction in this biochemical system. Our
mathematical formalism yields bounds on filter performance in cases important
to cellular function---like ultrasensitive response to stimuli. We highlight
features of the system relevant for optimizing filter efficiency, encoded in a
single, measurable, dimensionless parameter. Our theory, which describes noise
control in a large class of signal transduction networks, is also useful both
for the design of synthetic biochemical signaling pathways, and the
manipulation of pathways through experimental probes like oscillatory input.


  We study the competition and the evolution of nodes embedded in Euclidean
restricted spaces. The population evolves by a branching process in which new
nodes are generated when up to two new nodes are attached to the previous ones
at each time unit. The competition in the population is introduced by
considering the effect of overcrowding of nodes in the embedding space. The
branching process is suppressed if the newborn node is closer than a distance
$\xi$ of the previous nodes. This rule may be relevant to describe a
competition for resources, limiting the density of individuals and therefore
the total population. This results in an exponential growth in the initial
period, and, after some crossover time, approaching some limiting value. Our
results show that the competition among the nodes associated with geometric
restrictions can even, for certain conditions, lead the entire population to
extinction.


  We study the rock-paper-scissors game in structured populations, where the
invasion rates determine individual payoffs that govern the process of strategy
change. The traditional version of the game is recovered if the payoffs for
each potential invasion stem from a single pairwise interaction. However, the
transformation of invasion rates to payoffs also allows the usage of larger
interaction ranges. In addition to the traditional pairwise interaction, we
therefore consider simultaneous interactions with all nearest neighbors, as
well as with all nearest and next-nearest neighbors, thus effectively going
from single pair to group interactions in games of cyclic dominance. We show
that differences in the interaction range affect not only the stationary
fractions of strategies, but also their relations of dominance. The transition
from pairwise to group interactions can thus decelerate and even revert the
direction of the invasion between the competing strategies. Like in
evolutionary social dilemmas, in games of cyclic dominance too the indirect
multipoint interactions that are due to group interactions hence play a pivotal
role. Our results indicate that, in addition to the invasion rates, the
interaction range is at least as important for the maintenance of biodiversity
among cyclically competing strategies.


  Mechanical forces acting on cell adhesion receptor proteins regulate a range
of cellular functions by formation and rupture of non-covalent interactions
with ligands. Typically, force decreases the lifetimes of intact complexes
(slip-bonds), making the discovery that these lifetimes can also be prolonged
("catch-bonds"), a surprise. We created a microscopic analytic theory by
incorporating the structures of selectin and integrin receptors into a
conceptual framework based on the theory of stochastic equations, which
quantitatively explains a wide range of experimental data (including
catch-bonds at low forces and slip-bonds at high forces). Catch-bonds arise due
to force-induced remodeling of hydrogen bond networks, a finding that also
accounts for unbinding in structurally unrelated integrin-fibronectin and
actomyosin complexes. For the selectin family, remodeling of hydrogen bond
networks drives an allosteric transition resulting in the formation of maximum
number of hydrogen bonds determined only by the structure of the receptor and
is independent of the ligand. A similar transition allows us to predict the
increase in number of hydrogen bonds in a particular allosteric state of
$\alpha_5 \beta_1$ integrin-fibronectin complex, a conformation which is yet to
be crystallized. We also make a testable prediction that a single point
mutation (Tyr51Phe) in the ligand associated with selectin should dramatically
alter the nature of the catch-bond compared to the wild type. Our work suggests
that nature utilizes a ductile network of hydrogen bonds to engineer function
over a broad range of forces.


  Yeast cells grown in culture can spontaneously synchronize their respiration,
metabolism, gene expression and cell division. Such metabolic oscillations in
synchronized cultures reflect single-cell oscillations, but the relationship
between the oscillations in single cells and synchronized cultures is poorly
understood. To understand this relationship and the coordination between
metabolism and cell division, we collected and analyzed DNA-content,
gene-expression and physiological data, at hundreds of time-points, from
cultures metabolically-synchronized at different growth rates, carbon sources
and biomass densities. The data enabled us to extend and generalize an
ensemble-average-over-phases (EAP) model that connects the population-average
gene-expression of asynchronous cultures to the gene-expression dynamics in the
single-cells comprising the cultures. The extended model explains the
carbon-source specific growth-rate responses of hundreds of genes. Our data
demonstrate that for a given growth rate, the frequency of metabolic cycling in
synchronized cultures increases with the biomass density. This observation
underscores the difference between metabolic cycling in synchronized cultures
and in single cells and suggests entraining of the single-cell cycle by a
quorum-sensing mechanism. Constant levels of residual glucose during the
metabolic cycling of synchronized cultures indicate that storage carbohydrates
are required to fuel not only the G1/S transition of the division cycle but
also the metabolic cycle. Despite the large variation in profiled conditions
and in the time-scale of their dynamics, most genes preserve invariant dynamics
of coordination with each other and with the rate of oxygen consumption.
Similarly, the G1/S transition always occurs at the beginning, middle or end of
the high oxygen consumption phases, analogous to observations in human and
drosophila cells.


  The currently existing theory of fluorescence correlation spectroscopy(FCS)
is based on the linear fluctuation theory originally developed by Einstein,
Onsager, Lax, and others as a phenomenological approach to equilibrium
fluctuations in bulk solutions. For mesoscopic reaction-diffusion systems with
nonlinear chemical reactions among a small number of molecules, a situation
often encountered in single-cell biochemistry, it is expected that FCS time
correlation functions of a reaction-diffusion system can deviate from the
classic results of Elson and Magde. We first discuss this nonlinear effect for
reaction systems without diffusion. For nonlinear stochastic reaction-diffusion
systems here are no closed solutions; therefore, stochastic Monte-Carlo
simulations are carried out. We show that the deviation is small for a simple
bimolecular reaction; the most significant deviations occur when the number of
molecules is small and of the same order. Our results show that current linear
FCS theory could be adequate for measurements on biological systems that
contain many other sources of uncertainties. At the same time it provides a
framework for future measurements of nonlinear, fluctuating chemical reactions
with high-precision FCS. Extending Delbr\"uck-Gillespie's theory for stochastic
nonlinear reactions with rapidly stirring to reaction-diffusion systems
provides a mesoscopic model for chemical and biochemical reactions at
nanometric and mesoscopic level such as a single biological cell.


  We review the optimal control of systems modeling the dynamics of
tuberculosis. Time dependent control functions are introduced in the
mathematical models, representing strategies for the improvement of the
treatment and cure of active infectious and/or latent individuals. Optimal
control theory allows then to find the optimal way to implement the strategies,
minimizing the number of infectious and/or latent individuals and keeping the
cost of implementation as low as possible. An optimal control problem is
proposed and solved, illustrating the procedure. Simulations show an effective
reduction in the number of infectious individuals.


  Donor T cell mediated graft vs. host effects may result from the aggregate
alloreactivity to minor histocompatibility antigens (mHA) presented by the HLA
in each donor-recipient pair (DRP) undergoing stem cell transplantation (SCT).
Whole exome sequencing has demonstrated extensive nucleotide sequence variation
in HLA-matched DRP. Non-synonymous single nucleotide polymorphisms (nsSNPs) in
the GVH direction (polymorphisms present in recipient and absent in donor) were
identified in 4 HLA-matched related and 5 unrelated DRP. The nucleotide
sequence flanking each SNP was obtained utilizing the ANNOVAR software package.
All possible nonameric-peptides encoded by the non-synonymous SNP were then
interrogated in-silico for their likelihood to be presented by the HLA class I
molecules in individual DRP, using the Immune-Epitope Database (IEDB) SMM
algorithm. The IEDB-SMM algorithm predicted a median 18,396 peptides/DRP which
bound HLA with an IC50 of <500nM, and 2254 peptides/DRP with an IC50 of <50nM.
Unrelated donors generally had higher numbers of peptides presented by the HLA.
A similarly large library of presented peptides was identified when the data
was interrogated using the Net MHCPan algorithm. These peptides were uniformly
distributed in the various organ systems. The bioinformatic algorithm presented
here demonstrates that there may be a high level of minor histocompatibility
antigen variation in HLA-matched individuals, constituting an HLA-specific
alloreactivity potential. These data provide a possible explanation for how
relatively minor adjustments in GVHD prophylaxis yield relatively similar
outcomes in HLA matched and mismatched SCT recipients.


  Many common diseases have a complex genetic basis in which large numbers of
genetic variations combine with environmental and lifestyle factors to
determine risk. However, quantifying such polygenic effects and their
relationship to disease risk has been challenging. In order to address these
difficulties we developed a global measure of the information content of an
individual's genome relative to a reference population, which may be used to
assess differences in global genome structure between cases and appropriate
controls. Informally this measure, which we call relative genome information
(RGI), quantifies the relative "disorder" of an individual's genome. In order
to test its ability to predict disease risk we used RGI to compare single
nucleotide polymorphism genotypes from two independent samples of women with
early-onset breast cancer with three independent sets of controls. We found
that RGI was significantly elevated in both sets of breast cancer cases in
comparison with all three sets of controls, with disease risk rising sharply
with RGI (odds ratio greater than 12 for the highest percentile RGI).
Furthermore, we found that these differences are not due to associations with
common variants at a small number of disease-associated loci, but rather are
due to the combined associations of thousands of markers distributed throughout
the genome. Our results indicate that the information content of an
individual's genome may be used to measure the risk of a complex disease, and
suggest that early-onset breast cancer has a strongly polygenic basis.


  Understanding which phenotypic traits are consistently correlated throughout
evolution is a highly pertinent problem in modern evolutionary biology. Here,
we propose a multivariate phylogenetic latent liability model for assessing the
correlation between multiple types of data, while simultaneously controlling
for their unknown shared evolutionary history informed through molecular
sequences. The latent formulation enables us to consider in a single model
combinations of continuous traits, discrete binary traits and discrete traits
with multiple ordered and unordered states. Previous approaches have
entertained a single data type generally along a fixed history, precluding
estimation of correlation between traits and ignoring uncertainty in the
history. We implement our model in a Bayesian phylogenetic framework, and
discuss inference techniques for hypothesis testing. Finally, we showcase the
method through applications to columbine flower morphology, antibiotic
resistance in Salmonella and epitope evolution in influenza.


  The modelling of evolutionary game dynamics in finite populations requires
microscopic processes that determine how strategies spread. The exact details
of these processes are often chosen without much further consideration.
Different types of microscopic models, including in particular fitness-based
selection rules and imitation-based dynamics, are often used as if they were
interchangeable. We challenge this view and investigate how robust these
choices on the micro-level really are. Focusing on a key macroscopic
observable, the probability for a single mutant to take over a population of
wild-type individuals, we show that there is a unique pair of a fitness-based
process and an imitation process leading to identical outcomes for arbitrary
games and for all intensities of selection. This highlights the perils of
making arbitrary choices at the micro-level without regard of the consequences
at the macro-level.


  We propose a computational method to measure and visualize interrelationships
among any number of DNA sequences allowing, for example, the examination of
hundreds or thousands of complete mitochondrial genomes. An "image distance" is
computed for each pair of graphical representations of DNA sequences, and the
distances are visualized as a Molecular Distance Map: Each point on the map
represents a DNA sequence, and the spatial proximity between any two points
reflects the degree of structural similarity between the corresponding
sequences. The graphical representation of DNA sequences utilized, Chaos Game
Representation (CGR), is genome- and species-specific and can thus act as a
genomic signature. Consequently, Molecular Distance Maps could inform species
identification, taxonomic classifications and, to a certain extent,
evolutionary history. The image distance employed, Structural Dissimilarity
Index (DSSIM), implicitly compares the occurrences of oligomers of length up to
$k$ (herein $k=9$) in DNA sequences. We computed DSSIM distances for more than
5 million pairs of complete mitochondrial genomes, and used Multi-Dimensional
Scaling (MDS) to obtain Molecular Distance Maps that visually display the
sequence relatedness in various subsets, at different taxonomic levels. This
general-purpose method does not require DNA sequence homology and can thus be
used to compare similar or vastly different DNA sequences, genomic or
computer-generated, of the same or different lengths. We illustrate potential
uses of this approach by applying it to several taxonomic subsets: phylum
Vertebrata, (super)kingdom Protista, classes Amphibia-Insecta-Mammalia, class
Amphibia, and order Primates. This analysis of an extensive dataset confirms
that the oligomer composition of full mtDNA sequences can be a source of
taxonomic information.


  We propose an SEIR model for the populations and an SEI model for the vector
to describe the transmission dynamics of a four-strain model with both primary
and secondary dengue infections. In order to accomplish this, we propose and
obtain an analytic solution of a system of 47 coupled differential equations.
This would be the most complete epidemic model proposed to describe the dengue
epidemic.


  Nanopore sequencing of DNA is a single-molecule technique that may achieve
long reads, low cost, and high speed with minimal sample preparation and
instrumentation. Here, we build on recent progress with respect to nanopore
resolution and DNA control to interpret the procession of ion current levels
observed during the translocation of DNA through the pore MspA. As
approximately four nucleotides affect the ion current of each level, we
measured the ion current corresponding to all 256 four-nucleotide combinations
(quadromers). This quadromer map is highly predictive of ion current levels of
previously unmeasured sequences derived from the bacteriophage phi X 174
genome. Furthermore, we show nanopore sequencing reads of phi X 174 up to 4,500
bases in length that can be unambiguously aligned to the phi X 174 reference
genome, and demonstrate proof-of-concept utility with respect to hybrid genome
assembly and polymorphism detection. All methods and data are made fully
available.


  Actions such as sharing food and cooperating to reach a common goal have
played a fundamental role in the evolution of human societies. These good
actions may not maximise the actor's payoff, but they maximise the other's
payoff. Consequently, their existence is puzzling for evolutionary theories.
Why should you make an effort to help others, even when no reward seems to be
at stake? Indeed, experiments typically show that humans are heterogeneous:
some may help others, while others may not. With the aim of favouring the
emergence of 'successful cultures', a number of studies has recently
investigated what mechanisms promote the evolution of a particular good action.
But still little is known about if and how good actions can spread from person
to person. For instance, does being recipient of an altruistic act increase
your probability of being cooperative with others? Plato's quote, 'Good actions
give strength to ourselves and inspire good actions in others', suggests that
is possible. We have conducted an experiment on Amazon Mechanical Turk to test
this mechanism using economic games. We have measured willingness to be
cooperative through a standard Prisoner's dilemma and willingness to act
altruistically using a binary Dictator game. In the baseline treatments, the
endowments needed to play were given by the experimenters, as usual; in the
control treatments, they came from a good action made by someone else. Across
four different comparisons and a total of 572 subjects, we have never found a
significant increase of cooperation or altruism when the endowment came from a
good action. We conclude that good actions do not necessarily inspire good
actions in others, at least in the ideal scenario of a lab experiment with
anonymous subjects.


  Protein folding and design are major biophysical problems, the solution of
which would lead to important applications especially in medicine. Here a novel
protein model capable of simultaneously provide quantitative protein design and
folding is introduced. With computer simulations it is shown that, for a large
set of real protein structures, the model produces designed sequences with
similar physical properties to the corresponding natural occurring sequences.
The designed sequences are not yet fully realistic and require further
experimental testing. For an independent set of proteins, notoriously difficult
to fold, the correct folding of both the designed and the natural sequences is
also demonstrated. The folding properties are characterized by free energy
calculations. which not only are consistent among natural and designed
proteins, but we also show a remarkable precision when the folded structures
are compared to the experimentally determined ones. Ultimately, this novel
coarse-grained protein model is unique in the combination of its fundamental
three features: its simplicity, its ability to produce natural foldable
designed sequences, and its structure prediction precision. The latter
demonstrated by free energy calculations. It is also remarkable that low
frustration sequences can be obtained with such a simple and universal design
procedure, and that the folding of natural proteins shows funnelled free energy
landscapes without the need of any potentials based on the native structure.


  Motivation:
  FASTQ is a standard file format for DNA sequencing data which stores both
nucleotides and quality scores. A typical sequencing study can easily generate
hundreds of gigabytes of FASTQ files, while public archives such as ENA and
NCBI and large international collaborations such as the Cancer Genome Atlas can
accumulate many terabytes of data in this format. Text compression tools such
as gzip are often employed to reduce the storage burden, but have the
disadvantage that the data must be decompressed before it can be used.
  Here we present BEETL-fastq, a tool that not only compresses FASTQ-formatted
DNA reads more compactly than gzip, but also permits rapid search for $k$-mer
queries within the archived sequences. Importantly, the full FASTQ record of
each matching read or read pair is returned, allowing the search results to be
piped directly to any of the many standard tools that accept FASTQ data as
input.
  Results:
  We show that 6.6 terabytes of human reads in FASTQ format can be transformed
into 1.7 terabytes of indexed files, from where we can search for 1, 10, 100,
1000, a million of 30-mers in respectively 3, 8, 14, 45 and 567 seconds plus 20
ms per output read. Useful applications of the search capability are
highlighted, including the genotyping of structural variant breakpoints and "in
silico pull-down" experiments in which only the reads that cover a region of
interest are selectively extracted for the purposes of variant calling or
visualization.
  Availability:
  BEETL-fastq is part of the BEETL library, available as a github repository at
git@github.com:BEETL/BEETL.git.


  Evolutionary dynamics on graphs can lead to many interesting and
counterintuitive findings. We study the Moran process, a discrete time
birth-death process, that describes the invasion of a mutant type into a
population of wild-type individuals. Remarkably, the fixation probability of a
single mutant is the same on all regular networks. But non-regular networks can
increase or decrease the fixation probability. While the time until fixation
formally depends on the same transition probabilities as the fixation
probabilities, there is no obvious relation between them. For example, an
amplifier of selection, which increases the fixation probability and thus
decreases the number of mutations needed until one of them is successful, can
at the same time slow down the process of fixation. Based on small networks, we
show analytically that (i) the time to fixation can decrease when links are
removed from the network and (ii) the node providing the best starting
conditions in terms of the shortest fixation time depends on the fitness of the
mutant. Our results are obtained analytically on small networks, but numerical
simulations show that they are qualitatively valid even in much larger
populations.


  Using low molecular weight chitosan nanoparticles (CNPs) prepared by an ionic
gelation method, we report the effect of low-intensity pulsed ultrasound (US)
on cell viability and nanoparticle uptake in cultured murine pre-osteoblasts.
Particle size and zeta potential are measured using Dynamic Light Scattering
(DLS), and cell viability is evaluated using the MTS assay. Results show that
30 min delivery of CNPs at 0.5 mg/mL is able to prevent loss of cell viability
due to either serum starvation or subsequent exposure to US (1 W/cm2 or 2
W/cm2, up to 1 min). Additionally, flow cytometry data suggest that there is a
close association between cellular membrane integrity and the presence of CNPs
when US at 2 W/cm2 is administered.


  Phylogenetic analyses which include fossils or molecular sequences that are
sampled through time require models that allow one sample to be a direct
ancestor of another sample. As previously available phylogenetic inference
tools assume that all samples are tips, they do not allow for this possibility.
We have developed and implemented a Bayesian Markov Chain Monte Carlo (MCMC)
algorithm to infer what we call sampled ancestor trees, that is, trees in which
sampled individuals can be direct ancestors of other sampled individuals. We
use a family of birth-death models where individuals may remain in the tree
process after the sampling, in particular we extend the birth-death skyline
model [Stadler et al, 2013] to sampled ancestor trees. This method allows the
detection of sampled ancestors as well as estimation of the probability that an
individual will be removed from the process when it is sampled. We show that
sampled ancestor birth-death models where all samples come from different time
points are non-identifiable and thus require one parameter to be known in order
to infer other parameters. We apply this method to epidemiological data, where
the possibility of sampled ancestors enables us to identify individuals that
infected other individuals after being sampled and to infer fundamental
epidemiological parameters. We also apply the method to infer divergence times
and diversification rates when fossils are included among the species samples,
so that fossilisation events are modelled as a part of the tree branching
process. Such modelling has many advantages as argued in literature. The
sampler is available as an open-source BEAST2 package
(https://github.com/gavryushkina/sampled-ancestors).


  Rigid-body motions of protein secondary structure are often implicated in
mecha-nisms of protein function. In GPCRs, evidence suggests that intracellular
loop 2 (IL2) contains a segment characterized as a helix when the activated
receptor trig-gers downstream signaling. However, neither experimental nor
computational methods are readily available to assess quantitatively the degree
of collective mo-tions in such secondary structure motifs of proteins. Here we
describe a new el-ement of our N-body Information Theory (NbIT) framework to
address this problem. To this end we introduce total intercorrelation, a
measure in infor-mation theory that can be used to describe n-body correlated
dynamics between multivariate distributions, such as 3-dimensional atomic
fluctuations in simula-tions of proteins. We also define two additional
measures, the rigid-body fraction and correlation order, which can be
determined from the decomposition of the configurational entropy. Using these
measures, we analyze the dynamics of IL2 in microsecond Molecular Dynamics
simulations of the 5-HT2A receptor to demonstrate the powerful features of the
new analysis techniques in studying the collective motion dynamics of secondary
structure motifs. The analysis reveals an intriguing difference in the extent
of correlated motions in the helical segment of IL2 in the presence and absence
of bound 5-HT, the endogenous agonist that ac-tivates the receptor and triggers
downstream signaling, suggesting that IL2 rigid-body motions can display
distinct behaviors that may discriminate functional mechanism of GPCRs.


  The inference of the evolutionary history of a collection of organisms is a
problem of fundamental importance in evolutionary biology. The abundance of DNA
sequence data arising from genome sequencing projects has led to significant
challenges in the inference of these phylogenetic relationships. Among these
challenges is the inference of the evolutionary history of a collection of
species based on sequence information from several distinct genes sampled
throughout the genome. It is widely accepted that each individual gene has its
own phylogeny, which may not agree with the species tree. Many possible causes
of this gene tree incongruence are known. The best studied is incomplete
lineage sorting, which is commonly modeled by the coalescent process. Numerous
methods based on the coalescent process have been proposed for estimation of
the phylogenetic species tree given DNA sequence data. However, use of these
methods assumes that the phylogenetic species tree can be identified from DNA
sequence data at the leaves of the tree, although this has not been formally
established. We prove that the unrooted topology of the n-leaf phylogenetic
species tree is generically identifiable given observed data at the leaves of
the tree that are assumed to have arisen from the coalescent process under a
time-reversible substitution process with the possibility of site-specific rate
variation modeled by the discrete gamma distribution and a proportion of
invariable sites.


  The time variation of contacts in a networked system may fundamentally alter
the properties of spreading processes and affect the condition for large-scale
propagation, as encoded in the epidemic threshold. Despite the great interest
in the problem for the physics, applied mathematics, computer science and
epidemiology communities, a full theoretical understanding is still missing and
currently limited to the cases where the time-scale separation holds between
spreading and network dynamics or to specific temporal network models. We
consider a Markov chain description of the Susceptible-Infectious-Susceptible
process on an arbitrary temporal network. By adopting a multilayer perspective,
we develop a general analytical derivation of the epidemic threshold in terms
of the spectral radius of a matrix that encodes both network structure and
disease dynamics. The accuracy of the approach is confirmed on a set of
temporal models and empirical networks and against numerical results. In
addition, we explore how the threshold changes when varying the overall time of
observation of the temporal network, so as to provide insights on the optimal
time window for data collection of empirical temporal networked systems. Our
framework is both of fundamental and practical interest, as it offers novel
understanding of the interplay between temporal networks and spreading
dynamics.