  An important and well-studied problem in phylogenetics is to compute a
\emph{consensus tree} so as to summarize the common features within a
collection of rooted phylogenetic trees, all whose leaf-sets are bijectively
labeled by the same set~(X) of species. More recently, however, it has become
of interest to find a consensus for a collection of more general, rooted
directed acyclic graphs all of whose sink-sets are bijectively labeled by~(X),
so called rooted \emph{phylogenetic networks}. These networks are used to
analyse the evolution of species that cross with one another, such as plants
and viruses. In this paper, we introduce an algorithm for computing a consensus
for a collection of so-called 1-\emph{nested} phylogenetic networks. Our
approach builds on a previous result by Rosell\'o et al. that describes an
encoding for any 1-nested phylogenetic network in terms of a collection of
ordered pairs of subsets of (X).More specifically, we characterize those
collections of ordered pairs that arise as the encoding of some 1-nested
phylogenetic network, and then use this characterization to compute a
\emph{consensus network} for a collection of~$t$ 1-nested networks in
$O(t|X|^2+|X|^3)$ time. Applying our algorithm to a collection of phylogenetic
trees yields the well-known majority rule consensus tree. Our approach leads to
several new directions for futurework, and we expect that it should provide a
useful new tool to help understand complex evolutionary scenarios.


  Deterministic reaction networks (RNs) are tools to model diverse biological
phenomena characterized by particle systems, when there are abundant number of
particles. Examples include but are not limited to biochemistry, molecular
biology, genetics, epidemiology, and social sciences. In this chapter we
propose a new type of decomposition of RNs, called fiber decomposition. Using
this decomposition, we establish lifting of mass-action RNs preserving
stationary properties, including multistationarity and absolute concentration
robustness. Such lifting scheme is simple and explicit which imposes little
restriction on the reaction networks. We provide examples to illustrate how
this lifting can be used to construct RNs preserving certain dynamical
properties.


  The evolution and function of imitation have always been placed within the
confines of animal learning and associated with its crucial role in cultural
transmission and cultural evolution. Can imitation evolve as a form of
phenotypic plasticity in the absence of cultural transmission, in phenotypes
beyond behavior?
  We investigate a model in which imitation is unbundled from cultural
transmission: an organism's adult phenotype is plastically altered by its
experiences as a juvenile of other juveniles' genetically determined traits.
The only information transmitted between generations is genetic.
  We find that during a period of directional selection towards a phenotypic
optimum, natural selection favors modifiers which cause an organism to bias its
plastic phenotype in the direction opposite to the mean phenotype of the
population -- anti-imitation. As the population approaches the phenotypic
optimum and shifts into stabilizing selection, selection on the modifier
reverses and favors strong imitation of the population mean. Imitation can
evolve to overshoot the target and produce an evolutionary pathology where mean
fitness decreases. When purifying selection for an extreme phenotype is
modeled, only selection for anti-imitation occurs, even at a mutation-selection
balance.
  Imitation and anti-imitation emerge from these models in the absence of any
goal representation, cognitive understanding of its purpose, or discernment of
any kind. These theoretical outcomes are all novel evolutionary and biological
phenomena, and we discuss their implications.


  The emergence of novel pathogens and zoonotic diseases like the SARS-CoV-2
have underlined the need for developing novel diagnosis and intervention
pipelines that can learn rapidly from small amounts of labeled data. Combined
with technological advances in next-generation sequencing, metagenome-based
diagnostic tools hold much promise to revolutionize rapid point-of-care
diagnosis. However, there are significant challenges in developing such an
approach, the chief among which is to learn self-supervised representations
that can help detect novel pathogen signatures with very low amounts of labeled
data. This is particularly a difficult task given that closely related
pathogens can share more than 90% of their genome structure. In this work, we
address these challenges by proposing MG-Net, a self-supervised representation
learning framework that leverages multi-modal context using pseudo-imaging data
derived from clinical metagenome sequences. We show that the proposed framework
can learn robust representations from unlabeled data that can be used for
downstream tasks such as metagenome sequence classification with limited access
to labeled data. Extensive experiments show that the learned features
outperform current baseline metagenome representations, given only 1000 samples
per class.


  Forecasting new cases, hospitalizations, and disease-induced deaths is an
important part of infectious disease surveillance and helps guide health
officials in implementing effective countermeasures. For disease surveillance
in the U.S., the Centers for Disease Control and Prevention (CDC) combine more
than 65 individual forecasts of these numbers in an ensemble forecast at
national and state levels. We collected data on CDC ensemble forecasts of
COVID-19 fatalities in the United States, and compare them with easily
interpretable ``Euler'' forecasts serving as a model-free benchmark that is
only based on the local rate of change of the incidence curve. The term ``Euler
method'' is motivated by the eponymous numerical integration scheme that
calculates the value of a function at a future time step based on the current
rate of change. Our results show that CDC ensemble forecasts are not more
accurate than ``Euler'' forecasts on short-term forecasting horizons of one
week. However, CDC ensemble forecasts show a better performance on longer
forecasting horizons. Using the current rate of change in incidences as
estimates of future incidence changes is useful for epidemic forecasting on
short time horizons. An advantage of the proposed method over other forecasting
approaches is that it can be implemented with a very limited amount of work and
without relying on additional data (e.g., human mobility and contact patterns)
and high-performance computing systems.


  Applying a ML approach to the temporal variability of the Spike protein
sequence enables us to identify, classify and track emerging virus variants.
Our analysis is unbiased, in the sense that it does not require any prior
knowledge of the variant characteristics, and our results are validated by
other informed methods that define variants based on the complete genome.
Furthermore, correlating persistent variants of our approach to epidemiological
data, we discover that each new wave of the COVID-19 pandemic is driven and
dominated by a new emerging variant. Our results are therefore indispensable
for further studies on the evolution of SARS-CoV-2 and the prediction of
evolutionary patterns that determine current and future mutations of the Spike
proteins, as well as their diversification and persistence during the viral
spread. Moreover, our ML algorithm works as an efficient early warning system
for the emergence of new persistent variants that may pose a threat of
triggering a new wave of COVID-19. Capable of a timely identification of
potential new epidemiological threats when the variant only represents 1% of
the new sequences, our ML strategy is a crucial tool for decision makers to
define short and long term strategies to curb future outbreaks. The same
methodology can be applied to other viral diseases, influenza included, if
sufficient sequencing data is available.


  The spread of COVID-19 has been greatly impacted by regulatory policies and
behavior patterns that vary across counties, states, and countries.
Population-level dynamics of COVID-19 can generally be described using a set of
ordinary differential equations, but these deterministic equations are
insufficient for modeling the observed case rates, which can vary due to local
testing and case reporting policies and non-homogeneous behavior among
individuals. To assess the impact of population mobility on the spread of
COVID-19, we have developed a novel Bayesian time-varying coefficient
state-space model for infectious disease transmission. The foundation of this
model is a time-varying coefficient compartment model to recapitulate the
dynamics among susceptible, exposed, undetected infectious, detected
infectious, undetected removed, detected non-infectious, detected recovered,
and detected deceased individuals. The infectiousness and detection parameters
are modeled to vary by time, and the infectiousness component in the model
incorporates information on multiple sources of population mobility. Along with
this compartment model, a multiplicative process model is introduced to allow
for deviation from the deterministic dynamics. We apply this model to observed
COVID-19 cases and deaths in several US states and Colorado counties. We find
that population mobility measures are highly correlated with transmission rates
and can explain complicated temporal variation in infectiousness in these
regions. Additionally, the inferred connections between mobility and
epidemiological parameters, varying across locations, have revealed the
heterogeneous effects of different policies on the dynamics of COVID-19.


  We explore the commonalities between methods for assuring the security of
computer systems (cybersecurity) and the mechanisms that have evolved through
natural selection to protect vertebrates against pathogens, and how insights
derived from studying the evolution of natural defenses can inform the design
of more effective cybersecurity systems. More generally, security challenges
are crucial for the maintenance of a wide range of complex adaptive systems,
including financial systems, and again lessons learned from the study of the
evolution of natural defenses can provide guidance for the protection of such
systems.


  The acquisition of Antimicrobial Multidrug Resistance (AMR) in patients
admitted to the Intensive Care Units (ICU) is a major global concern. This
study analyses data in the form of multivariate time series (MTS) from 3476
patients recorded at the ICU of University Hospital of Fuenlabrada (Madrid)
from 2004 to 2020. 18\% of the patients acquired AMR during their stay in the
ICU. The goal of this paper is an early prediction of the development of AMR.
Towards that end, we leverage the time-series cluster kernel (TCK) to learn
similarities between MTS. To evaluate the effectiveness of TCK as a kernel, we
applied several dimensionality reduction techniques for visualization and
classification tasks. The experimental results show that TCK allows identifying
a group of patients that acquire the AMR during the first 48 hours of their ICU
stay, and it also provides good classification capabilities.


  The rich and varied ways that genetic material can be passed between species
has motivated extensive research into the theory of phylogenetic networks.
Features that align with biological processes, or with mathematical
tractability, have been used to define classes and prove results, with the
ultimate goal of developing the theoretical foundations for network
reconstruction methods. We have now reached the point that a collection of
recent results can be drawn together to make one class of network, the
so-called normal networks, a clear leader. We may be at the cusp of practical
inference with phylogenetic networks.


  In this paper, the SIR dynamics is extended by considering another
compartmental which represents hospitalization of the critical cases. So a
system of differential equations with four blocks is considered when there is
intensive care unit (ICU) to cure critical cases. Outgoing rate of survived
infected individuals is divided into $nI$ and $\frac{bI}{I+b}$. The second term
represents the rate of critical cases who enter ICUs. It is proved that there
are forward, backward and Hopf bifurcations in different regimes of parameters.


  Let $\mathcal{T}$ be an unrooted binary tree with $n$ distinctly labelled
