  Phylogenetic networks extend phylogenetic trees to allow for modeling
reticulate evolutionary processes such as hybridization. They take the shape of
a rooted, directed, acyclic graph, and when parameterized with evolutionary
parameters, such as divergence times and population sizes, they form a
generative process of molecular sequence evolution. Early work on computational
methods for phylogenetic network inference focused exclusively on reticulations
and sought networks with the fewest number of reticulations to fit the data. As
processes such as incomplete lineage sorting (ILS) could be at play
concurrently with hybridization, work in the last decade has shifted to
computational approaches for phylogenetic network inference in the presence of
ILS. In such a short period, significant advances have been made on developing
and implementing such computational approaches. In particular, parsimony,
likelihood, and Bayesian methods have been devised for estimating phylogenetic
networks and associated parameters using estimated gene trees as data. Use of
those inference methods has been augmented with statistical tests for specific
hypotheses of hybridization, like the D-statistic. Most recently, Bayesian
approaches for inferring phylogenetic networks directly from sequence data were
developed and implemented. In this chapter, we survey such advances and discuss
model assumptions as well as methods' strengths and limitations. We also
discuss parallel efforts in the population genetics community aimed at
inferring similar structures. Finally, we highlight major directions for future
research in this area.


  Time varying susceptibility of host at individual level due to waning and
boosting immunity is known to induce rich long-term behavior of disease
transmission dynamics. Meanwhile, the impact of the time varying heterogeneity
of host susceptibility on the shot-term behavior of epidemics is not
well-studied, even though the large amount of the available epidemiological
data are the short-term epidemics. Here we constructed a parsimonious
mathematical model describing the short-term transmission dynamics taking into
account natural-boosting immunity by reinfection, and obtained the explicit
solution for our model. We found that our system show "the delayed epidemic",
the epidemic takes off after negative slope of the epidemic curve at the
initial phase of epidemic, in addition to the common classification in the
standard SIR model, i.e., "no epidemic" as $\mathcal{R}_{0}\leq1$ or normal
epidemic as $\mathcal{R}_{0}>1$. Employing the explicit solution we derived the
condition for each classification.


  Studies of voluntary vaccination decisions by rational individuals predict
that the population will reach a Nash equilibrium with vaccination coverage
below the societal optimum. Human decision-making involves mechanisms in
addition to rational calculations of self-interest, such as imitation of
successful others. Previous research had shown that imitation alone cannot
achieve better results. Under realistic choices of the parameters it may lead
to equilibrium vaccination coverage even below the Nash equilibrium. However,
these findings rely on the widely accepted use of Fermi functions for modeling
the probabilities of switching to another strategy. We consider here a more
general functional form of the switching probabilities. It is consistent with
functions that give best fits for empirical data in a widely cited
psychological experiment and involves one additional parameter $\alpha$. This
parameter can be loosely interpreted as a degree of open-mindedness. We found
both by means of simulations and analytically that sufficiently high values of
$\alpha$ will drive the equilibrium vaccination coverage arbitrarily close to
the societal optimum.


  Auxin is considered one of the cardinal hormones in plant growth and
development. It regulates a wide range of processes throughout the plant.
Synthetic auxins exploit the auxin-signalling pathway and are valuable as
herbicidal agrochemicals. Currently, despite a diversity of chemical scaffolds
all synthetic auxins have a carboxylic acid as the active core group. By
applying bio-isosteric replacement we discovered that indole-3-tetrazole was
active by surface plasmon resonance (SPR) spectrometry, showing that the
tetrazole could initiate assembly of the TIR1 auxin co-receptor complex. We
then tested the tetrazole's efficacy in a range of whole plant physiological
assays and in protoplast reporter assays which all confirmed auxin activity,
albeit rather weak. We then tested indole-3-tetrazole against the AFB5
homologue of TIR1, finding that binding was selective against TIR1, absent with
AFB5. The kinetics of binding to TIR1 are contrasted to those for the herbicide
picloram, which shows the opposite receptor preference as it binds to AFB5 with
far greater affinity than to TIR1. The basis of the preference of
indole-3-tetrazole for TIR1 was revealed to be a single residue substitution
using molecular docking, and assays using tir1 and afb5 mutant lines confirmed
selectivity in vivo. Given the potential that a TIR1-selective auxin might have
for unmasking receptor-specific actions, we followed a rational design, lead
optimisation campaign and a set of chlorinated indole-3-tetrazoles was
synthesised. Improved affinity for TIR1 and the preference for binding to TIR1
was maintained for 4- and 6-chloroindole-3-tetrazoles, coupled with improved
efficacy in vivo. This work expands the range of auxin chemistry for the design
of receptor-selective synthetic auxins.


  Power-law-distributed species counts or clone counts arise in many biological
settings such as multispecies cell populations, population genetics, and
ecology. This empirical observation that the number of species $c_{k}$
represented by $k$ individuals scales as negative powers of $k$ is also
supported by a series of theoretical birth-death-immigration (BDI) models that
consistently predict many low-population species, a few intermediate-population
species, and very high-population species. However, we show how a simple global
population-dependent regulation in a neutral BDI model destroys the power law
distributions. Simulation of the regulated BDI model shows a high probability
of observing a high-population species that dominates the total population.
Further analysis reveals that the origin of this breakdown is associated with
the failure of a mean-field approximation for the expected species abundance
distribution. We find an accurate estimate for the expected distribution
$\langle c_k \rangle$ by mapping the problem to a lower-dimensional Moran
process, allowing us to also straightforwardly calculate the covariances
$\langle c_k c_\ell \rangle$. Finally, we exploit the concepts associated with
energy landscapes to explain the failure of the mean-field assumption by
identifying a phase transition in the quasi-steady-state species counts
triggered by a decreasing immigration rate.


  We present a new method that combines alchemical transformation with physical
pathway to accurately and efficiently compute the absolute binding free energy
of receptor-ligand complex. Currently, the double decoupling method (DDM) and
the potential of mean force approach (PMF) methods are widely used to compute
the absolute binding free energy of biomolecules. The DDM relies on
alchemically decoupling the ligand from its environments, which can be
computationally challenging for large ligands and charged ligands because of
the large magnitude of the decoupling free energies involved. On the other
hand, the PMF approach uses physical pathway to extract the ligand out of the
binding site, thus avoids the alchemical decoupling of the ligand. However, the
PMF method has its own drawback because of the reliance on a ligand
binding/unbinding pathway free of steric obstruction from the receptor atoms.
Therefore, in the presence of deeply buried ligand functional groups the
convergence of the PMF calculation can be very slow leading to large errors in
the computed binding free energy. Here we develop a new method called AlchemPMF
by combining alchemical transformation with physical pathway to overcome the
major drawback in the PMF method. We have tested the new approach on the
binding of a charged ligand to an allosteric site on HIV-1 Integrase. After 20
ns of simulation per umbrella sampling window, the new method yields absolute
binding free energies within ~1 kcal/mol from the experimental result, whereas
the standard PMF approach and the DDM calculations result in errors of ~5
kcal/mol and > 2 kcal/mol, respectively. Furthermore, the binding free energy
computed using the new method is associated with smaller statistical error
compared with those obtained from the existing methods.


  Purpose of review: The epithelial-Mesenchymal Transition (EMT) and the
generation of Cancer Stem Cells (CSC) are two fundamental aspects contributing
to tumor growth, acquisition of resistance to therapy, formation of metastases,
and tumor relapse. Recent experimental data identifying the circuits regulating
EMT and CSCs have driven the development of computational models capturing the
dynamics of these circuits and consequently various aspects of tumor
progression. Recent findings: We review the contribution made by these models
in a) recapitulating experimentally observed behavior, b) making experimentally
testable predictions, and c) driving emerging notions in the field, including
the emphasis on the aggressive potential of hybrid epithelial/mesenchymal (E/M)
phenotype(s). We discuss dynamical and statistical models at intracellular and
population level relating to dynamics of EMT and CSCs, and those focusing on
interconnections between these two processes. Summary: These models highlight
the insights gained via mathematical modeling approaches, and emphasizes that
the connections between hybrid E/M phenotype(s) and stemness can be explained
by analyzing underlying regulatory circuits. Such experimentally curated models
have the potential of serving as platforms for better therapeutic design
strategies.


  Cells are complex structures which require considerable amounts of
organization via transport of large intracellular cargo. While passive
diffusion is often sufficiently fast for the transport of smaller cargo, active
transport is necessary to organize large structures on short timescales. The
main mechanism of this transport is by cargo attachment to motors which walk in
a directed fashion along intracellular filaments. There are a number of models
which seek to describe the motion of motors with attached cargo, from detailed
microscopic to coarse phenomenological descriptions. We focus on the
intermediate-detailed discrete stochastic hopping models, and explore how cargo
transport changes depending on the number of motors, motor interaction, system
constraints and rate formulations which are derived from common thermodynamic
assumptions. We find that, despite obeying the same detailed balance
constraint, the choice of rate formulation considerably affects the
characteristics of the overall motion of the system, with one rate formulation
exhibiting novel behavior of loaded motor groups moving faster than a single
unloaded motor.


  The search for causative pharmacogenomic loci is being transformed by
integrative omics pipelines, but their outputs have only begun being applied to
test design. We assess the direction of the field in light of Biobanks/PheWAS,
omics atlases, and AI. We first assess the potential of recent epigenome and
spatial genome concepts, datasets, and methods to improve the functionality of
PIP-style pipelines. We then discuss new potential methods of genetic test
design on the basis of the outputs of such pipelines. We conclude with a vision
for a pharmacophenomic atlas, in which omics atlas data, PheWAS associations,
and biobank data would be used with AI to design thousands of genetic tests for
clinical deployment in an automated parallel process.


  We propose an experiment to demonstrate spontaneous ordering and symmetry
