  One of the most remarkable aspects of our species is that while we show
surprisingly little genetic diversity, we demonstrate astonishing amounts of
cultural diversity. Perhaps most impressive is the diversity of our
technologies, broadly defined as all the physical objects we produce and the
skills we use to produce them. Despite considerable focus on the evolution of
technology by social scientists and philosophers, there have been few attempts
to systematically quantify technological diversity and therefore the dynamics
of technological change remain poorly understood. Here we show a novel Bayesian
model for examining technological diversification adopted from paleontological
analysis of occurrence data. We use this framework to estimate the tempo of
diversification in American car and truck models produced between 1896 and 2014
and to test the relative importance of competition and extrinsic factors in
shaping changes in macroevolutionary rates. Our results identify a four-fold
decrease in the origination and extinction rates of car models and a negative
net diversification rate over the last thirty years. We also demonstrate that
competition played a more significant role in car model diversification than
either changes in oil prices or gross domestic product. Together our analyses
provide a set of tools that can enhance current research on technological and
cultural evolution by providing a flexible and quantitative framework for
exploring the dynamics of diversification.


  Bacterial growth depends crucially on metabolic fluxes, which are limited by
the cell's capacity to maintain metabolic enzymes. The necessary enzyme amount
per unit flux is a major determinant of metabolic strategies both in evolution
and bioengineering. It depends on enzyme parameters (such as kcat and KM
constants), but also on metabolite concentrations. Moreover, similar amounts of
different enzymes might incur different costs for the cell, depending on
enzyme-specific properties such as protein size and half-life. Here, we
developed enzyme cost minimization (ECM), a scalable method for computing
enzyme amounts that support a given metabolic flux at a minimal protein cost.
The complex interplay of enzyme and metabolite concentrations, e.g. through
thermodynamic driving forces and enzyme saturation, would make it hard to solve
this optimization problem directly. By treating enzyme cost as a function of
metabolite levels, we formulated ECM as a numerically tractable, convex
optimization problem. Its tiered approach allows for building models at
different levels of detail, depending on the amount of available data.
Validating our method with measured metabolite and protein levels in E. coli
central metabolism, we found typical prediction fold errors of 3.8 and 2.7,
respectively, for the two kinds of data. ECM can be used to predict enzyme
levels and protein cost in natural and engineered pathways, establishes a
direct connection between protein cost and thermodynamics, and provides a
physically plausible and computationally tractable way to include enzyme
kinetics into constraint-based metabolic models, where kinetics have usually
been ignored or oversimplified.


  Fluctuations in the physical properties of biological machines are
inextricably linked to their functions. Distributions of run-lengths and
velocities of processive molecular motors, like kinesin-1, are accessible
through single molecule techniques, yet there is lack a rigorous theoretical
model for these probabilities up to now. We derive exact analytic results for a
kinetic model to predict the resistive force ($F$) dependent velocity ($P(v)$)
and run-length ($P(n)$) distribution functions of generic finitely processive
molecular motors that take forward and backward steps on a track. Our theory
quantitatively explains the zero force kinesin-1 data for both $P(n)$ and
$P(v)$ using the detachment rate as the only parameter, thus allowing us to
obtain the variations of these quantities under load. At non-zero $F$, $P(v)$
is non-Gaussian, and is bimodal with peaks at positive and negative values of
$v$. The prediction that $P(v)$ is bimodal is a consequence of the discrete
step-size of kinesin-1, and remains even when the step-size distribution is
taken into account. Although the predictions are based on analyses of kinesin-1
data, our results are general and should hold for any processive motor, which
walks on a track by taking discrete steps.


  Although both RNA and proteins have densely packed native structures, chain
organizations of these two biopolymers are fundamentally different. Motivated
by the recent discoveries in chromatin folding that interphase chromosomes have
territorial organization with signatures pointing to metastability, we analyzed
the biomolecular structures deposited in the Protein Data Bank and found that
the intrachain contact probabilities, $P(s)$ as a function of the arc length
$s$, decay in power-law $\sim s^{-\gamma}$ over the intermediate range of $s$,
$10\lesssim s\lesssim 110$. We found that the contact probability scaling
exponent is $\gamma\approx 1.11$ for large RNA ($N>110$), $\gamma\approx 1.41$
for small sized RNA ($N<110$), and $\gamma\approx 1.65$ for proteins. Given
that Gaussian statistics is expected for a fully equilibrated chain in polymer
melts, the deviation of $\gamma$ value from $\gamma=1.5$ for the subchains of
large RNA in the native state suggests that the chain configuration of RNA is
not fully equilibrated. It is visually clear that folded structures of large
sized RNA ($N\gtrsim 110$) adopt crumpled structures, partitioned into modular
multi-domains assembled by proximal sequences along the chain, whereas the
polypeptide chain of folded proteins looks better mixed with the rest of the
structure. Our finding of $\gamma\approx 1$ for large RNA might be an
ineluctable consequence of the hierarchical ordering of the secondary to
tertiary elements in the folding process.


  To recognize pathogens, B and T lymphocytes are endowed with a wide
repertoire of receptors generated stochastically by V(D)J recombination.
Measuring and estimating the diversity of these receptors is of great
importance for understanding adaptive immunity. In this chapter we review
recent modeling approaches for analyzing receptor diversity from
high-throughput sequencing data. We first clarify the various existing notions
of diversity, with its many competing mathematical indices, and the different
biological levels at which it can be evaluated. We then describe inference
methods for characterizing the statistical diversity of receptors at different
stages of their history: generation, selection and somatic evolution. We
discuss the intrinsic difficulty of estimating the diversity of receptors
realized in a given individual from incomplete samples. Finally, we emphasize
the limitations of diversity defined at the level of receptor sequences, and
advocate the more relevant notion of functional diversity relative to the set
of recognized antigens.


  Mathematical models of population extinction have a variety of applications
in such areas as ecology, paleontology and conservation biology. Here we
propose and investigate two types of sub-exponential models of population
extinction.The results of performed analysis show that the principle of minimum
of information loss is the underlying law for the evolution of a broad class of
models of population extinction. Finally, we propose a possible application of
this modeling framework to mechanisms underlying time perception.


  Motor protein motion on biopolymers can be described by models related to the
totally asymmetric simple exclusion process (TASEP). Inspired by experiments on
the motion of kinesin-4 motors on antiparallel microtubule overlaps, we analyze
a model incorporating the TASEP on two antiparallel lanes with binding kinetics
and lane switching. We determine the steady-state motor density profiles using
phase plane analysis of the steady-state mean field equations and kinetic Monte
Carlo simulations. We focus on the the density-density phase plane, where we
find an analytic solution to the mean-field model. By studying the phase space
flows, we determine the model's fixed points and their changes with parameters.
Phases previously identified for the single-lane model occur for low switching
rate between lanes. We predict a new multiple coexistence phase due to
additional fixed points that appear as the switching rate increases: switching
moves motors from the higher-density to the lower-density lane, causing local
jamming and creating multiple domain walls. We determine the phase diagram of
the model for both symmetric and general boundary conditions.


  The amyloid $\beta$ peptide (A$\beta$42), whose aggregation is associated
with Alzheimer's disease, is an amphiphatic peptide with a high propensity to
self-assemble. A$\beta$42 has a net negative charge at physiological pH and
modulations of intermolecular electrostatic interactions can significantly
alter its aggregation behaviour. Variations in sequence and solution conditions
lead to varied macroscopic behaviour, often resulting in a number of different
mechanistic explanations for the aggregation of these closely related systems.
Here we alter the electrostatic interactions governing the fibril aggregation
kinetics by varying the ionic strength over an order of magnitude, which allows
us to sample the space of different reaction mechanisms, and develop a minimal
reaction network that explains the experimental kinetics under all the
different conditions. We find that an increase in the ionic strength leads to
an increased rate of surface catalysed nucleation over fragmentation and
eventually to a saturation of this nucleation process. More generally, this
reaction network connects previously separate systems, such as mutants of
A$\beta$42 and the wild type, on a continuous mechanistic landscape, thereby
providing a unified picture of the aggregation mechanism of A$\beta$42 and the
means of directly comparing the effects of intrinsic modifications of the
peptide to those of simple electrostatic shielding.


  A simplified model for the growth of a population is studied in which random
effects arise because reproducing individuals have a certain probability of
surviving until the next breeding season and hence contributing to the next
generation. The resulting Markov chain is that of a branching process with a
known generating function. For parameter values leading to non-extinction, an
approximating diffusion process is obtained for the population size.
  Results are obtained for the number of offspring $r_h$ and the initial
population size $N_0$ required to guarantee a given probabilty of survival. For
large probabilities of survival, increasing the initial population size from
$N_0=1$ to $N_0=2$ gives a very large decrease in required fecundity but
further increases in $N_0$ lead to much smaller decreases in $r_h$. For small
probabilities (< 0.2) of survival the decreases in required fecundity when
$N_0$ changes from 1 to 2 are very small. The calculations have relevance to
the survival of populations derived from colonizing individuals which could be
any of a variety of organisms.


  The Superior Colliculus, a laminar structure involved in the retinotopic
mapping of the visual field, plays a cardinal role in the several cortical and
subcortical loops of the saccadic system. Although the selection of saccade
targets has long been thought to be the sole product of cortical processes, a
growing body of evidence hints at the implication of the Superior Colliculus,
firstly by the lateral connections between the neurons of its maps, and
secondly by its interactions with the midbrain Basal Ganglia, already renowned
for their role in decision making. We propose a biomimetic population-coded
race model of selection based on a dynamic tecto-basal loop that reproduces the
observed ability of the Superior Colliculus to stochastically select between
similar stimuli, the accuracy of this selection depending on the
discriminability of the target and the distractors. Our model also offers an
explanation for the phenomenon of Remote Distractor Effect based on the lateral
connectivity within the Basal Ganglia circuitry rather than on lateral
inhibitions within the collicular maps. Finally, we propose a role for the
intermediate layers of the Superior Colliculus, as stochastic integrators
dynamically gated by the selective disinhibition of the Basal Ganglia channels
that is consistent with the recorded activity profiles of these neurons.


  Advances in science are being sought in newly available opportunities to
