  Proteins are only moderately stable. It has long been debated whether this
narrow range of stabilities is solely a result of neutral drift towards lower
stability or purifying selection against excess stability is also at work - for
which no experimental evidence was found so far. Here we show that mutations
outside the active site in the essential E. coli enzyme adenylate kinase result
in stability-dependent increase in substrate inhibition by AMP, thereby
impairing overall enzyme activity at high stability. Such inhibition caused
substantial fitness defects not only in the presence of excess substrate but
also under physiological conditions. In the latter case, substrate inhibition
caused differential accumulation of AMP in the stationary phase for the
inhibition prone mutants. Further, we show that changes in flux through Adk
could accurately describe the variation in fitness effects. Taken together,
these data suggest that selection against substrate inhibition and hence excess
stability may have resulted in a narrow range of optimal stability observed for
modern proteins.


  Biomolecular machines transduce free energy from one form to another to
fulfill many important roles inside cells, with dissipation required to achieve
directed progress. We investigate how to break time-reversal symmetry at a
given dissipation cost by using deterministic protocols to drive systems over
sawtooth potentials, which have frequently been used to model molecular
machines as ratchets. Time asymmetry increases for sawtooth potentials with
higher barriers and for driving potentials of intermediate width. For systems
driven over a sawtooth potential according to a protocol, we find that
symmetric sawtooths maximize time asymmetry, while earlier work examining
ratchet models of molecular machines required asymmetric sawtooth potentials to
achieve directed behavior. This distinction arises because deterministically
driven machines are externally provided with direction, while autonomous
machines must generate directed behavior.


  The basic model of purposeful kinesis was developed recently (Gorban and
Cabukoglu, Ecol. Complex. 2018, 33, 75 83) on the basis of the let well enough
alone idea. According to this model the diffusion drops while the reproduction
coefficient is increasing. That is, species prefer to stay in a good condition
and the population gives birth; otherwise, in the bad situation individuals
want to run away because of the fatal conditions. In this study, we analyse the
impact of the purposeful kinesis model on running waves. The running waves in
the population with kinesis are studied using numerical experiments. Both
monotonic and non monotonic (Allee effect) dependence of the reproduction
coefficient on the population density are studied. The possible benefits of the
purposeful kinesis are demonstrated: with the higher diffusion, while the
population without kinesis ends up with extinction, the population with kinesis
stays alive and has the running wave behavior. While the kinesis of the prey
population is decreasing, the wave amplitude gets smaller. On the other hand,
for the lower kinesis of predators they have a sharp increase.


  The computational problem of inferring the full haplotype of a cell starting
from read sequencing data is known as haplotype assembly, and consists in
assigning all heterozygous Single Nucleotide Polymorphisms (SNPs) to exactly
one of the two chromosomes. Indeed, the knowledge of complete haplotypes is
generally more informative than analyzing single SNPs and plays a fundamental
role in many medical applications. To reconstruct the two haplotypes, we
addressed the weighted Minimum Error Correction (wMEC) problem, which is a
successful approach for haplotype assembly. This NP-hard problem consists in
computing the two haplotypes that partition the sequencing reads into two
disjoint sub-sets, with the least number of corrections to the SNP values. To
this aim, we propose here GenHap, a novel computational method for haplotype
assembly based on Genetic Algorithms, yielding optimal solutions by means of a
global search process. In order to evaluate the effectiveness of our approach,
we run GenHap on two synthetic (yet realistic) datasets, based on the Roche/454
and PacBio RS II sequencing technologies. We compared the performance of GenHap
against HapCol, an efficient state-of-the-art algorithm for haplotype phasing.
Our results show that GenHap always obtains high accuracy solutions (in terms
of haplotype error rate), and is up to 4x faster than HapCol in the case of
Roche/454 instances and up to 20x faster when compared on the PacBio RS II
dataset. Finally, we assessed the performance of GenHap on two different real
datasets. Future-generation sequencing technologies, producing longer reads
with higher coverage, can highly benefit from GenHap, thanks to its capability
of efficiently solving large instances of the haplotype assembly problem.


  Molecular motion through pores plays a crucial role in various natural and
industrial processes. One of the most fascinating features of biological
channel-facilitated transport is a stochastic gating process, when the channels
dynamically fluctuate between several conformations during the translocation.
Although this phenomenon has been intensively investigated, many properties of
translocation in dynamically changing environment remain not well understood
microscopically. We developed a discrete-state stochastic framework to analyze
the molecular mechanisms of transport processes with stochastic gating by
explicitly calculating molecular fluxes through the pores. Two scenarios are
specifically investigated: 1) symmetry preserving stochastic gating with
free-energy changes, and 2) stochastic gating with symmetry changes but without
modifications in the overall particle-pore interactions. It is found that
stochastic gating can both accelerate or slow down the molecular translocation
depending on the specific parameters of the system. We argue that biological
systems might optimize their performance by utilizing conformational
fluctuations of channels. Our theoretical analysis clarifies physical-chemical
aspects of the molecular mechanisms of transport with stochastic gating.


  The issue of longevity has been long time recognized as one of key concepts
in demography. A particular aspect of longevity addressed in this paper is the
difference in average observed duration of life for female and male population,
called observed lifespan differential. Using the data from Human Mortality
Database, the dynamics of the observed lifespan differential is studied for a
large number of countries worldwide from 1960 to 2014. An interesting
phenomenon that the growing trend of the observed lifespan differential at the
beginning of the studied interval does not persist, i.e. that it reverts to
stagnation or even decline is revealed for a large majority of countries in the
dataset. In a number of case studies, a strong association of lifespan dynamic
with disruptive events such as wars, dissolutions or integrations of states or
policy measures is demonstrated. Finally, a novel method of calculating the
observed lifespan differential of a population from mortality indicators of its
subpopulations is introduced and applied to the analysis of the observed
lifespan differential in Israel from 1990 to 2000.


  Determination of \emph{optimal} arrangements of $N$ particles on a sphere is
a well-known problem in physics. A famous example of such is the Thomson
problem of finding equilibrium configurations of electrical charges on a
sphere. More recently however, similar problems involving other potentials and
non-spherical domains have arisen in biophysical systems. Many optimal
configurations have previously been computed, especially for the Thomson
problem, however few results exist for potentials that correspond to more
applied problems. Here we numerically compute optimal configurations
corresponding to the \emph{narrow escape} and \emph{narrow capture} problems in
biophysics. We provide comprehensive tables of global energy minima for
$N\leq120$ and local energy minima for $N\leq65$, and we exclude all saddle
points. Local minima up to $N=120$ are available online.


  Although algebraic graph theory based models have been widely applied in
physical modeling and molecular studies, they are typically incompetent in the
analysis and prediction of biomolecular properties when compared with other
quantitative approaches. There is a need to explore the capability and
limitation of algebraic graph theory for molecular and biomolecular modeling,
analysis, and prediction. In this work, we propose novel algebraic graph
learning (AGL) models that encode high-dimensional physical and biological
information into intrinsically low-dimensional representations. The proposed
AGL model introduces multiscale weighted colored subgraphs to describe crucial
molecular and biomolecular interactions via graph invariants associated with
the graph Laplacian, its pseudo-inverse, and adjacent matrix. Additionally, the
AGL models are incorporated with an advanced machine learning algorithm to
connect the low-dimensional graph representation of biomolecular structures
with their macroscopic properties. Three popular protein-ligand binding
affinity benchmarks, namely CASF-2007, CASF-2013, and CASF-2016, are employed
to validate the accuracy, robustness, and reliability of the present AGL model.
Numerical results indicate that the proposed AGL method outperforms the other
state-of-the-art methods in the binding affinity predictions of the
protein-ligand complexes.


  We consider the problem of estimating species trees from unrooted gene tree
topologies in the presence of incomplete lineage sorting, a common phenomenon
that creates gene tree heterogeneity in multilocus datasets. One popular class
of reconstruction methods in this setting is based on internode distances, i.e.
the average graph distance between pairs of species across gene trees. While
statistical consistency in the limit of large numbers of loci has been
established in some cases, little is known about the sample complexity of such
methods. Here we make progress on this question by deriving a lower bound on
the worst-case variance of internode distance which depends linearly on the
corresponding graph distance in the species tree. We also discuss some
algorithmic implications.


  By developing and leveraging an explicit molecular realisation of a
measurement-and-feedback-powered Szilard engine, we investigate the extraction
of work from complex environments by minimal machines with finite capacity for
memory and decision-making. Living systems perform inference to exploit complex
structure, or correlations, in their environment, but the physical limits and
underlying cost/benefit trade-offs involved in doing so remain unclear. To
probe these questions, we consider a minimal model for a structured environment
- a correlated sequence of molecules - and explore mechanisms based on extended
Szilard engines for extracting the work stored in these non-equilibrium
correlations. We consider systems limited to a single bit of memory making
binary 'choices' at each step. We demonstrate that increasingly complex
environments allow increasingly sophisticated inference strategies to extract
more energy than simpler alternatives, and argue that optimal design of such
machines should also consider the energy reserves required to ensure robustness
against fluctuations due to mistakes.


  Galled trees are widely studied as a recombination model in population
genetics. This class of phylogenetic networks is generalized into galled
networks by relaxing a structural condition. In this work, a linear recurrence
formula is given for counting 1-galled networks, which are galled networks
satisfying the condition that each reticulate node has only one leaf
descendant. Since every galled network consists of a set of 1-galled networks
stacked one on top of the other, a method is also presented to count and
enumerate galled networks.


  We propose an analysis for a novel experiment intended to locate the genetic