breaking of kinesin-driven microtubules confined to an optical trap.
Calculations involving the feasibility of such an experiment are first
performed which analyze the power needed to confine microtubules and address
heating concerns. We then present the results of first-principles simulations
of active microtubules confined in such a trap and analyze the types of motion
observed by the microtubules as well as the velocity of the surrounding fluid,
both near the trap and in the far-field. We find three distinct phases
characterized by breaking of distinct symmetries and also analyze the power
spectrum of the angular momenta of polymers to further quantify the differences
between these phases. Under the correct conditions, microtubules were found to
spontaneously align with one another and circle the trap in one direction.


  Inside cells of living organisms, actin filaments and microtubules
self-assemble and dissemble dynamically by incorporating actin or tubulin from
the cell plasma or releasing it into their tips' surroundings. Such
reaction-diffusion systems can show diffusion- or reaction-limited behaviour.
However, neither limit explains the experimental data: while the offset of the
linear relation between growth speed and bulk tubulin density contradicts the
diffusion limit, the surprisingly large variance of the growth speed rejects a
pure reaction limit. In this Letter, we accommodate both limits and use a
Doi-Peliti field-theory model to estimate how diffusive transport is perturbing
the chemical reactions at the filament tip. Furthermore, a crossover bulk
density is predicted at which the limiting process changes from chemical
reactions to diffusive transport. In addition, we explain and estimate larger
variances of the growth speed.


  A class of stochastic vector-borne infectious disease models is derived and
studied. The class type is determined by a general nonlinear incidence rate of
the disease. The disease spreads in a highly random environment with
variability from the disease transmission and natural death rates. Other
sources of variability include the random delays of disease incubation inside
the vector and the human being, and also the random delay due to the period of
effective acquired immunity against the disease. The basic reproduction number
and other threshold conditions for disease eradication are computed. The
qualitative behaviors of the disease dynamics are examined under the different
sources of variability in the system. A technique to classify the different
levels of the intensities of the noises in the system is presented, and used to
investigate the qualitative behaviors of the disease dynamics in the
infection-free steady state population under the different intensity levels of
the white noises in the system. Moreover, the possibility of population
extinction, whenever the intensities of the white noises in the system are high
is examined. Numerical simulation results are presented to elucidate the
theoretical results.


  Evolutionary game theory has been widely used to study the evolution of
cooperation in social dilemmas where imitation-led strategy updates are
typically assumed. However, results of recent behavioral experiments are not
compatible with the predictions based on the imitation dynamics, casting doubts
on the assumption of the imitation-led updates and calling for alternative
mechanisms of strategy updates. An aspiration-led update is often considered as
an alternative to the imitation-led one. While details of update rules can have
significant impacts on the evolutionary outcomes and many variations in
imitation-led updates are thus studied, however, few variations exist in
aspiration-led updates. We introduce an aspiration-led update mechanism
("Satisfied-Defect, Unsatisfied-Cooperate") that is psychologically intuitive
and can yield a behavior richer than the conventional aspiration-led update
does in Prisoner's Dilemma games. Using analytical and numerical methods, we
study and link the stochastic dynamics of it in well-mixed finite populations
and the deterministic dynamics of infinite populations.


  Biomolecular condensates undergirded by phase separations of proteins and
nucleic acids serve crucial biological functions. To gain physical insights
into their genetic basis, we study how liquid-liquid phase separation (LLPS) of
intrinsically disordered proteins (IDPs) depends on their sequence charge
patterns using a continuum Langevin chain model wherein each amino acid residue
is represented by a single bead. Charge patterns are characterized by the
`blockiness' measure $\kappa$ and the `sequence charge decoration' (SCD)
parameter. Consistent with random phase approximation (RPA) theory and lattice
simulations, LLPS propensity as characterized by critical temperature $T^*_{\rm
cr}$ increases with increasingly negative SCD for a set of sequences showing a
positive correlation between $\kappa$ and $-$SCD. Relative to RPA, the
simulated sequence-dependent variation in $T^*_{\rm cr}$ is often---though not
always---smaller, whereas the simulated critical volume fractions are higher.
However, for a set of sequences exhibiting an anti-correlation between $\kappa$
and $-$SCD, the simulated $T^*_{\rm cr}$'s are quite insensitive to either
parameters. Additionally, we find that blocky sequences that allow for strong
electrostatic repulsion can lead to coexistence curves with upward concavity as
stipulated by RPA, but the LLPS propensity of a strictly alternating charge
sequence was likely overestimated by RPA and lattice models because interchain
stabilization of this sequence requires spatial alignments that are difficult
to achieve in real space. These results help delineate the utility and
limitations of the charge pattern parameters and of RPA, pointing to further
efforts necessary for rationalizing the newly observed subtleties.


  Experiments reveal that DNA solvation dynamics (SD) is characterized by
multiple time scales ranging from a few ps to hundreds of ns and in some cases
even up to microseconds. The last part of decay is slow and is characterized by
a power law (PL). The microscopic origin of this PL is not yet clearly
understood. Here we present a theoretical study based on time dependent
statistical mechanics and computer simulations. Our investigations show that
the primary candidates responsible for this exotic nature of SD are the
counterions and ions from the buffer solution. We employ the model developed by
Oosawa for polyelectrolyte solution that includes effects of counterion
fluctuations to construct a frequency dependent dielectric function. We use it
in the continuum model of Bagchi, Fleming and Oxtoby only to find that it fails
to explain the slow PL decay of DNA solvation dynamics. We then extend the
model by employing the continuous time random walk technique developed by
Scher-Montroll-Lax. This approach can explain the long time PL decay, in terms
of the collective response of the counter ions. From MD simulations we find
frequent occurrence of random walk of tagged counter ions along the phosphate
backbone. We propose a generalized random walk model for counterion hopping and
carry out kinetic Monte Carlo simulations to show that the nonexponential
contribution to solvation dynamics can indeed arise from dynamics of such ions.
We also employ a Mode Coupling Theory analysis to understand the slow
relaxation that originates from ions in solution. Explicit evaluation suggests
that buffer ion contribution could explain logarithmic time dependence in the
ns time scale, but not a power law. From MD simulations we find log-normal
distributions of relaxation times of water dynamics inside the grooves. This is
responsible for the initial faster multiexponential decay of SD.


  With approximately half of the world's population at risk of contracting
dengue, this mosquito-borne disease is of global concern. International
travellers significantly contribute to dengue's rapid and large-scale spread by
importing the disease from endemic into non-endemic countries. To prevent
future outbreaks and dengue from establishing in non-endemic countries,
knowledge about the arrival time and location of infected travellers is
crucial. We propose a network model that predicts the monthly number of
dengue-infected air passengers arriving at any given airport. We consider
international air travel volumes to construct weighted networks, representing
passenger flows between airports. We further calculate the probability of
passengers, who travel through the international air transport network, being
infected with dengue. The probability of being infected depends on the
destination, duration and timing of travel. Our findings shed light onto dengue
importation routes and reveal country-specific reporting rates that have been
until now largely unknown. This paper provides important new knowledge about
the spreading dynamics of dengue that is highly beneficial for public health
authorities to strategically allocate the often limited resources to more
efficiently prevent the spread of dengue.


  We study the evolution of the population genealogy in the classic neutral
Moran Model of finite size and in discrete time. The stochastic transformations
that shape a Moran population can be realized directly on its genealogy and
give rise to a process with a state space consisting of the finite set of Yule
trees of a certain size. We derive a number of properties of this process, and
show that they are in agreement with existing results on the
infinite-population limit of the Moran Model. Most importantly, this process
admits time reversal, which gives rise to another tree-valued Markov Chain and
allows for a thorough investigation of the Most Recent Common Ancestor process.


  Accurate prediction of inter-residue contacts of a protein is important to
calcu- lating its tertiary structure. Analysis of co-evolutionary events among
residues has been proved effective to inferring inter-residue contacts. The
Markov ran- dom field (MRF) technique, although being widely used for contact
prediction, suffers from the following dilemma: the actual likelihood function
of MRF is accurate but time-consuming to calculate, in contrast, approximations
to the actual likelihood, say pseudo-likelihood, are efficient to calculate but
inaccu- rate. Thus, how to achieve both accuracy and efficiency simultaneously
remains a challenge. In this study, we present such an approach (called clmDCA)
for contact prediction. Unlike plmDCA using pseudo-likelihood, i.e., the
product of conditional probability of individual residues, our approach uses
composite- likelihood, i.e., the product of conditional probability of all
residue pairs. Com- posite likelihood has been theoretically proved as a better
approximation to the actual likelihood function than pseudo-likelihood.
Meanwhile, composite likelihood is still efficient to maximize, thus ensuring
the efficiency of clmDCA. We present comprehensive experiments on popular
benchmark datasets, includ- ing PSICOV dataset and CASP-11 dataset, to show
that: i) clmDCA alone outperforms the existing MRF-based approaches in
prediction accuracy. ii) When equipped with deep learning technique for
refinement, the prediction ac- curacy of clmDCA was further significantly
improved, suggesting the suitability of clmDCA for subsequent refinement
procedure. We further present successful application of the predicted contacts
to accurately build tertiary structures for proteins in the PSICOV dataset.
  Accessibility: The software clmDCA and a server are publicly accessible
through http://protein.ict.ac.cn/clmDCA/.


  The expansion of a population into new habitat is a transient process that
leaves its footprints in the genetic composition of the expanding population.
How the structure of the environment shapes the population front and the
evolutionary dynamics during such a range expansion is little understood. Here,
we investigate the evolutionary dynamics of populations consisting of many
selectively neutral genotypes expanding on curved surfaces. Using a combination
of individual-based off-lattice simulations, geometrical arguments, and
lattice-based stepping-stone simulations, we characterise the effect of
individual bumps on an otherwise flat surface. Compared to the case of a range
expansion on a flat surface, we observe a transient relative increase, followed
by a decrease, in neutral genetic diversity at the population front. In
addition, we find that individuals at the sides of the bump have a dramatically
increased expected number of descendants, while their neighbours closer to the
bump's centre are far less lucky. Both observations can be explained using an
analytical description of straight paths (geodesics) on the curved surface.
Complementing previous studies of heterogeneous flat environments, the findings
here build our understanding of how complex environments shape the evolutionary
dynamics of expanding populations.


  Urban areas with larger and more connected populations offer an auspicious