collect massive quantities of data about complex systems. While key advances
are being made in detailed mapping of systems, how to relate this data to
solving many of the challenges facing humanity is unclear. The questions we
often wish to address require identifying the impact of interventions on the
system and that impact is not apparent in the detailed data that is available.
Here we review key concepts and motivate a general framework for building
larger scale views of complex systems and for characterizing the importance of
information in physical, biological and social systems. We provide examples of
its application to evolutionary biology with relevance to ecology,
biodiversity, pandemics, and human lifespan, and in the context of social
systems with relevance to ethnic violence, global food prices, and stock market
panic. Framing scientific inquiry as an effort to determine what is important
and unimportant is a means for advancing our understanding and addressing many
practical concerns, such as economic development or treating disease.


  Network enrichment analysis is a powerful method, which allows to integrate
gene enrichment analysis with the information on relationships between genes
that is provided by gene networks. Existing tests for network enrichment
analysis deal only with undirected networks, they can be computationally slow
and are based on normality assumptions.
  We propose NEAT, a test for network enrichment analysis. The test is based on
the hypergeometric distribution, which naturally arises as the null
distribution in this context. NEAT can be applied not only to undirected, but
to directed and partially directed networks as well. Our simulations indicate
that NEAT is considerably faster than alternative resampling-based methods, and
that its capacity to detect enrichments is at least as good as the one of
alternative tests. We discuss applications of NEAT to network analyses in yeast
by testing for enrichment of the Environmental Stress Response target gene set
with GO Slim and KEGG functional gene sets, and also by testing for
associations between GO Slim categories themselves.
  NEAT is a flexible and efficient test for network enrichment analysis that
