  We developed simulation methodology to assess eventual therapeutic efficiency
of exogenous multiparametric changes in a four-component cellular system
described by the system of ordinary differential equations. The method is
numerically implemented to simulate the temporal behavior of a cellular system
of multiple myeloma cells. The problem is conceived as an inverse optimization
task where the alternative temporal changes of selected parameters of the
ordinary differential equations represent candidate solutions and the objective
function quantifies the goals of the therapy. The system under study consists
of two main cellular components, tumor cells and their cellular environment,
respectively. The subset of model parameters closely related to the environment
is substituted by exogenous time dependencies - therapeutic pulses combining
continuous functions and discrete parameters subordinated thereafter to the
optimization. Synergistic interaction of temporal parametric changes has been
observed and quantified whereby two or more dynamic parameters show effects
that absent if either parameter is stimulated alone. We expect that the
theoretical insight into unstable tumor growth provided by the sensitivity and
optimization studies could, eventually, help in designing combination
therapies.


  We explore a model of metapopulation genetics which is based on a more
ecologically motivated approach than is frequently used in population genetics.
The size of the population is regulated by competition between individuals,
rather than by artificially imposing a fixed population size. The increased
complexity of the model is managed by employing techniques often used in the
physical sciences, namely exploiting time-scale separation to eliminate fast
variables and then constructing an effective model from the slow modes.
Remarkably, an initial model with 2$\mathcal{D}$ variables, where $\mathcal{D}$
is the number of islands in the metapopulation, can be reduced to a model with
a single variable. We analyze this effective model and show that the
predictions for the probability of fixation of the alleles and the mean time to
fixation agree well with those found from numerical simulations of the original
model.


  Here we develop an effective approach to simplify two-time-scale Markov
chains with infinite state spaces by removal of states with fast leaving rates,
which improves the simplification method of finite Markov chains. We introduce
the concept of fast transition paths and show that the effective transitions of
the reduced chain can be represented as the superposition of the direct
transitions and the indirect transitions via all the fast transition paths.
Furthermore, we apply our simplification approach to the standard Markov model
of single-cell stochastic gene expression and provide a mathematical theory of
random gene expression bursts. We give the precise mathematical conditions for
the bursting kinetics of both mRNAs and proteins. It turns out that random
bursts exactly correspond to the fast transition paths of the Markov model.
This helps us gain a better understanding of the physics behind the bursting
kinetics as an emergent behavior from the fundamental multi-scale biochemical
reaction kinetics of stochastic gene expression.


  The emerging field of high-throughput compartmentalized in vitro evolution is
a promising new approach to protein engineering. In these experiments,
libraries of mutant genotypes are randomly distributed and expressed in
microscopic compartments - droplets of an emulsion. The selection of desirable
variants is performed according to the phenotype of each compartment. The
random partitioning leads to a fraction of compartments receiving more than one
genotype making the whole process a lab implementation of the group selection.
From a practical point of view (where efficient selection is typically sought),
it is important to know the impact of the increase in the mean occupancy of
compartments on the selection efficiency. We carried out a theoretical
investigation of this problem in the context of selection dynamics for an
infinite non-mutating subdivided population that randomly colonizes an infinite
number of patches (compartments) at each reproduction cycle. We derive here an
update equation for any distribution of phenotypes and any value of the mean
occupancy. Using this result, we demonstrate that, for the linear additive
fitness, the best genotype is still selected regardless of the mean occupancy.
Furthermore, the selection process is remarkably resilient to the presence of
multiple genotypes per compartments, and slows down approximately inversely
proportional to the mean occupancy at high values. We extend out results to
more general expressions that cover nonadditive and non-linear fitnesses, as
well non-Poissonian distribution among compartments. Our conclusions may also
apply to natural genetic compartmentalized replicators, such as viruses or
early trans-acting RNA replicators.


  In this paper we investigate the complexity of model selection and model
testing for dynamical systems with toric steady states. Such systems frequently
arise in the study of chemical reaction networks. We do this by formulating
these tasks as a constrained optimization problem in Euclidean space. This
optimization problem is known as a Euclidean distance problem; the complexity
of solving this problem is measured by an invariant called the Euclidean
distance (ED) degree. We determine closed-form expressions for the ED degree of
the steady states of several families of chemical reaction networks with toric
steady states and arbitrarily many reactions. To illustrate the utility of this
work we show how the ED degree can be used as a tool for estimating the
computational cost of solving the model testing and model selection problems.


  Modern time series gene expression and other omics data sets have enabled
unprecedented resolution of the dynamics of cellular processes such as cell
cycle and response to pharmaceutical compounds. In anticipation of the
proliferation of time series data sets in the near future, we use the Hopfield
model, a recurrent neural network based on spin glasses, to model the dynamics
of cell cycle in HeLa (human cervical cancer) and S. cerevisiae cells. We study
some of the rich dynamical properties of these cyclic Hopfield systems,
including the ability of populations of simulated cells to recreate
experimental expression data and the effects of noise on the dynamics. Next, we
use a genetic algorithm to identify sets of genes which, when selectively
inhibited by local external fields representing gene silencing compounds such
as kinase inhibitors, disrupt the encoded cell cycle. We find, for example,
that inhibiting the set of four kinases BRD4, MAPK1, NEK7, and YES1 in HeLa
cells causes simulated cells to accumulate in the M phase. Finally, we suggest
possible improvements and extensions to our model.


  Scaling laws in ecology, intended both as functional relationships among
ecologically-relevant quantities and the probability distributions that
characterize their occurrence, have long attracted the interest of empiricists
and theoreticians. Empirical evidence exists of power laws associated with the
number of species inhabiting an ecosystem, their abundances and traits.
Although their functional form appears to be ubiquitous, empirical scaling
exponents vary with ecosystem type and resource supply rate. The idea that
ecological scaling laws are linked had been entertained before, but the full
extent of macroecological pattern covariations, the role of the constraints
imposed by finite resource supply and a comprehensive empirical verification
are still unexplored. Here, we propose a theoretical scaling framework that
predicts the linkages of several macroecological patterns related to species'
abundances and body sizes. We show that such framework is consistent with the
stationary state statistics of a broad class of resource-limited community
dynamics models, regardless of parametrization and model assumptions. We verify
predicted theoretical covariations by contrasting empirical data and provide
testable hypotheses for yet unexplored patterns. We thus place the observed
variability of ecological scaling exponents into a coherent statistical
framework where patterns in ecology embed constrained fluctuations.


  Autoinducers are small signaling molecules that mediate intercellular
communication in microbial populations and trigger coordinated gene expression
via "quorum sensing". Elucidating the mechanisms that control autoinducer
production is, thus, pertinent to understanding collective microbial behavior,
such as virulence and bioluminescence. Recent experiments have shown a
heterogeneous promoter activity of autoinducer synthase genes, suggesting that
some of the isogenic cells in a population might produce autoinducers, whereas
others might not. However, the mechanism underlying this phenotypic
heterogeneity in quorum-sensing microbial populations has remained elusive. In
our theoretical model, cells synthesize and secrete autoinducers into the
environment, up-regulate their production in this self-shaped environment, and
non-producers replicate faster than producers. We show that the coupling
between ecological and population dynamics through quorum sensing can induce
phenotypic heterogeneity in microbial populations, suggesting an alternative
mechanism to stochastic gene expression in bistable gene regulatory circuits.


  We describe a continuous-time modelling framework for biological population
dynamics that accounts for demographic noise. In the spirit of the methodology
used by statistical physicists, transitions between the states of the system
are caused by individual events while the dynamics are described in terms of
the time-evolution of a probability density function. In general, the
application of the diffusion approximation still leaves a description that is
quite complex. However, in many biological applications one or more of the
processes happen slowly relative to the system's other processes, and the
dynamics can be approximated as occurring within a slow low-dimensional
subspace. We review these time-scale separation arguments and analyse the more
simple stochastic dynamics that result in a number of cases. We stress that it
is important to retain the demographic noise derived in this way, and emphasise
this point by showing that it can alter the direction of selection compared to
the prediction made from an analysis of the corresponding deterministic model.


  The symptoms of many infectious diseases influence their host to withdraw
from social activity limiting their own potential to spread. Successful
transmission therefore requires the onset of infectiousness to coincide with a
time when its host is socially active. Since social activity and infectiousness
are both temporal phenomena, we hypothesize that diseases are most pervasive
when these two processes are synchronized. We consider disease dynamics that
incorporate a behavioral response that effectively shortens the infectious
period of the disease. We apply this model to data collected from face-to-face
social interactions and look specifically at how the duration of the latent
period effects the reachability of the disease. We then simulate the spread of
the model disease on the network to test the robustness of our results.
Diseases with latent periods that synchronize with the temporal social behavior
of people, i.e. latent periods of 24 hours or 7 days, correspond to peaks in
the number of individuals who are potentially at risk of becoming infected. The
effect of this synchronization is present for a range of disease models with
realistic parameters. The relationship between the latent period of an
infectious disease and its pervasiveness is non-linear and depends strongly on
the social context in which the disease is spreading.


  Population genetics models typically consider a fixed population size and a
