  Extensive cooperation among biologically unrelated individuals is uniquely
human and much current research attempts to explain this fact. We draw upon
social, cultural, and psychological aspects of human uniqueness to present an
integrated theory of human cooperation that explains aspects of human
cooperation that are problematic for other theories (e.g., defector invasion
avoidance, preferential assortment to exclude free riders, and the second order
free rider problem). We propose that the evolution of human cooperative
behavior required (1) a capacity for self-sustained, self-referential thought
manifested as an integrated worldview, including a sense of identity and point
of view, and (2) the cultural formation of kinship-based social organizational
systems within which social identities can be established and transmitted
through enculturation. Human cooperative behavior arose, we argue, through the
acquisition of a culturally grounded social identity that included the
expectation of cooperation among kin. This identity is linked to basic survival
instincts by emotions that are mentally experienced as culture laden feelings.
As a consequence, individuals are motivated to cooperate with those perceived
culturally as kin, while deviations from expected social behavior are
experienced as threatening to one's social identity, leading to punishment of
those seen as violating cultural expectations regarding socially proper
behavior.


  As the number of solved protein structures increases, the opportunities for
meta-analysis of this dataset increase too. Protein structures are known to be
formed of domains; structural and functional subunits that are often repeated
across sets of proteins. These domains generally form compact, globular
regions, and are therefore often easily identifiable by inspection, yet the
problem of automatically fragmenting the protein into these compact
substructures remains computationally challenging. Existing domain
classification methods focus on finding subregions of protein structure that
are conserved, rather than finding a decomposition which spans the full protein
structure. However, such a decomposition would find ready application in
coarse-graining molecular dynamics, analysing the protein's topology, in de
novo protein design and in fitting electron microscopy maps. Here, we present a
tool for performing this modular decomposition using the Infomap community
detection algorithm. The protein structure is abstracted into a network in
which its amino acids are the nodes, and where the edges are generated using a
simple proximity test. Infomap can then be used to identify highly
intra-connected regions of the protein. We perform this decomposition
systematically across 4000 distinct protein structures, taken from the Protein
Data Bank. The decomposition obtained correlates well with existing PFAM
sequence classifications, but has the advantage of spanning the full protein,
with the potential for novel domains. The coarse-grained network formed by the
communities can also be used as a proxy for protein topology at the
single-chain level; we demonstrate that grouping these proteins by their
coarse-grained network results in a functionally significant classification.


  The phylogenetic effective sample size is a parameter that has as its goal
the quantification of the amount of independent signal in a phylogenetically
correlated sample. It was studied for Brownian motion and Ornstein-Uhlenbeck
models of trait evolution. Here, we study this composite parameter when the
trait is allowed to jump at speciation points of the phylogeny. Our numerical
study indicates that there is a non-trivial limit as the effect of jumps grows.
The limit depends on the value of the drift parameter of the Ornstein-Uhlenbeck
process.


  Motivation: Working with a large number of genomes simultaneously is of great
interest in genetic population and comparative genomics research. Bubbles
discovery in multi-genomes coloured de bruijn graph for de novo genome assembly
is a problem that can be translated to cycles enumeration in graph theory.
Cycle enumerations algorithms in big and complex de Bruijn graphs are time
consuming. Specialised fast algorithms for efficient bubble search are needed
for coloured de bruijn graph variant calling applications. In coloured de
Bruijn graphs, bubble paths coverages are used in downstream variants calling
analysis. Results: In this paper, we introduce a fast parallel graph search for
different K-mer cycle sizes. Coloured path coverages are used for SNP
prediction. The graph search method uses a combined multi-node and multi-core
design to speeds up cycles enumeration. The search algorithm uses an index
extracted from the raw assembly of a coloured de Bruijn graph stored in a hash
table. The index is distributed across different CPU-cores, in a shared memory
HPC compute node, to build undirected subgraphs then search independently and
simultaneously specific cycle sizes. This same index can also be split between
several HPC compute nodes to take advantage of as many CPU-cores available to
the user. The local neighbourhood parallel search approach reduces the graph's
complexity and facilitate cycles search of a multi-colour de Bruijn graph. The
search algorithm is incorporated into $\Pi$-cyc application and tested on a
number of Schizosaccharomyces Pombe genomes. Availability: $\Pi$-cyc is an
open-source software available at www.github.com/2kplus2P


  We analyze free energy estimators from simulation trials mimicking
single-molecule pulling experiments on a neck linker of a kinesin motor. For
that purpose, we have performed a version of steered molecular dynamics (SMD)
calculations. The sample trajectories have been analyzed to derive distribution
of work done on the system. In order to induce unfolding of the linker, we have
stretched the molecule at a constant pulling force and allowed for a subsequent
relaxation of its structure. The use of fluctuation relations (FR) relevant to
non-equilibrium systems subject to thermal fluctuations allows us to assess the
difference in free energy between stretched and relaxed conformations. To
further understand effects of potential mutations on elastic properties of the
linker, we have performed similar in silico studies on a structure formed of a
polyalanine sequence (Ala-only) and on three other structures, created by
substituting selected types of amino acid residues in the linker's sequence
with alanine (Ala) ones. The results of SMD simulations indicate a crucial role
played by the Asparagine (Asn) and Lysine (Lys) residues in controlling
stretching and relaxation properties of the linker domain of the motor.


  Investigating the emergence of a particular cell type is a recurring theme in
models of growing cellular populations. The evolution of resistance to therapy
is a classic example. Common questions are: when does the cell type first
occur, and via which sequence of steps is it most likely to emerge? For growing
populations, these questions can be formulated in a general framework of
branching processes spreading through a graph from a root to a target vertex.
Cells have a particular fitness value on each vertex and can transition along
edges at specific rates. Vertices represents cell states, say \mic{genotypes
}or physical locations, while possible transitions are acquiring a mutation or
cell migration. We focus on the setting where cells at the root vertex have the
highest fitness and transition rates are small. Simple formulas are derived for
the time to reach the target vertex and for the probability that it is reached
along a given path in the graph. We demonstrate our results on \mic{several
scenarios relevant to the emergence of drug resistance}, including: the
orderings of resistance-conferring mutations in bacteria and the impact of
imperfect drug penetration in cancer.


  Gene regulatory networks play a crucial role in controlling an organism's
biological processes, which is why there is significant interest in developing
computational methods that are able to extract their structure from
high-throughput genetic data. A typical approach consists of a series of
conditional independence tests on the covariance structure meant to
progressively reduce the space of possible causal models. We propose a novel
efficient Bayesian method for discovering the local causal relationships among
triplets of (normally distributed) variables. In our approach, we score the
patterns in the covariance matrix in one go and we incorporate the available
background knowledge in the form of priors over causal structures. Our method
is flexible in the sense that it allows for different types of causal
structures and assumptions. We apply the approach to the task of inferring gene
regulatory networks by learning regulatory relationships between gene
expression levels. We show that our algorithm produces stable and conservative
posterior probability estimates over local causal structures that can be used
to derive an honest ranking of the most meaningful regulatory relationships. We
demonstrate the stability and efficacy of our method both on simulated data and
on real-world data from an experiment on yeast.


  We study the probability of fixation in a stochastic two-species competition
model. By identifying a naturally occurring fast timescale, we derive an
approximation to the associated backward Kolmogorov equation that allows us to
obtain an explicit closed form solution for the probability of fixation of
either species. We use our result to study fitness tradeoff strategies and show
that, despite some tradeoffs having nearly negligible effects on the
corresponding deterministic dynamics, they can have large implications for the
outcome of the stochastic system.


  Our understanding of how chromosomes structurally organize and dynamically
interact has been revolutionized through the lens of long-chain polymer
physics. Major protein contributors to chromosome structure and dynamics are
condensin and cohesin that stochastically generate loops within and between
chains, and entrap proximal strands of sister chromatids. In this paper, we
explore the ability of transient, protein-mediated, gene-gene crosslinks to
induce clusters of genes, thereby dynamic architecture, within the highly
repeated ribosomal DNA that comprises the nucleolus of budding yeast. We
implement three approaches: live cell microscopy; computational modeling of the
full genome during G1 in budding yeast, exploring four decades of timescales
for transient crosslinks between 5k bp domains in the nucleolus on Chromosome
XII; and, temporal network models with automated community detection algorithms
applied to the full range of 4D modeling datasets. The data analysis tools
detect and track gene clusters, their size, number, persistence time, and their
plasticity. Of biological significance, our analysis reveals an optimal mean
crosslink lifetime that promotes pairwise and cluster gene interactions through
"flexible" clustering. In this state, large gene clusters self-assemble yet
frequently interact, marked by gene exchanges between clusters, which in turn
maximizes global gene interactions in the nucleolus. This regime stands between
two limiting cases each with far less global gene interactions: with shorter
crosslink lifetimes, "rigid" clustering emerges with clusters that interact
infrequently; with longer crosslink lifetimes, there is a dissolution of
clusters. These observations are compared with imaging experiments on a normal
yeast strain and two condensin-modified mutant cell strains, applying the same
image analysis pipeline to the experimental and simulated datasets.


  Metagenome assembly is the process of transforming a set of short,
