  In the study of the evolution of cooperation, many mechanisms have been
proposed to help overcome the self-interested cheating that is individually
optimal in the Prisoners' Dilemma game. These mechanisms include assortative or
networked social interactions, other-regarding preferences considering the
payoffs of others, reciprocity rules to establish cooperation as a social norm,
and multilevel selection involving simultaneous competition between individuals
favoring cheaters and competition between groups favoring cooperators. In this
paper, we build on recent work studying PDE replicator equations for multilevel
selection to understand how within-group mechanisms of assortment,
other-regarding preferences, and both direct and indirect reciprocity can help
to facilitate cooperation in concert with evolutionary competition between
groups. We consider a group-structured population in which interactions between
individuals consist of Prisoners' Dilemma games, and study the dynamics of
multilevel competition determined by the payoffs individuals receive when
interacting according to these within-group mechanisms. We find that the
presence of each of these mechanisms acts synergistically with multilevel
selection for the promotion of cooperation, decreasing the strength of
between-group competition required to sustain long-time cooperation and
increasing the collective payoff achieved by the population. However, we find
that only other-regarding preferences allows for the achievement of socially
optimal collective payoffs for Prisoners' Dilemma games in which average payoff
is maximized by an intermediate mix of cooperators and defectors. For the other
three mechanisms, the multilevel dynamics remain susceptible to a shadow of
lower-level selection, as the collective outcome fails to exceed the payoff of
the all-cooperator group.


  With the rapid growth in high-throughput biological sequencing technologies
and subsequently the amount of produced omics data, it is essential to develop
automated methods to annotate the functionality of unknown genes and proteins.
There are developed tools such as AHRD applying known proteins characterization
to annotate unknown ones. Some other algorithms such as eggNOG apply
orthologous groups of proteins to detect the most probable function. However,
while the available tools focus on the detection of the most similar
characterization, they are not able to generalize and integrate information
from multiple homologs while maintaining accuracy. Here, we devise AFDP, an
integrated approach for protein function prediction which benefits from the
combination of two available tools, AHRD and eggNOG, to predict the
functionality of novel proteins and produce more precise human readable
descriptions by applying our stCFExt algorithm. StCFExt creates function
descriptions applying available manually curated descriptions in swiss-prot.
Using a benchmark dataset we show that the annotations predicted by our
approach are more accurate than eggNOG and AHRD annotations.


  DNA and RNA quadruplexes are extensively studied for the key roles they are
suspected to play in the cellular regulation networks at both genomic and
transcriptomic levels. The reliable detection of quadruplexes in cells was and
remains a challenging task. Here, we describe the various strategies that have
been implemented over the past years to visualize functionally relevant DNA and
RNA quadruplexes in human cells, from immunodetection studies to the design and
use of quadruplex-specific turn-on fluorescent probes.


  In the modeling of complex biological systems, the use of power-law models
(such as S-systems and GMA systems) often provides a remarkable accuracy over
several orders of magnitude in concentrations, an unusually broad range not
fully understood at present. In order to provide additional insight in this
sense, this article is devoted to the renormalization group analysis of
reactions in fractal or self-similar media. In particular, the renormalization
group methodology is applied to the investigation of how rate-laws describing
such reactions are transformed when the geometric scale is changed. The precise
purpose of such analysis is to investigate whether or not power-law rate-laws
present some remarkable features accounting for the successes of power-law
modeling. As we shall see, according to the renormalization group point of view
the answer is positive, as far as power-laws are the critical solutions of the
renormalization group transformation, namely power-law rate-laws are the
renormalization group invariant solutions. Moreover, it is shown that these
results also imply invariance under the group of concentration scalings, thus
accounting for the reported power-law model accuracy over several orders of
magnitude in metabolite concentrations.


  Dendritic spines, the carriers of long-term memory, occupy a small fraction
of cortical space, and yet they are the major consumers of brain metabolic
energy. What fraction of this energy goes for synaptic plasticity, correlated
with learning and memory? It is estimated here based on neurophysiological and
proteomic data for rat brain that, depending on the level of protein
phosphorylation, the energy cost of synaptic plasticity constitutes a small
fraction of the energy used for fast excitatory synaptic transmission,
typically $4.0-11.2 \%$. Next, this study analyzes a metabolic cost of a new
learning and its memory trace in relation to the cost of prior memories, using
a class of cascade models of synaptic plasticity. It is argued that these
models must contain bidirectional cyclic motifs, related to protein
phosphorylation, to be compatible with basic thermodynamic principles. For most
investigated parameters longer memories generally require proportionally more
energy to store. The exception are the parameters controlling the speed of
molecular transitions (e.g. ATP driven phosphorylation rate), for which memory
lifetime per invested energy can increase progressively for longer memories.
Furthermore, in general, a memory trace decouples dynamically from a
corresponding synaptic metabolic rate such that the energy expended on a new
learning and its memory trace constitutes in most cases only a small fraction
of the baseline energy associated with prior memories. Taken together, these
empirical and theoretical results suggest a metabolic efficiency of
synaptically stored information.


  {\it Caenorhabditis elegans} nematode worms are the only animals with the
known detailed neural connectivity diagram, well characterized genomics, and
relatively simple quantifiable behavioral output. With this in mind, many
researchers view this animal as the best candidate for a systems biology
approach, where one can integrate molecular and cellular knowledge to gain
global understanding of worm's behavior. This work reviews some research in
this direction, emphasizing computational perspective, and points out some
successes and challenges to meet this lofty goal.


  The Peyrard-Bishop DNA model describes the molecular interactions with simple
potentials which allow efficient calculations of melting temperatures. However,
it is based on a Hamiltonian that does not consider the helical twist or any
other relevant molecular dimensions. Here, we start from a more realistic 3D
model and work out several approximations to arrive at a new non-linear 1D
Hamiltonian with a twist angle dependence. Our approximations were numerically
compared to full 3D calculations, and established its validity in the regime of
small angles. For long DNA sequences we obtain sharp, first-order-like melting,
transitions.


  Phylogenetic trees represent evolutionary relationships among sets of
organisms. Popular phylogenetic reconstruction approaches typically yield
hundreds to thousands of trees on a common leafset. Storing and sharing such
large collection of trees requires considerable amount of space and bandwidth.
Furthermore, the huge size of phylogenetic tree databases can make search and
retrieval operations time-consuming. Phylogenetic compression techniques are
specialized compression techniques that exploit redundant topological
information to achieve better compression of phylogenetic trees. Here, we
present EvoZip, a new approach for phylogenetic tree compression. On average,
EvoZip achieves 71.6% better compression and takes 80.71% less compression time
and 60.47% less decompression time than TreeZip, the current state-of-the-art
algorithm for phylogenetic tree compression. While EvoZip is based on TreeZip,
it betters TreeZip due to (a) an improved bipartition and support list encoding
scheme, (b) use of Deflate compression algorithm, and (c) use of an efficient
tree reconstruction algorithm. EvoZip is freely available online for use by the
scientific community.


  Peptide nucleic acids (PNAs) are artificial nucleic acids with a peptide
backbone instead of sugar phosphate backbone of DNA or RNA. Their resistance to
degradation, selectivity and greater binding affinity in comparison to usual
nucleic acids led to consideration of their great potential for different
applications. For example, they can be used in molecular diagnostics and
antisense therapeutics such as antimicrobial agents or gene regulatory tools.
On the other hand, large hydrophilic property of PNA molecules, which inhibit
them to cross cell membranes readily, is an obstacle for their delivery to
considered target cells and a limiting criterion for their applications.
Therefore, PNA delivery technologies have been developing to hurdle this
limitation. For example, addition of lysine residues, charged membrane
penetrating peptide sequences, PNAs conjugated with antibodies or steroids,
cationic liposomes as carriers of PNA conjugates, protective peptides and
technology of photochemical internalization (PCI) as well as by the recent
technology of nanoparticle-based delivery have been employed. In this article
we compared different delivery technologies which can be applicable to PNAs. As
a result nanoparticle-based delivery showed more advantages in comparison to
others and its application is growing fast. Keywords: Target cell delivery,
Peptide nucleic acid, PNA, Delivery.


  Ecological systems show a variety of characteristic patterns of biodiversity
in space and time. It is a challenge for theory to find models that can
reproduce and explain the observed patterns. Since the advent of island
biogeography these models revolve around speciation, dispersal, and extinction,
but they usually neglect trophic structure. Here, we propose and study a
spatially extended evolutionary food web model that allows us to study large
spatial systems with several trophic layers. Our computer simulations show that
the model gives rise simultaneously to several biodiversity patterns in space
and time, from species abundance distributions to the waxing and waning of
geographic ranges. We find that trophic position in the network plays a crucial
role when it comes to the time evolution of range sizes, because the trophic
context restricts the occurrence and survival of species especially on higher
trophic levels.


  Self-assembly and force generation are two central processes in biological
systems that usually are considered in separation. However, the signals that
activate non-muscle myosin II molecular motors simultaneously lead to
self-assembly into myosin II minifilaments as well as progression of the motor
heads through the crossbridge cycle. Here we investigate theoretically the
possible effects of coupling these two processes. Our assembly model, which
builds upon a consensus architecture of the minifilament, predicts a critical
aggregation concentration at which the assembly kinetics slows down
dramatically. The combined model predicts that increasing actin filament
concentration and force both lead to a decrease in the critical aggregation
concentration. We suggest that due to these effects, myosin II minifilaments in
a filamentous context might be in a critical state that reacts faster to
varying conditions than in solution. We finally compare our model to
experiments by simulating fluorescence recovery after photobleaching.


  We present a mathematical study of the emergence of phenotypic heterogeneity