aims to overcome some limitations of existing resampling-based tests. The
method is implemented in the R package neat, which can be freely downloaded
from CRAN (http://cran.r-project.org/package=neat).


  The pathogenesis of cancer in human is still poorly understood. With the
rapid development of high-throughput sequencing technologies, huge volumes of
cancer genomics data have been generated. Deciphering those data poses great
opportunities and challenges to computational biologists. One of such key
challenges is to distinguish driver mutations, genes as well as pathways from
passenger ones. Mutual exclusivity of gene mutations (each patient has no more
than one mutation in the gene set) has been observed in various cancer types
and thus has been used as an important property of a driver gene set or
pathway. In this article, we aim to review the recent development of
computational models and algorithms for discovering driver pathways or modules
in cancer with the focus on mutual exclusivity-based ones.


  A new method to identify all sufficiently long repeating substrings in one or
several symbol sequences is proposed. The method is based on a specific gauge
applied to symbol sequences that guarantees identification of the repeating
substrings. It allows the matching of substrings to contain a given level of
errors. The gauge is based on the development of a heavily sparse dictionary of
repeats, thus drastically accelerating the search procedure. Some genomic
applications illustrate the method.
  This paper is the extended and detailed version of the presentation at the
third International Conference on Algorithms for Computational Biology to be
held at Trujillo, Spain, June 21-22, 2016.


  Till now, in biological sciences, the term, transcription, mainly refers to
DNA to RNA transcription. But our recently published experimental findings
obtained from Plasmodium falciparum strongly suggest the existence of DNA to
DNA transcription in the genome of eukaryotic cells, which could shed some
light on the functions of certain noncoding DNA in the human and other
eukaryotic genomes.


  This is an analysis of how magnetic fields affect biological molecules and
cells. It was prompted by a series of prominent reports regarding magnetism in
biological systems. The first claims to have identified a protein complex that
acts like a compass needle to guide magnetic orientation in animals (Qin et
al., 2016). Two other articles report magnetic control of membrane conductance
by attaching ferritin to an ion channel protein and then tugging the ferritin
or heating it with a magnetic field (Stanley et al., 2015; Wheeler et al.,
2016). Here I argue that these claims conflict with basic laws of physics. The
discrepancies are large: from 5 to 10 log units. If the reported phenomena do
in fact occur, they must have causes entirely different from the ones proposed
by the authors. The paramagnetic nature of protein complexes is found to
seriously limit their utility for engineering magnetically sensitive cells.


  Flow diagrams are a common tool used to help build and interpret models of
dynamical systems, often in biological contexts such as consumer-resource
models and similar compartmental models. Typically, their usage is intuitive
and informal. Here, we present a formalized version of flow diagrams as a kind
of weighted directed graph which follow a strict grammar, which translate into
a system of ordinary differential equations (ODEs) by a single unambiguous
rule, and which have an equivalent representation as a relational database. (We
abbreviate this schema of "ODEs and formalized flow diagrams" as OFFl.) Drawing
a diagram within this strict grammar encourages a mental discipline on the part
of the modeler in which all dynamical processes of a system are thought of as
interactions between dynamical species that draw parcels from one or more
source species and deposit them into target species according to a set of
transformation rules. From these rules, the net rate of change for each species
can be derived. The modeling schema can therefore be understood as both an
epistemic and practical heuristic for modeling, serving both as an
organizational framework for the model building process and as a mechanism for
deriving ODEs. All steps of the schema beyond the initial scientific
(intuitive, creative) abstraction of natural observations into model variables
are algorithmic and easily carried out by a computer, thus enabling the future
development of a dedicated software implementation. Such tools would empower
the modeler to consider significantly more complex models than practical
limitations might have otherwise proscribed, since the modeling framework
itself manages that complexity on the modeler's behalf. In this report, we
describe the chief motivations for OFFl, outline its implementation, and
utilize a range of classic examples from ecology and epidemiology to showcase
its features.


  Tumours develop in an evolutionary process, in which the accumulation of
mutations produces subpopulations of cells with distinct mutational profiles,
called clones. This process leads to the genetic heterogeneity widely observed
in tumour sequencing data, but identifying the genotypes and frequencies of the
different clones is still a major challenge. Here, we present Cloe, a
phylogenetic latent feature model to deconvolute tumour sequencing data into a
set of related genotypes. Our approach extends latent feature models by placing
the features as nodes in a latent tree. The resulting model can capture both
the acquisition and the loss of mutations, as well as episodes of convergent
evolution. We establish the validity of Cloe on synthetic data and assess its
performance on controlled biological data, comparing our reconstructions to
those of several published state-of-the-art methods. We show that our method
provides highly accurate reconstructions and identifies the number of clones,
their genotypes and frequencies even at a modest sequencing depth. As a proof
of concept we apply our model to clinical data from three cases with chronic
lymphocytic leukaemia, and one case with acute myeloid leukaemia.


  Network alignment (NA) aims to find a node mapping between molecular networks
of different species that identifies topologically or functionally similar
network regions. Analogous to genomic sequence alignment, NA can be used to
transfer biological knowledge from well- to poorly-studied species between
aligned network regions. Pairwise NA (PNA) finds similar regions between two
networks while multiple NA (MNA) can align more than two networks. We focus on
MNA. Existing MNA methods aim to maximize total similarity over all aligned
nodes (node conservation). Then, they evaluate alignment quality by measuring
the amount of conserved edges, but only after the alignment is constructed.
Directly optimizing edge conservation during alignment construction in addition
to node conservation may result in superior alignments. Thus, we present a
novel MNA approach called multiMAGNA++ that can achieve this. Indeed,
multiMAGNA++ generally outperforms or is on par with the existing MNA methods,
while often completing faster than the existing methods. That is, multiMAGNA++
scales well to larger network data and can be parallelized effectively. During
method evaluation, we also introduce new MNA quality measures to allow for more
complete alignment characterization as well as more fair MNA method comparison
compared to using only the existing alignment quality measures.


  Motivation: High throughput DNA sequencing (HTS) technologies generate an
excessive number of small DNA segments -- called short reads -- that cause
significant computational burden. To analyze the entire genome, each of the
billions of short reads must be mapped to a reference genome based on the
similarity between a read and "candidate" locations in that reference genome.
The similarity measurement, called alignment, formulated as an approximate
string matching problem, is the computational bottleneck because: (1) it is
implemented using quadratic-time dynamic programming algorithms, and (2) the
majority of candidate locations in the reference genome do not align with a
given read due to high dissimilarity. Calculating the alignment of such
incorrect candidate locations consumes an overwhelming majority of a modern
read mapper's execution time. Therefore, it is crucial to develop a fast and
effective filter that can detect incorrect candidate locations and eliminate
them before invoking computationally costly alignment operations. Results: We
propose GateKeeper, a new hardware accelerator that functions as a
pre-alignment step that quickly filters out most incorrect candidate locations.
GateKeeper is the first design to accelerate pre-alignment using
Field-Programmable Gate Arrays (FPGAs), which can perform pre-alignment much
faster than software. GateKeeper can be integrated with any mapper that
performs sequence alignment for verification. When implemented on a single FPGA
chip, GateKeeper maintains high accuracy (on average >96%) while providing up
to 90-fold and 130-fold speedup over the state-of-the-art software
pre-alignment techniques, Adjacency Filter and Shifted Hamming Distance (SHD),
respectively. The addition of GateKeeper as a pre-alignment step can reduce the
verification time of the mrFAST mapper by a factor of 10. Availability:
https://github.com/BilkentCompGen/GateKeeper


  Two atmospheric responses to simulated astrophysical ionizing radiation
events significant to life on Earth are production of odd-nitrogen species,
especially NO2, and subsequent depletion of stratospheric ozone. Ozone
depletion increases incident short-wavelength ultraviolet radiation (UVB,
280-315 nm) and longer ( > 600 nm) wavelengths of photosynthetically available
radiation (PAR, 400 -700 nm). On the other hand, the NO2 haze decreases
atmospheric transmission in the long-wavelength UVA (315-400 nm) and short
wavelength PAR. Here we use the results of previous simulations of incident
spectral irradiance following an ionizing radiation event to predict changes in
Terran productivity focusing on photosynthesis of marine phytoplankton. The
prediction is based on a spectral model of photosynthetic response developed
for the dominant genera in central regions of the ocean (Synechococcus and
Prochlorococcus), and remote-sensing based observations of spectral water
transparency, temperature, wind speed and mixed layer depth. Predicted
productivity declined after a simulated ionizing event, but the effect
integrated over the water column was small. For integrations taking into
account the full depth range of PAR transmission (down to 0.1% of utilizable
PAR), the decrease was at most 2-3% (depending on strain), with larger effects
(5-7%) for integrations just to the depth of the surface mixed layer. The
deeper integrations were most affected by the decreased utilizable PAR at depth
due to the NO2 haze, whereas shallower integrations were most affected by the
increased surface UV.


  We present a semi-quantitative model of condensation of short nucleic acid
(NA) duplexes induced by tri-valent cobalt(III) hexammine (CoHex) ions. The
model is based on partitioning of bound counterion distribution around singleNA
duplex into "external" and "internal" ion binding shells distinguished by the
proximity to duplex helical axis. In the aggregated phase the shells overlap,
which leads to significantly increased attraction of CoHex ions in these
overlaps with the neighboring duplexes. The duplex aggregation free energy is
decomposed into attractive and repulsive components in such a way that they can
be represented by simple analytical expressions with parameters derived from
molecular dynamic (MD) simulations and numerical solutions of Poisson equation.
The short-range interactions described by the attractive term depend on the
fractions of bound ions in the overlapping shells and affinity of CoHex to the
"external" shell of nearly neutralized duplex. The repulsive components of the
free energy are duplex configurational entropy loss upon the aggregation and
the electrostatic repulsion of the duplexes that remains after neutralization
by bound CoHex ions. The estimates of the aggregation free energy are
consistent with the experimental range of NA duplex condensation propensities,
including the unusually poor condensation of RNA structures and subtle sequence
effects upon DNA condensation. The model predicts that, in contrast to DNA, RNA
duplexes may condense into tighter packed aggregates with a higher degree of
duplex neutralization. The model also predicts that longer NA fragments will
condense more readily than shorter ones. The ability of this model to explain
experimentally observed trends in NA condensation, lends support to proposed NA
condensation picture based on the multivalent "ion binding shells".


  Pathogens can spread epidemically through populations. Beneficial contagions,
such as viruses that enhance host survival or technological innovations that
improve quality of life, also have the potential to spread epidemically. How do
the dynamics of beneficial biological and social epidemics differ from those of
detrimental epidemics? We investigate this question using three theoretical
approaches. First, in the context of population genetics, we show that a
horizontally-transmissible element that increases fitness, such as viral DNA,
spreads superexponentially through a population, more quickly than a beneficial
mutation. Second, in the context of behavioral epidemiology, we show that
infections that cause increased connectivity lead to superexponential fixation
in the population. Third, in the context of dynamic social networks, we find
that preferences for increased global infection accelerate spread and produce
superexponential fixation, but preferences for local assortativity halt
epidemics by disconnecting the infected from the susceptible. We conclude that
the dynamics of beneficial biological and social epidemics are characterized by
the rapid spread of beneficial elements, which is facilitated in biological
systems by horizontal transmission and in social systems by active spreading
behavior of infected individuals.


  Motivation: The analysis of antibodies and T-cell receptors (TCRs)
concentrations in serum is a fundamental problem in immunoinformatics.
Repertoire construction is a preliminary step of analysis of clonal lineages,
understanding of immune response dynamics, population analysis of
immunoglobulin and TCR loci. Emergence of MiSeq Illumina sequencing machine in
2013 opened horizons of investigation of adaptive immune repertoires using
highly accurate reads. Reads produced by MiSeq are able to cover repertoires of
moderate size. At the same time, throughput of sequencing machines increases
from year to year. This will enable ultra deep scanning of adaptive immune
repertoires and analysis of their diversity. Such data requires both efficient
and highly accurate repertoire construction tools. In 2015 Safonova et al.
presented IgRepertoireConstructor, a tool for accurate construction of antibody
repertoire and immunoproteogenomics analysis. Unfortunately, proposed algorithm
was very time and memory consuming and could be a bottleneck of processing
large immunosequencing libraries. In this paper we overcome this challenge and
present IgReC, a novel algorithm for adaptive repertoire construction problem.
IgReC reconstructs a repertoire with high precision even if each input read
contains sequencing errors and performs well on contemporary datasets. Results
of computational experiments show that IgReC improves state-of-the-art in the
field. Availability: IgReC is an open source and freely available program
running on Linux platforms. The source code is available at GitHub:
yana-safonova.github.io/ig_repertoire_constructor. Contact:
safonova.yana@gmail.com


  Tumorigenic transformation of human epithelial cells in vitro has been
described experimentally as the potential result of a process known as
spontaneous immortalization. In this process a generic series of cell state
transitions occur in which normal epithelial cells acquire a senescent state,
later surpassed to attain first a mesenchymal state and then a final
mesenchymal stem-like phenotype, with potential tumorigenic behavior. In this
paper we integrate published data on the molecular components and interactions
that have been described as key regulators of such cell states and transitions.
A large network, that is provided, is constructed and then reduced with the aim
of recovering a minimal regulatory core incorporating the necessary and
sufficient restrictions to recover the observed cell states and their generic
progression patterns in epithelial-mesenchymal transition. Data is formalized
into logical regulatory rules that govern the dynamics of each of the network's
components as a function of the states of its regulators. The proposed core
gene regulatory network attains only three steady-state gene expression
configurations that correspond to the profiles characteristic of normal
epithelial, senescent, and mesenchymal stem-like cells. Interestingly,
epigenetic landscape analyses of the uncovered network shows that it also
recovers the generic time-ordered transitions documented during tumorigenic
transformation in vitro of epithelial cells. The latter strongly correlate with
the patterns observed during the progressive pathological description of
epithelial carcinogenesis in vivo.


  New long read sequencing technologies, like PacBio SMRT and Oxford NanoPore,
can produce sequencing reads up to 50,000 bp long but with an error rate of at
least 15%. Reducing the error rate is necessary for subsequent utilisation of
the reads in, e.g., de novo genome assembly. The error correction problem has
been tackled either by aligning the long reads against each other or by a
hybrid approach that uses the more accurate short reads produced by second
generation sequencing technologies to correct the long reads. We present an
error correction method that uses long reads only. The method consists of two
phases: first we use an iterative alignment-free correction method based on de
Bruijn graphs with increasing length of k-mers, and second, the corrected reads
are further polished using long-distance dependencies that are found using
multiple alignments. According to our experiments the proposed method is the
most accurate one relying on long reads only for read sets with high coverage.
Furthermore, when the coverage of the read set is at least 75x, the throughput
of the new method is at least 20% higher. LoRMA is freely available at
http://www.cs.helsinki.fi/u/lmsalmel/LoRMA/.


  High throughput mRNA expression profiling can be used to characterize the
response of cell culture models to perturbations such as pharmacologic
modulators and genetic perturbations. As profiling campaigns expand in scope,
it is important to homogenize, summarize, and analyze the resulting data in a
manner that captures significant biological signals in spite of various noise
sources such as batch effects and stochastic variation. We used the L1000
platform for large-scale profiling of 978 genes, chosen to be representative of
the genome as whole, across thousands of compound treatments. Here, a method is
described that uses deep learning techniques to convert the expression changes
of the landmark genes into a perturbation barcode that reveals important
features of the underlying data, performing better than the raw data in
revealing important biological insights. The barcode captures compound
structure and target information, in addition to predicting a compound's high
throughput screening promiscuity, to a higher degree than the original data
measurements, indicating that the approach uncovers underlying factors of the
expression data that are otherwise entangled or masked by noise. Furthermore,
we demonstrate that visualizations derived from the perturbation barcode can be
used to more sensitively assign functions to unknown compounds through a
guilt-by-association approach, which we use to predict and experimentally
validate the activity of compounds on the MAPK pathway. The demonstrated
application of deep metric learning to large-scale chemical genetics projects
highlights the utility of this and related approaches to the extraction of
insights and testable hypotheses from big, sometimes noisy data.


  Background. Large scale metagenomic projects aim to extract biodiversity
knowledge between different environmental conditions. Current methods for
comparing microbial communities face important limitations. Those based on
taxonomical or functional assignation rely on a small subset of the sequences
that can be associated to known organisms. On the other hand, de novo methods,
that compare the whole sets of sequences, either do not scale up on ambitious
metagenomic projects or do not provide precise and exhaustive results.
  Methods. These limitations motivated the development of a new de novo
metagenomic comparative method, called Simka. This method computes a large
collection of standard ecological distances by replacing species counts by
k-mer counts. Simka scales-up today's metagenomic projects thanks to a new
parallel k-mer counting strategy on multiple datasets.
  Results. Experiments on public Human Microbiome Project datasets demonstrate
that Simka captures the essential underlying biological structure. Simka was
able to compute in a few hours both qualitative and quantitative ecological
distances on hundreds of metagenomic samples (690 samples, 32 billions of
reads). We also demonstrate that analyzing metagenomes at the k-mer level is
highly correlated with extremely precise de novo comparison techniques which
rely on all-versus-all sequences alignment strategy or which are based on
taxonomic profiling.


  Ecological networks are theoretical abstractions that represent ecological
communities. These networks are usually defined as static entities, in which
the occurrence of a particular interaction between species is considered fixed
despite the intrinsic dynamics of ecological systems. However, empirical
analysis of the temporal variation of trophic interactions is constrained by
the lack of data with high spatial, temporal, and taxonomic resolution. Here,
we evaluate the spatiotemporal variability of multiple consumer-resource
interactions of large marine networks. The tropic interactions of all of the
analyzed networks had low temporal persistence, which was well described by a
common exponential decay in the rank-frequency relationship of
consumer-resource interactions. This common pattern of low temporal persistence
was evident despite the dissimilarities of environmental conditions among
sites. Between-site rank correlations of frequency of occurrence of
interactions ranged from 0.59 to 0.73. After removing the interactions with
<50% frequency, the between-site correlations decreased to values between 0.60
and 0.28, indicating that low-frequency interactions accounted for the apparent
similarities between sites. Our results showed that the communities studied
were characterized by few persistent interactions and a large number of
transient trophic interactions. We suggest that consumer-resource temporal
asynchrony in addition to varying local environmental conditions and
opportunistic foraging could be among the mechanisms generating the observed
rank-frequency relationship of trophic interactions. Therefore, our results
question the analysis of ecological communities as static and persistent
natural entities and stress the need for strengthening the analysis of temporal
variability in ecological networks and long-term studies.


  Motivation: As cancer researchers have come to appreciate the importance of
intratumor heterogeneity, much attention has focused on the challenges of
accurately profiling heterogeneity in individual patients. Experimental
technologies for directly profiling genomes of single cells are rapidly
improving, but they are still impractical for large-scale sampling. Bulk
genomic assays remain the standard for population-scale studies, but conflate
the influences of mixtures of genetically distinct tumor, stromal, and
infiltrating immune cells. Many computational approaches have been developed to
deconvolute these mixed samples and reconstruct the genomics of genetically
homogeneous clonal subpopulations. All such methods, however, are limited to
reconstructing only coarse approximations to a few major subpopulations. In
prior work, we showed that one can improve deconvolution of genomic data by
leveraging substructure in cellular mixtures through a strategy called
simplicial complex inference. This strategy, however, is also limited by the
difficulty of inferring mixture structure from sparse, noisy assays. Results:
We improve on past work by introducing enhancements to automate learning of
substructured genomic mixtures, with specific emphasis on genome-wide copy
number variation (CNV) data. We introduce methods for dimensionality estimation
to better decompose mixture model substructure; fuzzy clustering to better
identify substructure in sparse, noisy data; and automated model inference
methods for other key model parameters. We show that these improvements lead to
more accurate inference of cell populations and mixture proportions in
simulated scenarios. We further demonstrate their effectiveness in identifying
mixture substructure in real tumor CNV data. Availability: Source code is
available at http://www.cs.cmu.edu/~russells/software/WSCUnmix.zip


  The advent of next-generation sequencing (NGS) technologies enables
researchers to sequence complex microbial communities directly from
environment. Since assembly typically produces only genome fragments, also
known as contigs, instead of entire genome, it is crucial to group them into
operational taxonomic units (OTUs) for further taxonomic profiling and
down-streaming functional analysis. OTU clustering is also referred to as
binning. We present COCACOLA, a general framework automatically bin contigs
into OTUs based upon sequence composition and coverage across multiple samples.
  The effectiveness of COCACOLA is demonstrated in both simulated and real
datasets in comparison to state-of-art binning approaches such as CONCOCT,
GroopM, MaxBin and MetaBAT. The superior performance of COCACOLA relies on two
aspects. One is employing $L_{1}$ distance instead of Euclidean distance for
better taxonomic identification during initialization. More importantly,
COCACOLA takes advantage of both hard clustering and soft clustering by
sparsity regularization.
  In addition, the COCACOLA framework seamlessly embraces customized knowledge
to facilitate binning accuracy. In our study, we have investigated two types of
additional knowledge, the co-alignment to reference genomes and linkage of
contigs provided by paired-end reads, as well as the ensemble of both. We find
that both co-alignment and linkage information further improve binning in the
majority of cases. COCACOLA is scalable and faster than CONCOCT ,GroopM, MaxBin
and MetaBAT.
  The software is available at https://github.com/younglululu/COCACOLA


  A kinetic theory is developed for exonuclease-deficient DNA polymerases,
based on the experimental observation that the rates depend not only on the
newly incorporated nucleotide, but also on the previous one, leading to the
growth of Markovian DNA sequences from a Bernoullian template. The dependences
on nucleotide concentrations and template sequence are explicitly taken into
account. In this framework, the kinetic and thermodynamic properties of DNA
replication, in particular, the mean growth velocity, the error probability,
and the entropy production in terms of the rate constants and the
concentrations are calculated analytically. Theory is compared with numerical
simulations for the DNA polymerases of T7 viruses and human mitochondria.


  Kinetic theory and thermodynamics are applied to DNA polymerases with
exonuclease activity, taking into account the dependence of the rates on the
previously incorportated nucleotide. The replication fidelity is shown to
increase significantly thanks to this dependence at the basis of the mechanism
of exonuclease proofreading. In particular, this dependence can provide up to a
hundred-fold lowering of the error probability under physiological conditions.
Theory is compared with numerical simulations for the DNA polymerases of T7
viruses and human mitochondria.


  The reconstruction of phylogenetic trees from mixed populations has become
important in the study of cancer evolution, as sequencing is often performed on
bulk tumor tissue containing mixed populations of cells. Recent work has shown
how to reconstruct a perfect phylogeny tree from samples that contain mixtures
of two-state characters, where each character/locus is either mutated or not.
However, most cancers contain more complex mutations, such as copy-number
aberrations, that exhibit more than two states. We formulate the Multi-State
Perfect Phylogeny Mixture Deconvolution Problem of reconstructing a multi-state
perfect phylogeny tree given mixtures of the leaves of the tree. We
characterize the solutions of this problem as a restricted class of spanning
trees in a graph constructed from the input data, and prove that the problem is
NP-complete. We derive an algorithm to enumerate such trees in the important
special case of cladisitic characters, where the ordering of the states of each
character is given. We apply our algorithm to simulated data and to two cancer
datasets. On simulated data, we find that for a small number of samples, the
Multi-State Perfect Phylogeny Mixture Deconvolution Problem often has many
solutions, but that this ambiguity declines quickly as the number of samples
increases. On real data, we recover copy-neutral loss of heterozygosity,
single-copy amplification and single-copy deletion events, as well as their
interactions with single-nucleotide variants.


  Many systems in biology, physics and engineering can be described by systems
of ordinary differential equation containing many parameters. When studying the
dynamic behavior of these large, nonlinear systems, it is useful to identify
and characterize the steady-state solutions as the model parameters vary, a
technically challenging problem in a high-dimensional parameter landscape.
Rather than simply determining the number and stability of steady-states at
distinct points in parameter space, we decompose the parameter space into
finitely many regions, the steady-state solutions being consistent within each
distinct region. From a computational algebraic viewpoint, the boundary of
these regions is contained in the discriminant locus. We develop global and
local numerical algorithms for constructing the discriminant locus and
classifying the parameter landscape. We showcase our numerical approaches by
applying them to molecular and cell-network models.


  In epidemiological modelling, dynamics on networks, and in particular
adaptive and heterogeneous networks have recently received much interest. Here
we present a detailed analysis of a previously proposed model that combines
heterogeneity in the individuals with adaptive rewiring of the network
structure in response to a disease. We show that in this model qualitative
changes in the dynamics occur in two phase transitions. In a macroscopic
description one of these corresponds to a local bifurcation whereas the other
one corresponds to a non-local heteroclinic bifurcation. This model thus
provides a rare example of a system where a phase transition is caused by a
non-local bifurcation, while both micro- and macro-level dynamics are
accessible to mathematical analysis. The bifurcation points mark the onset of a
behaviour that we call network inoculation. In the respective parameter region
exposure of the system to a pathogen will lead to an outbreak that collapses,
but leaves the network in a configuration where the disease cannot reinvade,
despite every agent returning to the susceptible class. We argue that this
behaviour and the associated phase transitions can be expected to occur in a
wide class of models of sufficient complexity.


  Metagenomic binning is an essential task in analyzing metagenomic sequence
datasets. To analyze structure or function of microbial communities from
environmental samples, metagenomic sequence fragments are assigned to their
taxonomic origins. Although sequence alignment algorithms can readily be used
and usually provide high-resolution alignments and accurate binning results,
the computational cost of such alignment-based methods becomes prohibitive as
metagenomic datasets continue to grow. Alternative compositional-based methods,
which exploit sequence composition by profiling local short k-mers in
fragments, are often faster but less accurate than alignment-based methods.
Inspired by the success of linear error correcting codes in noisy channel
communication, we introduce Opal, a fast and accurate novel compositional-based
binning method. It incorporates ideas from Gallager's low-density parity-check
code to design a family of compact and discriminative locality-sensitive
hashing functions that encode long-range compositional dependencies in long
fragments. By incorporating the Gallager LSH functions as features in a simple
linear SVM, Opal provides fast, accurate and robust binning for datasets
consisting of a large number of species, even with mutations and sequencing
errors. Opal not only performs up to two orders of magnitude faster than BWA,
an alignment-based binning method, but also achieves improved binning accuracy
and robustness to sequencing errors. Opal also outperforms models built on
traditional k-mer profiles in terms of robustness and accuracy. Finally, we
demonstrate that we can effectively use Opal in the "coarse search" stage of a
compressive genomics pipeline to identify a much smaller candidate set of
taxonomic origins for a subsequent alignment-based method to analyze, thus
providing metagenomic binning with high scalability, high accuracy and high
resolution.


  Metagenomics has revealed hundreds of bacterial species in almost all
microbiota. In a few well-studied cases, bacterial communities have been
observed to coordinate their metabolic fluxes. In principle, bacteria can
divide tasks to reap the benefits of specialization, as in human economies.
However, the benefits and stability of an economy of bacterial specialists are
far from obvious. Here, we physically model the population dynamics of bacteria
that compete for steadily supplied resources. Importantly, we explicitly model
the metabolic fluxes yielding cellular biomass production under the constraint
of a limited enzyme budget. In our framework, we find that population dynamics
generally leads to the coexistence of different metabolic types, which satisfy
an extended competitive exclusion principle (even allowing for adaptive
mutation). We establish that these consortia act as cartels, whereby population
dynamics pins down resource concentrations at values for which no other
strategy can invade. Finally, we propose that at steady supply, cartels of
competing strategies automatically yield maximum biomass, thereby achieving a
collective optimum.


  The intrinsic stochasticity of gene expression is usually mitigated in higher
eukaryotes by post-transcriptional regulation channels that stabilise the
output layer, most notably protein levels. The discovery of small non-coding
RNAs (miRNAs) in specific motifs of the genetic regulatory network has led to
identifying noise buffering as the possible key function they exert in
regulation. Recent in vitro} and in silico studies have corroborated this
hypothesis. It is however also known that miRNA-mediated noise reduction is
hampered by transcriptional bursting in simple topologies. Here, using
stochastic simulations validated by analytical calculations based on van
Kampen's expansion, we revisit the noise-buffering capacity of the
miRNA-mediated Incoherent Feed Forward Loop (IFFL), a small module that is
widespread in the gene regulatory networks of higher eukaryotes, in order to
account for the effects of intermittency in the transcriptional activity of the
modulator gene. We show that bursting considerably alters the circuit's ability
to control static protein noise. By comparing with other regulatory
architectures, we find that direct transcriptional regulation significantly
outperforms the IFFL in a broad range of kinetic parameters. This suggests
that, under pulsatile inputs, static noise reduction may be less important than
dynamical aspects of noise and information processing in characterising the
performance of regulatory elements.


  Three-dimensional interphase organization of metazoan genomes has been linked