overlapping, and potentially erroneous DNA segments from environmental samples
into the accurate representation of the underlying microbiomes's genomes.
State-of-the-art tools require big shared memory machines and cannot handle
contemporary metagenome datasets that exceed Terabytes in size. In this paper,
we introduce the MetaHipMer pipeline, a high-quality and high-performance
metagenome assembler that employs an iterative de Bruijn graph approach.
MetaHipMer leverages a specialized scaffolding algorithm that produces long
scaffolds and accommodates the idiosyncrasies of metagenomes. MetaHipMer is
end-to-end parallelized using the Unified Parallel C language and therefore can
run seamlessly on shared and distributed-memory systems. Experimental results
show that MetaHipMer matches or outperforms the state-of-the-art tools in terms
of accuracy. Moreover, MetaHipMer scales efficiently to large concurrencies and
is able to assemble previously intractable grand challenge metagenomes. We
demonstrate the unprecedented capability of MetaHipMer by computing the first
full assembly of the Twitchell Wetlands dataset, consisting of 7.5 billion
reads - size 2.6 TBytes.


  The simplest model of DNA mechanics describes the double helix as a
continuous rod with twist and bend elasticity. Recent work has discussed the
relevance of a little-studied coupling $G$ between twisting and bending, known
to arise from the groove asymmetry of the DNA double helix. Here, the effect of
$G$ on the statistical mechanics of long DNA molecules subject to applied
forces and torques is investigated. We present a perturbative calculation of
the effective torsional stiffness $C_\text{eff}$ for small twist-bend coupling.
We find that the "bare" $G$ is "screened" by thermal fluctuations, in the sense
that the low-force, long-molecule effective free energy is that of a model with
$G=0$, but with long-wavelength bending and twisting rigidities that are
shifted by $G$-dependent amounts. Using results for torsional and bending
rigidities for freely-fluctuating DNA, we show how our perturbative results can
be extended to a non-perturbative regime. These results are in excellent
agreement with numerical calculations for Monte Carlo "triad" and molecular
dynamics "oxDNA" models, characterized by different degrees of coarse-graining,
validating the perturbative and non-perturbative analyses. While our theory is
in generally-good quantitative agreement with experiment, the predicted
torsional stiffness does systematically deviate from experimental data,
suggesting that there are as-yet-uncharacterized aspects of DNA
twisting-stretching mechanics relevant to low-force, long-molecule mechanical
response, which are not captured by widely-used coarse-grained models.


  The mechanism responsible for the compaction of the genomic DNA of bacteria
inside a structure called the nucleoid is a longstanding but still lively
debated question. Most puzzling is the fact that the nucleoid occupies only a
small fraction of the cell, although it is not separated from the rest of the
cytoplasm by any membrane and would occupy a volume about thousand times larger
outside from the cell. Here, by performing numerical simulations with
coarse-grained models, we elaborate on the conjecture that the formation of the
nucleoid may result from a segregative phase separation mechanism driven by the
demixing of the DNA coil and non-binding globular macromolecules present in the
cytoplasm, presumably functional ribosomes. Simulations performed with crowders
having spherical, dumbbell or octahedral geometry highlight the sensitive
dependence of the level of DNA compaction on the dissymmetry of DNA/DNA,
DNA/crowder, and crowder/crowder repulsive interactions, thereby supporting the
segregative phase separation scenario. Simulations also consistently predict
much stronger DNA compaction close to the jamming threshold. Moreover,
simulations performed with crowders of different sizes suggest that the final
density distribution of each species results from the competition between
thermodynamic forces and steric hindrance, so that bigger crowders are expelled
selectively from the nucleoid only at moderate total crowder concentrations.
This work leads to several predictions, which may eventually be tested
experimentally.


  Parametrized polynomial ordinary differential equation systems are broadly
used for modeling, specially in the study of biochemical reaction networks
under the assumption of mass-action kinetics. Understanding the qualitative
behavior of the solutions with respect to the parameter values gives rise to
complex problems within real algebraic geometry, concerning the study of the
signs of multivariate polynomials over the positive orthant. In this work we
provide further insight into the number of positive steady states of a
benchmark model, namely the two-site phosphorylation cycle. In particular, we
provide new conditions on the reaction rate constants for the existence of one
or three positive steady states, partially filling a gap left in previous
works.


  In this paper we show how to use one or more assembled or partially assembled
genome as the basis for a compressed full-text index of its readset.
Specifically, we build a labelled tree by taking the assembled genome as a
trunk and grafting onto it the reads that align to it, at the starting
positions of their alignments. Next, we compute the eXtended Burrows-Wheeler
Transform (XBWT) of the resulting labelled tree and build a compressed
full-text index on that. Although this index can occasionally return false
positives, it is usually much more compact than the alternatives. Following the
established practice for datasets with many repetitions, we compare different
full-text indices by looking at the number of runs in the transformed strings.
For a human Chr19 readset our preliminary experiments show that eliminating
separators characters from the EBWT reduces the number of runs by 19\%, from
220 million to 178 million, and using the XBWT reduces it by a further 15\%, to
150 million.


  Due to the low density of envelope (Env) spikes on the surface of HIV-1,
neutralizing IgG antibodies rarely bind bivalently using both antigen-binding
arms (Fabs) to crosslink between spikes (inter-spike crosslinking), instead
resorting to weaker monovalent binding that is more sensitive to Env mutations.
Synthetic antibodies designed to bivalently bind a single Env trimer
(intra-spike crosslinking) were previously shown to exhibit increased
neutralization potencies. In initial work, diFabs joined by varying lengths of
rigid double-stranded DNA (dsDNA) were considered. Anticipating future
experiments to improve synthetic antibodies, we investigate whether linkers
with different rigidities could enhance diFab potency by modeling DNA-Fabs
containing different combinations of rigid dsDNA and flexible single-stranded
DNA (ssDNA) and characterizing their neutralization potential. Model
predictions suggest that while a long flexible polymer may be capable of
bivalent binding, it exhibits weak neutralization due to the large loss in
entropic degrees of freedom when both Fabs are bound. In contrast, the
strongest neutralization potencies are predicted to require a rigid linker that
optimally spans the distance between two Fab binding sites on an Env trimer,
and avidity can be further boosted by incorporating more Fabs into these
constructs. These results inform the design of multivalent anti-HIV-1
therapeutics that utilize avidity effects to remain potent against HIV-1 in the
face of the rapid mutation of Env spikes.


  We study a general setting of neutral evolution in which the population is of
finite, constant size and can have spatial structure. Mutation leads to
different genetic types ("traits"), which can be discrete or continuous. Under
minimal assumptions, we show that the marginal trait distributions of the
evolutionary process, which specify the probability that any given individual
has a certain trait, all converge to the stationary distribution of the
mutation process. In particular, the stationary frequencies of traits in the
population are independent of its size, spatial structure, and evolutionary
update rule, and these frequencies can be calculated by evaluating a simple
stochastic process describing a population of size one (i.e. the mutation
process itself). We conclude by analyzing mixing times, which characterize
rates of convergence of the mutation process along the lineages, in terms of
demographic variables of the evolutionary process.


  Motivation: The ability to generate massive amounts of sequencing data
continues to overwhelm the processing capability of existing algorithms and
compute infrastructures. In this work, we explore the use of hardware/software
co-design and hardware acceleration to significantly reduce the execution time
of short sequence alignment, a crucial step in analyzing sequenced genomes. We
introduce Shouji, a highly-parallel and accurate pre-alignment filter that
remarkably reduces the need for computationally-costly dynamic programming
algorithms. The first key idea of our proposed pre-alignment filter is to
provide high filtering accuracy by correctly detecting all common subsequences
shared between two given sequences. The second key idea is to design a hardware
accelerator that adopts modern FPGA (Field-Programmable Gate Array)
architectures to further boost the performance of our algorithm.
  Results: Shouji significantly improves the accuracy of pre-alignment
filtering by up to two orders of magnitude compared to the state-of-the-art
pre-alignment filters, GateKeeper and SHD. Our FPGA-based accelerator is up to
three orders of magnitude faster than the equivalent CPU implementation of
Shouji. Using a single FPGA chip, we benchmark the benefits of integrating
Shouji with five state-of-the-art sequence aligners, designed for different
computing platforms. The addition of Shouji as a pre-alignment step reduces the
execution time of the five state-of-the-art sequence aligners by up to 18.8x.
Shouji can be adapted for any bioinformatics pipeline that performs sequence
alignment for verification. Unlike most existing methods that aim to accelerate
sequence alignment, Shouji does not sacrifice any of the aligner capabilities,
as it does not modify or replace the alignment step.
  Availability: https://github.com/CMU-SAFARI/Shouji


  Motivation: Despite being often perceived as the main contributors to cell
fate and physiology, genes alone cannot predict cellular phenotype. During the
process of gene expression, 95% of human genes can code for multiple proteins
due to alternative splicing. While most splice variants of a gene carry the
same function, variants within some key genes can have remarkably different
roles. To bridge the gap between genotype and phenotype, condition- and
tissue-specific models of metabolism have been constructed. However, current
metabolic models only include information at the gene level. Consequently, as
recently acknowledged by the scientific community, common situations where
changes in splice-isoform expression levels alter the metabolic outcome cannot
be modeled. Results: We here propose GEMsplice, the first method for the
incorporation of splice-isoform expression data into genome-scale metabolic
models. Using GEMsplice, we make full use of RNA-Seq quantitative expression
profiles to predict, for the first time, the effects of splice isoform-level
changes in the metabolism of 1455 patients with 31 different breast cancer
types. We validate GEMsplice by generating cancer-versus-normal predictions on
metabolic pathways, and by comparing with gene-level approaches and available
literature on pathways affected by breast cancer. GEMsplice is freely available
for academic use at https://github.com/GEMsplice/GEMsplice_code. Compared to
state-of-the-art methods, we anticipate that GEMsplice will enable for the
first time computational analyses at transcript level with splice-isoform
resolution.


  Comprehensive and confident identifications of metabolites and other