leaves. Deriving its name from the field of phylogenetics, a convex character
on $\mathcal{T}$ is simply a partition of the leaves such that the minimal
spanning subtrees induced by the blocks of the partition are mutually disjoint.
In earlier work Kelk and Stamoulis (Advances in Applied Mathematics 84 (2017),
pp. 34--46) defined $g_k(\mathcal{T})$ as the number of convex characters where
each block has at least $k$ leaves. Exact expressions were given for $g_1$ and
$g_2$, where the topology of $\mathcal{T}$ turns out to be irrelevant, and it
was noted that for $k \geq 3$ topological neutrality no longer holds. In this
article, for every $k \geq 3$ we describe tree topologies achieving the maximum
and minimum values of $g_k$ and determine corresponding expressions and
exponential bounds for $g_k$. Finally, we reflect briefly on possible
algorithmic applications of these results.


  Genomic selection (GS) is a technique that plant breeders use to select
individuals to mate and produce new generations of species. Allocation of
resources is a key factor in GS. At each selection cycle, breeders are facing
the choice of budget allocation to make crosses and produce the next generation
of breeding parents. Inspired by recent advances in reinforcement learning for
AI problems, we develop a reinforcement learning-based algorithm to
automatically learn to allocate limited resources across different generations
of breeding. We mathematically formulate the problem in the framework of Markov
Decision Process (MDP) by defining state and action spaces. To avoid the
explosion of the state space, an integer linear program is proposed that
quantifies the trade-off between resources and time. Finally, we propose a
value function approximation method to estimate the action-value function and
then develop a greedy policy improvement technique to find the optimal
resources. We demonstrate the effectiveness of the proposed method in enhancing
genetic gain using a case study with realistic data.


  The spread of an infectious disease depends on intrinsic properties of the
disease as well as the connectivity and actions of the population. This study
investigates the dynamics of an SIR type model which accounts for human
tendency to avoid infection while also maintaining preexisting, interpersonal
relationships. Specifically, we use a network model in which individuals
probabilistically deactivate connections to infected individuals and later
reconnect to the same individuals upon recovery. To analyze this network model,
a mean field approximation consisting of a system of fourteen ordinary
differential equations for the number of nodes and edges is developed. This
system of equations is closed using a moment closure approximation for the
number of triple links. By analyzing the differential equations, it is shown
that, in addition to force of infection and recovery rate, the probability of
deactivating edges and the average node degree of the underlying network
determine if an epidemic occurs.


  Fast Repetition Rate fluorometry (FRRf) generates high-resolution measures of
phytoplankton primary productivity as electron transport rates (ETRs). How ETRs
scale to corresponding inorganic carbon (C) uptake rates (the so-called
electron requirement for carbon fixation, e,C), inherently describes the extent
and effectiveness with which absorbed light energy drives C-fixation. However,
it remains unclear whether and how e,C follows predictable patterns for
oceanographic datasets spanning physically dynamic, and complex, environmental
gradients. We utilise a unique high-throughput approach, coupling ETRs and
14C-incubations to produce a semi-continuous dataset of e,C (n = 80),
predominantly from surface waters, along the Australian coast (Brisbane to the
Tasman Sea), including the East Australian Current (EAC). Environmental
conditions along this transect could be generally grouped into cooler, more
nutrient-rich waters dominated by larger size-fractionated Chl-a (>10 um)
versus warmer nutrient-poorer waters dominated by smaller size-fractionated
Chl-a (< 2 um). Whilst e,C was higher for warmer water samples, environmental
conditions alone explained less than 20% variance of e,C, and changes in
predominant size-fraction(s) distributions of Chl-a (biomass) failed to explain
variance of e,C. Instead, NPQNSV was a better predictor of e,C, explaining 55%
of observed variability. NPQNSV is a physiological descriptor that accounts for
changes in both long-term driven acclimation in non-radiative decay, and
quasi-instantaneous PSII downregulation, and thus may prove a useful predictor
of e,C across physically-dynamic regimes, provided the slope describing their
relationship is predictable.


  Nutritional regulation by ants emerges from a distributed process: food is
collected by a small fraction of workers, stored within the crops of
individuals, and spreads via local ant-to-ant interactions. The precise
individual level underpinnings of this collective regulation have remained
unclear mainly due to difficulties in measuring food within ants crops. Here we
image fluorescent liquid food in individually tagged Camponotus sanctus ants,
and track the real-time food flow from foragers to their gradually satiating
colonies. We show how the feedback between colony satiation level and food
inflow is mediated by individual crop loads; specifically, the crop loads of
recipient ants control food flow rates, while those of foragers regulate the
frequency of foraging trips. Interestingly, these effects do not rise from pure
physical limitations of crop capacity. Our findings suggest that the emergence
of food intake regulation does not require individual foragers to assess the
global state of the colony.


  We present several algorithms designed to learn a pattern of correspondence
between two data sets in situations where it is desirable to match elements
that exhibit a relationship belonging to a known parametric model. In the
motivating case study, the challenge is to better understand micro-RNA
regulation in the striatum of Huntington's disease model mice. The algorithms
unfold in two stages. First, an optimal transport plan P and an optimal affine
transformation are learned, using the Sinkhorn-Knopp algorithm and a mini-batch
gradient descent. Second, P is exploited to derive either several co-clusters
or several sets of matched elements. A simulation study illustrates how the
algorithms work and perform. The real data application further illustrates
their applicability and interest.


  In this paper we consider two aspects of the inverse problem of how to
construct merge trees realizing a given barcode. Much of our investigation
exploits a recently discovered connection between the symmetric group and
barcodes in general position, based on the simple observation that death order
is a permutation of birth order. The first important outcome of our study is a
clear combinatorial distinction between the space of phylogenetic trees (as
defined by Billera, Holmes and Vogtmann) and the space of merge trees. Generic
BHV trees on $n+1$ leaf nodes fall into $(2n-1)!!$ distinct strata, but the
analogous number for merge trees is equal to the number of maximal chains in
the lattice of partitions, i.e., $(n+1)!n!2^{-n}$. The second aspect of our
study is the derivation of precise formulas for the distribution of tree
realization numbers (the number of merge trees realizing a given barcode) when
we assume that barcodes are sampled using a uniform distribution on the
symmetric group. We are able to characterize some of the higher moments of this
distribution, thanks in part to a reformulation in terms of Dirichlet
convolution. This characterization provides a type of null hypothesis,
apparently different from the distributions observed in real neuron data and
opens the door to doing more precise science.


  We investigate a stochastic individual-based model for the population
dynamics of host-virus systems where the microbial hosts may transition into a
dormant state upon contact with virions, thus evading infection. Such a
contact-mediated defence mechanism was described in Bautista et al (2015) for
an archaeal host, while Jackson and Fineran (2019) and Meeske et al (2019)
describe a related, CRISPR-Cas induced, dormancy defense of bacterial hosts to
curb phage epidemics. We first analyse the effect of the dormancy-related model
parameters on the probability and time of invasion of a newly arriving virus
into a resident host population. Given successful invasion in the stochastic
system, we then show that the emergence (with high probability) of a persistent
virus infection ('epidemic') in a large host population can be determined by
the existence of a coexistence equilibrium for the dynamical system arising as
the deterministic many-particle limit of our model. This is an extension of a
dynamical system considered by Beretta and Kuang (1998) that is known to
exhibit a Hopf bifurcation, giving rise to a 'paradox of enrichment'. In our
system, we verify that the additional dormancy component can, at least for
certain parameter ranges, prevent the associated loss of stability. Finally, we
show that the presence of contact-mediated dormancy enables the host population
to attain higher equilibrium sizes - while still being able to avoid a
persistent epidemic - than host populations without this trait.


  Motivation: MicroRNAs (miRNAs) play pivotal roles in gene expression
regulation by binding to target sites of messenger RNAs (mRNAs). While
identifying functional targets of miRNAs is of utmost importance, their
prediction remains a great challenge. Previous computational algorithms have
major limitations. They use conservative candidate target site (CTS) selection
criteria mainly focusing on canonical site types, rely on laborious and
time-consuming manual feature extraction, and do not fully capitalize on the
information underlying miRNA-CTS interactions. Results: In this paper, we
introduce TargetNet, a novel deep learning-based algorithm for functional miRNA
target prediction. To address the limitations of previous approaches, TargetNet
has three key components: (1) relaxed CTS selection criteria accommodating
irregularities in the seed region, (2) a novel miRNA-CTS sequence encoding
scheme incorporating extended seed region alignments, and (3) a deep residual
network-based prediction model. The proposed model was trained with miRNA-CTS
pair datasets and evaluated with miRNA-mRNA pair datasets. TargetNet advances
the previous state-of-the-art algorithms used in functional miRNA target
classification. Furthermore, it demonstrates great potential for distinguishing
high-functional miRNA targets.


  Gene gain-loss-duplication models are commonly based on continuous-time
birth-death processes. Employed in a phylogenetic context, such models have
been increasingly popular in studies of gene content evolution across multiple
genomes. While the applications are becoming more varied and demanding,
bioinformatics methods for probabilistic inference on copy numbers (or
integer-valued evolutionary characters, in general) are scarce. We describe a
flexible probabilistic framework for phylogenetic gene-loss-duplication models.
The framework is based on a novel elementary representation by dependent random
variables with well-characterized conditional distributions: binomial, P\'olya
(negative binomial), and Poisson. The corresponding graphical model yields
exact numerical procedures for computing the likelihood and the posterior
distribution of ancestral copy numbers. The resulting algorithms take quadratic
time in the total number of copies. In addition, we show how the likelihood
gradient can be computed by a linear-time algorithm.


  Reference population databases are an essential tool in variant and gene