unique selection coefficient. However, population dynamics inherently generate
noise in numbers of individuals and selection acts on various components of the
individuals' fitness. In plant species with seed banks, the size of both the
above- and below-ground compartments present noise depending on seed production
and the state of the seed bank. We investigate if this noise has consequences
on 1)~the rate of genetic drift, and 2)~the efficacy of selection. We consider
four variants of two-allele Moran-type models defined by combinations of
presence and absence of noise in above-ground and seed bank compartments. Time
scale analysis and dimension reduction methods allow us to reduce the
corresponding Fokker-Planck equation to a one-dimensional diffusion
approximation of a Moran model. We first show that if the above-ground noise
classically affects the rate of genetic drift, below-ground noise reduces the
diversity storage effect of the seed bank. Second, we consider that selection
can act on four different components of the plant fitness: plant or seed death
rate, seed production or seed germination. Our striking result is that the
efficacy of selection for seed death rate or germination rate is reduced by
seed bank noise, whereas selection occurring on plant death rate or seed
production is not affected. We derive the expected site-frequency spectrum
reflecting this heterogeneity in selection efficacy between genes underpinning
different plant fitness components. Our results highlight the importance to
consider the effect of ecological noise to predict the impact of seed banks on
neutral and selective evolution.


  Identifying undocumented or potential future interactions among species is a
challenge facing modern ecologists. Recent link prediction methods rely on
trait data, however large species interaction databases are typically sparse
and covariates are limited to only a fraction of species. On the other hand,
evolutionary relationships, encoded as phylogenetic trees, can act as proxies
for underlying traits and historical patterns of parasite sharing among hosts.
We show that using a network-based conditional model, phylogenetic information
provides strong predictive power in a recently published global database of
host-parasite interactions. By scaling the phylogeny using an evolutionary
model, our method allows for biological interpretation often missing from
latent variable models. To further improve on the phylogeny-only model, we
combine a hierarchical Bayesian latent score framework for bipartite graphs
that accounts for the number of interactions per species with the host
dependence informed by phylogeny. Combining the two information sources yields
significant improvement in predictive accuracy over each of the submodels
alone. As many interaction networks are constructed from presence-only data, we
extend the model by integrating a correction mechanism for missing
interactions, which proves valuable in reducing uncertainty in unobserved
interactions.


  Many complex systems are modular. Such systems can be represented as
"component systems", i.e., sets of elementary components, such as LEGO bricks
in LEGO sets. The bricks found in a LEGO set reflect a target architecture,
which can be built following a set-specific list of instructions. In other
component systems, instead, the underlying functional design and constraints
are not obvious a priori, and their detection is often a challenge of both
scientific and practical importance, requiring a clear understanding of
component statistics. Importantly, some quantitative invariants appear to be
common to many component systems, most notably a common broad distribution of
component abundances, which often resembles the well-known Zipf's law. Such
"laws" affect in a general and non-trivial way the component statistics,
potentially hindering the identification of system-specific functional
constraints or generative processes. Here, we specifically focus on the
statistics of shared components, i.e., the distribution of the number of
components shared by different system-realizations, such as the common bricks
found in different LEGO sets. To account for the effects of component
heterogeneity, we consider a simple null model, which builds
system-realizations by random draws from a universe of possible components.
Under general assumptions on abundance heterogeneity, we provide analytical
estimates of component occurrence, which quantify exhaustively the statistics
of shared components. Surprisingly, this simple null model can positively
explain important features of empirical component-occurrence distributions
obtained from data on bacterial genomes, LEGO sets, and book chapters. Specific
architectural features and functional constraints can be detected from
occurrence patterns as deviations from these null predictions, as we show for
the illustrative case of the "core" genome in bacteria.


  A deep neural network based architecture was constructed to predict amino
acid side chain conformation with unprecedented accuracy. Amino acid side chain
conformation prediction is essential for protein homology modeling and protein
design. Current widely-adopted methods use physics-based energy functions to
evaluate side chain conformation. Here, using a deep neural network
architecture without physics-based assumptions, we have demonstrated that side
chain conformation prediction accuracy can be improved by more than 25%,
especially for aromatic residues compared with current standard methods. More
strikingly, the prediction method presented here is robust enough to identify
individual conformational outliers from high resolution structures in a protein
data bank without providing its structural factors. We envisage that our amino
acid side chain predictor could be used as a quality check step for future
protein structure model validation and many other potential applications such
as side chain assignment in Cryo-electron microscopy, crystallography model
auto-building, protein folding and small molecule ligand docking.


  We apply our statistically deterministic machine learning/clustering