chemicals in complex samples will revolutionize our understanding of the role
these chemically diverse molecules play in biological systems. Despite recent
advances, metabolomics studies still result in the detection of a
disproportionate number of features than cannot be confidently assigned to a
chemical structure. This inadequacy is driven by the single most significant
limitation in metabolomics: the reliance on reference libraries constructed by
analysis of authentic reference chemicals. To this end, we have developed the
in silico chemical library engine (ISiCLE), a high-performance
computing-friendly cheminformatics workflow for generating libraries of
chemical properties. In the instantiation described here, we predict probable
three-dimensional molecular conformers using chemical identifiers as input,
from which collision cross sections (CCS) are derived. The approach employs
state-of-the-art first-principles simulation, distinguished by use of molecular
dynamics, quantum chemistry, and ion mobility calculations to generate
structures and libraries, all without training data. Importantly, optimization
of ISiCLE included a refactoring of the popular MOBCAL code for
trajectory-based mobility calculations, improving its computational efficiency
by over two orders of magnitude. Calculated CCS values were validated against
1,983 experimentally-measured CCS values and compared to previously reported
CCS calculation approaches. An online database is introduced for sharing both
calculated and experimental CCS values (metabolomics.pnnl.gov), initially
including a CCS library with over 1 million entries. Finally, three successful
applications of molecule characterization using calculated CCS are described.
This work represents a promising method to address the limitations of small
molecule identification.


  We use the oxDNA coarse-grained model to provide a detailed characterization
of the fundamental structural properties of DNA origamis, focussing on
archetypal 2D and 3D origamis. The model reproduces well the characteristic
pattern of helix bending in a 2D origami, showing that it stems from the
intrinsic tendency of anti-parallel four-way junctions to splay apart, a
tendency that is enhanced both by less screened electrostatic interactions and
by increased thermal motion. We also compare to the structure of a 3D origami
whose structure has been determined by cryo-electron microscopy. The oxDNA
average structure has a root-mean-square deviation from the experimental
structure of 8.4 Angstrom, which is of the order of the experimental
resolution. These results illustrate that the oxDNA model is capable of
providing detailed and accurate insights into the structure of DNA origamis,
and has the potential to be used to routinely pre-screen putative origami
designs.


  Moran or Wright-Fisher processes are probably the most well known model to
study the evolution of a population under various effects. Our object of study
will be the Simpson index which measures the level of diversity of the
population, one of the key parameter for ecologists who study for example
forest dynamics. Following ecological motivations, we will consider here the
case where there are various species with fitness and immigration parameters
being random processes (and thus time evolving). To measure biodiversity,
ecologists generally use the Simpson index, who has no closed formula, except
in the neutral (no selection) case via a backward approach, and which is
difficult to evaluate even numerically when the population size is large. Our
approach relies on the large population limit in the "weak" selection case, and
thus to give a procedure which enable us to approximate, with controlled rate,
the expectation of the Simpson index at fixed time. Our approach will be
forward and valid for all time, which is the main difference with the
historical approach of Kingman, or Krone-Neuhauser. We will also study the long
time behaviour of the Wright-Fisher process in a simplified setting, allowing
us to get a full picture for the approximation of the expectation of the
Simpson index.


  We consider the problem of determining the time evolution of a trait
distribution in a mathematical model of non-uniform populations with parametric
heterogeneity. This means that we consider only heterogeneous populations in
which heterogeneity is described by an individual specific parameter that
differs in general from individual to individual, but does not change with time
for the whole lifespan of this individual. Such a restriction allows obtaining
a number of simple and yet important analytical results. In particular we show
that initial assumptions on time-dependent behavior of various characteristics,
such as the mean, variance, of coefficient of variation, restrict severely
possible choices for the exact form of the trait distribution. We illustrate
our findings by in-depth analysis of the variance evolution. We also reanalyze
a well known mathematical model for gypsy moth population showing that the
knowledge of how distributions evolve allows producing oscillatory behaviors
for highly heterogeneous populations.


  We introduce an R package, PCMBase, to rapidly calculate the likelihood for
multivariate phylogenetic comparative methods. The package is not specific to
particular models but offers the user the functionality to very easily
implement a wide range of models where the transition along a branch is
multivariate normal. We demonstrate the package's possibilities on the now
standard, multitrait Ornstein-Uhlenbeck process as well as the novel
multivariate punctuated equilibrium model. The package can handle trees of
various types (e.g. ultrametric, nonultrametric, polytomies, e.t.c.), as well
as measurement error, missing measurements or non-existing traits for some of
the species in the tree.


  Many theoretical models have been formulated to better understand the
coevolutionary patterns that emerge from antagonistic interactions. These
models usually assume that the attacks by the exploiters are random, so the
effect of victim selection by exploiters on coevolutionary patterns remains
unexplored. Here we analytically studied the payoff for predators and prey
under coevolution assuming that every individual predator can attack only a
small number of prey any given time, considering two scenarios: (i) predation
occurs at random; (ii) predators select prey according to phenotype matching.
We also develop an individual based model to verify the robustness of our
analytical prediction. We show that both scenarios result in well known similar
coevolutionary patterns if population sizes are sufficiently high: symmetrical
coevolutionary branching and symmetrical coevolutionary cycling (Red Queen
dynamics). However, for small population sizes, prey selection can cause
unexpected coevolutionary patterns. One is the breaking of symmetry of the
coevolutionary pattern, where the phenotypes evolve towards one of two
evolutionarily stable patterns. As population size increases, the phenotypes
oscillate between these two values in a novel form of Red Queen dynamics, the
episodic reversal between the two stable patterns. Thus, prey selection causes
prey phenotypes to evolve towards more extreme values, which reduces the
fitness of both predators and prey, increasing the likelihood of extinction.


  We present analysis of a novel tool for protein secondary structure
prediction using the recently-investigated Neural Machine Translation
framework. The tool provides a fast and accurate folding prediction based on
primary structure with subsecond prediction time even for batched inputs. We
hypothesize that Neural Machine Translation can improve upon current predictive
accuracy by better encoding complex relationships between nearby but
non-adjacent amino acids. We overview our modifications to the framework in
order to improve accuracy on protein sequences. We report 65.9% Q3 accuracy and
analyze the strengths and weaknesses of our predictive model.


  DNA is a flexible molecule, but the degree of its flexibility is subject to
debate. The commonly-accepted persistence length of $l_p \approx 500\,$\AA\ is
inconsistent with recent studies on short-chain DNA that show much greater
flexibility but do not probe its origin. We have performed X-ray and neutron
small-angle scattering on a short DNA sequence containing a strong nucleosome
positioning element, and analyzed the results using a modified Kratky-Porod
model to determine possible conformations. Our results support a hypothesis
from Crick and Klug in 1975 that some DNA sequences in solution can have sharp
kinks, potentially resolving the discrepancy. Our conclusions are supported by
measurements on a radiation-damaged sample, where single-strand breaks lead to
increased flexibility and by an analysis of data from another sequence, which
does not have kinks, but where our method can detect a locally enhanced
flexibility due to an $AT$-domain.


  Metabolic modelling has entered a mature phase with dozens of methods and
software implementations available to the practitioner and the theoretician. It
is not easy for a modeller to be able to see the wood (or the forest) for the
trees. Driven by this analogy, we here present a "forest" of principal methods
used for constraint-based modelling in systems biology. This provides a
tree-based view of methods available to prospective modellers, also available
in interactive version at http://modellingmetabolism.net, where it will be kept
updated with new methods after the publication of the present manuscript. Our
updated classification of existing methods and tools highlights the most
promising in the different branches, with the aim to develop a vision of how
existing methods could hybridise and become more complex. We then provide the
first hands-on tutorial for multi-objective optimisation of metabolic models in
R. We finally discuss the implementation of multi-view machine-learning
approaches in poly-omic integration. Throughout this work, we demonstrate the
optimisation of trade-offs between multiple metabolic objectives, with a focus
on omic data integration through machine learning. We anticipate that the
combination of a survey, a perspective on multi-view machine learning, and a
step-by-step R tutorial should be of interest for both the beginner and the
advanced user.


  Synthetic biology aims at designing modular genetic circuits that can be
assembled according to the desired function. When embedded in a cell, a circuit
module becomes a small subnetwork within a larger environmental network, and
its dynamics is therefore affected by potentially unknown interactions with the
environment. It is well-known that the presence of the environment not only
causes extrinsic noise but also memory effects, which means that the dynamics
of the subnetwork is affected by its past states via a memory function that is
characteristic of the environment. We study several generic scenarios for the
coupling between a small module and a larger environment, with the environment
consisting of a chain of mono-molecular reactions. By mapping the dynamics of
this coupled system onto random walks, we are able to give exact analytical
expressions for the arising memory functions. Hence, our results give insights
into the possible types of memory functions and thereby help to better predict
subnetwork dynamics.


  In extreme cold weather, living organisms produce Antifreeze Proteins (AFPs)
to counter the otherwise lethal intracellular formation of ice. Structures and
sequences of various AFPs exhibit a high degree of heterogeneity, consequently
the prediction of the AFPs is considered to be a challenging task. In this
research, we propose to handle this arduous manifold learning task using the
notion of localized processing. In particular an AFP sequence is segmented into
two sub-segments each of which is analyzed for amino acid and di-peptide
compositions. We propose to use only the most significant features using the
concept of information gain (IG) followed by a random forest classification
approach. The proposed RAFP-Pred achieved an excellent performance on a number
of standard datasets. We report a high Youden's index
(sensitivity+specificity-1) value of 0.75 on the standard independent test data
set outperforming the AFP-PseAAC, AFP\_PSSM, AFP-Pred and iAFP by a margin of
0.05, 0.06, 0.14 and 0.68 respectively. The verification rate on the UniProKB
dataset is found to be 83.19\% which is substantially superior to the 57.18\%
reported for the iAFP method.


  We study the effect of correlations in generation times on the dynamics of