to cellular identity. However, the principles governing 3D interphase genome
architecture and its faithful transmission through disruptive events of
cell-cycle, like mitosis, are not fully understood. By using Brownian dynamics
simulations of Drosophila chromosome 3R up to time-scales of minutes, we show
that chromatin binding profile of Polycomb-repressive-complex-1 robustly
predicts a sub-set of topologically associated domains (TADs), and inclusion of
other factors recapitulates the profile of all TADs, as observed
experimentally. Our simulations show that chromosome 3R attains interphase
organization from mitotic state by a two-step process in which formation of
local TADs is followed by long-range interactions. Our model also explains
statistical features and tracks the assembly kinetics of polycomb subnuclear
clusters. In conclusion, our approach can be used to predict structural and
kinetic features of 3D chromosome folding and its associated proteins in
biological relevant genomic and time scales.


  While metagenomics has emerged as a technology of choice for analyzing
bacterial populations, assembly of metagenomic data remains difficult thus
stifling biological discoveries. metaSPAdes is a new assembler that addresses
the challenge of metagenome analysis and capitalizes on computational ideas
that proved to be useful in assemblies of single cells and highly polymorphic
diploid genomes. We benchmark metaSPAdes against other state-of-the-art
metagenome assemblers across diverse da-tasets and demonstrate that it results
in high-quality assemblies.


  Efficient text indexing data structures have enabled large-scale genomic
