  The famous "two-fold cost of sex" is really the cost of anisogamy -- why
should females mate with males who do not contribute resources to offspring,
rather than isogamous partners who contribute equally? In typical anisogamous
populations, a single very fit male can have an enormous number of offspring,
far larger than is possible for any female or isogamous individual. If the
sexual selection on males aligns with the natural selection on females,
anisogamy thus allows much more rapid adaptation via super-successful males. We
show via simulations that this effect can be sufficient to overcome the
two-fold cost and maintain anisogamy against isogamy in populations adapting to
environmental change. The key quantity is the variance in male fitness -- if
this exceeds what is possible in an isogamous population, anisogamous
populations can win out in direct competition by adapting faster.


  We study the effect of cholesterol on the structure of
dipalmitoylphosphatidylcholine (DPPC) phospholipid bilayers. Using extensive
molecular dynamics computer simulations at atomistic resolution we observe and
quantify several structural changes upon increasing cholesterol content that
are collectively known as the condensing effect: a thickening of the bilayer,
an increase in lipid tail order, and a decrease in lateral area. We also
observe a change in leaflet interdigitation and a lack thereof in the
distributions of DPPC head group orientations. These results, obtained over a
wide range of cholesterol mole fractions, are then used to calibrate the
analysis of phospholipid properties in bilayers containing a single cholesterol
molecule per leaflet, which we perform in a spatially resolved way. We find
that a single cholesterol molecule affects phospholipids in its first and
second solvation shells, which puts the range of this interaction to be on the
order of one to two nanometers. We also observe a tendency of phospholipids to
orient their polar head groups toward the cholesterol, which provides
additional support for the umbrella model of bilayer organization.


  Traveling fronts describe the transition between two alternative states in a
great number of physical and biological systems. Examples include the spread of
beneficial mutations, chemical reactions, and the invasions by foreign species.
In homogeneous environments, the alternative states are separated by a smooth
front moving at a constant velocity. This simple picture can break down in
structured environments such as tissues, patchy landscapes, and microfluidic
devices. Habitat fragmentation can pin the front at a particular location or
lock invasion velocities into specific values. Locked velocities are not
sensitive to moderate changes in dispersal or growth and are determined by the
spatial and temporal periodicity of the environment. The synchronization with
the environment results in discontinuous fronts that propagate as periodic
pulses. We characterize the transition from continuous to locked invasions and
show that it is controlled by positive density-dependence in dispersal or
growth. We also demonstrate that velocity locking is robust to demographic and
environmental fluctuations and examine stochastic dynamics and evolution in
locked invasions.


  In this paper we present ACEMod, an agent-based modelling framework for
studying influenza epidemics in Australia. The simulator is designed to analyse
the spatiotemporal spread of contagion and influenza spatial synchrony across
the nation. The individual-based epidemiological model accounts for mobility
(worker and student commuting) patterns and human interactions derived from the
2006 Australian census and other national data sources. The high-precision
simulation comprises 19.8 million stochastically generated software agents and
traces the dynamics of influenza viral infection and transmission at several
scales. Using this approach, we are able to synthesise epidemics in Australia
with varying outbreak locations and severity. For each scenario, we investigate
the spatiotemporal profiles of these epidemics, both qualitatively and
quantitatively, via incidence curves, prevalence choropleths, and epidemic
synchrony. This analysis exemplifies the nature of influenza pandemics within
Australia and facilitates future planning of effective intervention, mitigation
and crisis management strategies.


  Single-cell experiments have revealed cell-to-cell variability in generation
times and growth rates for genetically identical cells. Theoretical models
relating the fluctuating generation times of single cells to the population
growth rate are usually based on the assumption that the generation times of
mother and daughter cells are uncorrelated. This assumption, however, is
inconsistent with the exponential growth of cell volume in time observed for
many cell types. Here we develop a more general and biologically relevant model
in which cells grow exponentially and generation times are correlated in a
manner which controls cell size. In addition to the fluctuating generation
times, we also allow the single-cell growth rates to fluctuate and account for
their correlations across the lineage tree. Surprisingly, we find that the
population growth rate only depends on the distribution of single-cell growth
rates and their correlations.


  In the present study, different photoperiods and nutritional conditions were
applied to a mixed wastewater-borne cyanobacterial culture in order to enhance
the intracellular accumulation of polyhydroxybutyrates (PHBs) and
carbohydrates. Two different experimental set-ups were used. In the first, the
culture was permanently exposed to illumination, while in the second; it was
submitted to light/dark alternation (12h cycles). In both cases, two different
nutritional regimes were also evaluated, N-limitation and P-limitation. Results
showed that the highest PHB concentration (104 mg L-1) was achieved under P
limited conditions and permanent illumination, whereas the highest carbohydrate
concentration (838 mg L-1) was obtained under N limited condition and
light/dark alternation. With regard to bioplastics and biofuel generation, this
study demonstrates that the accumulation of PHBs (bioplastics) and
carbohydrates (potential biofuel substrate) is favored in wastewater-borne
cyanobacteria under conditions where nutrients are limited.


  Direct cDNA preamplification protocols developed for single-cell RNA-seq have
enabled transcriptome profiling of precious clinical samples and rare cells
without sample pooling or RNA extraction. Currently, there is no algorithm
optimized to reveal and remove noisy transcripts in limiting-cell RNA-seq
(lcRNA-seq) data for downstream analyses. Herein, we present CLEAR, a workflow
that identifies reliably quantifiable transcripts in lcRNA-seq data for
differentially expressed gene (DEG) analysis. Libraries at three input amounts
of FACS-derived CD5+ and CD5- cells from a chronic lymphocytic leukemia patient
were used to develop CLEAR. When using CLEAR transcripts vs. using all
transcripts, downstream analyses revealed more shared transcripts across
different input RNA amounts, improved Principal Component Analysis (PCA)
separation, and yielded more DEGs between cell types. As proof-of-principle,
CLEAR was applied to an in-house lcRNA-seq dataset and two public datasets.
When imputation is used, CLEAR is also adaptable to large clinical studies and
for single cell analyses.


  We study the evolutionary dynamics of games under environmental feedback
using replicator equations for two interacting populations. One key feature is
to consider jointly the co-evolution of the dynamic payoff matrices and the
state of the environment: the payoff matrix varies with the changing
environment and at the same time, the state of the environment is affected
indirectly by the changing payoff matrix through the evolving population
profiles. For such co-evolutionary dynamics, we investigate whether convergence
will take place, and if so, how. In particular, we identify the scenarios where
oscillation offers the best predictions of long-run behavior by using
reversible system theory. The obtained results are useful to describe the
evolution of multi-community societies in which individuals' payoffs and
societal feedback interact.


  Public health services are constantly searching for new ways to reduce the
spread of infectious diseases, such as public vaccination of asymptomatic
individuals, quarantine (isolation) and treatment of symptomatic individuals.
Epidemic models have a long history of assisting in public health planning and
policy making. In this paper, we introduce epidemic models including variable
population size, degree-related imperfect vaccination and quarantine on
scale-free networks. More specifically, the models are formulated both on the
population with and without permanent natural immunity to infection, which
corresponds respectively to the
susceptible-vaccinated-infected-quarantined-recovered (SVIQR) model and the
susceptible-vaccinated-infected-quarantined (SVIQS) model. We develop different
mathematical methods and techniques to study the dynamics of two models,
including the basic reproduction number, the global stability of disease-free
and endemic equilibria. For the SVIQR model, we show that the system exhibits a
forward bifurcation. Meanwhile, the disease-free and unique endemic equilibria
are shown to be globally asymptotically stable by constructing suitable
Lyapunov functions. For the SVIQS model, conditions ensuring the occurrence of
multiple endemic equilibria are derived. Under certain conditions, this system
cannot undergo a backward bifurcation. The global asymptotical stability of
disease-free equilibrium, and the persistence of the disease are proved. The
endemic equilibrium is shown to be globally attractive by using monotone
iterative technique. Finally, stochastic network simulations yield quantitative
agreement with the deterministic mean-field approach.


  In this paper we analyze the length-spectrum of rainbows in RNA secondary
structures. A rainbow in a secondary structure is a maximal arc with respect to
the partial order induced by nesting. We show that there is a significant gap
in this length-spectrum. We shall prove that there asymptotically almost surely
exists a unique longest rainbow of length at least $n-O(n^{1/2})$ and that with
high probability any other rainbow has finite length. We show that the
distribution of the length of the longest rainbow converges to a discrete limit
law and that, for finite $k$, the distribution of rainbows of length $k$,
becomes for large $n$ a negative binomial distribution. We then put the results
of this paper into context, comparing the analytical results with those
observed in RNA minimum free energy structures, biological RNA structures and
relate our findings to the sparsification of folding algorithms.


  Although many phenotypic traits are determined by a large number of genetic
variants, how a polygenic trait adapts in response to the changes in the
environment is still poorly understood. Here we study the adaptation dynamics
of a polygenic trait that is determined by a finite number of genetic loci in
an infinitely large population which is evolving under stabilising selection
and recurrent mutations. We find that in a changing environment, modeled here
by a linearly moving phenotypic optimum, the mean trait also moves linearly
with time. But its speed is smaller than that of the phenotypic optimum when
the effect sizes of the genetic variants are small and approaches that of the
environmental change for larger effect sizes. Our study thus highlights the
influence of the genetic architecture of a polygenic trait on its adaptability.


  We review the standard lock and key (LK) model for binding small ligands to
larger adsorbent molecule. We discuss three levels of the traditional LK model
for binding. Within this model the binding constant or the Gibbs energy of the
binding process is related to the total interaction energy between the ligand
and the binding site of the adsorbent molecules. When solvent molecules are
present, which is the case in all binding processes in biochemistry, we find
that a major part of the Gibbs energy of binding could be due to interactions
mediated through the solvent molecules. This finding could have major
consequences to the applicability of the LK model in drug design, and perhaps
require a shift in the prevailing paradigm in this field of research.
  Keywords: Lock and key model, Binding constant, solvent effect on binding,