population growth of microorganisms. We show that any non-zero correlation that
is due to cell-size regulation, no matter how small, induces long-term
oscillations in the population growth rate. The population only reaches its
steady state when we include the often-neglected variability in the growth
rates of individual cells. We discover that the relaxation time scale of the
population to its steady state is determined by the distribution of single-cell
growth rates and is surprisingly independent of details of the division process
such as the noise in the timing of division and the mechanism of cell-size
regulation. We validate the predictions of our model using existing
experimental data and propose an experimental method to measure single-cell
growth variability by observing how long it takes for the population to reach
its steady state or balanced growth.


  In spite of the considerable progress towards reducing illiteracy rates, many
countries, including developed ones, have encountered difficulty achieving
further reduction in these rates. This is worrying because illiteracy has been
related to numerous health, social, and economic problems. Here, we show that
the spatial patterns of illiteracy in urban systems have several features
analogous to the spread of diseases such as dengue and obesity. Our results
reveal that illiteracy rates are spatially long-range correlated, displaying
non-trivial clustering structures characterized by percolation-like transitions
and fractality. These patterns can be described in the context of percolation
theory of long-range correlated systems at criticality. Together, these results
provide evidence that the illiteracy incidence can be related to a
transmissible process, in which the lack of access to minimal education
propagates in a population in a similar fashion to endemic diseases.


  How diffusion impacts on ecological dynamics under the Allee effect and
spatial constraints? That is the question we address. Employing a microscopic
minimal model in a metapopulation (without imposing nonlinear birth and death
rates) we evince --- both numerically and analitically --- the emergence of an
optimal diffusion that maximises the survival probability. Even though, at
first such result seems counter-intuitive, it has empirical support from recent
experiments with engineered bacteria. Moreover, we show that this optimal
diffusion disappears for loose spatial constraints.


  Accurate prediction of protein stability changes upon single-site variations
(DDG) is important for protein design, as well as our understanding of the
mechanism of genetic diseases. The performance of high-throughput computational
methods to this end is evaluated mostly based on the Pearson correlation
coefficient between predicted and observed data, assuming that the upper bound
would be 1 (perfect correlation). However, the performance of these predictors
can be limited by the distribution and noise of the experimental data. Here we
estimate, for the first time, a theoretical upper-bound to the DDG prediction
performances imposed by the intrinsic structure of currently available DDG
data. Given a set of measured DDG protein variations, the theoretically best
predictor is estimated based on its similarity to another set of experimentally
determined DDG values. We investigate the correlation between pairs of measured
DDG variations, where one is used as a predictor for the other. We analytically
derive an upper bound to the Pearson correlation as a function of the noise and
distribution of the DDG data. We also evaluate the available datasets to
highlight the effect of the noise in conjunction with DDG distribution. We
conclude that the upper bound is a function of both uncertainty and spread of
the DDG values, and that with current data the best performance should be
between 0.7-0.8, depending on the dataset used; higher Pearson correlations
might be indicative of overtraining. It also follows that comparisons of
predictors using different datasets are inherently misleading.


  In this study, we introduced a new unit, named "protein token", as a dynamic
protein structural unit for protein-protein interactions. Unlike the
conventional structural units, protein token is not based on the sequential or
spatial arrangement of residues, but comprises remote residues involved in
cooperative conformational changes during protein interactions. Application of
protein token on Ras GTPases revealed various tokens present in the
superfamily. Distinct token combinations were found in H-Ras interacting with
its various regulators and effectors, directing to a possible clue for the
multiplexer property of Ras superfamily. Thus, this protein token theory may
provide a new approach to study protein-protein interactions in broad
applications.


  Our paper investigates distributions of exposed and infectious time periods
in an epidemic model and how applying a disease control strategy affects the
model's accuracy. While ordinary differential equations are widely used for
their simplicity, they incorporate an exponential distribution for time spent
exposed or infectious. This allows for a high probability of unrealistically
short exposed and infectious time periods. We propose that caution must be
taken when applying intervention methods to basic models in order to avoid
inaccurate predictions. Delay differential equations, which use a delta
distribution for exposed and infectious periods, can provide better realism but
are more difficult to use and analyze. We introduce a multi-infected
compartment model to interpolate between an ODE model with exponential
distributions and a DDE model with delta distributions in order to investigate
the effect these distributions have on the dynamics of the system when an
intervention method is also included. Using steady state stability and
bifurcation analysis, this paper considers when simpler infectious disease
models can be used versus when more realistic time periods must be
incorporated. We find that the placement of control measures on subpopulations
and the length of the time delay impacts the accuracy of the simpler models.


  Extracting genetic information from a full range of sequencing data is
important for understanding diseases. We propose a novel method to effectively
explore the landscape of genetic mutations and aggregate them to predict cancer
type. We used multinomial logistic regression, nonsmooth non-negative matrix
factorization (nsNMF), and support vector machine (SVM) to utilize the full
range of sequencing data, aiming at better aggregating genetic mutations and
improving their power in predicting cancer types. Specifically, we introduced a
classifier to distinguish cancer types using somatic mutations obtained from
whole-exome sequencing data. Mutations were identified from multiple cancers
and scored using SIFT, PP2, and CADD, and grouped at the individual gene level.
The nsNMF was then applied to reduce dimensionality and to obtain coefficient
and basis matrices. A feature matrix was derived from the obtained matrices to
train a classifier for cancer type classification with the SVM model. We have
demonstrated that the classifier was able to distinguish the cancer types with
reasonable accuracy. In five-fold cross-validations using mutation counts as
features, the average prediction accuracy was 77.1% (SEM=0.1%), significantly
outperforming baselines and outperforming models using mutation scores as
features. Using the factor matrices derived from the nsNMF, we identified
multiple genes and pathways that are significantly associated with each cancer
type. This study presents a generic and complete pipeline to study the
associations between somatic mutations and cancers. The discovered genes and
pathways associated with each cancer type can lead to biological insights. The
proposed method can be adapted to other studies for disease classification and
pathway discovery.


  Symmetry breaking, a central principle of physics, has been hailed as the
driver of self-organization in biological systems in general and biogenesis of
cellular organelles in particular, but the molecular mechanisms of symmetry
breaking only begin to become understood. Centrioles, the structural cores of
centrosomes and cilia, must duplicate every cell cycle to ensure their faithful
inheritance through cellular divisions. Work in model organisms identified
conserved proteins required for centriole duplication and found that altering
their abundance affects centriole number. However, the biophysical principles
that ensure that, under physiological conditions, only a single procentriole is
produced on each mother centriole remain enigmatic. Here we propose a
mechanistic biophysical model for the initiation of procentriole formation in
mammalian cells. The model faithfully recapitulates the experimentally observed
transition from PLK4 uniformly distributed around the mother centriole, the
"ring", to a unique PLK4 focus, the "spot", that triggers the assembly of a new
procentriole. This symmetry breaking requires a dual positive feedback based on
autocatalytic activation of PLK4 and enhanced centriolar anchoring of PLK4-STIL
complexes by phosphorylated STIL. We find that, contrary to previous proposals,
in situ degradation of active PLK4 is insufficient to break symmetry. Instead,
the model predicts that competition between transient PLK4 activity maxima for
PLK4-STIL complexes explains both the instability of the PLK4 ring and
formation of the unique PLK4 spot. In the model, strong competition at
physiologically normal parameters robustly produces a single procentriole,
while increasing overexpression of PLK4 and STIL weakens the competition and
causes progressive addition of procentrioles in agreement with experimental
observations.


  We study a minimal model for the growth of a phenotypically heterogeneous
population of cells subject to a fluctuating environment in which they can
replicate (by exploiting available resources) and modify their phenotype within
a given landscape (thereby exploring novel configurations). The model displays
an exploration-exploitation trade-off whose specifics depend on the statistics
of the environment. Most notably, the phenotypic distribution corresponding to
maximum population fitness (i.e. growth rate) requires a non-zero exploration
rate when the magnitude of environmental fluctuations changes randomly over
time, while a purely exploitative strategy turns out to be optimal in two-state
environments, independently of the statistics of switching times. We obtain
analytical insight into the limiting cases of very fast and very slow
exploration rates by directly linking population growth to the features of the
environment.


  The flexibility of short DNA chains is investigated via computation of the
average correlation function between dimers which defines the persistence
length. Path integration techniques have been applied to confine the phase
space available to base pair fluctuations and derive the partition function.
The apparent persistence lengths of a set of short chains have been computed as
a function of the twist conformation both in the over-twisted and the untwisted
regimes, whereby the equilibrium twist is selected by free energy minimization.
The obtained values are significantly lower than those generally attributed to
kilo-base long DNA. This points to an intrinsic helix flexibility at short
length scales, arising from large fluctuational effects and local bending, in
line with recent experimental indications. The interplay between helical
untwisting and persistence length has been discussed for a heterogeneous
fragment by weighing the effects of the sequence specificities through the
non-linear stacking potential.


  The unprecedented size of the human population, along with its associated
economic activities, have an ever increasing impact on global environments.
Across the world, countries are concerned about the growing resource
consumption and the capacity of ecosystems to provide them. To effectively
conserve biodiversity, it is essential to make indicators and knowledge openly
available to decision-makers in ways that they can effectively use them. The
development and deployment of mechanisms to produce these indicators depend on
having access to trustworthy data from field surveys and automated sensors,
biological collections, molecular data, and historic academic literature. The
transformation of this raw data into synthesized information that is fit for
use requires going through many refinement steps. The methodologies and
techniques used to manage and analyze this data comprise an area often called
biodiversity informatics (or e-Biodiversity). Biodiversity data follows a life
cycle consisting of planning, collection, certification, description,
preservation, discovery, integration, and analysis. Researchers, whether
producers or consumers of biodiversity data, will likely perform activities
related to at least one of these steps. This article explores each stage of the
life cycle of biodiversity data, discussing its methodologies, tools, and
challenges.


  Evolutionary game dynamics in structured populations are strongly affected by