locus Xce (X-chromosome controlling element), which biases the stochastic
process of X-inactivation in the mouse. X-inactivation bias is a phenomenon
where cells in the embryo randomly choose one parental chromosome to
inactivate, but show an average bias towards one parental strain. Measurement
of allele-specific gene-expression through pyrosequencing was conducted on
mouse crosses of an uncharacterized parent with known carriers. Our Bayesian
analysis is suitable for this adaptive experimental design, accounting for the
biases and differences in precision among genes. Model identifiability is
facilitated by priors constrained to a manifold. We show that reparameterized
slice-sampling can suitably tackle a general class of constrained priors. We
demonstrate a physical model, based upon a "weighted-coin" hypothesis, that
predicts X-inactivation ratios in untested crosses. This model suggests that
Xce alleles differ due to a process known as copy number variation, where
stronger Xce alleles are shorter sequences.


  We explore simple models aimed at the study of social contagion, in which
contagion proceeds through two stages. When coupled with demographic turnover,
we show that two-stage contagion leads to nonlinear phenomena which are not
present in the basic `classical' models of mathematical epidemiology. These
include: bistability, critical transitions, endogenous oscillations, and
excitability, suggesting that contagion models with stages could account for
some aspects of the complex dynamics encountered in social life. These
phenomena, and the bifurcations involved, are studied by a combination of
analytical and numerical means.


  The Hassell model has been widely used as a general discrete-time population
dynamics model that describes both contest and scramble intraspecific
competition through a tunable exponent. Since the two types of competition
generally lead to different degrees of inequality in the resource distribution
among individuals, the exponent is expected to be related to this inequality.
However, among various first-principles derivations of this model, none is
consistent with this expectation. This paper explores whether a Hassell model
with an exponent related to inequality in resource allocation can be derived
from first principles. Indeed, such a Hassell model can be derived by assuming
random competition for resources among the individuals wherein each individual
can obtain only a fixed amount of resources at a time. Changing the size of the
resource unit alters the degree of inequality, and the exponent changes
accordingly. The Beverton-Holt and Ricker models can be regarded as special
cases of the derived Hassell model. Two additional Hassell models are derived
under some modified assumptions.


  With the rise of computers, simulation models have emerged beside the more
traditional statistical and mathematical models as a third pillar for
ecological analysis. Broadly speaking, a simulation model is an algorithm,
typically implemented as a computer program, which propagates the states of a
system forward. Unlike in a mathematical model, however, this propagation does
not employ the methods of calculus but rather a set of rules or formulae that
directly prescribe the next state. Such an algorithmic model specification is
particularly suited for describing systems that are difficult to capture or
analyze with differential equations such as: (a) systems that are highly
nonlinear or chaotic; (b) discrete systems, for example networks or groups of
distinct individuals; (c) systems that are stochastic; and (d) systems that are
too complex to be successfully treated with classical calculus. As these
situations are frequently encountered in ecology, simulation models are now
widely applied across the discipline. They have been instrumental in developing
new insights into classical questions of species' coexistence, community
assembly, population dynamics, biogeography, and many more. The methods for
this relatively young field are still being actively developed, and practical
work with simulation models requires ecologists to learn new skills such as
coding, sensitivity analysis, calibration, validation, and forecasting
uncertainties. Moreover, scientific inquiry with complex systems has led to
subtle changes to the philosophical and epistemological views regarding
simplicity, reductionism, and the relationship between prediction and
understanding.


  A key focus in current cancer research is the discovery of cancer biomarkers
that allow earlier detection with high accuracy and lower costs for both
patients and hospitals. Blood samples have long been used as a health status
indicator, but DNA methylation signatures in blood have not been fully
appreciated in cancer research. Historically, analysis of cancer has been
conducted directly with the patient's tumor or related tissues. Such analyses
allow physicians to diagnose a patient's health and cancer status; however,
physicians must observe certain symptoms that prompt them to use biopsies or
imaging to verify the diagnosis. This is a post-hoc approach. Our study will
focus on epigenetic information for cancer detection, specifically information
about DNA methylation in human peripheral blood samples in cancer discordant
monozygotic twin-pairs. This information might be able to help us detect cancer
much earlier, before the first symptom appears. Several other types of
epigenetic data can also be used, but here we demonstrate the potential of
blood DNA methylation data as a biomarker for pan-cancer using SAS 9.3 and SAS
EM. We report that 55 methylation CpG sites measurable in blood samples can be
used as biomarkers for early cancer detection and classification.


  In Costa Rica, the first known cases of Zika were reported in 2016. We looked
at the 2016-2017 Zika outbreak and explored the transmission dynamics using
weekly reported data. A nonlinear differential equation single-outbreak model
with sexual transmission, as well as host availability for vector-feeding was
used to estimate key parameters, fit the data and compute the {\it basic
reproductive number}, $\mathcal{R}_0$, distribution. Furthermore, a sensitivity
and elasticity analysis was computed based on the $\mathcal{R}_0$ parameters.


  Cyclization of DNA with sticky ends is commonly used to construct DNA
minicircles and to measure DNA bendability. The cyclization probability of
short DNA (< 150 bp) has a strong length dependence, but how it depends on the
rotational positioning of the sticky ends around the helical axis is less
clear. To shed light upon the determinants of the cyclization probability of
short DNA, we measured cyclization and decyclization rates of ~100-bp DNA with
sticky ends over two helical periods using single-molecule Fluorescence
Resonance Energy Transfer (FRET). The cyclization rate increases monotonically
with length, indicating no excess twisting, while the decyclization rate
oscillates with length, higher at half-integer helical turns and lower at
integer helical turns. The oscillation profile is kinetically and
thermodynamically consistent with a three-state cyclization model in which
sticky-ended short DNA first bends into a torsionally-relaxed teardrop, and
subsequently transitions to a more stable loop upon terminal base stacking. We
also show that the looping probability density (the J factor) extracted from
this study is in good agreement with the worm-like chain model near 100 bp. For
shorter DNA, we discuss various experimental factors that prevent an accurate
measurement of the J factor.


  Living systems are capable of locomotion, reconfiguration, and replication.
To perform these tasks, cells spatiotemporally coordinate the interactions of
force-generating, "active" molecules that create and manipulate non-equilibrium
structures and force fields that span up to millimeter length scales [1-3].
Experimental active matter systems of biological or synthetic molecules are
capable of spontaneously organizing into structures [4,5] and generating global
flows [6-9]. However, these experimental systems lack the spatiotemporal
control found in cells, limiting their utility for studying non-equilibrium
phenomena and bioinspired engineering. Here, we uncover non-equilibrium
phenomena and principles by optically controlling structures and fluid flow in
an engineered system of active biomolecules. Our engineered system consists of
purified microtubules and light-activatable motor proteins that crosslink and
organize microtubules into distinct structures upon illumination. We develop
basic operations, defined as sets of light patterns, to create, move, and merge
microtubule structures. By composing these basic operations, we are able to
create microtubule networks that span several hundred microns in length and
contract at speeds up to an order of magnitude faster than the speed of an
individual motor. We manipulate these contractile networks to generate and
sculpt persistent fluid flows. The principles of boundary-mediated control we
uncover may be used to study emergent cellular structures and forces and to
develop programmable active matter devices.


  Non-conding RNAs play a key role in the post-transcriptional regulation of
mRNA translation and turnover in eukaryotes. miRNAs, in particular, interact
with their target RNAs through protein-mediated, sequence-specific binding,
giving rise to extended and highly heterogeneous miRNA-RNA interaction
networks. Within such networks, competition to bind miRNAs can generate an
effective positive coupling between their targets. Competing endogenous RNAs
(ceRNAs) can in turn regulate each other through miRNA-mediated crosstalk.
Albeit potentially weak, ceRNA interactions can occur both dynamically,
affecting e.g. the regulatory clock, and at stationarity, in which case ceRNA
networks as a whole can be implicated in the composition of the cell's
proteome. Many features of ceRNA interactions, including the conditions under
which they become significant, can be unraveled by mathematical and in silico
models. We review the understanding of the ceRNA effect obtained within such
frameworks, focusing on the methods employed to quantify it, its role in the
processing of gene expression noise, and how network topology can determine its
reach.


  A strategy based on photo sequencing batch operation was used to select
cyanobacteria over unsettled green algae in a wastewater treatment system,
evaluating for the first time the effect of hydraulic regimes on nutritional
dynamics and microorganisms competition. During 30 days of operation, an
initial microalgae mixed consortia dominated by the green microalgae
Scenedesmus sp. was cultivated in two different photo-sequencing batch reactors
operated at hydraulic retention time (HRT) of 6 days (PSBR6) and 4 days (PSBR4)
at a theoretical solids retention time (SRT) of 10 d. Both reactors were
compared with a semi-continuous reactor (SC10) operated at 10 d of HRT and 10
days of SRT (used as a control). The results indicated that PSBR6 and PSBR4
decreased Scenedesmus sp. population by 88% and 48%, respectively. However,
only PSBR6 provided suitable conditions to select cyanobacteria from an initial
green algae dominated culture. These conditions included volumetric loads of
11.72 mg TN L-1d-1, 2.04 mg TP L-1d-1 and 53.31 mg TOC L-1d-1. The remaining
nutrients in the culture led also to a phosphorus limiting N:P ratio (34:1)
that improved the increase of cyanobacteria from an initial 2% until 70% of the
total population. In addition, PSBR6 reached a biomass production of 0.12 g
L-1d-1, while removing TN, TP and TOC by 58%, 83% and 85%, respectively.
Conversely, the application of higher nutrients loads caused by lower HRT
(PSBR4) led to an increase of only 13% of cyanobacteria while SC10 remained
with the same biomass composition during all the experimental time. Thus, this
study showed that the dominance of cyanobacteria in microalgal based wastewater
treatment systems can be controlled by the operational and nutritional
conditions. This knowledge could contribute to control microalgae contamination
from up-scaling cyanobacterial biomass production in wastewater treatment
systems.


  Biological organisms exhibit diverse strategies for adapting to varying
environments. For example, a population of organisms may express the same
phenotype in all environments (`unvarying strategy'), or follow environmental
cues and express alternative phenotypes to match the environment (`tracking
strategy'), or diversify into coexisting phenotypes to cope with environmental
uncertainty (`bet-hedging strategy'). We introduce a general framework for
studying how organisms respond to environmental variations, which models an
adaptation strategy by an abstract mapping from environmental cues to
phenotypic traits. Depending on the accuracy of environmental cues and the
strength of natural selection, we find different adaptation strategies
represented by mappings that maximize the longterm growth rate of a population.
The previously studied strategies emerge as special cases of our model: the
tracking strategy is favorable when environmental cues are accurate, whereas
when cues are noisy, organisms can either use an unvarying strategy or,
remarkably, use the uninformative cue as a source of randomness to bet-hedge.
Our model of the environment-to-phenotype mapping is based on a network with
hidden units; the performance of the strategies is shown to rely on having a
high-dimensional internal representation, which can even be random.


  We investigate the global behaviour of a SIRI epidemic model with distributed
delay and relapse. From the theory of functional differential equations with
delay, we prove that the solution of the system is unique, bounded, and
positive, for all time. The basic reproduction number $R_{0}$ for the model is
computed. By means of the direct Lyapunov method and LaSalle invariance
principle, we prove that the disease free equilibrium is globally
asymptotically stable when $R_{0} < 1$. Moreover, we show that there is a
unique endemic equilibrium, which is globally asymptotically stable, when
$R_{0} > 1$.


  In a language corpus, the probability that a word occurs $n$ times is often
proportional to $1/n^2$. Assigning rank, $s$, to words according to their
abundance, $\log s$ vs $\log n$ typically has a slope of minus one. That simple
Zipf's law pattern also arises in the population sizes of cities, the sizes of
corporations, and other patterns of abundance. By contrast, for the abundances
of different biological species, the probability of a population of size $n$ is
typically proportional to $1/n$, declining exponentially for larger $n$, the
log series pattern. This article shows that the differing patterns of Zipf's
law and the log series arise as the opposing endpoints of a more general
theory. The general theory follows from the generic form of all probability
patterns as a consequence of conserved average values and the associated
invariances of scale. To understand the common patterns of abundance, the
generic form of probability distributions plus the conserved average abundance
is sufficient. The general theory includes cases that are between the Zipf and
log series endpoints, providing a broad framework for analyzing widely observed
abundance patterns.


  The intricate three-dimensional geometries of protein tertiary structures
underlie protein function and emerge through a folding process from
one-dimensional chains of amino acids. The exact spatial sequence and
configuration of amino acids, the biochemical environment and the temporal
sequence of distinct interactions yield a complex folding process that cannot
yet be easily tracked for all proteins. To gain qualitative insights into the
fundamental mechanisms behind the folding dynamics and generic features of the
folded structure, we propose a simple model of structure formation that takes
into account only fundamental geometric constraints and otherwise assumes
randomly paired connections. We find that despite its simplicity, the model
results in a network ensemble consistent with key overall features of the
ensemble of Protein Residue Networks we obtained from more than 1000 biological
protein geometries as available through the Protein Data Base. Specifically,
the distribution of the number of interaction neighbors a unit (amino acid)
has, the scaling of the structure's spatial extent with chain length, the
eigenvalue spectrum and the scaling of the smallest relaxation time with chain
length are all consistent between model and real proteins. These results
indicate that geometric constraints alone may already account for a number of
generic features of protein tertiary structures.


  The urban transition, the increased ratio of urban to rural population
globally and within countries, is a hallmark of the 21st century. Our analysis
of publicly available data from the World Bank spanning several decades for
~195 countries show that across and within nations over time, per capita Gross
Domestic Product (GDP), energy use, and CO2 emissions are lowest in
predominantly rural countries (rural > urban pop.), increase rapidly across
urbanizing countries (rural urban pop.) and are highest in the most urban
countries (rural < urban pop.). These trends coincide with changes in
employment by sector and gender. Rural economies are based largely on
employment in the resource-extraction sector, which includes agriculture,
fisheries, forestry, and mining. In urbanizing nations, male employment is
predominantly in the industrial sector, including public utilities, while
female employment is higher in service-based than resource-based economies. In
the most urban nations, service economies predominate with some countries
employing 90% of women and 65% of men in the service sector. Our analysis shows
that per capita GDP, energy use, and CO2 emissions increase by over two orders
of magnitude from low-income, resource-based rural countries to high-income,
urbanized countries with predominantly service economies. Data from the U.S.
over the past 200 years illuminate a socio-metabolic urban transition similar
to that seen globally in recent decades across countries and through time. Our
study suggests that increased energy demand and climate consequences of burning
fossil fuels will continue to accompany a rapidly urbanizing planet posing
major challenges for global sustainability.


  Rule-based modelling allows to represent molecular interactions in a compact
and natural way. The underlying molecular dynamics, by the laws of stochastic
chemical kinetics, behaves as a continuous-time Markov chain. However, this
Markov chain enumerates all possible reaction mixtures, rendering the analysis
of the chain computationally demanding and often prohibitive in practice. We
here describe how it is possible to efficiently find a smaller, aggregate
chain, which preserves certain properties of the original one. Formal methods
and lumpability notions are used to define algorithms for automated and
efficient construction of such smaller chains (without ever constructing the
original ones). We here illustrate the method on an example and we discuss the
applicability of the method in the context of modelling large signalling
pathways.


  Motivation: With the development of droplet based systems, massive single
cell transcriptome data has become available, which enables analysis of
cellular and molecular processes at single cell resolution and is instrumental
to understanding many biological processes. While state-of-the-art clustering
methods have been applied to the data, they face challenges in the following
aspects: (1) the clustering quality still needs to be improved; (2) most models
need prior knowledge on number of clusters, which is not always available; (3)
there is a demand for faster computational speed. Results: We propose to tackle
these challenges with Parallel Split Merge Sampling on Dirichlet Process
Mixture Model (the Para-DPMM model). Unlike classic DPMM methods that perform
sampling on each single data point, the split merge mechanism samples on the
cluster level, which significantly improves convergence and optimality of the
result. The model is highly parallelized and can utilize the computing power of
high performance computing (HPC) clusters, enabling massive clustering on huge
datasets. Experiment results show the model outperforms current widely used
models in both clustering quality and computational speed. Availability: Source
code is publicly available on
https://github.com/tiehangd/Para_DPMM/tree/master/Para_DPMM_package


  Pro-social punishment and exclusion are common means to elevate the level of
cooperation among unrelated individuals. Indeed, it is worth pointing out that
the combined use of these two strategies is quite common across human
societies. However, it is still not known how a combined strategy where
punishment and exclusion are switched can promote cooperation from the
theoretical perspective. In this paper, we thus propose two different switching
strategies, namely peer switching that is based on peer punishment and peer
exclusion, and pool switching that is based on pool punishment and pool
exclusion. Individuals adopting the switching strategy will punish defectors
when their numbers are below a threshold and exclude them otherwise. We study
how the two switching strategies influence the evolutionary dynamics in the
public goods game. We show that an intermediate value of the threshold leads to
a stable coexistence of cooperators, defectors and players adopting the
switching strategy in a well-mixed population, and this regardless of whether
the pool-based or the peer-based switching strategy is introduced. Moreover, we
show that the pure exclusion strategy alone is able to evoke a limit cycle
attractor in the evolutionary dynamics, such that cooperation can coexist with
other strategies.


  Totally asymmetric simple exclusion process (TASEP) was originally introduced
as a model for the traffic-like collective movement of ribosomes on a messenger
RNA (mRNA) that serves as the track for the motor-like forward stepping of
individual ribosomes. In each step, a ribosome elongates a protein by a single
unit using the track also as a template for protein synthesis. But,
pre-fabricated, functionally competent, ribosomes are not available to begin
synthesis of protein; a subunit directionally scans the mRNA in search of the
pre-designated site where it is supposed to bind with the other subunit and
begin the synthesis of the corresponding protein. However, because of `leaky'
scanning, a fraction of the scanning subunits miss the target site and continue
their search beyond the first target. Sometimes such scanners successfully
identify the site that marks the site for initiation of the synthesis of a
different protein. In this paper, we develop an exclusion model, with three
interconvertible species of hard rods, to capture some of the key features of
these biological phenomena and study the effects of the interference of the
flow of the different species of rods on the same lattice. More specifically,
we identify the meantime for the initiation of protein synthesis as appropriate
mean {\it first-passage} time that we calculate analytically using the
formalism of backward master equations. In spite of the approximations made,
our analytical predictions are in reasonably good agreement with the numerical
data that we obtain by performing Monte Carlo simulations. We also compare our
results with a few experimental facts reported in the literature and propose
new experiments for testing some of our new quantitative predictions.


  Over the last two decades epidermal growth factor receptor (EGFR) kinase has
become an important target to treat non-small cell lung cancer (NSCLC).
Currently, three generations of EGFR kinase-targeted small molecule drugs have
been FDA approved. They nominally produce a response at the start of treatment
and lead to a substantial survival benefit for patients. However, long-term
treatment results in acquired drug resistance and further vulnerability to
NSCLC. Therefore, novel EGFR kinase inhibitors that specially overcome acquired
mutations are urgently needed. To this end, we carried out a comprehensive
study of different EGFR kinase mutants using a structural systems pharmacology
strategy. Our analysis shows that both wild-type and mutated structures exhibit
multiple conformational states that have not been observed in solved crystal
structures. We show that this conformational flexibility accommodates diverse
types of ligands with multiple types of binding modes. These results provide
insights for designing a new-generation of EGFR kinase inhibitor that combats
acquired drug-resistant mutations through a multi-conformation-based drug
design strategy.


  We investigate evolutionary dynamics of altruism with long-range interaction
on a cycle. The interaction between individuals is described by a simplified
version of the prisoner's dilemma (PD) game in which the payoffs are
parameterized by $c$, the cost of a cooperative action. In our model, the
probabilities of the game interaction and competition decay algebraically with
$r_{AB}$, the distance between two players $A$ and $B$, but with different
exponents: That is, the probability to play the PD game is proportional to
$r_{AB}^{-\alpha}$. If player $A$ is chosen for death, on the other hand, the
probability for $B$ to occupy the empty site is proportional to
$r_{AB}^{-\beta}$. In a limiting case of $\beta\to\infty$, where the
competition for an empty site occurs between its nearest neighbors only, we
analytically find the condition for the proliferation of altruism in terms of
$c_{th}$, a threshold of $c$ below which altruism prevails. For finite $\beta$,
we conjecture a formula for $c_{th}$ as a function of $\alpha$ and $\beta$. We
also propose a numerical method to locate $c_{th}$, according to which we
observe excellent agreement with the conjecture even when the selection
strength is of considerable magnitude.


  In this paper we consider an SIRS epidemic model under a general assumption
of density-dependent mortality. We prove the global stability of the
disease-free equilibrium and propose a Lyapunov function that allows to
demonstrate the global stability of the (unique) endemic state under broad
conditions.


  Direct comparison of three-dimensional (3D) objects is computationally
expensive due to the need for translation, rotation, and scaling of the objects
to evaluate their similarity. In applications of 3D object comparison, often
identifying specific local regions of objects is of particular interest. We
have recently developed a set of 2D moment invariants based on discrete
orthogonal Krawtchouk polynomials for comparison of local image patches. In
this work, we extend them to 3D and construct 3D Krawtchouk descriptors (3DKD)
that are invariant under translation, rotation, and scaling. The new
descriptors have the ability to extract local features of a 3D surface from any
region-of-interest. This property enables comparison of two arbitrary local
surface regions from different 3D objects. We present the new formulation of
3DKD and apply it to the local shape comparison of protein surfaces in order to
predict ligand molecules that bind to query proteins.


  In biodiversity and ecosystem functioning (BEF) research, the Loreau-Hector
(LH) statistical scheme is widely-used to partition the effect of biodiversity
on ecosystem properties into a "complementarity effect" and a "selection
effect". This selection effect was originally considered analogous to the
selection term in the Price equation from evolutionary biology. However, a key
paper published over thirteen years ago challenged this interpretation by
devising a new tripartite partitioning scheme that purportedly quantified the
role of selection in biodiversity experiments more accurately. This tripartite
method, as well as its recent spatiotemporal extension, were both developed as
an attempt to apply the Price equation in a BEF context. Here, we demonstrate
that the derivation of this tripartite method, as well as its spatiotemporal
extension, involve a set of incoherent and nonsensical mathematical arguments
driven largely by na\"ive visual analogies with the original Price equation,
that result in neither partitioning scheme quantifying any real property in the
natural world. Furthermore, we show that Loreau and Hector's original selection
effect always represented a true analog of the original Price selection term,
making the tripartite partitioning scheme a nonsensical solution to a
non-existent problem [...]


  Overpopulation and environmental degradation due to inadequate resource-use
are outcomes of human's ecosystem engineering that has profoundly modified the
world's landscape. Despite the age-old concern that unchecked population and
economic growth may be unsustainable, the prospect of societal collapse remains
contentious today. Contrasting with the usual approach to modeling human-nature
interactions, which are based on the Lotka-Volterra predator-prey model with
humans as the predators and nature as the prey, here we address this issue
using a discrete-time population dynamics model of ecosystem engineers. The
growth of the population of engineers is modeled by the Beverton-Holt equation
with a density-dependent carrying capacity that is proportional to the number
of usable habitats. These habitats (e.g., farms) are the products of the work
of the individuals on the virgin habitats (e.g., native forests), hence the
denomination engineers of ecosystems to those agents. The human-made habitats
decay into degraded habitats, which eventually regenerate into virgin habitats.
For slow regeneration resources, we find that the dynamics is dominated by
cycles of prosperity and collapse, in which the population reaches vanishing
small densities. However, increase of the efficiency of the engineers to
explore the resources eliminates the dangerous cyclical patterns of feast and
famine and leads to a stable equilibrium that balances population growth and
resource availability. This finding supports the viewpoint of growth optimists
that technological progress may avoid collapse.


  Ecological management and decision-making typically focus on uncertainty
about the future, but surprisingly little is known about how to account for
uncertainty of the present: that is, the realities of having only partial or
imperfect measurements. Our primary paradigms for handling decisions under
uncertainty -- the precautionary principle and optimal control -- have so far
given contradictory results. This paradox is best illustrated in the example of
fisheries management, where many ideas that guide thinking about ecological
decision making were first developed. We find that simplistic optimal control
approaches have repeatedly concluded that a manager should increase catch
quotas when faced with greater uncertainty about the fish biomass. Current best
practices take a more precautionary approach, decreasing catch quotas by a
fixed amount to account for uncertainty. Using comparisons to both simulated
and historical catch data, we find that neither approach is sufficient to avoid
stock collapses under moderate observational uncertainty. Using partially
observed Markov decision process (POMDP) methods, we demonstrate how this
paradox arises from flaws in the standard theory, which contributes to
over-exploitation of fisheries and increased probability of economic and
ecological collapse. In contrast, we find POMDP-based management avoids such
over-exploitation while also generating higher economic value. These results
have significant implications for how we handle uncertainty in both fisheries
and ecological management more generally.


  Reaction currents in chemical networks usually increase when increasing their
driving affinities. But far from equilibrium the opposite can also happen. We
find that such negative differential response (NDR) occurs in reaction schemes
of major biological relevance, namely, substrate inhibition and autocatalysis.
We do so by deriving the full counting statistics of two minimal representative
models using large deviation methods. We argue that NDR implies the existence
of optimal affinities that maximize the robustness against environmental and
intrinsic noise at intermediate values of dissipation. An analogous behavior is
found in dissipative self-assembly, for which we identify the optimal working
conditions set by NDR.


  Many instances of cellular signaling and transcriptional regulation involve
switch-like molecular responses to the presence or absence of input ligands. To
understand how these responses come about and how they can be harnessed, we
develop a statistical mechanical model to characterize the types of Boolean
logic that can arise from allosteric molecules following the
Monod-Wyman-Changeux (MWC) model. Building upon previous work, we show how an
allosteric molecule regulated by two inputs can elicit AND, OR, NAND and NOR
responses, but is unable to realize XOR or XNOR gates. Next, we demonstrate the
ability of an MWC molecule to perform ratiometric sensing - a response behavior
where activity depends monotonically on the ratio of ligand concentrations. We
then extend our analysis to more general schemes of combinatorial control
involving either additional binding sites for the two ligands or an additional
third ligand and show how these additions can cause a switch in the logic
behavior of the molecule. Overall, our results demonstrate the wide variety of
control schemes that biological systems can implement using simple mechanisms.


  The early fly embryo offers a relatively pure version of the problem of
spatial scaling in biological pattern formation. Within three hours, a
"blueprint" for the final segmented body plan of the animal is visible in
striped patterns of gene expression. We measure the positions of these stripes
in an ensemble of 100+ embryos from a laboratory strain of Drosophila
melanogaster, under controlled conditions. These embryos vary in length by only
4% (rms), yet stripes are positioned with 1% accuracy; precision and scaling of
the pattern are intertwined. We can see directly the variation of absolute
stripe positions with length, and the precision is so high as to exclude
alternatives, such as combinations of unscaled signals from the two ends of the
embryo.


  Biochemical molecules interact through modification and binding reactions,
giving raise to a combinatorial number of possible biochemical species. The
time-dependent evolution of concentrations of the species is commonly described
by a system of coupled ordinary differential equations (ODEs). However, the
analysis of such high-dimensional, non-linear system of equations is often
computationally expensive and even prohibitive in practice. The major challenge
towards reducing such models is providing the guarantees as to how the solution
of the reduced model relates to that of the original model, while avoiding to
solve the original model. In this paper, we have designed and tested an
approximation method for ODE models of biochemical reaction systems, in which
the guarantees are our major requirement. Borrowing from tropical analysis
techniques, dominance relations among terms of each species' ODE are exploited
to simplify the original model, by neglecting the dominated terms. As the
dominant subsystems can change during the system's dynamics, depending on which
species dominate the others, several possible modes exist. Thus, simpler models
consisting of only the dominant subsystems can be assembled into hybrid,
piecewise smooth models, which approximate the behavior of the initial system.
By combining the detection of dominated terms with symbolic bounds propagation,
we show how to approximate the original model by an assembly of simpler models,
consisting in ordinary differential equations that provide time-dependent lower
and upper bounds for the concentrations of the initial models species. Our
method provides sound interval bounds for the concentrations of the chemical
species, and hence can serve to evaluate the faithfulness of
tropicalization-based reduction heuristics for ODE models of biochemical
reduction systems. The method is tested on several case studies.


  The performance of most error-correction algorithms that operate on genomic
sequencer reads is dependent on the proper choice of its configuration
parameters, such as the value of k in k-mer based techniques. In this work, we
target the problem of finding the best values of these configuration parameters
to optimize error correction. We perform this in a data-driven manner, due to
the observation that different configuration parameters are optimal for
different datasets, i.e., from different instruments and organisms. We use
language modeling techniques from the Natural Language Processing (NLP) domain
in our algorithmic suite, Athena, to automatically tune the
performance-sensitive configuration parameters. Through the use of N-Gram and
Recurrent Neural Network (RNN) language modeling, we validate the intuition
that the EC performance can be computed quantitatively and efficiently using
the perplexity metric, prevalent in NLP. After training the language model, we
show that the perplexity metric calculated for runtime data has a strong
negative correlation with the correction of the erroneous NGS reads. Therefore,
we use the perplexity metric to guide a hill climbing-based search, converging
toward the best $k$-value. Our approach is suitable for both de novo and
comparative sequencing (resequencing), eliminating the need for a reference
genome to serve as the ground truth. This is important because the use of a
reference genome often carries forward the biases along the stages of the
pipeline.


  It is tempting to believe that we now own the genome. The ability to read and
re-write it at will has ushered in a stunning period in the history of science.
Nonetheless, there is an Achilles heel exposed by all of the genomic data that
has accrued: we still don't know how to interpret it. Many genes are subject to
sophisticated programs of transcriptional regulation, mediated by DNA sequences
that harbor binding sites for transcription factors which can up- or
down-regulate gene expression depending upon environmental conditions. This
gives rise to an input-output function describing how the level of expression
depends upon the parameters of the regulated gene { for instance, on the number
and type of binding sites in its regulatory sequence. In recent years, the
ability to make precision measurements of expression, coupled with the ability
to make increasingly sophisticated theoretical predictions, have enabled an
explicit dialogue between theory and experiment that holds the promise of
covering this genomic Achilles heel. The goal is to reach a predictive
understanding of transcriptional regulation that makes it possible to calculate
gene expression levels from DNA regulatory sequence. This review focuses on the
canonical simple repression motif to ask how well the models that have been
used to characterize it actually work. We consider a hierarchy of increasingly
sophisticated experiments in which the minimal parameter set learned at one
level is applied to make quantitative predictions at the next. We show that
these careful quantitative dissections provide a template for a predictive
understanding of the many more complex regulatory arrangements found across all
domains of life.


  Base pair mismatch can relieve mechanical stress in highly strained DNA
molecules, but how it affects their kinetic stability is not known. Using
single-molecule Fluorescence Resonance Energy Transfer (FRET), we measured the
lifetimes of tightly bent DNA loops with and without base pair mismatch.
Surprisingly, for loops captured by stackable sticky ends, the mismatch
decreased the loop lifetime despite reducing the overall bending stress, and
the decrease was largest when the mismatch was placed at the DNA midpoint.
These findings show that base pair mismatch transfers bending stress to the
opposite side of the loop through an allosteric mechanism known as cooperative
kinking. Based on this mechanism, we present a three-state model that explains
the apparent dichotomy between thermodynamic and kinetic stability of DNA
loops.


  The human adaptive immune response is known to weaken in advanced age,
resulting in increased severity of pathogen-born illness, poor vaccine
efficacy, and a higher prevalence of cancer in the elderly. Age-related erosion
of the T-cell compartment has been implicated as a likely cause, but the
underlying mechanisms driving this immunosenescence have not been
quantitatively modeled and systematically analyzed. T-cell receptor diversity,
or the extent of pathogen-derived antigen responsiveness of the T-cell pool, is
known to diminish with age, but inherent experimental difficulties preclude
accurate analysis on the full organismal level. In this paper, we formulate a
mechanistic mathematical model of T-cell population dynamics on the
immunoclonal subpopulation level, which provides quantitative estimates of
diversity. We define different estimates for diversity that depend on the
individual number of cells in a specific immunoclone. We show that diversity
decreases with age primarily due to diminished thymic output of new T-cells and
the resulting overall loss of small immunoclones.


  Motivated by the classical Susceptible-Infected-Recovered (SIR) epidemic
models proposed by Kermack and Mckendrick, we consider a class of stochastic
compartmental dynamical systems with a notion of partial ordering among the
compartments. We call such systems unidirectional Mass Transfer Models (MTMs).
We show that there is a natural way of interpreting a uni-directional MTM as a
Survival Dynamical System (SDS) that is described in terms of survival
functions instead of population counts. This SDS interpretation allows us to
employ tools from survival analysis to address various issues with data
collection and statistical inference of unidirectional MTMs. In particular, we
propose and numerically validate a statistical inference procedure based on
SDS-likelihoods. We use the SIR model as a running example throughout the paper
to illustrate the ideas.


  Living systems regulate many aspects of their behavior through periodic
oscillations of molecular concentrations, which function as `biochemical
clocks.' These clocks are intrinsically subject to thermal fluctuations, so
that the duration of a full oscillation cycle is random. Their success in
carrying out their biological function is thought to depend on the degree to
which these fluctuations in the cycle period can be suppressed. Biochemical
oscillators also require a constant supply of free energy in order to break
detailed balance and maintain their cyclic dynamics. For a given free energy
budget, the recently discovered `thermodynamic uncertainty relation' yields the
magnitude of period fluctuations in the most precise conceivable free-running
clock. In this paper, we show that computational models of real biochemical
clocks severely underperform this optimum, with fluctuations several orders of
magnitude larger than the theoretical minimum. We argue that this suboptimal
performance is due to the small number of internal states per molecule in these
models, combined with the high level of thermodynamic force required to
maintain the system in the oscillatory phase. We introduce a new model with a
tunable number of internal states per molecule, and confirm that it approaches
the optimal precision as this number increases.


  Communication signals often comprise an array of colors, lines, spots, notes
or odors that are arranged in complex patterns, melodies or blends. Receiver
perception is assumed to influence preference and thus the evolution of signal
design, but evolutionary biologists still struggle to understand how
perception, preference, and signal design are mechanistically linked. In
parallel, the field of empirical aesthetics aims to understand why people like
some designs more than others. The model of processing bias discussed here is
rooted in empirical aesthetics, which posits that preferences are influenced by
the emotional system as it monitors the dynamics of information processing, and
that attractive signals have either effective designs that maximize information
transmission, efficient designs that allow information processing at low
metabolic cost, or both. We refer to the causal link between preference and the
emotionally rewarding experience of effective and efficient information
processing as the processing bias, and we apply it to the evolutionary model of
sensory drive. A sensory drive model that incorporates processing bias
hypothesizes a causal chain of relationships between the environment,
perception, pleasure, preference, and ultimately the evolution of signal
design, from simple to complex.


  The more we learn about the cytoplasm of cells, the more we realise that the
cytoplasm is not uniform but instead is highly inhomogeneous. In any
inhomogeneous solution, there are concentration gradients, and particles move
either up or down these gradients due to a mechanism called diffusiophoresis. I
estimate that inside metabolically active cells, the dynamics of particles can
be strongly accelerated by diffusiophoresis, provided that they are at least
tens of nanometres across. The dynamics of smaller objects, such as single
proteins are largely unaffected.


  Deep Learning (DL) algorithms hold great promise for applications in the
field of computational biophysics. In fact, the vast amount of available
molecular structures, as well as their notable complexity, constitutes an ideal
context in which DL-based approaches can be profitably employed. To express the
full potential of these techniques, though, it is a prerequisite to express the
information contained in the molecule's atomic positions and distances in a set
of input quantities that the network can process. Many of the molecular
descriptors devised insofar are effective and manageable for relatively small
structures, but become complex and cumbersome for larger ones. Furthermore,
most of them are defined locally, a feature that could represent a limit for
those applications where global properties are of interest. Here, we build a
deep learning architecture capable of predicting non-trivial and intrinsically
global quantities, that is, the eigenvalues of a protein's lowest-energy
fluctuation modes. This application represents a first, relatively simple test
bed for the development of a neural network approach to the quantitative
analysis of protein structures, and demonstrates unexpected use in the
identification of mechanically relevant regions of the molecule.


  Exploring and understanding the protein-folding problem has been a
long-standing challenge in molecular biology. Here, using molecular dynamics
simulation, we reveal how parallel distributed adjacent planar peptide groups
of unfolded proteins fold reproducibly following explicit physical folding
codes in aqueous environments due to electrostatic attractions. Superfast
folding of protein is found to be powered by the contribution of the formation
of hydrogen bonds. Temperature-induced torsional waves propagating along
unfolded proteins break the parallel distributed state of specific amino acids,
inferred as the beginning of folding. Electric charge and rotational resistance
differences among neighboring side-chains are used to decipher the physical
folding codes by means of which precise secondary structures develop. We
present a powerful method of decoding amino acid sequences to predict native
structures of proteins. The method is verified by comparing the results
available from experiments in the literature.


  Proteins are essential for maintaining life. For example, knowing the
structure of a protein, cell regulatory mechanisms of organisms can be modeled,
supporting the development of disease treatments or the understanding of
relationships between protein structures and food attributes. However,
discovering the structure of a protein can be a difficult and expensive task,
since it is hard to explore the large search to predict even a small protein.
Template-based methods (coarse-grained, homology, threading etc) depend on
Prior Knowledge (PK) of proteins determined using other methods as X-Ray
Crystallography or Nuclear Magnetic Resonance. On the other hand, template-free
methods (full-atom and ab initio) rely on atoms physical-chemical properties to
predict protein structures. In comparison with other approaches, the Estimation
of Distribution Algorithms (EDAs) can require significant less PK, suggesting
that it could be adequate for proteins of low-level of PK. Finding an EDA able
to handle both prediction quality and computational time is a difficult task,
since they are strong inversely correlated. We developed an EDA specific for
the ab initio Protein Structure Prediction (PSP) problem using full-atom
representation. We developed one univariate and two bivariate probabilistic
models in order to design a proper EDA for PSP. The bivariate models make
relationships between dihedral angles $\phi$ and $\psi$ within an amino acid.
Furthermore, we compared the proposed EDA with other approaches from the
literature. We noticed that even a relatively simple algorithm such as Random
Walk can find the correct solution, but it would require a large amount of
prior knowledge (biased prediction). On the other hand, our EDA was able to
correctly predict with no prior knowledge at all, characterizing such a
prediction as pure ab initio.


  We translate a coagulation-framentation model, describing the dynamics of
animal group size distributions, into a model for the population distribution
and associate the \blue{nonlinear} evolution equation with a Markov jump
process of a type introduced in classic work of H.~McKean. In particular this
formalizes a model suggested by H.-S. Niwa [J.~Theo.~Biol.~224 (2003)] with
simple coagulation and fragmentation rates. Based on the jump process, we
develop a numerical scheme that allows us to approximate the equilibrium for
the Niwa model, validated by comparison to analytical results by Degond et al.
[J.~Nonlinear Sci.~27 (2017)], and study the population and size distributions
for more complicated rates. Furthermore, the simulations are used to describe
statistical properties of the underlying jump process. We additionally discuss
the relation of the jump process to models expressed in stochastic differential
equations and demonstrate that such a connection is justified in the case of
nearest-neighbour interactions, as opposed to global interactions as in the
Niwa model.


  Bacterial biofilms, surface-attached communities of cells, are in some
respects similar to colloidal solids; both are densely packed with non-zero
yield stresses. However, unlike non-living materials, bacteria reproduce and
die, breaking mechanical equilibrium and inducing collective dynamic responses.
We report experiments and theory investigating the motion of immotile Vibrio
cholerae, which can kill each other and reproduce in biofilms. We vary
viscosity by using bacterial variants that secrete different amounts of
extracellular matrix polymers, but are otherwise identical. Unlike
thermally-driven diffusion, in which diffusivity decreases with increased
viscosity, we find that cellular motion mediated by death and reproduction is
independent of viscosity over timescales relevant to bacterial reproduction. To
understand this surprising result, we use two separate modeling approaches.
First we perform explicitly mechanical simulations of one-dimensional chains of
Voigt-Kelvin elements that can die and reproduce. Next, we perform an
independent statistical approach, modeling Brownian motion with the classic
Langevin equation under an effective temperature that depends on cellular
division rate. The diffusion of cells in both approaches agrees quite well,
supporting a kinetic interpretation for the effective temperature used here and
developed in previous work. As the viscoelastic behavior of biofilms is
believed to play a large role in their anomalous biological properties, such as
antibiotic resistance, the independence of cellular diffusive motion ---
important for biofilm growth and remodeling --- on viscoelastic properties
likely holds ecological, medical, and industrial relevance.


  Life forms exhibit such a degree of exquisite organization that it seems
impossible that they could have developed out of a process of trial and error,
as intimated by the theory of Darwinian evolution. In this general public paper
I discuss how differential reproduction rates work in producing an exceedingly
high degree of improbability, and the conceptual tools of the theory of
evolution help us to predict, to some degree, the course of evolution -- as it
is routinely done, e.g., in the process leading to the yearly influenza
vaccines.


  The Red Queen Hypothesis, which suggests that continual evolution can result
from solely biotic interactions, has been studied in macroevolutionary and
microevolutionary contexts. While the latter has been effective in describing
examples in which evolution does not cease, describing which properties lead to
continual evolution or to stasis remains a major challenge. In many contexts it
is unclear which assumptions are necessary for continual evolution, and whether
described behavior is robust under perturbations. Our aim here is to prove
continual evolution under minimal conditions and in a general framework, thus
automatically obtaining robustness. We show that the combination of a fast
positive and a slow negative feedback causes continual evolution with a single
evolving trait, provided the ecological timescale is sufficiently separated
from the timescales of mutations and negative feedback. Our approach and
results form a next step towards a deeper understanding of the evolutionary
dynamics resulting from biotic interactions.


  Time is arguably the key limiting factor for interstellar exploration. At
high speeds, flyby missions to nearby stars by laser propelled wafersats taking
50-100 years would be feasible. Directed energy launch systems could accelerate
on the other side also crafts weighing several tons to cruising speeds of the
order of 1000\,km/s (c/300). At these speeds, superconducting magnetic sails
would be able to decelerate the craft by transferring kinetic energy to the
protons of the interstellar medium. A tantalizing perspective, which would
allow interstellar probes to stop whenever time is not a limiting factor. Prime
candidates are in this respect Genesis probes, that is missions aiming to offer
terrestrial life new evolutionary pathways on potentially habitable but
hitherto barren exoplanets.
  Genesis missions raise important ethical issues, in particular with regard to
planetary protection. Here we argue that exoplanetary and planetary protection
differ qualitatively as a result of the vastly different cruising times for
payload delivering probes, which are of the order of millennia for interstellar
probes, but only of years for solar system bodies. Furthermore we point out
that our galaxy may harbor a large number of habitable exoplanets, M-dwarf
planets, which could be sterile due to the presence of massive primordial
oxygen atmospheres. We believe that the prospect terrestrial life has in our
galaxy would shift on a fundamental level in case that the existence of this
type of habitable but sterile oxygen planets will be corroborated by future
research. It may also explain why our sun is not a M dwarf, the most common
star type, but a medium-sized G-class star.


  Understanding whether preferences are sensitive to the frame has been a major
topic of debate in the last decades. For example, several works have explored
whether the dictator game in the give frame gives rise to a different rate of
pro-sociality than the same game in the take frame, leading to mixed results.
Here we contribute to this debate with two experiments. In Study 1 ($N=567$) we
implement an extreme dictator game in which the dictator either gets \$0.50 and
the recipient gets nothing, or the opposite (i.e., the recipient gets \$0.50
and the dictator gets nothing). We experimentally manipulate the words
describing the available actions using six terms, from very negative (e.g.,
stealing) to very positive (e.g., donating) connotations. We find that the rate
of pro-sociality is affected by the words used to describe the available
actions. In Study 2 ($N=221$) we ask brand new participants to rate each of the
words used in Study 1 from ``extremely wrong'' to ``extremely right'' . We find
that these moral judgments explain the framing effect in Study 1. In sum, our
studies provide evidence that framing effects in an extreme Dictator game can
be generated using morally loaded language.


  The V(D)J recombination forms the immunoglobulin genes by joining the
variable (V), diversity (D), and joining (J) germline genes. Since variations
in germline genes have been linked to various diseases, personalized
immunogenomics aims at finding alleles of germline genes across various
patients. Although recent studies described algorithms for de novo inference of
V and J genes from immunosequencing data, they stopped short of solving a more
difficult problem of reconstructing D genes that form the highly divergent CDR3
regions and provide the most important contribution to the antigen binding. We
present the IgScout algorithm for de novo D gene reconstruction and apply it to
reveal new alleles of human D genes and previously unknown D genes in camel, an
important model organism in immunology. We further analyze non-canonical V(DD)J
recombination that results in unusually long tandem CDR3s and thus expands the
diversity of the antibody repertoires. We demonstrate that tandem CDR3s
represent a consistent and functional feature of all analyzed immunosequencing
datasets, reveal ultra-long tandem CDR3s, and shed light on the mechanism
responsible for their formation.


  Since the largest 2014-2016 Ebola virus disease outbreak in West Africa,
understanding of Ebola virus infection has improved, notably the involvement of
innate immune mediators. Amongst them, collectins are important players in the
antiviral innate immune defense. A screening of Ebola glycoprotein
(GP)-collectins interactions revealed the specific interaction of human
surfactant protein D (hSP-D), a lectin expressed in lung and liver, two
compartments where Ebola was found in vivo. Further analyses have demonstrated
an involvement of hSP-D in the enhancement of virus infection in several in
vitro models. Similar effects were observed for porcine SP-D (pSP-D). In
addition, both hSP-D and pSP-D interacted with Reston virus (RESTV) GP and
enhanced pseudoviral infection in pulmonary cells. Thus, our study reveals a
novel partner of Ebola GP that may participate to enhance viral spread.


  In this paper, we introduce a general framework for co-infection as
cooperative SIR dynamics. We first solve analytically CGCG model [1] and then
the generalized model in symmetric scenarios. We calculate transition points,
order parameter, i.e. total number of infected hosts. Also we show analytically
there is a saddle-node bifurcation for two cooperative SIR dynamics and the
transition is hybrid. Moreover, we investigate where symmetric solution is
stable for initial fluctuations. Then we study asymmetric cases of parameters.
The more asymmetry, for the primary and secondary infection rates of one
pathogen in comparison to the other pathogen, can lead to the less infected
hosts, the higher epidemic threshold and continuous transitions. Our model and
results for co-infection in combination with super-infection [2] can open a
road to model disease ecology.


  Predator-prey networks originating from different aqueous and terrestrial
environments are compared to assess if the difference in environments of these
networks produce any significant difference in the structure of such
predator-prey networks. Spectral graph theory is used firstly to discriminate
between the structure of such predator-prey networks originating from aqueous
and terrestrial environments and secondly to establish that the difference
observed in the structure of networks originating from these two environments
are precisely due to the way edges are oriented in these networks and are not a
property of random networks.We use random projections in $\mathbb{R^2}$ and
$\mathbb{R^3}$ of weighted spectral distribution (WSD) of the networks
belonging to the two classes viz. aqueous and terrestrial to differentiate
between the structure of these networks. The spectral theory of graph
non-randomness and relative non-randomness is used to establish the deviation
of structure of these networks from having a topology similar to random
networks.We thus establish the absence of a universal structural pattern across
predator-prey networks originating from different environments.


  We consider the adaptive learning rule of Harley (1981) for behavior
selection in symmetric conflict games in large populations. The rule uses
organisms' past, accumulated rewards as the predictor for the future behavior,
and can be traced in many life forms from bacteria to humans. We derive a
partial differential equation (PDE) that describes the stochastic learning in a
population of agents. The equation has simple structure of the `conservation of
mass'-type equation in the space of stimuli to engage in a particular type of
behavior. We analyze the solutions of the PDE model for typical 2x2 games. It
is found that in games with small residual stimuli, adaptive learning rules
with faster memory decay have an evolutionary advantage.


  Linear mixed models (LMMs) are widely used for heritability estimation in
genome-wide association studies (GWAS). In standard approaches to heritability
estimation with LMMs, a genetic relationship matrix (GRM) must be specified. In
GWAS, the GRM is frequently a correlation matrix estimated from the study
population's genotypes, which corresponds to a normalized Euclidean distance
kernel. In this paper, we show that reliance on the Euclidean distance kernel
contributes to several unresolved modeling inconsistencies in heritability
estimation for GWAS. These inconsistencies can cause biased heritability
estimates in the presence of linkage disequilibrium (LD), depending on the
distribution of causal variants. We show that these biases can be resolved (at
least at the modeling level) if one adopts a Mahalanobis distance-based GRM for
LMM analysis. Additionally, we propose a new definition of partitioned
heritability -- the heritability attributable to a subset of genes or single
nucleotide polymorphisms (SNPs) -- using the Mahalanobis GRM, and show that it
inherits many of the nice consistency properties identified in our original
analysis. Partitioned heritability is a relatively new area for GWAS analysis,
where inconsistency issues related to LD have previously been known to be
especially pernicious.


  A phylogenetic tree is an acyclic graph with distinctly labeled leaves, whose
internal edges have a positive weight. Given a set of n leaves, the collection
of all phylogenetic trees with this leaf set can be assembled into a metric
cube complex known as phylogenetic tree space, or Billera-Holmes-Vogtmann tree
space. In this largely combinatorial paper, we show that the isometry group of
this space is the symmetric group on n elements. This fact is relevant to
distance-based analyses of phylogenetic tree sets.


  Collapse of the polypeptide backbone is an integral part of protein folding.
Using polyglycine as a probe, we explore the nonequilibrium pathways of protein
collapse in water. We find that the collapse depends on the competition between
hydration effects and intra-peptide interactions. Once intra-peptide van der
Waal interactions dominate, the chain collapses along a nonequilibrium pathway
characterized by formation of pearl-necklace-like local clusters as
intermediates that eventually coagulate into a single globule. By describing
this coarsening through the contact probability as a function of distance along
the chain, we extract a time-dependent length scale that grows in linear
fashion. The collapse dynamics is characterized by a dynamical critical
exponent $z=0.5$ that is much smaller than the values of $z=1-2$ reported for
non-biological polymers. This difference in the exponents is explained by the
instantaneous formation of intra-chain hydrogen bonds and local ordering that
may be correlated with the observed fast folding times in proteins.


  Calcein leakage assays are a standard experimental set-up for probing the
extent of damage induced by external agents on synthetic lipid vesicles. The
fluorescence signal associated with calcein release from liposomes is the
signature of vesicle disruption, transient pore formation or vesicle fusion.
This type of assay is widely used to test the membrane disruptive effect of
biological macromolecules, such as proteins, antimicrobial peptides and RNA and
is also used on synthetic nanoparticles with a polymer, metal or oxide core.
Little is known about the effect that calcein and other fluorescent dyes may
have on the properties of lipid bilayers, potentially altering their structure
and permeability. Here we develop a coarse-grained model of calcein that is
compatible with the Martini force field for lipids. We validate the model by
comparing its dimerization free energy, aggregation behavior at different
concentrations and interaction with a
1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane to those
obtained at atomistic resolution. Our coarse-grained description of calcein
makes it suitable for the simulation of large calcein-filled liposomes and of
their interactions with external agents, allowing for a direct comparison
between simulations and experimental liposome leakage assays.


  There is great interest in predicting rare and extreme events in complex
systems, and in particular, understanding the role of network topology in
facilitating such events. In this work, we show that degree dispersion -- the
fact that the number of local connections in networks varies broadly --
increases the probability of large, rare fluctuations in population networks
generically. We perform explicit calculations for two canonical and distinct
classes of rare events: network extinction and switching. When the distance to
threshold is held constant, and hence stochastic effects are fairly compared
among networks, we show that there is a universal, exponential increase in the
rate of rare events proportional to the variance of a network's degree
distribution over its mean squared.


  Whether or not cooperation is favored over defection in evolutionary games
can be assigned by structure coefficients for any arrangement of cooperators
and defectors on any network modeled as a regular graph. We study how these
structure coefficients relate to a scaling of dilemma strength in social
dilemma games. It is shown that some graphs permit certain arrangements of
cooperators and defectors to possess particularly large structure coefficients.
Moreover, these large coefficients imply particularly large sections of a
bounded parameter plane spanned by scaling gamble-intending and risk-averting
dilemma strength.


  Different evolutionary models are known to make disparate predictions for the
success of an invading mutant in some situations. For example, some
evolutionary mechanics lead to amplification of selection in structured
populations, while others suppress it. Here, we use computer simulations to
study evolutionary populations moved by flows, and show how the speed of this
motion impacts the fixation probability of an invading mutant. Flows of
different speeds interpolate between evolutionary dynamics on fixed
heterogeneous graphs and in well-stirred populations. We find that the motion
has an active role in amplifying or suppressing selection, accomplished by
fragmenting and reconnecting the interaction graph. While increasing flow
speeds suppress selection for most evolutionary models, we identify
characteristic responses to flow for the different update rules we test. We
suggest these responses as a potential aid for choosing the most suitable
update rule for a given biological system.


  The art of viticulture and the quest for making wines has a long tradition
and it just started recently that mathematicians entered this field with their
main contribution of modelling alcoholic fermentation. These models consist of
systems of ordinary differential equations that describe the kinetics of the
bio-chemical reactions occurring in the fermentation process. The aim of this
paper is to present a new model of wine fermentation that accurately describes
the yeast dying component, the presence of glucose transporters, and the
formation of aromas and acids. Therefore the new model could become a valuable
tool to predict the taste of the wine and provide the starting point for an
emerging control technology that aims at improving the quality of the wine by
steering a well-behaved fermentation process that is also energetically more
efficient. Results of numerical simulations are presented that successfully
confirm the validity of the proposed model by comparison with real data.


  Seasonal influenza is a sometimes surprisingly impactful disease, causing
thousands of deaths per year along with much additional morbidity. Timely
knowledge of the outbreak state is valuable for managing an effective response.
The current state of the art is to gather this knowledge using in-person
patient contact. While accurate, this is time-consuming and expensive. This has
motivated inquiry into new approaches using internet activity traces, based on
the theory that lay observations of health status lead to informative features
in internet data.
  These approaches risk being deceived by activity traces having a
coincidental, rather than informative, relationship to disease incidence; to
our knowledge, this risk has not yet been quantitatively explored. We evaluated
both simulated and real activity traces of varying deceptiveness for influenza
incidence estimation using linear regression.
  We found that deceptiveness knowledge does reduce error in such estimates,
that it may help automatically-selected features perform as well or better than
features that require human curation, and that a semantic distance measure
derived from the Wikipedia article category tree serves as a useful proxy for
deceptiveness. This suggests that disease incidence estimation models should
incorporate not only data about how internet features map to incidence but also
additional data to estimate feature deceptiveness. By doing so, we may gain one
more step along the path to accurate, reliable disease incidence estimation
using internet data. This capability would improve public health by decreasing
the cost and increasing the timeliness of such estimates.


  The generalized Fibonacci recurrence $g_n=g_{n-k}+g_{n-m}$ was recently used
to demonstrate the theoretically optimal nature of limited senescence in
morphologically symmetrically dividing bacteria. Here, we study this recurrence
from a more abstract viewpoint, as a general model for asymmetric branching,
and interpret solutions for different initial conditions in terms of
branching-related quantities. We provide a compact diagrammatic representation
for the evolution of this process which leads to an explicit binomial identity
for the sums of elements lying on the diagonals $kx+my=n$ in Pascal's triangle
$\mathbb N_0\times \mathbb N_0\ni(x,y)\mapsto {x+y\choose x}$, previously
sought by Dickinson [Dic50], Raab [Raa63], and Green [Gre68].


  We review research papers which use game theory to model the decision making
of individuals during an epidemic, attempting to classify the literature and
identify the emerging trends in this field. We show that the literature can be
classified based on (i) type of population modelling (compartmental or
network-based), (ii) frequency of the game (non-iterative or iterative), and
(iii) type of strategy adoption (self-evaluation or imitation). We highlight
that the choice of model depends on many factors such as the type of immunity
the disease confers, the type of immunity the vaccine confers, and size of
population and level of mixing therein. We show that while early studies used
compartmental modelling with self-evaluation based strategy adoption, the
recent trend is to use network-based modelling with imitation-based strategy
adoption. Our review indicates that game theory continues to be an effective
tool to model intervention (vaccination or social distancing) decision-making
by individuals.


  Turkey hosts almost 3.5M refugees and has to face a humanitarian emergency of
unprecedented levels. We use mobile phone data to map the mobility patterns of
both Turkish and Syrian refugees, and use these patterns to build data-driven
computational models for quantifying the risk of epidemics spreading for
measles -- a disease having a satisfactory immunization coverage in Turkey but
not in Syria, due to the recent civil war -- while accounting for hypothetical
policies to integrate the refugees with the Turkish population. Our results
provide quantitative evidence that policies to enhance social integration
between refugees and the hosting population would reduce the transmission
potential of measles by almost 50%, preventing the onset of widespread large
epidemics in the country. Our results suggest that social segregation does not
hamper but rather boosts potential outbreaks of measles to a greater extent in
Syrian refugees but also in Turkish citizens, although to a lesser extent. This
is due to the fact that the high immunization coverage of Turkish citizens can
shield Syrian refugees from getting exposed to the infection and this in turn
reduces potential sources of infection and spillover of cases among Turkish
citizens as well, in a virtuous cycle reminiscent of herd immunity.


  More than 30 years ago, we published a paper entitled as abnormal chromatin
configuration and oncogenesis, which proposed the first hypothesis that links
oncogenesis to abnormal three-dimensional (3D) genome structure. Recently, many
studies have demonstrated that the 3D genome structure plays a major role in
oncogenesis, which strongly supports our hypothesis. In this paper, further
thoughts about our hypothesis is presented.


  Given a gene tree topology and a species tree topology, a coalescent history
represents a possible mapping of the list of gene tree coalescences to
associated branches of a species tree on which those coalescences take place.
Enumerative properties of coalescent histories have been of interest in the
analysis of relationships between gene trees and species trees. The simplest
enumerative result identifies a bijection between coalescent histories for a
matching caterpillar gene tree and species tree with monotonic paths that do
not cross the diagonal of a square lattice, establishing that the associated
number of coalescent histories for $n$-taxon matching caterpillar trees ($n
\geqslant 2$) is the Catalan number $C_{n-1} = \frac{1}{n} {2n-2 \choose n-1}$.
Here, we show that a similar bijection applies for \emph{non-matching}
caterpillars, connecting coalescent histories for a non-matching caterpillar
gene tree and species tree to a class of \emph{roadblocked} monotonic paths.
The result provides a simplified algorithm for enumerating coalescent histories
in the non-matching caterpillar case. It enables a rapid proof of a known
result that given a caterpillar species tree, no non-matching caterpillar gene
tree has a number of coalescent histories exceeding that of the matching gene
tree. Additional results on coalescent histories can be obtained by a bijection
between permissible roadblocked monotonic paths and Dyck paths. We study the
number of coalescent histories for non-matching caterpillar gene trees that
differ from the species tree by nearest-neighbor-interchange and
subtree-prune-and-regraft moves, characterizing the non-matching caterpillar
with the largest number of coalescent histories. We discuss the implications of
the results for the study of the combinatorics of gene trees and species trees.


  Can prelife proceed without cell division? A recently proposed mechanism
suggests that transient compartmentalization could have preceded cell division
in prebiotic scenarios. Here, we study transient compartmentalization dynamics
in the presence of mutations and noise in replication, as both can be
detrimental the survival of compartments. Our study comprises situations where
compartments contain uncoupled autocatalytic reactions feeding on a common
resource, and systems based on RNA molecules copied by replicases, following a
recent experimental study.
  Using the theory of branching processes, we show analytically that two
regimes are possible. In the diffusion-limited regime, replication is
asynchronous which leads to a large variability in the composition of
compartments. In contrast, in a replication-limited regime, the growth is
synchronous and thus the compositional variability is low. Typically, simple
autocatalysts are in the former regime, while polymeric replicators can access
the latter.
  For deterministic growth dynamics, we introduce mutations that turn
functional replicators into parasites. We derive the phase boundary separating
coexistence or parasite dominance as a function of relative growth, inoculation
size and mutation rate. We show that transient compartmentalization allows
coexistence beyond the classical error threshold, above which the parasite
dominates. Our findings invite to revisit major prebiotic transitions, notably
the transitions towards cooperation, complex polymers and cell division.


  In this paper, a mathematical model of pneumococcal pneumonia with time
delays is proposed. The stability theory of delay differential equations is
used to analyze the model. The results show that the disease-free equilibrium
is asymptotically stable if the control reproduction ratio R0 is less than
unity and unstable otherwise. The stability of equilibria with delays shows
that the endemic equilibrium is locally stable without delays and stable if the
delays are under conditions. The existence of Hopf-bifurcation is investigated
and transversality conditions proved. The model results suggest that as the
respective delays exceed some critical value past the endemic equilibrium, the
system loses stability through the process of local birth or death of
oscillations. Further, a decrease or an increase in the delays leads to
asymptotic stability or instability of the endemic equilibrium respectively.
The analytical results, are supported by numerical simulations.
  Keywords: Time delay, Pneumococcal pneumonia, Vaccination, Stability,
Hopf-bifurcation


  The study of human genes and diseases is very rewarding and can lead to
improvements in healthcare, disease diagnostics and drug discovery. In this
paper, we further our previous study on gene disease relationship specifically
with the multifunctional genes. We investigate the multifunctional gene disease
relationship based on the published molecular function annotations of genes
from the Gene Ontology which is the most comprehensive source on gene
functions.


  The stochastic extinction and stability in the mean of a family of SEIRS
malaria models with a general nonlinear incidence rate is presented. The
dynamics is driven by independent white noise processes from the disease
transmission and natural death rates. The basic reproduction number
$R^{*}_{0}$, the expected survival probability of the plasmodium
$E(e^{-(\mu_{v}T_{1}+\mu T_{2})})$, and other threshold values are calculated.
A sample Lyapunov exponential analysis for the system is utilized to obtain
extinction results. Moreover, the rate of extinction of malaria is estimated,
and innovative local Martingale and Lyapunov functional techniques are applied
to establish the strong persistence, and asymptotic stability in the mean of
the malaria-free steady population. %The extinction of malaria depends on
$R^{*}_{0}$, and $E(e^{-(\mu_{v}T_{1}+\mu T_{2})})$. Moreover, for either
$R^{*}_{0}<1$, or $E(e^{-(\mu_{v}T_{1}+\mu T_{2})})<\frac{1}{R^{*}_{0}}$,
whenever $R^{*}_{0}\geq 1$, respectively, extinction of malaria occurs.
Furthermore, the robustness of these threshold conditions to the intensity of
noise from the disease transmission rate is exhibited. Numerical simulation
results are presented.


  This paper investigates the deterministic extinction and permanence of a
family of SEIRS malaria models with multiple random delays, and with a general
nonlinear incidence rate. The conditions for the extinction and permanence of
the disease are presented. For the persistence of disease, improved analytical
techniques are employed to obtain eventual lower bounds for the states of the
system, while extinction results are obtained via analyzing the Lyapunov
exponent. The basic reproduction number $R^{*}_{0}$ is calculated, and for
$R^{*}_{0}<1$, extinction of disease occurs, while for $R^{*}_{0}\geq 1$, the
disease persists. But if $R^{*}_{0}\geq 1$, and the expected survival
probability of the plasmodium $E(e^{-(\mu_{v}T_{1}+\mu
T_{2})})<\frac{1}{R^{*}_{0}}$, then extinction occurs. Results are interpreted,
and numerical simulation examples are presented to (a) compare and evaluate the
impacts of the intensity of the incidence rate, and the intensity of response
to malaria control, on the persistence of malaria in the human population, and
(2) to justify the asymptotic stability of the disease-free equilibrium.


  Realistic evolutionary fitness landscapes are notoriously difficult to
construct. A recent cutting-edge model of virus assembly consists of a
dodecahedral capsid with $12$ corresponding packaging signals in three affinity
bands. This whole genome/phenotype space consisting of $3^{12}$ genomes has
been explored via computationally expensive stochastic assembly models, giving
a fitness landscape in terms of the assembly efficiency. Using latest
machine-learning techniques by establishing a neural network, we show that the
intensive computation can be short-circuited in a matter of minutes to
astounding accuracy.


  This study presents a family of stochastic models for the dynamics of
influenza in a closed human population. We consider treatment for the disease
in the form of vaccination, and incorporate the periods of effectiveness of the
vaccine and infectiousness for the individuals in the population. Our model is
a SVIR model, with trinomial transition probabilities, where all individuals
who recover from the disease acquire permanent natural immunity against the
strain of the disease. Special SVIR models in the family are presented, based
on the structure of the probability of getting infection and vaccination at any
instant. The methods of maximum likelihood, and expectation maximization are
derived for the parameters of the chain. Moreover, estimators for some
epidemiological assessment parameters, such as the basic reproduction number
are computed. Numerical simulation examples are presented for the model.


  Background: Although Alzheimer's disease (AD) is a central nervous system
disease and type 2 diabetes mellitus (T2DM) is a metabolic disorder, an
increasing number of genetic epidemiological studies show clear link between AD
and T2DM. The current approach to uncovering the shared pathways between AD and
T2DM involves association analysis; however, such analyses lack power to
discover the mechanisms of the diseases. Methods: We develop novel statistical
methods to shift the current paradigm of genetic analysis from association
analysis to deep causal inference for uncovering the shared mechanisms between
AD and T2DM, and develop pipelines to infer multilevel omics causal networks
which lead to shifting the current paradigm of genetic analysis from genetic
analysis alone to integrated causal genomic, epigenomic, transcriptional and
phenotypic data analysis. To discover common causal paths from genetic variants
to AD and T2DM, we also develop algorithms that can automatically search the
causal paths from genetic variants to diseases and Results: The proposed
methods and algorithms are applied to ROSMAP dataset with 432 individuals who
simultaneously had genotype, RNA-seq, DNA methylation and some phenotypes. We
construct multi-omics causal networks and identify 13 shared causal genes, 16
shared causal pathways between AD and T2DM, and 754 gene expression and 101
gene methylation nodes that were connected to both AD and T2DM in multi-omics
causal networks. Conclusions: The results of application of the proposed
pipelines for identifying causal paths to real data analysis of AD and T2DM
provided strong evidence to support the link between AD and T2DM and unraveled
causal mechanism to explain this link.


  A feed-forward loop (FFL) is a common gene-regulatory motif in which usually
the upstream regulator is a protein, a transcription factor, that regulates the
expression of the target protein in two parallel pathways. Here, we study a
distinct sRNA-driven FFL (sFFL) discovered recently in Salmonella enterica.
Unlike previously studied transcriptional FFLs (tFFL) and sRNA-mediated FFLs
(smFFL), here the upstream regulator is an sRNA that activates the target
protein and its transcriptional activator. Such sFFL has not been subjected to
rigorous analysis. We, therefore, set out to understand two aspects. First is a
quantitative comparison of the regulatory response of sFFL with tFFL and smFFL
using a differential equation framework. Since the process of gene expression
is inherently stochastic, the second objective is to find how the noise affects
the functionality of sFFL. We find the response of sFFLto be stronger, faster,
and more sensitive to the initial concentration of the upstream regulator than
tFFL and smFFL. Further, a generating function based analysis and stochastic
simulations lead to a non-trivial prediction that an optimal noise filtration
can be attained depending on the synthesis rate of the sRNA and the degradation
rate of the transcriptional activator. We conclude that in sFFL, sRNA plays a
critical role not only in driving a rapid and strong response, but also a
reliable response that depends critically on its concentration. Given the
advantages of sFFL brought out in this work, it should not be surprising if
future work reveals their employment in different biological contexts.