hydrophilic effect, molecular recognition.


  Tree-grass coexistence in savanna ecosystems depends strongly on
environmental disturbances out of which crucial is fire. Most modeling attempts
in the literature lack stochastic approach to fire occurrences which is
essential to reflect their unpredictability. Existing models that actually
include stochasticity of fire are usually analyzed only numerically. We
introduce new minimalistic model of tree-grass coexistence where fires occur
according to stochastic process. We use the tools of linear semigroup theory to
provide more careful mathematical analysis of the model. Essentially we show
that there exists a unique stationary distribution of tree and grass biomasses.


  Motivation: Despite its great success in various physical modeling,
differential geometry (DG) has rarely been devised as a versatile tool for
analyzing large, diverse and complex molecular and biomolecular datasets due to
the limited understanding of its potential power in dimensionality reduction
and its ability to encode essential chemical and biological information in
differentiable manifolds. Results: We put forward a differential geometry based
geometric learning (DG-GL) hypothesis that the intrinsic physics of
three-dimensional (3D) molecular structures lies on a family of low-dimensional
manifolds embedded in a high-dimensional data space. We encode crucial
chemical, physical and biological information into 2D element interactive
manifolds, extracted from a high-dimensional structural data space via a
multiscale discrete-to-continuum mapping using differentiable density
estimators. Differential geometry apparatuses are utilized to construct element
interactive curvatures %Gaussian curvature, mean curvature, minimum curvature
and maximum curvature in analytical forms for certain analytically
differentiable density estimators. These low-dimensional differential geometry
representations are paired with a robust machine learning algorithm to showcase
their descriptive and predictive powers for large, diverse and complex
molecular and biomolecular datasets. Extensive numerical experiments are
carried out to demonstrated that the proposed DG-GL strategy outperforms other
advanced methods in the predictions of drug discovery related protein-ligand
binding affinity, drug toxicity, and molecular solvation free energy.


  Stochastic epidemic models, generally more realistic than deterministic
counterparts, have often been seen too complex for rigorous mathematical
analysis because of level of details it requires to comprehensively capture the
dynamics of diseases. This problem further becomes intense when complexity of
diseasees increases as in the case of vector-borne diseases (VBD). The VBDs are
human illnesses caused by pathogens transmitted among humans by intermediate
species, which are primarily arthropods. In this study, a stochastic VBD model
is developed and novel mathematical methods are described and evaluated to
systematically analyze the model and understand its complex dynamics. The VBD
model incorporates some relevant features of the VBD transmission process
including demographical, ecological and social mechanisms. The analysis is
based on dimensional reductions and model simplications via scaling limit
theorems. The results suggest that the dynamics of the stochastic VBD depends
on a threshold quantity R_0, the initial size of infectives, and the type of
scaling in terms of host population size. The quantity R_0 for deterministic
counterpart of the model, interpreted as threshold condition for infection
persistence as is mentioned in the literature for many infectious disease
models, can be computed. Different scalings yield different approximations of
the model, and in particular, if vectors have much faster dynamics, the effect
of the vector dynamics on the host population averages out, which largely
reduces the dimension of the model.


  The present study investigates different strategies for the treatment of a
mixture of digestate from an anaerobic digester diluted and secondary effluent
from a high rate algal pond. To this aim, the performance of two
photo-sequencing batch reactors (PSBRs) operated at high nutrients loading
rates and different solids retention times (SRTs) were compared with a
semi-continuous photobioreactor (SC). Performances were evaluated in terms of
wastewater treatment, biomass composition and biopolymers accumulation during
30 days of operation. PSBRs were operated at a hydraulic retention time (HRT)
of 2 days and SRTs of 10 and 5 days (PSBR2-10 and PSBR2-5, respectively),
whereas the semi-continuous reactor was operated at a coupled HRT/SRT of 10
days (SC10-10). Results showed that PSBR2-5 achieved the highest removal rates
in terms of TN (6.7 mg L-1 d-1), TP (0.31 mg L-1 d-1), TOC (29.32 mg L-1 d-1)
and TIC (3.91 mg L-1 d-1). These results were in general 3-6 times higher than
the removal rates obtained in the SC10-10 (TN 29.74 mg L-1 d-1, TP 0.96 mg L-1
d-1, TOC 29.32 mg L-1 d-1 and TIC 3.91 mg L-1 d-1). Furthermore, both PSBRs
were able to produce biomass up to 0.09 g L-1 d-1, more than twofold the
biomass produced by the semi-continuous reactor (0.04 g L-1 d-1), and achieved
a biomass settleability of 86-92%. This study also demonstrated that the
microbial composition could be controlled by the nutrients loads, since the
three reactors were dominated by different species depending on the nutritional
conditions. Concerning biopolymers accumulation, carbohydrates concentration
achieved similar values in the three reactors (11%), whereas <0.5 % of
polyhydrohybutyrates (PHB) was produced. These low values in biopolymers
production could be related to the lack of microorganisms as cyanobacteria that
are able to accumulate carbohydrates/PHB.


  The evolutionary dynamics of molecular populations are strongly dependent on
the structure of genotype spaces. The map between genotype and phenotype
determines how easily genotype spaces can be navigated and the accessibility of
evolutionary innovations. In particular, the size of neutral networks
corresponding to specific phenotypes and its statistical counterpart, the
distribution of phenotype abundance, have been studied through multiple
computationally tractable genotype-phenotype maps. In this work, we test a
theory that predicts the abundance of a phenotype and the corresponding
asymptotic distribution (given the compositional variability of its genotypes)
through the exact enumeration of several GP maps. Our theory predicts with high
accuracy phenotype abundance, and our results show that, in navigable genotype
spaces ---characterised by the presence of large neutral networks---, phenotype
abundance converges to a log-normal distribution.


  Reconciling gene trees with a species tree is a fundamental problem to
understand the evolution of gene families. Many existing approaches reconcile
each gene tree independently. However, it is well-known that the evolution of
gene families is interconnected. In this paper, we extend a previous approach
to reconcile a set of gene trees with a species tree based on segmental
macro-evolutionary events, where segmental duplication events and losses are
associated with cost $\delta$ and $\lambda$, respectively. We show that the
problem is polynomial-time solvable when $\delta \leq \lambda$ (via
LCA-mapping), while if $\delta > \lambda$ the problem is NP-hard, even when
$\lambda = 0$ and a single gene tree is given, solving a long standing open
problem on the complexity of the reconciliation problem. On the positive side,
we give a fixed-parameter algorithm for the problem, where the parameters are
$\delta/\lambda$ and the number $d$ of segmental duplications, of time
complexity $O(\lceil \frac{\delta}{\lambda} \rceil^{d} \cdot n \cdot
\frac{\delta}{\lambda})$. Finally, we demonstrate the usefulness of this
algorithm on two previously studied real datasets: we first show that our
method can be used to confirm or refute hypothetical segmental duplications on
a set of 16 eukaryotes, then show how we can detect whole genome duplications
in yeast genomes.


  We examine salient trends of influenza pandemics in Australia, a rapidly
urbanizing nation. To do so, we implement state-of-the-art influenza
transmission and progression models within a large-scale stochastic computer
simulation, generated using comprehensive Australian census datasets from 2006,
2011, and 2016. Our results offer a simulation-based investigation of a
population's sensitivity to pandemics across multiple historical time points,
and highlight three significant trends in pandemic patterns over the years:
increased peak prevalence, faster spreading rates, and decreasing
spatiotemporal bimodality. We attribute these pandemic trends to increases in
two key quantities indicative of urbanization: population fraction residing in
major cities, and international air traffic. In addition, we identify features
of the pandemic's geographic spread that we attribute to changes in the
commuter mobility network. The generic nature of our model and the ubiquity of
urbanization trends around the world make it likely for our results to be
applicable in other rapidly urbanizing nations.


  Membrane transporters contribute to the regulation of the internal
environment of cells by translocating substrates across cell membranes. Like
all physical systems, the behaviour of membrane transporters is constrained by
the laws of thermodynamics. However, many mathematical models of transporters,
especially those incorporated into whole-cell models, are not thermodynamically
consistent, leading to unrealistic behaviour. In this paper we use a
physics-based modelling framework, in which the transfer of energy is
explicitly accounted for, to develop thermodynamically consistent models of
transporters. We then apply this methodology to model two specific
transporters: the cardiac sarcoplasmic/endoplasmic Ca$^{2+}$ ATPase (SERCA) and
the cardiac Na$^+$/K$^+$ ATPase.


  Antibodies are a critical part of the immune system, having the function of