updating rules. Previous studies usually focus on imitation-based rules, which
rely on payoff information of social peers. Recent behavioral experiments
suggest that whether individuals use such social information for strategy
updating may be crucial to the outcomes of social interactions. This hints at
the importance of considering updating rules without dependence on social
peers' payoff information, which, however, is rarely investigated. Here, we
study aspiration-based self-evaluation rules, with which individuals
self-assess the performance of strategies by comparing own payoffs with an
imaginary value they aspire, called the aspiration level. We explore the fate
of strategies on population structures represented by graphs or networks. Under
weak selection, we analytically derive the condition for strategy dominance,
which is found to coincide with the classical condition of risk-dominance. This
condition holds for all networks and all distributions of aspiration levels,
and for individualized ways of self-evaluation. Our condition can be
intuitively interpreted: one strategy prevails over the other if the strategy
brings more satisfaction to individuals than the other does. Our work thus
sheds light on the intrinsic difference between evolutionary dynamics induced
by aspiration-based and imitation-based rules.


  Ciliopathies are a class of human diseases marked by dysfunction of the
cellular organelle, cilia. While many of the molecular components that make up
cilia have been identified and studied, comparatively little is understood
about the transcriptional regulation of genes encoding these components. The
conserved transcription factor Regulatory Factor X (RFX)/DAF-19, which acts
through binding to the cis-regulatory motif known as X-box, has been shown to
regulate ciliary genes in many animals from Caenorhabditis elegans to humans.
However, accumulating evidence suggests that RFX is unable to initiate
transcription on its own. Therefore, other factors and cis-regulatory elements
are likely required. One such element, a DNA motif called the C-box, has
recently been identified in C. elegans. It is still unclear if the X-box and
C-boxes are the only regulatory elements involved and how they interact. To
this end, I analyzed the transcriptional regulation of dyf-5, the C. elegans
ortholog of the human ciliopathy gene Male-Associated Kinase (MAK). Using
computational methods, I was able to confirm the presence of the previously
reported X-box and C-boxes as well as identifying an additional C-box. By
sequentially mutating each of the identified motifs, I identified the role each
potential motif plays in transcriptional regulation of dyf-5. My results showed
that only the X-box and the three C-boxes are necessary and are sufficient to
drive transcription, with the X-box and the centre C-box being the major
contributors and the other two C-boxes enhancing expression. This study
advances the knowledge of gene regulation in general and will further our
understanding of ciliopathies and the mutations that cause them.


  Experimental evolution has yielded surprising insights into human history and
evolution by shedding light on the roles of chance and contingency in history
and evolution, and on the deep evolutionary roots of cooperation, conflict and
kin discrimination. We argue that an interesting research direction would be to
develop computational and experimental systems for studying evolutionary
processes that involve multiple layers of inheritance (such as genes,
epigenetic inheritance, language, and culture) and feedbacks (such as
gene-culture coevolution and mate choice) as well as open-ended niche
construction---all of which are important in human history and evolution. Such
systems would also be a clear way to motivate evolution and computation to
scholars and students across diverse cultural and socioeconomic backgrounds, as
well as to scholars and students in the social sciences and humanities. In
principle, computational models of cultural evolution could be compared to
data, given that large-scale datasets already exist for tracking cultural
change in real-time. Thus, experimental evolution, as a laboratory and
computational science, is poised to grow as an educational tool for people to
question and study where we come from, why we believe what we believe, and
where we as a species may be headed.


  Non-invasive light scattering methods provide data on biological
macromolecules (i.e. proteins, nucleic acids, as well as assemblies and larger
entities composed of them) that are complementary with those of size exclusion
chromatography, gel electrophoresis, analytical ultracentrifugation and mass
spectrometry methods. Static light scattering measurements are useful to
determine the mass of macromolecules and to monitor aggregation phenomena.
Dynamic light scattering measurements are suitable for the quality control and
to assess sample homogeneity, to determine particle size, examine the effect of
physical and chemical treatments, probe the binding of ligands, and study
interactions between macromolecules.


  Fishing activities have broad impacts that affect, although not exclusively,
the targeted stocks. These impacts affect predators and prey of the harvested
species, as well as the whole ecosystem it inhabits. Ecosystem models can be
used to study the interactions that occur within a system, including those
between different organisms and those between fisheries and targeted species.
Trophic web models like Ecopath with Ecosim (EwE) can handle fishing fleets as
a top predator, with top-down impact on harvested organisms. The aim of this
study was to better understand the Icelandic marine ecosystem and the
interactions within. This was done by constructing an EwE model of Icelandic
waters. The model was run from 1984 to 2013 and was fitted to time series of
biomass estimates, landings data and mean annual temperature. The final model
was chosen by selecting the model with the lowest Akaike information criterion.
A skill assessment was performed using the Pearson's correlation coefficient,
the coefficient of determination, the modelling efficiency and the reliability
index to evaluate the model performance. The model performed satisfactorily
when simulating previously estimated biomass and known landings. Most of the
groups with time series were estimated to have top-down control over their
prey. These are harvested species with direct and/or indirect links to lower
trophic levels and future fishing policies should take this into account. This
model could be used as a tool to investigate how such policies could impact the
marine ecosystem in Icelandic waters.


  In this work we have theoretically investigated how the action potential
generation and its associated intrinsic properties are affected in presence of
ion channel blockers by adapting Gillepie's stochastic simulation technique on
a very basic neuron of Hodgkin-Huxley type. With a simple extension of the
Hodgkin-Huxley Markov model we have mainly investigated three types of drug
blocking mechanisms and showed that the major experimental and physiological
observations such as ionic currents, spiking frequency trends, change in action
potential shape and duration, altered gating dynamics etc due to the presence
of ion channel blockers can be well reproduced. The nature of action potential
termination process in presence of sodium and potassium channel blockers are
distinct and physiologically very different from each other. Channel blockers
have distinct signatures on ionic currents. In presence of only sodium channel
blockers the frequency of action potential generation falls off exponentially
with increasing drug affinity, whereas in contrast, for only potassium channel
blockers initially an enhanced spiking activity of action potential is found
followed by a gradual decrease of the spiking frequency as the drug affinity
increases. In case of dual type blockers with equal sodium and potassium
channel binding affinity, the spiking frequency passes through maxima and
minima due to the competition between channel number fluctuation and overall
sodium and potassium conductances. We have found that sodium channel blockers
shorten the duration of action potential while the potassium channel blockers
delay it. We have also shown how the ion channel blockers alter the gating
dynamics. Some experimental results of ion channel blocking in diverse systems
have been validated through our site selected binding scheme.


  Biological control, the use of predators and pathogens to control target
pests, is a promising alternative to chemical control. It is hypothesized that
the introduced predators efficacy can be boosted by providing them with an
additional food source. The current literature \cite{SP11} claims that if the
additional food is of sufficient constant quantity and quality then pest
eradication is possible in \emph{finite} time. The purpose of the current
manuscript is to show that to the contrary, pest eradication is \emph{not
possible} in finite time, for \emph{any} quantity and quality of constant
additional food. However for a density dependent quantity of additional food,
we show pest eradication in finite time is indeed possible. Our results have
large scale implications for the effective design of biological control methods
involving additional food.


  In its simplest form, the competitive exclusion principle states that a
number of species competing for a smaller number of resources cannot coexist.
However, it has been observed empirically that in some settings it is possible
to have coexistence. One example is Hutchinson's `paradox of the plankton'.
This is an instance where a large number of phytoplankton species coexist while
competing for a very limited number of resources. Both experimental and
theoretical studies have shown that temporal fluctuations of the environment
can facilitate coexistence for competing species. Hutchinson conjectured that
one can get coexistence because nonequilibrium conditions would make it
possible for different species to be favored by the environment at different
times.
  In this paper we show in various settings how a variable (stochastic)