in vascularised tumours. Our study is based on formal asymptotic analysis and
numerical simulations of a system of non-local parabolic equations that
describes the phenotypic evolution of tumour cells and their nonlinear dynamic
interactions with the oxygen, which is released from the intratumoural vascular
network. Numerical simulations are carried out both in the case of arbitrary
distributions of intratumour blood vessels and in the case where the
intratumoural vascular network is reconstructed from clinical images obtained
using dynamic optical coherence tomography. The results obtained support a more
in-depth theoretical understanding of the eco-evolutionary process which
underpins the emergence of phenotypic heterogeneity in vascularised tumours. In
particular, our results offer a theoretical basis for empirical evidence
indicating that the phenotypic properties of cancer cells in vascularised
tumours vary with the distance from the blood vessels, and establish a relation
between the degree of tumour tissue vascularisation and the level of
intratumour phenotypic heterogeneity.


  Motivation: Identifying interaction clusters of large gene regulatory
networks (GRNs) is critical for its further investigation, while this task is
very challenging, attributed to data noise in experiment data, large scale of
GRNs, and inconsistency between gene expression profiles and function modules,
etc. It is promising to semi-supervise this process by prior information, but
shortage of prior information sometimes make it very challenging. Meanwhile, it
is also annoying, and sometimes impossible to discovery gold standard for
evaluation of clustering results.\\ Results: With assistance of an online
enrichment tool, this research proposes a semi-supervised hierarchical
clustering method via deconvolved correlation matrix~(SHC-DC) to discover
interaction clusters of large-scale GRNs. Three benchmark networks including a
\emph{Ecoli} network and two \emph{Yeast} networks are employed to test
semi-supervision scheme of the proposed method. Then, SHC-DC is utilized to
cluster genes in sleep study. Results demonstrates it can find interaction
modules that are generally enriched in various signal pathways. Besides the
significant influence on blood level of interleukins, impact of sleep on
important pathways mediated by them is also validated by the discovered
interaction modules.


  The prisoner's dilemma game is the most known contribution of game theory