algorithm *K-means (recently developed in https://ssrn.com/abstract=2908286) to
10,656 published exome samples for 32 cancer types. A majority of cancer types
exhibit mutation clustering structure. Our results are in-sample stable. They
are also out-of-sample stable when applied to 1,389 published genome samples
across 14 cancer types. In contrast, we find in- and out-of-sample
instabilities in cancer signatures extracted from exome samples via nonnegative
matrix factorization (NMF), a computationally costly and non-deterministic
method. Extracting stable mutation structures from exome data could have
important implications for speed and cost, which are critical for early-stage
cancer diagnostics such as novel blood-test methods currently in development.


  In response to cold stress, Drosophila Melanogaster increase their expression
of Frost, a candidate gene involved in cold response. The direct role of Frost
in cold tolerance is yet to be determined, and its importance for survival in
cold environments has been questioned. In this study, I attempt to better
understand the molecular machinery of Frost by selecting for its protein in fly
lysate knowing only its RNA has been found. Detection of Frost expression and
subsequently studying its protein will hopefully lead to enhanced comprehension
of cold responses in insects, which is the long-term goal of this research. I
predict that Frost will be expressed in flies that undergo cold stress at 0oC
for 2 hours before recovering at 220C for 3 hours before lysing. Two Western
blots were executed using fly lysate that underwent the above treatment, whose
antibodies were specific to Frost. No bands were seen in any of the lanes
containing treated samples in either of the Immunoblots, indicating that Frost
was either not expressed or not detected in fly lysate.


  Inside cells, cargos such as vesicles and organelles are transported by
molecular motors to their correct locations via active motion on cytoskeletal
tracks and passive, Brownian diffusion. During the transportation of cargos,
motor-cargo complexes (MCC) navigate the confining and crowded environment of
the cytoskeletal network and other macromolecules. Motivated by this, we study
a minimal two-state model of motor-driven cargo transport in confinement and
predict transport properties that can be tested in experiments. We assume that
the motion of the MCC is directly affected by the entropic barrier due to
confinement if it is in the passive, unbound state, but not in the active,
bound state where it moves with a constant bound velocity. We construct a
lattice model based on a Fokker Planck description of the two-state system,
study it using a kinetic Monte Carlo method and compare our numerical results
with analytical expressions for a mean field limit. We find that the effect of
confinement strongly depends on the bound velocity and the binding kinetics of
the MCC. Confinement effectively reduces the effective diffusivity and average
velocity, except when it results in an enhanced average binding rate and
thereby leads to a larger average velocity than when unconfined.


  Biologically functional liquid-liquid phase separation of intrinsically
disordered proteins (IDPs) is driven by interactions encoded by their amino
acid sequences. Little is currently known about the molecular recognition
mechanisms for distributing different IDP sequences into various cellular
membraneless compartments. Pertinent physics was addressed recently by applying
random-phase-approximation (RPA) polymer theory to electrostatics, which is a
major energetic component governing IDP phase properties. RPA accounts for
charge patterns and thus has advantages over Flory-Huggins and Overbeek-Voorn
mean-field theories. To make progress toward deciphering the phase behaviors of
multiple IDP sequences, the RPA formulation for one IDP species plus solvent is
hereby extended to treat polyampholyte solutions containing two IDP species.
The new formulation generally allows for binary coexistence of two phases, each
containing a different set of volume fractions $(\phi_1,\phi_2)$ for the two
different IDP sequences. The asymmetry between the two predicted coexisting
phases with regard to their $\phi_1/\phi_2$ ratios for the two sequences
increases with increasing mismatch between their charge patterns. This finding
points to a multivalent, stochastic, "fuzzy" mode of molecular recognition that
helps populate various IDP sequences differentially into separate phase
compartments. An intuitive illustration of this trend is provided by
Flory-Huggins models, whereby a hypothetical case of ternary coexistence is
also explored. Augmentations of the present RPA theory with a relative
permittivity $\epsilon_{\rm r}(\phi)$ that depends on IDP volume fraction
$\phi=\phi_1+\phi_2$ lead to higher propensities to phase separate, in line
with the case with one IDP species we studied previously. ...


  We studied intermediate filaments (IFs) in the retina of the Pied flycatcher
(Ficedula hypoleuca) in the foveolar zone. Single IFs span M\"uller cells (MC)
lengthwise; cylindrical bundles of IFs (IFBs) appear inside the cone inner
segment (CIS) at the outer limiting membrane (OLM) level. IFBs adjoin the cone
cytoplasmatic membrane, following lengthwise regularly spaced, forming a
skeleton of the CIS, located above the OLM. IFBs follow along the cone outer
segment (COS), with single IFs separating from the IFB, touching and entering
in-between the light-sensitive disks of the cone membrane. We propose a
mechanism of exciton transfer from the inner retinal surface to the visual
pigments in the photoreceptor cells. This includes excitation transfer in
donor-acceptor systems, from the IF donors to the rhodopsin acceptors, with
theoretic efficiency over 80%. This explains high image contrast in fovea and
foveola in daylight, while the classical mechanism that describes M\"uller
cells as optical lightguides operates in night vision, with loss of resolution
traded for sensitivity. Our theory receives strong confirmation in morphology
and function of the cones and pigment cells. In daylight the lateral surface of
the photosensor disks is blocked from the (scattered or oblique) light by the
pigment cells. Thus the light energy can only get to the cone via intermediate
filaments that absorb photons in the M\"uller cell endfeet and conduct excitons
to the cone. Thus, the disks are consumed at their lateral surfaces, moving to
the apex of the cone, with new disks produced below. An alternative hypothesis
of direct light passing through the cone with its organelles and hitting the
lowest disk contradicts morphological evidence, as thus all of the other disks
would have no useful function in daylight vision.


  A minimal, one-parameter model of the excitonic Hamiltonian for the light
harvesting complex of purple bacteria (LH2) based on ab initio calculation of
the excitonic energies is proposed. The set of input parameters contains
positions of atoms in 27 bacteriochlorophyll pigments only. The excitation
energies, transition dipole and quadrupole moments of the bacteriochlorophyll
units were calculated by advanced multiconfigurational multireference method.
In the work we investigate influence of the dipole-quadrupole interaction on
the particularities of the absorption spectra and circular dichroism and on the
wave-functions localization. We demonstrate that although the spectra changes
are small the dipole-quadrupole interaction term leads, however, to almost
complete separation of the system onto two subsystems: the small B800 ring of
pigments and the large B850 ring. The latter observation is seemed essential
for formulation of energy transport models in light harvesting complexes.


  Sensitivity analysis of biochemical reactions aims at quantifying the
dependence of the reaction dynamics on the reaction rates. The computation of
the parameter sensitivities, however, poses many computational challenges when
taking stochastic noise into account. This paper proposes a new finite
difference method for efficiently computing sensitivities of biochemical
reactions. We employ propensity bounds of reactions to couple the simulation of
the nominal and perturbed processes. The exactness of the simulation is
reserved by applying the rejection-based mechanism. For each simulation step,
the nominal and perturbed processes under our coupling strategy are
synchronized and often jump together, increasing their positive correlation and
hence reducing the variance of the estimator. The distinctive feature of our
approach in comparison with existing coupling approaches is that it only needs
to maintain a single data structure storing propensity bounds of reactions
during the simulation of the nominal and perturbed processes. Our approach
allows to computing sensitivities of many reaction rates simultaneously.
Moreover, the data structure does not require to be updated frequently, hence
improving the computational cost. This feature is especially useful when
applied to large reaction networks. We benchmark our method on biological
reaction models to prove its applicability and efficiency.


  Fluctuations in biodiversity, large and small, are pervasive in the fossil
record, yet we do not understand the processes generating them. Here we extend
theory from non-equilibrium statistical physics to describe the previously
unaccounted for fat-tailed form of fluctuations in marine invertebrate richness
through the Phanerozoic. Using this theory, known as superstatistics, we show
that the simple fact of heterogeneous origination and extinction rates between
clades and conserved rates within clades is sufficient to account for this
fat-tailed form. We identify orders and the families they subsume as the
taxonomic level at which clades experience inter-clade heterogeneity and within
clade homogeneity of rates. Following superstatistics we would posit that
orders and families are subsystems in local statistical equilibrium while the
entire system is not in equilibrium. The separation of timescales between
background origination and extinction within clades compared to the origin of
major ecological and evolutionary innovations leading to new clades allows
within-clade dynamics to reach equilibrium, while between-clade diversification
is non-equilibrial. This between clade non-equilibrium accounts for the
fat-tailed nature of the system as a whole. The distribution of shifts in
diversification dynamics across orders and families is consistent with niche
conservatism and pulsed exploration of adaptive landscapes by higher taxa.
Compared to other approaches that have used simple birth-death processes,
simple equilibrial dynamics, or non-linear theories from complexity science,
superstatistics is superior in its ability to account for both small and
extreme fluctuations in the richness of fossil taxa. Its success opens up new
research directions to better understand the evolutionary processes leading to
stasis in an adaptive landscape interrupted by innovations that lead to novel
forms.


  A model of bit-strings, that uses the technique of multi-spin coding, was
previously used to study the time evolution of B-cell clone repertoire, in a
paper by Lagreca, Almeida and Santos. In this work we extend that simplified
model to include independently the role of the populations of antibodies, in
the control of the immune response, producing mechanisms of differentiation and
regulation in a more complete way. Although the antibodies have the same
molecular shape of the B-cells receptors (BCR), they should present a different
time evolution and thus should be treated separately. We have also studied a
possible model for the network immune memory, suggesting a random memory
regeneration, which is self-perpetuating.


  One important feature of the mammalian immune system is the highly specific
binding of antigens to antibodies. Antibodies generated in response to one
infection may also provide some level of cross immunity to other infections.
One model to describe this cross immunity is the notion of antigenic space,
which assigns each antibody and each virus a point in $\mathbb{R}^n$. Past
studies have suggested the dimensionality of antigenic space, $n$, may be
small. In this study we show that data from hemagglutination assays suggest a
high dimensional random walk (or self avoiding random random walk). The
discrepancy between our result and prior studies is due to the fact that random
walks can appear low dimensional according to a variety of analyses. including
principal component analysis (PCA) and multidimensional scaling (MDS).


  Many simple RNA viruses enclose their genetic material by a protein shell
called the capsid. While the capsid structures are well characterized for most
viruses, the structure of RNA inside the shells and the factors contributing to
it remain poorly understood. We study the impact of base pairing on the
conformations of RNA and find that it undergoes a swollen coil to globule
continuous transition as a function of the strength of the pairing interaction.
We also observe a first order transition and kink profile as a function of RNA
length. All these transitions could explain the different RNA profiles observed
inside viral shells.


  Chromosome conformation capture experiments such as Hi-C are used to map the
three-dimensional spatial organization of genomes. One specific feature of the
3D organization is known as topologically associating domains (TADs), which are
densely interacting, contiguous chromatin regions playing important roles in
regulating gene expression. A few algorithms have been proposed to detect TADs.
In particular, the structure of Hi-C data naturally inspires application of
community detection methods. However, one of the drawbacks of community
detection is that most methods take exchangeability of the nodes in the network
for granted; whereas the nodes in this case, i.e. the positions on the
chromosomes, are not exchangeable. We propose a network model for detecting
TADs using Hi-C data that takes into account this non-exchangeability. In
addition, our model explicitly makes use of cell-type specific CTCF binding
sites as biological covariates and can be used to identify conserved TADs
across multiple cell types. The model leads to a likelihood objective that can
be efficiently optimized via relaxation. We also prove that when suitably
initialized, this model finds the underlying TAD structure with high
probability. Using simulated data, we show the advantages of our method and the
caveats of popular community detection methods, such as spectral clustering, in
this application. Applying our method to real Hi-C data, we demonstrate the
domains identified have desirable epigenetic features and compare them across
different cell types.


  There does not exist a general positive correlation between important
life-supporting properties and the entropy production rate. The simple reason
is that nondissipative and time-symmetric kinetic aspects are also relevant for
establishing optimal functioning. In fact those aspects are even crucial in the
nonlinear regimes around equilibrium where we find biological processing on
mesoscopic scales. We make these claims specific via examples of molecular
motors, of circadian cycles and of sensory adaptation, whose performance in
some regimes is indeed spoiled by increasing the dissipated power. We use the
relation between dissipation and the amount of time-reversal breaking to keep
the discussion quantitative also in effective models where the physical entropy
production is not clearly identifiable.


  In this chapter we review concepts and theories of polymer dynamics. We think
of it as an introduction to the topic for scientists specializing in other
subfields of statistical mechanics and condensed matter theory, so, for the
readers reference, we start with a short review of the equilibrium static
properties of polymer systems. Most attention is paid to the dynamics of
unentangled polymer systems, where apart from classical Rouse and Zimm models
we review some recent scaling and analytical generalizations. The dynamics of
systems with entanglements is also briefly reviewed. Special attention is paid
to the discussion of comparatively weakly understood topological states of
polymer systems and possible approaches to the description of their dynamics.


  The Shapley value, a solution concept from cooperative game theory, has
recently been considered for both unrooted and rooted phylogenetic trees. Here,
we focus on the Shapley value of unrooted trees and first revisit the so-called
split counts of a phylogenetic tree and the Shapley transformation matrix that
allows for the calculation of the Shapley value from the edge lengths of a
tree. We show that non-isomorphic trees may have permutation-equivalent Shapley
transformation matrices and permutation-equivalent null spaces. This implies
that estimating the split counts associated with a tree or the Shapley values
of its leaves does not suffice to reconstruct the correct tree topology. We
then turn to the use of the Shapley value as a prioritization criterion in
biodiversity conservation and compare it to a greedy solution concept. Here, we
show that for certain phylogenetic trees, the Shapley value may fail as a
prioritization criterion, meaning that the diversity spanned by the top $k$
species (ranked by their Shapley values) cannot approximate the total diversity
of all $n$ species.


  We describe the use of the Frechet mean and variance in the
Billera-Holmes-Vogtmann (BHV) treespace to summarize and explore the diversity
of a set of phylogenetic trees. We show that the Frechet mean is comparable to
other summary methods, and, despite its stickiness property, is more likely to
be binary than the majority-rules consensus tree. We show that the Frechet
variance is faster and more precise than commonly used variance measures. The
Frechet mean and variance are more theoretically justified, and more robust,
than previous estimates of this type, and can be estimated reasonably
efficiently, providing a foundation for building more advanced statistical
methods and leading to applications such as mean hypothesis testing.


  Knowledge of regional net primary productivity (NPP) is important for the
systematic understanding of the global carbon cycle. In this study,
multi-source data were employed to conduct a 33-year regional NPP study in
southwest China, at a 1-km scale. A multi-sensor fusion framework was applied
to obtain a new normalized difference vegetation index (NDVI) time series from
1982 to 2014, combining the respective advantages of the different remote
sensing datasets. As another key parameter for NPP modeling, the total solar
radiation was calculated by the improved Yang hybrid model (YHM), using
meteorological station data. The verification described in this paper proved
the feasibility of all the applied data processes, and a greatly improved
accuracy was obtained for the NPP calculated with the final processed NDVI. The
spatio-temporal analysis results indicated that 68.07% of the study area showed
an increasing NPP trend over the past three decades. Significant heterogeneity
was found in the correlation between NPP and precipitation at a monthly scale,
specifically, the negative correlation in the growing season and the positive
correlation in the dry season. The lagged positive correlation in the growing
season and no lag in the dry season indicated the important impact of
precipitation on NPP.


  Recent developments in automated tracking allow uninterrupted,
high-resolution recording of animal trajectories, sometimes coupled with the
identification of stereotyped changes of body pose or other behaviors of
interest. Analysis and interpretation of such data represents a challenge: the
timing of animal behaviors may be stochastic and modulated by kinematic
variables, by the interaction with the environment or with the conspecifics
within the animal group, and dependent on internal cognitive or behavioral
state of the individual. Existing models for collective motion typically fail
to incorporate the discrete, stochastic, and internal-state-dependent aspects
of behavior, while models focusing on individual animal behavior typically
ignore the spatial aspects of the problem. Here we propose a probabilistic
modeling framework to address this gap. Each animal can switch stochastically
between different behavioral states, with each state resulting in a possibly
different law of motion through space. Switching rates for behavioral
transitions can depend in a very general way, which we seek to identify from
data, on the effects of the environment as well as the interaction between the
animals. We represent the switching dynamics as a Generalized Linear Model and
show that: (i) forward simulation of multiple interacting animals is possible
using a variant of the Gillespie's Stochastic Simulation Algorithm; (ii)
formulated properly, the maximum likelihood inference of switching rate
functions is tractably solvable by gradient descent; (iii) model selection can
be used to identify factors that modulate behavioral state switching and to
appropriately adjust model complexity to data. To illustrate our framework, we
apply it to two synthetic models of animal motion and to real zebrafish
tracking data.


  Since the pioneering works of Pethig, Grant and Wuthrich on protein hydration
layer, many studies have been devoted to find out if there are any general and
universal characteristic features that can distinguish water molecules inside
the protein hydration layer from bulk. Given that the surface itself varies
from protein to protein, and that each surface facing the water is
heterogeneous, search for universal features has been elusive. Here, we perform
atomistic molecular dynamics simulation in order to propose and demonstrate
that such defining characteristics can emerge if we look not at average
properties but the distribution of relaxation times. We present results of
calculations of distributions of residence times and rotational relaxation
times for four different protein-water systems, and compare them with the same
quantities in the bulk. The distributions in the hydration layer is unusually
broad and log-normal in nature, due to the simultaneous presence of peptide
backbones that form weak hydrogen bonds, hydrophobic amino acid side chains
that form no hydrogen bond and charged polar groups that form strong hydrogen
bond with the surrounding water molecules. The broad distribution is
responsible for the non-exponential dielectric response and also agrees with
large specific heat of the hydration water. Our calculations reveal that while
the average time constant is just about 2-3 times larger than that of bulk
water, it provides a poor representation of the real behaviour. In particular,
the average leads to the erroneous conclusion that water in the hydration layer
is bulk-like. However, the observed and calculated lower value of static
dielectric constant of hydration layer remained difficult to reconcile with the
broad distribution observed in dynamical properties. We offer a plausible
explanation of these unique properties.


  In this work, we propose a novel centrality metric, referred to as star
centrality, which incorporates information from the closed neighborhood of a
node, rather than solely from the node itself, when calculating its topological
importance. More specifically, we focus on degree centrality and show that in
the complex protein-protein interaction networks it is a naive metric that can
lead to misclassifying protein importance. For our extension of degree
centrality when considering stars, we derive its computational complexity,
provide a mathematical formulation, and propose two approximation algorithms
that are shown to be efficient in practice. We portray the success of this new
metric in protein-protein interaction networks when predicting protein
essentiality in several organisms, including the well-studied Saccharomyces
cerevisiae, Helicobacter pylori, and Caenorhabditis elegans, where star
centrality is shown to significantly outperform other nodal centrality metrics
at detecting essential proteins. We also analyze the average and worst case
performance of the two approximation algorithms in practice, and show that they
are viable options for computing star centrality in very large-scale
protein-protein interaction networks, such as the human proteome, where exact
methodologies are bound to be time and memory intensive.


  We study the phenomenon of cyclic dominance in the paradigmatic
Rock--Paper--Scissors model, as occurring in both stochastic individual-based
models of finite populations and in the deterministic replicator equations. The
mean-field replicator equations are valid in the limit of large populations
and, in the presence of mutation and unbalanced payoffs, they exhibit an
attracting limit cycle. The period of this cycle depends on the rate of
mutation; specifically, the period grows logarithmically as the mutation rate
tends to zero. We find that this behaviour is not reproduced in stochastic
simulations with a fixed finite population size. Instead, demographic noise
present in the individual-based model dramatically slows down the progress of
the limit cycle, with the typical period growing as the reciprocal of the
mutation rate. Here we develop a theory that explains these scaling regimes and
delineates them in terms of population size and mutation rate. We identify a
further intermediate regime in which we construct a stochastic differential
equation model describing the transition between stochastically-dominated and
mean-field behaviour.


  The structure of ecological interactions is commonly understood through
analyses of interaction networks. However, these analyses may be sensitive to
sampling biases in both the interactors (the nodes of the network) and
interactions (the links between nodes). These issues may affect the accuracy of
empirically constructed ecological networks. We explore the properties of
sampled ecological networks by simulating large-scale ecological networks with
predetermined topologies, and sampling them with different mathematical
procedures. Several types of modular networks were generated, intended to
represent a wide variety of communities that vary in size and types of
ecological interactions. We sampled these networks with different sampling
designs that may be encountered in field experiments. The observed networks
generated by each sampling process were analyzed with respect to number and
size of components. We show that the sampling effort needed to estimate
underlying network properties depends both on the sampling design and on
network topology. Networks with random or scale-free modules require more
complete sampling compared to networks whose modules are nested or bipartite.
Overall, the structure of nested modules was the easiest to detect, regardless
of sampling design. Sampling according to species degree was consistently found
to be the most accurate strategy to estimate network structure. Conversely,
sampling according to module results in an accurate view of certain modules,
but fails to provide a global picture of the underlying network. We recommend
that these findings are incorporated into the design of projects aiming to
characterize large networks of species interactions in the field, to reduce
sampling biases. The software scripts developed to construct and sample
networks are provided for further explorations of network structure and
comparisons to real interaction data.


  In this paper we propose some mathematical models for the transmission of
dengue using a system of reaction-diffusion equations. The mosquitoes are
divided into infected, uninfected and aquatic subpopulations, while the humans,
which are divided into susceptible, infected and recovered, are considered
homogeneously distributed in space and with a constant total population. We
find Lie point symmetries of the models and we study theirs temporal dynamics,
which provides us the regions of stability and instability, depending on the
values of the basic offspring and the basic reproduction numbers. Also, we
calculate the possible values of the wave speed for the mosquitoes invasion and
dengue spread and compare them with those found in the literature.


  In this paper we study some deterministic mathematical models that seek to
explain the expansion of zika virus, as a viral epidemic, using published data
for Brazil. SIR type models are proposed and validated using the epidemic data
found, considering several aspects in the spread of the disease. Finally, we
confirmed that the crucial epidemic parameter such as $R_0$ is consistent with
those previously reported in the literature for other areas. We also explored
variations of the parameters within Brazil for different federal entities. We
concluded that a parsimonious model that includes both human and vector
populations best describe the epidemic parameters.


  Motivation: Recent advances in sequencing technologies promise ultra-long
reads of $\sim$100 kilo bases (kb) in average, full-length mRNA or cDNA reads
in high throughput and genomic contigs over 100 mega bases (Mb) in length.
Existing alignment programs are unable or inefficient to process such data at
scale, which presses for the development of new alignment algorithms.
  Results: Minimap2 is a general-purpose alignment program to map DNA or long
mRNA sequences against a large reference database. It works with accurate short
reads of $\ge$100bp in length, $\ge$1kb genomic reads at error rate $\sim$15%,
full-length noisy Direct RNA or cDNA reads, and assembly contigs or closely
related full chromosomes of hundreds of megabases in length. Minimap2 does
split-read alignment, employs concave gap cost for long insertions and
deletions (INDELs) and introduces new heuristics to reduce spurious alignments.
It is 3-4 times faster than mainstream short-read mappers at comparable
accuracy and $\ge$30 times faster at higher accuracy for both genomic and mRNA
reads, surpassing most aligners specialized in one type of alignment.
  Availability and implementation: https://github.com/lh3/minimap2
  Contact: hengli@broadinstitute.org


  Drug repositioning, that is finding new uses for existing drugs to treat more
patients. Cumulative studies demonstrate that the mature miRNAs as well as
their precursors can be targeted by small molecular drugs. At the same time,
human diseases result from the disordered interplay of tissue- and cell
lineage-specific processes. However, few computational researches predict
drug-disease potential relationships based on miRNA data and tissue
specificity. Therefore, based on miRNA data and the tissue specificity of
diseases, we propose a new method named as miTS to predict the potential
treatments for diseases. Firstly, based on miRNAs data, target genes and
information of FDA approved drugs, we evaluate the relationships between miRNAs
and drugs in the tissue-specific PPI network. Then, we construct a tripartite
network: drug-miRNA-disease Finally, we obtain the potential drug-disease
associations based on the tripartite network. In this paper, we take breast
cancer as case study and focus on the top-30 predicted drugs. 25 of them
(83.3%) are found having known connections with breast cancer in CTD benchmark
and the other 5 drugs are potential drugs for breast cancer. We further
evaluate the 5 newly predicted drugs from clinical records, literature mining,
KEGG pathways enrichment analysis and overlapping genes between enriched
pathways. For each of the 5 new drugs, strongly supported evidences can be
found in three or more aspects. In particular, Regorafenib has 15 overlapping
KEGG pathways with breast cancer and their p-values are all very small. In
addition, whether in the literature curation or clinical validation,
Regorafenib has a strong correlation with breast cancer. All the facts show
that Regorafenib is likely to be a truly effective drug, worthy of our further
study. It further follows that our method miTS is effective and practical for
predicting new drug indications.


  In this paper we study a class of stochastic individual-based models that
describe the evolution of haploid populations where each individual is
characterised by a phenotype and a genotype. The phenotype of an individual
determines its natural birth- and death rates as well as the competition
kernel, $c(x,y)$ which describes the induced death rate that an individual of
type $x$ experiences due to the presence of an individual or type $y$. When a
new individual is born, with a small probability a mutation occurs, i.e. the
offspring has different genotype as the parent. The novel aspect of the models
we study is that an individual with a given genotype may express a certain set
of different phenotypes, and during its lifetime it may switch between
different phenotypes, with rates that are much larger then the mutation rates
and that, moreover, may depend on the state of the entire population. The
evolution of the population is described by a continuous-time, measure-valued
Markov process. In our last paper [4], such a model was proposed to describe
tumor evolution under immunotherapy. In the present paper we consider a large
class of models which comprises the example studied in [4] and analyse their
scaling limits as the population size tends to infinity and the mutation rate
tends to zero. Under suitable assumptions, we prove convergence to a Markov
jump process that is a generalisation of the polymorphic evolution sequence
(PES) as analysed by Champagnat and M\'el\'eard [8, 10].


  We present a model of the dynamical transition of atomic displacements in
proteins. Increased mean-square displacement at higher temperatures is caused
by softening of the vibrational force constant by electrostatic and van der
Waals forces from the protein-water thermal bath. Vibrational softening passes
through a nonergodic dynamical transition when the relaxation time of the
force-force correlation function enters, with increasing temperature, the
instrumental observation window. Two crossover temperatures are identified. The
lower crossover, presently connected to the glass transition, is related to the
dynamical unfreezing of rotations of water molecules within nanodomains
polarized by charged surface residues of the protein. The higher crossover
temperature, usually assigned to the dynamical transition, marks the onset of
water translations. All crossovers are ergodicity breaking transitions
depending on the corresponding observation windows. Allowing stretched
exponential relaxation of the protein-water thermal bath significantly improves
the theory-experiment agreement when applied to solid protein samples studied
by Moessbauer spectroscopy.


  The analysis of equilibrium points in random games has been of great interest
in evolutionary game theory, with important implications for understanding of
complexity in a dynamical system, such as its behavioural, cultural or
biological diversity. The analysis so far has focused on random games of
independent payoff entries. In this paper, we overcome this restrictive
assumption by considering multi-player two-strategy evolutionary games where
the payoff matrix entries are correlated random variables. Using techniques
from the random polynomial theory we establish a closed formula for the mean
numbers of internal (stable) equilibria. We then characterise the asymptotic
behaviour of this important quantity for large group sizes and study the effect
of the correlation. Our results show that decreasing the correlation among
payoffs (namely, of a strategist for different group compositions) leads to
larger mean numbers of (stable) equilibrium points, suggesting that the system
or population behavioural diversity can be promoted by increasing independence
of the payoff entries. Numerical results are provided to support the obtained
analytical results.


  The exchange of diffusive metabolites is known to control the spatial
patterns formed by microbial populations, as revealed by recent studies in the
laboratory. However, the matrices used, such as agarose pads, lack the
structured geometry of many natural microbial habitats, including in the soil
or on the surfaces of plants or animals. Here we address the important question
of how such geometry may control diffusive exchanges and microbial interaction.
We model mathematically mutualistic interactions within a minimal unit of
structure: two growing reservoirs linked by a diffusive channel through which
metabolites are exchanged. The model is applied to study a synthetic mutualism,
experimentally parameterised on a model algal-bacterial co-culture. Analytical
and numerical solutions of the model predict conditions for the successful
establishment of remote mutualisms, and how this depends, often
counterintutively, on diffusion geometry. We connect our findings to
understanding complex behaviour in synthetic and naturally occurring microbial
communities.


  We introduce and study a set of training-free methods of
information-theoretic and algorithmic complexity nature applied to DNA
sequences to identify their potential capabilities to determine nucleosomal
binding sites. We test our measures on well-studied genomic sequences of
different sizes drawn from different sources. The measures reveal the known in
vivo versus in vitro predictive discrepancies and uncover their potential to
pinpoint (high) nucleosome occupancy. We explore different possible signals
within and beyond the nucleosome length and find that complexity indices are
informative of nucleosome occupancy. We compare against the gold standard
(Kaplan model) and find similar and complementary results with the main
difference that our sequence complexity approach. For example, for high
occupancy, complexity-based scores outperform the Kaplan model for predicting
binding representing a significant advancement in predicting the highest
nucleosome occupancy following a training-free approach.


  Many proteoforms - arising from alternative splicing, post-translational
modifications (PTMs), or paralogous genes - have distinct biological functions,
such as histone PTM proteoforms. However, their quantification by existing
bottom-up mass-spectrometry (MS) methods is undermined by peptide-specific
biases. To avoid these biases, we developed and implemented a first-principles
model (HIquant) for quantifying proteoform stoichiometries. We characterized
when MS data allow inferring proteoform stoichiometries by HIquant, derived an
algorithm for optimal inference, and demonstrated experimentally high accuracy
in quantifying fractional PTM occupancy without using external standards, even
in the challenging case of the histone modification code.
  HIquant server is implemented at:
https://web.northeastern.edu/slavov/2014_HIquant/


  Retroactivity arises when the coupling of a molecular network $\mathcal{U}$
to a downstream network $\mathcal{D}$ results in signal propagation back from
$\mathcal{D}$ to $\mathcal{U}$. The phenomenon represents a breakdown in
modularity of biochemical circuits and hampers the rational design of complex
functional networks. Considering simple models of signal-transduction
architectures, we demonstrate the strong dependence of retroactivity on the
properties of the upstream system, and explore the cost and efficacy of
fuel-consuming insulating motifs that can mitigate retroactive effects. We find
that simple insulating motifs can suppress retroactivity at a low fuel cost by
coupling only weakly to the upstream system $\mathcal{U}$. However, this design
approach reduces the signalling network's robustness to perturbations from leak
reactions, and potentially compromises its ability to respond to
rapidly-varying signals.


  Stochastic partial differential equations (SPDE) on graphs were introduced by
Cerrai and Freidlin [Ann. Inst. Henri Poincar\'e Probab. Stat. 53 (2017)
865-899]. This class of stochastic equations in infinite dimensions provides a
minimal framework for the study of the effective dynamics of much more complex
systems. However, how they emerge from microscopic individual-based models is
still poorly understood, partly due to complications near vertex singularities.
In this work, motivated by the study of the dynamics and the genealogies of
expanding populations in spatially structured environments, we obtain a new
class of SPDE on graphs of Wright-Fisher type which have nontrivial boundary
conditions on the vertex set. We show that these SPDE arise as scaling limits
of suitably defined biased voter models (BVM), which extends the scaling limits
of Durrett and Fan [Ann. Appl. Probab. 26 (2016) 456-490]. We further obtain a
convergent simulation scheme for each of these SPDE in terms of a system of
It\^o SDEs, which is useful when the size of the BVM is too large for
stochastic simulations. These give the first rigorous connection between SPDE
on graphs and more discrete models, specifically, interacting particle systems
and interacting SDEs. Uniform heat kernel estimates for symmetric random walks
approximating diffusions on graphs are the keys to our proofs.


  Protein-protein interactions (PPIs) perform important roles on biological
functions. Researches of mutants on protein interactions can further understand
PPIs. In the past, many researchers have developed databases that stored
mutants on protein interactions, which are old and not updated till now. To
address the issue, we developed a kinetic and thermodynamic database of mutant
protein interactions (dbMPIKT) that can be freely accessible at our website.
This database contains 5291 mutants that integrated data from previous
databases and data from literatures for nearly three years. Furthermore, the
data were analyzed, involving mutation number, mutation type, protein pair
source and network map construction. On the whole, the database provides new
data to further improve the study on PPIs. Website:
http://210.45.212.128/lqy/index.php


  We show that higher order epistasis has a substantial impact on evolutionary
dynamics by analyzing peaks in the fitness landscapes. There are 193,270,310
fitness graphs, or cube orientations, for 4-locus systems, distributed on 511,
863 isomorphism classes. We identify all fitness graphs with 6 or more peaks.
81 percent of them imply 4-way epistasis, whereas 9 percent of all 4-locus
fitness graphs imply 4-way epistasis. Fitness graphs are useful in that they
reflect the entire collection of fitness landscapes rather than focusing on a
particular model. Our results depend on a characterization of fitness graphs
that imply $n$-way epistasis. The characterization is expressed in terms of a
partition property that can be derived from Hall's marriage theorem for
bipartite graphs. A similar partition condition holds for any partial order.
The result answers an open problem posed at a conference on interactions
between algebra and the sciences at the Max Planck institute.


  Under supercoiling constraints, naked DNA, such as a large part of bacterial
DNA, folds into braided structures called plectonemes. The double-helix can
also undergo local structural transitions, leading to the formation of
denaturation bubbles and other alternative structures. Various polymer models
have been developed to capture these properties, with Monte-Carlo (MC)
approaches dedicated to the inference of thermodynamic properties. In this
chapter, we explain how to perform such Monte-Carlo simulations, following two
objectives. On one hand, we present the self-avoiding supercoiled Worm-Like
Chain (ssWLC) model, which is known to capture the folding properties of
supercoiled DNA, and provide a detailed explanation of a standard MC simulation
method. On the other hand, we explain how to extend this ssWLC model to include
structural transitions of the helix.


  The generalized distance matrix of a graph is the matrix whose entries depend
only on the pairwise distances between vertices, and the generalized distance
spectrum is the set of eigenvalues of this matrix. This framework generalizes
many of the commonly studied spectra of graphs. We show that for a large class
of graphs these eigenvalues can be computed explicitly. We also present the
applications of our results to competition models in ecology and rapidly mixing
Markov Chains.


  Protein-protein interaction (PPI) networks are the backbone of all processes
in living cells. In this work we relate conservation, essentiality and
functional repertoire of a gene to the connectivity $k$ (i.e., the number of
interaction links) of the corresponding protein in the PPI network. On a set of
42 bacterial genomes of different sizes, and with reasonably separated
evolutionary trajectories, we investigate three issues: i) whether the
distribution of connectivities changes between PPI subnetworks of essential and
nonessential genes; ii) how gene conservation, measured both by the
evolutionary retention index (ERI) and by evolutionary pressures, is related to
the the connectivity of the corresponding protein; iii) how PPI connectivities
are modulated by evolutionary and functional relationships, as represented by
the Clusters of Orthologous Genes (COGs). We show that conservation,
essentiality and functional specialisation of genes constrain the connectivity
of the corresponding proteins in bacterial PPI networks. In particular, we
isolate a core of highly connected proteins (with connectivities $k\ge40$),
which is ubiquitous among the species considered here -- though mostly visible
in the degree distributions of bacteria with small genomes (less than 1000
genes). The genes that belong to this highly connected core are conserved,
essential and, in most cases, belong to the COG cluster J, related to ribosomal
functions and to the processing of genetic information.


  We consider the problem of inferring an ancestral state from observations at