environment enables a set of competing species limited by a smaller number of
resources or other density dependent factors to coexist. If the environmental
fluctuations are modeled by white noise, and the per-capita growth rates of the
competitors depend linearly on the resources, we prove that there is
competitive exclusion. However, if either the dependence between the growth
rates and the resources is not linear or the white noise term is nonlinear we
show that coexistence on fewer resources than species is possible. Even more
surprisingly, if the temporal environmental variation comes from switching the
environment at random times between a finite number of possible states, it is
possible for all species to coexist even if the growth rates depend linearly on
the resources. We show in an example (a variant of which first appeared in
Benaim and Lobry '16) that, contrary to Hutchinson's explanation, one can
switch between two environments in which the same species is favored and still
get coexistence.


  Our view of the universe of genomic regions harboring various types of
candidate human-specific regulatory sequences (HSRS) has been markedly expanded
in recent years. To infer the evolutionary origins of loci harboring HSRS,
analyses of conservations patterns of 59,732 loci in Modern Humans, Chimpanzee,
Bonobo, Gorilla, Orangutan, Gibbon, and Rhesus genomes have been performed. Two
major evolutionary pathways have been identified comprising thousands of
sequences that were either inherited from extinct common ancestors (ECAs) or
created de novo in humans after human/chimpanzee split. Thousands of HSRS
appear inherited from ECAs yet bypassed genomes of our closest evolutionary
relatives, presumably due to the incomplete lineage sorting and/or
species-specific loss or regulatory DNA. The bypassing pattern is prominent for
HSRS associated with development and functions of human brain. Common genomic
loci that may contributed to speciation during evolution of Great Apes comprise
248 insertions sites of African Great Ape-specific retrovirus PtERV1 (45.9%; p
= 1.03E-44) intersecting regions harboring 442 HSRS, which are enriched for
HSRS associated with human-specific (HS) changes of gene expression in cerebral
organoids. Among non-human primates (NHP), most significant fractions of
candidate HSRS associated with HS expression changes in both excitatory neurons
(347 loci; 67%) and radial glia (683 loci; 72%) are highly conserved in Gorilla
genome. Modern Humans acquired unique combinations of regulatory sequences
highly conserved in distinct species of six NHP separated by 30 million years
of evolution. Concurrently, this unique mosaic of regulatory sequences
inherited from ECAs was supplemented with 12,486 created de novo HSRS. These
observations support the model of complex continuous speciation process during
evolution of Great Apes that is not likely to occur as an instantaneous event.


  In ecology, species can mitigate their extinction risks in uncertain
environments by diversifying individual phenotypes. This observation is
quantified by the theory of bet-hedging, which provides a reason for the degree
of phenotypic diversity observed even in clonal populations. The theory of
bet-hedging in well-mixed populations is rather well developed. However, many
species underwent range expansions during their evolutionary history, and the
importance of phenotypic diversity in such scenarios still needs to be
understood. In this paper, we develop a theory of bet-hedging for populations
colonizing new, unknown environments that fluctuate either in space or time. In
this case, we find that bet-hedging is a more favorable strategy than in
well-mixed populations. For slow rates of variation, temporal and spatial
fluctuations lead to different outcomes. In spatially fluctuating environments,
bet-hedging is favored compared to temporally fluctuating environments. In the
limit of frequent environmental variation, no opportunity for bet-hedging
exists, regardless of the nature of the environmental fluctuations. For the
same model, bet-hedging is never an advantageous strategy in the well-mixed
case, supporting the view that range expansions strongly promote
diversification. These conclusions are robust against stochasticity induced by
finite population sizes. Our findings shed light on the importance of
phenotypic heterogeneity in range expansions, paving the way to novel
approaches to understand how biodiversity emerges and is maintained.


  Several experiments show that the three dimensional (3D) organization of
chromosomes affects genetic processes such as transcription and gene
regulation. To better understand this connection, researchers developed the
Hi-C method that is able to detect the pairwise physical contacts of all
chromosomal loci. The Hi-C data show that chromosomes are composed of 3D
compartments that range over a variety of scales. However, it is challenging to
systematically detect these cross-scale structures. Most studies have therefore
designed methods for specific scales to study foremost topologically associated
domains (TADs) and A/B compartments. To go beyond this limitation, we tailor a
network community detection method that finds communities in compact fractal
globule polymer systems. Our method allows us to continuously scan through all
scales with a single resolution parameter. We found: (i) polymer segments
belonging to the same 3D community do not have to be in consecutive order along
the polymer chain. In other words, several TADs may belong to the same 3D
community. (ii) CTCF proteins---a loop-stabilizing protein that is ascribed a
big role in TAD formation---are well correlated with community borders only at
one level of organization. (iii) TADs and A/B compartments are traditionally
treated as two weakly related 3D structures and detected with different
algorithms. With our method, we detect both by simply adjusting the resolution
parameter. We therefore argue that they represent two specific levels of a
continuous spectrum 3D communities, rather than seeing them as different
structural entities.


  Dengue is a vector-borne viral disease increasing dramatically over the past
years due to improvement in human mobility. The movement of host individuals
between and within the patches are captured via a residence-time matrix. A
system of ordinary differential equations (ODEs) modeling the spatial spread of
disease among the multiple patches is used to create a system of stochastic
differential equations (SDEs). Numerical solutions of the system of SDEs are
compared with the deterministic solutions obtained via ODEs.


  Computational identification of promoters is notoriously difficult as human
genes often have unique promoter sequences that provide regulation of
transcription and interaction with transcription initiation complex. While
there are many attempts to develop computational promoter identification
methods, we have no reliable tool to analyze long genomic sequences. In this
work we further develop our deep learning approach that was relatively
successful to discriminate short promoter and non-promoter sequences. Instead
of focusing on the classification accuracy, in this work we predict the exact
positions of the TSS inside the genomic sequences testing every possible
location. We studied human promoters to find effective regions for
discrimination and built corresponding deep learning models. These models use
adaptively constructed negative set which iteratively improves the models
discriminative ability. The developed promoter identification models
significantly outperform the previously developed promoter prediction programs
by considerably reducing the number of false positive predictions. The best
model we have built has recall 0.76, precision 0.77 and MCC 0.76, while the
next best tool FPROM achieved precision 0.48 and MCC 0.60 for the recall of
0.75. Our method is available at http://www.cbrc.kaust.edu.sa/PromID/.


  We state and prove a theorem on the partitioning of a randomly selected large
population into stationary and non-stationary components by using a property of
stationary population identity. Applications of this theorem for practical
purposes is summarized at the end.


  We present a novel method for real-time monitoring and kinetic analysis of
fluorescence in situ hybridization (FISH). We implement the method using a
vertical microfluidic probe containing a microstructure designed for rapid
switching between a probe solution and a non-fluorescent imaging buffer. The
FISH signal is monitored in real time during the imaging buffer wash, during
which signal associated with unbound probes is removed. We provide a
theoretical description of the method as well as a demonstration of its
applicability using a model system of centromeric probes (Cen17). We
demonstrate the applicability of the method for the characterization of FISH
kinetics under conditions of varying probe concentration, destabilizing agent
(formamide) content, volume exclusion agent (dextran sulfate) content, and
ionic strength. We show that our method can be used to investigate the effect
of each of these variables and provide insight into processes affecting in situ
hybridization, facilitating the design of new assays.


  We investigate dynamics of a kinetic model of inhibitory autoregulation as
exemplified when a protein inhibits its own production by interfering with its
messenger RNA, known in molecular biology as translational autoregulation. We
first show how linear models without feedback set the stage with a
nonequilibrium steady state that constitutes the target of the regulation.
However, regulation in the simple linear model is far from optimal. The
negative feedback mechanism whereby the protein "jams" the mRNA greatly
enhances the effectiveness of the control, with response to perturbation that
is targeted, rapid, and metabolically efficient. Understanding the full
dynamics of the system phase space is essential to understanding the
autoregulation process.


  Changes in lattice structure across sub-regions of protein crystals are
challenging to assess when relying on whole crystal measurements. Because of
this difficulty, macromolecular structure determination from protein micro and
nano crystals requires assumptions of bulk crystallinity and domain block
substructure. To evaluate the fidelity of these assumptions in protein
nanocrystals we map lattice structure across micron size areas of cryogenically
preserved three-dimensional peptide crystals using a nano-focused electron
beam. This approach produces diffraction from as few as 1,500 molecules in a
crystal, is sensitive to crystal thickness and three-dimensional lattice
orientation. Real-space maps reconstructed from unsupervised classification of
diffraction patterns across a crystal reveal regions of crystal order/disorder
and three-dimensional lattice reorientation on a 20nm scale. The lattice
nano-ripples observed in micron-sized macromolecular crystals provide a direct
view of their plasticity. Knowledge of these features is a first step to
understanding crystalline macromolecular self-assembly and improving the
determination of structures from protein nano and microcrystals from single or
serial crystal diffraction.


  Even though the persistence length $L_P$ of double-stranded DNA plays a
pivotal role in cell biology and nanotechnologies, its dependence on ionic
strength $I$ lacks a consensual description. Using a high-throughput
single-molecule technique and statistical physics modeling, we measure $L_P$ in
presence of monovalent (Li$^+$, Na$^+$, K$^+$) and divalent (Mg$^{2+}$,
Ca$^{2+}$) metallic and alkyl ammonium ions, over a large range 0.5 mM $\leq
I\leq 5$ M. We show that linear Debye-H\"uckel-type theories do not describe
even part of these data. By contrast, the Netz-Orland and Trizac-Shen formulas,
two approximate theories including non-linear electrostatic effects and the
finite DNA radius, fit our data with divalent and monovalent ions,
respectively, over the whole $I$ range. Furthermore the metallic ion type does
not influence $L_P(I)$, in contrast to alkyl ammonium monovalent ions at high
$I$.


  Evolutionary and ecosystem dynamics are often treated as different processes
--operating at separate timescales-- even if evidence reveals that rapid
evolutionary changes can feed back into ecological interactions. A recent
long-term field experiment has explicitly shown that communities of competing
plant species can experience very fast phenotypic diversification, and that
this gives rise to enhanced complementarity in resource exploitation and to
enlarged ecosystem-level productivity. Here, we build on progress made in
recent years in the integration of eco-evolutionary dynamics, and present a
computational approach aimed at describing these empirical findings in detail.
In particular we model a community of organisms of different but similar
species evolving in time through mechanisms of birth, competition, sexual
reproduction, descent with modification, and death. Based on simple rules, this
model provides a rationalization for the emergence of rapid phenotypic
diversification in species-rich communities. Furthermore, it also leads to
non-trivial predictions about long-term phenotypic change and ecological
interactions. Our results illustrate that the presence of highly specialized,
non-competing species leads to very stable communities and reveals that
phenotypically equivalent species occupying the same niche may emerge and
coexist for very long times. Thus, the framework presented here provides a
simple approach --complementing existing theories, but specifically devised to
account for the specificities of the recent empirical findings for plant
communities-- to explain the collective emergence of diversification at a
community level, and paves the way to further scrutinize the intimate
entanglement of ecological and evolutionary processes, especially in
species-rich communities.


  Motivated by the occurrence of a moderately nearby supernova near the
beginning of the Pleistocene, possibly as part of a long-term series beginning
in the Miocene, we investigate whether nitrate rainout resulting from the
atmospheric ionization of enhanced cosmic ray flux could have, through its
fertilizer effect, initiated carbon dioxide drawdown. Such a drawdown could
possibly reduce the greenhouse effect and induce the climate change that led to
the Pleistocene glaciations. We estimate that the nitrogen flux enhancement
onto the surface from an event at 50 pc would be of order 10%, probably too
small for dramatic changes. We estimate deposition of iron (another potential
fertilizer) and find it is also too small to be significant. There are also
competing effects of opposite sign, including muon irradiation and reduction in
photosynthetic yield caused by UV increase from stratospheric ozone layer
depletion, leading to an ambiguous result. However, if the atmospheric
ionization induces a large increase in the frequency of lightning, as argued
elsewhere, the amount of nitrate synthesis should be much larger, dominate over
the other effects, and induce the climate change. More work needs to be done to
clarify effects on lightning frequency.


  Background: High-throughput techniques bring novel tools but also statistical
challenges to genomic research. Identifying genes with differential expression
between different species is an effective way to discover evolutionarily
conserved transcriptional responses. To remove systematic variation between
different species for a fair comparison, the normalization procedure serves as
a crucial pre-processing step that adjusts for the varying sample sequencing
depths and other confounding technical effects.
  Results: In this paper, we propose a scale based normalization (SCBN) method
by taking into account the available knowledge of conserved orthologous genes
and hypothesis testing framework. Considering the different gene lengths and
unmapped genes between different species, we formulate the problem from the
perspective of hypothesis testing and search for the optimal scaling factor
that minimizes the deviation between the empirical and nominal type I errors.
Conclusions: Simulation studies show that the proposed method performs
significantly better than the existing competitor in a wide range of settings.
An RNA-seq dataset of different species is also analyzed and it coincides with
the conclusion that the proposed method outperforms the existing method. For
practical applications, we have also developed an R package named "SCBN" and
the software is available at
http://www.bioconductor.org/packages/devel/bioc/html/SCBN.html.


  Cell-cell communication is often achieved by secreted signaling molecules
that bind membrane-bound receptors. A common class of such receptors are
G-protein coupled receptors, where extracellular binding induces changes on the
membrane affinity near the receptor for certain cytosolic proteins, effectively
altering their chemical potential. We analyze the minimum-dissipation schedules
for dynamically changing chemical potential to induce steady-state changes in
protein copy-number distributions, and illustrate with analytic solutions for
linear chemical reaction networks. Protocols that change chemical potential on
biologically relevant timescales are experimentally accessible using
optogenetic manipulations, and our framework provides non-trivial predictions
about functional dynamical cell-cell interactions.


  We study the folding of RNA secondary structures with quenched sequence
randomness by means of the constrained annealing method. A thermodynamic phase
transition is induced by including the conformational weight of loop
structures. In addition to the expected melting at high temperature, a cold
melting transition appears. Our results suggest that the cold denaturation of
RNA found experimentally is, in fact, a continuous phase transition triggered
by quenched sequence disorder. We calculate both hot and cold melting critical
temperatures for the competing energy scenario between favorable and
unfavorable base pairs and present a phase diagram as a function of the loop
exponent and temperature.


  The rate of biological evolution depends on the fixation probability and on
the fixation time of new mutants. Intensive research has focused on identifying
population structures that augment the fixation probability of advantageous
mutants. But these `amplifiers of natural selection' typically increase
fixation time. Here we study population structures that achieve a trade-off
between high fixation probability and short fixation time. First, we show that
no amplifiers can have asymptotically lower absorption time than the well-mixed
population. Then we design population structures that substantially augment the
fixation probability with just a minor increase in fixation time. Finally, we
show that those structures enable higher effective rate of evolution than the
well-mixed population provided that the rate of generating advantageous mutants
is relatively low. Our work sheds light on how population structure affects the
rate of evolution. Moreover, our structures could be useful for lab-based,
medical or industrial applications of evolutionary optimization.


  We show that discrete distributions on the $d$-dimensional non-negative
integer lattice can be approximated arbitrarily well via the marginals of
stationary distributions for various classes of stochastic chemical reaction
networks. We begin by providing a class of detailed balanced networks and prove
that they can approximate any discrete distribution to any desired accuracy.
However, these detailed balanced constructions rely on the ability to
initialize a system precisely, and are therefore susceptible to perturbations
in the initial conditions. We therefore provide another construction based on
the ability to approximate point mass distributions and prove that this
construction is capable of approximating arbitrary discrete distributions for
any choice of initial condition. In particular, the developed models are
ergodic, so their limit distributions are robust to a finite number of
perturbations over time in the counts of molecules.


  Living cells communicate information about physiological conditions by
producing signaling molecules in a specific timed manner. Different conditions
can result in the same total amount of a signaling molecule, differing only in
the pattern of the molecular concentration over time. Such temporally coded
information can be completely invisible to even state-of-the-art molecular
sensors with high chemical specificity that respond only to the total amount of
the signaling molecule. Here, we demonstrate design principles for circuits
with temporal specificity, that is, molecular circuits that respond to specific
temporal patterns in a molecular concentration. We consider pulsatile patterns
in a molecular concentration characterized by three fundamental temporal
features - time period, duty fraction and number of pulses. We develop circuits
that respond to each one of these features while being insensitive to the
others. We demonstrate our design principles using abstract Chemical Reaction
Networks and with explicit simulations of DNA strand displacement reactions. In
this way, our work develops building blocks for temporal pattern recognition
through molecular computation.


  Purpose: Pulsed-laser irradiation of light-absorbing gold nanoparticles
(AuNPs) attached to cells transiently increases cell membrane permeability for
targeted molecule delivery. Here, we targeted EGFR on the ovarian carcinoma
cell line OVCAR-3 with AuNPs. In order to optimize membrane permeability and to
demonstrate molecule delivery into adherent OVCAR-3 cells, we systematically
investigated different experimental conditions. Materials and methods: AuNPs
(30 nm) were functionalized by conjugation of the antibody cetuximab against
EGFR. Selective binding of the particles was demonstrated by silver staining,
multiphoton imaging, and fluorescence-lifetime imaging. After laser
irradiation, membrane permeability of OVCAR-3 cells was studied under different
conditions of AuNP concentration, cell-incubation medium, and cell - AuNP
incubation time. Membrane permeability and cell viability were evaluated by
flow cytometry, measuring propidium iodide and fluorescein isothiocyanate -
dextran uptake. Results: Adherently growing OVCAR-3 cells can be effectively
targeted with EGFR-AuNP. Laser irradiation led to successful permeabilization,
and 150 kDa dextran was successfully delivered into cells with about 70%
efficiency. Conclusion: Antibody-targeted and laser-irradiated AuNPs can be
used to deliver molecules into adherent cells. Efficacy depends not only on
laser parameters but also on AuNP:cell ratio, cell-incubation medium, and cell
- AuNP incubation time.


  Plant diseases are responsible for substantial and sometimes devastating
economic and societal costs and thus are a major limiting factor for stable and
sustainable agricultural production. Diseases of crops are particular crippling
in developing countries that are heavily dependent on agriculture for food
security and income. Various techniques have been developed to reduce the
negative impact of plant diseases and eliminate the associated parasites, but
the success of these approaches strongly depends on population awareness and
the degree of engagement with disease control and prevention programs. In this
paper we derive and analyse a mathematical model of mosaic disease of {\it
Jatropha curcas}, an important biofuel plant, with particular emphasis on the
effects of interventions in the form of nutrients and insecticides, whose use
depends on the level of population awareness. Two contributions to disease
awareness are considered in the model: global awareness campaigns, and
awareness from observing infected plants. All steady states of the model are
found, and their stability is analysed in terms of system parameters. We
identify parameter regions associated with eradication of disease, stable
endemic infection, and periodic oscillations in the level of infection.
Analytical results are supported by numerical simulations that illustrate the
behaviour of the model in different dynamical regimes. Implications of
theoretical results for practical implementation of disease control are
discussed.


  Transcriptional repressor CTCF is an important regulator of chromatin 3D
structure, facilitating the formation of topologically associating domains
(TADs). However, its direct effects on gene regulation is less well understood.
Here, we utilize previously published ChIP-seq and RNA-seq data to investigate
the effects of CTCF on alternative splicing of genes with CTCF sites. We
compared the amount of RNA-seq signals in exons upstream and downstream of
binding sites following auxin-induced degradation of CTCF in mouse embryonic
stem cells. We found that changes in gene expression following CTCF depletion
were significant, with a general increase in the presence of upstream exons. We
infer that a possible mechanism by which CTCF binding contributes to
alternative splicing is by causing pauses in the transcription mechanism during
which splicing elements are able to concurrently act on upstream exons already
transcribed into RNA.


  Using the formalism of information theory, we analyze the mechanism of
information transduction in a simple one-step signaling cascade S$\rightarrow$X
representing the gene regulatory network. Approximating the signaling channel
to be Gaussian, we describe the dynamics using Langevin equations. Upon
discretization, we calculate the associated second moments for linear and
nonlinear regulation of the output by the input, which follows the birth-death
process. While mutual information between the input and the output
characterizes the channel capacity, the Fano factor of the output gives a clear
idea of how internal and external fluctuations assemble at the output level. To
quantify the contribution of the present state of the input to predict the
future output, transfer entropy is computed. We find that higher amount of
transfer entropy is accompanied by the greater magnitude of external
fluctuations (quantified by the Fano factor of the output) propagation from the
input to the output. We notice that low input population characterized by the
number of signaling molecules S, which fluctuates in a relatively slower
fashion compared to its downstream (target) species X, is maximally able to
predict (as quantified by transfer entropy) the future state of the output. Our
computations also reveal that with increased linear nature of the input-output
interaction, all three metrics of mutual information, Fano factor and, transfer
entropy achieve relatively larger magnitudes.


  Extinction is the ultimate absorbing state of any stochastic birth-death
process, hence the time to extinction is an important characteristic of any
natural population. Here we consider logistic and logistic-like systems under
the combined effect of demographic and bounded environmental stochasticity.
Three phases are identified: an inactive phase where the mean time to
extinction $T$ increases logarithmically with the initial population size, an
active phase where $T$ grows exponentially with the carrying capacity $N$, and
temporal Griffiths phase, with power-law relationship between $T$ and $N$. The
system supports an exponential phase only when the noise is bounded, in which
case the continuum (diffusion) approximation breaks down within the Griffiths
phase. This breakdown is associated with a crossover between qualitatively
different survival statistics and decline modes. To study the power-law phase
we present a new WKB scheme which is applicable both in the diffusive and in
the non-diffusive regime.


  Understanding the mechanisms responsible for the formation and growth of FtsZ
polymers and their subsequent formation of the $Z$-ring is important for
gaining insight into the cell division in prokaryotic cells. In this work, we
present a minimal stochastic model that qualitatively reproduces {\it in vitro}
observations of polymerization, formation of dynamic contractile ring that is
stable for a long time and depolymerization shown by FtsZ polymer filaments. In
this stochastic model, we explore different mechanisms for ring breaking and
hydrolysis. In addition to hydrolysis, which is known to regulate the dynamics
of other tubulin polymers like microtubules, we find that the presence of the
ring allows for an additional mechanism for regulating the dynamics of FtsZ
polymers. Ring breaking dynamics in the presence of hydrolysis naturally induce
rescue and catastrophe events in this model irrespective of the mechanism of
hydrolysis.


  In this work we propose a method to compute continuous embeddings for kmers
from raw RNA-seq data, without the need for alignment to a reference genome.
The approach uses an RNN to transform kmers of the RNA-seq reads into a 2
dimensional representation that is used to predict abundance of each kmer. We
report that our model captures information of both DNA sequence similarity as
well as DNA sequence abundance in the embedding latent space, that we call the
Latent Transcriptome. We confirm the quality of these vectors by comparing them
to known gene sub-structures and report that the latent space recovers exon
information from raw RNA-Seq data from acute myeloid leukemia patients.
Furthermore we show that this latent space allows the detection of genomic
abnormalities such as translocations as well as patient-specific mutations,
making this representation space both useful for visualization as well as
analysis.


  Temporal resolution of visual information processing is thought to be an
important factor in predator-prey interactions, shaped in the course of
evolution by animals' ecology. Here I show that light can be considered to have
a dual role of a source of information, which guides motor actions, and an
environmental feedback for those actions. I consequently show how temporal
perception might depend on behavioral adaptations realized by the nervous
system. I propose an underlying mechanism of synaptic clock, with every synapse
having its characteristic time unit, determined by the persistence of memory
traces of synaptic inputs, which is used by the synapse to tell time. The
present theory offers a testable framework, which may account for numerous
experimental findings, including the interspecies variation in temporal
resolution and the properties of subjective time perception, specifically the
variable speed of perceived time passage, depending on emotional and
attentional states or tasks performed.


  We answer two questions raised by Bryant, Francis and Steel in their work on
consensus methods in phylogenetics. Consensus methods apply to every practical
instance where it is desired to aggregate a set of given phylogenetic trees
(say, gene evolution trees) into a resulting, "consensus" tree (say, a species
tree). Various stability criteria have been explored in this context, seeking
to model desirable consistency properties of consensus methods as the
experimental data are updated (e.g., more taxa, or more trees, are mapped).
However, such stability conditions can be incompatible with some basic
regularity properties that are widely accepted to be essential in any
meaningful consensus method. Here, we prove that such an incompatibility does
arise in the case of extension stability on binary trees and in the case of
associative stability. Our methods combine general theoretical considerations
with the use of computer programs tailored to the given stability requirements.


  We propose a model of rumor spreading in which susceptible, but skeptically
oriented individuals may oppose the rumor. Resistance may be implemented either
by skeptical activists trying to convince spreaders to stop their activity,
becoming stiflers or, passively (non-reactive) as a consequence, for example,
of fact-checking. Interestingly, these two mechanisms, when combined, are
similar to the (assumed) spreading of a fictitious zombie outbreak, where
survivors actively target infected people. We analyse the well-mixed
(mean-field) description and obtain the conditions for rumors (zombies) to
spread through the whole population. The results show that when the skepticism
is strong enough, the model predicts the coexistence of two fixed points (such
bistability may be related to polarized situations), with the fate of rumors
depending on the initial exposure to it.


  Spatio temporal data are more ubiquitous and richer than even before and the
availability of such data poses great challenges in data analytics. Ecological
facilitation, the positive effect of density of individuals on the individual's
survival across a stress gradient, is a complex phenomenon. A large number of
tree individuals coupled with soil moisture, temperature, and water stress data
across a long temporal period were followed. Data driven analysis in the
absence of hypothesis was performed. Information theoretic analysis of multiple
statistical models was employed in order to quantify the best data-driven index
of vegetation density and spatial scale of interactions. Sequentially, tree
survival was quantified as a function of the size of the individual, vegetation
density, and time at the optimal spatial interaction scale. Land surface
temperature and soil moisture were also statistically explained by tree size,
density, and time. Results indicated that in space both facilitation and
competition coexist in the same ecosystem and the sign and magnitude of this
depend on the spatial scale. Overall, within the optimal data driven spatial
scale, tree survival was best explained by the interaction between density and
year, sifting overall from facilitation to competition through time. However,
small sized trees were always facilitated by increased densities, while large
sized trees had either negative or no density effects. Tree size was more
important predictor than density in survival and this has implications for
nature based solutions: maintaining large tree individuals or planting species
that can become large-sized can safeguard against tree less areas by promoting
survival at long time periods through harsh environmental conditions. Large
trees had also a significant effect in moderating land surface temperature.


  Computation biology helps to understand all processes in organisms from
interaction of molecules to complex functions of whole organs. Therefore, there
is a need for mathematical methods and models that deliver logical explanations
in a reasonable time. For the last few years there has been a growing interest
in biological theory connected to finite fields: the algebraic modeling tools
used up to now are based on Gr\"obner bases or Boolean group. Let $n$ variables
representing gene products, changing over the time on $p$ values. A Polynomial
dynamical system (PDS) is a function which has several components, each one is
a polynom with $n$ variables and coefficient in the finite field $Z/pZ$ that
model the evolution of gene products. We propose herein a method using
algebraic separators, which are special polynomials abundantly studied in
effective Galois theory. This approach avoids heavy calculations and provides a
first Polynomial model in linear time.


  We analyze the fate of dynamical systems that consist of two kind of
processes. The first type is supposed to perform a certain function by
processing information at a required high accuracy, which is, however, limited
to less than 100 percent, while the second process serves to maintain the
required precision. Both processes are assumed to be subject to a trade-off
between cost and precision, where the cost have to be paid from renewable but
limited resources. In a discrete map we pursue the time evolution of errors and
determine the conditions under which the fate of the system is either a stable
performance at the desired accuracy, or a deterioration. Deterioration may be
realized either as an accumulation of errors or a decline of resources when
they are all absorbed for maintenance. We point to possible implications for
living organisms and their perspectives to avoid an accumulation of errors in
the course of time.


  Blue mussels (Mytilus edulis) are an important keystone species that have
been declining in the Gulf of Maine. This could be attributed to a variety of
complex factors such as indirect effects due to invasion by epibionts, which
remains unexplored mathematically. Based on classical optimal foraging theory
and anti-fouling defense mechanisms of mussels, we derive an ODE model for
crab-mussel interactions in the presence of an invasive epibiont, Didemnum
vexillum. The dynamical analysis leads to results on stability, global
boundedness and bifurcations of the model. Next, via optimal control methods we
predict various ecological outcomes. Our results have key implications for
preserving mussel populations in the advent of invasion by non-native
epibionts. In particular they help us understand the changing dynamics of local
predator-prey communities, due to indirect effects that epibionts confer.


  In this work we combing models of disease dynamics and economic production,
and we show the potential implications of this for demonstrating the importance
of savings for buffering an economy during the period of an epidemic.
  Finding an explicit function that relates poverty and the production of a
community is an almost impossible task because of the number of variables and
parameters that should be taken into account. However, studying the dynamics of
an endemic disease in a region that affects its population, and therefore its
ability to work, is an honest approach to understanding this function. We
propose a model, perhaps the simplest, that couples two dynamics, the dynamics
of an endemic disease and the dynamics of a closed economy of products and
goods that the community produces in the epidemic period.
  Some of the results of this study are expected and known in the literature
but some others are not. We highlight three of them: the interdependence that
exists between health and the product of that economy. The majority of the
known results only show the dependence in a single direction. Another of the
most important goals in this work is to show the existence of a poverty
Malthusian trap in economies with low levels of capital investment due to
pandemic disease. And another is to show that there is an optimal level of
savings that maximizes the \emph{per capita} product of the economy, which is
Machiavellian, but has important implications for market efficiency. The model
tells us why saving is important in an aggregate economy despite being affected
by an epidemic.