sequence analysis and are used to help solve problems ranging from assembly to
read mapping. However, these data structures typically assume that the
underlying reference text is static and will not change over the course of the
queries being made. Some progress has been made in exploring how certain text
indices, like the suffix array, may be updated, rather than rebuilt from
scratch, when the underlying reference changes. Yet, these update operations
can be complex in practice, difficult to implement, and give fairly pessimistic
worst-case bounds. We present a novel data structure, SkipPatch, for
maintaining a k-mer-based index over a dynamically changing genome. SkipPatch
pairs a hash-based k-mer index with an indexable skip list that is used to
efficiently maintain the set of edits that have been applied to the original
genome. SkipPatch is practically fast, significantly outperforming the dynamic
extended suffix array in terms of update and query speed.


  Protein-protein interaction (PPI) networks, providing a comprehensive
landscape of protein interacting patterns, enable us to explore biological
processes and cellular components at multiple resolutions. For a biological
process, a number of proteins need to work together to perform the job.
Proteins densely interact with each other, forming large molecular machines or
cellular building blocks. Identification of such densely interconnected
clusters or protein complexes from PPI networks enables us to obtain a better
understanding of the hierarchy and organization of biological processes and
cellular components. Most existing methods apply efficient graph clustering
algorithms on PPI networks, often failing to detect possible densely connected
subgraphs and overlapped subgraphs. Besides clustering-based methods, dense
subgraph enumeration methods have also been used, which aim to find all densely
connected protein sets. However, such methods are not practically tractable
even on a small yeast PPI network, due to high computational complexity. In
this paper, we introduce a novel approximate algorithm to efficiently enumerate
putative protein complexes from biological networks. The key insight of our
algorithm is that we do not need to enumerate all dense subgraphs. Instead we
only need to find a small subset of subgraphs that cover as many proteins as
possible. The problem is formulated as finding a diverse set of dense
subgraphs, where we develop highly effective pruning techniques to guarantee
efficiency. To handle large networks, we take a divide-and-conquer approach to
speed up the algorithm in a distributed manner. By comparing with existing
clustering and dense subgraph-based algorithms on several human and yeast PPI
networks, we demonstrate that our method can detect more putative protein
complexes and achieves better prediction accuracy.


  Motivation: De novo transcriptome assembly of non-model organisms is the
first major step for many RNA-seq analysis tasks. Current methods for de novo
assembly often report a large number of contiguous sequences (contigs), which
may be fractured and incomplete sequences instead of full-length transcripts.
Dealing with a large number of such contigs can slow and complicate downstream
analysis.
  Results :We present a method for clustering contigs from de novo
transcriptome assemblies based upon the relationships exposed by multi-mapping
sequencing fragments. Specifically, we cast the problem of clustering contigs
as one of clustering a sparse graph that is induced by equivalence classes of
fragments that map to subsets of the transcriptome. Leveraging recent
developments in efficient read mapping and transcript quantification, we have
developed RapClust, a tool implementing this approach that is capable of
accurately clustering most large de novo transcriptomes in a matter of minutes,
while simultaneously providing accurate estimates of expression for the
resulting clusters. We compare RapClust against a number of tools commonly used
for de novo transcriptome clustering. Using de novo assemblies of organisms for
which reference genomes are available, we assess the accuracy of these
different methods in terms of the quality of the resulting clusterings, and the
concordance of differential expression tests with those based on ground truth
clusters. We find that RapClust produces clusters of comparable or better
quality than existing state-of-the-art approaches, and does so substantially
faster. RapClust also confers a large benefit in terms of space usage, as it
produces only succinct intermediate files - usually on the order of a few
megabytes - even when processing hundreds of millions of reads.


  Recent research in both the experimental and mathematical communities has
focused on biochemical interaction systems that satisfy an "absolute
concentration robustness" (ACR) property. The ACR property was first discovered
experimentally when, in a number of different systems, the concentrations of
key system components at equilibrium were observed to be robust to the total
concentration levels of the system. Followup mathematical work focused on
deterministic models of biochemical systems and demonstrated how chemical
reaction network theory can be utilized to explain this robustness. Later
mathematical work focused on the behavior of this same class of reaction
networks, though under the assumption that the dynamics were stochastic. Under
the stochastic assumption, it was proven that the system will undergo an
extinction event with a probability of one so long as the system is
conservative, showing starkly different long-time behavior than in the
deterministic setting. Here we consider a general class of stochastic models
that intersects with the class of ACR systems studied previously. We consider a
specific system scaling over compact time intervals and prove that in a limit
of this scaling the distribution of the abundances of the ACR species converges
to a certain product-form Poisson distribution whose mean is the ACR value of
the deterministic model. This result is in agreement with recent conjectures
pertaining to the behavior of ACR networks endowed with stochastic kinetics,
and helps to resolve the conflicting theoretical results pertaining to
deterministic and stochastic models in this setting.


  Pentameric ligand-gated ion channels (pLGICs) of the Cys-loop superfamily are
important neuroreceptors that mediate fast synaptic transmission. They are
activated by the binding of a neurotransmitter, but the details of this process
are still not fully understood. As a prototypical pLGIC, here we choose the
insect resistance to dieldrin (RDL) receptor, involved in the resistance to
insecticides, and investigate the binding of the neurotransmitter GABA to its
extracellular domain at the atomistic level. We achieve this by means of
$\mu$-sec funnel-metadynamics simulations, which efficiently enhance the
sampling of bound and unbound states by using a funnel-shaped restraining
potential to limit the exploration in the solvent. We reveal the sequence of
events in the binding process, from the capture of GABA from the solvent to its
pinning between the charged residues Arg111 and Glu204 in the binding pocket.
We characterize the associated free energy landscapes in the wild-type RDL
receptor and in two mutant forms, where the key residues Arg111 and Glu204 are
mutated to Ala. Experimentally these mutations produce non-functional channels,
which is reflected in the reduced ligand binding affinities, due to the loss of
essential interactions. We also analyze the dynamical behaviour of the crucial
loop C, whose opening allows the access of GABA to the binding site, while its
closure locks the ligand into the protein. The RDL receptor shares structural
and functional features with other pLGICs, hence our work outlines a valuable
protocol to study the binding of ligands to pLGICs beyond conventional docking
and molecular dynamics techniques.


  Motivation: High-throughput experimental techniques have been producing more