environment for contagion processes such as the spread of pathogens. Empirical
evidence reveals a systematic increase in the rates of certain sexually
transmitted diseases (STDs) with larger urban population size. However, the
main drivers of these systemic infection patterns are still not well
understood, and rampant urbanization rates worldwide makes it critical to
advance our understanding on this front. Using confirmed-cases data for three
STDs in US metropolitan areas, we investigate the scaling patterns of
infectious disease incidence in urban areas. The most salient features of these
patterns are that, on average, the incidence of infectious diseases that
transmit with less ease-- either because of a lower inherent transmissibility
or due to a less suitable environment for transmission-- scale more steeply
with population size, are less predictable across time and more variable across
cities of similar size. These features are explained, first, using a simple
mathematical model of contagion, and then through the lens of a new theory of
urban scaling. These theoretical frameworks help us reveal the links between
the factors that determine the transmissibility of infectious diseases and the
properties of their scaling patterns across cities.


  If a system undergoing phase transitions exhibits some characteristics of
both first and second order, it is said to be of 'mixed order' or to display
the Thouless effect. Such a transition is present in a simple model of a
dynamic social network, in which $N_{I/E}$ extreme introverts/extroverts always
cut/add random links. In particular, simulations showed that $\left\langle
f\right\rangle $, the average fraction of cross-links between the two groups
(which serves as an 'order parameter' here), jumps dramatically when $\Delta
\equiv N_{I}-N_{E}$ crosses the 'critical point' $\Delta _{c}=0$, as in typical
first order transitions. Yet, at criticality, there is no phase co-existence,
but the fluctuations of $f$ are much larger than in typical second order
transitions. Indeed, it was conjectured that, in the thermodynamic limit, both
the jump and the fluctuations become maximal, so that the system is said to
display an 'extreme Thouless effect.' While earlier theories are partially
successful, we provide a mean-field like approach that accounts for all known
simulation data and validates the conjecture. Moreover, for the critical system
$N_{I}=N_{E}=L$, an analytic expression for the mesa-like stationary
distribution, $P\left( f\right) $, shows that it is essentially flat in a range
$\left[ f_{0},1-f_{0}\right] $, with $f_0 \ll 1$. Numerical evaluations of
$f_{0}$ provides excellent agreement with simulation data for $L\lesssim 2000$.
For large $L$, we find $f_{0}\rightarrow \sqrt{\left( \ln L^2 \right) /L}$ ,
though this behavior begins to set in only for $L>10^{100}$. For accessible
values of $L$, we provide a transcendental equation for an approximate $f_{0}$
which is better than $\sim$1% down to $L=100$. We conjecture how this approach
might be used to attack other systems displaying an extreme Thouless effect.


  The Bond Graph approach and the Chemical Reaction Network approach to
modelling biomolecular systems developed independently. This paper brings
together the two approaches by providing a bond graph interpretation of the
chemical reaction network concept of complexes. Both closed and open systems
are discussed.
  The method is illustrated using a simple enzyme-catalysed reaction and a
trans-membrane transporter.


  Clostridium difficile infections (CDIs) affect patients in hospitals and in
the community, but the relative importance of transmission in each setting is
unknown. We developed a mathematical model of C. difficile transmission in a
hospital and surrounding community that included infants, adults, and
transmission from animal reservoirs. We assessed the role of these transmission
routes in maintaining disease and evaluated the recommended classification
system for hospital and community-acquired CDIs. The reproduction number in the
hospital was <1 (range: 0.16-0.46) for all scenarios. Outside the hospital, the
reproduction number was >1 for nearly all scenarios without transmission from
animal reservoirs (range: 1.0-1.34). However, the reproduction number for the
human population was <1 if a minority (>3.5-26.0%) of human exposures
originated from animal reservoirs. Symptomatic adults accounted for <10%
transmission in the community. Under conservative assumptions, infants
accounted for 17% of community transmission. An estimated 33-40% of
community-acquired cases were reported but 28-39% of these reported cases were
misclassified as hospital-acquired by recommended definitions. Transmission
could be plausibly sustained by asymptomatically colonized adults and infants
in the community or exposure to animal reservoirs, but not hospital
transmission alone. Underreporting of community-onset cases and systematic
misclassification underplays the role of community transmission.


  Alcohol is widely believed to have an effect on diabetes, often considered
beneficial in small amounts but detrimental in excess. The reasons are not
fully known but questions have been asked about the stability of insulin
oligomers in the presence of ethanol. We compute the free energy surface (FES)
and the pathway of insulin dimer dissociation in water and in 5% and 10%
water-ethanol mixture. We find that in the presence of ethanol the barrier
energy of dissociation reaction decreases by about 40% even in 5% water-ethanol
solution. In addition, ethanol induces a significant change in the reaction
pathway. We obtain estimates of the rate of reaction and binding energy for all
the three systems and those agree well with the previous experimental results
for the insulin dimer dissociation in water. The computed FES in water exhibits
ruggedness due to the existence of a number of intermediate states surrounded
by high and broad transition state region. However, the presence of ethanol
smoothens out the ruggedness. We extracted the conformations of the
intermediate states along the minimum energy pathway in all the three systems
and analyzed the change in microscopic structures in the presence of ethanol.
Interestingly, we discover a stable intermediate state in water-ethanol
mixtures where the monomers are separated (center-to-center) by about 3 nm and
the contact order parameter is close to zero. This intermediate is stabilized
by the distribution of ethanol and water molecules at the interface and which,
significantly, serves to reduce the dissociation rate constant .The solvation
of the two monomers during the dissociation and proteins' departure from native
state configuration are analyzed to obtain insight into the dimer dissociation
processes.


  The basic reproduction number ($R_0$) is a threshold parameter for disease
extinction or survival in isolated populations. However no human population is
fully isolated from other human or animal populations. We use compartmental
models to derive simple rules for the basic reproduction number for populations
with local person-to-person transmission and exposure from some other source:
either a reservoir exposure or imported cases. We introduce the idea of a
reservoir-driven or importation-driven disease: diseases that would become
extinct in the population of interest without reservoir exposure or imported
cases (since $R_0<1$), but nevertheless may be sufficiently transmissible that
many or most infections are acquired from humans in that population. We show
that in the simplest case, $R_0<1$ if and only if the proportion of infections
acquired from the external source exceeds the disease prevalence and explore
how population heterogeneity and the interactions of multiple strains affect
this rule. We apply these rules in two cases studies of Clostridium difficile
infection and colonisation: C. difficile in the hospital setting accounting for
imported cases, and C. difficile in the general human population accounting for
exposure to animal reservoirs. We demonstrate that even the hospital-adapted,
highly-transmissible NAP1/RT027 strain of C. difficile had a reproduction
number <1 in a landmark study of hospitalised patients and therefore was
sustained by colonised and infected admissions to the study hospital. We argue
that C. difficile should be considered reservoir-driven if as little as 13.0%
of transmission can be attributed to animal reservoirs.


  DNA read mapping is a computationally expensive bioinformatics task, required
for genome assembly and consensus polishing. It requires to find the
best-fitting location for each DNA read on a long reference sequence. A novel
resistive approximate similarity search accelerator, RASSA, exploits charge
distribution and parallel in-memory processing to reflect a mismatch count
between DNA sequences. RASSA implementation of DNA long read pre-alignment
outperforms the state-of-art solution, minimap2, by 16-77x with comparable
accuracy and provides two orders of magnitude higher throughput than
GateKeeper, a short-read pre-alignment hardware architecture implemented in
FPGA.


  Hydrogen bonding between nucleobases produces diverse DNA structural motifs,
including canonical duplexes, guanine (G) quadruplexes and cytosine (C)
i-motifs. Incorporating metal-mediated base pairs into nucleic acid structures
can introduce new functionalities and enhanced stabilities. Here we
demonstrate, using mass spectrometry (MS), ion mobility spectrometry (IMS) and
fluorescence resonance energy transfer (FRET), that parallel-stranded
structures consisting of up to 20 G-Ag(I)-G contiguous base pairs are formed
when natural DNA sequences are mixed with silver cations in aqueous solution.
FRET indicates that duplexes formed by poly(cytosine) strands with 20
contiguous C-Ag(I)-C base pairs are also parallel. Silver-mediated G-duplexes
form preferentially over G-quadruplexes, and the ability of Ag+ to convert
G-quadruplexes into silver-paired duplexes may provide a new route to
manipulating these biologically relevant structures. IMS indicates that
G-duplexes are linear and more rigid than B-DNA. DFT calculations were used to
propose structures compatible with the IMS experiments. Such inexpensive,
defect-free and soluble DNA-based nanowires open new directions in the design
of novel metal-mediated DNA nanotechnology.


  Diffusion in cell biology is important and complicated. Diffusing particles
must contend with a complex environment as they make their way through the
cell. We analyze a particular type of complexity that arises when diffusing
particles reversibly bind to elastically tethered binding partners. Using
asymptotic analysis, we derive effective equations for the transport of both
single and multiple particles in the presence of such elastic tethers. We show
that for the case of linear elasticity and simple binding kinetics, the elastic
tethers have a weak hindering effect on particle motion when only one particle
is present, while, remarkably, strongly enhancing particle motion when multiple
particles are present. We give a physical interpretation of this result that
suggests a similar effect may be present in other biological settings.


  Topological entanglements severely interfere with important biological