the leaves of a tree, assuming the state evolves along the tree according to a
two-state symmetric Markov process. We establish a general branching rate
condition under which maximum parsimony, a common reconstruction method
requiring only the knowledge of the tree, succeeds better than random guessing
uniformly in the depth of the tree. We thereby generalize previous results of
(Zhang et al., 2010) and (Gascuel and Steel, 2010). Our results apply to both
deterministic and i.i.d. edge weights.


  Background: While RSV circulation results in high burden of hospitalization,
particularly among infants, young children and the elderly, little is known
about the role of different age groups in propagating annual RSV epidemics in
the community.
  Methods: During a communicable disease outbreak, some subpopulations may play
a disproportionate role during the outbreak's ascent due to increased
susceptibility and/or contact rates. Such subpopulations can be identified by
considering the proportion that cases in a subpopulation represent among all
cases in the population occurring before (Bp) and after the epidemic peak (Ap)
to calculate the subpopulation's relative risk, RR=Bp/Ap. We estimated RR for
several age groups using data on RSV hospitalizations in the US between
2001-2012 from the Healthcare Cost and Utilization Project (HCUP).
  Results: Children aged 3-4y and 5-6y each had the highest RR estimate for
5/11 seasons in the data, with RSV hospitalization rates in infants being
generally higher during seasons when children aged 5-6y had the highest RR
estimates. Children aged 2y had the highest RR estimate during one season. RR
estimates in infants and individuals aged 11y and older were mostly lower than
in children aged 1-10y.
  Conclusions: The RR estimates suggest that preschool and young school-age