and more protein-protein interaction (PPI) data. PPI network alignment greatly
benefits the understanding of evolutionary relationship among species, helps
identify conserved sub-networks and provides extra information for functional
annotations. Although a few methods have been developed for multiple PPI
network alignment, the alignment quality is still far away from perfect and
thus, new network alignment methods are needed. Result: In this paper, we
present a novel method, denoted as ConvexAlign, for joint alignment of multiple
PPI networks by convex optimization of a scoring function composed of sequence
similarity, topological score and interaction conservation score. In contrast
to existing methods that generate multiple alignments in a greedy or
progressive manner, our convex method optimizes alignments globally and
enforces consistency among all pairwise alignments, resulting in much better
alignment quality. Tested on both synthetic and real data, our experimental
results show that ConvexAlign outperforms several popular methods in producing
functionally coherent alignments. ConvexAlign even has a larger advantage over
the others in aligning real PPI networks. ConvexAlign also finds a few
conserved complexes among 5 species which cannot be detected by the other
methods.


  Cyclic adenosine monophosphate (cAMP) is considered a model system for signal
transduction, the mechanism by which cells exchange chemical messages. Our
previous work calculated the Shannon capacity of a single cAMP receptor;
however, a typical cell may have thousands of receptors operating in parallel.
In this paper, we calculate the capacity of a cAMP signal transduction system
with an arbitrary number of independent, indistinguishable receptors. By
leveraging prior results on feedback capacity for a single receptor, we show
(somewhat unexpectedly) that the capacity is achieved by an IID input
distribution, and that the capacity for n receptors is n times the capacity for
a single receptor.


  In cell proliferation, stem cell differentiation, chemoresistance and tissue
organization, the ubiquitous role of YAP/TAZ continues to impact our
fundamental understanding in numerous physiological and disease systems.
YAP/TAZ is an important signaling nexus integrating diverse mechanical and
biochemical signals, such as ECM stiffness, adhesion ligand density, or
cell-cell contacts, and thus strongly influences cell fate. Recent studies show
that YAP/TAZ mechanical sensing is dependent on RhoA-regulated stress fibers.
However, current understanding of YAP/TAZ still remains limited due to the
unknown interaction between the canonical Hippo pathway and cell tension. To
identify the roles of key signaling molecules in mechanical signal sensing and
transduction, we present a novel computational model of the YAP/TAZ signaling
pathway. This model converts ECM mechanical properties to biochemical signals
via adhesion, and integrates intracellular signaling cascades associated with
cytoskeleton dynamics. Adhesion molecules, such as FAK, are predicted to rescue
YAP/TAZ activity in soft environments via the RhoA pathway. We found that
changes of molecule concentrations result in different pattern of YAP/TAZ
stiffness response. We also investigate the sensitivity of YAP/TAZ activity to
ECM stiffness. In addition, the model shows that the unresolved synergistic
effect of YAP/TAZ activity between the mechanosensing and the Hippo pathways
can be explained by the interaction of LIMK and LATS. Overall, our model
provides a novel platform for studying YAP/TAZ activity in the context of
integrating different signaling pathways. This platform can be used to gain new
fundamental insights into roles of key molecular and mechanical regulators on
development, tissue engineering or tumor progression.


  A central concern of network operators is to estimate the probability of an
incident that affects a significant part and thus may yield to a breakdown. We
answer this question by modeling how a failure of either a node or an edge will
affect the rest of the network using percolation theory. Our model is general
in the sense that it only needs two inputs: the topology of the network and the
chances of failure of its components. These chances may vary to represent
different types of edges having different tendencies to fail. We illustrate the
approach by an example, for which we can even obtain closed form expressions
for the likelihood of an outbreak remaining bounded or spreading unlimitedly.


  The string graph for a collection of next-generation reads is a lossless data
representation that is fundamental for de novo assemblers based on the
overlap-layout-consensus paradigm. In this paper, we explore a novel approach
to compute the string graph, based on the FM-index and Burrows-Wheeler
Transform. We describe a simple algorithm that uses only the FM-index
representation of the collection of reads to construct the string graph,
without accessing the input reads. Our algorithm has been integrated into the
SGA assembler as a standalone module to construct the string graph.
  The new integrated assembler has been assessed on a standard benchmark,
showing that FSG is significantly faster than SGA while maintaining a moderate
use of main memory, and showing practical advantages in running FSG on multiple
threads.


  Circadian (~24hr) clocks are self-sustained endogenous oscillators with which
organisms keep track of daily and seasonal time. Circadian clocks frequently
rely on interlocked transcriptional- translational feedback loops to generate
rhythms that are robust against intrinsic and extrinsic perturbations. To
investigate the dynamics and mechanisms of the intracellular feedback loops in
circadian clocks, a number of mathematical models have been developed. The
majority of the models use Hill functions to describe transcriptional
repression in a way that is similar to the Goodwin model. Recently, a new class
of models with protein sequestration-based repression has been introduced.
Here, we discuss how this new class of models differs dramatically from those
based on Hill-type repression in several fundamental aspects: conditions for
rhythm generation, robust network designs and the periods of coupled
oscillators. Consistently, these fundamental properties of circadian clocks
also differ among Neurospora, Drosophila, and mammals depending on their key
transcriptional repression mechanisms (Hill-type repression or protein
sequestration). Based on both theoretical and experimental studies, this review
highlights the importance of careful modeling of transcriptional repression
mechanisms in molecular circadian clocks.


  We study a continuous time model for the frequency distribution of an
infinitely large asexual population in which both beneficial and deleterious
mutations occur and the fitness is additive. When beneficial mutations are
ignored, the exact solution for the frequency distribution is known to be a
Poisson distribution. Here we include beneficial mutations and obtain exact
expressions for the frequency distribution at all times using an eigenfunction
expansion method. We find that the stationary distribution is non-Poissonian
and related to the Bessel function of the first kind. We also provide suitable
approximations for the stationary distribution and the time to relax to the
steady state. Our exact results, especially at mutation-selection equilibrium,
can be useful in developing semi-deterministic approaches to understand
stochastic evolution.


  We show that some natural output conventions for error-free computation in
chemical reaction networks (CRN) lead to a common level of computational
expressivity. Our main results are that the standard consensus-based output
convention have equivalent computational power to (1) existence-based and (2)
democracy-based output conventions. The CRNs using the former output convention
have only "yes" voters, with the interpretation that the CRN's output is yes if
any voters are present and no otherwise. The CRNs using the latter output
convention define output by majority vote among "yes" and "no" voters.
  Both results are proven via a generalized framework that simultaneously
captures several definitions, directly inspired by a Petri net result of
Esparza, Ganty, Leroux, and Majumder [CONCUR 2015]. These results support the
thesis that the computational expressivity of error-free CRNs is intrinsic, not
sensitive to arbitrary definitional choices.


  In the face of serious infectious diseases, governments endeavour to
implement containment measures such as public vaccination at a macroscopic
level. Meanwhile, individuals tend to protect themselves by avoiding contacts
with infections at a microscopic level. However, a comprehensive understanding
of how such combined strategy influences epidemic dynamics is still lacking. We
study a susceptible-infected-susceptible epidemic model with imperfect
vaccination on dynamic contact networks, where the macroscopic intervention is
represented by random vaccination of the population and the microscopic
protection is characterised by susceptible individuals rewiring contacts from
infective neighbours. In particular, the model is formulated both in
populations without and then with demographic effects. Using the pairwise
approximation and the probability generating function approach, we investigate
both dynamics of the epidemic and the underlying network. For populations
without demography, the emerging degree correlations, bistable states, and
oscillations demonstrate the combined effects of the public vaccination program
and individual protective behavior. Compared to either strategy in isolation,
the combination of public vaccination and individual protection is more
effective in preventing and controlling the spread of infectious diseases by
increasing both the invasion threshold and the persistence threshold. For
populations with additional demographic factors, the integration between
vaccination intervention and individual rewiring may promote epidemic spreading
due to the birth effect. Moreover, the degree distributions of both networks in
the steady state is closely related to the degree distribution of newborns,
which leads to uncorrelated connectivity. All the results demonstrate the
importance of both local protection and global intervention, as well as the
demographic effects.


  The influence of hypertension on low-density lipoproteins intake into the
arterial wall is an important factor for understanding mechanisms of
atherosclerosis. It has been experimentally observed that the increased
pressure leads to the higher level of the LDL inside the wall. In this paper we
attempt to construct a model of the LDL transport which reproduces
quantitatively experimental outcomes. We supplement the well known four-layer
arterial wall model to include two pressure induced effects: the compression of
the intima tissue and the increase of the fraction of leaky junctions. We
demonstrate that such model can reach the very good agreement with experimental
data.


  The Zika virus is the primary suspect in the large increase in microcephaly
cases in 2015-6 in Brazil, however its role is unconfirmed despite individual
cases of viral infections found in neural tissue. Here we consider the
alternative that the insecticide pyriproxyfen, used in Brazilian drinking water
for mosquito control, may actually be the cause. Pyriproxifen is an analog of
juvenile hormone, which corresponds in mammals to regulatory molecules
including retinoic acid, a vitamin A metabolite, with which it has
cross-reactivity and whose application during development causes microcephaly.
Methoprene, another juvenile hormone analog approved as an insecticide has
metabolites that bind to the retinoid X receptor, and causes developmental
disorders in mammals. Isotretinoin is another example of a retinoid causing
microcephaly in human babies via activation of the retinoid X receptor.
Moreover, tests of pyriproxyfen by the manufacturer, Sumitomo, widely quoted as
giving no evidence for developmental toxicity, actually found some evidence for
such an effect, including low brain mass and arhinencephaly--incomplete
formation of the anterior cerebral hemispheres--in rat pups. Finally, the
pyriproxyfen use in Brazil is unprecedented--it has never before been applied
to a water supply on such a scale. Claims that it is not being used in Recife,
the epicenter of microcephaly cases, do not distinguish the metropolitan area
of Recife, where it is widely used, and the municipality, where it is not.
Given this combination of information we strongly recommend that the use of
pyriproxyfen in Brazil be suspended pending further investigation.


  We develop a continuum theory for equilibrium elasticity of a network of