processes. For this reason, genomes must be kept unknotted and unlinked during
most of a cell cycle. Type II Topoisomerase (TopoII) enzymes play an important
role in this process but the precise mechanisms yielding systematic
disentanglement of DNA in vivo are not clear. Here we report computational
evidence that Structural Maintenance of Chromosomes (SMC) proteins -- such as
cohesins and condensins -- can cooperate with TopoII to establish a synergistic
mechanism to resolve topological entanglements. SMC-driven loop extrusion (or
diffusion) induces the spatial localisation of essential crossings in turn
catalysing the simplification of knots and links by TopoII enzymes even in
crowded and confined conditions. The mechanism we uncover is universal in that
it does not qualitatively depend on the specific substrate, whether DNA or
chromatin, or on SMC processivity; we thus argue that this synergy may be at
work across organisms and throughout the cell cycle.


  Arboviral infections, especially dengue, continue to cause significant health
burden in their endemic regions. One of the strategies to tackle these
infections is to replace the main vector agent, Ae. aegypti, with the ones
incapable of transmitting the virus. Wolbachia, an intracellular bacterium, has
shown promise in achieving this goal. However, key factors such as imperfect
maternal transmission, loss of Wolbachia infection, reduced reproductive
capacity and shortened life-span affect the dynamics of Wolbachia in different
forms in the Ae. aegypti population. In this study, we developed a Wolbachia
transmission dynamic model adjusting for imperfect maternal transmission and
loss of Wolbachia infection. The invasive reproductive number that determines
the likelihood of replacement of the Wolbachia-uninfected (WU) population is
derived and with it, we established the local and global stability of the
equilibrium points. This analysis clearly shows that cytoplasmic
incompatibility (CI) does not guarantee establishment of the Wolbachia-infected
(WI) mosquitoes as imperfect maternal transmission and loss of Wolbachia
infection could outweigh the gains from CI. Optimal release programs depending
on the level of imperfect maternal transmission and loss of Wolbachia infection
are shown. Hence, it is left to decision makers to either aim for replacement
or co-existence of both populations.


  Identifying significant subsets of the genes, gene shaving is an essential
and challenging issue for biomedical research for a huge number of genes and
the complex nature of biological networks,. Since positive definite kernel
based methods on genomic information can improve the prediction of diseases, in
this paper we proposed a new method, "kernel gene shaving (kernel canonical
correlation analysis (kernel CCA) based gene shaving). This problem is
addressed using the influence function of the kernel CCA. To investigate the
performance of the proposed method in a comparison of three popular gene
selection methods (T-test, SAM and LIMMA), we were used extensive simulated and
real microarray gene expression datasets. The performance measures AUC was
computed for each of the methods. The achievement of the proposed method has
improved than the three well-known gene selection methods. In real data
analysis, the proposed method identified a subsets of $210$ genes out of $2000$
genes. The network of these genes has significantly more interactions than
expected, which indicates that they may function in a concerted effort on colon
cancer.


  According to the evolutionary game theory principle, a strategy representing
a higher payoff can spread among competitors. But there are cases when a player
consistently overestimates or underestimates her own payoff, which undermines
proper comparison. Interestingly, both underconfident and overconfident
individuals are capable of elevating the cooperation level significantly. While
former players stimulate a local coordination of strategies, the presence of
overconfident individuals enhances the spatial reciprocity mechanism. In both
cases the propagations of competing strategies are influenced in a biased way
resulting in a cooperation supporting environment. These effects are strongly
related to the nonlinear character of invasion probabilities which is a common
and frequently observed feature of microscopic dynamics.


  Data commons collate data with cloud computing infrastructure and commonly
used software services, tools and applications to create biomedical resources
for the large-scale management, analysis, harmonization, and sharing of
biomedical data. Over the past few years, data commons have been used to
analyze, harmonize and share large scale genomics datasets. Data ecosystems can
be built by interoperating multiple data commons. It can be quite labor
intensive to curate, import and analyze the data in a data commons. Data lakes
provide an alternative to data commons and simply provide access to data, with
the data curation and analysis deferred until later and delegated to those that
access the data. We review software platforms for managing, analyzing and
sharing genomic data, with an emphasis on data commons, but also covering data
ecosystems and data lakes.


  Cell-to-cell heterogeneity drives a range of (patho)physiologically important