children have the leading relative roles during RSV epidemics. We hope that
those results will aid in the design of RSV vaccination policies.


  The phylogenetic tree space introduced by Billera, Holmes, and Vogtmann (BHV
tree space) is a CAT(0) continuous space that represents trees with edge
weights with an intrinsic geodesic distance measure. The geodesic distance
measure unique to BHV tree space is well known to be computable in polynomial
time, which makes it a potentially powerful tool for optimization problems in
phylogenetics and phylogenomics. Specifically, there is significant interest in
comparing and combining phylogenetic trees. For example, BHV tree space has
been shown to be potentially useful in tree summary and consensus methods,
which require combining trees with different number of leaves. Yet an open
problem is to transition between BHV tree spaces of different maximal
dimension, where each maximal dimension corresponds to the complete set of
edge-weighted trees with a fixed number of leaves. We show a combinatorial
method to transition between copies of BHV tree spaces in which trees with
different numbers of taxa can be studied, derived from its topological
structure and geometric properties. This method removes obstacles for embedding
problems such as supertree and consensus methods in the BHV treespace
framework.


  Due to the asymmetric nature of the nucleotides, the extant informational
biomolecule, DNA, is constrained to replicate unidirectionally on a template.
As a product of molecular evolution that sought to maximize replicative
potential, DNA's unidirectional replication poses a mystery since symmetric
bidirectional self-replicators obviously would replicate faster than
unidirectional self-replicators and hence would have been evolutionarily more
successful. Here we carefully examine the physico-chemical requirements for
evolutionarily successful primordial self-replicators and theoretically show
that at low monomer concentrations that possibly prevailed in the primordial
oceans, asymmetric unidirectional self-replicators would have an evolutionary
advantage over bidirectional self-replicators. The competing requirements of
low and high kinetic barriers for formation and long lifetime of inter-strand
bonds respectively are simultaneously satisfied through asymmetric kinetic
influence of inter-strand bonds, resulting in evolutionarily successful
unidirectional self-replicators.


  Starting from the Fisher-Kolmogorov-Petrovskii-Piskunov equation (FKPP) we