crosslinked semiflexible filaments, spanning the full range between flexible
entropy-driven chains to stiff athermal rods. We choose the 3-chain
constitutive model over several plausible candidates, and derive analytical
expressions for the elastic energy at arbitrary strain, with the corresponding
stress-strain relationship. The theory fits well to a wide range of
experimental data on simple shear in different filament networks,
quantitatively matching the differential shear modulus variation with stress,
with only two adjustable parameters (which represent the filament stiffness and
pre-tension in the network, respectively). The general theory accurately
describes the crossover between the positive and negative Poynting effect
(normal stress on imposed shear) on increasing the stiffness of filaments
forming the network.


  The ongoing Zika virus (ZIKV) epidemic poses a major global public health
emergency. It is known that ZIKV is spread by \textit{Aedes} mosquitoes, recent
studies show that ZIKV can also be transmitted via sexual contact and cases of
sexually transmitted ZIKV have been confirmed in the U.S., France, and Italy.
How sexual transmission affects the spread and control of ZIKV infection is not
well-understood. We presented a mathematical model to investigate the impact of
mosquito-borne and sexual transmission on spread and control of ZIKV and used
the model to fit the ZIKV data in Brazil, Colombia, and El Salvador. Based on
the estimated parameter values, we calculated the median and confidence
interval of the basic reproduction number R0=2.055 (95% CI: 0.523-6.300), in
which the distribution of the percentage of contribution by sexual transmission
is 3.044 (95% CI: 0.123-45.73). Our study indicates that R0 is most sensitive
to the biting rate and mortality rate of mosquitoes while sexual transmission
increases the risk of infection and epidemic size and prolongs the outbreak. In
order to prevent and control the transmission of ZIKV, it must be treated as
not only a mosquito-borne disease but also a sexually transmitted disease.


  Duality plays an important role in population genetics. It can relate results
from forwards-in-time models of allele frequency evolution with those of
backwards-in-time genealogical models; a well known example is the duality
between the Wright-Fisher diffusion for genetic drift and its genealogical
counterpart, the coalescent. There have been a number of articles extending
this relationship to include other evolutionary processes such as mutation and
selection, but little has been explored for models also incorporating crossover
recombination. Here, we derive from first principles a new genealogical process
which is dual to a Wright-Fisher diffusion model of drift, mutation, and
recombination. Our approach is based on expressing a putative duality
relationship between two models via their infinitesimal generators, and then
seeking an appropriate test function to ensure the validity of the duality
equation. This approach is quite general, and we use it to find dualities for
several important variants, including both a discrete L-locus model of a gene
and a continuous model in which mutation and recombination events are scattered
along the gene according to continuous distributions. As an application of our
results, we derive a series expansion for the transition function of the
diffusion. Finally, we study in further detail the case in which mutation is
absent. Then the dual process describes the dispersal of ancestral genetic
material across the ancestors of a sample. The stationary distribution of this
process is of particular interest; we show how duality relates this
distribution to haplotype fixation probabilities. We develop an efficient
method for computing such probabilities in multilocus models.


  Dietrich and Haider (2014) justify their integrative framework for creativity
founded on evolutionary theory and prediction research on the grounds that
"theories and approaches guiding empirical research on creativity have not been
supported by the neuroimaging evidence". Although this justification is
controversial, the general direction holds promise. This commentary clarifies
points of disagreement and unresolved issues, and addresses mis-applications of
evolutionary theory that lead the authors to adopt a Darwinian (versus
Lamarckian) approach. To say that creativity is Darwinian is not to say that it
consists of variation plus selection---in the everyday sense of the term---as
the authors imply; it is to say that evolution is occurring because selection
is affecting the distribution of randomly generated heritable variation across
generations. In creative thought the distribution of variants is not key, i.e.,
one is not inclined toward idea A because 60% of one's candidate ideas are
variants of A while only 40% are variants of B; one is inclined toward
whichever seems best. The authors concede that creative variation is partly
directed; however, the greater the extent to which variants are generated
non-randomly, the greater the extent to which the distribution of variants can
reflect not selection but the initial generation bias. Since each thought in a
creative process can alter the selective criteria against which the next is
evaluated, there is no demarcation into generations as assumed in a Darwinian
model. We address the authors' claim that reduced variability and individuality
are more characteristic of Lamarckism than Darwinian evolution, and note that a
Lamarckian approach to creativity has addressed the challenge of modeling the
emergent features associated with insight.


  When evolutionary games are contested in structured populations, the degree
of each player in the network plays an important role. If they exist, hubs
often determine the fate of the population in remarkable ways. Recent research
based on optimal percolation in random networks has shown, however, that the
degree is neither the sole nor the best predictor of influence in complex
networks. Low-degree nodes may also be optimal influencers if they are
hierarchically linked to hubs. Taking this into account leads to the formalism
of collective influence in complex networks, which as we show here, has
far-reaching implications for the favorable resolution of social dilemmas. In
particular, there exists an optimal hierarchical depth for the determination of
collective influence that we use to describe the potency of players for passing
their strategies, which depends on the strength of the social dilemma.
Interestingly, the degree, which corresponds to the baseline depth zero, is
optimal only when the temptation to defect is small. Our research reveals that
evolutionary success stories are related to spreading processes which are
rooted in favorable hierarchical structures that extend beyond local
neighborhoods.


  The detection of rare variants is important for understanding the genetic
heterogeneity in mixed samples. Recently, next-generation sequencing (NGS)
technologies have enabled the identification of single nucleotide variants
(SNVs) in mixed samples with high resolution. Yet, the noise inherent in the
biological processes involved in next-generation sequencing necessitates the
use of statistical methods to identify true rare variants. We propose a novel
Bayesian statistical model and a variational expectation-maximization (EM)
algorithm to estimate non-reference allele frequency (NRAF) and identify SNVs
in heterogeneous cell populations. We demonstrate that our variational EM
algorithm has comparable sensitivity and specificity compared with a Markov
Chain Monte Carlo (MCMC) sampling inference algorithm, and is more
computationally efficient on tests of low coverage ($27\times$ and $298\times$)
data. Furthermore, we show that our model with a variational EM inference
algorithm has higher specificity than many state-of-the-art algorithms. In an
analysis of a directed evolution longitudinal yeast data set, we are able to
identify a time-series trend in non-reference allele frequency and detect novel
variants that have not yet been reported. Our model also detects the emergence
of a beneficial variant earlier than was previously shown, and a pair of
concomitant variants.


  Cell membranes are out of thermodynamic equilibrium notably because of
membrane recycling, i.e. active exchange of material with the cytosol. We
propose an analytically tractable model of biomembrane predicting the effects
of recycling on the size of protein nanodomains. It includes a short-range
attraction between proteins and a weaker long-range repulsion which ensures the
existence of so-called cluster phases at equilibrium, where monomeric proteins
coexist with finite-size domains. Our main finding is that when taking
recycling into account, the typical cluster size increases logarithmically with
the recycling rate. Using physically realistic model parameters, the predicted
two-fold increase due to recycling in living cells is very likely
experimentally measurable with the help of super-resolution microscopy.


  We studied the structuredness in a chloroplast genome of Siberian larch. The
clusters in 63-dimensional space were identified with elastic map technique,
where the objects to be clusterized are the different fragments of the genome.
A seven-cluster structure in the distribution of those fragments reported
previously has been found. Unlike the previous results, we have found the
drastically other composition of the clusters comprising the fragments
extracted from coding and non-coding regions of the genome.


  Contacts between individuals play an important role in determining how
infectious diseases spread. Various methods to gather data on such contacts
co-exist, from surveys to wearable sensors. Comparisons of data obtained by
different methods in the same context are however scarce, in particular with
respect to their use in data-driven models of spreading processes. Here, we use
a combined data set describing contacts registered by sensors and friendship
relations in the same population to address this issue in a case study. We
investigate if the use of the friendship network is equivalent to a sampling
procedure performed on the sensor contact network with respect to the outcome
of simulations of spreading processes: such an equivalence might indeed give
hints on ways to compensate for the incompleteness of contact data deduced from
surveys. We show that this is indeed the case for these data, for a
specifically designed sampling procedure, in which respondents report their
neighbors with a probability depending on their contact time. We study the
impact of this specific sampling procedure on several data sets, discuss
limitations of our approach and its possible applications in the use of data
sets of various origins in data-driven simulations of epidemic processes.


  Considering agricultural landscapes as networks can provide information about
spatial connectivity relevant for a wide range of applications including
pollination, pest management, and ecology. Global agricultural networks are
well-described by power law rank-size distributions. However, regional analyses
capture only a subset of the total global network. Most analyses are regional.
In this paper, we seek to address the following questions: Does the globally
observed scale-free property of agricultural networks hold over smaller spatial
domains? Can similar properties be observed at kilometer to meter scales? We
analyze 9 intensively cultivated Landsat scenes on 5 continents with a wide
range of vegetation distributions. We find that networks of vegetation fraction
within the domain of each of these Landsat scenes exhibit substantial
variability - but still possess similar scaling properties to the global
distribution of agriculture. We also find similar results using a 39 km2 IKONOS
image. To illustrate an application of spatial network analysis, we show an
example of network disruption. We compare two networks with similar rank-size
distributions that behave differently when nodes are progressively removed. We
suggest that treating agricultural land cover as spatial networks can provide a
straightforward way of characterizing the connectivity of complex spatial
distributions of agriculture across a wide range of landscapes and at spatial
scales relevant for practical agricultural applications.


  When controlling an emerging outbreak of an infectious disease it is