phenomena, such as cell fate and chemotherapeutic resistance. The role of
metabolism, and particularly mitochondria, is increasingly being recognised as
an important explanatory factor in cell-to-cell heterogeneity. Most eukaryotic
cells possess a population of mitochondria, in the sense that mitochondrial DNA
(mtDNA) is held in multiple copies per cell, where the sequence of each
molecule can vary. Hence intra-cellular mitochondrial heterogeneity is
possible, which can induce inter-cellular mitochondrial heterogeneity, and may
drive aspects of cellular noise. In this review, we discuss sources of
mitochondrial heterogeneity (variations between mitochondria in the same cell,
and mitochondrial variations between supposedly identical cells) from both
genetic and non-genetic perspectives, and mitochondrial genotype-phenotype
links. We discuss the apparent homeostasis of mtDNA copy number, the
observation of pervasive intra-cellular mtDNA mutation (we term
`microheteroplasmy') and developments in the understanding of inter-cellular
mtDNA mutation (`macroheteroplasmy'). We point to the relationship between
mitochondrial supercomplexes, cristal structure, pH and cardiolipin as a
potential amplifier of the mitochondrial genotype-phenotype link. We also
discuss mitochondrial membrane potential and networks as sources of
mitochondrial heterogeneity, and their influence upon the mitochondrial genome.
Finally, we revisit the idea of mitochondrial complementation as a means of
dampening mitochondrial genotype-phenotype links in light of recent
experimental developments. The diverse sources of mitochondrial heterogeneity,
as well as their increasingly recognised role in contributing to cellular
heterogeneity, highlights the need for future single-cell mitochondrial
measurements in the context of cellular noise studies.


  Mitochondrial DNA (mtDNA) mutations cause severe congenital diseases but may
also be associated with healthy aging. MtDNA is stochastically replicated and
degraded, and exists within organelles which undergo dynamic fusion and
fission. The role of the resulting mitochondrial networks in the time evolution
of the cellular proportion of mutated mtDNA molecules (heteroplasmy), and
cell-to-cell variability in heteroplasmy (heteroplasmy variance), remains
incompletely understood. Heteroplasmy variance is particularly important since
it modulates the number of pathological cells in a tissue. Here, we provide the
first wide-reaching theoretical framework which bridges mitochondrial network
and genetic states. We show that, under a range of conditions, the (genetic)
rate of increase in heteroplasmy variance and de novo mutation are
proportionally modulated by the (physical) fraction of unfused mitochondria,
independently of the absolute fission-fusion rate. In the context of selective
fusion, we show that intermediate fusion/fission ratios are optimal for the
clearance of mtDNA mutants. Our findings imply that modulating network state,
mitophagy rate and copy number to slow down heteroplasmy dynamics when mean
heteroplasmy is low could have therapeutic advantages for mitochondrial disease
and healthy aging.


  How individuals adapt their behavior in cultural evolution remains elusive.
Theoretical studies have shown that the update rules chosen to model individual
decision making can dramatically modify the evolutionary outcome of the
population as a whole. This hints at the complexities of considering the
personality of individuals in a population, where each one uses its own rule.
Here, we investigate whether and how heterogeneity in the rules of behavior
update alters the evolutionary outcome. We assume that individuals update
behaviors by aspiration-based self-evaluation and they do so in their own ways.
Under weak selection, we analytically reveal a simple property that holds for
any two-strategy multi-player games in well-mixed populations and on regular
graphs: the evolutionary outcome in a population with heterogeneous update
rules is the weighted average of the outcomes in the corresponding homogeneous
populations, and the associated weights are the frequencies of each update rule
in the heterogeneous population. Beyond weak selection, we show that this
property holds for public goods games. Our finding implies that heterogeneous
aspiration dynamics is additive. This additivity greatly reduces the complexity
induced by the underlying individual heterogeneity. Our work thus provides an
efficient method to calculate evolutionary outcomes under heterogeneous update
rules.


  Background: In the marine environment, where there are few absolute physical
barriers, contemporary contact between previously isolated species can occur
across great distances, and in some cases, may be inter-oceanic. [..] in the
minke whale species complex [...] migrations [..] have been documented and
fertile hybrids and back-crossed individuals between both species have also
been identified. However, it is not known whether this represents a
contemporary event, potentially driven by ecosystem changes in the Antarctic,
or a sporadic occurrence happening over an evolutionary time-scale. We
successfully used whole genome resequencing to identify a panel of diagnostic
SNPs which now enable us address this evolutionary question.
  Results: A large number of SNPs displaying fixed or nearly fixed allele
frequency differences among the minke whale species were identified from the
sequence data. Five panels of putatively diagnostic markers were established on
a genotyping platform for validation of allele frequencies; two panels (26 and
24 SNPs) separating the two species of minke whale, and three panels (22, 23,
and 24 SNPs) differentiating the three subspecies of common minke whale. The
panels were validated against a set of reference samples, demonstrating the
ability to accurately identify back-crossed whales up to three generations.
  Conclusions: This work has resulted in the development of a panel of novel
diagnostic genetic markers to address inter-oceanic and global contact among
the genetically isolated minke whale species and sub-species. These markers,
including a globally relevant genetic reference data set for this species
complex, are now openly available for researchers [..]. The approach used here,
combining whole genome resequencing and high-throughput genotyping, represents
a universal approach to develop similar tools for other species and population
complexes.


  We devise an approach for targeted molecular design, a problem of interest in
computational drug discovery: given a target protein site, we wish to generate
a chemical with both high binding affinity to the target and satisfactory
pharmacological properties. This problem is made difficult by the enormity and
discreteness of the space of potential therapeutics, as well as the
graph-structured nature of biomolecular surface sites. Using a dataset of
protein-ligand complexes, we surmount these issues by extracting a signature of
the target site with a graph convolutional network and by encoding the discrete
chemical into a continuous latent vector space. The latter embedding permits
gradient-based optimization in molecular space, which we perform using learned
differentiable models of binding affinity and other pharmacological properties.
We show that our approach is able to efficiently optimize these multiple
objectives and discover new molecules with potentially useful binding
properties, validated via docking methods.


  Proteins must fold quickly to acquire their biologically functional
three-dimensional native structures. Hence, these are mainly stabilized by
local contacts, while intricate topologies such as knots are rare. Here, we
reveal the existence of specific patterns adopted by protein sequences and
structures to deal with backbone self-entanglement. A large scale analysis of
the Protein Data Bank shows that loops significantly intertwined with another
chain portion are typically closed by weakly bound amino acids. Why is this
energetic frustration maintained? A possible picture is that entangled loops
are formed only toward the end of the folding process to avoid kinetic traps.
Consistently, these loops are more frequently found to be wrapped around a
portion of the chain on their N-terminal side, the one translated earlier at
the ribosome. Finally, these motifs are less abundant in natural native states
than in simulated protein-like structures, yet they appear in 32% of proteins,
which in some cases display an amazingly complex intertwining.


  The two most fundamental processes describing change in biology, development
and evolu-tion, occur over drastically different timescales, difficult to
reconcile within a unified framework. Development involves temporal sequences
of cell states controlled by hierarchies of regulatory structures. It occurs
over the lifetime of a single individual, and is associated to the gene
expression level change of a given genotype. Evolution, by contrast entails
genotypic change through the acquisition/loss of genes and changes in the
network topology of interactions among genes. It involves the emergence of new,
environmentally selected phenotypes over the lifetimes of many individuals.
Here we present a model of regulatory network evolution that accounts for both
timescales. We extend the framework of Boolean models of gene regulatory
networks (GRN)-currently only applicable to describing development to include
evolutionary processes. As opposed to one-to-one maps to specific attractors,
we identify the phenotypes of the cells as the relevant macrostates of the GRN.
A phenotype may now correspond to multiple attractors, and its formal
definition no longer requires a fixed size for the genotype. This opens the
possibility for a quantitative study of the phenotypic change of a genotype,
which is itself changing over evolutionary timescales. We show how the
realization of specific phenotypes can be controlled by gene duplication events
(used here as an archetypal evolutionary event able to change the genotype),
and how successive events of gene duplication lead to new regulatory structures
via selection. At the same time, we show that our generalized framework does
not inhibit network controllability and the possibility for network control
theory to describe epigenetic signaling during development.


  Motivation: P values derived from the null hypothesis significance testing
framework are strongly affected by sample size, and are known to be
irreproducible in underpowered studies, yet no suitable replacement has been
proposed. Results: Here we present implementations of non-parametric
standardized median effect size estimates, dNEF, for high-throughput sequencing
datasets. Case studies are shown for transcriptome and tag-sequencing datasets.
The dNEF measure is shown to be more reproducible and robust than P values and
requires sample sizes as small as 3 to reproducibly identify differentially
abundant features. Availability: Source code and binaries freely available at:
https://bioconductor.org/packages/ALDEx2.html , omicplotR, and
https://github.com/ggloor/CoDaSeq .


  The aim of this work is to investigate the optimal vaccine sharing between
two SIR centers in the presence of migration fluxes of susceptibles and
infected individuals during the mumps outbreak. Optimality of the vaccine
allocation means the minimization of the total number of lost working days
during the whole period of epidemic outbreak $[0,t_f]$, which can be described
by the functional $Q=\int_0^{t_f}I(t){\rm d}t$ where $I(t)$ stands for the
number of infectives at time $t$. We explain the behavior of the optimal
allocation, which depends on the model parameters and the amount of available
vaccine.


  This work investigates the emergence of oscillations in one of the simplest
cellular signaling networks exhibiting oscillations, namely, the dual-site
phosphorylation and dephosphorylation network (futile cycle), in which the
mechanism for phosphorylation is processive while the one for dephosphorylation
is distributive (or vice-versa). The fact that this network yields oscillations
was shown recently by Suwanmajo and Krishnan. Our results, which significantly
extend their analyses, are as follows. First, in the three-dimensional space of
total amounts, the border between systems with a stable versus unstable steady
state is a surface defined by the vanishing of a single Hurwitz determinant.
Second, this surface consists generically of simple Hopf bifurcations. Next,
simulations suggest that when the steady state is unstable, oscillations are
the norm. Finally, the emergence of oscillations via a Hopf bifurcation is
enabled by the catalytic and association constants of the distributive part of
the mechanism: if these rate constants satisfy two inequalities, then the
system generically admits a Hopf bifurcation. Our proofs are enabled by the
Routh-Hurwitz criterion, a Hopf-bifurcation criterion due to Yang, and a
monomial parametrization of steady states.


  This is an epidemiological SIRV model based study that is designed to analyze
the impact of vaccination in containing infection spread, in a 4-tiered
population compartment comprised of susceptible, infected, recovered and
vaccinated agents. While many models assume a lifelong protection through
vaccination, we focus on the impact of waning immunization due to conversion of
vaccinated and recovered agents back to susceptible ones. Two asymptotic states
exist, the "disease-free equilibrium" and the "endemic equilibrium"; we express
the transitions between these states as function of the vaccination and
conversion rates using the basic reproduction number as a descriptor. We find
that the vaccination of newborns and adults have different consequences in
controlling epidemics. We also find that a decaying disease protection within
the recovered sub-population is not sufficient to trigger an epidemic at the
linear level. Our simulations focus on parameter sets that could model a
disease with waning immunization like pertussis. For a diffusively coupled
population, a transition to the endemic state can be initiated via the
propagation of a traveling infection wave, described successfully within a
Fisher-Kolmogorov framework.


  The aim of this review is to present and analyze the probabilistic models of
mathematical phylogenetics which have been intensively used in recent years in
biology as the cornerstone of attempts to infer and reconstruct the ancestral
relationships between species. We outline the development of theoretical
phylogenetics, from the earliest studies based on morphological characters,
through to the use of molecular data in a wide variety of forms. We bring the
lens of mathematical physics to bear on the formulation of theoretical models,
focussing on the applicability of many methods from the toolkit of that
tradition -- techniques of groups and representations to guide model
specification and to exploit the multilinear setting of the models in the
presence of underlying symmetries; extensions to coalgebraic properties of the
generators associated to rate matrices underlying the models, in relation to
the graphical structures (trees and networks) which form the search space for
inferring evolutionary trees. Aspects presented, include relating model classes
to relevant matrix Lie algebras, as well as manipulations with group characters
to enumerate various natural polynomial invariants, for identifying robust,
low-parameter quantities for use in inference. Above all, we wish to emphasize
the many features of multipartite entanglement which are shared between
descriptions of quantum states on the physics side, and the multi-way tensor
probability arrays arising in phylogenetics. In some instances, well-known
objects such as the Cayley hyperdeterminant (the `tangle') can be directly
imported into the formalism -- for models with binary character traits, and
triplets of taxa. In other cases new objects appear, such as the remarkable
quintic `squangle' invariants for quartet tree discrimination and DNA data,
with their own unique interpretation in the phylogenetic modeling context.


  Microtubules are biological protein polymers with critical and diverse
functions. Their structures share some similarities with photosynthetic antenna
complexes, particularly in the ordered arrangement of photoactive molecules
with large transition dipole moments. As the role of photoexcitations in
microtubules remains an open question, here we analyze tryptophan molecules,
the amino acid building block of microtubules with the largest transition
dipole strength. By taking their positions and dipole orientations from
realistic models capable of reproducing tubulin experimental spectra, and using
a Hamiltonian widely employed in quantum optics to describe light-matter
interactions, we show that such molecules arranged in their native microtubule
configuration exhibit a superradiant ground state, which represents an
excitation fully extended on the chromophore lattice. We also show that such a
superradiant ground state emerges due to supertransfer coupling between the
ground states of smaller blocks of the microtubule. In the dynamics we find
that the spreading of excitation is ballistic in the absence of external
sources of disorder and strongly dependent on initial conditions. The velocity
of photoexcitation spreading is shown to be enhanced by the supertransfer
effect with respect to the velocity one would expect from the strength of the
nearest-neighbor coupling between tryptophan molecules in the microtubule.
Finally, such structures are shown to have an enhanced robustness to static
disorder when compared to geometries that include only short-range
interactions. These cooperative effects (superradiance and supertransfer) may
induce ultra-efficient photoexcitation absorption and could enhance excitonic
energy transfer in microtubules over long distances under physiological
conditions.


  Liquid cellular compartments spatially segregate from the cytoplasm and can
regulate aberrant protein aggregation, a process linked to several medical
conditions, including Alzheimer's and Parkinson's diseases. Yet the mechanisms
by which these droplet-like compartments affect protein aggregation remain
unknown. Here, we combine kinetic theory of protein aggregation and
liquid-liquid phase separation to study the spatial control of irreversible
protein aggregation in the presence of liquid compartments. We find that, even
for weak interactions between the compartment constituents and the aggregating
monomers, aggregates are strongly enriched inside the liquid compartment
relative to the surrounding cytoplasm. We show that this enrichment is caused
by a positive feedback mechanism of aggregate nucleation and growth which is
mediated by a flux maintaining the phase equilibrium between the compartment
and the cytoplasm. Our model predicts that the compartment volume that
maximizes aggregate enrichment in the compartment is determined by the reaction
orders of aggregate nucleation. The underlying mechanism of aggregate
enrichment could be used to confine cytotoxic protein aggregates inside
droplet-like compartments suggesting potential new avenues against aberrant
protein aggregation. Our findings could also represent a common mechanism for
the spatial control of irreversible chemical reactions in general.


  We investigate a time-dependent spatial vector-host epidemic model with
non-coincident domains for the vector and host populations. The host population
resides in small non-overlapping sub-regions, while the vector population
resides throughout a much larger region. The dynamics of the populations are
modeled by a reaction-diffusion-advection compartmental system of partial
differential equations. The disease is transmitted through vector and host
populations in criss-cross fashion. We establish global well-posedness and
uniform a prior bounds as well as the long-term behavior. The model is applied
to simulate the outbreak of bluetongue disease in sheep transmitted by midges
infected with bluetongue virus. We show that the long-range directed movement
of the midge population, due to wind-aided movement, enhances the transmission
of the disease to sheep in distant sites.


  Effects of changes in functional composition of soil communities on nutrient
cycling are still not well understood. Models simulating community dynamics
overcome the technical challenges of conducting species removal experiments in
the field. However, to date, available soil food web models do not adequately
represent the organic matter processing chain which is key for soil dynamics.
Here, we present a new model of soil food web dynamics accounting for
allometric scaling of metabolic rate, ontogeny of organic matter, and explicit
representation of nitrogen and carbon flows. We use this model to investigate
what traits are best predictors of species effects on community productivity
and on nutrient cycling. To do so, we removed 161 tropho-species (groups of
functionally identical species) one at a time from 48 forest soil food webs,
and simulated their dynamics until equilibrium. Simulations revealed that
combinations of traits better determine removal effects than single ones. The
smallest species are the most competitive ones, but carnivores of various body
masses presenting the highest connectivity and resource similarity could be key
stone species in the regulation of competitive forces. Despite this, most
removals had low effects, suggesting functional redundancy provides a high
resistance of soil food webs to single tropho-species extinction. We also
highlight for the first time that food web structure and soil fertility can
drastically change species effects in an unpredictable way. Moreover, the
exclusion of detritus and stoichiometric constraints in past studies lead to
underestimations of indirect effects and retroactions. While additional work is
needed to incorporate complementarity between detritivores, it is essential to
take into account these mechanisms in models in order to improve the
understanding of soil food web functioning.


  Since 1990 Jeanne Louise Calment has held the record for human longevity. She
was born on 21 February 1875, became the longest-lived human on 12 May 1990
when she was 115.21 and died on 4 August 1997 at age 122.45 years. In this
chapter, we use data available on 25 September 2017 on people who reached age
110, supercentenarians, to address the following questions: 1. How likely is it
that a person has reached age 122.45? 2. How unlikely is it that Calment's
record has not yet been broken? 3. How soon might it be broken? Assuming a
constant annual probability of death of 50% after age 110, we found that the
probability that a person who survived to age 110 would have lived to 122.45 by
25 September 2017 is 17.1%. Furthermore, we calculated that there was only a
20.3% chance that Calment's record would have been broken after 1997 but before
2017. Finally, we estimated that there is less than a 50% chance that someone
will surpass Calment's lifespan before 2045. Jeanne Louise Calment's record is
exceptional but not impossible. It does not provide evidence that her lifespan
is an upper limit to human lifespans.


  The goal of this dissertation is to study the sequence polymorphism in
retrotransposable elements of Entamoeba histolytica. The Quasispecies theory, a
concept of equilibrium (stationary), has been used to understand the behaviour
of these elements. Two datasets of retrotransposons of Entamoeba histolytica
have been used. We present results from both datasets of retrotransposons
(SINE1s) of E. histolytica. We have calculated the mutation rate of EhSINE1s
for both datasets and drawn a phylogenetic tree for newly determined EhSINE1s
(dataset II). We have also discussed the variation in lengths of EhSINE1s for
both datasets. Using the quasispecies model, we have shown how sequences of
SINE1s vary within the population. The outputs of the quasispecies model are
discussed in the presence and the absence of back mutation by taking different
values of fitness. From our study of Non-long terminal repeat retrotransposons
(LINEs and their non-autonomous partner's SINEs) of Entamoeba histolytica, we
can conclude that an active EhSINE can generate very similar copies of itself
by retrotransposition. Due to this reason, it increases mutations which give
the result of sequence polymorphism. We have concluded that the mutation rate
of SINE is very high. This high mutation rate provides an idea for the
existence of SINEs, which may affect the genetic analysis of EhSINE1
ancestries, and calculation of phylogenetic distances.


  In nature and society problems arise when different interests are difficult
to reconcile, which are modeled in game theory. While most applications assume
uncorrelated games, a more detailed modeling is necessary to consider the
correlations that influence the decisions of the players. The current theory
for correlated games, however, enforces the players to obey the instructions
from a third party or "correlation device" to reach equilibrium, but this
cannot be achieved for all initial correlations. We extend here the existing
framework of correlated games and find that there are other interesting and
previously unknown Nash equilibria that make use of correlations to obtain the
best payoff. This is achieved by allowing the players the freedom to follow or
not to follow the suggestions of the correlation device. By assigning
independent probabilities to follow every possible suggestion, the players
engage in a response game that turns out to have a rich structure of Nash
equilibria that goes beyond the correlated equilibrium and mixed-strategy
solutions. We determine the Nash equilibria for all possible correlated
Snowdrift games, which we find to be describable by Ising Models in thermal
equilibrium. We believe that our approach paves the way to a study of
correlations in games that uncovers the existence of interesting underlying
interaction mechanisms, without compromising the independence of the players.


  This paper investigates the stochastic permanence of malaria and the
existence of a stationary distribution for the stochastic process describing
the disease dynamics over sufficiently longtime. The malaria system is highly
random with fluctuations from the disease transmission and natural deathrates,
which are expressed as independent white noise processes in a family of
stochastic differential equation epidemic models. Other sources of variability
in the malaria dynamics are the random incubation and naturally acquired
immunity periods of malaria. Improved analytical techniques and local
martingale characterizations are applied to describe the character of the
sample paths of the solution process of the system in the neighborhood of an
endemic equilibrium. Emphasis of this study is laid on examination of the
impacts of (1) the sources of variability- disease transmission and natural
death rates, and (2) the intensities of the white noise processes in the system
on the stochastic permanence of malaria, and also on the existence of the
stationary distribution for the solution process over sufficiently long time.
Numerical simulation examples are presented to illuminate the persistence and
stochastic permanence of malaria, and also to numerically approximate the
stationary distribution of the states of the solution process.


  A comparative stochastic and deterministic study of a family of SEIRS
epidemic dynamic models for malaria is presented. The family type is determined
by the qualitative behavior of the nonlinear incidence rates of the disease.
Furthermore, the malaria models exhibit three random delays:- two of the delays
represent the incubation periods of the disease inside the vector and human
hosts, whereas the third delay is the period of effective natural immunity
against the disease. The stochastic malaria models are improved by including
the random environmental fluctuations in the disease transmission and natural
death rates of humans. Insights about the effects of the delays and the noises
on the malaria dynamics are gained via comparative analyses of the family of
stochastic and deterministic models, and further critical examination of the
significance of the intensities of the white noises in the system on (1) the
existence and stability of the equilibria, and also on (2) the eradication and
persistence of malaria in the human population. The basic reproduction numbers
and other threshold values for malaria in the stochastic and deterministic
settings are determined and compared for the cases of constant or random delays
in the system. Numerical simulation results are presented.


  The ICTV develops, refines and maintains a universal virus taxonomy; Order is
the highest taxon in the branching hierarchy of recognised viral taxa.
Historically, ICTV (sub)committees have classified viruses on the basis of
morphological characteristics and various other attributes. Today, virtually
all new viral genomes are assembled from metagenomic datasets and are not
linked directly to biological agents. Thus, placing a virus into a taxonomic
scheme solely from primary genome structure is an increasingly important
problem. Various simple descriptive statistics of a viral genome sequence have
been used successfully for virus classification. Here, we use the NCBI's viral
and viroid reference sequence collection (RefSeq) and a common experimental
framework to compare the performance of different genome sequence-derived
features and classifiers in the task of assigning a virus to one of seven ICTV
Orders. The nucleotide-, word-, and compression-based features we consider
include genome length, the k-mer Natural Vector (k = 1, ..., 6) and its
derivatives, return time distribution, and general-purpose and DNA-specific
compression ratios; the classifiers used are the k-NN and SVM. The combination
of genome length and k-NN has the worst, yet still respectable, performance
(mean error rate of 0.137); the best performance is achieved using 4-mer counts
and SVM (mean error rate of 0.006). We investigate the main causes of
misclassification, explore which viruses are more difficult to classify, and
use the best performing combination to predict the Orders of 1,834 unclassified
viruses. A subsequent version of RefSeq assigned Orders to 17 of these
previously unlabelled viruses. Since 16 of our predictions match these
assignments, our approach could aid virologists dealing with viruses that are
known only from sequence data.


  Clusters of IP3 receptor channels in the membranes of the endoplasmic
reticulum (ER) of many non-excitable cells release calcium ions in a
cooperative manner giving rise to dynamical patterns such as Ca2+ puffs, waves,
and oscillations that occur on multiple spatial and temporal scales. We
introduce a minimal yet descriptive reaction-diffusion model of IP3 receptors
for a saturating concentration of IP3 using a principled reduction of a
detailed Markov chain description of individual channels. A dynamical systems
analysis reveals the possibility of excitable, bistable and oscillatory
dynamics of this model that correspond to three types of observed patterns of
calcium release -- puffs, waves, and oscillations respectively. We explain the
emergence of these patterns via a bifurcation analysis of a coupled two-cluster
model, compute the phase diagram and quantify the speed of the waves and period
of oscillations in terms of system parameters. We connect the termination of
large-scale Ca2+ release events to IP3 unbinding or stochasticity.


  We study the dynamics of predator-prey systems where prey are confined to a
single region of space and where predators move randomly according to a
power-law (L\'evy) dispersal kernel. Site fidelity, an important feature of
animal behaviour, is incorporated in the model through a stochastic resetting
dynamics of the predators to the prey patch. We solve in the long time limit
the rate equations of Lotka-Volterra type that describe the evolution of the
two species densities. Fixing the demographic parameters and the L\'evy
exponent, the total population of predators can be maximized for a certain
value of the resetting rate. This optimal value achieves a compromise between
over-exploitation and under-utilization of the habitat. Similarly, at fixed
resetting rate, there exists a L\'evy exponent which is optimal regarding
predator abundance. These findings are supported by 2D stochastic simulations
and show that the combined effects of diffusion and resetting can broadly
extend the region of species coexistence in ecosystems facing resources
scarcity.


  The growth of a population divided among spatial sites, with migration
between the sites, is sometimes modelled by a product of random matrices, with
each diagonal elements representing the growth rate in a given time period, and
off-diagonal elements the migration rate. The randomness of the matrices then
represents stochasticity of environmental conditions. We consider the case
where the off-diagonal elements are small, representing a situation where
migration has been introduced into an otherwise sessile meta-population. We
examine the asymptotic behaviour of the long-term growth rate. When there is a
single site with the highest growth rate, under the assumption of Gaussian log
growth rates at the individual sites (or having Gaussian-like tails) we show
that the behavior near zero is like a power of $\epsilon$, and derive upper and
lower bounds for the power in terms of the difference in the growth rates and
the distance between the sites. In particular, when the difference in mean log
growth rate between two sites is sufficiently small, or the variance of the
difference between the sites sufficiently large, migration will always be
favored by natural selection, in the sense that introducing a small amount of
migration will increase the growth rate of the population relative to the
zero-migration case.


  We follow up on a companion work that considered growth rates of populations
growing at different sites, with different randomly varying growth rates at
each site, in the limit as migration between sites goes to 0. We extend this
work here to the special case where the maximum average log growth rate is
achieved at two different sites. The primary motivation is to cover the case
where `sites' are understood as age classes for the same individuals. The
theory then calculates the effect on growth rate of introducing a rare delay in
development, a diapause, into an otherwise fixed-length semelparous life
history.
  Whereas the increase in stochastic growth rate due to rare migrations was
found to grow as a power of the migration rate, we show that under quite
general conditions that in the diapause model --- or in the migration model
with two or more sites having equal individual stochastic growth rates --- the
increase in stochastic growth rate due to diapause at rate $\epsilon$ behaves
like $(\log \epsilon^{-1})^{-1}$ as $\epsilon\downarrow 0$. In particular, this
implies that a small random disruption to the deterministic life history will
always be favored by natural selection, in the sense that it will increase the
stochastic growth rate relative to the zero-delay deterministic life history.


  Quadratic programming (QP) is a common and important constrained optimization
problem. Here, we derive a surprising duality between constrained optimization
with inequality constraints -- of which QP is a special case -- and consumer
resource models describing ecological dynamics. Combining this duality with a
recent `cavity solution', we analyze high-dimensional, random QP where the
optimization function and constraints are drawn randomly. Our theory shows
remarkable agreement with numerics and points to a deep connection between
optimization, dynamical systems, and ecology.


  What happens when a new social convention replaces an old one? While the
possible forces favoring norm change - such as institutions or committed
activists - have been identified since a long time, little is known about how a
population adopts a new convention, due to the difficulties of finding
representative data. Here we address this issue by looking at changes occurred
to 2,541 orthographic and lexical norms in English and Spanish through the
analysis of a large corpora of books published between the years 1800 and 2008.
We detect three markedly distinct patterns in the data, depending on whether
the behavioral change results from the action of a formal institution, an
informal authority or a spontaneous process of unregulated evolution. We
propose a simple evolutionary model able to capture all the observed behaviors
and we show that it reproduces quantitatively the empirical data. This work
identifies general mechanisms of norm change and we anticipate that it will be
of interest to researchers investigating the cultural evolution of language
and, more broadly, human collective behavior.


  Debye-Waller factor, a measure of X-ray attenuation, can be experimentally
observed in protein X-ray crystallography. Previous theoretical models have
made strong inroads in the analysis of B-factors by linearly fitting protein
B-factors from experimental data. However, the blind prediction of B-factors
for unknown proteins is an unsolved problem. This work integrates machine
learning and advanced graph theory, namely, multiscale weighted colored graphs
(MWCGs), to blindly predict B-factors of unknown proteins. MWCGs are local
features that measure the intrinsic flexibility due to a protein structure.
Global features that connect the B-factors of different proteins, e.g., the
resolution of X-ray crystallography, are introduced to enable the cross-protein
B-factor predictions. Several machine learning approaches, including ensemble
methods and deep learning, are considered in the present work. The proposed
method is validated with hundreds of thousands of experimental B-factors.
Extensive numerical results indicate that the blind B-factor predictions
obtained from the present method are more accurate than the least squares
fittings using traditional methods.


  Pairwise models are used widely to model epidemic spread on networks. These
include the modelling of susceptible-infected-removed (SIR) epidemics on
regular networks and extensions to SIS dynamics and contact tracing on more
exotic networks exhibiting degree heterogeneity, directed and/or weighted links
and clustering. However, extra features of the disease dynamics or of the
network lead to an increase in system size and analytical tractability becomes
problematic. Various `closures' can be used to keep the system tractable.
Focusing on SIR epidemics on regular but clustered networks, we show that even
for the most complex closure we can determine the epidemic threshold as an
asymptotic expansion in terms of the clustering coefficient.We do this by
exploiting the presence of a system of fast variables, specified by the
correlation structure of the epidemic, whose steady state determines the
epidemic threshold. While we do not find the steady state analytically, we
create an elegant asymptotic expansion of it. We validate this new threshold by
comparing it to the numerical solution of the full system and find excellent
agreement over a wide range of values of the clustering coefficient,
transmission rate and average degree of the network. The technique carries over
to pairwise models with other closures [1] and we note that the epidemic
threshold will be model dependent. This emphasises the importance of model
choice when dealing with realistic outbreaks.


  We show that the expected size of the maximum agreement subtree of two
$n$-leaf trees, uniformly random among all trees with the shape, is
$\Theta(\sqrt{n})$. To derive the lower bound, we prove a global structural
result on a decomposition of rooted binary trees into subgroups of leaves
called blobs. To obtain the upper bound, we generalize a first moment argument
for random tree distributions that are exchangeable and not necessarily
sampling consistent.


  Large parallel gains in the development of both computational resources as
well as sampling methods have now made it possible to simulate dissociation
events in ligand-protein complexes with all--atom resolution. Such encouraging
progress, together with the inherent spatiotemporal resolution associated with
molecular simulations, has left their use for investigating dissociation
processes brimming with potential, both in rational drug design, where it can
be an invaluable tool for determining the mechanistic driving forces behind
dissociation rate constants, as well as in force-field development, where it
can provide a catalog of transient molecular structures on which to refine
force-fields. Although much progress has been made in making force-fields more
accurate, reducing their error for transient structures along a transition path
could yet prove to be a critical development helping to make kinetic
predictions much more accurate. In what follows we will provide a
state-of-the-art compilation of the molecular dynamics (MD) methods used to
investigate the kinetics and mechanisms of ligand-protein dissociation
processes. Due to the timescales of such processes being slower than what is
accessible using straightforward MD simulations, several ingenious schemes are
being devised at a rapid rate to overcome this obstacle. Here we provide an
up-to-date compendium of such methods and their achievements/shortcomings in
extracting mechanistic insight into ligand-protein dissociation. We conclude
with a critical and provocative appraisal attempting to answer the title of
this review.


  Reaction networks are mathematical models of interacting chemical species
that are primarily used in biochemistry. There are two modeling regimes that
are typically used, one of which is deterministic and one that is stochastic.
In particular, the deterministic model consists of an autonomous system of
differential equations, whereas the stochastic system is a continuous time
Markov chain. Connections between the two modeling regimes have been studied
since the seminal paper by Kurtz (1972), where the deterministic model is shown
to be a limit of a properly rescaled stochastic model over compact time
intervals. Further, more recent studies have connected the long-term behaviors
of the two models when the reaction network satisfies certain graphical
properties, such as weak reversibility and a deficiency of zero.
  These connections have led some to conjecture a link between the long-term
behavior of the two models exists, in some sense. In particular, one is tempted
to believe that positive recurrence of all states for the stochastic model
implies the existence of positive equilibria in the deterministic setting, and
that boundary equilibria of the deterministic model imply the occurrence of an
extinction event in the stochastic setting. We prove in this paper that these
implications do not hold in general, even if restricting the analysis to
networks that are bimolecular and that conserve the total mass. In particular,
we disprove the implications in the special case of models that have absolute
concentration robustness, thus answering in the negative a conjecture stated in
the literature in 2014.


  Many plants harbor complex mechanisms that promote outcrossing and efficient
pollen transfer. These include floral adaptations as well as genetic
mechanisms, such as molecular self-incompatibility (SI) systems. The
maintenance of such systems over long evolutionary timescales suggests that
outcrossing is favorable over a broad range of conditions. Conversely, SI has
repeatedly been lost, often in association with transitions to
self-fertilization (selfing). This transition is favored when the short-term
advantages of selfing outweigh the costs, primarily inbreeding depression. The
transition to selfing is expected to have major effects on population genetic
variation and adaptive potential, as well as on genome evolution. In the
Brassicaceae, many studies on the population genetic, gene regulatory and
genomic effects of selfing have centered on the model plant Arabidopsis
thaliana and the crucifer genus Capsella. The accumulation of population
genomics datasets have allowed detailed investigation of where, when and how
the transition to selfing occurred. Future studies will take advantage of the
development of population genetics theory on the impact of selfing, especially
regarding positive selection. Furthermore, investigation of systems including
recent transitions to selfing, mixed mating populations and/or multiple
independent replicates of the same transition will facilitate dissecting the
effects of mating system variation from processes driven by demography.


  Transcription is the first step of gene expression, in which a particular
segment of DNA is copied to RNA by the enzyme RNA polymerase (RNAP). Despite
many details of the complex interactions between DNA and RNA synthesis
disclosed experimentally, much of physical behavior of transcription remains
largely unknown. Understanding torsional mechanics of DNA and RNAP together
with its nascent RNA and RNA-bound proteins in transcription maybe the first
step towards deciphering the mechanism of gene expression. In this study, based
on the balance between viscous drag on RNA synthesis and torque resulted from
untranscribed supercoiled DNA template, a simple model is presented to describe
mechanical properties of transcription. With this model, the rotation and
supercoiling density of the untranscribed DNA template are discussed in detail.
Two particular cases of transcription are considered, transcription with
constant velocity and transcription with torque dependent velocity. Our results
show that, during the initial stage of transcription, rotation originated from
the transcribed part of DNA template is mainly released by the rotation of RNAP
synthesis. During the intermediate stage, the rotation is usually released by
both the supercoiling of the untranscribed part of DNA template and the
rotation of RNAP synthesis, with proportion depending on the friction
coefficient in environment and the length of nascent RNA. However, with the
approaching to the upper limit of twisting of the untranscribed DNA template,
the rotation resulted from transcription will then be mainly released by the
rotation of RNAP synthesis.


  Introduction: To leverage functionality and clinical relevance into
understanding systems biology, one needs to understand the pathway of the
genetic effects on risk factors/disease through intermediate molecular levels,
such as metabolomics. Systems approaches integrate multi-omic information to
find pathways to disease endpoints and make optimal inference decisions.
Method: Here, we introduce a multi-stage approach to integrate causal networks
in observational studies and GWAS to facilitate mechanistic understanding
through identification of pathways from the genome to risk factors/disease via
metabolomics. The pathways in causal networks reveal the underlying
relationships behind observations, which do not play a significant role in more
traditional correlative analyses, where one variable at a time is considered.
Results: We identified a causal network over the metabolomic level using the
genome directed acyclic graph (G-DAG), to systematically assess whether
variations in the genome lead to variations in triglyceride levels as a risk
factor of cardiovascular disease. We found LRRC46 and LRRC69 harboring
loss-of-function mutations have significant effect on two metabolites with
direct effects on triglyceride levels. We also found pathways of FAM198B and
C6orf25 to triglycerides through indirect paths from metabolites. Conclusion:
Integrating causal networks with GWAS facilitates mechanistic understanding in
comparison to one-variable-at-a-time approaches due to accounting for
relationships among components at intermediate molecular levels. This approach
is complementary to experimental studies to identify efficacious targets in the
age of big data sets.


  This study assessed the effects of migration on the dynamics of a species
population. It was considered that the species in its natural state and without
the presence of migration exhibited Allee effect. This work also considered
migration as a density-dependent function, which, from a maximum rate,
decreases to a minimum of zero when the population reaches its carrying
capacity.


  Learning from a partner who collects higher payoff is a frequently used
working hypothesis in evolutionary game theory. One of the alternative
dynamical rules is when the focal player prefers to follow the strategy choice
of the majority in the local neighborhood, which is often called as
conformity-driven strategy update. In this work we assume that both strategy
learning methods are present and compete for space within the framework of a
coevolutionary model. Our results reveal that the presence of payoff-driven
strategy learning method becomes exclusive for high succer's payoff and/or high
temptation values that represent a snowdrift game dilemma situation. In
general, however, the competition of the mentioned strategy learning methods
could be useful to enlarge the parameter space where only cooperators prevail.
This success of cooperation is based on the enforced coordination of cooperator
players which reveals the benefit of the latter strategy. Interestingly, the
payoff-based and the conformity-based cooperator players can form an effective
alliance against defectors that can also extend the parameter space of full
cooperator solution in the stag-hunt game region. Our work highlights that the
coevolution of strategies and individual features such as learning method can
provide novel type of pattern formation mechanism that cannot be observed in a
static model, hence remains hidden in traditional models.


  Repeated exposure to low-level blast may initiate a range of adverse health
problem such as traumatic brain injury (TBI). Although many studies
successfully identified genes associated with TBI, yet the cellular mechanisms
underpinning TBI are not fully elucidated. In this study, we investigated
underlying relationship among genes through constructing transcript Bayesian
networks using RNA-seq data. The data for pre- and post-blast transcripts,
which were collected on 33 individuals in Army training program, combined with
our system approach provide unique opportunity to investigate the effect of
blast-wave exposure on gene-gene interactions. Digging into the networks, we
identified four subnetworks related to immune system and inflammatory process
that are disrupted due to the exposure. Among genes with relatively high fold
change in their transcript expression level, ATP6V1G1, B2M, BCL2A1, PELI,
S100A8, TRIM58 and ZNF654 showed major impact on the dysregulation of the
gene-gene interactions. This study reveals how repeated exposures to traumatic
conditions increase the level of fold change of transcript expression and
hypothesizes new targets for further experimental studies.


  Stochastic simulation can make the molecular processes of cellular control
more vivid than the traditional differential-equation approach by generating
typical system histories instead of just statistical measures such as the mean
and variance of a population. Simple simulations are now easy for students to
construct from scratch, that is, without recourse to black-box packages. In
some cases, their results can also be compared directly to single-molecule
experimental data. After introducing the stochastic simulation algorithm, this
article gives two case studies, involving gene expression and error correction,
respectively. Code samples and resulting animations showing results are given
in the online supplements.


  Invasive species have been recognized as a leading threat to biodiversity. In
particular, lakes are especially affected by species invasions because they are
closed systems sensitive to disruption. Accurately controlling the spread of
invasive species requires solving a complex spatial-dynamic optimization
problem. In this work we propose a novel framework for determining the optimal
management strategy to maximize the value of a lake system net of damages from
invasive species, including an endogenous diffusion mechanism for the spread of
invasive species through boaters' trips between lakes. The proposed method
includes a combined global iterative process which determines the optimal
number of trips to each lake in each season and the spatial-dynamic optimal
boat ramp fee.


  Resonance energy transfer has become an indispensable experimental tool for
single-molecule and single-cell biophysics. Its physical underpinnings,
however, are subtle: It involves a discrete jump of excitation from one
molecule to another, and so we regard it as a strongly quantum-mechanical
process. And yet, its kinetics differ from what many of us were taught about
two-state quantum systems; quantum superpositions of the states do not seem to
arise; and so on. Although J. R. Oppenheimer and T. F\"orster navigated these
subtleties successfully, it remains hard to find an elementary derivation in
modern language. The key step involves acknowledging quantum decoherence.
Appreciating that aspect can be helpful when we attempt to extend our
understanding to situations where F\"orster's original analysis is not
applicable.


  The chapter reviews the syntax to store machine-readable annotations and
describes the mapping between rule-based modelling entities (e.g., agents and
rules) and these annotations. In particular, we review an annotation framework
and the associated guidelines for annotating rule-based models of molecular
interactions, encoded in the commonly used Kappa and BioNetGen languages, and
present prototypes that can be used to extract and query the annotations. An
ontology is used to annotate models and facilitate their description.


  If a cultural feature is transmitted over generations and exposed to
stochastic selection when spreading in a population, its evolution may be
governed by statistical laws and be partly predictable, as in the case of
genetic evolution. Music exhibits steady changes of styles over time, with new
characteristics developing from traditions. Recent studies have found trends in
the evolution of music styles, but little is known about their relations to the
evolution theory. Here we analyze Western classical music data and find
statistical evolutionary laws. For example, distributions of the frequencies of
some rare musical events (e.g. dissonant intervals) exhibit steady increase in
the mean and standard deviation as well as constancy of their ratio. We then
study an evolutionary model where creators learn their data-generation models
from past data and generate new data that will be socially selected by
evaluators according to the content dissimilarity (novelty) and style
conformity (typicality) with respect to the past data. The model reproduces the
observed statistical laws and can make non-trivial predictions for the
evolution of independent musical features. In addition, the same model with
different parameterization can predict the evolution of Japanese enka music,
which is developed in a different society and has a qualitatively different
tendency of evolution. Our results suggest that the evolution of musical styles
can partly be explained and predicted by the evolutionary model incorporating
statistical learning, which can be important for other cultures and future
music technologies.


  Motile organisms often use finite spatial perception of their surroundings to
navigate and search their habitats. Yet standard models of search are usually
based on purely local sensory information. To model how a finite perceptual
horizon affects ecological search, we propose a framework for optimal
navigation that combines concepts from random walks and optimal control theory.
We show that, while local strategies are optimal on asymptotically long and
short search times, finite perception yields faster convergence and increased
search efficiency over transient time scales relevant in biological systems.
The benefit of the finite horizon can be maintained by the searchers tuning
their response sensitivity to the length scale of the stimulant in the
environment, and is enhanced when the agents interact as a result of increased
consensus within subpopulations. Our framework sheds light on the role of
spatial perception and transients in search movement and collective sensing of
the environment.


  Relative fluctuations of observables in discrete stochastic systems are
bounded at all times by the mean dynamical activity in the system, quantified
by the mean number of jumps. This constitutes a kinetic uncertainty relation
that is fundamentally different from the thermodynamic uncertainty relation
recently discussed in the literature. The thermodynamic constraint is more
relevant close to equilibrium while the kinetic constraint is the limiting
factor of the precision of a observables in regimes far from equilibrium. This
is visualized for paradigmatic simple systems and with an example of molecular
motor dynamics. Our approach is based on the recent fluctuation response
inequality by Dechant and Sasa [arXiv:1804.08250] and can be applied to generic
Markov jump systems, which describe a wide class of phenomena and observables,
including the irreversible predator-prey dynamics that we use as an
illustration.