model the dynamic of a diffusive system with two mutually communicating
identical patches and isolated of the remaining matrix. For this system we find
the minimal size of each fragment in the explicit form and compare with the
explicit results for similar problems found in the literature. From this
comparison emerges an unexpected result that for a same set of the parameters,
the isolated system studied in this work with size L, can be better or worst
than the non isolated systems with the same size L, uniquely depending on the
parameter $a_{0}$ (internal conditions of the patches). Due to the fact that
this result is unexpected we propose a experimental verification.


  Despite high vaccine coverage, pertussis has re-emerged as a public health
concern in many countries. One hypothesis posed for re-emergence is the waning
of immunity. In some disease systems, the process of waning immunity can be
non-linear, involving a complex relationship between the duration of immunity
and subsequent boosting of immunity through asymptomatic re-exposure.
  We present and analyse a model of infectious disease transmission to examine
the interplay between infection and immunity. By allowing the duration of
infection-acquired immunity to differ from that of vaccine-acquired immunity,
we explore the impact of the difference in durations on long-term disease
patterns and prevalence of infection.
  Our model demonstrates that vaccination may induce cyclic behaviour, and its
ability to reduce the infection prevalence increases with both the duration of
infection-acquired immunity and duration of vaccine-acquired immunity. We find
that increasing vaccine coverage, while capable of leading to an increase in
overall transmission, always results in a reduction in prevalence of primary
infections, with epidemic cycles characterised by a longer interepidemic period
and taller peaks.
  Our results show that the epidemiological patterns of an infectious disease
may change considerably when the duration of vaccine-acquired immunity differs
from that of infection-acquired immunity. Our study highlights that for any
particular disease and associated vaccine, a detailed understanding of the
duration of protection and how that duration is influenced by infection
prevalence is important as we seek to optimise vaccination strategies.


  HIV/AIDS spread depends upon complex patterns of interaction among various
sub-sets emerging at population level. This added complexity makes it difficult
to study and model AIDS and its dynamics. AIDS is therefore a natural candidate
to be modeled using agent-based modeling, a paradigm well-known for modeling
Complex Adaptive Systems (CAS). While agent-based models are also well-known to
effectively model CAS, often times models can tend to be ambiguous and the use
of purely text-based specifications (such as ODD) can make models difficult to
be replicated. Previous work has shown how formal specification may be used in
conjunction with agent-based modeling to develop models of various CAS.
However, to the best of our knowledge, no such model has been developed in
conjunction with AIDS. In this paper, we present a Formal Agent-Based
Simulation modeling framework (FABS-AIDS) for an AIDS-based CAS. FABS-AIDS
employs the use of a formal specification model in conjunction with an
agent-based model to reduce ambiguity as well as improve clarity in the model
definition. The proposed model demonstrates the effectiveness of using formal
specification in conjunction with agent-based simulation for developing models
of CAS in general and, social network-based agent-based models, in particular.


  Markov state models (MSMs) have been widely used to analyze computer
simulations of various biomolecular systems. They can capture conformational
transitions much slower than an average or maximal length of a single molecular
dynamics (MD) trajectory from the set of trajectories used to build the MSM. A
rule of thumb claiming that the slowest implicit timescale captured by an MSM
should be comparable by the order of magnitude to the aggregate duration of all
MD trajectories used to build this MSM has been known in the field. However,
this rule have never been formally proved. In this work, we present analytical
results for the slowest timescale in several types of MSMs, supporting the
above rule. We conclude that the slowest implicit timescale equals the product
of the aggregate sampling and four factors that quantify: (1) how much
statistics on the conformational transitions corresponding to the longest
implicit timescale is available, (2) how good the sampling of the destination
Markov state is, (3) the gain in statistics from using a sliding window for
counting transitions between Markov states, and (4) a bias in the estimate of
the implicit timescale arising from finite sampling of the conformational
transitions. We demonstrate that in many practically important cases all these
four factors are on the order of unity, and we analyze possible scenarios that
could lead to their significant deviation from unity. Overall, we provide for
the first time analytical results on the slowest timescales captured by MSMs.
These results can guide further practical applications of MSMs to biomolecular
dynamics and allow for higher computational efficiency of simulations.


  We develop and analyze an ordinary differential equation model to assess the