into social sciences. Here we describe new implications of this game for
transactional and transformative leadership. While the autocratic
(Stackelberg's) leadership is inefficient for this game, we discuss a
Pareto-optimal scenario, where the leader L commits to react probabilistically
to pure strategies of the follower F, which is free to make the first move.
Offering F to resolve the dilemma, L is able to get a larger average pay-off.
The exploitation can be stabilized via repeated interaction of L and F, and
turns to be more stable than the egalitarian regime, where the pay-offs of L
and F are equal. The total (summary) pay-off of the exploiting regime is never
larger than in the egalitarian case. We discuss applications of this solution
to a soft method of fighting corruption and to modeling the Machiavellian
leadership. Whenever the defection benefit is large, the optimal strategies of
F are mixed, while the summary pay-off is maximal. One mechanism for sustaining
this solution is that L recognizes intentions of F.


  Auto-regulatory feedback loops are one of the most common network motifs. A
wide variety of stochastic models have been constructed to understand how the
fluctuations in protein numbers in these loops are influenced by the kinetic
parameters of the main biochemical steps. These models differ according to (i)
which sub-cellular processes are explicitly modelled; (ii) the modelling
methodology employed (discrete, continuous or hybrid); (iii) whether they can
be analytically solved for the steady-state distribution of protein numbers. We
discuss the assumptions and properties of the main models in the literature,
summarize our current understanding of the relationship between them and
highlight some of the insights gained through modelling.


  Several orders of magnitude typically separate the contour length of
eukaryotic chromosomes and the size of the nucleus where they are confined. The
ensuing topological constraints can slow down the relaxation dynamics of
genomic filaments to the point that mammalian chromosomes are never in
equilibrium over a cell's lifetime. In this opinion article, we revisit these
out-of-equilibrium effects and discuss how their inclusion in physical models
can enhance the spatial reconstructions of interphase eukaryotic genomes from
phenomenological constraints collected during interphase.


  Motivation: We introduce SneakySnake, a highly parallel and highly accurate
pre-alignment filter that remarkably reduces the need for computationally
costly sequence alignment. The key idea of SneakySnake is to reduce the
approximate string matching (ASM) problem to the single net routing (SNR)
problem in VLSI chip layout. In the SNR problem, we are interested in finding
the optimal path that connects two terminals with the least routing cost on a
special grid layout that contains obstacles. The SneakySnake algorithm quickly
solves the SNR problem and uses the found optimal path to decide whether or not
performing sequence alignment is necessary. Reducing the ASM problem into SNR
also makes SneakySnake efficient to implement on CPUs, GPUs, and FPGAs.
Results: SneakySnake significantly improves the accuracy of pre-alignment
filtering by up to four orders of magnitude compared to the state-of-the-art
pre-alignment filters, Shouji, GateKeeper, and SHD. For short sequences,
SneakySnake accelerates Edlib (state-of-the-art implementation of Myers's
bit-vector algorithm) and Parasail (state-of-the-art sequence aligner with a
configurable scoring function), by up to 37.7x and 43.9x (>12x on average),
respectively, with its CPU implementation, and by up to 413x and 689x (>400x on
average), respectively, with FPGA and GPU acceleration. For long sequences, the
CPU implementation of SneakySnake accelerates Parasail and KSW2 (sequence
aligner of minimap2) by up to 979x (276.9x on average) and 91.7x (31.7x on
average), respectively. As SneakySnake does not replace sequence alignment,
users can still obtain all capabilities (e.g., configurable scoring functions)
of the aligner of their choice, unlike existing acceleration efforts that
sacrifice some aligner capabilities. Availability:
https://github.com/CMU-SAFARI/SneakySnake


  A key step in the origin of life is the emergence of a primitive metabolism.
This requires the formation of a subset of chemical reactions that is both
self-sustaining and collectively autocatalytic. A generic theory to study such
processes (called 'RAF theory') has provided a precise and computationally
effective way to address these questions, both on simulated data and in
laboratory studies. One of the classic applications of this theory (arising
from Stuart Kauffman's pioneering work in the 1980s) involves networks of
polymers under cleavage and ligation reactions; in the first part of this
paper, we provide the first exact description of the number of such reactions
under various model assumptions. Conclusions from earlier studies relied on
either approximations or asymptotic counting, and we show that the exact counts
lead to similar (though not always identical) asymptotic results. In the second
part of the paper, we solve some questions posed in more recent papers
concerning the computational complexity of some key questions in RAF theory. In
particular, although there is a fast algorithm to determine whether or not a
catalytic reaction network contains a subset that is both self-sustaining and
autocatalytic (and, if so, find one), determining whether or not sets exist
that satisfy certain additional constraints exist turns out to be NP-complete.


  Histopathology is a reflection of the molecular changes and provides
prognostic phenotypes representing the disease progression. In this study, we
introduced feature scores generated from hematoxylin and eosin histology images
based on deep learning (DL) models developed for prostate pathology. We
demonstrated that these feature scores were significantly prognostic for time
to event endpoints (biochemical recurrence and cancer-specific survival) and
had simultaneously molecular biologic associations to relevant genomic
alterations and molecular subtypes using already trained DL models that were
not previously exposed to the datasets of the current study. Further, we
discussed the potential of such feature scores to improve the current tumor
grading system and the challenges that are associated with tumor heterogeneity
and the development of prognostic models from histology images. Our findings
uncover the potential of feature scores from histology images as digital
biomarkers in precision medicine and as an expanding utility for digital
pathology.


  The phase separation of oppositely-charged polyelectrolytes in solution is of
current interest . In this work we study the driving force for polyelectrolyte
complexation using molecular dynamics simulations. We calculate the potential
of mean force between poly(lysine) and poly(glutamate) oligomers using three
different forcefields, an atomistic force field and two coarse-grained force
fields. There is qualitative agreement between all forcefields, i.e., the sign
and magnitude of the free energy and the nature of the driving force are
similar, which suggests that the molecular nature of water does not play a
significant role. For fully charged peptides, we find that the driving force
for association is entropic in all cases when small ions either neutralize the
poly-ions, or are in excess. The removal of all counterions switches the
driving force, making complexation energetic. This suggests that the entropy of
complexation is dominated by the counterions. When only 6 residues of a 11-mer
are charged, however, the driving force is enthalpic in salt-free conditions.
The simulations shed insight into the mechanism of complex coacervation and the
importance of realistic models for the polyions.


  The curvature sensitive localization of proteins on membranes is vital for
many cell biological processes. Coarse-grained models are routinely employed to
study the curvature sensing phenomena and membrane morphology at the length
scale of few micrometers. Two prevalent phenomenological models exist for
modeling experimental observations of curvature sensing, (1) the spontaneous
curvature model and (2) the curvature mismatch model, which differ in their
treatment of the change in elastic energy due to the binding of proteins on the
membrane. In this work, the prediction of sensing and generation behaviour, by
these two models, are investigated using analytical calculations as well as
Dynamic Triangulation Monte Carlo simulations of quasi-spherical vesicles.
While the spontaneous curvature model yields a monotonically decreasing sensing
curve as a function of vesicle radius, the curvature mismatch model results in
a non-monotonic sensing curve. We highlight the main differences in the
interpretation of the protein-related parameters in the two models. We further
propose that the spontaneous curvature model is appropriate for modeling
peripheral proteins employing the hydrophobic insertion mechanism, with minimal
modification of membrane rigidity, while the curvature mismatch model is
appropriate for modeling curvature generation using scaffolding mechanism where
there is significant stiffening of the membrane due to protein binding.


  microRNAs (miRNAs) regulate gene expression at post-transcriptional level by
repressing target RNA molecules. Competition to bind miRNAs tends in turn to
correlate their targets, establishing effective RNA-RNA interactions that can
influence expression levels, buffer fluctuations and promote signal
propagation. Such a potential has been characterized mathematically for small
motifs both at steady state and \red{away from stationarity}. Experimental
evidence, on the other hand, suggests that competing endogenous RNA (ceRNA)
crosstalk is rather weak. Extended miRNA-RNA networks could however favour the
integration of many crosstalk interactions, leading to significant large-scale
effects in spite of the weakness of individual links. To clarify the extent to
which crosstalk is sustained by the miRNA interactome, we have studied its
emergent systemic features in silico in large-scale miRNA-RNA network
reconstructions. We show that, although generically weak, system-level
crosstalk patterns (i) are enhanced by transcriptional heterogeneities, (ii)
can achieve high-intensity even for RNAs that are not co-regulated, (iii) are
robust to variability in transcription rates, and (iv) are significantly
non-local, i.e. correlate weakly with miRNA-RNA interaction parameters.
Furthermore, RNA levels are generically more stable when crosstalk is
strongest. As some of these features appear to be encoded in the network's
topology, crosstalk may functionally be favoured by natural selection. These
results suggest that, besides their repressive role, miRNAs mediate a weak but
resilient and context-independent network of cross-regulatory interactions that
interconnect the transcriptome, stabilize expression levels and support
system-level responses.


  Spatial structure is one of the simplest and most studied ecological factors
that affect the evolution of cooperation. It has been shown that spatial
reciprocity promotes cooperation due to the formation of cooperative clusters,
which provide mutual support against defectors. The usual assumption is that of
constant population size, where no density-related effect is possible. Here, we
extend the investigation of density-related effects on the evolution of
cooperation. We integrate evolutionary game theory to the Lattice
Lotka-Volterra Model. In our model, the birth rate depends on the local density
and on the payoff accumulated in the interactions. We characterize the
evolution of cooperation in terms of the coexistence and the extinction of the
types. The main result is that cooperation is most favoured at moderate levels
of the death rate. Interestingly, defectors are extinct at values of the death
rate that are lower than that at which cooperators are extinct. When death
knocks the door, defectors are the first to perish, whereas cooperators stand
longer due to mutual support.


  Living systems evolve one mutation at a time, but a single mutation can alter
the effect of subsequent mutations. The underlying mechanistic determinants of
such epistasis are unclear. Here, we demonstrate that the physical dynamics of
a biological system can generically constrain epistasis. We analyze models and
experimental data on proteins and regulatory networks. In each, we find that if
the long-time physical dynamics is dominated by a slow, collective mode, then
the dimensionality of mutational effects is reduced. Consequently, epistatic
coefficients for different combinations of mutations are no longer independent,
even if individually strong. Such epistasis can be summarized as resulting from
a global non-linearity applied to an underlying linear trait, i.e., as global
epistasis. This constraint, in turn, reduces the ruggedness of the
sequence-to-function map. By providing a generic mechanistic origin for
experimentally observed global epistasis, our work suggests that slow
collective physical modes can make biological systems evolvable.


  During the last two decades, ecological speciation has been a major research
theme in evolutionary biology. Ecological speciation occurs when reproductive
isolation between populations evolves as a result of niche differentiation.
Phytophagous insects represent model systems for the study of this evolutionary
process. The host-plants on which these insects feed and often spend parts of
their life cycle constitute ideal agents of divergent selection for these
organisms. Adaptation to feeding on different host-plant species can
potentially lead to ecological specialization of populations and subsequent
speciation. This process is thought to have given birth to the astonishing
diversity of phytophagous insects and is often put forward in macroevolutionary
scenarios of insect diversification. Consequently, numerous phylogenetic
studies on phytophagous insects have aimed at testing whether speciation driven
by host-plant adaptation is the main pathway for the diversification of the
groups under investigation. The increasing availability of comprehensive and
well-resolved phylogenies and the recent developments in phylogenetic
comparative methods are offering an unprecedented opportunity to test
hypotheses on insect diversification at a macroevolutionary scale, in a robust
phylogenetic framework. Our purpose here is to review the contribution of
phylogenetic analyses to investigate the importance of plant-mediated
speciation in the diversification of phytophagous insects and to present
suggestions for future developments in this field.


  Gene interaction graphs aim to capture various relationships between genes
and represent decades of biology research. When trying to make predictions from
genomic data, those graphs could be used to overcome the curse of
dimensionality by making machine learning models sparser and more consistent
with biological common knowledge. In this work, we focus on assessing whether
those graphs capture dependencies seen in gene expression data better than
random. We formulate a condition that graphs should satisfy to provide a good
prior knowledge and propose to test it using a `Single Gene Inference' (SGI)
task. We compare random graphs with seven major gene interaction graphs
published by different research groups, aiming to measure the true benefit of
using biologically relevant graphs in this context. Our analysis finds that
dependencies can be captured almost as well at random which suggests that, in
terms of gene expression levels, the relevant information about the state of
the cell is spread across many genes.


  Recent work has found that the behavior of an individual can be altered when
infected by a parasite. Here we explore the question: under what conditions, in
principle, can a general parasitic infection control system-wide social
behaviors? We analyze fixed points and hysteresis effects under the Master
Equation, with transitions between two behaviors given two different
subpopulations, healthy vs. parasitically-infected, within a population which
is kept fixed overall. The key model choices are: (i) the internal opinion of
infected humans may differ from that of the healthy population, (ii) the extent
that interaction drives behavioral changes may also differ, and (iii) indirect
interactions are most important. We find that the socio-configuration can be
controlled by the parasitically-infected population, under some conditions,
even if the healthy population is the majority and of opposite opinion.


  The ability of homologous chromosomes (or selected chromosomal loci) to pair
specifically in the apparent absence of DNA breakage and recombination
represents a prominent feature of eukaryotic biology. The mechanism of homology
recognition at the basis of such recombination-independent pairing has remained
elusive. A number of studies have supported the idea that sequence homology can
be sensed between intact DNA double helices in vivo. In particular, recent
analyses of the two silencing phenomena in fungi, known as repeat-induced point
mutation (RIP) and meiotic silencing by unpaired DNA (MSUD), have provided
genetic evidence for the existence of the direct homologous dsDNA-dsDNA
pairing. Both RIP and MSUD likely rely on the same search strategy, by which
dsDNA segments are matched as arrays of interspersed base-pair triplets. This
process is general and very efficient, yet it proceeds normally without the
RecA/Rad51/Dmc1 proteins. Further studies of RIP and MSUD may yield surprising
insights into the function of DNA in the cell.


  Many complex systems experience damage accumulation which leads to aging,
manifest as an increasing probability of system collapse with time. This
naturally raises the question of how to maximize health and longevity in an
aging system at minimal cost of maintenance and intervention. Here, we pose
this question in the context of a simple interdependent network model of aging
in complex systems, and use both optimal control theory and reinforcement
learning alongside a combination of analysis and simulation to determine
optimal maintenance protocols. These protocols may motivate the rational design
of strategies for promoting longevity in aging complex systems with potential
applications in therapeutic schedules and engineered system maintenance.


  Genetic data are often used to infer demographic history and changes or
detect genes under selection. Inferential methods are commonly based on models
making various strong assumptions: demography and population structures are
supposed \textit{a priori} known, the evolution of the genetic composition of a
population does not affect demography nor population structure, and there is no
selection nor interaction between and within genetic strains. In this paper, we
present a stochastic birth-death model with competitive interactions and
asexual reproduction. We develop an inferential procedure for ecological,
demographic and genetic parameters. We first show how genetic diversity and
genealogies are related to birth and death rates, and to how individuals
compete within and between strains. {This leads us to propose an original model
of phylogenies, with trait structure and interactions, that allows multiple
merging}. Second, we develop an Approximate Bayesian Computation framework to
use our model for analyzing genetic data. We apply our procedure to simulated
data from a toy model, and to real data by analyzing the genetic diversity of
microsatellites on Y-chromosomes sampled from Central Asia human populations in
order to test whether different social organizations show significantly
different fertility.


  Histone variants regulate replication, transcription, DNA damage repair, and
chromosome segregation. Though widely accepted as a paradigm, it has not been
rigorously demonstrated that histone variants encode unique mechanical
properties. Here, we present a new theoretical approach called Minimal Cylinder
Analysis (MCA) to determine the Young's modulus of nucleosomes from all-atom
Molecular Dynamics (MD) simulations. Recently, we validated this computational
analysis against in vitro single-molecule nanoindentation of histone variant
nucleosomes. In this report, we further extend MCA to study the biophysical
properties of hybrid nucleosomes that are known to exist in human cancer cells
and contain H3 histone variants CENP-A and H3.3. We investigate the mechanism
by which the elasticity of these heterotypic nucleosomes augments cryptic
binding surfaces. Further, we derive biological predictions that might arise
when such heterotypic nucleosomes take over large parts of the genome.


  In the dominant livestock systems of Sahelian countries herds have to move
across territories. Their mobility is often a source of conflict with farmers
in the areas crossed, and helps spread diseases such as Rift Valley Fever.
Knowledge of the routes followed by herds is therefore core to guiding the
implementation of preventive and control measures for transboundary animal
diseases, land use planning and conflict management. However, the lack of
quantitative data on livestock movements, together with the high temporal and
spatial variability of herd movements, has so far hampered the production of
fine resolution maps of animal movements. This paper proposes a general
framework for mapping potential paths for livestock movements and identifying
areas of high animal passage potential for those movements. The method consists
in combining the information contained in livestock mobility networks with
landscape connectivity, based on different mobility conductance layers. We
illustrate our approach with a livestock mobility network in Senegal and
Mauritania in the 2014 dry and wet seasons.


  DNA methylation is an important biological mechanism to regulate gene
expression and control cell development. Mechanistic modeling has become a
popular approach to enhance our understanding of the dynamics of methylation
pattern formation in living cells. Recent findings suggest that the methylation
state of a cytosine base can be influenced by its DNA neighborhood. Therefore,
it is necessary to generalize existing mathematical models that consider only
one cytosine and its partner on the opposite DNA-strand (CpG), in order to
include such neighborhood dependencies. One approach is to describe the system
as a stochastic automata network (SAN) with functional transitions. We show
that single-CpG models can successfully be generalized to multiple CpGs using
the SAN description and verify the results by comparing them to results from
extensive Monte-Carlo simulations.


  In this article, we investigate different parsimony-based approaches towards
finding recombination breakpoints in a multiple sequence alignment. This
recombination detection task is crucial in order to avoid errors in
evolutionary analyses caused by mixing together portions of sequences which had
a different evolution history. Following an overview of the field of
recombination detection, we formulate four computational problems for this task
with different objective functions. The four problems aim to minimize (1) the
total homoplasy of all blocks (2) the maximum homoplasy per block (3) the total
homoplasy ratio of all blocks and (4) the maximum homoplasy ratio per block. We
describe algorithms for each of these problems, which are fixed-parameter
tractable (FPT) when the characters are binary. We have implemented and tested
the algorithms on simulated data, showing that minimizing the total homoplasy
gives, in most cases, the most accurate results. Our implementation and
experimental data have been made publicly available. Finally, we also consider
the problem of combining blocks into non-contiguous blocks consisting of at
most p contiguous parts. Fixing the homoplasy h of each block to 0, we show
that this problem is NP-hard when p >= 3, but polynomial-time solvable for p =
2. Furthermore, the problem is FPT with parameter h for binary characters when
p = 2. A number of interesting problems remain open.


  DNA methylation (meDNA) is a suspected modulator of alternative splicing,
while splicing in turn is involved in tumour formations nearly as frequently as
DNA mutations. Yet, the impact of meDNA on tumorigenesis via its effect on
splicing has not been thoroughly explored. Here, we find that HCT116 colon
carcinoma cells inactivated for the DNA methylases DNMT1 and DNMT3b undergo a
partial epithelial to mesenchymal transition (EMT) associated with alternative
splicing of the CD44 transmembrane receptor. The skipping of CD44 variant exons
is in part explained by altered expression or splicing of splicing and
chromatin factors. A direct effect of meDNA on alternative splicing was
sustained by transient depletion of DNMT1 and the methyl-binding genes MBD1,
MBD2, and MBD3. Yet, local changes in intragenic meDNA also altered recruitment
of MBD1 protein and of the chromatin factor HP1$\gamma$ known to alter
transcriptional pausing and alternative splicing decisions. We further tested
if meDNA level has sufficiently strong direct impact on the outcome of
alternative splicing to have a predictive value in the MCF10A model for breast
cancer progression and in patients with acute lymphoblastic leukemia (B ALL).
We found that a small number of differentially spliced genes mostly involved in
splicing and signal transduction is systematically correlated with local meDNA.
Altogether, our observations suggest that, although DNA methylation has
multiple avenues to alternative splicing, its indirect effect may be also
mediated through alternative splicing isoforms of these sensors of meDNA.


  Intrinsically disordered proteins (IDPs) are important for biological
functions. In contrast to folded proteins, molecular recognition among certain
IDPs is "fuzzy" in that their binding and/or phase separation are
stochastically governed by the interacting IDPs' amino acid sequences while
their assembled conformations remain largely disordered. To help elucidate a
basic aspect of this fascinating yet poorly understood phenomenon, the binding
of a homo- or hetero-dimeric pair of polyampholytic IDPs is modeled statistical
mechanically using cluster expansion. We find that the binding affinities of
binary fuzzy complexes in the model correlate strongly with a newly derived
simple "jSCD" parameter readily calculable from the pair of IDPs' sequence
charge patterns. Predictions by our analytical theory are in essential
agreement with coarse-grained explicit-chain simulations. This computationally
efficient theoretical framework is expected to be broadly applicable to
rationalizing and predicting sequence-specific IDP-IDP polyelectrostatic
interactions.


  Nature provides us with abundant examples of how large numbers of individuals
can make decisions without the coordination of a central authority. Social
insects, birds, fishes, and many other living collectives, rely on simple
interaction mechanisms to do so. They individually gather information from the
environment; small bits of a much larger picture that are then shared locally
among the members of the collective and processed together to output a commonly
agreed choice. Throughout evolution, Nature found solutions to collective
decision-making problems that are intriguing to engineers for their robustness
to malfunctioning or lost individuals, their flexibility in face of dynamic
environments, and their ability to scale with large numbers of members. In the
last decades, whereas biologists amassed large amounts of experimental
evidence, engineers took inspiration from these and other examples to design
distributed algorithms that, while maintaining the same properties of their
natural counterparts, come with guarantees on their performance in the form of
predictive mathematical models. In this paper, we review the fundamental
processes that lead to a collective decision. We discuss examples of collective
decisions in biological systems and show how similar processes can be
engineered to design artificial ones. During this journey, we review a
framework to design distributed decision-making algorithms that are modular,
can be instantiated and extended in different ways, and are supported by a suit
of predictive mathematical models.


  Tumor metastasis is one of the main factors responsible for the high fatality
rate of cancer. Metastasis can occur after malignant cells transition from the
epithelial phenotype to the mesenchymal phenotype. This transformation allows
cells to migrate via the circulatory system and subsequently settle in distant
organs after undergoing the reverse transition from the mesenchymal to the
epithelial phenotypes. The core gene regulatory network controlling these
transitions consists of a system made up of coupled SNAIL/miRNA-34 and
ZEB1/miRNA-200 subsystems. In this work, we formulate a mathematical model of
the core regulatory motif and analyze its long-term behavior. We start by
developing a detailed reaction network with 24 state variables. Assuming fast
promoter and mRNA kinetics, we then show how to reduce our model to a monotone
four-dimensional system. For the reduced system, monotone dynamical systems
theory can be used to prove generic convergence to the set of equilibria for
all bounded trajectories. The theory does not apply to the full model, which is
not monotone, but we briefly discuss results for singularly-perturbed monotone
systems that provide a tool to extend convergence results from reduced to full
systems, under appropriate time separation assumptions.


  Tiered graph autoencoders provide the architecture and mechanisms for
learning tiered latent representations and latent spaces for molecular graphs
that explicitly represent and utilize groups (e.g., functional groups). This
enables the utilization and exploration of tiered molecular latent spaces,
either individually - the node (atom) tier, the group tier, or the graph
(molecule) tier - or jointly, as well as navigation across the tiers. In this
paper, we discuss the use of tiered graph autoencoders together with graph
prediction for molecular graphs. We show features of molecular graphs used, and
groups in molecular graphs identified for some sample molecules. We briefly
review graph prediction and the QM9 dataset for background information, and
discuss the use of tiered graph embeddings for graph prediction, particularly
weighted group pooling. We find that functional groups and ring groups
effectively capture and represent the chemical essence of molecular graphs
(structures). Further, tiered graph autoencoders and graph prediction together
provide effective, efficient and interpretable deep learning for molecular
graphs, with the former providing unsupervised, transferable learning and the
latter providing supervised, task-optimized learning.


  Epidemiological contact network models have emerged as an important tool in
understanding and predicting the spread of infectious disease, due to their
capacity to engage individual heterogeneity that may underlie essential
dynamics of a particular host-pathogen system. Just as fundamental are the
changes that real-world contact networks undergo over time, both independently
of and in response to pathogen spreading. These dynamics play a central role in
determining whether a disease will die out or become an epidemic within a
population, known as the epidemic threshold. In this paper, we provide an
overview of methods to predict the epidemic threshold for temporal contact
network models, and discuss areas that remain unexplored.


  Implicit in the RNA world hypothesis is that prebiotic RNA synthesis, despite
occurring in an environment without biochemical catalysts, produced the long
RNA polymers which are essential to the formation of life. In order to
investigate the prebiotic formation of long RNA polymers, we consider a general
solution of functionally identical monomer units that are capable of bonding to
form linear polymers by a step-growth process. Under the assumptions that (1)
the solution is well-mixed and (2) bonding/unbonding rates are independent of
polymerization state, the concentration of each length of polymer follows the
geometric Flory-Schulz distribution. We consider the rate dynamics that produce
this equilibrium; connect the rate dynamics, Gibbs free energy of bond
formation, and the bonding probability; solve the dynamics in closed form for
the representative special case of a Flory-Schulz initial condition; and
demonstrate the effects of imposing a maximum polymer length. Afterwards, we
derive a lower bound on the error introduced by truncation and compare this
lower bound to the actual error found in our simulation. Finally, we suggest
methods to connect these theoretical predictions to experimental results.


  Selection in a time-periodic environment is modeled via the continuous-time
two-player replicator dynamics, which for symmetric pay-offs reduces to the
Fisher equation of mathematical genetics. For a sufficiently rapid and cyclic
[fine-grained] environment, the time-averaged population frequencies are shown
to obey a replicator dynamics with a non-linear fitness that is induced by
environmental changes. The non-linear terms in the fitness emerge due to
populations tracking their time-dependent environment. These terms can induce a
stable polymorphism, though they do not spoil the polymorphism that exists
already without them. In this sense polymorphic populations are more robust
with respect to their time-dependent environments. The overall fitness of the
problem is still given by its time-averaged value, but the emergence of
polymorphism during genetic selection can be accompanied by decreasing mean
fitness of the population. The impact of the uncovered polymorphism scenario on
the models of diversity is examplified via the rock-paper-scissors dynamics,
and also via the prisoner's dilemma in a time-periodic environment.


  Mathematical epidemiological models have a broad use, including both
qualitative and quantitative applications. With the increasing availability of
data, large-scale quantitative disease spread models can nowadays be
formulated. Such models have a great potential, e.g., in risk assessments in
public health. Their main challenge is model parameterization given
surveillance data, a problem which often limits their practical usage.
  We offer a solution to this problem by developing a Bayesian methodology
suitable to epidemiological models driven by network data. The greatest
difficulty in obtaining a concentrated parameter posterior is the quality of
surveillance data; disease measurements are often scarce and carry little
information about the parameters. The often overlooked problem of the model's
identifiability therefore needs to be addressed, and we do so using a hierarchy
of increasingly realistic known truth experiments.
  Our proposed Bayesian approach performs convincingly across all our synthetic
tests. From pathogen measurements of shiga toxin-producing Escherichia coli
O157 in Swedish cattle, we are able to produce an accurate statistical model of
first-principles confronted with data. Within this model we explore the
potential of a Bayesian public health framework by assessing the efficiency of
disease detection and -intervention scenarios.


  Protein identification and profiling is critical for the advancement of cell
and molecular biology as well as medical diagnostics. Although mass
spectrometry and protein microarrays are commonly used for protein
identification, both methods require extensive experimental steps and long data
analysis times. Here we present a microfluidic top down proteomics platform
giving multidimensional read outs of the essential amino acids of proteins. We
obtain hydrodynamic radius and fluorescence signals relating to the content of
tryptophans, tyrosines and lysines of proteins using a combination of
diffusional sizing of proteins, label-free detection and on-chip labelling of
proteins with a latent fluorophore in the solution phase. We thereby achieve
identification of proteins on a single microfluidic chip by separating and
mapping proteins in multidimensional space based on their characteristic
physical parameters. Our results have significant implications in the
development of easy and rapid platforms to use for native protein
identification in clinical and laboratory settings.


  The diverse repertoire of T-cell receptors (TCR) plays a key role in the
adaptive immune response to infections. Previous studies show that secondary
responses to the yellow fever vaccine - the model for acute infection in humans
- are weaker than primary ones, but only quantitative measurements can describe
the concentration changes and lineage fates for distinct T-cell clones in vivo
over time. Using TCR alpha and beta repertoire sequencing for T-cell subsets,
as well as single-cell RNAseq and TCRseq, we track the concentrations and
phenotypes of individual T-cell clones in response to primary and secondary
yellow fever immunization showing their large diversity. We confirm the
secondary response is an order of magnitude weaker, albeit $\sim10$ days faster
than the primary one. Estimating the fraction of the T-cell response directed
against the single immunodominant epitope, we identify the sequence features of
TCRs that define the high precursor frequency of the two major TCR motifs
specific for this particular epitope. We also show the consistency of clonal
expansion dynamics between bulk alpha and beta repertoires, using a new
methodology to reconstruct alpha-beta pairings from clonal trajectories.


  Counting how many particles pass through a specific space within a specific
time is an interesting question in applied physics and social science. Here a
logistic model is developed to estimate the total number of flowing particles.
This model sheds light on a collective contribution of particle growth rate and
transient probability within a specific space in particle counting. This model
may offer a basic concept to understand transport dynamics of flowing
particles.


  Thyroid cancer is the commonest endocrine malignancy. Mutation in the BRAF
serine/threonine kinase is the most frequent genetic alteration in thyroid
cancer. Target therapy for advanced and poorly differentiated thyroid
carcinomas include BRAF pathway inhibitors. Here, we evaluated the role of
MAP3K8 expression as a potential driver of resistance to BRAF inhibition in
thyroid cancer. By analyzing Gene Expression Omnibus data repository, across
all thyroid cancer histotypes, we found that MAP3K8 is up-regulated in poorly
differentiated thyroid carcinomas and its expression is related to a stem cell
like phenotype and a poorer prognosis and survival. Taken together these data
unravel a novel mechanism for thyroid cancer progression and chemo-resistance
and confirm previous results obtained in cultured thyroid cancer stem cells


  We present a theoretical formalism to study steady state information
transmission in type 1 coherent feed-forward loop motif with an additive signal
integration mechanism. Our construct allows a two-step cascade to be slowly
transformed into a bifurcation network via a feed-forward loop which is a
prominent network motif. Utilizing a Gaussian framework, we show that the
feed-forward loop motif harnesses the maximum amount of Shannon mutual
information fractions constructed between the input signaling species and the
final gene product and also between the two downstream gene products. Our
findings may explain why a feed-forward loop is naturally selected by the force
of evolution as a motif in bacterial transcription control network while the
two-step cascade and bifurcation network are not.


  Reaction networks are commonly used within the mathematical biology and
mathematical chemistry communities to model the dynamics of interacting
species. These models differ from the typical graphs found in random graph
theory since their vertices are constructed from elementary building blocks,
i.e., the species. In this paper, we consider these networks in an Erd\H
os-R\'enyi framework and, under suitable assumptions, derive a threshold
function for the network to have a deficiency of zero, which is a property of
great interest in the reaction network community. Specifically, if the number
of species is denoted by $n$ and if the edge probability is denote by $p_n$,
then we prove that the probability of a random binary network being deficiency
zero converges to 1 if $\frac{p_n}{r(n)}\to 0$, as $n \to \infty$, and
converges to 0 if $\frac{p_n}{r(n)}\to \infty$, as $n \to \infty$, where
$r(n)=\frac{1}{n^3}$.


  Mountain ecosystems are sensitive indicators of climate change. Long-term
studies may be extremely useful in assessing the responses of high-elevation
ecosystems to climate change and other anthropogenic drivers. Mountain research
sites within the LTER (Long-Term Ecosystem Research) network are representative
of various types of ecosystems and span a wide bioclimatic and elevational
range. Here, we present a synthesis and a review of the main results from
long-term ecological studies in mountain ecosystems at 20 LTER sites in Italy,
Switzerland and Austria. We analyzed a set of key climate parameters, such as
temperature and snow cover duration, in relation to vascular species
composition, plant traits, abundance patterns, pedoclimate, nutrient dynamics
in soils and water, phenology and composition of freshwater biota. The overall
results highlight the rapid response of mountain ecosystems to climate change.
As temperatures increased, vegetation cover in alpine and subalpine summits
increased as well. Years with limited snow cover duration caused an increase in
soil temperature and microbial biomass during the growing season. Effects on
freshwater ecosystems were observed, in terms of increases in solutes,
decreases in nitrates and changes in plankton phenology and benthos
communities. This work highlights the importance of comparing and integrating
long-term ecological data collected in different ecosystems, for a more
comprehensive overview of the ecological effects of climate change.
Nevertheless, there is a need for i) adopting co-located monitoring site
networks to improve our ability to obtain sound results from cross-site
analysis, ii) carrying out further studies, with fine spatial and temporal
resolutions to improve understanding of responses to extreme events, and iii)
increasing comparability and standardizing protocols across networks to clarify
local from global patterns.


  The Loreau and Hector (2019) Comment on our paper (Pillai and Gouhier, 2019)
failed to address the two core elements of our critique, both the circularity
of the BEF research program, in general, and the mathematical flaws of the
Loreau-Hector partitioning scheme, in particular. Loreau and Hector avoided
dealing with the first part of our critique by arguing against a non-existent
claim that all biodiversity effects could be reduced to coexistence, while the
mathematical flaws in the Loreau-Hector partitioning method that we described
in the second part of our critique were ignored altogether. Here, we address
these misconceptions and demonstrate that all of the claims that were made in
our original paper hold. We conclude that (i) BEF studies need to adopt
baselines that account for coexistence in order to avoid overestimating the
effects of biodiversity and (ii) the LH partitioning method should not be used
unless the linearity of the abundance-ecosystem functioning relationship in
monocultures can be verified for all species.


  Large bottlebrush complexes formed from the polysaccharide hyaluronan (HA)
and the proteoglycan aggrecan contribute to cartilage compression resistance
and are necessary for healthy joint function. A variety of mechanical forces
act on these complexes in the cartilage extracellular matrix, motivating the
need for a quantitative description which links their structure and mechanical
response. Studies using electron microscopy have imaged the HA-aggrecan brush
but require adsorption to a surface, dramatically altering the complex from its
native conformation. We use magnetic tweezers force spectroscopy to measure
changes in extension and mechanical response of an HA chain as aggrecan
monomers bind and form a bottlebrush. This technique directly measures changes
undergone by a single complex with time and under varying solution conditions.
Upon addition of aggrecan, we find a large swelling effect manifests when the
HA chain is under very low external tension (i.e. stretching forces less than
~1 pN). We use models of force-extension behavior to show that repulsion
between the aggrecans induces an internal tension in the HA chain. Through
reference to theories of bottlebrush polymer behavior, we demonstrate that the
experimental values of internal tension are consistent with a polydisperse
aggrecan population, likely caused by varying degrees of glycosylation. By
enzymatically deglycosylating aggrecan, we show that aggrecan glycosylation is
the structural feature which causes HA stiffening. We then construct a simple
stochastic binding model to show that variable glycosylation leads to a wide
distribution of internal tensions in HA, causing variations in the mechanics at
much longer length-scales. Our results provide a mechanistic picture of how
flexibility and size of HA and aggrecan lead to the brush architecture and
mechanical properties of this important component of cartilage.


  To synthesize peptides alongside the RNAs making the so-called RNA world,
some genetic coding involving RNA had to develop. Herein, it is proposed that
the first real-coding setup was a direct one, made up of continuous
poly-tRNA-like molecules, with each tRNA-like moiety carrying, beyond and near
its 5 prime or 3 prime end, a trinucleotide site for specific amino acid
binding: the sequence and continuity of the tRNA moieties of a particular
poly-tRNA would ensure the sequence and continuity of the amino acids of the
corresponding peptide or small protein. In parallel with these particular
entities, and enhancing their peptide-forming function, a proto-ribosome and
primitive amino acid-activation system would develop. At some stage, one
critical innovation would be the appearance of RNA fragments that could tighten
several adjacent tRNA moieties together on a particular poly-tRNA molecule, by
pairing with the second trinucleotide sequence (identical to the first one
carrying the specific amino acid-binding site) situated at, or close to, the
middle of each tRNA moiety (i.e., the present anticodon site). These RNA
fragments, acting as authentic co-ribozymes in the peptide-synthesizing
apparatus, would constitute the ancestors of the present mRNAs. Later, on these
mRNA-like guiding fragments, free tRNA forms would be additionally used, first
keeping their amino acid-binding sites, then losing them in favor of a specific
amino acid attachment at a CCA arm at their 3 prime end. Finally, these latter
mechanisms would progressively prevail, leading to the modern and universal
indirect genetic coding system. Experimental and theoretical arguments are
presented and discussed in favor of such a scenario for the origin and
evolution of genetic coding.


  We demonstrate that the issues described in the Wagg et al. (2019) Comment on
our paper (Pillai and Gouhier, 2019) are all due to misunderstandings about the
implications of pairwise effects, the nature of the null baseline in both our
framework and in the Loreau-Hector (LH) partitioning scheme (i.e., the midpoint
of the monocultures), and the impact of nonlinearity on the LH partitioning
results. Specifically, we show that (i) pairwise effects can be computed over
any time horizon and thus do not imply stable coexistence, (ii) the midpoint of
the monocultures corresponds to a neutral community so coexistence was always
part of the LH baseline, and (iii) contrary to what Wagg et al. suggested,
generalized diversity-interaction models do not account for (and may in fact
exacerbate) the problem of nonlinearity in monocultures, which inflates the LH
net biodiversity effect and generates incorrect estimates of selection and
complementarity. Hence, all of our original claims about the triviality
inherent in biodiversity-ecosystem functioning research and the issues with the
LH partitioning scheme hold.


  A ubiquitous way that cells share information is by exchanging molecules.
Yet, the fundamental ways that this information exchange is influenced by
intracellular dynamics remain unclear. Here we use information theory to
investigate a simple model of two interacting cells with internal feedback. We
show that cell-to-cell molecule exchange induces a collective two-cell critical
point and that the mutual information between the cells peaks at this critical
point. Information can remain large far from the critical point on a manifold
of cellular states, but scales logarithmically with the correlation time of the
system, resulting in an information-correlation time tradeoff. This tradeoff is
strictly imposed, suggesting the correlation time as a proxy for the mutual
information.


  Many proteins in cells are capable of sensing and responding to piconewton
scale forces, a regime in which conformational changes are small but
significant for biological processes. In order to efficiently and effectively
sample the response of these proteins to small forces, enhanced sampling
techniques will be required. In this work, we derive, implement, and evaluate
an efficient method to simultaneously sample the result of applying any
constant pulling force within a specified range to a molecular system of
interest. We start from Simulated Tempering in Force, whereby force is applied
as a linear bias on a collective variable to the system's Hamiltonian, and the
coefficient is taken as a continuous auxiliary degree of freedom. We derive a
formula for an average collective-variable-dependent force, which depends on a
set of weights, learned on-the-fly throughout a simulation, that reflect the
limit where force varies infinitely quickly. These weights can then be used to
retroactively compute averages of any observable at any force within the
specified range. This technique is based on recent work deriving similar
equations for Infinite Switch Simulated Tempering in Temperature, that showed
the infinite switch limit is the most efficient for sampling. Here, we
demonstrate that our method accurately and simultaneously samples molecular
systems at all forces within a user defined force range, and show how it can
serve as an enhanced sampling tool for cases where the pulling direction
destabilizes states of low free-energy at zero-force. This method is
implemented in, and will be freely-distributed with, the PLUMED open-source
sampling library, and hence can be readily applied to problems using a wide
range of molecular dynamics software packages.


  Chromosome organisation is increasingly recognised as an essential component
of genome regulation, cell fate and cell health. Within the realm of
transposable elements (TEs) however, the spatial information of how genomes are
folded is still only rarely integrated in experimental studies or accounted for
in modelling. Here, we propose a new predictive modelling framework for the
study of the integration patterns of TEs based on extensions of widely employed
polymer models for genome organisation. Whilst polymer physics is recognised as
an important tool to understand the mechanisms of genome folding, we now show
that it can also offer orthogonal and generic insights into the integration and
distribution profiles (or "topography") of TEs across organisms. Here, we
present polymer physics arguments and molecular dynamics simulations on TEs
inserting into heterogeneously flexible polymers and show with a simple model
that polymer folding and local flexibility affects TE integration patterns. The
preliminary discussion presented herein lay the foundations for a large-scale
analysis of TE integration dynamics and topography as a function of the
three-dimensional host genome.


  Understanding the spatial organisation of the genome in the cell nucleus is
one of the current grand challenges in biophysics. Certain biochemical -- or
epigenetic -- marks that are deposited along the genome are thought to play an
important, yet poorly understood, role in determining genome organisation and
cell identity. The physical principles underlying the interplay between
epigenetic dynamics and genome folding remain elusive. Here we propose and
study a theory that assumes a coupling between epigenetic mark and genome
densities, and which can be applied at the scale of the whole nucleus. We show
that equilibrium models are not compatible with experiments and a qualitative
agreement is recovered by accounting for non-equilibrium processes which can
stabilise microphase separated epigenomic domains. We finally discuss the
potential biophysical origin of these terms.


  Mass-spectrometry technologies are widely used in the fields of ionomics and
metabolomics to simultaneously profile at the genome scale intracellular
concentrations of e.g. amino acids or elements. Short profiles of molecular or
sub-molecular features are intrinsically non-Gaussian and may reveal patterns
of correlations that reflect the system nature of the cell biochemistry and
biology. Here we introduce two profile similarity measures that enforce
information from the empirical covariance matrix of the data, the Mahalanobis
cosine and the hybrid-Mahalanobis cosine. We evaluate the performance of these
similarity measures in the task of inferring and integrating genetic networks
from omics data by analysing experimental datasets derived from the ionome and
the metabolome of the model organism S. cerevisiae, and several large curated
databases of genetic annotations. The proposed covariance-based similarity
measures can in general recover known and predicted associations between genes
better than the commonly used Pearson's correlation and the standard cosine
similarity. The choice of which of the two measures to recommend depends upon
whether the focus is on extracting genetic associations at a global or local
genetic network scale.


  The quality of biomolecular simulations critically depends on the accuracy of
the force field used to calculate the potential energy of the molecular
configurations. Currently, most simulations employ non-polarisable force
fields, which describe electrostatic interactions as the sum of Coulombic
interactions between fixed atomic charges. Polarization of these charge
distributions is incorporated only in a mean-field manner. In the past decade,
extensive efforts have been devoted to developing simple, efficient, and yet
generally applicable polarisable force fields for biomolecular simulations. In
this review, we summarise the latest developments in accounting for key
biomolecular interactions with polarisable force fields and applications to
address challenging biological questions. In the end, we provide an outlook for
future development in polarisable force fields.


  Several theoretical frameworks have been proposed to explain observed
biodiversity patterns, ranging from the classical niche-based theories, mainly
employing a continuous formalism, to neutral theories, based on statistical
mechanics of discrete communities. Differences in the descriptions of
biodiversity can arise due to the discrete or continuous nature of the
underlying models and the way internal or external perturbations appear in
their formulations. Here, we trace the effects of stochastic population
dynamics on biodiversity, from the scale of the individuals to the community
and based on both discrete and continuous representations of the system, by
consistently using measures of community diversity like the species abundance
distribution and the rank abundance curve and applying them to both discrete
and continuous populations. A novel measure, the community abundance
distribution, is introduced to facilitate the comparison across different
levels of description, from microscopic to macroscopic. Using a simple birth
and death process and an interacting population model, we highlight
discrepancies in their discrete and continuous distributions and discuss
relevant implications for the analysis of rare species and extinction dynamics.
Quantitative consideration of these issues is useful for better understanding
of the contributions of non-neutral processes and the mathematical
approximations to various measures of biodiversity.


  Why are living systems complex? Why does the biosphere contain living beings
with complexity features beyond those of the simplest replicators? What kind of
evolutionary pressures result in more complex life forms? These are key
questions that pervade the problem of how complexity arises in evolution. One
particular way of tackling this is grounded in an algorithmic description of
life: living organisms can be seen as systems that extract and process
information from their surroundings in order to reduce uncertainty. Here we
take this computational approach using a simple bit string model of coevolving
agents and their parasites. While agents try to predict their worlds, parasites
do the same with their hosts. The result of this process is that, in order to
escape their parasites, the host agents expand their computational complexity
despite the cost of maintaining it. This, in turn, is followed by increasingly
complex parasitic counterparts. Such arms races display several qualitative
phases, from monotonous to punctuated evolution or even ecological collapse.
Our minimal model illustrates the relevance of parasites in providing an active
mechanism for expanding living complexity beyond simple replicators, suggesting
that parasitic agents are likely to be a major evolutionary driver for
biological complexity.


  How should dispersal strategies be chosen to increase the likelihood of
survival of a species? We obtain the answer for the spatially extended versions
of three well-known models of two competing species with unequal diffusivities.
Though identical at the mean-field level, the three models exhibit drastically
different behaviour leading to different optimal strategies for survival, with
or without a selective advantage for one species. With conserved total particle
number, dispersal has no effect on survival probability. With a fluctuating
number, faster dispersal is advantageous if intra-species competition is
present, while moving slower is the optimal strategy for the disadvantaged
species if there is no intra-species competition: it is imperative to include
fluctuations to properly formulate survival strategies.


  We consider a model of extracellular signal-regulated kinase (ERK) regulation
by dual-site phosphorylation and dephosphorylation, which exhibits bistability
and oscillations, but loses these properties in the limit in which the
mechanisms underlying phosphorylation and dephosphorylation become processive.
Our results suggest that anywhere along the way to becoming processive, the
model remains bistable and oscillatory. More precisely, in simplified versions
of the model, precursors to bistability and oscillations (specifically,
multistationarity and Hopf bifurcations, respectively) exist at all
"processivity levels". Finally, we investigate whether bistability and
oscillations can exist together.


  One algorithm to predict protein structure is the residual dipolar coupling
based residue assembly and filter tool (REDCRAFT). This algorithm exploits an
exponential reduction of the search space of all possible structures to find a
structure that best fits a set of experimental residual dipolar couplings.
However, the minimum amount of data required to successfully determine a
protein's structure using REDCRAFT has not been previously investigated. Here
we explore the effect of reducing the amount of data used to fold proteins. Our
goal is to reduce experimental data collection times while retaining the
accuracy levels previously achieved with larger amounts of data. We also
investigate incorporating a priori secondary structure information into
REDCRAFT to improve its structure prediction ability.


  The multiple-biomarker classifier problem and its assessment are reviewed
against the background of some fundamental principles from the field of
statistical pattern recognition, machine learning, or the recently so-called
"data science". A narrow reading of that literature has led many authors to
neglect the contribution to the total uncertainty of performance assessment
from the finite training sample. Yet the latter is a fundamental indicator of
the stability of a classifier; thus its neglect may be contributing to the
problematic status of many studies. A three-level strategy is proposed for
moving forward in this field. The lowest level is that of construction, where
candidate features are selected and the choice of classifier architecture is
made. At that point, the effective dimensionality of the classifier is
estimated and used to size the next level of analysis, a pilot study on
previously unseen cases. The total (training and testing) uncertainty resulting
from the pilot study is, in turn, used to size the highest level of analysis, a
pivotal study with a target level of uncertainty. Some resources available in
the literature for implementing this approach are reviewed. Although the
concepts explained in the present article may be fundamental and
straightforward for many researchers in the machine learning community they are
subtle for many practitioners, for whom we provided a general advice for the
best practice in \cite{Shi2010MAQCII} and elaborate here in the present paper.


  Locusts are significant agricultural pests. Under favorable environmental
conditions flightless juveniles may aggregate into coherent, aligned swarms
referred to as hopper bands. These bands are often observed as a propagating
wave having a dense front with rapidly decreasing density in the wake. A
tantalizing and common observation is that these fronts slow and steepen in the
presence of green vegetation. This suggests the collective motion of the band
is mediated by resource consumption. Our goal is to model and quantify this
effect. We focus on the Australian plague locust, for which excellent field and
experimental data is available. Exploiting the alignment of locusts in hopper
bands, we concentrate solely on the density variation perpendicular to the
front. We develop two models in tandem; an agent-based model that tracks the
position of individuals and a partial differential equation model that
describes locust density. In both these models, locust are either stationary
(and feeding) or moving. Resources decrease with feeding. The rate at which
locusts transition between moving and stationary (and vice versa) is enhanced
(diminished) by resource abundance. This effect proves essential to the
formation, shape, and speed of locust hopper bands in our models. From the
biological literature we estimate ranges for the ten input parameters of our
models. Sobol sensitivity analysis yields insight into how the band's
collective characteristics vary with changes in the input parameters. By
examining 4.4 million parameter combinations, we identify biologically
consistent parameters that reproduce field observations. We thus demonstrate
that resource-dependent behavior can explain the density distribution observed
in locust hopper bands. This work suggests that feeding behaviors should be an
intrinsic part of future modeling efforts.


  Single-cell RNA-Sequencing (scRNA-Seq) has undergone major technological
advances in recent years, enabling the conception of various organism-level
cell atlassing projects. With increasing numbers of datasets being deposited in
public archives, there is a need to address the challenges of enabling the
reproducibility of such data sets. Here, we describe guidelines for a minimum
set of metadata to sufficiently describe scRNA-Seq experiments, ensuring
reproducibility of data analyses.


  It remains a challenging task to generate a vast variety of novel compounds
with desirable pharmacological properties. In this work, a generative network
complex (GNC) is proposed as a new platform for designing novel compounds,
predicting their physical and chemical properties, and selecting potential drug
candidates that fulfill various druggable criteria such as binding affinity,
solubility, partition coefficient, etc. We combine a SMILES string generator,
which consists of an encoder, a drug-property controlled or regulated latent
space, and a decoder, with verification deep neural networks, a target-specific
three-dimensional (3D) pose generator, and mathematical deep learning networks
to generate new compounds, predict their drug properties, construct 3D poses
associated with target proteins, and reevaluate druggability, respectively. New
compounds were generated in the latent space by either randomized output,
controlled output, or optimized output. In our demonstration, 2.08 million and
2.8 million novel compounds are generated respectively for Cathepsin S and BACE
targets. These new compounds are very different from the seeds and cover a
larger chemical space. For potentially active compounds, their 3D poses are
generated using a state-of-the-art method. The resulting 3D complexes are
further evaluated for druggability by a championing deep learning algorithm
based on algebraic topology, differential geometry, and algebraic graph
theories. Performed on supercomputers, the whole process took less than one
week. Therefore, our GNC is an efficient new paradigm for discovering new drug
candidates.


  Vector or pest control is essential to reduce the risk of vector-borne
diseases or crop losses. Among the available biological control tools, the
Sterile Insect Technique (SIT) is one of the most promising. However,
SIT-control campaigns must be carefully planned in advance in order to render
desirable outcomes. In this paper, we design SIT-control intervention programs
that can avoid the real-time monitoring of the wild population and require to
mass-rear a minimal overall number of sterile insects, in order to induce a
local elimination of the wild population in the shortest time. Continuous-time
release programs are obtained by applying an optimal control approach, and then
laying the groundwork of more practical SIT-control programs consisting of
periodic impulsive releases.


  Constructing gene regulatory networks is a critical step in revealing disease
mechanisms from transcriptomic data. In this work, we present NO-BEARS, a novel
algorithm for estimating gene regulatory networks. The NO-BEARS algorithm is
built on the basis of the NOTEARS algorithm with two improvements. First, we
propose a new constraint and its fast approximation to reduce the computational
cost of the NO-TEARS algorithm. Next, we introduce a polynomial regression loss
to handle non-linearity in gene expressions. Our implementation utilizes modern
GPU computation that can decrease the time of hours-long CPU computation to
seconds. Using synthetic data, we demonstrate improved performance, both in
processing time and accuracy, on inferring gene regulatory networks from gene
expression data.


  Motivation: Since long non-coding RNAs (lncRNAs) have involved in a wide
range of functions in cellular and developmental processes, an increasing
number of methods have been proposed for distinguishing lncRNAs from coding
RNAs. However, most of the existing methods are designed for lncRNAs in animal
systems, and only a few methods focus on the plant lncRNA identification.
Different from lncRNAs in animal systems, plant lncRNAs have distinct
characteristics. It is desirable to develop a computational method for accurate
and robust identification of plant lncRNAs. Results: Herein, we present a plant
lncRNA identification method ItLnc-BXE, which utilizes multiple features and
the ensemble learning strategy. First, a diversity of lncRNA features is
collected and filtered by feature selection to represent RNA transcripts. Then,
several base learners are trained and further combined into a single
meta-learner by ensemble learning, and thus an ItLnc-BXE model is constructed.
ItLnc-BXE models are evaluated on datasets of six plant species, the results
show that ItLnc-BXE outperforms other state-of-the-art plant lncRNA
identification methods, achieving better and robust performances (AUC>95.91%).
We also perform some experiments about cross-species lncRNA identification, and
the results indicate that dicots-based and monocots-based models can be used to
accurately identify lncRNAs in lower plant species, such as mosses and algae.
Availability: source codes are available at
https://github.com/BioMedicalBigDataMiningLab/ItLnc-BXE. Contact:
zhangwen@mail.hzau.edu.cn (or) zhangwen@whu.edu.cn Supplementary information:
Supplementary data are available at Bioinformatics online.


  Single-cell RNA sequencing data have complex features such as dropout events,
over-dispersion, and high-magnitude outliers, resulting in complicated
probability distributions of mRNA abundances that are statistically
characterized in terms of a zero-inflated negative binomial (ZINB) model. Here
we provide a mesoscopic kinetic foundation of the widely used ZINB model based
on the biochemical reaction kinetics underlying transcription. Using multiscale
modeling and simplification techniques, we show that the ZINB distribution of
mRNA abundance and the phenomenon of transcriptional bursting naturally emerge
from a three-state stochastic transcription model. We further reveal a
nontrivial quantitative relation between dropout events and transcriptional
bursting, which provides novel insights into how and to what extent the burst
size and burst frequency could reduce the dropout rate. Three different
biophysical origins of over-dispersion are also clarified at the single-cell
level.


  Despite many advances in computational modeling of protein structures, these
methods have not been widely utilized by experimental structural biologists.
Two major obstacles are preventing the transition from a purely-experimental to
a purely-computational mode of protein structure determination. The first
problem is that most computational methods need a large library of computed
structures that span a large variety of protein fold families, while structural
genomics initiatives have slowed in their ability to provide novel protein
folds in recent years. The second problem is an unwillingness to trust
computational models that have no experimental backing. In this paper we test a
potential solution to these problems that we have called Probability Density
Profile Analysis (PDPA) that utilizes unassigned residual dipolar coupling data
that are relatively cheap to acquire from NMR experiments.


  Lung cancer is the leading cause of the largest number of deaths worldwide
and lung adenocarcinoma (LUAD) is the most common form of lung cancer. In this
study, we carried out an integrated meta-analysis of the mutations including
single-nucleotide variations (SNVs), the copy number variations (CNVs), RNA-seq
and clinical data of patients with LUAD downloaded from The Cancer Genome Atlas
(TCGA). We integrated significant SNV and CNV genes, differentially expressed
genes (DEGs) and the DEGs in active subnetworks to construct a prognosis
signature. Cox proportional hazards model (LOOCV) with Lasso penalty was used
to identify the best gene signature among different gene categories. The
patients in both training and test data were clustered into high-risk and
low-risk groups by using risk scores of the patients calculated based on
selected gene signature. We generated a 12-gene signature (DEPTOR, ZBTB16,
BCHE, MGLL, MASP2, TNNI2, RAPGEF3, SGK2, MYO1A, CYP24A1, PODXL2, CCNA1) for
overall survival prediction. The survival time of high-risk and low-risk groups
was significantly different. This 12-gene signature could predict prognosis and
they are potential predictors for the survival of the patients with LUAD.


  Resonance assignment is a critical first step in the investigation of protein
structures using NMR spectroscopy. The development of assignment methods that
require less experimental data is possible with prior knowledge of the
macromolecular structure. Automated methods of performing the task of resonance
assignment can significantly reduce the financial cost and time requirement for
protein structure determination. Such methods can also be beneficial in
validating a protein's solution state structure. Here we present a new approach
to the assignment problem. Our approach uses only RDC data to assign backbone
resonances. It provides simultaneous order tensor estimation and assignment.
Our approach compares independent order tensor estimates to determine when the
correct order tensor has been found. We demonstrate the algorithm's viability
using simulated data from the protein domain 1A1Z.


  It is well known that excessive harvesting or hunting has driven species to
extinction both on local and global scales. This leads to one of the
fundamental problems of conservation ecology: how should we harvest a
population so that economic gain is maximized, while also ensuring that the
species is safe from extinction? We study an ecosystem of interacting species
that are influenced by random environmental fluctuations. At any point in time,
we can either harvest or seed (repopulate) species. Harvesting brings an
economic gain while seeding incurs a cost. The problem is to find the optimal
harvesting-seeding strategy that maximizes the expected total income from
harvesting minus the cost one has to pay for the seeding of various species. We
consider what happens when one, or both, of the seeding and harvesting rates
are bounded. The focus of this paper is the analysis of these three novel
settings: bounded seeding and infinite harvesting, bounded seeding and bounded
harvesting, and infinite seeding and bounded harvesting.
  We prove analytical results and develop numerical approximation methods. By
implementing these approximations, we are able to gain qualitative information
about how to best harvest and seed species. We are able to show that in the
single species setting there are thresholds $0<L_1<L_2<\infty$ such that: 1) if
the population size is `low', so that it lies in $(0, L_1]$, there is seeding
using the maximal seeding rate; 2) if the population size `moderate', so that
it lies in $(L_1,L_2)$, there is no harvesting or seeding; 3) if the population
size is `high', so that it lies in the interval $[L_2, \infty)$, there is
harvesting using the maximal harvesting rate. Once we have a system with at
least two species, numerical experiments show that constant threshold
strategies are not optimal anymore.


  The amount of sequence data obtained from ancient samples has dramatically
expanded in the last decade, and so have the types of questions that can now be
addressed using ancient DNA. In the field of human history, while ancient DNA
has provided answers to long-standing debates about major movements of people,
it has also recently begun to inform on other important facets of the human
experience. The field is now moving from not only focusing on large-scale
supra-regional studies to also taking a more local perspective, shedding light
on socioeconomic processes, inheritance rules, marriage practices and
technological diffusion. In this review, we summarize recent studies showcasing
these types of insights, focusing on the methods used to infer sociocultural
aspects of human behaviour. This work often involves working across disciplines
that have, until recently, evolved in separation. We argue that
multidisciplinary dialogue is crucial for a more integrated and richer
reconstruction of human history, as it can yield extraordinary insights about
past societies, reproductive behaviours and even lifestyle habits that would
not have been possible to obtain otherwise.


  Recently, molecular fingerprints extracted from three-dimensional (3D)
structures using advanced mathematics, such as algebraic topology, differential
geometry, and graph theory have been paired with efficient machine learning,
especially deep learning algorithms to outperform other methods in drug
discovery applications and competitions. This raises the question of whether
classical 2D fingerprints are still valuable in computer-aided drug discovery.
This work considers 23 datasets associated with four typical problems, namely
protein-ligand binding, toxicity, solubility and partition coefficient to
assess the performance of eight 2D fingerprints. Advanced machine learning
algorithms including random forest, gradient boosted decision tree, single-task
deep neural network and multitask deep neural network are employed to construct
efficient 2D-fingerprint based models. Additionally, appropriate consensus
models are built to further enhance the performance of 2D-fingerprintbased
methods. It is demonstrated that 2D-fingerprint-based models perform as well as
the state-of-the-art 3D structure-based models for the predictions of toxicity,
solubility, partition coefficient and protein-ligand binding affinity based on
only ligand information. However, 3D structure-based models outperform 2D
fingerprint-based methods in complex-based protein-ligand binding affinity
predictions.


  The chemical master equation (CME) is a fundamental description of
interacting molecules commonly used to model chemical kinetics and noisy gene
regulatory networks. Exact time-dependent solutions of the CME -- which
typically consists of infinitely many coupled differential equations -- are
rare, and are valuable for numerical benchmarking and getting intuition for the
behavior of more complicated systems. Jahnke and Huisinga's landmark
calculation of the exact time-dependent solution of the CME for monomolecular
reaction systems is one of the most general analytic results known; however, it
is hard to generalize, because it relies crucially on properties of
monomolecular reactions. In this paper, we rederive Jahnke and Huisinga's
result on the time-dependent probability distribution and moments of
monomolecular reaction systems using the Doi-Peliti path integral approach,
which reduces solving the CME to evaluating many integrals. While the
Doi-Peliti approach is less intuitive, it is also more mechanical, and hence
easier to generalize. To illustrate how the Doi-Peliti approach can go beyond
the method of Jahnke and Huisinga, we also find an explicit and exact
time-dependent solution to a problem involving an autocatalytic reaction that
Jahnke and Huisinga identified as not solvable using their method. We also find
a formal exact time-dependent solution for any CME whose list of reactions
involves only zero and first order reactions, which may be the most general
result currently known.


  With recent advances in sequencing technologies, large amounts of epigenomic
data have become available and computational methods are contributing
significantly to the progress of epigenetic research. As an orthogonal approach
to methods based on machine learning, mechanistic modeling aims at a
description of the mechanisms underlying epigenetic changes. Here, we propose
an efficient method for parameter estimation for stochastic models that
describe the dynamics of DNA methylation patterns over time. Our method is
based on the Generalized Method of Moments (GMM) and gives results with an
accuracy similar to that of maximum likelihood-based estimation approaches.
However, in contrast to the latter, the GMM still allows an efficient and
accurate calibration of parameters even if the complexity of the model is
increased by considering longer methylation patterns. We show the usefulness of
our method by applying it to hairpin bisulfite sequencing data from mouse ESCs
for varying pattern lengths.


  The traditional focus of physiological and functional genomic research is on
molecular processes that play out within a single body. In contrast, when
social interactions occur, molecular and behavioral responses in interacting
individuals can lead to physiological processes that are distributed across
multiple individuals. In eusocial insect colonies, such multi-body processes
are tightly integrated, involving social communication mechanisms that regulate
the physiology of colony members. As a result, conserved physiological
mechanisms, for example related to pheromone detection and neural signaling
pathways, are deployed in novel contexts and regulate emergent colony traits
during the evolutionary origin and elaboration of social complexity. Here we
review conceptual frameworks for organismal and colony physiology, and
highlight functional genomic, physiological, and behavioral research exploring
how colony-level traits arise from physical and chemical interactions among
nestmates. We highlight mechanistic work exploring how colony traits arise from
physical and chemical interactions among physiologically-specialized nestmates
of various developmental stages. We consider similarities and differences
between organismal and colony physiology, and make specific predictions based
on a decentralized perspective on the function and evolution of colony traits.
Integrated models of colony physiological function will be useful to address
fundamental questions related to the evolution and ecology of collective
behavior in natural systems.


  A labeled gene tree topology that is more probable than the labeled gene tree
topology matching a species tree is called \textit{anomalous}. Species trees
that can generate such anomalous gene trees are said to be in the
\textit{anomaly zone}. Here, probabilities of \textit{unranked} and
\textit{ranked} gene tree topologies under the multispecies coalescent are
considered. A ranked tree depicts not only the topological relationship among
gene lineages, as an unranked tree does, but also the sequence in which the
lineages coalesce. In this article, we study how the parameters of a species
tree simulated under a constant rate birth-death process can affect the
probability that the species tree lies in the anomaly zone. We find that with
more than five taxa, it is possible for species trees have both AGTs and ARGTs.
The probability of being in either type of anomaly zones increases with more
taxa. The probability of AGTs also increases with higher speciation rates. We
observe that the probabilities of unranked anomaly zones are higher and grow
much faster than those of ranked anomaly zones as the speciation rate
increases. Our simulation shows that the most probable ranked gene tree is
likely to have the same unranked topology as the species tree. We design the
software {\it PRANC} which computes probabilities of ranked gene tree
topologies given a species tree under the coalescent model.


  This paper is dedicated to the study of the interaction between dynamical
systems and percolation models, with views towards the study of viral
infections whose virus mutate with time. Recall that r-bootstrap percolation
describes a deterministic process where vertices of a graph are infected once r
neighbors of it are infected. We generalize this by introducing F(t)-bootstrap
percolation, a time-dependent process where the number of neighbouring vertices
which need to be infected for a disease to be transmitted is determined by a
percolation function F(t) at each time t. After studying some of the basic
properties of the model, we consider smallest percolating sets and construct a
polynomial-timed algorithm to find one smallest minimal percolating set on
finite trees for certain F(t)-bootstrap percolation models.