essential to know the key epidemiological parameters, such as the basic
reproduction number $R_0$ and the control effort required to prevent a large
outbreak. These parameters are estimated from the observed incidence of new
cases and information about the infectious contact structures of the population
in which the disease spreads. However, the relevant infectious contact
structures for new, emerging infections are often unknown or hard to obtain.
Here we show that for many common true underlying heterogeneous contact
structures, the simplification to neglect such structures and instead assume
that all contacts are made homogeneously in the whole population, results in
conservative estimates for $R_0$ and the required control effort. This means
that robust control policies can be planned during the early stages of an
outbreak, using such conservative estimates of the required control effort.


  Recently, a long-term model of HIV infection dynamics was developed to
describe the entire time course of the disease. It consists of a large system
of ODEs with many parameters, and is expensive to simulate. In the current
paper, this model is analyzed by determining all infection-free steady states
and studying the local stability properties of the unique biologically-relevant
equilibrium. Active subspace methods are then used to perform a global
sensitivity analysis and study the dependence of an infected individual's
T-cell count on the parameter space. Building on these results, a
global-in-time approximation of the T-cell count is created by constructing
dynamic active subspaces and reduced order models are generated, thereby
allowing for inexpensive computation.


  Existing mathematical models of the shoot apical meristem (SAM) explain
nucleation and confinement of a stem cell domain by Turing's mechanism,
assuming that the diffusion coefficients of the activator (WUSCHEL) and
inhibitor (CLAVATA) are significantly different. As there is no evidence for
this assumption of differential diffusivity, we recently proposed a new
mechanism based on a bistable switch model of the SAM. Here we study the
bistable-switch mechanism in detail, demonstrating that it can be understood as
localized switches of WUSHEL activity in individual cells driven by a
non-uniform field of a peptide hormone. By comparing domain formation by Turing
and bistable-switch mechanisms on a cell network, we show that the latter does
not require the assumptions needed by the former, which are not supported by
biological evidences.


  A simple model to predict the directionality of vibrational energy flow at
molecular level is presented. This model is based on a vibrational energy
propagation analysis using ab intio molecular dynamics and the Fukui function
and local softness reactivity indexes derived from DFT. By using this simple
conceptual model we are giving a cogent rationale to previous theoretical and
experimental reports of asymmetrical vibrational energy diffusion in proteins
and energetic materials. We proposed here three basic rules for the Vibrational
Energy Relaxation in molecules: (i) the vibrational energy flows form a soft
site to a hard site but not in the opposite direction, (ii) when vibrational
energy is injected directly to a hard site, it is trapped in this site, and
(iii) if the vibrational energy is pump in a polarizable site and two sites
with different softness are available, the energy will propagate to the softest
one.


  We propose a frame-based representation of k-mers for detecting sequencing
errors and rare variants in next generation sequencing data obtained from
populations of closely related genomes. Frames are sets of non-orthogonal basis
functions, traditionally used in signal processing for noise removal. We define
a frame for genomes and sequenced reads to consist of discrete spatial signals
of every k-mer of a given size. We show that each k-mer in the sequenced data
can be projected onto multiple frames and these projections are maximized for
spatial signals corresponding to the k-mer's substrings. Our proposed
classifier, MultiRes, is trained on the projections of k-mers as features used
for marking k-mers as erroneous or true variations in the genome. We evaluate
MultiRes on simulated and real viral population datasets and compare it to
other error correction methods known in the literature. MultiRes has 4 to 500
times less false positives k-mer predictions compared to other methods,
essential for accurate estimation of viral population diversity and their
de-novo assembly. It has high recall of the true k-mers, comparable to other
error correction methods. MultiRes also has greater than 95% recall for
detecting single nucleotide polymorphisms (SNPs), fewer false positive SNPs,
while detecting higher number of rare variants compared to other variant
calling methods for viral populations. The software is freely available from
the GitHub link (https://github.com/raunaq-m/MultiRes).


  Size varies. Small things are typically more frequent than large things. The
logarithm of frequency often declines linearly with the logarithm of size. That
power law relation forms one of the common patterns of nature. Why does the
complexity of nature reduce to such a simple pattern? Why do things as
different as tree size and enzyme rate follow similarly simple patterns? Here I
analyze such patterns by their invariant properties. For example, a common
pattern should not change when adding a constant value to all observations.
That shift is essentially the renumbering of the points on a ruler without
changing the metric information provided by the ruler. A ruler is shift
invariant only when its scale is properly calibrated to the pattern being
measured. Stretch invariance corresponds to the conservation of the total
amount of something, such as the total biomass and consequently the average
size. Rotational invariance corresponds to pattern that does not depend on the
order in which underlying processes occur, for example, a scale that additively
combines the component processes leading to observed values. I use tree size as
an example to illustrate how the key invariances shape pattern. A simple
interpretation of common pattern follows. That simple interpretation connects
the normal distribution to a wide variety of other common patterns through the
transformations of scale set by the fundamental invariances.


  Over the past decade, several targeted therapies (e.g. imatinib, dasatinib,
nilotinib) have been developed to treat Chronic Myeloid Leukemia (CML). Despite
an initial response to therapy, drug resistance remains a problem for some CML
patients. Recent studies have shown that resistance mutations that preexist
treatment can be detected in a substan- tial number of patients, and that this
may be associated with eventual treatment failure. One proposed method to
extend treatment efficacy is to use a combination of multiple targeted
therapies. However, the design of such combination therapies (timing, sequence,
etc.) remains an open challenge. In this work we mathematically model the
dynamics of CML response to combination therapy and analyze the impact of
combination treatment schedules on treatment efficacy in patients with
preexisting resistance. We then propose an optimization problem to find the
best schedule of multiple therapies based on the evolution of CML according to
our ordinary differential equation model. This resulting optimiza- tion problem
is nontrivial due to the presence of ordinary different equation constraints
and integer variables. Our model also incorporates realistic drug toxicity
constraints by tracking the dynamics of patient neutrophil counts in response
to therapy. Using realis- tic parameter estimates, we determine optimal
combination strategies that maximize time until treatment failure.


  The mutant allele frequencies in oncogenes peak around 0.40 and rapidly
decrease. In this article, we explain why this is the case. Invoking a key
result from mathematical analysis in our model, namely, the inverse function
theorem, we estimate the selection pressures of the mutations as a function of
germline allele frequencies. Under complete dominance of oncogenic mutations,
this selection function is expected to be linearly correlated with the
distribution of the mutant alleles. We demonstrate that this is the case by
investigating the allele frequencies of mutations in oncogenes across various
cancer types, validating our model for mean effective selection. Consistent
with the population genetics model of fitness, the selection function fits a
gamma distribution curve that accurately describes the trend of the mutant
allele frequencies. While existing equations for selection explain evolution at
low allele frequencies, our equations are general formulas for natural
selection under complete dominance operating at all frequencies. We show that
selection exhibits linear behavior at all times, favoring dominant alleles with
respect to the change in recessive allele frequency. Also, these equations
show, selection behaves like power-law against the recessive alleles at low
dominant allele frequency.


  Deterministically growing (wild-type) populations which seed stochastically
developing mutant clones have found an expanding number of applications from
microbial populations to cancer. The special case of exponential wild-type
population growth, usually termed the Luria-Delbr\"uck or Lea-Coulson model, is
often assumed but seldom realistic. In this article we generalise this model to
different types of wild-type population growth, with mutants evolving as a
birth-death branching process. Our focus is on the size distribution of clones
- that is the number of progeny of a founder mutant - which can be mapped to
the total number of mutants. Exact expressions are derived for exponential,
power-law and logistic population growth. Additionally for a large class of
population growth we prove that the long time limit of the clone size
distribution has a general two-parameter form, whose tail decays as a
power-law. Considering metastases in cancer as the mutant clones, upon
analysing a data-set of their size distribution, we indeed find that a
power-law tail is more likely than an exponential one.


  The $\text{NJ}_{st}$ method was proposed by Liu and Yu to infer a species
tree topology from unrooted topological gene trees. While its statistical
consistency under the multispecies coalescent model was established only for a
4-taxon tree, simulations demonstrated its good performance on gene trees
inferred from sequences for many taxa. Here we prove the statistical
consistency of the method for an arbitrarily large species tree. Our approach
connects $\text{NJ}_{st}$ to a generalization of the STAR method of Liu, Pearl
and Edwards, and a previous theoretical analysis of it. We further show
$\text{NJ}_{st}$ utilizes only the distribution of splits in the gene trees,
and not their individual topologies. Finally, we discuss how multiple samples
per taxon per gene should be handled for statistical consistency.


  In this paper, we focus on model reduction of biomolecular systems with
multiple time-scales, modeled using the Linear Noise Approximation. Considering
systems where the Linear Noise Approximation can be written in singular
perturbation form, with $\epsilon$ as the singular perturbation parameter, we
obtain a reduced order model that approximates the slow variable dynamics of
the original system. In particular, we show that, on a finite time-interval,
the first and second moments of the reduced system are within an
$O(\epsilon)$-neighborhood of the first and second moments of the slow variable
dynamics of the original system. The approach is illustrated on an example of a
biomolecular system that exhibits time-scale separation.


  We propose a physically-realisable biochemical device that is coupled to a
biochemical reservoir of mutual information, fuel molecules and a chemical
bath. Mutual information allows work to be done on the bath even when the fuel
molecules appear to be in equilibrium; alternatively, mutual information can be
created by driving from the fuel or the bath. The system exhibits diverse
behaviour, including a regime in which the information, despite increasing
during the reaction, enhances the extracted work. We further demonstrate that a
modified device can function without the need for external manipulation,
eliminating the need for a complex and potentially costly control.