potential effectiveness of infecting mosquitoes with the Wolbachia bacteria to
control the ongoing mosquito-borne epidemics, such as dengue fever,
chikungunya, and Zika. Wolbachia is a natural parasitic microbe that stops the
proliferation of the harmful viruses inside the mosquito and reduces disease
transmission. It is difficult to sustain an infection of the maternal
transmitted Wolbachia in a wild mosquito population because of the reduced
fitness of the Wolbachia-infected mosquitoes and cytoplasmic incompatibility
limiting maternal transmission. The infection will only persist if the fraction
of the infected mosquitoes exceeds a minimum threshold. Our two-sex mosquito
model captures the complex transmission-cycle by accounting for heterosexual
transmission, multiple pregnant states for female mosquitoes, and the
aquatic-life stage. We identify important dimensionless numbers and analyze the
critical threshold condition for obtaining a sustained Wolbachia infection in
the natural population. This threshold effect is characterized by a backward
bifurcation with three coexisting equilibria of the system of differential
equations: a stable disease-free equilibrium, an unstable
intermediate-infection endemic equilibrium and a stable high-infection endemic
equilibrium. We perform sensitivity analysis on epidemiological and
environmental parameters to determine their relative importance to Wolbachia
transmission and prevalence. We also compare the effectiveness of different
integrated mitigation strategies and observe that the most efficient approach
to establish the Wolbachia infection is to first reduce the natural mosquitoes
and then release both infected males and pregnant females. The initial
reduction of natural population could be accomplished by either residual
spraying or ovitraps.


  The main result of this article is the determination of the minimal size for
the general case of problems with two identical patches. This solution is
presented in the explicit form, which allows to recuperate all the cases found
in the literature as particular cases, namely, one isolated fragment, one
single fragment communicating with its neighborhood, a system with two
identical fragments isolated from the matrix but mutually communicating and a
system of two identical fragments inserted in a homogeneous matrix. It is also
addressed the new problem of a single fragment communicating with the matrix,
with different life difficulty of each side. As application, it is found that
the internal condition $a_{0}$ can set which system is the worst to life. This
prediction confirms and extends the prediction already found in the literature
between isolated and non-isolated systems.


  Empirical observations show that ecological communities can have a huge
number of coexisting species, also with few or limited number of resources.
These ecosystems are characterized by multiple type of interactions, in
particular displaying cooperative behaviors. However, standard modeling of
population dynamics based on Lotka-Volterra type of equations predicts that
ecosystem stability should decrease as the number of species in the community
increases and that cooperative systems are less stable than communities with
only competitive and/or exploitative interactions. Here we propose a stochastic
model of population dynamics, which includes exploitative interactions as well
as cooperative interactions induced by cross-feeding. The model is exactly
solved and we obtain results for relevant macro-ecological patterns, such as
species abundance distributions and correlation functions. In the large system
size limit, any number of species can coexist for a very general class of
interaction networks and stability increases as the number of species grows.
For pure mutualistic/commensalistic interactions we determine the topological
properties of the network that guarantee species coexistence. We also show that
the stationary state is globally stable and that inferring species interactions
through species abundance correlation analysis may be misleading. Our
theoretical approach thus show that appropriate models of cooperation naturally
leads to a solution of the long-standing question about complexity-stability
paradox and on how highly biodiverse communities can coexist.


  Recent technological changes have increased connectivity between individuals
around the world leading to higher frequency interactions between members of
communities that would be otherwise distant and disconnected. This paper
examines a model of opinion dynamics in interacting communities and studies how
increasing interaction frequency affects the ability for communities to retain
distinct identities versus falling into consensus or polarized states in which
community identity is lost. We also study the effect (if any) of opinion noise
related to a tendency for individuals to assert their individuality in
homogenous populations. Our work builds on a model we developed previously [11]
where the dynamics of opinion change is based on individual interactions that
seek to minimize some energy potential based on the differences between
opinions across the population.


  Consanguinity or inter-cousin marriage is a phenomenon quite prevalent in
certain regions around the globe. Consanguineous parents have a higher risk of
having offspring with congenital disorders. It is difficult to model large
scale consanguineous parental populations because of disparate cultural issues
unique to regions and cultures across the globe. Although consanguinity, as a
social problem has been studied previously, consanguinity from a biological
perspective has yet to be modeled. Discrete Event System Specification (DEVS)
formalism is a powerful modeling formalism for the study of intricate details
of real-world complex systems. In this article, we develop a DEVS model to get
an insight into the role of consanguineous marriages in the evolution of
congenital disorders in a population. As proof-of-concept, we develop a
consanguinity simulation model in Simio simulation software. Our results show
the effectiveness of DEVS in the modeling of consanguinity effects in causing
congenital defects.


  Motivation: The rapid growth of diverse biological data allows us to consider
interactions between a variety of objects, such as genes, chemicals, molecular
signatures, diseases, pathways and environmental exposures. Often, any pair of
objects--such as a gene and a disease--can be related in different ways, for
example, directly via gene-disease associations or indirectly via functional
annotations, chemicals and pathways. Different ways of relating these objects
carry different semantic meanings. However, traditional methods disregard these
semantics and thus cannot fully exploit their value in data modeling.
  Results: We present Medusa, an approach to detect size-k modules of objects
that, taken together, appear most significant to another set of objects. Medusa
operates on large-scale collections of heterogeneous data sets and explicitly
distinguishes between diverse data semantics. It advances research along two
dimensions: it builds on collective matrix factorization to derive different
semantics, and it formulates the growing of the modules as a submodular
optimization program. Medusa is flexible in choosing or combining semantic
meanings and provides theoretical guarantees about detection quality. In a
systematic study on 310 complex diseases, we show the effectiveness of Medusa
in associating genes with diseases and detecting disease modules. We
demonstrate that in predicting gene-disease associations Medusa compares
favorably to methods that ignore diverse semantic meanings. We find that the
utility of different semantics depends on disease categories and that, overall,
Medusa recovers disease modules more accurately when combining different
semantics.


  Ecosystems display a complex spatial organization. Ecologists have long tried
to characterize them by looking at how different measures of biodiversity
change across spatial scales. Ecological neutral theory has provided simple
predictions accounting for general empirical patterns in communities of
competing species. However, while neutral theory in well-mixed ecosystems is
mathematically well understood, spatial models still present several open
problems, limiting the quantitative understanding of spatial biodiversity. In
this review, we discuss the state of the art in spatial neutral theory. We
emphasize the connection between spatial ecological models and the physics of
non-equilibrium phase transitions and how concepts developed in statistical
physics translate in population dynamics, and vice versa. We focus on
non-trivial scaling laws arising at the critical dimension D = 2 of spatial
neutral models, and their relevance for biological populations inhabiting
two-dimensional environments. We conclude by discussing models incorporating
non-neutral effects in the form of spatial and temporal disorder, and analyze
how their predictions deviate from those of purely neutral theories.


  The shape of allele-frequency clines maintained by migration-selection
balance depends not only on the properties of migration and selection, but also
on the dominance relations among alleles and on linkage to other loci under
selection. We investigate a two-locus model in which two diallelic, recombining
loci are subject to selection caused by an abrupt environmental change. The
habitat is one-dimensional and unbounded, selection at each locus is modeled by
step functions such that in one region one allele at each locus is advantageous
and in the other deleterious. We admit an environmentally independent,
intermediate degree of dominance at both loci, including complete dominance.
First, we derive an explicit expression for the single-locus cline with
dominance, thus generalizing classical results by Haldane (1948). We show that
the slope of the cline in the center (at the step) or, equivalently, the width
of the cline, is independent of the degree of dominance. Second, under the
assumption of strong recombination relative to selection and migration, the
first-order approximations of the allele-frequency clines at each of the loci
and of the linkage disequilibrium are derived. This may be interpreted as the
quasi-linkage-equilibrium approximation of the two-locus cline. Explicit
asymptotic expressions for the clines are deduced as $x\to\pm\infty$. For
equivalent loci, explicit expressions for the whole clines are derived. The
influence of dominance and of linkage on the slope of the cline in the center
and on a global measure of steepness are investigated. This global measure
reflects the influence of dominance. Finally, the accuracy of the
approximations and the dependence of the shape of the two-locus cline on the
full range of recombination rates is explored by numerical integration of the
underlying system of partial differential equations.


  We derive an analytical approximation for making quantitative predictions for
ecological communities as a function of the mean intensity of the
inter-specific competition and the species richness. This method, with only a
fraction of the model parameters (carrying capacities and competition
coefficients), is able to predict accurately empirical measurements covering a
wide variety of taxa (algae, plants, protozoa).


  Surveillance data serving for epidemic alert systems are typically fully
aggregated in space. However, epidemics may be spatially heterogeneous,
undergoing distinct dynamics in distinct regions of the surveillance area. We
unveil this in retrospective analyses by classifying incidence time series. We
use Pearson correlation to quantify the similarity between local time series
and then classify them using modularity maximization. The surveillance area is
thus divided into regions with different incidence patterns. We analyzed 31
years of data on influenza-like-illness from the French system Sentinelles and
found spatial heterogeneity in 19/31 influenza seasons. However, distinct
epidemic regions could be identified only 4-5 weeks after the nationwide alert.
The impact of spatial heterogeneity on influenza epidemiology was complex.
First, when the nationwide alert was triggered, 32-41% of the administrative
regions were experiencing an epidemic, while the others were not. Second, the
nationwide alert was timely for the whole surveillance area, but, subsequently,
regions experienced distinct epidemic dynamics. Third, the epidemic dynamics
were homogeneous in space. Spatial heterogeneity analyses can provide the
timing of the epidemic peak and finish, in various regions, to tailor disease
monitoring and control.


  In evolutionary processes, population structure has a substantial effect on