directly neutralising or tagging undesirable objects (the antigens) for future
destruction. Being able to predict which amino acids belong to the paratope,
the region on the antibody which binds to the antigen, can facilitate antibody
design and contribute to the development of personalised medicine. The
suitability of deep neural networks has recently been confirmed for this task,
with Parapred outperforming all prior physical models. Our contribution is
twofold: first, we significantly outperform the computational efficiency of
Parapred by leveraging \`a trous convolutions and self-attention. Secondly, we
implement cross-modal attention by allowing the antibody residues to attend
over antigen residues. This leads to new state-of-the-art results on this task,
along with insightful interpretations.


  In this paper we analyze the length-spectrum of blocks in
$\gamma$-structures. $\gamma$-structures are a class of RNA pseudoknot
structures that plays a key role in the context of polynomial time RNA folding.
A $\gamma$-structure is constructed by nesting and concatenating specific
building components having topological genus at most $\gamma$. A block is a
substructure enclosed by crossing maximal arcs with respect to the partial
order induced by nesting. We show that, in uniformly generated
$\gamma$-structures, there is a significant gap in this length-spectrum, i.e.,
there asymptotically almost surely exists a unique longest block of length at
least $n-O(n^{1/2})$ and that with high probability any other block has finite
length. For fixed $\gamma$, we prove that the length of the longest block
converges to a discrete limit law, and that the distribution of short blocks of
given length tends to a negative binomial distribution in the limit of long
sequences. We refine this analysis to the length spectrum of blocks of specific
pseudoknot types, such as H-type and kissing hairpins. Our results generalize
the rainbow spectrum on secondary structures by the first and third authors and
are being put into context with the structural prediction of long non-coding
RNAs.


  We define a general class of models representing natural selection between
two alleles. The population size and spatial structure are arbitrary, but
fixed. Genetics can be haploid, diploid, or otherwise; reproduction can be
asexual or sexual. Biological events (e.g.~births, deaths, mating, dispersal)
depend in arbitrary fashion on the current population state. Our formalism is
based on the idea of genetic sites. Each genetic site resides at a particular
locus and houses a single allele. Each individual contains a number of sites
equal to its ploidy (one for haploids, two for diploids, etc.). Selection
occurs via replacement events, in which alleles in some sites are replaced by
copies of others. Replacement events depend stochastically on the population
state, leading to a Markov chain representation of natural selection. Within
this formalism, we define reproductive value, fitness, neutral drift, and
fixation probability, and prove relationships among them. We identify four
criteria for evaluating which allele is selected and show that these become
equivalent in the limit of low mutation. We then formalize the method of weak
selection. The power of our formalism is illustrated with applications to
evolutionary games on graphs and to selection in a haplodiploid population.


  Chlamydia trachomatis (Ct) is the most common bacterial sexually transmitted
infection (STI) in the United States and is major cause of infertility, pelvic
inflammatory disease, and ectopic pregnancy among women. Despite decades of
screening women for Ct, rates continue to increase among them in high prevalent
areas such as New Orleans. A pilot study in New Orleans found approximately 11%
of 14-24 year old of African Americans (AAs) were infected with Ct.
  Our goal is to mathematically model the impact of different interventions for
AA men resident in New Orleans on the general rate of Ct among women resident
at the same region. We create and analyze mathematical models such as
multi-risk and continuous-risk compartmental models and agent-based network
model to first help understand the spread of Ct and second evaluate and
estimate behavioral and biomedical interventions including condom-use,
screening, partner notification, social friend notification, and rescreening.
Our compartmental models predict the Ct prevalence is a function of the number
of partners for a person, and quantify how this distribution changes as a
function of condom-use. We also observe that although increased Ct screening
and rescreening, and treating partners of infected people will reduce the
prevalence, these mitigations alone are not sufficient to control the epidemic.
A combination of both sexual partner and social friend notification is needed
to mitigate Ct.


  Understanding cell identity is an important task in many biomedical areas.
Expression patterns of specific marker genes have been used to characterize
some limited cell types, but exclusive markers are not available for many cell
types. A second approach is to use machine learning to discriminate cell types
based on the whole gene expression profiles (GEPs). The accuracies of simple
classification algorithms such as linear discriminators or support vector
machines are limited due to the complexity of biological systems. We used deep
neural networks to analyze 1040 GEPs from 16 different human tissues and cell
types. After comparing different architectures, we identified a specific
structure of deep autoencoders that can encode a GEP into a vector of 30
numeric values, which we call the cell identity code (CIC). The original GEP
can be reproduced from the CIC with an accuracy comparable to technical
replicates of the same experiment. Although we use an unsupervised approach to
train the autoencoder, we show different values of the CIC are connected to
different biological aspects of the cell, such as different pathways or
biological processes. This network can use CIC to reproduce the GEP of the cell
types it has never seen during the training. It also can resist some noise in
the measurement of the GEP. Furthermore, we introduce classifier autoencoder,
an architecture that can accurately identify cell type based on the GEP or the
CIC.


  The log-det distance between two aligned DNA sequences was introduced as a
tool for statistically consistent inference of a gene tree under simple
non-mixture models of sequence evolution. Here we prove that the log-det
distance, coupled with a distance-based tree construction method, also permits
consistent inference of species trees under mixture models appropriate to
aligned genomic-scale sequences data. Data may include sites from many genetic
loci, which evolved on different gene trees due to incomplete lineage sorting
on an ultrametric species tree, with different time-reversible substitution
processes. The simplicity and speed of distance-based inference suggests
log-det based methods should serve as benchmarks for judging more elaborate and
computationally-intensive species trees inference methods.


  We analyze the fluctuation of the number of individuals when two competing
species, beginning with a few initial individuals, are submitted to a logistic
growth. We show that when the total number of individuals reaches the carrying
capacity, the number of each species is subject to giant fluctuations (variance
$\sim$ mean 2) if the two species have similar growth rate. We show that the
deterministic logistic equation can be used only when the growth rates are
significantly different, otherwise such growth has to be investigated by
stochastic processes tools. These results generalize to a wide class of growth
law.


  In the present work, we review the fundamental methods which have been
developed in the last few years for classifying into families and clans the
distribution of amino acids in protein databases. This is done through
functions of random variables, the Entropy Measures of probabilities of
occurrence of the amino acids. An intensive study of the Pfam databases is
presented with restrictions to families which could be represented by
rectangular arrays of amino acids with m rows (protein domains) and n columns
(amino acids). This work is also an invitation to scientific research groups
worldwide to undertake the statistical analysis with different numbers of rows
and columns since we believe in the mathematical characterization of the
distribution of amino acids as a fundamental insight on the determination of
protein structure and evolution.


  Despite the advent of several novel, synthetic gene drive mechanisms and
their potential to one-day control a number of devastating diseases, among
other applications, practical use of these systems remains contentious and
risky. In particular, there is little in the way of empirical evidence of the
long-term robustness of these synthetic systems against mutational breakdown.
Rather, most existing systems are either known or predicted to be susceptible
to rapid inactivation, though methodological designs continue to be refined.
Here we evaluate a currently existing synthetic, underdominance-based gene
drive system 200+ generations after it was first established in a laboratory
colony of Drosophila melanogaster. Not only do we find that the system is still
functioning as designed, we also show evidence that disruptions to the genetic
construct are highly likely to be removed by natural selection, contributing to
the system's robust, long-term stability. This stability appears to be a result
of a fundamental relationship between ribosomal proteins (a novel target of the
system) and natural cellular defenses that protect against cancer development.
As far as we are aware, this is the longest continually functioning synthetic
gene drive system thus verified, making it highly appropriate for additional
research into its eventual suitability for field trials. Due to inherent
properties of this gene drive, it is also likely to be adaptable for use in
many different species. The insect lines established and used to test this
system have been deposited at a Drosophila stock center, and are available to
labs for further, independent testing.


  In this paper, a simple dynamical model in which fractal networks are formed
by self-organized critical (SOC) dynamics is proposed; the proposed model
consists of growth and collapse processes. It has been shown that SOC dynamics
are realized by the combined processes in the model. Thus, the distributions of
the cluster size and collapse size follow a power-law function in the
stationary state. Moreover, through SOC dynamics, the networks become fractal
in nature. The criticality of SOC dynamics is the same as the universality
class of mean-field theory. The model explains the possibility that the fractal
nature in complex networks emerges by SOC dynamics in a manner similar to the
case with fractal objects embedded in a Euclidean space.


  The Devil Facial Tumour Disease (DFTD), a unique case of a transmissible
cancer, had a devastating effect on its host, the Tasmanian Devil. Current
estimates of its density are at roughly 20% of the pre-disease state, and
single-population epidemiological models have predicted the likely extinction
of the host. Here we take advantage of extensive surveys across Tasmania
providing data on the spatial and temporal spread of DFTD, and investigate the
dynamics of this host-pathogen system using a spatial metapopulation model. We
first confirm a most likely origin of DFTD in the north-east corner of the
island, and then use the inferred dynamics to predict the fate of the species.
We find that our medium-term predictions match additional data not used for
fitting, and that on the long-term, Tasmanian Devils are predicted to coexist
with the tumour. The key process allowing persistence is the repeated
reinvasion of extinct patches from neighbouring areas where the disease has
flared up and died out, resulting in a dynamic equilibrium with different
levels of spatial heterogeneity. However, this dynamic equilibrium is predicted
to keep this apex predator at about 9 % of its original density, with possible
dramatic effects on the Tasmanian ecosystem.


  An adaptive agent predicting the future state of an environment must weigh
trust in new observations against prior experiences. In this light, we propose
a view of the adaptive immune system as a dynamic Bayesian machinery that
updates its memory repertoire by balancing evidence from new pathogen
encounters against past experience of infection to predict and prepare for
future threats. This framework links the observed initial rapid increase of the
memory pool early in life followed by a mid-life plateau to the ease of
learning salient features of sparse environments. We also derive a modulated
memory pool update rule in agreement with current vaccine response experiments.
Our results suggest that pathogenic environments are sparse and that memory
repertoires significantly decrease infection costs even with moderate sampling.
The predicted optimal update scheme maps onto commonly considered competitive
dynamics for antigen receptors.


  The minimal number of rooted subtree prune and regraft (rSPR) operations
needed to transform one phylogenetic tree into another one induces a metric on
phylogenetic trees - the rSPR-distance. The rSPR-distance between two
phylogenetic trees $T$ and $T'$ can be characterised by a maximum agreement
forest; a forest with a minimum number of components that covers both $T$ and
$T'$. The rSPR operation has recently been generalised to phylogenetic networks
with, among others, the subnetwork prune and regraft (SNPR) operation. Here, we
introduce maximum agreement graphs as an explicit representations of
differences of two phylogenetic networks, thus generalising maximum agreement
forests. We show that maximum agreement graphs induce a metric on phylogenetic
networks - the agreement distance. While this metric does not characterise the
distances induced by SNPR and other generalisations of rSPR, we prove that it
still bounds these distances with constant factors.


  One interesting question in love relationships is: finally, what and when is
the end of this love relationship? Using a prey-predator
Verhulst-Lotka-Volterra (VLV) model we imply cooperation and competition
tendency between people in order to describe a "love dilemma game". We select
the most simple but immediately most complex case for studying the set of
nonlinear differential equations, i.e. that implying three persons, being at
the same time prey and predator. We describe four different scenarios in such a
love game containing either a one-way love or a love triangle. Our results show
that it is hard to love more than one person simultaneously. Moreover, to love
several people simultaneously is an unstable state. We find some condition in
which persons tend to have a friendly relationship and love someone in spite of
their antagonistic interaction. We demonstrate the dynamics by displaying flow
diagrams.


  Population annealing is a powerful tool for large-scale Monte Carlo
simulations. We adapt this method to molecular dynamics simulations and
demonstrate its excellent accelerating effect by simulating the folding of a
short peptide commonly used to gauge the performance of algorithms. The method
is compared to the well established parallel tempering approach and is found to
yield similar performance for the same computational resources. In contrast to
other methods, however, population annealing scales to a nearly arbitrary
number of parallel processors and it is thus a unique tool that enables
molecular dynamics to tap into the massively parallel computing power available
in supercomputers that is so much needed for a range of difficult computational
problems.


  We review algorithms for protein design in general. Although these algorithms
have a rich combinatorial, geometric, and mathematical structure, they are
almost never covered in computer science classes. Furthermore, many of these
algorithms admit provable guarantees of accuracy, soundness, complexity,
completeness, optimality, and approximation bounds. The algorithms represent a
delicate and beautiful balance between discrete and continuous computation and
modeling, analogous to that which is seen in robotics, computational geometry,
and other fields in computational science. Finally, computer scientists may be
unaware of the almost direct impact of these algorithms for predicting and
introducing molecular therapies that have gone in a short time from mathematics
to algorithms to software to predictions to preclinical testing to clinical
trials. Indeed, the overarching goal of these algorithms is to enable the
development of new therapeutics that might be impossible or too expensive to
discover using experimental methods. Thus the potential impact of these
algorithms on individual, community, and global health has the potential to be
quite significant.


  The epidemic threshold is probably the most studied quantity in the modelling
of epidemics on networks. For a large class of networks and dynamics the
epidemic threshold is well studied and understood. However, it is less so for
clustered networks where theoretical results are mostly limited to idealised
networks. In this paper we focus on a class of models known as pairwise models
where, to our knowledge, no analytical result for the epidemic threshold
exists. We show that by exploiting the presence of fast variables and using
some standard techniques from perturbation theory we are able to obtain the
epidemic threshold analytically. We validate this new threshold by comparing it
to the numerical solution of the full system. The agreement is found to be
excellent over a wide range of values of the clustering coefficient,
transmission rate and average degree of the network. Interestingly, we find
that the analytical form of $R_0$ depends on the choice of closure,
highlighting the importance of model choice when dealing with real-world
epidemics. Nevertheless, we ex- pect that our method will extend to other
systems in which fast variables are present.


  We consider the task of estimating a high-dimensional directed acyclic graph,
given observations from a linear structural equation model with arbitrary noise
distribution. By exploiting properties of common random graphs, we develop a
new algorithm that requires conditioning only on small sets of variables. The
proposed algorithm, which is essentially a modified version of the
PC-Algorithm, offers significant gains in both computational complexity and
estimation accuracy. In particular, it results in more efficient and accurate
estimation in large networks containing hub nodes, which are common in
biological systems. We prove the consistency of the proposed algorithm, and
show that it also requires a less stringent faithfulness assumption than the
PC-Algorithm. Simulations in low and high-dimensional settings are used to
illustrate these findings. An application to gene expression data suggests that
the proposed algorithm can identify a greater number of clinically relevant
genes than current methods.


  Stochastic formalisms are necessary to describe the behavior of many
biological systems. However, there remains a lack of numerical methods
available to calculate the stationary probability distributions of stochastic
reaction networks. We have previously development a numerical approach to
calculate stationary probability distributions of stochastic networks, named ZI
closure scheme. In this work, we present a free applications based on
ZI-closure scheme, called ZICS.


  The emerging era of personalized medicine relies on medical decisions,
practices, and products being tailored to the individual patient. Point-of-care
systems, at the heart of this model, play two important roles. First, they are
required for identifying subjects for optimal therapies based on their genetic
make-up and epigenetic profile. Second, they will be used for assessing the
progression of such therapies. Central to this vision is designing systems
that, with minimal user-intervention, can transduce complex signals from
biosystems in complement with clinical information to inform medical decision
within point-of-care settings. To reach our ultimate goal of developing
point-of-care systems and realizing personalized medicine, we are taking a
multistep systems-level approach towards understanding cellular processes and
biomolecular profiles, to quantify disease states and external interventions.


  Open-ended evolution requires unbounded possibilities that evolving entities
can explore. The cardinality of a set of those possibilities thus has a
significant implication for the open-endedness of evolution. We propose that
facilitating formation of higher-order entities is a generalizable, effective
way to cause a "cardinality leap" in the set of possibilities that promotes
open-endedness. We demonstrate this idea with a simple, proof-of-concept toy
model called "Hash Chemistry" that uses a hash function as a fitness evaluator
of evolving entities of any size/order. Simulation results showed that the
cumulative number of unique replicating entities that appeared in evolution
increased almost linearly along time without an apparent bound, demonstrating
the effectiveness of the proposed cardinality leap. It was also observed that
the number of individual entities involved in a single replication event
gradually increased over time, indicating evolutionary appearance of
higher-order entities. Moreover, these behaviors were not observed in control
experiments in which fitness evaluators were replaced by random number
generators. This strongly suggests that the dynamics observed in Hash Chemistry
were indeed evolutionary behaviors driven by selection and adaptation taking
place at multiple scales.


  Chemical evolution is essential in understanding the origins of life. We
present a theory for the evolution of molecule masses and show that small
molecules grow by random diffusion and large molecules by a preferential
attachment process leading eventually to life's molecules. It reproduces
correctly the distribution of molecules found via mass spectroscopy for the
Murchison meteorite and estimates the start of chemical evolution back to 12.8
billion years following the birth of stars and supernovae. From the Frontier
mass between the random and preferential attachment dynamics the birth time of
molecule families can be estimated. Amino acids emerge about 165 million years
after chemical elements emerge in stars. Using the scaling of reaction rates
with the distance of the molecules in space we recover correctly the few days
emergence time of amino acids in the Miller-Urey experiment. The distribution
of interstellar and extragalactic molecules are both consistent with the
evolutionary mass distribution, and their age is estimated to 108 and 65
million years after the start of evolution. From the model, we can determine
the number of different molecule compositions at the time of the emergence of
Earth to be 1.6 million and the number of molecule compositions in interstellar
space to a mere 719 species.


  We study the effect of biological confounders on the model selection problem
between Kingman coalescents with population growth, and Xi-coalescents
involving simultaneous multiple mergers. We use a low dimensional,
computationally tractable summary statistic, dubbed the singleton-tail
statistic, to carry out approximate likelihood ratio tests between these model
classes. The singleton-tail statistic has been shown to distinguish between
them with high power in the simple setting of neutrally evolving, panmictic
populations without recombination. We extend this work by showing that cryptic
recombination and selection do not diminish the power of the test, but that
misspecifying population structure does. Furthermore, we demonstrate that the
singleton-tail statistic can also solve the more challenging model selection
problem between multiple mergers due to selective sweeps, and multiple mergers
due to high fecundity with moderate power of up to 30%.


  A central goal of protein-folding theory is to predict the stochastic
dynamics of transition paths --- the rare trajectories that transit between the
folded and unfolded ensembles --- using only thermodynamic information, such as
a low-dimensional equilibrium free-energy landscape. However, commonly used
one-dimensional landscapes typically fall short of this aim, because an
empirical coordinate-dependent diffusion coefficient has to be fit to
transition-path trajectory data in order to reproduce the transition-path
dynamics. We show that an alternative, first-principles free-energy landscape
predicts transition-path statistics that agree well with simulations and
single-molecule experiments without requiring dynamical data as an input. This
'topological configuration' model assumes that distinct, native-like
substructures assemble on a timescale that is slower than native-contact
formation but faster than the folding of the entire protein. Using only
equilibrium simulation data to determine the free energies of these
coarse-grained intermediate states, we predict a broad distribution of
transition-path transit times that agrees well with the transition-path
durations observed in simulations. We further show that both the distribution
of finite-time displacements on a one-dimensional order parameter and the
ensemble of transition-path trajectories generated by the model are consistent
with the simulated transition paths. These results indicate that a landscape
based on transient folding intermediates, which are often hidden by
one-dimensional projections, can form the basis of a predictive model of
protein-folding transition-path dynamics.


  We study the challenges of applying deep learning to gene expression data. We
find experimentally that there exists non-linear signal in the data, however is
it not discovered automatically given the noise and low numbers of samples used
in most research. We discuss how gene interaction graphs (same pathway,
protein-protein, co-expression, or research paper text association) can be used
to impose a bias on a deep model similar to the spatial bias imposed by
convolutions on an image. We explore the usage of Graph Convolutional Neural
Networks coupled with dropout and gene embeddings to utilize the graph
information. We find this approach provides an advantage for particular tasks
in a low data regime but is very dependent on the quality of the graph used. We
conclude that more work should be done in this direction. We design experiments
that show why existing methods fail to capture signal that is present in the
data when features are added which clearly isolates the problem that needs to
be addressed.


  The physical mechanism behind the spontaneous formation of metachronal waves
in microtubule arrays in a low Reynolds number fluid has been of interest for
the past several years, yet is still not well understood. We present a model
implementing the hydrodynamic coupling hypothesis from first principles, and
use this model to simulate kinesin-driven microtubule arrays and observe their
emergent behavior. The results of simulations are compared to known
experimental observations by Sanchez et al. By varying parameters, we determine
regimes in which the metachronal wave phenomenon emerges, and categorize other
types of possible microtubule motion outside these regimes.


  We use plasmon rulers to follow the conformational dynamics of a single
protein for up to 24 h at a video rate. The plasmon ruler consists of two gold
nanospheres connected by a single protein linker. In our experiment, we follow
the dynamics of the molecular chaperone heat shock protein 90, which is known
to show open and closed conformations. Our measurements confirm the previously
known conformational dynamics with transition times in the second to minute
time scale and reveals new dynamics on the time scale of minutes to hours.
Plasmon rulers thus extend the observation bandwidth 3/4 orders of magnitude
with respect to single-molecule fluorescence resonance energy transfer and
enable the study of molecular dynamics with unprecedented precision.


  In recent studies, non-coding protein RNAs have been identified as microRNA
that can be used as biomarkers for early diagnosis and treatment of cancer,
that decrease mortality in cancer. A microRNA may target hundreds or thousands
of genes and a gene may regulate several microRNAs, so determining which
microRNA is associated with which cancer is a big challenge. Many computational
methods have been performed to detect micoRNAs association with cancer, but
more effort is needed with higher accuracy. Increasing research has shown that
relationship between microRNAs and TFs play a significant role in the diagnosis
of cancer. Therefore, we developed a new computational framework (CAMIRADA) to
identify cancer-related microRNA based on the relationship between microRNAs
and disease genes (DG) in the protein network, the functional relationships
between microRNAs and Transcription Factors (TF) on the co-expression network,
and the relationship between microRNAs and the Differential Expression Gene
(DEG) on co-expression network. The CAMIRADA was applied to assess breast
cancer data from two HMDD and miR2Disease databases. In this study, the AUC for
the 65 microRNA of the top of the list was 0.95, which was more accurate than
the similar methods used to detect microRNA associated with the cancer artery.


  To explore how the topology of interaction networks determines the robustness
of dynamical systems, we study the antisymmetric Lotka-Volterra equation
(ALVE). The ALVE is the replicator equation of zero-sum games in evolutionary
game theory, in which the strengths of pairwise interactions between strategies
are defined by an antisymmetric matrix such that typically some strategies go
extinct over time. Here we show that there also exist topologically robust
zero-sum games, such as the rock-paper-scissors game, for which all strategies
coexist for all choices of interaction strengths. We refer to such zero-sum
games as coexistence networks and construct coexistence networks with an
arbitrary number of strategies. By mapping the long-time dynamics of the ALVE
to the algebra of antisymmetric matrices, we identify simple graph-theoretical
rules by which coexistence networks are constructed. Examples are
triangulations of cycles characterized by the golden ratio $\varphi =
1.6180...$, cycles with complete subnetworks, and non-Hamiltonian networks. In
graph-theoretical terms, we extend the concept of a Pfaffian orientation from
even-sized to odd-sized networks. Our results show that the topology of
interaction networks alone can determine the long-time behavior of nonlinear
dynamical systems, and may help to identify robust network motifs arising, for
example, in ecology.


  Protein gamma-turn prediction is useful in protein function studies and
experimental design. Several methods for gamma-turn prediction have been
developed, but the results were unsatisfactory with Matthew correlation
coefficients (MCC) around 0.2-0.4. One reason for the low prediction accuracy
is the limited capacity of the methods; in particular, the traditional
machine-learning methods like SVM may not extract high-level features well to
distinguish between turn or non-turn. Hence, it is worthwhile exploring new
machine-learning methods for the prediction. A cutting-edge deep neural
network, named Capsule Network (CapsuleNet), provides a new opportunity for
gamma-turn prediction. Even when the number of input samples is relatively
small, the capsules from CapsuleNet are very effective to extract high-level
features for classification tasks. Here, we propose a deep inception capsule
network for gamma-turn prediction. Its performance on the gamma-turn benchmark
GT320 achieved an MCC of 0.45, which significantly outperformed the previous
best method with an MCC of 0.38. This is the first gamma-turn prediction method
utilizing deep neural networks. Also, to our knowledge, it is the first
published bioinformatics application utilizing capsule network, which will
provides a useful example for the community.


  Cell competition is recognized to be instrumental to the dynamics and
structure of the tumor-host interface in invasive cancers. In mild competition
scenarios, the healthy tissue and cancer cells can coexist. When the
competition is aggressive, competitive cells, the so called super-competitors,
expand by killing other cells. Novel cytotoxic drugs and molecularly targeted
drugs are commonly administered as part of cancer therapy. Both types of drugs
are susceptible to various mechanisms of drug resistance, obstructing or
preventing a successful outcome. In this paper, we develop a cancer growth
model that accounts for the competition between cancer cells and healthy cells.
The model incorporates resistance to both cytotoxic and targeted drugs. In both
cases, the level of drug resistance is assumed to be a continuous variable
ranging from fully-sensitive to fully-resistant. Using our model we demonstrate
that when the competition is moderate, therapies using both drugs are more
effective compared with single drug therapies. However, when cancer cells are
highly competitive, targeted drugs become more effective. In this case,
therapies that are initiated with a targeted drug and are exposed to it for a
sufficiently long time are shown to have better outcomes. The results of the
study stress the importance of adjusting the therapy to the pre-treatment
resistance levels. We conclude with a study of the spatiotemporal propagation
of drug resistance in a competitive setting, verifying that the same
conclusions hold in the spatially heterogeneous case.


  Motivation: Drug discovery demands rapid quantification of compound-protein
interaction (CPI). However, there is a lack of methods that can predict
compound-protein affinity from sequences alone with high applicability,
accuracy, and interpretability.
  Results: We present a seamless integration of domain knowledges and
learning-based approaches. Under novel representations of
structurally-annotated protein sequences, a semi-supervised deep learning model
that unifies recurrent and convolutional neural networks has been proposed to
exploit both unlabeled and labeled data, for jointly encoding molecular
representations and predicting affinities. Our representations and models
outperform conventional options in achieving relative error in IC$_{50}$ within
5-fold for test cases and 20-fold for protein classes not included for
training. Performances for new protein classes with few labeled data are
further improved by transfer learning. Furthermore, separate and joint
attention mechanisms are developed and embedded to our model to add to its
interpretability, as illustrated in case studies for predicting and explaining
selective drug-target interactions. Lastly, alternative representations using
protein sequences or compound graphs and a unified RNN/GCNN-CNN model using
graph CNN (GCNN) are also explored to reveal algorithmic challenges ahead.
  Availability: Data and source codes are available at
https://github.com/Shen-Lab/DeepAffinity
  Supplementary Information: Supplementary data are available at
http://shen-lab.github.io/deep-affinity-bioinf18-supp-rev.pdf


  Multidrug resistance consists of a series of genetic and epigenetic
alternations that involve multifactorial and complex processes, which are a
challenge to successful cancer treatments. Accompanied by advances in
biotechnology and high-dimensional data analysis techniques that are bringing
in new opportunities in modeling biological systems with continuous phenotypic
structured models, we study a cancer cell population model that considers a
multi-dimensional continuous resistance trait to multiple drugs to investigate
multidrug resistance. We compare our continuous resistance trait model with
classical models that assume a discrete resistance state and classify the cases
when the continuum and discrete models yield different dynamical patterns in
the emerging heterogeneity in response to drugs. We also compute the maximal
fitness resistance trait for various continuum models and study the effect of
epimutations. Finally, we demonstrate how our approach can be used to study
tumor growth regarding the turnover rate and the proliferating fraction, and
show that a continuous resistance level may result in a different dynamics when
compared with the predictions of other discrete models.


  Rooted phylogenetic networks are rooted acyclic digraphs. They are used to
model complex evolution where hybridization, recombination and other
reticulation events play important roles. A rigorous definition of network
compression is introduced on the basis of the recent studies of the
relationships between cluster, tree and rooted phylogenetic network. The
concept reveals another interesting connection between the two well-studied
network classes|tree-child networks and reticulation-visible networks|and
enables us to define a new class of networks for which the cluster containment
problem has a linear-time algorithm.


  1. Theoretical models pertaining to feedbacks between ecological and
evolutionary processes are prevalent in multiple biological fields. An
integrative overview is currently lacking, due to little crosstalk between the
fields and the use of different methodological approaches.
  2. Here we review a wide range of models of eco-evolutionary feedbacks and
highlight their underlying assumptions. We discuss models where feedbacks occur
both within and between hierarchical levels of ecosystems, including
populations, communities, and abiotic environments, and consider feedbacks
across spatial scales.
  3. Identifying the commonalities among feedback models, and the underlying
assumptions, helps us better understand the mechanistic basis of
eco-evolutionary feedbacks. Eco-evolutionary feedbacks can be readily modelled
by coupling demographic and evolutionary formalisms. We provide an overview of
these approaches and suggest future integrative modelling avenues.
  4. Our overview highlights that eco-evolutionary feedbacks have been
incorporated in theoretical work for nearly a century. Yet, this work does not
always include the notion of rapid evolution or concurrent ecological and
evolutionary time scales. We discuss the importance of density- and
frequency-dependent selection for feedbacks, as well as the importance of
dispersal as a central linking trait between ecology and evolution in a spatial
context.


  The recent development of metagenomic sequencing makes it possible to
sequence microbial genomes including viruses in an environmental sample.
Identifying viral sequences from metagenomic data is critical for downstream
virus analyses. The existing reference-based and gene homology-based methods
are not efficient in identifying unknown viruses or short viral sequences. Here
we have developed a reference-free and alignment-free machine learning method,
DeepVirFinder, for predicting viral sequences in metagenomic data using deep
learning techniques. DeepVirFinder was trained based on a large number of viral
sequences discovered before May 2015. Evaluated on the sequences after that
date, DeepVirFinder outperformed the state-of-the-art method VirFinder at all
contig lengths. Enlarging the training data by adding millions of purified
viral sequences from environmental metavirome samples significantly improves
the accuracy for predicting under-represented viruses. Applying DeepVirFinder
to real human gut metagenomic samples from patients with colorectal carcinoma
(CRC) identified 51,138 viral sequences belonging to 175 bins. Ten bins were
associated with the cancer status, indicating their potential use for
non-invasive diagnosis of CRC. In summary, DeepVirFinder greatly improved the
precision and recall rates of viral identification, and it will significantly
accelerate the discovery rate of viruses.


  Molecular motors walk along filaments until they detach stochastically with a
force-dependent unbinding rate. Here, we show that this unbinding rate can be
obtained from the analysis of experimental data of molecular motors moving in
stationary optical traps. Two complementary methods are presented, based on the
analysis of the distribution for the unbinding forces and of the motor's force
traces. In the first method, analytically derived force distributions for slip
bonds, slip-ideal bonds, and catch bonds are used to fit the cumulative
distributions of the unbinding forces. The second method is based on the
statistical analysis of the observed force traces. We validate both methods
with stochastic simulations and apply them to experimental data for kinesin-1.


  The burnt-bridges ratchet (BBR) mechanism is a model for biased molecular
motion whereby the construct destroys track binding sites as it progresses, and
therefore acts as a diffusing forager, seeking new substrate sites. Using Monte
Carlo simulations that implement the Gillespie algorithm, we investigate the
kinetic characteristics of simple polyvalent BBRs as they move on tracks of
increasing width. We find that as the track width is increased the BBRs remain
nearly ballistic for considerable track widths proportional to the span (leg
length) of the polyvalent walker, before transitioning to near-conventional
diffusion on two-dimensional tracks. We find there exists a tradeoff in BBR
track association time and superdiffusivity in the BBR design parameter space
of span, polyvalency and track width. Furthermore, we develop an analytical
model to describe the ensembleaverage motion on the track and find it is in
good agreement with our Gillespie simulation results. This work offers insights
into design criteria for de novo BBRs and their associated tracks, where
experimentalists seek to optimize directionality and track association time.


  This article reports prototype web apps that use commodity, open-source
technologies for augmented and virtual reality to provide immersive,
interactive human-computer interfaces for chemistry, structural biology and
related disciplines. The examples, which run in any standard web browser and
are accessible at
https://lucianoabriata.altervista.org/jsinscience/arjs/armodeling/ together
with demo videos, showcase applications that could go well beyond pedagogy,
i.e. advancing actual utility in research settings: molecular visualization at
atomistic and coarse-grained levels in interactive immersive 3D, coarse-grained
modeling of molecular physics and chemistry, and on-the-fly calculation of
experimental observables and overlay onto experimental data. From this
playground, I depict perspectives on how these emerging technologies might
couple in the future to neural network-based quantum mechanical calculations,
advanced forms of human-computer interaction such as speech-based
communication, and sockets for concurrent collaboration through the internet
-all technologies that are today maturing in web browsers- to deliver the next
generation of tools for truly interactive, immersive molecular modeling that
can streamline human thought and intent with the numerical processing power of
computers.


  A two-dimensional (2D) model is developed to describe the growth of spotty
patterns (domains) in the lipid membranes which result from the coagulation of
ion channels. It is assumed that the ion channels can coagulate when activated,
e.g. the neurotransmitter may act as the surface-active substance for some ion
channels in the lipid matrix or the activation may result in increase of
polarity of the ion channels that give rise to their coagulation process due to
the Brownian diffusive motion. The results show that the domain radius rD of
the ion channels scales with time as rD~time^0.4.The slowing and subsequent
saturation of the domain radius at longer time periods are in qualitative
agreement with experiments.


  Mutation rate is a key determinant of the pace as well as outcome of
evolution, and variability in this rate has been shown in different scenarios
to play a key role in evolutionary adaptation and resistance evolution under
stress caused by selective pressure. Here we investigate the dynamics of
resistance fixation in a bacterial population with variable mutation rates and
show that evolutionary outcomes are most sensitive to mutation rate variations
when the population is subject to environmental and demographic conditions that
suppress the evolutionary advantage of high-fitness subpopulations. By directly
mapping a biophysical fitness function to the system-level dynamics of the
population we show that both low and very high, but not intermediate, levels of
stress in the form of an antibiotic result in a disproportionate effect of
hypermutation on resistance fixation. We demonstrate how this behavior is
directly tied to the extent of genetic hitchhiking in the system, the
propagation of high-mutation rate cells through association with high-fitness
mutations. Our results indicate a substantial role for mutation rate
flexibility in the evolution of antibiotic resistance under conditions that
present a weak advantage over wildtype to resistant cells.


  The likelihood ratio statistic, with its asymptotic $\chi^2$ distribution at
regular model points, is often used for hypothesis testing. At model
singularities and boundaries, however, the asymptotic distribution may not be
$\chi^2$, as highlighted by recent work of Drton. Indeed, poor behavior of a
$\chi^2$ for testing near singularities and boundaries is apparent in
simulations, and can lead to conservative or anti-conservative tests. Here we
develop a new distribution designed for use in hypothesis testing near
singularities and boundaries, which asymptotically agrees with that of the
likelihood ratio statistic. For two example trinomial models, arising in the
context of inference of evolutionary trees, we show the new distributions
outperform a $\chi^2$.


  We present a novel approach to the study of epidemics on networks as
thermodynamic phenomena, considering the thermodynamic efficiency of
contagions, considered as distributed computational processes. Modelling SIS
dynamics on a contact network statistical-mechanically, we follow the Maximum
Entropy principle to obtain steady state distributions and derive, under
certain assumptions, relevant thermodynamic quantities both analytically and
numerically. In particular, we obtain closed form solutions for some cases,
while interpreting key epidemic variables, such as the reproductive ratio $R_0$
of a SIS model, in a statistical mechanical setting. On the other hand, we
consider configuration and free entropy, as well as the Fisher Information, in
the epidemiological context. This allowed us to identify criticality and
distinct phases of epidemic processes. For each of the considered thermodynamic
quantities, we compare the analytical solutions informed by the Maximum Entropy
principle with the numerical estimates for SIS epidemics simulated on
Watts-Strogatz random graphs.


  Non-invasive prenatal testing or NIPT is currently among the top researched
topic in obstetric care. While the performance of the current state-of-the-art
NIPT solutions achieve high sensitivity and specificity, they still struggle
with a considerable number of samples that cannot be concluded with certainty.
Such uninformative results are often subject to repeated blood sampling and
re-analysis, usually after two weeks, and this period may cause a stress to the
future mothers as well as increase the overall cost of the test. We propose a
supplementary method to traditional z-scores to reduce the number of such
uninformative calls. The method is based on a novel analysis of the length
profile of circulating cell free DNA which compares the change in such profiles
when random-based and length-based elimination of some fragments is performed.
The proposed method is not as accurate as the standard z-score; however, our
results suggest that combination of these two independent methods correctly
resolves a substantial portion of healthy samples with an uninformative result.
Additionally, we discuss how the proposed method can be used to identify
maternal aberrations, thus reducing the risk of false positive and false
negative calls.
  Keywords: Next-generation sequencing, Cell-free DNA, Uninformative result,
Method, Trisomy, Prenatal testing


  We study hairpin folding dynamics by means of extensive computer simulations,
with particular attention paid to the influence of helicity on the folding time
$\tau$. We find that the dynamical exponent $\alpha$ of the anomalous scaling
$\tau \sim N^\alpha$ for a hairpin with length N changes from 1.6 ($1+\nu$) to
1.2 ($2\nu$) in three dimensions, when duplex helicity is removed. The relation
$\alpha = 2\nu$ in rotationless hairpin folding is further verified in two
dimensions ($\nu = 0.75$), and for a ghost-chain ($\nu = 0.5$). This, to our
knowledge, is the first observation of the theoretical lower bound on $\alpha$,
which was predicted earlier on the basis of energy conservation for polymer
translocation through a pore. Our findings suggest that the folding dynamics in
long helical chains is governed by the duplex dynamics, contrasting the earlier
understanding based on the stem-flower picture of unpaired segments. We propose
a scaling argument for $\alpha = 1+\nu$ in helical chains, assuming that duplex
relaxation required for orientational positioning of the next pair of bases is
the rate-limiting process.


  Muller's ratchet describes the irreversible accumulation of deleterious
mutations in asexual populations. In well-mixed populations the speed of
fitness decline is exponentially small in the population size, and any positive
rate of beneficial mutations is sufficient to reverse the ratchet in large
populations. The behavior is fundamentally different in populations with
spatial structure, because the speed of the ratchet remains nonzero in the
infinite size limit when the deleterious mutation rate exceeds a critical
value. Based on the relation between the spatial ratchet and directed
percolation, we develop a scaling theory incorporating both deleterious and
beneficial mutations. The theory is verified by extensive simulations in one
and two dimensions.


  We describe an "embarrassingly parallel" method for Bayesian phylogenetic
inference, annealed Sequential Monte Carlo, based on recent advances in the
Sequential Monte Carlo literature such as adaptive determination of annealing
parameters. The algorithm provides an approximate posterior distribution over
trees and evolutionary parameters as well as an unbiased estimator for the
marginal likelihood. This unbiasedness property can be used for the purpose of
testing the correctness of posterior simulation software. We evaluate the
performance of phylogenetic annealed Sequential Monte Carlo by reviewing and
comparing with other computational Bayesian phylogenetic methods, in
particular, different marginal likelihood estimation methods. Unlike previous
Sequential Monte Carlo methods in phylogenetics, our annealed method can
utilize standard Markov chain Monte Carlo tree moves and hence benefit from the
large inventory of such moves available in the literature. Consequently, the
annealed Sequential Monte Carlo method should be relatively easy to incorporate
into existing phylogenetic software packages based on Markov chain Monte Carlo
algorithms. We illustrate our method using simulation studies and real data
analysis.


  Acute respiratory infections have epidemic and pandemic potential and thus
are being studied worldwide, albeit in many different contexts and study
formats. Predicting infection from symptom data is critical, though using
symptom data from varied studies in aggregate is challenging because the data
is collected in different ways. Accordingly, different symptom profiles could
be more predictive in certain studies, or even symptoms of the same name could
have different meanings in different contexts. We assess state-of-the-art
transfer learning methods for improving prediction of infection from symptom
data in multiple types of health care data ranging from clinical, to home-visit
as well as crowdsourced studies. We show interesting characteristics regarding
six different study types and their feature domains. Further, we demonstrate
that it is possible to use data collected from one study to predict infection
in another, at close to or better than using a single dataset for prediction on
itself. We also investigate in which conditions specific transfer learning and
domain adaptation methods may perform better on symptom data. This work has the
potential for broad applicability as we show how it is possible to transfer
learning from one public health study design to another, and data collected
from one study may be used for prediction of labels for another, even collected
through different study designs, populations and contexts.


  Diagnosis and risk stratification of cancer and many other diseases require
the detection of genomic breakpoints as a prerequisite of calling copy number
alterations (CNA). This, however, is still challenging and requires
time-consuming manual curation. As deep-learning methods outperformed classical
state-of-the-art algorithms in various domains and have also been successfully
applied to life science problems including medicine and biology, we here
propose Deep SNP, a novel Deep Neural Network to learn from genomic data.
Specifically, we used a manually curated dataset from 12 genomic single
nucleotide polymorphism array (SNPa) profiles as truth-set and aimed at
predicting the presence or absence of genomic breakpoints, an indicator of
structural chromosomal variations, in windows of 40,000 probes. We compare our
results with well-known neural network models as well as Rawcopy though this
tool is designed to predict breakpoints and in addition genomic segments with
high sensitivity. We show, that Deep SNP is capable of successfully predicting
the presence or absence of a breakpoint in large genomic windows and
outperforms state-of-the-art neural network models. Qualitative examples
suggest that integration of a localization unit may enable breakpoint detection
and prediction of genomic segments, even if the breakpoint coordinates were not
provided for network training. These results warrant further evaluation of
DeepSNP for breakpoint localization and subsequent calling of genomic segments.


  Taylor's Law (TL) relates the variance to the mean of a random variable via
power law. In ecology it applies to populationsand it is a common empirical
pattern shared among different ecosystems. Measurements give power law exponent
to be between 1 and 2, and more often to cluster around 2, whereas theoretical
models predict TL exponent can assume any real value. In this paper, adopting
the framework of multiplicative growth models in a Markovian environment, we
investigate the possibility of evolutionary strategies to be responsible for TL
exponent to be in a finite range. We implement three different strategies the
individuals can follow and for each strategy set two different optimization
investment objectives. In all the studied cases we find TL exponent can assume
any real value due to the existence of regions of the model parameters in which
the exponent can diverge. Furthermore, under natural hypothesis on the dynamics
of the environment, the shapes of these regions do not depend on different
strategies adopted and nor on the optimization objective. Thus the introduction
of strategies dose not affect the range of TL exponent in the model. In our
theoretical framework rare events are shaping the value of the TL exponent,
suggesting, as hinted by previous works, that empirical values may be a
statistical artifact following from under sampling.


  Markov State Models (MSMs) are a powerful framework to reproduce the
long-time conformational dynamics of biomolecules using a set of short
Molecular Dynamics (MD) simulations. However, precise kinetics predictions of
MSMs heavily rely on the features selected to describe the system. Despite the
importance of feature selection for large system, determining an optimal set of
features remains a difficult unsolved problem. Here, we introduce an automatic
approach to optimize feature selection based on genetic algorithms (GA), which
adaptively evolves the most fitted solution according to natural selection
laws. The power of the GA-based method is illustrated on long atomistic folding
simulations of four proteins, varying in length from 28 to 80 residues. Due to
the diversity of tested proteins, we expect that our method will be extensible
to other proteins and drive MSM building to a more objective protocol.


  Diffuse pollution in rural areas due to agricultural run-off is a widespread
and difficult problem to address due to the vast areas affected. Drainage
channels do receive these polluted waters, but its introduction in the
conventional treatment network is unfeasible. Within this context,
microalgae-based treatment systems could be used as alternative treatment
plants. In this study, a new design of semi-closed (hybrid) tubular horizontal
photobioreactor (HTH-PBR) with low energy requirements has been evaluated for
microalgae cultivation at full-scale (8.5 m3), using agricultural runoff as
feedstock. This novel system was tested in batch and continuous mode during a
total of 4 and 135 days. Considering a full-scale application in an
agricultural context, a batch test was carried out to evaluate the performance
of the system. An increase of 22% in the biomass concentration in 4 days was
registered, and all nutrients were consumed during the first two days. In the
continuous experiment (December- April), a productivity between 2-14 g TSS/m3d
was reached in winter, whereas values up to 76.4 g TSS/m3d were reached at the
end of the study in spring, despite the low nutrients concentration in the
feedstock. The elimination of emerging contaminants was also evaluated,
obtaining the highest removals for the fragrances tonalide and galaxolide (73%
and 68%), and the antiinflamatory diclofenac (61%).


  The design of novel proteins has many applications but remains an attritional
process with success in isolated cases. Meanwhile, deep learning technologies
have exploded in popularity in recent years and are increasingly applicable to
biology due to the rise in available data. We attempt to link protein design
and deep learning by using variational autoencoders to generate protein
sequences conditioned on desired properties. Potential copper and calcium
binding sites are added to non-metal binding proteins without human
intervention and compared to a hidden Markov model. In another use case, a
grammar of protein structures is developed and used to produce sequences for a
novel protein topology. One candidate structure is found to be stable by
molecular dynamics simulation. The ability of our model to confine the vast
search space of protein sequences and to scale easily has the potential to
assist in a variety of protein design tasks.


  Probabilistic inference provides a language for describing how organisms may
learn from and adapt to their environment. The computations needed to implement
probabilistic inference often require specific representations, akin to having
the suitable data structures for implementing certain algorithms in computer
programming. Yet it is unclear how such representations can be instantiated in
the stochastic, parallel-running biochemical machinery found in cells (such as
single-celled organisms). Here, we show how representations for supporting
inference in Markov models can be embodied in cellular circuits, by combining a
concentration-dependent scheme for encoding probabilities with a mechanism for
directional counting. We show how the logic of protein production and
degradation constrains the computation we set out to implement. We argue that
this process by which an abstract computation is shaped by its biochemical
realization strikes a compromise between "rationalistic" information-processing
perspectives and alternative approaches that emphasize embodiment.


  Stochastic quasi-cycles for a two species model of the excitatory-inhibitory
type, arranged on a triangular loop, are studied. By increasing the strength of
the inter-nodes coupling, one moves the system towards the Hopf bifurcation and
the amplitude of the stochastic oscillations are consequently magnified. When
the system is instead constrained to evolve on specific manifolds, selected so
as to return a constant rate of deterministic damping for the perturbations,
the observed amplification correlates with the degree of non normal reactivity,
here quantified by the numerical abscissa. The thermodynamics of the reactive
loop is also investigated and the degree of inherent reactivity shown to
facilitate the out-of-equilibrium exploration of the available phase space.


  The Protein Data Bank (PDB) contains the atomic structures of over 105
biomolecules with better than 2.8A resolution. The listing of the identities
and coordinates of the atoms comprising each macromolecule permits an analysis
of the slow-time vibrational response of these large systems to minor
perturbations. 3D video animations of individual modes of oscillation
demonstrate how regions interdigitate to create cohesive collective motions,
providing a comprehensive framework for and familiarity with the overall 3D
architecture. Furthermore, the isolation and representation of the softest,
slowest deformation coordinates provide opportunities for the development of
me- chanical models of enzyme function. The eigenvector decomposition,
therefore, must be accurate, reliable as well as rapid to be generally reported
upon. We obtain the eigenmodes of a 1.2A 34kDa PDB entry using either
exclusively heavy atoms or partly or fully reduced atomic sets; Cartesian or
internal coordinates; interatomic force fields derived either from a full
Cartesian potential, a reduced atomic potential or a Gaussian
distance-dependent potential; and independently devel- oped software. These
varied technologies are similar in that each maintains proper stereochemistry
either by use of dihedral degrees of freedom which freezes bond lengths and
bond angles, or by use of a full atomic potential that includes realistic bond
length and angle restraints. We find that the shapes of the slowest
eigenvectors are nearly identical, not merely similar.


  Bazhin has analyzed ATP coupling in terms of quasiequilibrium states where
fast reactions have reached an approximate steady state while slow reactions
have not yet reached equilibrium. After an expository introduction to the
relevant aspects of reaction network theory, we review his work and explain the
role of emergent conserved quantities in coupling. These are quantities, left
unchanged by fast reactions, whose conservation forces exergonic processes such
as ATP hydrolysis to drive desired endergonic processes.


  Plasmids are autonomously replicating genetic elements in bacteria. At cell
division plasmids are distributed among the two daughter cells. This gene
transfer from one generation to the next is called vertical gene transfer. We
study the dynamics of a bacterial population carrying plasmids and are in
particular interested in the long-time distribution of plasmids. Starting with
a model for a bacterial population structured by the discrete number of
plasmids, we proceed to the continuum limit in order to derive a continuous
model. The model incorporates plasmid reproduction, division and death of
bacteria, and distribution of plasmids at cell division. It is a hyperbolic
integro-differential equation and a so-called growth-fragmentation-death model.
As we are interested in the long-time distribution of plasmids we study the
associated eigenproblem and show existence of eigensolutions. The stability of
this solution is studied by analyzing the spectrum of the integro-differential
operator given by the eigenproblem. By relating the spectrum with the spectrum
of an integral operator we find a simple real dominating eigenvalue with a
non-negative corresponding eigenfunction. Moreover, we describe an iterative
method for the numerical construction of the eigenfunction.


  The T-cell (TCR) repertoire relies on the diversity of receptors composed of
two chains, called $\alpha$ and $\beta$, to recognize pathogens. Using results
of high throughput sequencing and computational chain-pairing experiments of
human TCR repertoires, we quantitively characterize the $\alpha\beta$
generation process. We estimate the probabilities of a rescue recombination of
the $\beta$ chain on the second chromosome upon failure or success on the first
chromosome. Unlike $\beta$ chains, $\alpha$ chains recombine simultaneously on
both chromosomes, resulting in correlated statistics of the two genes which we
predict using a mechanistic model. We find that $\sim 28 \%$ of cells express
both $\alpha$ chains. We report that clones sharing the same $\beta$ chain but
different $\alpha$ chains are overrepresented, suggesting that they respond to
common immune challenges. Altogether, our statistical analysis gives a complete
quantitative mechanistic picture that results in the observed correlations in
the generative process. We learn that the probability to generate any
TCR$\alpha\beta$ is lower than $10^{-12}$ and estimate the generation diversity
and sharing properties of the $\alpha\beta$ TCR repertoire.


  Hierarchy is an efficient way for a group to organize, but often goes along
with inequality that benefits leaders. To control despotic behaviour, followers
can assess leaders decisions by aggregating their own and their neighbours
experience, and in response challenge despotic leaders. But in hierarchical
social networks, this interactional justice can be limited by (i) the high
influence of a small clique who are treated better, and (ii) the low
connectedness of followers. Here we study how the connectedness of a social
network affects the co-evolution of despotism in leaders and tolerance to
despotism in followers. We simulate the evolution of a population of agents,
where the influence of an agent is its number of social links. Whether a leader
remains in power is controlled by the overall satisfaction of group members, as
determined by their joint assessment of the leaders behaviour. We demonstrate
that centralization of a social network around a highly influential clique
greatly increases the level of despotism. This is because the clique is more
satisfied, and their higher influence spreads their positive opinion of the
leader throughout the network. Finally, our results suggest that increasing the
connectedness of followers limits despotism while maintaining hierarchy.


  An analytical expression is derived for the transition path time distribution
for a one-dimensional particle crossing of a parabolic barrier. Two cases are
analyzed: (i) A non-Markovian process described by a generalized Langevin
equation with a power-law memory kernel and (ii) a Markovian process with a
noise violating the fluctuation-dissipation theorem, modeling the stochastic
dynamics generated by active forces. In the case (i) we show that the anomalous
dynamics strongly affecting the short time behavior of the distributions, but
this happens only for very rare events not influencing the overall statistics.
At long times the decay is always exponential, in disagreement with a recent
study suggesting a stretched exponential decay. In the case (ii) the active
forces do not substantially modify the short time behavior of the distribution,
but lead to an overall decrease of the average transition path time. These
findings offer some novel insights, useful for the analysis of experiments of
transition path times in (bio)molecular systems.


  Metabolism and evolution are closely connected: if a mutation incurs extra
energetic costs for an organism, there is a baseline selective disadvantage
that may or may not be compensated for by other adaptive effects. A
long-standing, but to date unproven, hypothesis is that this disadvantage is
equal to the fractional cost relative to the total resting metabolic
expenditure. This hypothesis has found a recent resurgence as a powerful tool
for quantitatively understanding the strength of selection among different
classes of organisms. Our work explores the validity of the hypothesis from
first principles through a generalized metabolic growth model, versions of
which have been successful in describing organismal growth from single cells to
higher animals. We build a mathematical framework to calculate how
perturbations in maintenance and synthesis costs translate into contributions
to the selection coefficient. This allows us to show that the hypothesis is an
approximation to the actual baseline selection coefficient. Moreover we can
derive the correct prefactor in its functional form, as well as analytical
bounds on its accuracy for any given realization of the model. We illustrate
our general framework using a special case of the growth model, which we show
provides a quantitative description of metabolic expenditures in data collected
from a wide array of unicellular organisms (both prokaryotes and eukaryotes).
In all these cases we demonstrate that the hypothesis is an excellent
approximation, allowing estimates of baseline selection coefficients to within
15% of their actual values. Even in a broader biological parameter range,
covering growth data from multicellular organisms, the hypothesis continues to
work well, always within an order of magnitude of the correct result. Our work
thus justifies its use as a versatile tool, setting the stage for its wider
deployment.


  The analysis of the demographic transition of the past century and a half,
using both empirical data and mathematical models, has rendered a wealth of
well-established facts, including the dramatic increases in life expectancy.
Despite these insights, such analyses have also occasionally triggered debates
which spill over many disciplines, from genetics, to biology, or demography.
Perhaps the hottest discussion is happening around the question of maximum
human lifespan, which --besides its fascinating historical and philosophical
interest-- poses urgent pragmatic warnings on a number of issues in public and
private decision-making. In this paper, we add to the controversy some results
which, based on purely statistical grounds, suggest that the maximum human
lifespan is not fixed, or has not reached yet a plateau. Quite the contrary,
analysis on reliable data for over 150 years in more than 20 industrialized
countries point at a sustained increase in the maximum age at death.
Furthermore, were this trend to continue, a limitless lifespan could be
achieved by 2102. Finally, we quantify the dependence of increases in the
maximum lifespan on socio-economic factors. Our analysis indicates that in some
countries the observed rising patterns can only be sustained by progressively
larger increases in GDP, setting the problem of longevity in a context of
diminishing returns.


  We introduce a new method for transforming chemical systems into desired
logical operators (e.g. NAND-gates) or similar signal-manipulation components.
The method is based upon open-loop dynamic regulation, where external
conditions such as feed-rate, lighting conditions, etc. are modulated according
to a prescribed temporal sequence that is independent of the input to the
network. The method is first introduced using a simple didactic model. We then
show its application in transforming a well-stirred cubic autocatalytic
reaction (often referred to as the Selkov-Gray-Scott model) into a logical
NAND-gate. We also comment on the applicability of the method to biological and
other systems.


  The emerging field at the intersection of quantitative biology, network
modeling, and control theory has enjoyed significant progress in recent years.
This Special Issue brings together a selection of papers on complementary
approaches to observe, identify, and control biological and biologically
inspired networks. These approaches advance the state of the art in the field
by addressing challenges common to many such networks, including high
dimensionality, strong nonlinearity, uncertainty, and limited opportunities for
observation and intervention. Because these challenges are not unique to
biological systems, it is expected that many of the results presented in these
contributions will also find applications in other domains, including physical,
social, and technological networks.


  We consider linear elimination of variables in steady state equations of a
chemical reaction network. Particular subsets of variables corresponding to
sets of so-called reactant-noninteracting species, are introduced. The steady
state equations for the variables in such a set, taken together with potential
linear conservation laws in the variables, define a linear system of equations.
We give conditions that guarantee that the solution to this system is
nonnegative, provided it is unique. The results are framed in terms of spanning
forests of a particular multidigraph derived from the reaction network and
thereby conditions for uniqueness and nonnegativity of a solution are derived
by means of the multidigraph. Though our motivation comes from applications in
systems biology, the results have general applicability in applied sciences.


  Segmental duplications (SDs), or low-copy repeats (LCR), are segments of DNA
greater than 1 Kbp with high sequence identity that are copied to other regions
of the genome. SDs are among the most important sources of evolution, a common
cause of genomic structural variation, and several are associated with diseases
of genomic origin. Despite their functional importance, SDs present one of the
major hurdles for de novo genome assembly due to the ambiguity they cause in
building and traversing both state-of-the-art overlap-layout-consensus and de
Bruijn graphs. This causes SD regions to be misassembled, collapsed into a
unique representation, or completely missing from assembled reference genomes
for various organisms. In turn, this missing or incorrect information limits
our ability to fully understand the evolution and the architecture of the
genomes. Despite the essential need to accurately characterize SDs in
assemblies, there is only one tool that has been developed for this purpose,
called Whole Genome Assembly Comparison (WGAC). WGAC is comprised of several
steps that employ different tools and custom scripts, which makes it difficult
and time consuming to use. Thus there is still a need for algorithms to
characterize within-assembly SDs quickly, accurately, and in a user friendly
manner.
  Here we introduce a SEgmental Duplication Evaluation Framework (SEDEF) to
rapidly detect SDs through sophisticated filtering strategies based on Jaccard
similarity and local chaining. We show that SEDEF accurately detects SDs while
maintaining substantial speed up over WGAC that translates into practical run
times of minutes instead of weeks. Notably, our algorithm captures up to 25%
pairwise error between segments, where previous studies focused on only 10%,
allowing us to more deeply track the evolutionary history of the genome.
  SEDEF is available at https://github.com/vpc-ccg/sedef