interpretation. Their use guides the identification of pathogenic variants
amidst the sea of benign variation present in every human genome, and supports
the discovery of new disease-gene relationships. The Genome Aggregation
Database (gnomAD) is currently the largest and most widely used publicly
available collection of population variation from harmonized sequencing data.
The data is available through the online gnomAD browser
(https://gnomad.broadinstitute.org/) that enables rapid and intuitive variant
analysis. This review provides guidance on the content of the gnomAD browser,
and its usage for variant and gene interpretation. We introduce key features
including allele frequency, per-base expression levels, constraint scores, and
variant co-occurrence, alongside guidance on how to use these in analysis, with
a focus on the interpretation of candidate variants and novel genes in rare
disease.


  Trajectories of endosomes inside living eukaryotic cells are highly
heterogeneous in space and time and diffuse anomalously due to a combination of
viscoelasticity, caging, aggregation and active transport. Some of the
trajectories display switching between persistent and anti-persistent motion
while others jiggle around in one position for the whole measurement time. By
splitting the ensemble of endosome trajectories into slow moving sub-diffusive
and fast moving super-diffusive endosomes, we analyzed them separately. The
mean squared displacements and velocity auto-correlation functions confirm the
effectiveness of the splitting methods. Applying the local analysis, we show
that both ensembles are characterized by a spectrum of local anomalous
exponents and local generalized diffusion coefficients. Slow and fast endsomes
have exponential distributions of local anomalous exponents and power law
distributions of generalized diffusion coefficients. This suggests that
heterogeneous fractional Brownian motion is an appropriate model for both fast
and slow moving endosomes. This article is part of a Special Issue entitled:
"Recent Advances In Single-Particle Tracking: Experiment and Analysis" edited
by Janusz Szwabi\'nski and Aleksander Weron.


  In contexts ranging from embryonic development to bacterial ecology, cell
populations migrate chemotactically along self-generated chemical gradients,
often forming a propagating front. Here, we theoretically show that the
stability of such chemotactic fronts to morphological perturbations is
determined by limitations in the ability of individual cells to sense and
thereby respond to the chemical gradient. Specifically, cells at bulging parts
of a front are exposed to a smaller gradient, which slows them down and
promotes stability, but they also respond more strongly to the gradient, which
speeds them up and promotes instability. We predict that this competition leads
to chemotactic fingering when sensing is limited at too low chemical
concentrations. Guided by this finding and by experimental data on E. coli
chemotaxis, we suggest that the cells' sensory machinery might have evolved to
avoid these limitations and ensure stable front propagation. Finally, as
sensing of any stimuli is necessarily limited in living and active matter in
general, the principle of sensing-induced stability may operate in other types
of directed migration such as durotaxis, electrotaxis, and phototaxis.


  Cells sense their surrounding by employing intracellular signaling pathways
that transmit hormonal signals from the cell membrane to the nucleus.
TGF-$\beta$/SMAD signaling encodes various cell fates, controls tissue
homeostasis and is deregulated in diseases such as cancer. The pathway shows
strong heterogeneity at the single-cell level, but quantitative insights into
mechanisms underlying fluctuations at various time scales are still missing,
partly due to inefficiency in the calibration of stochastic models that
mechanistically describe signaling processes. In this work we analyze
single-cell TGF-$\beta$/SMAD signaling and show that it exhibits temporal
stochastic bursts which are dose-dependent and whose number and magnitude
correlate with cell migration. We propose a stochastic modeling approach to
mechanistically describe these pathway fluctuations with high computational
efficiency. Employing high-order numerical integration and fitting to burst
statistics we enable efficient quantitative parameter estimation and
discriminate models that assume noise in different reactions at the receptor
level. This modeling approach suggests that stochasticity in the
internalization of TGF-$\beta$ receptors into endosomes plays a key role in the
observed temporal bursting. Further, the model predicts the single-cell
dynamics of TGF-$\beta$/SMAD signaling in untested conditions, e.g.,
successfully reflects memory effects of signaling noise and cellular
sensitivity towards repeated stimulation. Taken together, our computational
framework based on burst analysis, noise modeling and path computation scheme
is a suitable tool for the data-based modeling of complex signaling pathways,
capable of identifying the source of temporal noise.


  Proteins tend to bury hydrophobic residues inside their core during the
folding process to provide stability to the protein structure and to prevent
aggregation. Nevertheless, proteins do expose some 'sticky' hydrophobic
residues to the solvent. These residues can play an important functional role,
for example in protein-protein and membrane interactions. Here, we investigate
how hydrophobic protein surfaces are by providing three measures for surface
hydrophobicity: the total hydrophobic surface area, the relative hydrophobic
surface area, and - using our MolPatch method - the largest hydrophobic patch.
Secondly, we analyse how difficult it is to predict these measures from
sequence: by adapting solvent accessibility predictions from NetSurfP2.0, we
obtain well-performing prediction methods for the THSA and RHSA, while
predicting LHP is more difficult. Finally, we analyse implications of exposed
hydrophobic surfaces: we show that hydrophobic proteins typically have low
expression, suggesting cells avoid an overabundance of sticky proteins.


  We apply optimal control theory to a generalized SEIR-type model. The
proposed system has three controls, representing social distancing, preventive
means, and treatment measures to combat the spread of the COVID-19 pandemic. We
analyze such optimal control problem with respect to real data transmission in
Italy. Our results show the appropriateness of the model, in particular with
respect to the number of quarantined/hospitalized (confirmed and infected) and
recovered individuals. Considering the Pontryagin controls, we show how in a
perfect world one could have drastically diminish the number of susceptible,
exposed, infected, quarantined/hospitalized, and death individuals, by
increasing the population of insusceptible/protected.


  Combining different modalities of data from human tissues has been critical
in advancing biomedical research and personalised medical care. In this study,
we leverage a graph embedding model (i.e VGAE) to perform link prediction on
tissue-specific Gene-Gene Interaction (GGI) networks. Through ablation
experiments, we prove that the combination of multiple biological modalities
(i.e multi-omics) leads to powerful embeddings and better link prediction
performances. Our evaluation shows that the integration of gene methylation
profiles and RNA-sequencing data significantly improves the link prediction
performance. Overall, the combination of RNA-sequencing and gene methylation
data leads to a link prediction accuracy of 71% on GGI networks. By harnessing
graph representation learning on multi-omics data, our work brings novel
insights to the current literature on multi-omics integration in
bioinformatics.


  In order to target threatening pathogens, the adaptive immune system performs
a continuous reorganization of its lymphocyte repertoire. Following an immune
challenge, the B cell repertoire can evolve cells of increased specificity for
the encountered strain. This process of affinity maturation generates a memory
pool whose diversity and size remain difficult to predict. We assume that the
immune system follows a strategy that maximizes the long-term immune coverage
and minimizes the short-term metabolic costs associated with affinity
maturation. This strategy is defined as an optimal decision process on a finite
dimensional phenotypic space, where a pre-existing population of naive cells is
sequentially challenged with a neutrally evolving strain. We unveil a trade-off
between immune protection against future strains and the necessary
reorganization of the repertoire. This plasticity of the repertoire drives the
emergence of distinct regimes for the size and diversity of the memory pool,
depending on the density of naive cells and on the mutation rate of the strain.
The model predicts power-law distributions of clonotype sizes observed in data,
and rationalizes antigenic imprinting as a strategy to minimize metabolic costs
while keeping good immune protection against future strains.


  The spread of COVID-19 has been thwarted in most countries through
non-pharmaceutical interventions. In particular, the most effective measures in
this direction have been the stay-at-home and closure strategies of businesses
and schools. However, population-wide lockdowns are far from being optimal
carrying heavy economic consequences. Therefore, there is nowadays a strong
interest in designing more efficient restrictions. In this work, starting from
a recent kinetic-type model which takes into account the heterogeneity
described by the social contact of individuals, we analyze the effects of
introducing an optimal control strategy into the system, to limit selectively
the mean number of contacts and reduce consequently the number of infected
cases. Thanks to a data-driven approach, we show that this new mathematical
model permits to assess the effects of the social limitations. Finally, using
the model introduced here and starting from the available data, we show the
effectivity of the proposed selective measures to dampen the epidemic trends.


  Protein-RNA interactions are of vital importance to a variety of cellular
activities. Both experimental and computational techniques have been developed
to study the interactions. Due to the limitation of the previous database,
especially the lack of protein structure data, most of the existing
computational methods rely heavily on the sequence data, with only a small
portion of the methods utilizing the structural information. Recently,
AlphaFold has revolutionized the entire protein and biology field. Foreseeably,
the protein-RNA interaction prediction will also be promoted significantly in
the upcoming years. In this work, we give a thorough review of this field,
surveying both the binding site and binding preference prediction problems and
covering the commonly used datasets, features, and models. We also point out
the potential challenges and opportunities in this field. This survey
summarizes the development of the RBP-RNA interaction field in the past and
foresees its future development in the post-AlphaFold era.


  A delayed, discrete-time, prey-predator model with Allee effects imposed on
prey and predator populations is defined, and dynamics of the system is
characterized computationally. The parametric conditions for local asymptotic
stability of equilibria are investigated with several illustrating examples. In
addition, several periodic solutions were achieved. Furthermore, high degree of
fractality of the prey-predator population trajectories is observed when the
system is sufficiently delayed. This system has turned out to be a cooperative
system, which is depicted using fractal dimension.


  The paper is an introduction to the concept of vaccine efficacy and its
practical applications. We provide illustrative examples aimed at a wider
audience.


  Research shows that gene duplication followed by either repurposing or
removal of duplicated genes is an important contributor to evolution of gene
and protein interaction networks. We aim to identify which characteristics of a
network can arise through this process, and which must have been produced in a
different way. To model the network evolution, we postulate vertex duplication
and edge deletion as evolutionary operations on graphs. Using the novel concept
of an ancestrally distinguished subgraph, we show how features of present-day
networks require certain features of their ancestors. In particular,
ancestrally distinguished subgraphs cannot be introduced by vertex duplication.
Additionally, if vertex duplication and edge deletion are the only evolutionary
mechanisms, then a graph's ancestrally distinguished subgraphs must be
contained in all of the graph's ancestors. We analyze two experimentally
derived genetic networks and show that our results accurately predict lack of
large ancestrally distinguished subgraphs, despite this feature being
statistically improbable in associated random networks. This observation is
consistent with the hypothesis that these networks evolved primarily via vertex
duplication. The tools we provide open the door for analysing ancestral
networks using current networks. Our results apply to edge-labeled (e.g.
signed) graphs which are either undirected or directed.


  Geometric deep learning (GDL), which is based on neural network architectures
that incorporate and process symmetry information, has emerged as a recent
paradigm in artificial intelligence. GDL bears particular promise in molecular
modeling applications, in which various molecular representations with
different symmetry properties and levels of abstraction exist. This review
provides a structured and harmonized overview of molecular GDL, highlighting
its applications in drug discovery, chemical synthesis prediction, and quantum
chemistry. Emphasis is placed on the relevance of the learned molecular
features and their complementarity to well-established molecular descriptors.
This review provides an overview of current challenges and opportunities, and
presents a forecast of the future of GDL for molecular sciences.


  The COVID-19 pandemic is the most significant global crisis since World War
II that affected almost all the countries of our planet. To control the
COVID-19 pandemic outbreak, it is necessary to understand how the virus is
transmitted to a susceptible individual and eventually spread in the community.
The primary transmission pathway of COVID-19 is human-to-human transmission
through infectious droplets. However, a recent study by Greenhalgh et al.
(Lancet: 397:1603-1605, 2021) demonstrates 10 scientific reasons behind the
airborne transmission of SARS-COV-2. In the present study, we introduce a novel
mathematical model of COVID-19 that considers the transmission of free viruses
in the air besides the transmission of direct contact with an infected person.
The basic reproduction number of the epidemic model is calculated using the
next-generation operator method and observed that it depends on both the
transmission rate of direct contact and free virus contact. The local and
global stability of disease-free equilibrium (DFE) is well established.
Analytically it is found that there is a forward bifurcation between the DFE
and an endemic equilibrium using central manifold theory. Next, we used the
nonlinear least-squares technique to identify the best-fitted parameter values
in the model from the observed COVID-19 mortality data of two major districts
of India. Using estimated parameters for Bangalore urban and Chennai, different
control scenarios for mitigation of the disease are investigated. Results
indicate that when a vaccine crisis is there, the public health authorities
should prefer to vaccinate the susceptible people compared to the recovered
persons who are now healthy. Along with face mask use, treatment of
hospitalized patients and vaccination of susceptibles, immigration should be
allowed in a supervised manner so that economy of the overall society remains
healthy.


  Furin is a serine protease that takes part in the processing and activation
of the host cell pre-proteins. The enzyme also plays an important role in the
activation of several viruses like the newly emerging SARS-CoV-2 virus that
causes COVID-19 disease with a high rate of virulence and mortality. Unlike
viral enzymes, furin owns a constant sequence and active site characteristics
and seems to be a better target for drug design for COVID-19 treatment.
Considering furin active site as receptor and some approved drugs from
different classes including antiviral, antibiotics, and anti protozoa/anti
parasites with suspected beneficial effects on COVID-19, as ligands we have
carried out docking experiments in HEX software to pickup those capable to bind
furin active site with high affinity and suggest them as probable candidates
for clinical trials assessments. Our docking experiments show that saquinavir,
nelfinavir, and atazanavir with cumulative inhibitory effects of 2.52, 2.16,
and 2.13 respectively seem to be the best candidates for furin inhibition even
in severe cases of COVID-19 as adjuvant therapy, while clarithromycin,
niclosamide, and erythromycin with cumulative inhibitory indices of 1.97, 1.90,
and 1.84 respectively with lower side effects than antiviral drugs could be
suggested as prophylaxes for the first stage of COVID-19 as a promising treat.


  Bam-readcount is a utility for generating low-level information about
sequencing data at specific nucleotide positions. Originally designed to help
filter genomic mutation calls, the metrics it outputs are useful as input for
variant detection tools and for resolving ambiguity between variant callers .
In addition, it has found broad applicability in diverse fields including tumor
evolution, single-cell genomics, climate change ecology, and tracking community
spread of SARS-CoV-2. Here we report on the release of version 1.0 of this
tool, which adds CRAM support, among other improvements. It is released under a
permissive MIT license and available at
https://github.com/genome/bam-readcount.


  The emerging field of Molecular Communication (MC) aims to characterize
biological-based signaling environment through information that are encoded
into molecules. Since the birth of this field, a number of different
applications and biological systems have been characterized using MC theory.
This study proposes a new application and direction for MC, focusing on the
digestive system, where we characterize and model the starch and glucose
propagation along the small intestine. Based on the advection-diffusion and
reaction mechanisms, we define a channel capacity for the small intestine
digestive tract that is dependent on the starch to glucose conversion, velocity
flow within the tract, viscosity of the digest product, and length of the tract
and position of the receivers for glucose absorption. The numerical results
from the derived channel capacity model shows that the SI digestive capacity
depends both on physiological factors of the digestive system and type of
consumed food, where the digestive capacity is greater for shorter gastric
emptying time, low viscosity of the digest product and efficient enzyme
activity. We believe that our digital MC model of the digestive tract and lead
to personalized diet for each individual, which can potentially avoid a number
of different diseases (e.g., celiac disease).


  In this paper, we conduct theoretical analyses on inferring the structure of
gene regulatory networks. Depending on the experimental method and data type,
the inference problem is classified into 20 different scenarios. For each
scenario, we discuss the problem that with enough data, under what assumptions,
what can be inferred about the structure. For scenarios that have been covered
in the literature, we provide a brief review. For scenarios that have not been
covered in literature, if the structure can be inferred, we propose new
mathematical inference methods and evaluate them on simulated data. Otherwise,
we prove that the structure cannot be inferred.


  Recently a technique based on the interaction between adhesion proteins
extracted from Streptococcus pyogenes, known as SpyRing, has been widely used
to improve the thermal resilience of enzymes, the assembly of biostructures,
cancer cell recognition and other fields. In SpyRing, the two termini of the
target enzyme are respectively linked to the peptide SpyTag and its protein
partner SpyCatcher. SpyTag spontaneously reacts with SpyCatcher to form an
isopeptide bond, with which the target enzyme forms a close ring structure. It
was believed that the covalent cyclization of protein skeleton caused by
SpyRing reduces the conformational entropy of biological structure and improves
its rigidity, thus improving the thermal resilience of the target enzyme.
However, the effects of SpyTag/ SpyCatcher interaction with this enzyme are
poorly understood, and their regulation of enzyme properties remains unclear.
Here, for simplicity, we took the single domain enzyme lichenase from Bacillus
subtilis 168 as an example, studied the interface interactions in the SpyRing
system by molecular dynamics simulations, and examined the effects of the
changes of electrostatic interaction and van der Waals interaction on the
thermal resilience of target enzyme. The simulations showed that the interface
between SpyTag/SpyCatcher and lichenase is different from that found by
geometric matching method and highlighted key mutations that affect the
intensity of interactions at the interface and might have effect on the thermal
resilience of the enzyme. Our calculations provided new insights into the
rational designs in the SpyRing.


  The physical interpretation of the functioning of the adaptive immune system,
which has been thoroughly characterized on genetic and molecular levels,
provides a unique opportunity to define an adaptive self-organizing biological
system in its entirety. This paper describes a configuration space model of
immune function, where directed chemical potentials of the system constitute a
space of interactions. In the physical sense, the humoral adaptive immune
system adjusts the chemical potential of all available antigenic molecules by
tuning the chemical potential and organizing the network hierarchy of its
sensor-effector molecules, antibodies. The coupling of sensors and effectors
allows the system to adjust the thermodynamic activity of antigens and
antibodies, while network organization helps minimize chemical potentials and
maximize diversity. Mathematically the system couples the variance of Gaussian
distributed interaction energies in its interaction space to the exponentially
distributed chemical potentials of its effector molecules to maintain its
stationary state. This process creates a scale-free network in interaction
space, where absolute thermodynamic activity corresponds to node degree. In the
thermodynamic interpretation, the system is an ensemble carrying out {mu}N
work, adjusting chemical potentials according to the changes in the chemical
potentials of the surroundings. The validity of the model is supported by
identifying an interaction flexibility index, the corresponding variables in
thermodynamics and network science, and by confirming its applicability to the
humoral immune system. Overall, this statistical thermodynamics model of
adaptive immunity describes how adaptive biological self-organization arises
from the maintenance of a scale-free, directed interaction network with fractal
topology.


  The stationary asymptotic properties of the diffusion limit of a multi-type
branching process with neutral mutations are studied. For the critical and
subcritical processes the interesting limits are those of quasi-stationary
distributions conditioned on non-extinction. Pedagogical derivations are given
for known results that the limiting distributions for supercritical and
critical processes are found to collapse onto rays aligned with stationary
eigenvectors of the mutation rate matrix, in agreement with discrete multi-type
branching processes. For the sub-critical process the previously unsolved
quasi-stationary distribution is obtained to first order in the overall
mutation rate, which is assumed to be small. The sampling distribution over
allele types for a sample of given finite size is found to agree to first order
in mutation rates with the analogous sampling distribution for a Wright-Fisher
diffusion with constant population size.


  Genome-scale screening experiments in cancer produce long lists of candidate
genes that require extensive interpretation for biological insight and
prioritization for follow-up studies. Interrogation of gene lists frequently
represents a significant and time-consuming undertaking, in which experimental
biologists typically combine results from a variety of bioinformatics resources
in an attempt to portray and understand cancer relevance. As a means to
simplify and strengthen the support for this endeavor, we have developed
oncoEnrichR, a flexible bioinformatics tool that allows cancer researchers to
comprehensively interrogate a given gene list along multiple facets of cancer
relevance. oncoEnrichR differs from general gene set analysis frameworks
through the integration of an extensive set of prior knowledge specifically
relevant for cancer, including ranked gene-tumor type associations,
literature-supported proto-oncogene and tumor suppressor gene annotations,
target druggability data, regulatory interactions, synthetic lethality
predictions, as well as prognostic associations, gene aberrations, and
co-expression patterns across tumor types. The software produces a structured
and user-friendly analysis report as its main output, where versions of all
underlying data resources are explicitly logged, the latter being a critical
component for reproducible science. We demonstrate the usefulness of
oncoEnrichR through interrogation of two candidate lists from proteomic and
CRISPR screens. oncoEnrichR is freely available as a web-based workflow hosted
by the Galaxy platform (https://oncotools.elixir.no), and can also be accessed
as a stand-alone R package (https://github.com/sigven/oncoEnrichR).


  Several external sources may change the river's geomorphology, hydrodynamic
characteristics, and water quality. Therefore, having a comprehensive
assessment of the rivers' hydrodynamic would bring the preliminary information
to reduce the human-made structures' short- and long-term effects on a
different aspect of a region. An automated GIS toolbox is developed in the
present study, which interpolates the hydrodynamic data from the HEC-RAS
hydrodynamic model. A suitable fish habitat model is designed in the GIS
toolbox, which uses the hydrodynamic interpolation result to assess the
suitable habitat for aquatic species. In the study, different fish species'
physical properties were considered to be used in conjunction with the
hydrodynamic model to generate the habitat suitability raster datasets. This
study provides the designers and scholars with a graphical interface of the GIS
platform to model the potential migratory maps of the fish species under
different hydrodynamic conditions.


  Genes have specific functional roles, however, since they are dependent on
each other, they can play a structural role within a network structure of their
interactions. In this study, we analyze the structure of the gene interaction
network and detect the most contributing genes through the random matrix
theory. Specifically, we compare the interaction network of essential and
nonessential genes of the yeast Saccharomyces cerevisiae. Most remarkably, this
well-established combined framework by measuring the node participation ratio
$(NPR)$ index helps detect important genes, which control the insightful
structural patterns in the underlying networks. Results indicate that the
essential genes have higher values of $NPR$ rather than the nonessential ones
which means that they have the most contribution to the network structure. It
is worth mentioning that among all essential genes, the $NPR$ value of 5
significant ones is considerably higher than the other essential genes, and
also the same is for 15 significant nonessential genes compared to the others.
Thus, the significant essential genes strongly manage their network structure,
while the significant nonessential genes, besides their global contributions,
have weak effects on their network structure. Most strikingly, these genes
existing in a limited number of structural patterns are responsible for the
specific bioprocesses which are the signature of their networks.


  We introduce a new model for plant metapopulations with a seed bank
component, living in a fragmented environment in which local extinction events
are frequent. This model is an intermediate between population dynamics models
with a seed bank component, based on the classical Wright-Fisher model, and
Stochastic Patch Occupancy Models (SPOMs) used in metapopulation ecology. Its
main feature is the use of "ghost" individuals, which can reproduce but with a
very strong selective disadvantage against "real" individuals, to artificially
ensure a constant population size. We show the existence of an extinction
threshold above which persistence of the subpopulation of "real" individuals is
not possible, and investigate how the seed bank characteristics affect this
extinction threshold. We also show the convergence of the model to a SPOM under
an appropriate scaling, bridging the gap between individual-based models and
occupancy models.


  Pathogen introduction in plant communities can cause serious impact and
biodiversity losses that may take long time to manage and restore. Effective
control of epidemic spreading in the wild is a problem of paramount importance,
because of its implications in conservation and potential economic losses.
Understanding the mechanisms that hinder pathogen's propagation is, therefore,
crucial. Usual modelization approaches in epidemic spreading are based in
compartmentalized models, without keeping track of pathogen's concentrations
during spreading. In this contribution we present and fully analyze a dynamical
model for plant epidemic spreading based on pathogen's abundances. The model,
which is defined on top of network substrates, is amenable to a deep
mathematical analysis in the absence of a limit in the amount of pathogen a
plant can tolerate before dying. In the presence of such death threshold, we
observe that the fraction of dead plants peaks at intermediate values of
network's connectivity, and mortality decreases for large average degrees. We
discuss the implications of our results as mechanisms to halt infection
propagation.


  The received wisdom on how activity affects energy expenditure is that the
more activity is undertaken, the more calories will have been burned by the end
of the day. Yet traditional hunter-gatherers, who lead physically hard lives,
burn no more calories each day than western populations living in labour-saving
environments. Indeed, there is now a wealth of data, both for humans and other
animals, demonstrating that long-term lifestyle changes involving increases in
exercise or other physical activities do not result in commensurate increases
in daily energy expenditure (DEE). This is because humans and other animals
exhibit a degree of energy compensation at the organismal level, ameliorating
some of the increases in DEE that would occur from the increased activity by
decreasing the energy expended on other biological processes. And energy
compensation can be sizable, reaching many hundreds of calories in humans. But
the processes that are downregulated in the long-term to achieve energy
compensation are far from clear, particularly in humans. We do not know how
energy compensation is achieved. My review here of the literature on relevant
exercise intervention studies, for both humans and other species, indicates
conflict regarding the role that basal metabolic rate (BMR) or low level
activity such as fidgeting play, if any, particularly once changes in body
composition are factored out. In situations where BMR and low-level activity
are not major components of energy compensation, what then drives it? I discuss
how changes in mitochondrial efficiency and changes in circadian fluctuations
in BMR may contribute to our understanding of energy management. Currently
unexplored, these mechanisms and others may provide important insights into the
mystery of how energy compensation is achieved.


  In this study, we classify the COVID-19 anomalous diffusion in two categories
of countries based on the mean squared displacement (MSD) of daily new cases,
which includes the top four countries and four randomly selected countries in
terms of the total cases. The COVID-19 diffusion is a stochastic process, and
the daily new cases are regarded as the displacements of diffusive particles.
The diffusion environment of COVID-19 in each country is heterogeneous, in
which the underlying dynamic process is anomalous diffusion. The calculated MSD
is a power law function of time, and the power law exponent is not a constant
but varies with time. The power law exponents are estimated by using the
bi-exponential model and the long short-term memory network (LSTM). The
bi-exponential model frequently use in magnetic resonance imaging (MRI) can
quantify the power law exponent and make an easy prediction. The LSTM network
has much better accuracy than the bi-exponential model in predicting the power
law exponent. The LSTM network is more flexible and preferred to predict the
power law exponent, which is independent on the unique mathematical formula.
The diffusion process of COVID-19 can be classified based on the power law
exponent. More specific evaluation and suggestion can be proposed and submitted
to the government in order to control the COVID-19 diffusion.


  The next-generation sequencers such as Illumina and SOLiD platforms generate
a large amount of data, commonly above 10 Gigabytes of text files.
Particularly, the SOLiD platform allows the sequencing of multiple samples in a
single run, called multiplex run, through a tagging system called Barcode. This
feature requires a computational process for separation of the data sample
because the sequencer provides a mixture of all samples in a single output.
This process must be secure to avoid any harm that may scramble further
analysis. In this context, realized the need to develop a probabilistic model
capable of assigning a degree of confidence in the marking system used in
multiplex sequencing. The results confirmed the adequacy of the model obtained,
which allows, among other things, to guide a process of filtering the data and
evaluation of the sequencing protocol used.


  Coupled chemical interactions in a well-mixed solution are commonly
formalized as chemical reaction networks (CRNs). However, despite the
widespread use of CRNs in the natural sciences, the range of computational
behaviors exhibited by CRNs is not well understood. Here we study the following
problem: what functions $f:\mathbb{R}^k \to \mathbb{R}$ can be computed by a
CRN, in which the CRN eventually produces the correct amount of the "output"
molecule, no matter the rate at which reactions proceed? This captures a
previously unexplored, but very natural class of computations: for example, the
reaction $X_1 + X_2 \to Y$ can be thought to compute the function $y =
\min(x_1, x_2)$. Such a CRN is robust in the sense that it is correct no matter
the kinetic model of chemistry, so long as it respects the stoichiometric
constraints.
  We develop a reachability relation based on "what could happen" if reaction
rates can vary arbitrarily over time. We define *stable computation*
analogously to probability 1 computation in distributed computing, and connect
it with a seemingly stronger notion of rate-independent computation based on
convergence under a wide class of generalized rate laws. We also consider the
"dual-rail representation" that can represent negative values as the difference
of two concentrations and allows the composition of CRN modules. We prove that
a function is rate-independently computable if and only if it is piecewise
linear (with rational coefficients) and continuous (dual-rail representation),
or non-negative with discontinuities occurring only when some inputs switch
from zero to positive (direct representation). The many contexts where
continuous piecewise linear functions are powerful targets for implementation,
combined with the systematic construction we develop for computing these
functions, demonstrate the potential of rate-independent chemical computation.


  Efficient testing and vaccination protocols are critical aspects of epidemic
management. To study the optimal allocation of limited testing and vaccination
resources in a heterogeneous contact network of interacting susceptible,
recovered, and infected individuals, we present a degree-based testing and
vaccination model for which we use control-theoretic methods to derive optimal
testing and vaccination policies. Within our framework, we find that optimal
intervention policies first target high-degree nodes before shifting to
lower-degree nodes in a time-dependent manner. Using such optimal policies, it
is possible to delay outbreaks and reduce incidence rates to a greater extent
than uniform and reinforcement-learning-based interventions, particularly on
certain scale-free networks.


  The twin crises of climate change and biodiversity loss define a strong need
for functional diversity monitoring. While the availability of high-quality
ecological monitoring data is increasing, the quantification of functional
diversity so far requires the identification of species traits, for which data
is harder to obtain. However, the traits that are relevant for the ecological
function of a species also shape its performance in the environment and hence
should be reflected indirectly in its spatio-temporal distribution. Thus it may
be possible to reconstruct these traits from a sufficiently extensive
monitoring dataset. Here we use diffusion maps, a deterministic and de-facto
parameter-free analysis method, to reconstruct a proxy representation of the
species' traits directly from monitoring data and use it to estimate functional
diversity. We demonstrate this approach both with simulated data and real-world
phytoplankton monitoring data from the Baltic sea. We anticipate that wider
application of this approach to existing data could greatly advance the
analysis of changes in functional biodiversity.


  Cell division, aging, and stress recovery triggers spatial reorganization of
cellular components in the cytoplasm, including membrane bound organelles, with
molecular changes in their compositions and structures. However, it is not
clear how these events are coordinated and how they integrate with regulation
of molecular crowding. We use the budding yeast Saccharomyces cerevisiae as a
model system to study these questions using recent progress in optical
fluorescence microscopy and crowding sensing probe technology. We used a
F\"{o}rster Resonance Energy Transfer (FRET) based sensor, illuminated by
confocal microscopy for high throughput analyses and Slimfield microscopy for
single-molecule resolution, to quantify molecular crowding. We determine
crowding in response to cellular growth of both mother and daughter cells, in
addition to osmotic stress, and reveal hot spots of crowding across the bud
neck in the burgeoning daughter cell. This crowding might be rationalized by
the packing of inherited material, like the vacuole, from mother cells. We
discuss recent advances in understanding the role of crowding in cellular
regulation and key current challenges and conclude by presenting our recent
advances in optimizing FRET-based measurements of crowding whilst
simultaneously imaging a third color, which can be used as a marker that labels
organelle membranes. Our approaches can be combined with synchronised cell
populations to increase experimental throughput and correlate molecular
crowding information with different stages in the cell cycle.


  The Coronavirus 2019 (Covid-19) pandemic caused by the SARS-CoV-2 virus
represents an unprecedented crisis for our planet. It is a bane of the \"uber
connected world that we live in that this virus has affected almost all
countries and caused mortality and economic upheaval at a scale whose effects
are going to be felt for generations to come. While we can all be buoyed at the
pace at which vaccines have been developed and brought to market, there are
still challenges ahead for all countries to get their populations vaccinated
equitably and effectively. This paper provides an overview of ongoing
immunization efforts in various countries. In this early draft, we have
identified a few key factors that we use to review different countries' current
COVID-19 immunization strategies and their strengths and draw conclusions so
that policymakers worldwide can learn from them. Our paper focuses on processes
related to vaccine approval, allocation and prioritization, distribution
strategies, population to vaccine ratio, vaccination governance, accessibility
and use of digital solutions, and government policies. The statistics and
numbers are dated as per the draft date [June 24th, 2021].


  The second wave of Covid-19 that started in mid-February 2021 in Mumbai is
now subsiding. Increasingly the focus amongst the policy makers and general
public is on the potential third wave. Due to uncertainties regarding emergence
of new variants and reinfections, instead of projecting our best guess
scenario, in this report we conduct an extensive scenario analysis for Mumbai
and track peak fatalities in the coming months in each of these scenarios. Our
key conclusions are - As per our model, about 80% of Mumbai population has been
exposed to Covid-19 by June 1, 2021. Under the assumption that all who are
exposed have immunity against further infection, it is unlikely that Mumbai
would see a large third wave. It is the reinfections that may lead to a large
wave. Reinfections could occur because of declining antibodies amongst the
infected as well as by variants that can break through the immunity provided by
prior infections. Even under a reasonably pessimistic scenario we observe the
resulting peak to be no larger than that under the second wave. We further
observe that under the scenario where the reinfections are mild so that they
affect the fatality figures negligibly, where the new variants (beyond the
existing delta variant) have a mild impact, as the city opens up, we observe a
small wave in the coming months. However, if by then the vaccine coverage is
extensive, this wave will be barely noticeable. We also plot $R_t$, the
infection growth rate at time $t$, and highlight some interesting observations.


  Homochirality, i.e. the dominance across all living matter of one enantiomer
over the other among chiral molecules, is thought to be a key step in the
emergence of life. Building on ideas put forward by Frank and many others, we
proposed recently one such mechanism in G. Laurent et al., PNAS (2021) based on
the properties of large out of equilibrium chemical networks. We showed that in
such networks, a phase transition towards an homochiral state is likely to
occur as the number of chiral species in the system becomes large or as the
amount of free energy injected into the system increases. This paper aims at
clarifying some important points in that scenario, not covered by our previous
work. We first analyze the various conventions used to measure chirality,
introduce the notion of chiral symmetry of a network, and study its
implications regarding the relative chiral signs adopted by different groups of
molecules. We then propose a generalization of Frank's model for large chemical
networks, which we characterize completely using methods of random matrices.
This analysis can be extended to sparse networks, which shows that the
emergence of homochirality is a robust transition.


  Cortical actin networks are highly dynamic and play critical roles in shaping
the mechanical properties of cells. The actin cytoskeleton undergoes
significant reorganization over the course of the cell cycle, when cortical
actin transitions between open patched meshworks, homogeneous distributions,
and aligned bundles. Several types of myosin motor proteins, characterized by
different kinetic parameters, have been involved in this reorganization of
actin filaments. Given the limitations in studying the interactions of actin
with myosin in vivo, we propose stochastic agent-based model simulations and
develop a set of data analysis measures to assess how myosin motor proteins
mediate various actin organizations. In particular, we identify individual
motor parameters, such as motor binding rate and step size, that generate actin
networks with different levels of contractility and different patterns of
myosin motor localization. In simulations where two motor populations with
distinct kinetic parameters interact with the same actin network, we find that
motors may act in a complementary way, by tuning the actin network
organization, or in an antagonistic way, where one motor emerges as dominant.
This modeling and data analysis framework also uncovers parameter regimes where
spatial segregation between motor populations is achieved. By allowing for
changes in kinetic rates during the actin-myosin dynamic simulations, our work
suggests that certain actin-myosin organizations may require additional
regulation beyond mediation by motor proteins in order to reconfigure the
cytoskeleton network on experimentally-observed timescales.


  Conformational dynamics of biomolecules are of fundamental importance for
their function. Single-molecule F\"orster Resonance Energy Transfer (smFRET) is
a powerful approach to inform on the structure and the dynamics of labeled
molecules. If the dynamics occur on the sub-millisecond timescale, capturing
and quantifying conformational dynamics can be challenging by intensity-based
smFRET. Multiparameter fluorescence detection (MFD) addresses this challenge by
simultaneously registering intensities and fluorescence lifetimes. Together,
the mean donor fluorescence lifetime and the fluorescence intensities inform on
the variance, and the mean FRET efficiency tells the conformational dynamics.
Here, we present a general framework that relates average fluorescence
lifetimes and intensities in smFRET counting histograms. Using this framework,
we show how to compute parametric relations (FRET-lines) of these observables
that facilitate a graphical interpretation of experimental data, can be used to
test models, identify conformational states, resolve the connectivity of
states, and can be applied to unstructured systems to infer properties of
polymer chains or study fast protein folding. To simplify the graphical
analysis of complex kinetic networks, we derive a moment-based representation
of the experimental data and show how to decouple the motion of the
fluorescence labels from the conformational dynamics of the biomolecule.


  Single-molecule F\"orster Resonance Energy Transfer (smFRET) is ideally
suited to resolve the dynamics of biomolecules. A significant challenge to date
is capturing and quantifying the exchange between multiple conformational
states, mainly when these dynamics occur on the sub-millisecond timescale. Many
methods for the quantitative analysis are challenged if more than two states
are involved, and the appropriate choice of the number of states in the kinetic
network is difficult. An additional complication arises if dynamically active
molecules coexist with solely static molecules in similar conformational
states. To address these problems, we developed an integrative analysis
framework that combines the information from FRET-lines, time-correlated single
photon counting, and fluorescence correlation spectroscopy. While individually,
these methodologies provide ambiguous results for quantitatively characterize
the dynamics in complex kinetic networks, the global analysis enables accurate
determination of the number of states, their kinetic connectivity, and the
transition rate constants. To challenge the potential of smFRET to study
multi-state kinetic networks, we benchmark the framework using synthetic data
of three-state systems with different kinetic connectivity and interconversion
rates.


  By circumventing the optical diffraction limit, super-resolved fluorescence
microscopies enable the study of larger cellular structures and molecular
assemblies. However, fluorescence nanoscopy currently lacks the spatiotemporal
resolution to resolve distances on the size of individual molecules and reveal
the conformational fine structure and dynamics of molecular complexes. Here we
establish FRET nanoscopy by combining colocalization STED microscopy with
multiparameter FRET spectroscopy. We simultaneously localize donor and acceptor
dyes of single FRET pairs with nanometer resolution and quantitatively measure
intramolecular distances with sub-nanometer precision over a large dynamic
range. While FRET provides isotropic 3D distance information, colocalization
measures the projected distance onto the image plane. The combined information
allows us to directly determine its 3D orientation using Pythagoras's theorem.
Studying two DNA model systems and the human guanylate binding protein hGBP1,
we demonstrate that FRET nanoscopy unravels the interplay between their spatial
organization and local molecular conformation in a complex environment.


  Global change is reshaping ecosystems and societies. Strategic choices that
were best yesterday may be sub-optimal tomorrow; and environmental conditions
that were once taken for granted may soon cease to exist. In this setting, how
people choose behavioral strategies has important consequences for
environmental dynamics. Economic and evolutionary theories make similar
predictions for strategic behavior in a static environment, even though one
approach assumes perfect rationality and the other assumes no cognition
whatsoever; but predictions differ in a dynamic environment. Here we explore a
middle ground between economic rationality and evolutionary myopia. Starting
from a population of myopic agents, we study the emergence of a new type that
forms environmental forecasts when making strategic decisions. We show that
forecasting types can have an advantage in changing environments, even when the
act of forecasting is costly. Forecasting types can invade but not overtake the
population, producing a stable coexistence with myopic types. Moreover,
forecasters provide a public good by reducing the amplitude of environmental
oscillations and increasing mean payoff to forecasting and myopic types alike.
We interpret our results for understanding the evolution of different modes of
decision-making. And we discuss implications for the management of
environmental systems of great societal importance.


  A novel model is here introduced for the SOC change index defined as the
normalized difference between the actual Soil Organic Carbon and the value
assumed at an initial reference year. It is tailored on the RothC carbon model
dynamics and assumes as baseline the value of the SOC equilibrium under
constant environmental conditions. A sensitivity analysis is performed to
evaluate the response of the model to changes of temperature, Net Primary
Production (NPP), and land use soil class (forest, grassland, arable). A
non-standard monthly time-stepping procedure has been proposed to approximate
the SOC change index in the Alta Murgia National Park, a protected area in the
Italian Apulia region, selected as test site. In the case of arable class, the
SOC change index exhibits a negative trend which can be inverted by a suitable
organic fertilization program here proposed.


  We investigate the macroscopic time evolution and stationary states of a mean
field generalized contact process in $\mathbb{R}^d$. The model is described by
a coupled set of nonlinear integral-differential equations. It was inspired by
a model of neurons with discrete voltages evolving by a stochastic integrate
and fire mechanism. We obtain a complete solution in the spatially uniform case
and partial solutions in the general case. The system has one or more fixed
points and also traveling wave solutions.


  The discovery of disease biomarkers from gene expression data has been
greatly advanced by feature selection (FS) methods, especially using ensemble
FS (EFS) strategies with perturbation at the data level (i.e., homogeneous,
Hom-EFS) or method level (i.e., heterogeneous, Het-EFS). Here we proposed a
Hybrid EFS (Hyb-EFS) design that explores both types of perturbation to improve
the stability and the predictive power of candidate biomarkers. With this,
Hyb-EFS aims to disrupt associations of good performance with a single dataset,
single algorithm, or a specific combination of both, which is particularly
interesting for better reproducibility of genomic biomarkers. We investigated
the adequacy of our approach for microarray data related to four types of
cancer, carrying out an extensive comparison with other ensemble and single FS
approaches. Five FS methods were used in our experiments: Wx, Symmetrical
Uncertainty (SU), Gain Ratio (GR), Characteristic Direction (GeoDE), and
ReliefF. We observed that the Hyb-EFS and Het-EFS approaches attenuated the
large performance variation observed for most single FS and Hom-EFS across
distinct datasets. Also, the Hyb-EFS improved upon the stability of the Het-EFS
within our domain. Comparing the Hyb-EFS and Het-EFS composed of the
top-performing selectors (Wx, GR, and SU), our hybrid approach surpassed the
equivalent heterogeneous design and the best Hom-EFS (Hom-Wx). Interestingly,
the rankings produced by our Hyb-EFS reached greater biological plausibility,
with a notably high enrichment for cancer-related genes and pathways. Thus, our
experiments suggest the potential of the proposed Hybrid EFS design in
discovering candidate biomarkers from microarray data. Finally, we provide an
open-source framework to support similar analyses in other domains, both as a
user-friendly application and a plain Python package.


  The prevalence of multicellular organisms is due in part to their ability to
form complex structures. How cells pack in these structures is a fundamental
biophysical issue, underlying their functional properties. However, much
remains unknown about how cell packing geometries arise, and how they are
affected by random noise during growth - especially absent developmental
programs. Here, we quantify the statistics of cellular neighborhoods of two
different multicellular eukaryotes: lab-evolved "snowflake" yeast and the green
alga $Volvox~carteri$. We find that despite large differences in cellular
organization, the free space associated with individual cells in both organisms
closely fits a modified gamma distribution, consistent with maximum entropy
predictions originally developed for granular materials. This 'entropic'
cellular packing ensures a degree of predictability despite noise, facilitating
parent-offspring fidelity even in the absence of developmental regulation.
Together with simulations of diverse growth morphologies, these results suggest
that gamma-distributed cell neighborhood sizes are a general feature of
multicellularity, arising from conserved statistics of cellular packing.


  A deterministic two-species predator-prey model with prey herd behavior is
considered incorporating mutual interference and the effect of fear. We provide
guidelines to the dynamical analysis of biologically feasible equilibrium
points. We give conditions for the existence of some local and global
bifurcations at the coexistence equilibrium. We also show that fear can induce
extinction of the prey population from a coexistence zone in finite time. Our
numerical simulations reveal that varying the strength of fear of predators
with suitable choice of parameters can stabilize and destabilize the
coexistence equilibrium solutions of the model. Additionally, we discuss the
outcome of introducing a constant harvesting effort to the predator population
in terms of changing the dynamics of the system, in particular, from finite
time extinction to stable coexistence. Furthermore, we perform extensive
numerical experiments to visualize the dynamical behavior of the model and
substantiate the results we obtained.


  An information processing circuit must be flexible to perform multiple tasks:
biochemical information processing circuits in living systems can have multiple
functions, requiring them to perform distinct computational tasks in different
scenarios. Here, we describe a tunable dynamical system that performs different
logical operations (such as AND, OR, and XOR) depending on which basin of
attraction the dynamics resides in. We analyze the robustness of these
dynamically switchable logic gates, and characterize the trade-off between
reliability and efficiency in their function.


  Motivated by recent epidemic outbreaks, including those of COVID-19, we solve
the canonical problem of calculating the dynamics and likelihood of extensive
outbreaks in a population within a large class of stochastic epidemic models
with demographic noise, including the Susceptible-Infected-Recovered (SIR)
model and its general extensions. In the limit of large populations, we compute
the probability distribution for all extensive outbreaks, including those that
entail unusually large or small (extreme) proportions of the population
infected. Our approach reveals that, unlike other well-known examples of rare
events occurring in discrete-state stochastic systems, the statistics of
extreme outbreaks emanate from a full continuum of Hamiltonian paths, each
satisfying unique boundary conditions with a conserved probability flux.


  Global strategies to contain a pandemic, such as social distancing and
protective measures, are designed to reduce the overall transmission rate
between individuals. Despite such measures, essential institutions, including
hospitals, schools, and food producing plants, remain focal points of local
outbreaks. Here we develop a model for the stochastic outbreak dynamics in such
local communities. We derive analytical expressions for the probability of
containment of the outbreak, which is complementary to the probability of
seeding a deterministically growing epidemic. This probability depends on the
statistics of the intra-community contact network and the initial conditions,
in particular, on the contact degree of patient zero. Based on this model, we
suggest surveillance protocols by which individuals are tested proportionally
to their degree in the contact network. We characterize the efficacy of
contact-based protocols as a function of the epidemiological and the contact
network parameters, and show numerically that such protocols outperform random
testing.


  The technology to generate Spatially Resolved Transcriptomics (SRT) data is
rapidly being improved and applied to investigate a variety of biological
tissues. The ability to interrogate how spatially localised gene expression can
lend new insight to different tissue development is critical, but the
appropriate tools to analyse this data are still emerging. This chapter reviews
available packages and pipelines for the analysis of different SRT datasets
with a focus on identifying spatially variable genes (SVGs) alongside other
aims, while discussing the importance of and challenges in establishing a
standardised 'ground truth' in the biological data for benchmarking.


  In several previous works, I presented the mirror symmetry in the set of
protein amino acids, expressed through the number of atoms. Here, however, the
same thing is shown but over the number of nucleons and molecules mass.
Compared to the previous version of the paper, minimal changes have been made,
and Display 2 as well as Figures 3 and 4 have been added.


  The transcription factors, such as activators and repressors, can interact
with the promoter of gene either in a competitive or non-competitive way. In
this paper, we construct a stochastic model with non-competitive
transcriptional regulatory architecture and develop an analytical theory that
re-establishes the experimental results with an improved data fitting. The
analytical expressions in the theory allow us to study the nature of the system
corresponding to any of its parameters, and hence enable us to find out the
factors that govern the regulation of gene expression for that architecture. We
notice that, along with transcriptional reinitiation and repressors, there are
other parameters that can control the noisiness of this network. We also
observe that, the Fano factor (at mRNA level) varies from sub-Poissonian regime
to superPoissonian regime. In addition to the aforementioned properties, we
observe some anomalous characteristics of the Fano factor (at mRNA level) and
that of the variance of protein at lower activator concentrations in presence
of repressor molecules. This model is useful to understand the architecture of
interactions which may buffer the stochasticity inherent to gene transcription.


  Inference of network-like evolutionary relationships between species from
genomic data must address the interwoven signals from both gene flow and
incomplete lineage sorting. The heavy computational demands of standard
approaches to this problem severely limit the size of datasets that may be
analyzed, in both the number of species and the number of genetic loci. Here we
provide a theoretical pointer to more efficient methods, by showing that logDet
distances computed from genomic-scale sequences retain sufficient information
to recover network relationships in the level-1 ultrametric case. This result
is obtained under the Network Multispecies Coalescent model combined with a
mixture of General Time-Reversible sequence evolution models across individual
gene trees, but does not depend on partitioning sequences by genes. Thus under
standard stochastic models statistically justifiable inference of network
relationships from sequences can be accomplished without consideration of
individual genes or gene trees.


  We investigate melting transition of DNA sequences embedded in a Langevin
fluctuation-dissipation thermal bath. Torsional effects are considered by a
twist angle $\varphi$ between neighboring base pairs stacked along the molecule
backbone. Our simulation results show that the increase of twist angle
translates linearly the melting temperature with a positive slope. After the so
called equilibrium angle $\varphi_\mathrm{eq}$, the DNA chain becomes very
rigid against opening and accordingly very high temperatures are required to
initiate the melting process. In such cases however, the biofunctionality of
DNA is destroyed before so that the observed in our model melting process
becomes biologically irrelevant. We believe that the outcome of this survey
would deeper understanding of the interplay between DNA twisting and melting
transition for precise control of DNA behavior.


  We develop theoretical equivalences between stochastic and deterministic
models for populations of individual cells stratified by age. Specifically, we
develop a hierarchical system of equations describing the full dynamics of an
age-structured multi-stage Markov process for approximating cell cycle time
distributions. We further demonstrate that the resulting mean behaviour is
equivalent, over large timescales, to the classical McKendrick-von Foerster
integro-partial differential equation. We conclude by extending this framework
to a spatial context, facilitating the modelling of travelling wave phenomena
and cell-mediated pattern formation. More generally, this methodology may be
extended to myriad reaction-diffusion processes for which the age of
individuals is relevant to the dynamics.


  The rapidly emerging field of deep learning-based computational pathology has
demonstrated promise in developing objective prognostic models from histology
whole slide images. However, most prognostic models are either based on
histology or genomics alone and do not address how histology and genomics can
be integrated to develop joint image-omic prognostic models. Additionally
identifying explainable morphological and molecular descriptors from these
models that govern such prognosis is of interest. We used multimodal deep
learning to integrate gigapixel whole slide pathology images, RNA-seq
abundance, copy number variation, and mutation data from 5,720 patients across
14 major cancer types. Our interpretable, weakly-supervised, multimodal deep
learning algorithm is able to fuse these heterogeneous modalities for
predicting outcomes and discover prognostic features from these modalities that
corroborate with poor and favorable outcomes via multimodal interpretability.
We compared our model with unimodal deep learning models trained on histology
slides and molecular profiles alone, and demonstrate performance increase in
risk stratification on 9 out of 14 cancers. In addition, we analyze morphologic
and molecular markers responsible for prognostic predictions across all cancer
types. All analyzed data, including morphological and molecular correlates of
patient prognosis across the 14 cancer types at a disease and patient level are
presented in an interactive open-access database
(http://pancancer.mahmoodlab.org) to allow for further exploration and
prognostic biomarker discovery. To validate that these model explanations are
prognostic, we further analyzed high attention morphological regions in WSIs,
which indicates that tumor-infiltrating lymphocyte presence corroborates with
favorable cancer prognosis on 9 out of 14 cancer types studied.


  A zombie apocalypse is one pandemic that would likely be worse than anything
humanity has ever seen. However, despite the mechanisms for zombie uprisings in
pop culture, it is unknown whether zombies, from an evolutionary point of view,
can actually rise from the dead. To provide insight into this unknown, we
created a mathematical model that predicts the trajectory of human and zombie
populations during a zombie apocalypse. We parameterized our model according to
the demographics of the US, the zombie literature, and then conducted an
evolutionary invasion analysis to determine conditions that permit the
evolution of zombies. Our results indicate a zombie invasion is theoretically
possible, provided there is a sufficiently large ratio of transmission rate to
the zombie death rate. While achieving this ratio is uncommon in nature, the
existence of zombie ant fungus illustrates it is possible and thereby suggests
that a zombie apocalypse among humans could occur.


  Novel diagnostic and therapeutic radiopharmaceuticals are increasingly
becoming a central part of personalized medicine. Continued innovation in the
development of new radiopharmaceuticals is key to sustained growth and
advancement of precision medicine. Artificial intelligence (AI) has been used
in multiple fields of medicine to develop and validate better tools for patient
diagnosis and therapy, including in radiopharmaceutical design. In this review,
we first discuss common in silico approaches and focus on their utility and
challenges in radiopharmaceutical development. Next, we discuss the practical
applications of in silico modeling in design of radiopharmaceuticals in various
diseases.


  The spread of pandemic viruses and invasive species can be catastrophic for
human societies and natural ecosystems. SARS-CoV-2 demonstrated that the speed
of our response is critical, as each day of delay permitted exponential growth
and dispersion of the virus. Here we propose a global Nucleic Acid Observatory
(NAO) to monitor the relative frequency of everything biological through
comprehensive metagenomic sequencing of waterways and wastewater. By searching
for divergences from historical baseline frequencies at sites throughout the
world, NAO could detect any virus or invasive organism undergoing exponential
growth whose nucleic acids end up in the water, even those previously unknown
to science. Continuously monitoring nucleic acid diversity would provide us
with universal early warning, obviate subtle bioweapons, and generate a wealth
of sequence data sufficient to transform ecology, microbiology, and
conservation. We call for the immediate construction of a global NAO to defend
and illuminate planetary health.


  While many good textbooks are available on Protein Structure, Molecular
Simulations, Thermodynamics and Bioinformatics methods in general, there is no
good introductory level book for the field of Structural Bioinformatics. This
book aims to give an introduction into Structural Bioinformatics, which is
where the previous topics meet to explore three dimensional protein structures
through computational analysis. We provide an overview of existing
computational techniques, to validate, simulate, predict and analyse protein
structures. More importantly, it will aim to provide practical knowledge about
how and when to use such techniques. We will consider proteins from three major
vantage points: Protein structure quantification, Protein structure prediction,
and Protein simulation & dynamics.
  The main emphasis of this work is to provide a background on experimental
techniques for protein structure determination. The focus is set on X-ray
crystallography and Nuclear Magnetic Resonance spectroscopy (NMR), which are by
far the main methods used to determine the structure of soluble proteins. We
will also introduce cryogenic Electron Microscopy (cryo-EM) and electron
diffraction which are more suited to analyze membrane proteins and larger
protein complexes. At the end, more qualitative techniques are summarized that
are used to obtain insight on the overall structure and dynamics of proteins.
Note that this introduction to protein structure determination aims at
familiarizing the reader to different experimental techniques, their benefits
and bottlenecks, but that a thorough mathematical and technical description of
the concept is beyond the scope of this work.


  The spreading of bacterial populations is central to processes in
agriculture, the environment, and medicine. However, existing models of
spreading typically focus on cells in unconfined settings--despite the fact
that many bacteria inhabit complex and crowded environments, such as soils,
sediments, and biological tissues/gels, in which solid obstacles confine the
cells and thereby strongly regulate population spreading. Here, we develop an
extended version of the classic Keller-Segel model of bacterial spreading that
incorporates the influence of confinement in promoting both cell-solid and
cell-cell collisions. Numerical simulations of this extended model demonstrate
how confinement fundamentally alters the dynamics and morphology of spreading
bacterial populations, in good agreement with recent experimental results. In
particular, with increasing confinement, we find that cell-cell collisions
increasingly hinder the initial formation and the long-time propagation speed
of chemotactic pulses. Moreover, also with increasing confinement, we find that
cellular growth and division plays an increasingly dominant role in driving
population spreading--eventually leading to a transition from chemotactic
spreading to growth-driven spreading via a slower, jammed front. This work thus
provides a theoretical foundation for further investigations of the influence
of confinement on bacterial spreading. More broadly, these results help to
provide a framework to predict and control the dynamics of bacterial
populations in complex and crowded environments.


  Summary: We present several recent improvements to minimap2, a versatile
pairwise aligner for nucleotide sequences. Now minimap2 v2.22 can more
accurately map long reads to highly repetitive regions and align through
insertions or deletions up to 100kb by default, addressing major weakness in
minimap2 v2.18 or earlier.
  Availability and implementation: https://github.com/lh3/minimap2


  The ancestral sequence reconstruction problem is the inference, back in time,
of the properties of common sequence ancestors from measured properties of
contemporary populations. Standard algorithms for this problem assume
independent (factorized) evolution of the characters of the sequences, which is
generally wrong (e.g. proteins and genome sequences). In this work, we have
studied this problem for sequences described by global co-evolutionary models,
which reproduce the global pattern of cooperative interactions between the
elements that compose it. For this, we first modeled the temporal evolution of
correlated real valued characters by a multivariate Ornstein-Uhlenbeck process
on a finite tree. This represents sequences as Gaussian vectors evolving in a
quadratic potential, who describe selection forces acting on the evolving
entities. Under a Bayesian framework, we developed a reconstruction algorithm
for these sequences and obtained an analytical expression to quantify the
quality of our estimation. We extend this formalism to discrete valued
sequences by applying our method to a Potts model. We showed that for both
continuous and discrete configurations, there is a wide range of parameters
where, to properly reconstruct the ancestral sequences, intra-species
correlations must be taken into account. We also demonstrated that, for
sequences with discrete elements, our reconstruction algorithm outperforms
traditional schemes based on independent site approximations.