natural selection. Here, we analyze how motion of individuals affects constant
selection in structured populations. Motion is relevant because it leads to
changes in the distribution of types as mutations march toward fixation or
extinction. We describe motion as the swapping of individuals on graphs, and
more generally as the shuffling of individuals between reproductive updates.
Beginning with a one-dimensional graph, the cycle, we prove that motion
suppresses natural selection for death-birth updating or for any process that
combines birth-death and death-birth updating. If the rule is purely
birth-death updating, no change in fixation probability appears in the presence
of motion. We further investigate how motion affects evolution on the square
lattice and weighted graphs. In the case of weighted graphs we find that motion
can be either an amplifier or a suppressor of natural selection. In some cases,
whether it is one or the other can be a function of the relative reproductive
rate, indicating that motion is a subtle and complex attribute of evolving
populations. As a first step towards understanding less restricted types of
motion in evolutionary graph theory, we consider a similar rule on dynamic
graphs induced by a spatial flow and find qualitatively similar results
indicating that continuous motion also suppresses natural selection.


  Non protein coding regions of the human genome contain many complex patterns
which regulate the cellular activity. Studying the human genome is limited by
the lack of understanding of its features and their complex interactions.
However, recent advances in AI research have enabled automatically learning
representations of high dimensional complex data without feature engineering,
using deep neural networks. Therefore, in this paper, we demonstrate that a
convolutional neural network can learn a representation of DNA sequence without
specifying any motifs or patterns, such that it becomes capable of predicting
whether a DNA sequence is natural or artificial. The trained model could
distinguish scrambled vs real DNA sequences for scrambling lengths of 2 bp, 10
bp, 50 bp and even 100 bp, with a significantly higher accuracy than linear
SVMs. With this study, we have discovered that regions of non protein coding
DNA might have meaningful interactions at even longer than 100 bp distances
even though they do not code proteins.


  Artificial intelligence (AI) tools are gaining more and more ground each year
in bioinformatics. Learning algorithms can be taught easily by using the
existing enormous biological databases, and the resulting models can be used
for the high-quality classification of novel, un-categorized data in numerous
areas, including biological sequence analysis. Here we introduce SECLAF, an
artificial neural-net based biological sequence classifier framework, which
uses the Tensorflow library of Google, Inc. By applying SECLAF for
residue-sequences, we have reported (Methods (2017),
https://doi.org/10.1016/j.ymeth.2017.06.034) the most accurate multi-label
protein classifier to date (UniProt --into 698 classes-- AUC 99.99\%; Gene
Ontology --into 983 classes-- AUC 99.45\%). Our framework SECLAF can be applied
for other sequence classification tasks, as we describe in the present
contribution.
  Availability and implementation: The program SECLAF is implemented in Python,
and is available for download, with example datasets at the website
https://pitgroup.org/seclaf/. For Gene Ontology and UniProt based
classifications a webserver is also available at the address above.


  Motivation: Seed filtering is critical in DNA read mapping, a process where
billions of DNA fragments (reads) sampled from a donor are mapped onto a
reference genome to identify genomic variants of the donor. Read mappers 1)
quickly generate possible mapping locations (i.e., seeds) for each read, 2)
extract reference sequences at each of the mapping locations, and then 3) check
similarity between each read and its associated reference sequences with a
computationally expensive dynamic programming algorithm (alignment) to
determine the origin of the read. Location filters come into play before
alignment, discarding seed locations that alignment would have deemed a poor
match. The ideal location filter would discard all poor matching locations
prior to alignment such that there is no wasted computation on poor alignments.
  Results: We propose a novel filtering algorithm, GRIM-Filter, optimized to
exploit emerging 3D-stacked memory systems that integrate computation within a
stacked logic layer, enabling processing-in-memory (PIM). GRIM-Filter quickly
filters locations by 1) introducing a new representation of coarse-grained
segments of the reference genome and 2) using massively-parallel in-memory
operations to identify read presence within each coarse-grained segment. Our
evaluations show that for 5% error acceptance rates, GRIM-Filter eliminates
5.59x-6.41x more false negatives and exhibits end-to-end speedups of
1.81x-3.65x compared to mappers employing the best previous filtering
algorithm.


  Graphlets are small connected induced subgraphs of a larger graph $G$.
Graphlets are now commonly used to quantify local and global topology of
networks in the field. Methods exist to exhaustively enumerate all graphlets
(and their orbits) in large networks as efficiently as possible using orbit
counting equations. However, the number of graphlets in $G$ is exponential in
both the number of nodes and edges in $G$. Enumerating them all is already
unacceptably expensive on existing large networks, and the problem will only
get worse as networks continue to grow in size and density. Here we introduce
an efficient method designed to aid statistical sampling of graphlets up to
size $k=8$ from a large network. We define graphettes as the generalization of
graphlets allowing for disconnected graphlets. Given a particular (undirected)
graphette $g$, we introduce the idea of the canonical graphette $\mathcal K(g)$
as a representative member of the isomorphism group $Iso(g)$ of $g$. We compute
the mapping $\mathcal K$, in the form of a lookup table, from all
$2^{k(k-1)/2}$ undirected graphettes $g$ of size $k\le 8$ to their canonical
representatives $\mathcal K(g)$, as well as the permutation that transforms $g$
to $\mathcal K(g)$. We also compute all automorphism orbits for each canonical
graphette. Thus, given any $k\le 8$ nodes in a graph $G$, we can in constant
time infer which graphette it is, as well as which orbit each of the $k$ nodes
belongs to. Sampling a large number $N$ of such $k$-sets of nodes provides an
approximation of both the distribution of graphlets and orbits across $G$, and
the orbit degree vector at each node.


  Biomolecular folding, at least in simple systems, can be described as a two
state transition in a free energy landscape with two deep wells separated by a
high barrier. Transition paths are the short part of the trajectories that
cross the barrier. Average transition path times and, recently, their full
probability distribution have been measured for several biomolecular systems,
e.g. in the folding of nucleic acids or proteins. Motivated by these
experiments, we have calculated the full transition path time distribution for
a single stochastic particle crossing a parabolic barrier, focusing on the
underdamped regime. Our analysis thus includes inertial terms, which were
neglected in previous studies. These terms influence the short time scale
dynamics of a stochastic system, and can be of experimental relevance in view
of the short duration of transition paths. We derive the full transition path
time distribution in the underdamped case and discuss the similarities and
differences with the high friction (overdamped) limit.


  Genetically encoded regulatory circuits control biological function. A major
focus of systems biology is to understand these circuits by establishing the
relationship between specific structures and functions. Of special interest are
multifunctional circuits that are capable of performing distinct behaviors
without changing their topology. A particularly simple example of such a system
is the AC-DC circuit. Found in multiple regulatory processes, this circuit
consists of three genes connected in a combination of a toggle switch and a
repressilator. Using dynamical system theory we analyze the available dynamical
regimes to show that the AC-DC can exhibit both oscillations and bistability.
We found that both dynamical regimes can coexist robustly in the same region of
parameter space, generating novel emergent behaviors not available to the
individual subnetwork components. We demonstrate that the AC-DC circuit
provides a mechanism to rapidly switch between oscillations and steady
expression and, in the presence of noise, the multi-functionality of the
circuit offers the possibility to control the coherence of oscillations.
Additionally, we provide evidence that the availability of a bistable
oscillatory regime allows the AC-DC circuit to behave as an excitable system
capable of stochastic pulses and spatial signal propagation. Taken together
these results reveal how a system as simple as the AC-DC circuit can produce
multiple complex dynamical behaviors in the same parameter region and is well
suited for the construction of multifunctional synthetic genetic circuits.
Likewise, the analysis reveals the potential of this circuit to facilitate the
evolution of distinct patterning mechanisms.


  Actin filaments are critical components of the eukaryotic cytoskeleton,
playing important roles in a number of cellular functions, such as cell
migration, organelle transport, and mechanosensation. They are helical polymers
with a well-defined polarity, composed of globular monomers that bind
nucleotides in one of three hydrolysis states (ATP, ADP-Pi, or ADP). Mean-field
models of the dynamics of actin polymerization have succeeded in, among other
things, determining the nucleotide profile of an average filament and resolving
the mechanisms of accessory proteins, however these models require numerical
solution of a high-dimensional system of nonlinear ODE's. By truncating a set
of recursion equations, the Brooks-Carlsson model reduces dimensionality to 11,
but it remains nonlinear and does not admit an analytical solution, hence,
significantly hindering understanding of its resulting dynamics. In this work,
by taking advantage of the fast timescales of the hydrolysis states of the
filament tips, we propose two model reduction schemes that achieve low
dimensionality and linearity. We provide an exact solution of the resulting
linear equations and use it to shed light on the dynamical behaviors of the
full BC model, highlighting the relative ordering of the timescales of various
collective processes, and explaining some unusual dependence of the
steady-state behavior on initial conditions.


  Template directed replication of nucleic acids is at the essence of all
living beings and a major milestone for any origin of life scenario. We here
present an idealized model of prebiotic sequence replication, where binary
polymers act as templates for their autocatalytic replication, thereby serving
as each others reactants and products in an intertwined molecular ecology. Our
model demonstrates how autocatalysis alters the qualitative and quantitative
system dynamics in counter-intuitive ways. Most notably, numerical simulations
reveal a very strong intrinsic selection mechanism that favours the appearance
of a few population structures with highly ordered and repetitive sequence
patterns when starting from a pool of monomers. We demonstrate both
analytically and through simulation how this "selection of the dullest" is
caused by continued symmetry breaking through random fluctuations in the
transient dynamics that are amplified by autocatalysis and eventually propagate
to the population level. The impact of these observations on related prebiotic
mathematical models is discussed.


  Hydroxymethylcytosine (5hmC) methylation is known to be a possible epigenetic
mark impacting genome stability. In this paper, we address the existing 5hmC
measure $\Delta \beta$ and discuss its properties both analytically and
empirically on real data. Then we introduce several alternative
hydroxymethylation measures and compare their properties with those of $\Delta
\beta$. All results will be illustrated by means of real data analyses.